考虑光压力和热载荷的太阳帆薄膜力学行为分析

基于冯卡门板理论描述薄膜结构大变形,耦合光压力与热载荷,构建10μm量级薄膜非线性动力学主控方程。采用伽辽金离散方法和龙格-库塔推进方法数值求解方程时域解,并引入最大李雅普诺夫指数判别薄膜混沌运动。分析表明,法向光压力远大于切向光压力并抑制薄膜动力学颤振,仅计入光压力时,太阳帆膜呈静变形特征。当计入的光辐射热载荷达到临界值后,考虑切向光压力与法向光压力的薄膜都呈现动力学混沌运动。综上,光压力抑制太阳帆薄膜颤振,热载荷导致薄膜进入混沌运动。     

斜激波总压损失率极小值理论解与物理意义

以法向马赫数作为激波强度表征量,对斜激波关系式进行重新推导,得到了穿过激波总压损失率极小值的理论解。控制方程表达式为激波角对物面角的线性函数。依据斜激波总压损失率极小值解析公式,首先,绘制了针对超声速流总压损失率应用的楔形角-激波角-马赫数的斜激波效率图。其次,通过生成斜激波三维总压损失率等值线图,呈现了总压损失率在楔形角-激波角-特征马赫数空间上的分布规律。此外,利用斜激波效率图,揭示了等总压损失率条件下马赫数与激波角的对称双解现象。     

进口导流叶片对转子叶片颤振特性的影响

上游静子叶片的尾迹扰动是下游转子叶片发生强迫响应问题的主要原因,但通常认为不会对颤振特性产生影响。运用CFD技术求解非定常流场,用能量法对典型涡轮风扇NASA Rotor 67进行了气动弹性分析。结果表明:不加进口导流叶片时其叶片最有可能在第二阶模态、60°叶间振动相角的情况下发生颤振;加上IGV(进口导流叶片)后,IGV的尾迹会显著改变转子叶片的颤振特性,当转子叶片的最危险颤振模态频率接近IGV尾迹的扰动频率时,即使单独转子叶片十分稳定,但在尾迹的激励下颤振仍然可能发生;当该频率远离IGV尾迹扰动频率时,尾迹扰动同样会使转子叶片气动阻尼降低。该结果意味着尽管IGV的引入不会明显改变压气机性能,但有可能会对其颤振特性带来明显的负面影响,需要在涡轮机设计中予以考虑。     

斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。     

内埋式弹舱舱门气动载荷计算分析研究

采用中心有限体积法进行空间离散和四步Runge-Kutta显式推进方法进行时间推进数值求解了非定常欧拉方程;采用Delaunay方法生成内埋式弹舱舱盖打开时的非结构网格.数值模拟了内埋式轰炸机弹舱单开舱门在分别打开8°和45°时的复杂流场.通过对上述复杂流场的分析,发现作用于舱门的气动载荷具有使其趋向于关闭的特性,且有随打开角度的增大而减小的趋势.     

复合材料超临界机翼跨音速颤振实验与数值计算研究

颤振课题是飞机设计过程中常常遇到的一个关键技术问题。以支线客机ARJ21超临界机翼颤振特性研究为背景,在俄罗斯TsAGI的T-106风洞中完成了该复合材料机翼跨音速颤振实验,基于N-S方程和无限插值方法（TFI）生成三维贴体运动网格对超临界机翼跨音速颤振进行了并行计算。结果表明：复合材料的超临界机翼在跨音速区域具有跨音速颤振＂凹坑＂现象;与风洞实验结果相比,有较好的一致性,为使用超临界机翼的运输类飞机跨音速颤振特性预计提供了一定的参考。     

基于非结构运动对接网格的旋翼前飞流场数值模拟

发展了一种基于线性弹簧理论的非结构运动对接网格方法,通过求解非定常Euler方程模拟了旋翼的前飞流场。利用运动区网格的转动模拟旋翼的旋转运动,并运用线性弹簧动态网格技术处理旋翼的挥舞和变距运动,方程求解中采用一种精度较高的插值和重构方法。算例表明,本文的方法能够有效地模拟旋翼的前飞流场。     

两种跨声速气动弹性问题分析研究

在非结构运动网格基础上,采用中心有限体积法进行空间离散和双时间方法进行时间推进求解非定常欧拉方程.通过与气动力方程的联立求解,在时域内用四步龙格-库塔方法求解结构运动方程.分析和研究了二维嗡鸣和三维机翼颤振这两种跨声速非线性气动弹性问题.二维嗡鸣问题的研究考虑了翼面-舵面系统的缝隙间网格运动、缝隙对嗡鸣的影响和扰流片对嗡鸣的抑制.耦合多自由度Lagrange结构运动方程数值模拟了三维机翼的颤振问题.通过跨声速标模算例AGARD445.6机翼的颤振计算,计算的颤振临界速度与实验值有5%左右的误差,验证本方法的正确性.由于本方法是在对外形具有良好普适性的非结构动网格基础上完成的,具有良好的工程实践价值.     

典型超声速/高超声速动导数计算方法研究

以当地流活塞理论为基础完成超声速/高超声速动导数计算公式推导,发展了一种高效的动导数计算方法.选取1个超声速、2个高超声速动导数算例标模,计算了超声速/高超声速下飞行器动导数变化规律.研究结果表明:BFM超声速流动算例的俯仰组合动导数绝对值随马赫数的增大而减小,计算重心位置后移亦会导致俯仰组合动导数绝对值减小,纵向动稳定性降低;在0° ～20°迎角范围内,高超声速流动中的尖锥和钝头旋成体俯仰组合动导数绝对值随迎角的增大而增大,纵向动稳定性增强;所提方法计算时间约为双时间动导数计算方法的1/27.     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。