悬停共轴双旋翼干扰流动数值模拟

基于非结构动态嵌套网格方法,通过求解三维非定常Euler方程对悬停共轴双旋翼直升机的复杂流场进行了数值模拟研究.双旋翼之间的干扰使得两者的拉力性能下降,上旋翼对下旋翼的干扰程度远大于下旋翼对上旋翼的影响程度.下旋翼处的下洗速度比单旋翼的要大得多,而上旋翼的下洗速度比单旋翼的略大一些,几乎相同.桨叶总星巨角增加,则下洗速...     

气体润滑弹性波箔片轴承动态特性分析

通过摄动方法并结合气体润滑压力控制Reynolds方程和气膜压力、气膜厚度及箔片变形等变量的Taylor展开式,推导出弹性波箔片轴承气膜动力特性系数计算微分方程组.同时,考虑气体可压缩性等因素,建立气膜厚度、箔片变形和气膜压力三者之间关系的数学模型方程组.运用Newton-Raphson迭代法和有限差分法耦合求解气膜动...     

冲击＋同向对流＋气膜冷却的三维壁温计算与分析

采用Fluent 6.3软件对冲击＋同向对流＋气膜冷却结构的壁温进行了流固耦合传热计算,获得了该复合冷却结构的壁温和冷却效果分布规律,并将计算结果与试验结果进行了对比.结果表明：冲击＋同向对流＋气膜冷却具有壁温分布均匀、冷却效果好等优点.计算结果与试验结果的分布趋势一致,最大相对误差为3.6%.     

结构隔绝内-外空间的流-热-固耦合仿真应用

为验证和指导高速飞行器的防隔热设计,准确地模拟气动热产生的热量穿透防隔热材料进而影响舱内温度空间分布和时间变化的过程,研究了一种同时求解机体外流场及气动热、机体结构传热及舱内流场温度场仿真计算方法,其中的传热方式包括热传导、热对流及热辐射。采用两套计算模型、两种求解器、一个数据交换文件的计算结构,构建了一种针对流场-热-结构的多场耦合分析方法,实现了对固体隔绝内外流场温度动态变化问题的仿真分析。最后通过计算实例验证了整套计算方法,得到的飞行器舱内温度变化特性能够用于指导高速飞行器的防隔热设计。     

高超声速复杂气动问题数值方法研究进展

高超声速流场具有复杂流动特征,其中真实气体效应、磁流体干扰效应和力热结构耦合效应等对气动力分析产生了重要影响。将流体力学研究扩展到分子动力学、电磁流体力学以及流固耦合等交叉学科领域,这给数值模拟方法带来了巨大挑战。针对高超声速气动力/热分析的热点问题,重点关注高温效应与低密度流动效应、磁流体干扰效应和力热结构耦合效应等,结合算例分析了相应的数值求解技术;在气动热方面主要比较了3类求解方法(纯工程方法、纯数值方法和基于Prandtl边界层理论的方法),并给出了相应算例;对于气动力/热/结构耦合问题,从耦合模型及耦合计算方法两方面开展了分析。最后指出了高超声速复杂气动问题数值求解技术未来需重点关注的几个方面。     

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

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

基于退化模型的失效机理一致性检验方法

有效的加速试验必须保证产品的失效机理不变。针对加速退化试验(ADT)可靠性建模的特点,提出了一种基于退化模型的失效机理一致性检验方法。以基于随机过程和基于退化轨迹拟合的退化模型为例,通过加速系数的2种等效定义分别建立了模型参数与失效机理不变的内在联系,进一步推导出其参数在不同应力下须满足的关系。在合理假定检验样本服从正态分布的基础上,利用t统计量对样本的一致性进行检验。通过2个ADT实例验证了所提检验方法的合理性与可行性。     

基于波动耦合的复杂结构动力响应分析与研究

有限元法(FEM)适合低频分析,统计能量分析(SEA)方法适合高频分析,而一类较宽频带(中频段)的声振响应问题不适合单独用FEM或SEA方法解决。基于波动耦合的混合FE-SEA方法兼顾了FEM与SEA方法的优点,可预测中频段含一定不确定性的声振系统的稳态响应。本文阐释了混合FE-SEA方法的理论,分别通过典型结构的蒙特卡罗仿真与复杂系统的噪声试验两案例对该方法进行了验证,其中梁-板典型结构能量响应计算值与蒙特卡罗仿真值具有较好的一致性;复杂飞行器系统的振动响应与试验结果较为吻合,舱内噪声声压级(SPL)误差小于3dB,证明了混合FESEA方法的有效性,解决了FEM与SEA方法在中频复杂声振响应中的局限性这一问题。     

直升机的气动弹性问题

直升机的气动弹性问题与固定翼飞机不同,不仅要考虑单片桨叶,更要将旋翼视为一个整体,考虑其动态入流、尾迹影响以及旋翼与机身之间的相互耦合等。就单片桨叶而言,在结构动力学上,需要考虑离心力场、几何非线性以及桨叶的非线性挥舞-摆振-扭转耦合;在气动力上,需要考虑动态入流以及桨尖处可能的失速效应,本质上属于非线性气动弹性力学范畴。由于旋翼气动力通常是以周期形式通过旋翼轴传给机身,并引起机身振动,而机身运动又通过改变桨叶根部形态反过来影响旋翼的气动弹性特性,这种旋翼/机身耦合问题,也是近年来直升机气动弹性问题研究中的重要方向和热点之一。此外,随着旋翼流场数值分析方法的日趋成熟,采用动态重叠网格或滑移网格方法来实现桨叶运动,并通过动网格技术来实现桨叶的弹性变形,从而实现弹性旋翼流场的数值模拟,目前正呈现出勃勃生机,成为直升机气动弹性研究的又一重要方向和热点。随着各种新构型直升机的相继出现,如倾转旋翼机、前行桨叶概念旋翼(ABC)直升机和复合式直升机等,也带来了新的气动弹性问题。不断发现问题、解决问题,推动本学科持续发展,永远是气动弹性工作者终身奋斗的目标。     

预置接触区长度系数的Spirac切齿调整计算方法

基于参考点处瞬时接触椭圆大体决定接触区长度的方法,提出一种切齿调整计算方法.根据奥利康刀倾半展成法(Spirac)切齿原理,建立摆线齿准双曲面齿轮的数学模型,得到了大轮齿面参考点处瞬时接触椭圆参数;在满足参考点位置、压力角和螺旋角的基础上,将刀倾角也作为迭代变量,并将接触区长度系数等于预置值作为新增的迭代条件,迭代求解切齿调整参数.最后以一对摆线齿准双曲面齿轮副为例,对比该方法与Spirac切齿调整计算结果,验证了该方法的准确性,并对不同预置接触区长度系数得到的齿面进行齿面接触分析.结果表明,该方法解决了反复调整刀倾角的问题,并通过精确控制大轮齿面接触区长度系数等于设计值,实现了对大轮齿面接触区长度的预控.