静压气体轴承气膜力及其与转子耦合动力学特性研究

针对静压气体轴承,进行三维实体建模,采用供气孔区域非结构化网格和非供气孔区域结构化网格相结合的网格划分方法;采用基于有限体积法的计算流体动力学(CFD)商业软件对三维稳态可压缩气体Navier-Stokes(N-S)方程进行求解;根据计算结果,通过数据拟合获得了考虑转子偏心和转速的静压气体轴承气膜支承力模型.基于有限元法建立了气体轴承-转子系统动力学模型,采用Newmark逐步积分法求解了系统的临界转速和不平衡响应.在此基础上进行实验测试,验证了数值仿真结果.研究结果表明:低速、小偏心下,气膜主支承力随偏心呈近似线性变化;高速、大偏心下,气膜主支承力急剧增大,气体轴承的动压效应显著增强;气膜x,y向耦合支承力随转速和偏心呈非线性变化;转子系统一、二阶临界转速对当前结构刚度变化比较敏感,而三阶临界转速对此不敏感.因此,气体轴承气膜支承力的非线性特性及其与转子耦合动力学特性较为复杂,在对气体轴承进行结构设计时,应充分考虑其与转子的耦合,合理设计工作转速范围.      

喷管内缝隙流场与热环境数值模拟研究

用有限元素法数值模拟了固体火箭发动机喷管内缝隙的流场和热环境。采用非结构的高分辨率的有限元离散格式,计算了不同攻角和不同Ｍａ数的缝隙内压力分布和气动加热,分析了缝隙内流场和热环境特性,数值结果与风洞实验结果符合很好。     

Langmuir双探针测量Hall推力器羽流特性

应用Langmu ir双探针测量了距离Hall推力器出口30~100 cm的羽流特性。实验中真空舱压力保持在5×10-4~2×10-3Pa之间。实验证明,电子温度和数密度随着角度和轴向距离的变化而变化,但羽流中远场的电子温度变化幅度很小,基本稳定在0.455 eV左右;电子数密度随着与轴向距离和中心轴夹角的增加而减小。实验结果可为高性能的Hall推力器设计提供参考依据,并为等离子体羽流的数值模拟提供输入条件和验证依据。     

弹塑性梯度材料热机耦合分析的有限元方法

陶瓷/金属功能梯度材料(FGMs)的组份及性能梯度分布,在航空航天高温热结构中有重要应用。材料性能在空间上的连续变化给这类结构强度设计的有限元分析造成了一定困难,针对上述问题,提出一种热弹塑性梯度有限元方法。该方法通过编写有限元软件ABAQUS的材料用户子程序(UMAT)来实现,在UMAT中将FGMs的弹塑性本构模型参数定义为空间坐标的函数;用热弹塑性梯度有限元方法对陶瓷-FGMs-金属典型的三明治板结构进行热机械分析。结果表明:FGMs组分布函数对其热应力分布有很大的影响,特别是在三明治结构的界面上,影响更为突出,这和现有文献的实验结果相符。     

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

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

搅拌头几何特征对搅拌摩擦焊试板温度场的影响

通过实验测量了不同类型搅拌头条件下铝合金焊接试板特征点的温度曲线,分析了搅拌头轴肩尺寸和搅拌针形状对焊接试板温度场的影响。结果表明,轴肩是搅拌摩擦焊接热输入的主要来源,搅拌头轴肩尺寸越大,产生的焊接热量越大,对应试板测温点的温度越高。搅拌头针形对焊接初始阶段试板的温度有明显影响,稳定焊接阶段,试板温度与搅拌针的几何特征及其作用下焊缝金属的塑性流动有关。采用圆台形搅拌针时焊接试板温度最高,采用圆柱形搅拌针时试板温度次之,采用螺纹形搅拌针时试板温度最低。通过对采用圆柱形搅拌针时垂直于焊缝方向上焊接试板温度数据的回归分析,得到了焊接试板宽度方向温度分布的二次解析式。     

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

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

面向对地成像观测任务的高空飞艇应急调度

针对应急条件下高空飞艇(HAA)对地成像观测任务调度问题进行研究,分析了问题中的主要约束条件,建立了以任务收益(TB)和巡航距离为优化目标的约束满足问题(CSP)模型。考虑飞艇侦察载荷具有侧摆观测能力,在构建视场范围约束模型和分辨率约束模型的基础上,对成像观测任务进行合成。提出了元任务与合成任务的概念,给出了任务合成的步骤与方法。将HAA应急调度问题转换为车辆路径问题(VRP),并进一步分解为任务排序主问题和路径选择子问题,分别应用改进粒子群(IPSO)算法和关键节点搜索(KNS)算法求解。详细介绍了算法中的编码、解码和移动等操作,以及采用的混沌初始化和禁忌搜索(TS)策略。通过仿真实验,对文中所提方法的有效性进行了验证。     

基于V型槽的空气/煤油燃气发生器火焰稳定方法研究

为实现空气/煤油燃气发生器的可靠点火和稳定燃烧,采用试验与数值仿真相结合的方法对V槽火焰稳定器应用于燃气发生器时阻塞比和安装位置的选取进行了研究。发现不同阻塞比的火焰稳定器安装于燃气发生器不同位置时局部油气比可以是富油或贫油,甚至出现火焰稳定器一侧为贫油,另一侧为富油的状态。从接近当量比侧向贫油侧的热量和质量传递可以使得贫油侧的火焰稳定性能提高。钝体稳定的空气/煤油部分预混火焰存在“熄火-再点火”过程,其熄火频率为30Hz,高于贫油预混火焰的10Hz。阻塞比约0.6的V槽火焰稳定性能较好,获得了空气/煤油燃气发生器采用钝体稳定火焰时稳定器安装的特征长度范围。      

CFD计算结果的静态及动态可视化显示方法

针对某翼型设计的 CFD计算结果,采用色调连续变化的图像映射方法,基于划分的计算网格,提出了三角单元结构的抽取和将 Gouraud-Shading与质心坐标插值法相结合的扫描线转换算法,并根据 HSV颜色模型的特点建立了可视化标量与颜色值之间的两种线性对应关系,在此基础上实现了翼型设计的 CFD计算结果中流体流函数和旋涡强度二维标量场的静态和动态可视化显示,实例证明了算法的有效性。