Fatigue-crack-initiation numerical modelling of a Ni-base powder metallurgy alloy

A simplified three-dimensional numerical model was presented to simulate the micro-crack nucleation and growth to some predefined dimension (approximately 0.38mm) on the throat surface of a Ni-base powder metallurgy (PM) specimen.The numerical simulation of micro-crack initiation was based on the Tanaka-Mura micro-crack initiation models,where individual grains of the mesoscopic model were simulated using the Voronoi tessellation.Four improvements were made in the model.(1) Considering crack initiation along with 12 principal slip systems on octahedral slip planes of face centered cubic (FCC) crystal in three-dimensional (3-D) models.(2) Considering the effect of secondary phase precipitate due to impinging slip and dislocation pileup.(3) The Tanaka-Mura theory of fatigue-crack-initiation from notches was applied to simulate the crack initiation from another crack tip.(4) The coalescence of random initiated micro-cracks was simulated once they intersected with each other and a macro-crack was finally formed.The calculated results were in good agreement with the experimental data which verified the rationality of the simulation model.The applicability of the proposed model for treating fatigue-crack-initiation life in engineering structures was preliminarily achieved.      

基于GasTurb/MATLAB的航空发动机部件级模型研究

基于航空发动机总体性能分析软件GasTurb及其源代码,利用动态链接库技术实现在控制系统开发平台MATLAB下直接调用GasTurb部件级动态模型,并在Simulink下建立了包含涡喷、涡扇、涡轴、涡桨在内的22种发动机类型的部件级模型库,实现了两者之间的无缝衔接.应用实例与仿真结果验证了所建模型的有效性与高度可定制性.消除了从总体性能分析阶段过渡到控制系统开发阶段存在的交互障碍,为发动机控制和故障诊断研究提供了一种灵活的仿真平台.      

飞机推出和跑道掉头滑行行为计算方法

针对飞机推出和跑道掉头两种地面滑行特殊行为的位置和姿态计算问题进行了研究。从机场塔台管制运行模拟仿真的角度出发,提出了非完整运动学约束下飞机推出滑行和跑道掉头运动计算方法。依据飞机推出实际和非完整运动路径跟踪模型,设计出适合仿真推演计算的飞机推出运动求解方法。根据飞机跑道掉头的几何约束和转弯角运动约束,设计出开环控制下完整刻画飞机运动的求解方法。与动力学求解模型相比,该方法具有计算简便,适合于视景仿真驱动的优点,并通过数学仿真验证了方法的有效性。     

温度冲击下变截面圆孔装药的热力耦合分析

为研究温度冲击载荷对某改性双基推进剂装药结构完整性的影响,基于热力耦合基本理论,开展了不同温度下改性双基推进剂力学松弛试验,得到了不同温度下松弛模量的6阶Prony级数及以293.15 K为参考温度的时温等效W.L.F(Williams-Landel-Ferry)方程。利用Ansys有限元软件对由333.15 K降低至218.15 K的温度冲击条件下装药结构完整性进行分析,结果表明:装药的温度和应变在温度冲击前期变化较剧烈,至12 000 s时基本达到平衡,装药前端面两个拐角处在温度冲击前期发生破坏;为验证仿真计算结果,开展了由333.15 K降低至218.15 K的温度冲击试验。在温度冲击过程中,药柱沿轴向和径向出现裂纹,起裂点为推进剂装药前端两个拐角。温度冲击试验结果验证了仿真计算的准确性,说明该计算过程可用于装药结构完整性分析。      

基于大涡模拟的离心式喷嘴雾化过程研究

为实现离心式喷嘴雾化过程的精确数值仿真,探究喷嘴内部流动特性与外部液膜破碎形式,采用基于大涡模拟的仿真方法,对一种典型的四进口离心式喷嘴进行研究,仿真结果揭示了喷嘴内部相界面的振荡现象与外部液膜的破碎细节,并通过耦合流体体积法（VOF）与离散相模型（DPM）,获得液滴粒径的空间分布特征。研究结果表明：在液体填充过程中,喷嘴内的气液相界面存在波动与褶皱,形状并不稳定,内部的空气芯直径呈现正弦模式的振荡变化,喷嘴出口液膜厚度沿周向分布不均,这些因素导致出口附近的液膜表面出现扰动。在不同的进口条件下,不稳定性导致液膜表面上的扰动波形式不同。进口压力为0.3MPa时,液膜破碎由开尔文-亥姆霍兹（K-H）不稳定性产生的轴向正弦波所导致,产生沿周向分布的环形液带;在0.7MPa下,液膜表面开始出现由瑞利-泰勒（R-T）不稳定性引发的周向扰动波;随着压力增加至1.1MPa,液膜的破碎则由R-T不稳定性主导,产生沿轴向分布的液带结构,随后在气动力与表面张力的作用下破碎成液滴。二次雾化破碎后,喷嘴外部截面内的粒径呈“单谷”分布,液滴平均粒径计算结果与实验的最大相对误差为5.1%,与实验数据吻合度较高。     

正十烷简化机理的构建及其在发动机燃烧数值模拟中的应用

发动机燃烧室中燃料的能量释放与燃烧特性对于发动机设计具有重要作用,为了预测发动机点火包线和贫/富油极限等关键性能,迫切需要发展航空燃料及其典型组分的高精度化学动力学模型。本文针对燃料典型组分正十烷,采用自主开发的机理生成程序ReaxGen构建了其燃烧详细机理（1499种组分、5713步反应）。为了验证机理的合理性与可靠性,在当量比Φ=0.5-2.0,压力P=1-80 atm的宽工况条件下进行了点火延迟模拟验证,结果表明本文提出的正十烷详细机理在较宽的温度、压力和当量比条件下具有较高的模拟精度。为获得适用于发动机燃烧模拟的高精度简化机理,本文基于误差传播的直接关系图方法简化了正十烷燃烧详细机理,得到包含709种组分、2793步反应的正十烷半详细机理。进一步在高温范围（1000-1500 K）,采用路径通量分析方法简化得到含77种组分、359个反应的骨架机理。获得的骨架机理能够合理描述正十烷在高温下的燃烧特性,且该骨架机理尺度规模可用于基于火焰面模型的燃烧数值模拟。基于此高精度的骨架机理模型,结合火焰面生成流形湍流燃烧模型,采用大涡模拟方法进行了航空发动机环形燃烧室单头部扇形的燃烧模拟,初步获得了非稳态流场结构,其中温度模拟结果与实验值基本符合。     

基于流固热耦合的非金属迷宫密封泄漏特性与公式构造

建立了考虑齿变形的非金属迷宫密封泄漏特性流固热耦合数值求解模型,在验证求解模型准确性的基础上,研究了聚醚醚酮（PEEK）、填充30%碳纤维增强聚醚醚酮（PEEK-CA30）和铝合金三种材料迷宫密封在不同压比、温度下的流场特性、结构力学特性与泄漏特性,并基于Vermes公式构造了考虑齿变形的非金属迷宫密封泄漏量理论公式。研究结果表明:建立的迷宫密封流固热耦合模型可以准确计算非金属密封齿的变形量和变形后的泄漏量。在研究的三种材料中,采用PEEK-CA30材料密封齿的变形量相对较小,占密封齿径向长度的0.35%~0.67%,其密封性能较好,相比于未考虑齿变形密封的泄漏量增加1.2%~6.8%。当温度高于500 K,压比大于3时,采用铝合金材料密封齿的最大等效应力达到材料的屈服极限而引起密封件失效。所构造的泄漏量理论公式能够准确预测考虑齿变形的非金属迷宫密封泄漏量,为非金属迷宫密封泄漏特性分析提供理论依据。      

刷式密封流动传热特性数值方法

开展了刷式密封流动传热特性数值方法研究,分别建立了刷式密封多孔介质、稳态实体与瞬态流固热耦合求解模型,设计搭建了刷式密封泄漏流动特性实验装置,在实验验证数值方法准确性基础上,对比分析了3种数值方法的差异性,研究了刷式密封流动传热特性,揭示了刷式密封的封严与传热机理。研究结果表明:在研究工况下,刷式密封多孔介质、稳态实体、瞬态流固热耦合模型泄漏量计算值与实验值的对比误差分别为9.8%~17.1%、8.1%~10%、6.92%~9.01%。刷式密封多孔介质模型计算速度较快,但需通过实验修正孔隙率,湍流模型对稳态实体模型流动传热特性结果影响较大,瞬态流固热耦合模型考虑了流场、刷丝及摩擦热三者间相互耦合作用,计算精度较高,但所需计算时间较长;同一压比下刷丝束温度从上游至下游逐渐增加,刷丝束最高温度随压比的增加而增大。气流流经刷丝间隙形成的节流效应致使泄漏气流能量耗散是刷式密封封严的主要原因,泄漏气流与刷丝表面间的对流换热是刷式密封摩擦热耗散的主要形式。      

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

低压天然气喷射的船用二冲程低速双燃料发动机爆震研究

为了解决大缸径船用预混天然气双燃料发动机的爆震问题,并拓展天然气的稀燃边界。基于三维数值模拟的方法对大缸径船机进行了仿真模拟。分析了大缸径双燃料发动机爆震的特点,并对缸内涡流强度和废气再循环（EGR,Exhaust Gas Re-circulation）率对爆震的影响进行了研究。研究结果表明：大缸径预混天然气发动机的爆震位置往往发生在气缸边缘,火焰面的传播过程是引起缸内爆震的主要因素。随着缸内涡流从无到有的增强,缸内的爆震强度随之增强;当涡流到达一定程度后,随着涡流的增强,缸内的爆震强度反而降低;缸内加入EGR可以提高天然气当量比的同时减少爆震的强度,可以拓展天然气的稀燃边界。