多种群遗传算法在PBX本构模型参数识别中的应用

利用多种群并行结构对标准遗传算法SGA进行并行化处理,引入移民算子和精华种群形成多种群遗传算法MPGA,并设计了自适应交叉和变异概率对算法的收敛速度进行改进。结合ABAQUS软件和改进的多种群遗传算法,建立了材料本构模型参数识别方法。采用该方法对PBX炸药黏弹性损伤本构模型参数进行了模拟识别,并同基于标准遗传算法的参数识别方法进行了比较。结果证明,基于改进多种群遗传算法IMPGA的方法对克服算法未成熟收敛有显著的效果,识别结果更稳定。同时该方法的收敛速度更快,寻优能力更强,适合复杂非线性问题的优化,此方法可以被应用到其他材料本构模型的参数识别中。     

喷管对水下爆轰气泡形态和压力特征影响的实验研究

针对具有不同类型喷管的爆轰管在水下爆轰中形成的燃气射流问题,搭建了爆轰实验平台,研究了单次爆轰过程中尾部喷管对水下气泡形态与压力特征的影响。采用数字粒子图像测速技术对高速摄影机拍摄得到的气泡脉动图片进行流场可视化分析,得到各喷管工况下的气泡速度场。为了确认爆轰管内是否形成稳定爆轰波,并观察爆轰波在气液两相界面上的透反射特性,爆轰管尾部安装有2个动态压力传感器,同时在距离喷管一定距离处设置一个水下爆炸传感器,以监测水中传播的压力波。结果表明：扩张喷管工况下的气泡脉动过程与直喷管工况基本一致,但扩张喷管提高了燃气射流速度,气泡膨胀体积更大;因为燃气射流的持续性,收敛喷管工况下的气泡脉动过程具有明显差异,气泡膨胀体积较小,但气泡二次脉动时长相较于一次脉动时长衰减更小;扩张喷管提高了气泡脉动强度,扩张喷管工况下的气泡脉动压力与透射冲击波压力远大于直喷管工况下的气泡脉动压力与透射冲击波压力;收敛喷管工况下的气泡脉动压力与透射冲击波压力都较小,但收敛喷管燃气射流的持续性减缓了气泡脉动压力的衰减速度。相比于直喷管,扩张喷管工况下的气泡脉动时间、气泡脉动压力与透射冲击波压力都更大。收敛喷管工况下的气泡...     

水下爆炸气泡脉动测量及分析

水下爆炸气泡破坏效应是水中兵器的重要毁伤模式之一。为研究水下爆炸气泡脉动现象,建立了小当量水下爆炸实验系统,并进行了爆炸当量分别为0.125g、1.0g、3.375g和8g TNT的水下爆炸实验。采用球形PETN装药并中心起爆,产生球形对称的气泡和冲击波载荷,并利用高速摄像系统记录水下爆炸气泡脉动过程,以及布置压力传感器测量水中冲击波压力。实验获得了清晰的水下爆炸气泡脉动过程图像,得到了冲击波和气泡脉动压力曲线。对数字化图像进行判读,得到气泡脉动直径和周期。另外根据冲击波曲线测量了气泡脉动周期,对比分析了气泡脉动相关参数。结果表明,高速摄像数据测量的气泡直径与经验公式较接近,高速摄像测量的气泡周期与冲击波曲线测量的气泡脉动周期以及经验公式结果具有较好的一致性。本文提出的实验技术安全、经济、可靠,气泡脉动参数判读精确,满足水下爆炸气泡脉动研究需求。     

近场水下爆炸瞬态强非线性流固耦合无网格数值模拟研究

近场水下爆炸涉及多相流体的掺杂耦合以及结构的大变形、损伤和断裂等瞬态强非线性现象, 传统的网格算法在模拟近场水下爆炸时面临结构网格畸变、多相界面捕捉精度不足等难题, 鉴于此, 本文建立了完全无网格的近场水下爆炸冲击波和气泡全物理过程瞬态强非线性流固耦合动力学模型. 流体采用基于黎曼求解器的光滑粒子流体动力学(SPH)方法求解, 结构采用重构核粒子法(RKPM)求解, 并基于法向通量边界条件实现流固耦合. 为提高SPH对流场间断的求解精度, 引入黎曼问题思想并结合MUSCL重构算法, 为解决流场粒子体积变化剧烈导致的精度下降问题, 应用了自适应粒子分割与合并方法. 为模拟水下爆炸对结构造成的损伤断裂, 基于退化实体几何表述, 采用Lemaitre损伤算法, 建立了RKPM壳结构断裂损伤模型. 依据所建立的SPH-RKPM流固耦合模型, 对近场水下爆炸冲击波传播、气泡脉动与射流以及结构毁伤进行了模拟, 将得到的冲击波载荷、气泡演化以及结构响应与实验值和其他数值解对比, 验证了当前建立的SPH-RKPM流固耦合模型的有效性和精度, 并给出了水下爆炸载荷特性及其对结构的流固耦合毁伤机制与规律, 旨在为近场水下爆炸载荷预报提供理论和基础性技术支撑, 为毁伤威力评估和舰船防护结构设计提供参考.      

遗传算法改进及其在岩土参数反分析中的应用

本文的主要目的是开发基于实数编码的杂交遗传算法来识别土体的本构参数。该杂交遗传算法在经典遗传算法框架下开发,融合两个新开发的交叉算子,形成了一个新的杂交策略。为了保持种群的多样性,在算法中采用了一个动态随机变异算子。另外,为了提高算法收敛性,采用了一个基于混沌的局部搜索技术。分别基于室内试验和现场试验,通过识别土的本构参数来测试新算法的搜索能力和搜索效率。为了测试新开发算法的突出表现,特选用5种经典的随机类算法（遗传算法、粒子群算法、模拟退火算法、差分算法和蜂巢算法）,分析同样的案例进行比较。结果表明,在收敛速度和最优解的准确度方面,新改进的算法可以很好地处理岩土工程的参数反演。     

装药量及水深对水下爆炸气泡动态特性的影响

以气泡体积加速度模型为基础研究水下爆炸气泡运动的初始条件,采用MSC.DYTRAN 非线性 有限元软件,结合开发的定义流场初始条件与边界条件的子程序,研究水下爆炸气泡运动特性,包括气泡的脉 动、坍塌以及射流等运动特性,并将气泡脉动体积计算结果与实验及边界积分方法计算结果进行对比,验证了 有限元模型的正确性与有效性。以此为基础,得到初始水深、装药量与气泡的脉动体积、最大半径、周期以及 射流速度之间的关系,计算结果与经验公式具有较好的一致性。得到一些有规律性的曲线,可为相关水下爆 炸气泡动态特性研究提供参考。     

连续体结构非概率可靠性拓扑优化

基于非概率可靠性 指标的定义，考虑材料、几何及荷载大小的不确定性，提出以结构体积最小化为目标、具有 位移非概率可靠性约束的三维连续体拓扑优化数学模型. 采用目标性能方法对优化模型进行 转换，给出目标性能值的伴随法灵敏度分析算法，利用数学规划法实现优化问题的求解. 数 值算例验证了所提出优化模型的正确性及算法的有效性，并指出相对于确定性优化而言，非 概率可靠性拓扑优化能够给出在考虑不确定参数和荷载条件下更合理的材料分布.     

遗传算法改进及其在岩土参数反分析中的应用

本文的主要目的是开发基于实数编码的杂交遗传算法来识别土体的本构参数。该杂交遗传算法在经典遗传算法框架下开发,融合两个新开发的交叉算子,形成了一个新的杂交策略。为了保持种群的多样性,在算法中采用了一个动态随机变异算子。另外,为了提高算法收敛性,采用了一个基于混沌的局部搜索技术。分别基于室内试验和现场试验,通过识别土的本构参数来测试新算法的搜索能力和搜索效率。为了测试新开发算法的突出表现,特选用5种经典的随机类算法(遗传算法、粒子群算法、模拟退火算法、差分算法和蜂巢算法),分析同样的案例进行比较。结果表明,在收敛速度和最优解的准确度方面,新改进的算法可以很好地处理岩土工程的参数反演。     

基于改进遗传算法的桥梁结构损伤识别应用研究

提出了一种基于残余力向量和改进遗传算法的桥梁结构损伤识别方法。在无噪声的情况下,使用任意一阶模态数据,残余力向量法都能够对损伤进行准确定位。但是,振动测试数据中往往包含噪声,导致运用残余力向量法进行损伤识别完全不可行。考虑到这个问题,在常规模态分析的基础上,以节点的残余力向量构造目标函数,提出了一种用于遗传搜索优化的目标函数形式。利用改进遗传算法重点进行了噪声条件下的结构损伤定位和定量研究,并对遗传算法的参数选取问题进行了深入探讨。遗传算法的主要改进包括:采用浮点编码、采用基于标准化几何分布排名的选择策略、采用最优保存策略、采用算术交叉算子、采用自适应变异算子。最后,本文用一个连续梁桥模型进行了数值模拟,验证了所提出方法的有效性,并对方法应用中存在的一些问题进行了深入探讨。     

可压缩流场中气泡脉动数值模拟

在应用边界元方法对气泡动力学的研究中, 绝大多数模型是建立在不压缩势流理论基础之上, 针对可压缩流场中气泡运动特性的研究很少. 从波动方程出发, 分别在气泡运动前期和后期对波动方程进行简化, 得到气泡运动局部和全局简化方程, 采用双渐进方法对简化方程进行匹配, 提出了考虑流场可压缩性的非球状气泡运动模型. 该模型的计算结果与Prospertti 等的解析结果吻合很好, 气泡脉动最大半径和内部最大压力随气泡脉动逐渐减小. 基于该模型对比了自由场中药包爆炸考虑可压缩性与不考虑可压缩性的计算结果, 发现考虑可压缩性气泡射流速度较小, 随后基于该模型计算了刚性边界下气泡的运动特性.     

Attenuation analysis of hydraulic transients with laminar-turbulent flow alternations

An improved compound mathematical model is established to simulate the attenuation of hydraulic transients with laminar-turbulent alternations,which usually occur when the pipeline flow velocity fluctuates near the critical velocity.The laminar friction resistance and the turbulent friction resistance are considered respectively in this model by applying different resistance schemes to the characteristics method of fluid transient analysis.The hydraulic transients of the valve closing process are simulated using the model.A more reasonable attenuation of hydraulic transients is obtained.The accurate attenuation is more distinct than that obtained from the traditional mathematical model.The research shows that the hydraulic transient is a type of energy waves,and its attenuation is caused by the friction resistance.The laminar friction resistance is more important than the turbulent friction resistance if the flow velocity is smaller than the critical velocity.Otherwise the turbulent friction resistance is more important.The laminar friction resistance is important in the attenuation of hydraulic transients for the closing process.Thus,it is significant to consider the different resistances separately to obtain more accurate attenuation of hydraulic transients.     

Numerical computation of fluid pressure transients in pumping installations with air entrainment

In pumping installations such as sewage pumping stations, where gas content and air entrainment exist, the computation of fluid pressure transients in the pipelines becomes grossly inaccurate when constant wave speed and constant friction are assumed. A numerical model and computational procedure have been developed here to better compute the fluid pressure transient in a pipeline by including the effects of air entrainment and gas evolution characteristics of the transported fluid. Free and dissolved gases in the fluid and cavitation at the fluid vapour pressure are modelled. Numerical experiments show that entrained, entrapped or released gases amplify the pressure peak, increase surge damping and produce asymmetric pressure surges. The transient pressure shows a longer period for down-surge and a shorter period for up-surge. The up-surge is considerably amplified and the down-surge marginally reduced when compared with the gas-free case. These observations are consistent with the experimental observations of other investigators. Numerical experiments also show that the use of a variable loss factor in the pressure transient analysis produces marginally higher maximum and lower minimum pressure transients when compared with the constant-loss-factor model for pipelines where the pressures are above the fluid vapour pressure.     

Double method of characteristics to analyze hydraulic-thermal transients of pipeline flow

ＩｎｔｒｏｄｕｃｔｉｏｎＭｅｔｈｏｄｏｆｃｈａｒａｃｔｅｒｉｓｔｉｃｓ (ＭＯＣ)ｈａｓｓｔｉｌｌｂｅｅｎｔｈｅｐｏｗｅｒｆｕｌｔｏｏｌｔｏａｎａｌｙｚｅｈｙｄｒａｕｌｉｃｔｒａｎｓｉｅｎｔｏｆｐｉｐｅｌｉｎｅｆｌｏｗｓｉｎｃｅｔｈｅ 1 960ｓ.Ｓｅｌｄｏｍｈａｓｔｈｅｔｈｅｒｍａｌｔｒａｎｓｉｅｎｔｂｅｅｎｔａｋｅｎｉｎｔｏａｃｃｏｕｎｔ.Ｈｅａｔｄｕｅｔｏｆｒｉｃｔｉｏｎｖａｒｉｅｓｗｈｅｎｆｌｏｗｒａｔｅｃｈａｎｇｅｄｂｅｃａｕｓｅｏｆｖａｒｉｏｕｓｋｉｎｄｓｏｆｒｅａｓｏｎｓ,ｗｈｉｃｈｃａｕｓ…     

黏弹性管道水柱分离弥合水锤有限体积法模型

基于有限体积法二阶Godunov求解格式对黏弹性输水管道中水柱分离弥合现象进行建模和模拟研究. 在传统的弹性管道模型基础上考虑管道黏弹性效应的影响. 在瞬变流控制方程中引入管道黏弹性项和动态摩阻项, 采用有限体积法进行求解, 考虑压力修正系数来模拟自由气体对计算单元的影响, 同时为避免数值模拟结果产生虚假震荡引入斜率限制器MINMOD函数; 通过虚拟单元法进行边界构建, 实现了计算区域的统一计算. 将所建模型计算结果与已有模型结果、试验结果进行对比, 并对影响模型的各参数进行敏感性分析. 结果表明, 本文模型能够准确模拟出纯水锤、水柱分离弥合水锤两种情况下的瞬态压力变化, 均能与试验数据高度吻合; 与传统的特征线方法相比, 当库朗数C_r小于1时, 有限体积法二阶Godunov格式计算结果更准确、稳定; 在压力波动的衰减过程中, 黏弹性效应相比于管道摩阻起主导作用; 与弹性管道模型相比, 考虑管道黏弹性效应后可显著提高模拟结果的准确性, 尤其是压力波峰值的相对误差明显降低.      

基于瞬变流和遗传算法的管道泄漏辨识

输水系统管道的泄漏影响瞬变水击波波形的畸变和衰减特性,对泄漏的定位和定量是泄漏辨识研究的重要环节.考虑水力瞬变中非恒定摩阻的影响,将泄漏参数反映在建立的管道瞬变流时域数学模型中,给出离散网格和求解方法,分析不同泄漏参数对瞬变水击波的影响,用遗传算法结合瞬变流时域模型的特征线法最小化管道测点压力实测和计算响应值的均方误差来辨识泄漏参数,其中实测信号首先经过预滤波和小波去噪处理.通过模型试验对瞬变流数学模型和泄漏反问题分析模型作了率定和验证,结果表明此法能有效的辨识出小泄漏参数,并且检测信号主动,阀门扰动可操作性强.     

Modèle multi-bulles pour la cavitation

In this study we propose new multi-bubble model for cavitation, in which, to simulate the interactions within a cloud of cavitation at the initial stage, the dynamic behaviour of two nonidentical bubbles localised in a volume of control is studied. The presence of two bubbles introduces an instability in which the exchange of volume seems an additional degree of freedom. Depending on the conditions of expansion, the small bubble can disappear or not. If the small bubble disappears, the volume of control is readjusted to introduce a new small bubble and to continue calculation in a new sequence. The model makes it possible for many small bubbles to disappear as in the appearance of cavitation, which is at the origin of certain phenomena observed in the zone of the appearance, such as emission of the noise. The model reveals especially the pressure rather like a result than a datum.The comparison of the size of the bubbles and the pressure varying in time, obtained with the model are coherent with the measurements taken by Ohl [Phys. Fluids 14 (10) (2002) 3512–3521]. To cite this article: M. Adama Maiga, D. Buisine, C. R. Mecanique 337 (2009).     

MODELLING OF RELEASE OF GAS FROM HIGH-PRESSURE PIPELINES

The problem investigated is the break of a high-pressure pipeline carrying natural single-phase gas which may condensate (retrograde) when the pressure drops. Single-phase non-ideal gas is assumed using a general- ized equation of state. Taking advantage of the choked massflow condition, the break is split into a pipe flow problem and a dispersion flow problem, both solved using a finite difference control volume scheme. The transient flow field from the pipeline break location is expanded analytically, using an approximation of the governing equations, until ambient pressure is reached and matched to the corresponding gas dispersion flow field using as subgrid model a jet box with a time-varying equivalent nozzle area as an internal boundary of the dispersion domain. The turbulence models used for the pipe and dispersion flow fields are an empirical model of Reichard and the k–ϵ model for buoyant flow respectively. The pipe flow simulations indicate that the flow from the pipeline might include dispersed condensate which will affect quantitatively the mass flow rate from the pipeline and qualitatively the gas dispersion if the condensate rains out. The transient dispersion simulation shows that an entrainment flow field develops and mixes supersaturated gas with ambient warmer air to an unsaturated mixture. Because of the inertia of the ambient air, it takes time to develop the entrainment flow field. As a consequence of this and the decay of the mass flow with time, the lower flammability limit of the gas–air mixture reaches its most remote downstream position relatively early in the simulation (about 15 s) and withdraws closer to the break location.     

Water hammer with column separation: A historical review

Column separation refers to the breaking of liquid columns in fully filled pipelines. This may occur in a water-hammer event when the pressure in a pipeline drops to the vapor pressure at specific locations such as closed ends, high points or knees (changes in pipe slope). The liquid columns are separated by a vapor cavity that grows and diminishes according to the dynamics of the system. The collision of two liquid columns, or of one liquid column with a closed end, may cause a large and nearly instantaneous rise in pressure. This pressure rise travels through the entire pipeline and forms a severe load for hydraulic machinery, individual pipes and supporting structures. The situation is even worse: in one water-hammer event many repetitions of cavity formation and collapse may occur.This paper reviews water hammer with column separation from the discovery of the phenomenon in the late 19th century, the recognition of its danger in the 1930s, the development of numerical methods in the 1960s and 1970s, to the standard models used in commercial software packages in the late 20th century. A comprehensive survey of laboratory tests and field measurements is given. The review focuses on transient vaporous cavitation. Gaseous cavitation and steam condensation are beyond the scope of the paper.     

Approximate Solutions for Non-stationary Gas Permeability Tests

Full and approximate models describing changes in gas pressure for non-stationary single- and double-reservoir methods of testing of gas permeability of porous materials were formulated. The models include effects of compressibility and slip of gas. For derivation of the approximate models it was assumed that a transient mass flux of gas is constant along the tested sample (spatially homogeneous). The assumption allowed finding analytical solutions for evolution of reservoir pressures as well as the transient spatial distributions of pore pressure in tested samples. The full model was solved with help of the finite element method and the results were compared with the derived approximations. The comparison shows very good agreement of predictions of both models. Parametric studies have given the evaluation of the role of selected material parameters (permeability, porosity and slip factor) and parameters of the experimental system (initial pressure, volume of reservoirs), especially important from the point of view of identification procedures of the material parameters.     

超声波作用下气泡的非线性振动

Flynn方程在引入泡内气体的压缩性修正后能够更好地描述超声波作用下液体内气泡的非线性振动。我们以此为基础通过数值分析得到了不同初始半径的气泡在同一声场中的振动位移时间图像和相轨迹,发现气泡的运动行为对其初始状态有很强的依赖关系,在频率为26.5 kHz,幅度为1.35 atm的声波作用下,初始半径小于1 μm的气泡作小幅受迫振动,大于200 μm的气泡作小幅准本征振动,不具备空化气泡特征。声场中的小气泡对声波强度变化的反应更为激烈,当增加声强度时,可将一定范围内的小幅受迫振动气泡转化为空化气泡,并且,当驱动声波压力幅值增加时,初始半径越小的气泡的最大位移增加幅度越大。驱动声波频率同样影响气泡的振动。随着声波频率的升高,空化气泡的初始半径取值区间越来越小,空化振荡也越来越弱。本文还通过高速摄影系统对换能器作用于未除气的自来水所引起的变化进行了实验研究,结果表明,超声波作用下液体内的气泡场是一个混合场,场内气泡尺寸呈一定的分布状态,不仅有空化气泡,还有毫米级的大气泡。气泡的运动行为直接影响空化效果。超声空化场是一个复杂的物理场,场内除了有气泡的振动外,还有气泡间的相互吸引、碰撞和结合。