基于响应面近似技术的平流层飞艇协同设计优化

针对平流层飞艇多学科设计优化问题,引入协同优化算法,分析基于响应面近似模型的协同优化方法CO/RSA的执行步骤和特点,并介绍响应面近似模型的构建方法;建立平流层飞艇各学科模型,并对学科间耦合关系进行分析;运用建立的气动/推进、结构、能源三个耦合子系统的学科分析模型和系统优化模型,以平流层飞艇总质量最小为优化目标,综合考虑浮重、推阻和能源平衡,采用CO/RSA方法对平流层飞艇进行设计优化.结果表明:所建立的平流层飞艇的MDO模型是合理的,CO/RSA算法应用于飞艇总体设计优化是有效的,研究结果可为平流层飞艇方案论证和方案设计提供参考.     

旋转爆震发动机爆震波建立过程实验研究

为了深入研究旋转爆震发动机爆震波建立过程及形成机理,采用小能量火花单次点火的方式进行了一系列旋转爆震发动机起爆实验。发动机采用环缝-喷孔对撞式掺混方式,燃料为H2,氧化剂为空气,实验成功起爆旋转爆震波,并连续旋转稳定传播,爆震波传播频率为5.09~6.45k Hz,传播速度为1286~1644.8m/s。在发动机稳定工作过程中,集气腔与燃烧室相互影响,二者处于平稳的动态平衡。其次,通过对旋转爆震波起爆过程详细分析发现,点火形成的初始火焰在环形燃烧室经历一个类似DDT的火焰发展过程,成功转变为爆震波,且从点火到爆震波建立之间的火焰发展传播过程和时间间隔均表现出很强的随机性。此外,为验证小能量火花点火的可靠性,还进行了小能量点火重复性实验,发现在稳定工况条件下采用小能量点火成功率最高可达100%,各组旋转爆震波传播速度在1440m/s附近波动。     

俯仰姿态保持模式下飞机结冰边界保护方法

在研究飞机结冰机翼和平尾失速机理的基础上,以飞机迎角作为关键参数,飞机俯仰指令作为指令参数,提出一种基于飞机自驾仪的结冰后边界保护方法.通过引入铰链力矩检测模块,提前告警飞机失速,为边界保护提供了裕度.建立飞机纵向动力学方程,针对俯仰姿态保持（PAH）模式下机翼失速进行了仿真计算.结果表明:结冰条件下,通过铰链力矩的检测,飞机能在失速迎角为2°之前,进入边界保护模式,增加飞行安全的裕度.通过控制指令的限制,自动驾驶模式下飞机迎角能保持在失速迎角之内.      

基于自然激励技术的颤振边界预测

为了预测紊流激励条件下机翼的颤振边界,基于自然激励技术提取紊流响应的自由衰减信号,采用矩阵束方法识别模态参数,最后通过Z-W(Zimmerman-Weissenburger)方法计算稳定性判据,拟合判据变化曲线并外推颤振边界.对平板机翼模型进行了数值仿真分析,对单独机翼模型风洞颤振试验数据进行了计算.结果表明:采用自然激励技术与矩阵束方法能够较准确地识别紊流激励响应的模态参数,频率识别误差小于6％,阻尼比识别误差小于30％,结合Z-W方法能够在较低风速较早地预测颤振边界,有助于提高试验的安全性.     

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

基于功能危险性评估的民用飞机驾驶舱门设计探讨

FAA分别于2002年和2008年颁布25-106及25-127修正案,增加驾驶舱门安保事项以提高驾驶舱的安全水平。修正案要求加固驾驶舱门,驾驶舱门必须设计成可以阻止未被授权人员进入驾驶舱并且阻碍危险物体射入以保证机组人员的安全。驾驶舱门必须能检测驾驶舱失压并可以释放门或者卸压板平衡压力差,还应有措施使飞机机组成员在该舱门卡住的情况下能直接从驾驶舱进入客舱。因此,在驾驶舱门上设有一个可以快速释放的紧急出口。根据飞机驾驶舱门结构及功能的主要特点,对驾驶舱门功能故障进行探讨,同时对功能故障引发的危险性等级进行划分,以确定安全指标,进而辅助驾驶舱门设计工作。     

亚声速进气道设计敏感参数仿真分析

对不同气动布局的两种亚声速进气道建模,并针对两种进气道的敏感参数进行分析。通过对背部双发布局飞机进气道的隔道高度和喉道设计马赫数参数的研究可知,隔道高度增大或者喉道设计马赫数减小,使得进气道出口总压恢复系数增大,增量约0.5%。而对背部单发复杂形状进口(M型进口)布局飞机进气道的喉道形状研究显示,喉道形状越接近进口形状,进气道性能越好,较其他喉道形状进气道,M型喉道进气道出口总压恢复系数最高,最大增量约为1.5%。     

Numerical analysis and optimization of boundary layer suction on airfoils

Numerical approach of hybrid laminar flow control(HLFC) is investigated for the suction hole with a width between 0.5 mm and 7 mm. The accuracy of Menter and Langtry’s transition model applied for simulating the flow with boundary layer suction is validated. The experiment data are compared with the computational results. The solutions show that this transition model can predict the transition position with suction control accurately. A well designed laminar airfoil is selected in the present research. For suction control with a single hole, the physical mechanism of suction control, including the impact of suction coefficient and the width and position of the suction hole on control results, is analyzed. The single hole simulation results indicate that it is favorable for transition delay and drag reduction to increase the suction coefficient and set the hole position closer to the trailing edge properly. The modified radial basis function(RBF) neural network and the modified differential evolution algorithm are used to optimize the design for suction control with three holes. The design variables are suction coefficient, hole width, hole position and hole spacing. The optimization target is to obtain the minimum drag coefficient. After optimization,the transition delay can be up to 17% and the aerodynamic drag coefficient can decrease by 12.1%.     

二元双喉道射流推力矢量喷管流动参数影响的数值研究

采用数值模拟方法研究了不同流动参数对二元双喉道射流推力矢量喷管(Dual-throat fluidic Thrust-vectoring Nozzle,DTN)内流特性和推力矢量控制效果的影响。结果表明,DTN在非推力矢量情况下,NPR在3~4范围时,推力系数较大,达到0.968,而流量系数较小,仅为0.93;NPR再增大,推力系数迅速下降。在推力矢量情况下,落压比一定时,随着次流流量比的增加,推力矢量角增加,而流量系数、推力系数、推力矢量效率减小;次流流量比一定时,随着落压比的增加,推力矢量角减小,系统推力系数先增加后减小,流量系数略微增加。     

民用飞机壁板蒙皮及长桁布置结构优化设计

机身壁板是飞机结构设计的重要承载组件,轻量化、高效率、共通性设计及优化是民机设计关注的重点。首先提出一种耦合ABAQUS的Buckle分析及ISIGHT优化的设计方法,利用自编子程序获取ABAQUS屈曲特征值,将特征值输入ISIGHT中计算临界屈曲载荷,同步更新变量参数及ABAQUS文件并提交计算,迭代分析直至优化流程结束。采用上述方法考虑轴向压缩载荷情况,以壁板整体重量最小为优化目标,疲劳应力值为约束条件,对单曲度金属机身壁板的蒙皮厚度,长桁数量及长桁截面厚度等几何参数进行优化。在满足壁板结构承载能力及总重量最小条件下,综合考虑结构载重比,临界应力及壁板加筋比,对比分析出一组最优参数,并与工程算法结果对比吻合程度较好,两者相对误差为3.73%。该优化思路实现FEA平台与优化工作一体化,可用于复合材料壁板设计及结构件减重优化工作,一定程度上可缩短零组件设计周期。