悬停状态直升机桨叶扭转分布的优化数值计算

建立了一套基于高精度计算流体力学(CFD)技术和代理模型优化算法的旋翼气动外形设计方法。在该方法中,旋翼流场气动性能的计算采用了基于Navier-Stokes/Euler方程的CFD方法,并根据流场特点、精度和效率的要求采用Baldwin-Lomax(B-L)湍流模型,通量计算采用Roe-MUSCL格式进行。为提高网格生成质量和便于流场控制方程的求解,将流场分成两个区域,即围绕旋翼的黏性区和无黏的背景网格区。其中,桨叶网格使用了基于二维翼型网格的参数化方法生成,数值计算结果表明该方法有效地提高了网格生成质量及效率。在参考旋翼流场及桨叶细节流动分析的基础上给出设计变量及范围,有效减小了优化问题的规模;为满足优化和机理分析的需要,将基于置换遗传算法优化的拉丁超立方(PermGA LHS)方法和径向基函数(RBF)的代理模型优化方法引入到桨叶外形的优化设计中。首先以Helishape 7A旋翼为算例,检验了数值模拟方法的准确性。然后,应用所建立的优化方法针对旋翼负扭转分布进行了优化计算,结果表明优化后的旋翼悬停气动性能比优化前有了明显提高。     

计入畸变修正的旋翼尾迹前飞状态稳定性分析

现有的旋翼尾迹稳定性分析方法在涡线过近处会得出不现实的过大发散率结果,因而存在一定的局限性。为了消除这种前飞状态下尾迹高度畸变对其稳定性分析所造成的不合理干扰,提出了一个旋翼尾迹在正弦扰动下,结合涡线畸变修正的线性化稳定性分析方法。在该方法中,将尾迹涡线离散为一系列涡元,并对其端点施加沿尾迹涡线分布的正弦扰动;考虑了桨尖涡的自诱导、互诱导以及与桨叶的干扰作用,并创新地引入在伪涡核对稳定性分析中,因尾迹高度畸变产生的不合理干扰而进行的修正。以H-34型直升机旋翼模型为例,在前飞状态下,重点计算分析了尾迹向后飘移、涡核大小和前飞速度等对旋翼尾迹稳定性的影响。研究结果表明:引入伪涡核模型可以有效地修正在涡线过近处所产生的不现实稳定性分析结果;前飞速度、涡核大小等因素会在一定条件下有限地影响尾迹的不稳定性,但都不改变前飞状态旋翼尾迹的本质不稳定性。     

复合材料机身壁板的纵向连接设计与失效分析

按照连接设计准则及机身压差载荷水平,开展了复合材料机身壁板的纵向连接设计研究。为提高壁板多钉连接结构分析精度及设计效率,发展了一种基于Fastener单元的钉群载荷计算方法,在此基础上结合单钉失效分析模型,提出了一种壁板多钉连接区的失效评估方法。首先,通过与试验数据对比,验证了采用Fastener单元求解钉群载荷的可行性;然后运用Fastener单元分析壁板连接结构的钉载分配;最后基于钉载分析结果,对局部危险区域采用单钉模型进行失效载荷计算并进而评估壁板连接区的失效载荷。本方法特别适用于快速、有效地校核多钉连接区的连接强度。     

双层壳型冲击/气膜综合流量系数实验研究

为了研究新型高效涡轮铸冷叶片内部冷却通道中的流动特性,以铸冷叶片典型结构——双层壳型冲击/气膜复合冷却结构的平壁模型为研究对象,实验研究了该种冷却结构的综合流量系数Cd随冲击雷诺数Red(600～2800),主次流吹风比M(0.7～7),冲击间距比H/d(0.5～1.25),冲击孔和气膜孔间距比P/d(0,3,4)等参数的变化规律。研究结果表明:(1)Cd随着冲击雷诺数Red,吹风比M的增加而增大,在吹风比小于3工况下更为显著;(2)在相同的气动参数条件下,综合流量系数Cd随着孔排距比P/d的增加而增大,随着冲击间距比H/d的减小而显著降低。     

预旋进气小尺寸涡轮叶片冷却的流场研究

为了了解和掌握一种具有直通式冷气预旋进气系统的小型燃气轮机涡轮转子叶片的流场,在旋转雷诺数Reθ=4.66×106和冷却空气的无量纲质量流量Cw=1750时改变预旋角θ的大小,使其在15°～90°变化,通过数值研究得到了预旋角对涡轮盘腔、连管和涡轮叶片内冷却空气的流动以及叶栅通道中燃气的流动的影响。结果表明:(1)预旋角的变化会改变涡轮盘腔、连管和涡轮叶片冷气进口附近局部区域的流场,但是对涡轮叶片内其它区域和叶栅通道中的流动基本没有影响。(2)随着预旋角的增大,涡轮盘腔内预旋进气冷气射流的轴向穿透深度先增大后减小;当θ<45°时冷却空气沿外围屏流向转盘接收孔,而当θ>45°时冷却空气沿内围屏流向转盘接收孔;气流的周向速度随着预旋角的增大而减小。(3)垂直进气时连管内存在多个回流区和很大的涡流,流动损失较大,而采用预旋进气能够减弱或消除这些流动结构,存在最优预旋角θopt,θopt≈45°,此时连管的有效流通面积最大。     

基于Matlab／Simulink的航空发动机部件级建模与分析

借助Matlab搭建完整的双轴涡轮风扇发动机模型,对各部件进行封装,利用Simulink显示模块观察仿真结果。在动态仿真时,引入容积动力效应,避免反复迭代过程,提高模型实时性。与GSP仿真软件结果进行对比,结果表明：所建立的模型具有较高精度,可以准确捕捉航空发动机的动态响应和稳态响应。     

圆转矩形喷管出口宽高比对射流雷诺剪应力的影响

采用热线风速仪,利用单斜丝,对宽高比W/H分别为1,4,8,12,16的5个圆转矩形收敛喷管和一个轴对称喷管的射流对称面上雷诺剪应力分布特性进行了实验研究.研究发现:在喷口下游不同截面上,射流宽、窄对称面上的雷诺剪应力沿径向均先缓慢增大,到达射流边界后迅速减小,射流边界逐渐沿径向外移.矩形喷管射流相比轴对称射流具有较强的旋流,雷诺剪应力较大,且随着宽高比增大,旋流强度增大,剪应力也逐渐提高,导致了射流与外流掺混增强.宽高比大于8以后,增大幅度逐渐减小.射流宽、窄对称面上的分布规律相同.     

基于实测的古斯塔夫飓风非平稳特性分析

根据美国路易斯安娜州立大学飓风研究中心对古斯塔夫(Gustavo)飓风的实测资料,首先采用传统方法统计出不同高度处强风的基本参数特性,然后运用轮次法对不同高度、时距和起始点的风速样本进行非平稳特性检验,最后基于希尔伯特-黄变换(HHT)方法对于呈现非平稳特性的风速序列进行经验模式分解(EMD)和时-频-谱联合特征分析.分析发现,飓风的湍流强度和积分尺度较良态风场相差较大,并且其平均风速也极不稳定;通过轮次法对不同高度和时段的风速时程检验发现,非平稳特性与高度和统计时段大小密切相关,相应的10min时距的风速样本也存在非平稳特性;HHT对非平稳风速序列分析结果表明,飓风的非平稳特性主要是由时变的平均风速引起,其能量和频率随着时间的变化明显,并且主要能量集中在0.05~0.2Hz频段内,但EMD分解和Hilbert谱结果表明高频能量同样不能忽略.最终得出的结论为:基于传统的分析方法可以很好地获得基于平稳假定的飓风统计参数,但不能获取其非平稳信息,而轮次法和HHT方法是识别非平稳和分析其时频谱联合特性的良好方法.     

行星齿轮非线性振动系统参数稳定域的计算方法

研究了行星齿轮非线性传动系统参数稳定域计算的一般方法.该方法通过选取合理的失稳阀值,根据考查参数域内系统的运动状态选取合适的数值积分时间段,以循环套嵌的手段计算考查参数在各自范围内不同组合下的系统位移响应最大值,比较失稳阀值以判稳,参数稳定域的图形输出等5个步骤完成对行星轮系参数稳定域的计算.最后,以四自由度行星轮系纯扭转非线性振动模型为例,以行星轮输入转速、系统的齿侧间隙以及齿轮副的啮合阻尼系数为考查参数,分别计算得到了系统的单参数稳定域、双参数稳定域以及三参数稳定域,为行星轮系的设计取值提供了重要参考.      

可调静子旋转凸台位置对端区流动的影响

对包含凸台的可调静子模型进行简化,研究了不同凸台位置对端区流动的影响.通过实验测量与数值模拟对比可以发现:可调静子端区流动损失主要是由叶片气动负荷引起的泄漏及气流绕过凸台造成的掺混2个部分组成;通过合理优化设计,将凸台置于气动负荷高的叶片前缘并覆盖尽可能多的前缘间隙,能够有效降低损失,改善气动性能.