基于人机闭环稳定性的舵机速率限制边界

舵机速率限制是造成电传操纵飞机人机耦合的主要原因。利用描述函数法对舵机速率限制非线性进行建模,分析舵机速率限制非线性、人机耦合发生频率及人机闭环稳定性间的关系,并基于此提出舵机速率边界的确定方法;以典型放宽静稳定性飞机为例,基于最优McRuer驾驶员模型,确定人机闭环稳定性所需的最小舵机偏转速率;基于开环起始点(OLOP)准则对所确定的速率限制边界进行验证。结果表明:本文提出的舵机速率边界的确定方法最小成本地避免了人机耦合;所确定的舵机速率限制边界与OLOP准则边界对应的舵机速率基本吻合,即所建立的舵机速率限制边界确定方法合理。     

基于改进牛顿法的航空发动机起动建模技术

为进行涡扇发动机起动过程性能仿真研究,在发动机慢车以上部件级模型的基础上,通过指数外推法获得低转速部件特性,在建立起动附件模型的同时,修正起动过程各部件总压恢复系数,从而建立起包含地面起动过程的全状态性能模型。在求解发动机各部件共同工作非线性方程组时,采用基于传统牛顿拉夫森迭代求解方法的改进方法,即在第1次完成中心差分得到Jacobian矩阵后,使用修正项迭代更新上一步的Jacobian矩阵。结果表明:该建模方法可以有效地建立起动过程性能模型,并大幅改善模型的执行效率、提高模型仿真求解速度。     

固定几何气动矢量喷管气动性能数值仿真

为掌握固定几何气动矢量喷管气动性能,通过CFD数值模拟的方法,研究了主流落压比、扩张段二次流落压比、扩张段二次流角度和引射对固定几何气动矢量喷管轴向推力系数的影响;主流落压比、扩张段二次流落压比和扩张段二次流角度对矢量角的影响;主流落压比、喉道二次流落压比和喉道二次流角度对喉道控制率的影响。结果表明:随主流落压比增大轴向推力系数增大,矢量角减小,喉道控制率减小;随扩张段二次流落压比增大推力系数减小,矢量角增大;随喉道二次流落压比增大,喉道控制率增大;随扩张段二次流角增大轴向推力系数减小,矢量角略有减小;随喉道二次流角增大喉道控制率增大;随引射方式增加喷管推力系数增大。     

复合固体推进剂细观结构建模及脱黏过程数值模拟

为研究复合固体推进剂损伤演化规律,基于分子动力学颗粒填充算法构建了HTPB(hydroxyl terminated polybutadiene)推进剂细观结构模型,通过在AP(ammonium perchlorate)颗粒/HTPB基体界面处引入黏接接触替代传统的黏接单元,并基于Hooke Jeeves的参数优化算法反演得到颗粒/基体界面处内聚力模型参数,利用双线性和自定义指数型损伤内聚力模型模拟了AP颗粒和HTPB基体黏接界面处损伤的萌生、发展、聚合直至宏观裂纹破坏的过程。通过数值仿真与实验结果对比发现,指数型损伤内聚力模型比双线性模型能更准确描述推进剂单轴拉伸过程中颗粒与HTPB基体界面间脱黏过程。最后对比了多阶段加载实验结果与仿真结果曲线,发现两者变化趋势基本一致,最大偏差仅为10%,验证了所建细观模型的可靠性及反演所得界面参数的准确性。      

基于正激波位置计算的进气道/发动机系统实时仿真模型

研究了超声速条件下,基于进气道内正激波位置计算的进气道/发动机系统综合仿真问题。提出了将超声速进气道分为内、外两部分独立计算的方案,外部通过斜激波计算公式获得进气道内部的边界条件,对进气道内部采用准一维CFD模型模拟背压对正激波位置的影响。相比二维CFD计算,在保证参数精度要求的前提下有效缩短了计算时间,提高了仿真实时性。进一步将进气道模型与发动机部件级模型匹配运行,构建了进气道与发动机综合实时仿真模型。进气道/发动机综合模型的仿真表明,其能够准确模拟进气道正激波位置对涡扇发动机的动稳态影响,当来流发生扰动时,通过喷管喉道面积和压缩部件导叶角的快速调节,能够抑制正激波位置的变化,并使推力下降量减少50%。      

叶片丢失后发动机整机响应模拟试验与仿真

依据大涵道比涡扇发动机结构设计了叶片丢失试验台,开展了一系列模拟叶片丢失试验,并采用显式有限元方法进行了数值仿真,研究了发动机叶片丢失后整机结构响应与载荷传递规律。结果表明:丢失叶片与机匣存在叶尖与叶身两次撞击,对应的加速度曲线存在两个响应峰值;转子转速越高,加速度响应幅值越大。叶片飞断后转子不平衡载荷传递路径为前轴承支承-中轴承支承-中介支板-机匣结构;叶片撞击机匣导致的冲击载荷则由风扇机匣向后传递,最终传给吊装结构;发动机承受的载荷是由不平衡和冲击影响耦合得到,其中冲击载荷为主要部分。该研究为掌握真实发动机叶片丢失下整机响应规律提供了试验模拟方法与数值仿真分析工具。      

基于CFD和CSM耦合的通用静气弹分析方法

提出了一种适用于有限元精细化建模的流固耦合插值点选择方法,通过RBF(径向基函数)方法实现流固耦合面的数据交换,实现了基于CFD/CSM(computational fluid dynamics/computational structural mechanics)耦合的通用非线性静气弹分析方法。以HIRENASD(high Reynolds number aero-structural dynamics)风洞试验模型为验证对象,数值结果很好地与风洞试验结构变形、气动压力分布吻合,验证了所发展非线性CFD/CSM耦合静气弹求解器的精度。详细研究了HIRENASD模型在大迎角(AOA)流动下的静气动弹性特性,以及该模型弹性变形对机翼气动特性影响规律。研究表明:HIRENASD弹性模型变形后其升力小于刚性模型;在小迎角范围内刚性、弹性模型升力差随迎角增大呈线性增长;当迎角大于4°后,升力差先减小后基本保持不变,呈非线性关系。      

热冲击后气膜冷却叶片结构强度流固热耦合仿真

研究了冷气流量对气孔周围热应力的影响,为气膜冷却叶片可靠性设计提供参考。改变气孔的孔径,并建立有限元模型,结合有限元/边界元理论,通过流固热三场耦合技术获得热冲击后的叶片最大温度、温度不均衡程度及最大热应力。研究表明:增加冷气量有利于改善叶片冷却效率降低叶片温度,但也会使叶片温度不均衡程度增加,加剧尾缘气孔内的热应力载荷;增加前缘气孔直径可提升66%的平均冷却效率,有利于减缓气孔内的热应力,增加尾缘气孔的直径对冷却效率及热应力的影响均较小。此外,数值计算结果与试验及解析解较为吻合,对于气膜冷却叶片结构设计具有参考价值。      

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符;反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。      

大涵道比风扇/增压级进气畸变数值模拟研究

为探究进口总压畸变下大涵道比风扇/增压级内部流场的主要特征,基于三维彻体力模型的思想,开发了一套能够实现风扇/增压级内外涵联算功能的三维数值计算程序。利用该程序模拟了某大涵道比风扇/增压级在周向总压畸变进气下的三维流场。计算结果表明：大涵道比风扇单转子不同叶高处的畸变传递特征差异较大,转子出口总压畸变强度由叶根到叶尖逐渐降低,在叶尖处衰减为最小值1.5%;在转子出口相应诱导出的总温畸变强度由叶根到叶尖逐渐升高,在叶尖处达到最大值1.4%;进口周向总压畸变导致风扇转子总压比下降0.5%,而风扇转子出口形成的总压总温复合畸变导致增压级总压比下降2%;总压畸变在增压级中呈逐级衰减趋势,而高温畸变区的周向范围在逐级增加。