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

基于光滑粒子流体动力学方法（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°非线性加剧。     

中应变率下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推进剂的压缩变形,可为低温中应变率下固体火箭发动机药柱的结构完整性分析提供理论基础。     

低温高应变率条件下HTPB推进剂拉伸力学性能研究

为研究固体推进剂在低温高应变率条件下的力学性能,通过单轴拉伸实验和扫描电镜(SEM)断面观察,分别获取了HTPB推进剂在温度范围为-40~25℃及0.4~14.29s-1应变率下的应力-应变曲线和拉伸断面形貌。结果表明,HTPB推进剂的力学性能具有明显的温度和应变率效应。随温度降低和应变率升高,应力-应变曲线特性变得更加复杂,断面形貌基本上呈现"脱湿"越困难、颗粒断裂越明显的规律,低温和高应变率的"耦合"作用使得推进剂的损伤变得更加严重。初始弹性模量E和最大拉伸应力σm随温度的降低和应变率的升高而逐渐增加,且均与应变率具有相对较好的线性对数关系。低温和高应变率的"耦合"作用,使得-40℃及14.29s-1条件下的初始弹性模量和最大拉伸应力分别为25℃及0.4s-1条件下数值的1.6倍和3.2倍。与模量和强度相比,应变的规律性较复杂,其值随温度的升高而增加,且在常温下随应变率的升高而增加,但在低温下随应变率的升高而降低。通过双因素方差分析表明,低温高应变率条件下,温度和应变率均对最大拉伸应力有更显著的影响,同时温度较应变率对最大拉伸应力影响更加明显,而对模量则较弱。基于时温等效原理,得到了低温高应变率条件下HTPB推进剂的拉伸力学性能主曲线,该主曲线较大地拓宽了对推进剂力学性能的预测范围。     

固体火箭发动机喉衬用轴编C/C复合材料的烧蚀及热结构特性研究进展

固体火箭发动机喉衬用轴编C/C复合材料的工作环境面临高温、高压、高速燃气流和大量凝聚相颗粒的烧蚀和冲刷,对材料的抗烧蚀性能和热结构特性要求十分严格。因此,从烧蚀实验和热结构特性实验研究、热结构特性预测与气体-颗粒两相流数值模拟三个方面,论述了轴编C/C复合材料的烧蚀及热结构特性研究进展。总结讨论了实现真实烧蚀工作环境的模拟和影响烧蚀实验参数的控制是高温烧蚀实验的难点,对于铝颗粒添加下工况的烧蚀实验和在细观尺度下热结构特性参数的测定实验是重点;提出从实验件类型、实验系统设计和对比有无铝颗粒添加下的工况进行烧蚀实验;提出采用一种热稳定性材料取代界面的实验方案进行热结构特性参数的测定实验。在热结构特性研究的细观尺度方面,组分材料之间的界面对热结构特性的影响有待深入研究,提出在代表性体积单元模型的基础上引入温度的周期性边界条件来实现热结构参数的预测。在气粒两相流数值模拟方面,发动机内不同相之间相互耦合作用以及对轴编C/C复合材料的机械侵蚀是数值模拟研究的难点,提出使用SDPH-FVM耦合的方法去解决内流场燃烧流动的问题,进一步可用来揭示内流场燃烧流动机理。     

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

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

基于SIFT算法的含碳颗粒凝胶推进剂雾化场速度分析

雾化是凝胶推进技术的关键问题之一,雾化过程中的液膜、液丝及液滴等的速度分布、液膜厚度等参数对于判断雾化效果、揭示雾化机理具有重要作用。为了对含碳颗粒凝胶推进剂雾化场的速度进行定量分析,提出了一种基于SIFT关键点匹配的雾化场速度计算方法,并以雾化场速度分析为基础,提出了一种新的液膜厚度估算方法。结果表明:雾化场速度随射流速度的增大而增大,随撞击角度的增大而减小;雾化场平均速度与射流速度的比值vato/vjet在0.6~0.9,可用于表征雾化效果,vato/vjet越小,雾化效果越好;射流撞击形成液膜的厚度在0.04~0.13mm,液膜厚度随着射流速度及撞击角度的增大而减小;凝胶推进剂的雾化效果随着碳颗粒浓度的增大而降低,随着碳颗粒粒径的增大而略有改善。     

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

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

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

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

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

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