基于约束预测控制的飞行失控包线保护研究

风切变等大气扰动现象严重威胁飞机起降飞行安全。在受到随机风切变干扰下,离线设计的最优控制器无法保证飞行状态始终处在正常包线内。在大型民机线性变参数模型的基础上,从预防飞行失控的角度研究风切变下飞机进场着陆的预测控制方法。将飞行失控包线量化为状态约束,在此基础上设计在线模型预测控制器。仿真结果表明,利用预测控制的滚动优化策略,可使飞机在进场下滑阶段遭遇风切变时进行良好的下滑道捕获和跟踪,且将飞机状态始终约束在安全包线内。     

穿越微下冲气流的飞翼布局无人机控制方法

微下冲气流是最危险的低空风切变形式,为在起降阶段安全穿越该气流,飞翼布局的无人机控制律应具有快速响应能力和良好的鲁棒性。针对大展弦比飞翼布局无人机舵面附加升力大和低速状态俯仰操纵效能低的特点,提出了舵面附加升力和机体气动力相结合的复合控制方案,改进了以输出误差为参考量的非线性指令分配策略,设计了基于迎角保护的指令分配策略。将风干扰和模型的不确定性视为未知扰动,采用自抗扰控制(ADRC)理论设计飞翼布局无人机非线性控制律,使之对风干扰和模型的不确定性进行估计补偿。仿真结果表明,复合控制与ADRC相结合的方法加速了航迹倾角的单位阶跃响应速度,使上升时间缩短了64%,同时能够实现对风干扰的有效观测和补偿,使高度损失低于2m;能够在风切变中有效保护迎角,使其维持在5.5°以内。因此,该方法能够为飞翼布局无人机安全平稳地穿越微下冲气流提供一种参考方案。     

连续尺度变换方程在稳态及非稳态流动模拟中的应用

为了提高涡粘性假设的湍流模型对于非稳态流动的求解精度,同时兼顾其对于稳态流动的求解性能,将雷诺应力项与连续变换方程(CSSE)结合而形成新的应力项,使其根据流场尺度、网格尺度及Kolmogorov尺度来自动调节当地的应力雷诺应力模化水平,避免网格因素在流场模拟中产生不利影响,改正了混合RANS/LES方法的速度型偏离对数率问题;同时,该方法并未引入显式亚格子模型(SGS),因此回避了亚格子系数确定对于流场模拟精度产生的影响,改善了湍流模型对于流动不稳定性的辨识精度。在湍流平板算例中,CSSE方法计算的边界层速度型精度与雷诺平均方法(RANS)相当,而对于圆柱尾迹的模拟则证明了CSSE方法具有混合RANS/LES方法的优点,即能够准确模拟流动的不稳定性特征。     

航空铆钉的动态力学性能测试

在飞机鸟撞的课题研究中,数值模拟占据着日益重要的地位,而数值模拟的精度在很大程度上依赖于铆钉力学参数的准确性,但是目前航空铆钉的动态力学性能参数还很匮乏。本文测定了航空铝合金铆钉在不同加载速度下的剪切和拉伸力学性能。设计了基于分离式Hopkinson拉杆系统的铆钉动态性能测试装置,并利用此装置对7种不同型号的航空铆钉进行了动态剪切和拉伸试验,得到了每种型号铆钉的动态力学性能。利用电子万能试验机对7种铆钉进行了两种应变率下的准静态剪切和拉伸试验,并且和动态试验结果进行对比。讨论了加载速度、加载形式、铆钉直径以及铆钉形式对铆钉力学性能的影响,并利用扫描电子显微镜(SEM)观察了铆钉剪切破坏以及拉伸破坏的断口形貌,分析了其破坏形式。本文的研究结果对于航空铝合金铆钉在工程中的应用、尤其是在抗冲击领域的应用具有指导意义。     

柔性热膜剪应力传感器水下测量温度修正

柔性热膜微传感器应用于流体壁面剪应力测量时,流体温度对传感器输出有很大影响,因此有必要对工作在非标定温度下的传感器输出信号进行修正补偿。重点研究了水下测量时传感器输出与水温的关系以及温度修正方法。通过分析传感器过热比和流体物理性质与温度的相关性,建立了流体温度与传感器过热比和标定系数的函数关系。在此基础上,提出了一种用于水下剪应力测量的温度修正方法,可以有效减小水温对传感器输出的影响。经实验验证,该方法可以使工作在25℃、28℃水温环境下的传感器输出值与其20℃标定值的相对误差从23．7％和37．1％回落到0．82％和0．83％。     

剪切层与边界层组合流动的线性不稳定性分析

对不可压缩剪切层与边界层的组合流动完成了线性稳定性研究.组合流动的数学模型为Blasius边界层相似解与双曲正切函数的叠加,采用整体数值方法求解组合流动的稳定性方程,并验证了程序的准确性及网格无关性.研究给出组合流动的不稳定模态的辨识,即边界层模态和剪切层模态.在此基础上研究了剪切层对边界层模态不稳定性的影响以及壁面对剪切层模态的影响.由于剪切层的存在,使边界层模态中性曲线向左上方平移,临界雷诺数减小.此外,边界层模态不稳定性得到增强或抑制的影响,取决于扰动频率以及剪切层速度比的变化.组合流动中壁面边界层促使剪切层不稳定性减弱,主要表现在低频区域;而在高频区域,剪切层不稳定性几乎不受壁面边界层的影响.      

高度欠膨胀声速射流的自相似特性

构建了高压气体射流模型和高分辨率的计算网格,对喷压比（NPR）为5.60、流动雷诺数Re为105量级的高度欠膨胀射流进行了三维大涡模拟（LES）计算.讨论了时均的射流近场结构,发现大涡模拟成功捕捉到了高度欠膨胀射流近场的典型波系结构,并与文献结果吻合较好.研究着重定量考察流向速度和氮气质量分数的分布规律,以揭示高度欠膨胀射流的自相似特性.结果表明:在射流核心区之后,不同流向位置上流向速度和氮气质量分数沿径向的分布呈现出与亚声速射流类似的自相似特性.但射流流场开始呈现自相似的位置仍然为超声速,这是高度欠膨胀射流自相似特性的独有特征.提出了高斯拟合公式,在靠近和远离射流中心线的区域内均可较好地表征射流的自相似特性.此外,还考察了射流剪切层的发展特征,比较了计算得到的射流剪切层增长率与先前实验测量结果的差异,并分析了其中可能的原因.      

边界层厚度对腔体气动声学特性影响数值模拟

为了研究来流边界层厚度对开式腔体气动声学特性的影响,基于分离涡模拟方法,计算了来流马赫数为2.0条件下,不同来流边界层厚度与腔体深度比时,长深比为5.88的腔体流动特性,得到了该腔体声压级的频谱特性.计算结果表明:随着来流边界层厚度增加,形成的剪切层稳定性增强,失稳后上下摆动幅度减少,失稳生成的大尺度涡与超声速主流的相互作用减弱,使得大尺度涡发展到腔体后缘时所具有的平动动能和转动动能降低.大尺度涡撞击腔体后缘在腔体内形成的气动噪声的声压级降低,最大减小幅度达7.5dB.同时各阶模态的频率也发生偏移,偏移值在100Hz左右.基于新的假设重新推导了Rossiter公式,明确了经验常数的物理意义,并以此解释了频率偏移现象.      

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。     

在平板不同位置开孔的剪切稳定性分析及应力集中影响

为满足管线通过、减重、维修等需求,通常需要在结构上开孔。基于Patran和Nastran软件,针对剪切载荷下的四边简支开孔矩形平板,对其在不同开孔位置下的稳定性及口边应力集中情况进行仿真分析,并整理仿真分析结果,得到不同开孔位置对矩形平板剪切稳定性和口边应力集中的影响,供结构开孔设计使用。