基于光滑粒子流体动力学方法的液滴蒸发问题数值模拟研究

基于光滑粒子流体动力学方法（Smoothed Particle Hydrodynamics，SPH），开展了SPH新算法在蒸发燃烧领域的研究。建立了适用于SPH方法的蒸发数值模型，推导了基于傅立叶热传导公式和菲克扩散定律的SPH离散方程；借鉴VOF方法（Volume of Fluid）的思想，提出了SPH粒子的液相质量分数的概念，以有效表征蒸发过程中的相变问题。采用SPH方法对高温环境中单个液滴的蒸发过程进行数值模拟，结果符合D2定律，与理论模型相一致；在强迫对流环境中，液滴的蒸发过程受到对流作用及表面张力的影响，蒸发速率加快；进一步对双液滴在静止、对流环境中的蒸发过程进行数值模拟研究。结果表明，液滴的间距、滴径对多个液滴的蒸发过程影响至关重要，液滴间距至少在两倍的液滴直径以上，相互之间的影响才可以近似忽略。通过本文研究，拓宽了SPH方法在蒸发相变领域的应用范围，研究结果也能够为进一步的燃烧问题研究奠定基础。     

面向一体化计算的整体式液体冲压发动机建模及性能评估

为发展冲压发动机性能工程预估方法,建立了面向一体化计算的整体式液体冲压发动机性能计算模型,并使其兼容基团贡献算法,提高了拓展性。利用模型分析了冲压发动机在攻角0°～6°、高度0～18km、马赫数2.0～3.5、余气系数1.0～2.9范围工况多维度连续变化下,比冲、推力系数和燃油质量流量的速度-高度特性、高度-节流特性、节流-速度特性以及攻角特性。研究结果表明:性能计算模型可在一体化计算条件下、基团贡献算法允许范围内,不依赖试验数据对冲压发动机性能预估,计算结果与技术参考值相比相对误差均小于14%。推力系数和比冲具有基本一致的速度-高度特性和相似的高度-节流特性,受燃气组分影响,推力系数和比冲的节流-速度特性差异明显。燃油质量流量的变化规律不同于比冲和推力系数,在进入平流层后呈现折缓趋势,而呈现连续性,速度越大、高度越低、余气系数越小,燃油质量流量越高,反之则越小。引入攻角以后,攻角越大,比冲和推力系数越低,进气道起动马赫数越高;攻角-起动马赫数曲线小范围内近似线性,攻角超过5.6°非线性加剧。     

运输环境中某战术导弹结构振动特性分析

应用三维有限单元和导弹实体简化模型对一种战术导弹的结构振动特性进行分析，并动用ＡＮＳＹＳ有限元分析程序，计算得到了“自由－固支”条件下，导弹竖立状态的固有模态，取得了与试验和数据基本一致的结果，证明了所建模型物正确性，进而对运输环境中导弹水平支承条件下的固有频率和振型进行了分析计算，所得结果为导弹的使用可靠性分析提供了依据。     

温度参数下轴编C/C复合材料的超高温热结构性能分析

为了研究轴编碳/碳（C/C）复合材料的超高温热结构性能,开展了2800℃超高温载荷下复合材料的拉伸实验,表征了复合材料组分材料的高温失效形貌;提出了立方体单胞和正六棱柱单胞的温度周期性边界条件的算法和预测复合材料等效热结构参数的方法,在材料属性各向同性正六棱柱单胞模型上验证了方法的正确性;最后建立了轴编C/C复合材料的代表性体积单元,计算了复合材料随温度载荷变化下的等效热结构参数,与实验数据进行了对比分析。分析表明：开展的轴编C/C复合材料的超高温拉伸实验,可以得到复合材料在2800℃载荷下复合材料的等效刚度值,同时分析了轴向试件和径向试件的高温失效机理,径向试件的纤维束拉伸强度对复合材料的径向拉伸强度起到决定性的作用。在应用位移以及温度周期性边界条件的基础上,提出的方法可以得到随温度载荷变化下复合材料的等效热结构参数,得出复合材料在横向满足近似各向同性的特性,预测结果表明只改变温度参数下复合材料的等效热导率不会发生改变。提出的方法同样也适用于其他编织复合材料的热结构性能研究。     

二氧化氮催化硝化甘油分解机理的数值计算

采用密度泛函理论（DFT）在B3LYP/6-31+G（d,p）水平上模拟计算了NO₂对硝化甘油（NG）催化分解的反应机理。通过与O—NO₂均裂、HNO₂分子内消去反应对比发现,NO₂分子对NG分解反应有显著的催化效应;NG的催化分解途径有α-H夺氢催化和β-H夺氢催化,计算得出2种途径的反应能垒分别为105.337 k J/mol和124.381 k J/mol,说明α-H夺氢催化更容易发生。     

凝胶推进剂锥形管道流动特性数值分析

对凝胶推进剂在锥形圆管中的流动过程进行了数值计算,研究了锥形圆管收敛角对轴向速度、平均表观粘性、压降的影响关系,得到了凝胶推进荆在锥形圆管中的轴向速度与平均表观粘性分布.计算结果表明,随着收敛角的增大,出、入口截面平均表观粘性降低幅度不断增大,出口截面平均表观粘性不断减小至近牛顿粘性水平ηoo;并且,粘性的减小是以增大压降需求为前提的,当角度改变达某定值之后,角度的改变引起粘性的变化将不再显著.结果表明,圆管收敛角是影响粘性变化的一个重要参数,粘性变化与压降需求之间存在最佳结合点.     

残余应变对锥柱形膜片多次翻转行为的影响研究

弧长法能跟踪复杂的非线性平衡路径全过程,但计算受初始设定载荷影响较大,而针对这一局限设计的改进算法能更好地求解膜片翻转过程。考虑新型Ti3012合金、纯钛、理想材料,基于改进弧长法研究残余应变对锥柱形膜片多次翻转行为的影响,得出新型Ti3012合金膜片的正翻与回翻过程均没有纯钛膜片好,且存在严重偏心现象,其98%正向翻转状态时回翻过程中不仅偏心量达到60mm,而且出现褶皱等严重失效行为;理想材料在不同正向翻转状态时反向回翻,膜片均存在严重的偏心问题,最大偏心量大于80mm;纯钛膜片仅98%正向翻转状态时回翻过程中偏心量稍大,约为10mm。无残余应变这一材料特性降低了膜片结构刚度,同时伴有初始弹性模量减小的性质,这些因素均不利于膜片多次翻转。     

幂律型流体雾化SPH方法数值分析

幂律型流体的雾化过程存在复杂的界面运动,用传统网格法很难精确追踪运动界面.为研究幂律型流体的雾化特性,运用SPH方法对典型的双股幂律型流体撞击雾化问题进行三维数值分析.根据文献实验,采用SPH方法模拟获得与实验条件相应的数值结果,对比后表明,二者雾化角相当吻合,数值结果还成功捕捉到液膜向液丝的破碎过程及网状的液丝分布状态,验证了方法的有效性.分析了射流速度和撞击角度对雾化角的影响,得到雾化角随着射流速度和撞击角度的增加而增大;对雾化后的速度场进行数值分析后表明,撞击点附近,惯性力的作用使速度场变化剧烈;在流体远离撞击点的过程中,粘性耗散作用使速度场趋于稳定,但速度大小小于初始撞击速度.     

中应变率下HTPB推进剂压缩力学性能和本构模型研究

为研究固体推进剂在中应变率条件下的压缩力学性能,在高应变率液压伺服试验机上开展了单轴压缩实验,并获取了温度范围为-40~25℃及0.40~85.71s-1应变率下HTPB推进剂的应力-应变曲线。结果表明,本文的实验方法是有效的,温度和应变率对HTPB推进剂的压缩力学性能影响显著。随温度降低和应变率升高,应力-应变曲线特性变得更加复杂,并与准静态下的应力-应变曲线特性有明显区别。压缩模量E和压缩应力σ0.17随温度的降低和应变率的升高而逐渐增加,且均与应变率具有相对较好的线性双对数关系。在低温和较高应变率的双重作用下,-40℃,85.71s-1条件下的压缩模量E和压缩应力σ0.17分别为25℃,0.40s-1条件下数值的10.64倍和4.25倍。基于时温等效原理,得到了HTPB推进剂的压缩力学性能主曲线,该主曲线能够对低温较宽应变率范围内推进剂的压缩力学性能进行预测。在朱-王-唐非线性粘弹性本构模型的基础上,构建了考虑温度和应变率效应的固体推进剂中应变率压缩本构模型,并采用遗传算法拟合了本构参数。通过不同温度和应变率下预测结果与实验数据的比较,验证了模型的有效性。所建模型能够较好地描述0.17应变以内HTPB推进剂的压缩变形,可为低温中应变率下固体火箭发动机药柱的结构完整性分析提供理论基础。     

压电阻抗技术监测固体推进剂老化研究

为了进一步研究表面粘贴式的主动传感技术——压电阻抗技术(EMI, Electro-Mechanical Impedance)监测固体推进剂的老化,构建了基于线粘弹杆结构的一维机电耦合模型,根据粘弹性波的理论与边界条件进行分析,建立一种线粘弹结构动态模量与其机械阻抗的表征关系,并对构建的机电耦合阻抗模型进行了数值计算及试验验证。针对HTPB固体推进剂开展高温热加速老化试验及压电主动激励试验,根据监测所得导纳频谱并通过提取结构机械阻抗进行分析,结果表明:压电陶瓷片与推进剂结构在高频机电耦合运动时(300kHz,280kHz)有明显的共振现象,结构固有频率高于机电耦合共振频率,并且在200kHz～400kHz频段与700kHz～900kHz频段内,结构机械阻抗频谱峰值会因固体推进剂的热老化时间增长而降低,且与热老化时间之间满足线性关系。由此可见,压电阻抗法能够与固体推进剂在力学性能上建立表征关系,通过得到固体推进剂结构机械阻抗能够监测其老化损伤。