基于Fibonacci搜索方法的航空发动机性能寻优

介绍了Fibonacci搜索方法应用于发动机性能优化中的原理,建立带两个出口面积、气流马赫数可变的引射器的涡喷发动机的总体性能计算模型,应用Fibonacci搜索方法在优化参数区间内对该发动机进行了最大推力、最低油耗、最低涡轮前温度模式的寻优.举例说明了Fibonacci搜索方法应用的具体过程.给出了不同高度、马赫数下的寻优结果.结果表明:最优目标性能满足的约束边界与飞行高度、马赫数密切相关.不同的飞行高度、马赫数对应了不同的最优控制策略.但不论飞行任务如何,根据该方法进行的优化控制所获得的目标性能是最佳的.该优化方法不仅适用于该类涡喷发动机,也为其他类型航空发动机的优化控制打下一定的基础.      

引入压力梯度的针对角区分离流动的Spalart-Allmaras模型的改进

针对Spalart-Allmaras（S-A）模型在角区分离计算中的问题,将无量纲的压力梯度引入其涡黏性输运方程的生成项,得到了改进的S-A模型.通过对两套含角区分离的低速压气机叶栅进行验证计算发现:与实验结果相比,原始S-A模型所得的分离区偏大,分离区内壁面压力偏低;而改进模型得到了与实验一致的分离区尺寸以及吸力面、压力面压力系数分布等结果.针对S-A模型涡黏性生成项和耗散项的分析表明:引入的无量纲压力梯度有效的识别了角区分离,在分离区内改变了涡黏性的生成、耗散关系,增大了涡黏性,从而缩小了计算所得分离区,同时在主流区保留了原始S-A模型的计算结果,进而带来了良好的改进效果.      

300M钢激光熔覆耐磨防腐自润滑涂层温度场数值模拟研究

为了提高飞机起落架减震支柱300M超高强钢的抗磨性能,突破由于温度梯度过大诱发的激光熔覆耐磨防腐自润滑涂层裂纹等技术瓶颈,基于ANSYS软件的生死单元法编制热循环程序,考虑自润滑相和耐磨相材料热物性参数随温度的变化、相变潜热、激光熔覆过程中与外界的换热、激光熔覆功率、激光熔覆扫描速率等因素对激光熔覆过程温度场、熔池、温度梯度的影响,建立300M超高强钢激光熔覆耐磨防腐自润滑涂层温度仿真模型。结果表明:基体的熔化需要激光、熔化的粉末等综合作用使传导到该区域的有效能量达到熔化的临界值,熔化高度增加率随激光功率的增加先降低后增加,熔化高度减少率随激光扫描速率的增加先变小后变大;由于激光熔覆的不同区域温度、冷却速率差异等综合因素影响,熔覆层熔池的纵截面为勺状熔池;伴随激光功率的增加,由于熔覆层不同区域对能量输入产生的温度响应速率差异导致Z方向温度梯度值增加,且最大冷却速率增加;伴随激光扫描速率的增加,激光输入能量降低,降低了高温区域温度及激光的快速局部加热等综合影响,Z方向温度梯度值降低。通过调控激光参数,在保持熔覆层结合强度的条件下,基体熔凝区能被控制到最小,并能够降低温度梯度。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

TGO非均匀增长对热障涂层应力演化和破坏机理的影响

采用内聚力模型和热生长氧化层(TGO)非均匀增长子程序,数值模拟了在热循环载荷作用下热障涂层(TBC)内部应力演化规律和开裂行为。涂层失效过程首先是源自陶瓷层(TC)内近波峰位置的拉伸和切应力共同主导的陶瓷层Ⅰ、陶瓷层Ⅱ混合型裂纹;随着循环数增加,则转向由TC内近波峰位置的切应力主导的Ⅱ型裂纹和波峰波谷中间的涂层厚度方向拉伸应力主导的Ⅰ型裂纹。整体非均匀增长和波谷均匀增长模式下的最大拉伸应力经过一定循环数后几乎不再随循环数而增加;而在波峰均匀增长和整体均匀增长模式下,最大拉伸应力则会随着循环数增加持续增长。整体非均匀增长、波谷非均匀增长模式下,20个循环后最大切应力出现在近波峰位置,分别为-16241 MPa和-15428 MPa;而整体波峰均匀增长和整体均匀增长模式下,最大切应力为-11382 MPa和-11198 MPa。对于波谷均匀增长和整体非均匀增长模式,在9个循环后出现界面裂纹。而对于波峰均匀增长和整体非均匀增长模式,在第17个循环出现界面裂纹。     

激光冲击强化对钛合金高温微动疲劳寿命的影响

为了解高温工作环境下激光冲击强化工艺(LSP)对钛合金材料微动疲劳寿命的影响,开展了强化前后TC11钛合金在室温、300°C和500°C下的微动疲劳试验并测试了试验件表层的残余应力及硬度。结果表明:随着温度的升高,激光冲击强化对TC11钛合金微动疲劳寿命的提高倍数逐渐减小。在轴向载荷为400MPa,法向载荷为65.5MPa时,经激光冲击强化后TC11钛合金试验件在室温、300°C和500°C下的微动疲劳寿命分别为强化前的5.5倍、3.5倍和1.7倍;强化后试验件表层的残余应力会在高温下发生松弛,且松弛程度会随温度的升高而增大,这是激光冲击强化效果随温度升高而逐渐弱化的主要原因。     

AZ31B镁合金双道次热变形过程流变应力特征模拟研究

建立了一类镁合金双道次热变形过程的二维元胞自动机模型。考虑到镁合金具有六重对称的密排六方晶体结构,模型采用了六边形元胞网格。另外,通过计算位错密度演变和设定晶粒形核长大规则,模型将热变形过程中动态再结晶、静态回复、静态再结晶、亚动态再结晶和晶粒长大等单个物理过程耦合到一起。为了验证模型,将其应用到AZ31B镁合金单道次和双道次热压缩过程。探讨了不同变形温度、应变速率、道次间隔和预应变对应力-应变曲线的影响。将计算结果和实验结果进行了对比,吻合较好。     

Mesh-dependence of Material with Softening Behavior

The issue of mesh-dependence emerges when the conventional continuum damage model is applied to handling the softening behavior. In order to circumvent the mesh-dependence, the non-local theory is introduced into the conventional damage model and the finite element formulas are derived for two-dimensional gradient-enhanced damage model. A new element is proposed in which the basic unknown quantities are displacement, non-local equivalent strain and the gradient of non-local equivalent strain. The element and constitutive equation proposed in this article are added to the finite element software ABAQUS through user subroutine UEL. Numerical results show that the gradient-enhanced damage model can eliminate the mesh-dependence and is effective for dealing with the issue of softening behavior.     

新电位移边界条件下的压电介质断裂问题

应用边界元法，并结合一种特殊的复变函数基本解，研究了压电介质的断裂问题。由于该基本解完全遵循本构方程没有对电位移边界作假设，所以结论同非穿透性裂纹假设的结果有较大差异。通过计算表明，外加力载荷不但会引起裂纹尖端应力奇异性，也会引起电场奇异性，同样，外加电压也会引起裂纹尖端应力奇异性。     

终端区点融合进近程序设计方法研究

流量的增长导致终端区拥挤和航班延误,增大了管制压力和公司成本.以长沙黄花国际机场36号跑道为例,空城受限,4条航线汇聚一点,梯度曲折多变,航班排序难,指挥复杂度高.结合我国飞行程序设计和管制间隔要求,根据点融合的基本概念,设计跑道点融合进近程序,确定相关参数.对程序进行容量和下降梯度评估表明,静态容量可达到63架次/h,是当前运行最大流量的3.7倍;下降梯度得到明显改善,为连续下降实施奠定了基础.