一种简化算法的航空发动机全状态数学模型

通过一种简化的燃气空间模型对某型双转子混排加力涡扇发动机各部件进行了建模和仿真计算,并给出了详细的计算流程以及数学公式.在仿真过程中结合实际工程应用经验,对涡轮处采用一个局部的小循环迭代以减小涡轮部件内燃气空间计算的时间步长,提高计算精度.与其他数值方法相比,该简化算法避免了求解非线性方程组时大范围的循环迭代,提高了全状态数学模型的计算速度.通过仿真结果与发动机实际试车数据的对比,证明该方法有足够的计算精度.      

运用SIMPLE方法计算大曲率弯道中的湍流流动

在三维任意贴体非正交曲线坐标系下应用同位网格技术,发展了适于求解复杂流动的压力修正算法,采用速度的协变物理分量为动量方程中速度的求解变量,湍流模型选用标准的κ-ε模型.计算表明,方法具有良好的收敛性、准确性和稳定性,计算值与试验数值吻合良好,而且能预示通道内复杂的流动现象.     

三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

预拉伸条件下铝合金焊接残余应力的数值模拟

利用ANSYS有限元分析软件,通过APDL语言编程,实现了不同预拉伸应力条件下铝合金平板对接焊过程中,温度场及应力场的模拟分析.计算结果表明,预拉伸焊接法可以有效地控制焊接残余应力,随着预拉伸应力的增大,其焊后残余应力值逐渐减小.当预拉伸应力σP从0增加到90%σ0.2时,纵向残余应力降低了85.6%.模拟分析结果与实验测试结果基本吻合.     

改进的电热除冰系统仿真

建立了电热除冰系统矩形微段的二维简化物理模型,采用焓法模型建立相变导热问题的微分控制方程,用控制容积法对该模型的控制方程及其边界、初始条件进行离散,使用块修正技术对原有的离散方程组的求解方法进行了改进,得到各个时刻的温度分布及冰的融解情况.计算结果与文献实验结果吻合,验证了改进的算法和程序的正确性.并计算分析了时间步长的选择、网格的划分对计算结果的影响.      

微型涡喷发动机燃烧室设计研究

设计了微型环型燃烧室,并对其燃烧过程进行数值模拟,同时根据模拟结果对燃烧室进行了优化设计及二次计算.计算中入口边界参数采用等熵数值计算结果,喷雾性能参数通过试验获得,考虑了湍流化学反应、热辐射、液滴的蒸发及二次雾化等情况,计算结果可以较为准确地反映燃烧室的实际燃烧情况.对优化设计前后的模拟结果进行比较分析,结果表明优化后的燃烧室设计合理,燃料可以充分燃烧,各项性能参数基本符合设计要求,能够为燃烧室的进一步优化设计提供改进依据.      

一种TVD方法在三维非结构网格中的应用

以三维非结构网格的显式有限体积法为基础,采用了一种TVD方法求解三维Euler方程,使用Roe通量来计算网格单元边界处的守恒量通量.为了验证方法的可行性,用该方法模拟三维爆炸问题,得出的结论是我们的方法可行,稳定且有效.     

机械连接件中螺栓孔角裂纹应力强度因子数值计算(英文)

The three-dimensional finite element method is used to solve the problem of the quarter-elliptical comer crack of the bolt-hole in mechanical joints being subjected to remote tension. The square-root stress singularity around the corner crack front is simulated using the collapsed 20-node quarter point singular elements. The contact interaction between the bolt and the hole boundary is considered in the finite element analysis. The stress intensity factors （SIFs） along the crack front are evaluated by using the displacement correlation technique. The effects of the amount of clearance between the hole and the bolt on the SIFs are investigated. The numerical results indicate that the SIF for mode I decrease with the decreases in clearance, and in the cases of clearance being present, the corner crack is in a mix-mode, even if mode I loading is dominant.     

炮弹飞行速度等对核电站建筑的毁伤研究

对核电站设施进行了分析,研究了炮弹对核电站安全产生最大威胁的情况,并对此过程进行了数值模拟;应用有限元建模技术,实现了对核电站安全壳的有限元建模.运用ANSYS-DYNA有限元软件平台进行了数值仿真,得出炮弹对核电站的侵彻效果.通过仿真验证,较准确地评估了炮弹对核电站构成的威胁,可供核电站选址和防护设计时参考.     

高超声速圆锥边界层转捩数值研究

以类似于处理湍流脉动的方式来模型化非湍流脉动,考虑湍流的间歇性特征及引起超声速流动失稳机制,研究包括转捩起始位置及转捩区范围的高超声速流动转捩问题,发展可反映扰动模态和可压缩性影响的两方程SST转捩/湍流模式.针对尖锥模型,在马赫数6、攻角O°～2°时,采用第二模态(Second-mode)机制,预测并比较完全层流态、完全湍流态,以及自然转捩发展过程中的表面热流、表面摩阻等气动特性.