二维叶栅气固耦合颤振分析

引入叶片结构动力学方程,通过数值求解气动/结构耦合方程体系,研究了气固耦合条件下叶栅内的二维非定常粘性流动,并初步进行了叶栅的颤振分析。针对两种叶型分别计算了不同折合频率下的流场,通过比较振动叶栅相关参数随时间变化的曲线,表明折合频率是决定叶片颤振与否的重要因素。      

风扇/压气机叶型厚度对颤振特性的影响

采用计算流体力学与结构动力学相结合的方法,数值模拟了三维振荡叶栅非定常黏性流场;通过对叶片表面非定常气动力及其所做非定常气动功的计算分析,采用能量法对叶片颤振与否进行预估判断。针对三维直叶片和风扇转子叶片,通过调整叶片的相对厚度,研究了叶片厚度变化对风扇/压气机颤振特性的影响。此外,通过对振荡叶栅非定常流场结构的研究,发现了叶片吸力面气流分离与叶片振动之间的耦合关系。本文的研究在风扇/压气机设计中,可用于评估最大相对厚度等叶型结构设计参数对气弹稳定性的影响,对叶轮机颤振机理研究具有一定的参考意义。     

基于CFD技术的高效叶栅耦合颤振分析方法

为了考虑叶轮机叶片结构与流体之间的耦合效应,同时提高叶轮机颤振数值研究的效率,发展了一种基于非定常气动力降阶模型(ROM)的叶栅耦合颤振分析方法。该方法运用时域计算流体力学(CFD)技术计算少数几个叶片的非定常气动力,通过系统辨识及一些假设构建整个叶栅振动的非定常气动力降阶模型,并在状态空间耦合叶栅结构动力学方程建立叶栅气动弹性方程,采用特征值和时域仿真分析该系统稳定性。运用该降阶耦合方法对STCF4(Standard Test Configuration 4)以及NASA Rotor67叶栅系统的稳定性进行了计算。通过与直接计算流体力学/计算结构动力学(CFD/CSD)耦合方法和非耦合方法计算结果的比较验证了该方法的准确性,且该降价耦合方法的计算效率相对于直接CFD/CSD耦合方法提高了1~2个量级,为叶轮机气动弹性参数研究、失谐研究以及多模态耦合计算等提供了便利。     

非失速二维振荡叶栅非定常流动数值模拟研究

进一步完善了非定常雷诺平均N-S方程求解程序,借助快速网格生成方法,对周期性俯仰运动下NACA0012叶栅的非失速绕流问题进行了数值模拟研究,并分析了振荡叶栅的迟滞效应以及俯仰运动各个参数对迟滞效应的影响.数值模拟结果与实验结果对比表明,本文发展的数值模拟方法是可靠的;结果还表明,随着振幅,折合频率的增加迟滞效应增强,平均攻角对迟滞效应无影响;保持无量纲参数--振幅和折合频率的乘积--恒定时,迟滞效应相同.     

叶片失调及盘片耦合对叶栅颤振边界的影响

 本文用特征值法研完了叶片频率失调及盘片耦合对叶栅颤振边界的影响。计算结果表明:失调可延迟颤振边界,在一定的失调幅度范围内,随着失调幅度的增大,颤振边界向后移;不同的失调阶次对颤振边界的影响不同,影响的大小取决于颤振时的振动阶次。此外,计及盘片耦合效应后,失调叶栅的颤振边界前移。     

基于CFD/CSD时域耦合方法的多通道叶栅颤振分析

为了对叶轮机多通道叶栅颤振进行研究分析,基于时域计算流体力学(CFD)方法,耦合简化的结构动力学方程,运用杂交的线性多步方法实现了对多通道叶栅的时域耦合颤振仿真.该耦合方法通过分析叶片的动力响应特性,获得叶栅的气动弹性稳定性特征.以标准叶栅颤振模型——STCF4 (Standard Test Configuration 4)为算例,计算结果表明该耦合方法与经典非耦合方法在对失稳区域的预测结果吻合得很好,证实了该耦合方法的正确性.运用该耦合方法计算了不同通道数下的自由响应,对比发现真实的叶轮机颤振与外流机翼颤振失稳的不同之处在于叶轮机叶排颤振失稳在很大程度上是一个多模态同时失稳问题;且由于出现不稳定的分支频率和阻尼都接近,叶片的时域响应会出现典型的葫芦节型发散.     

带扰流片的二维翼型非定常计算及分析

通过求解二维非定常可压N-S方程,研究了带扰流片的NACA0012翼型绕流问题.首先计算了扰流片固定在三个不同张开角度时的流动,气动力系数变化规律和已有实验数据对比,结论一致.并着重描述了扰流片张开时,背风区涡形成、发展和脱落过程,分析了非定常气动力产生的机理.为了进一步探讨振动扰流片对主翼的影响,采用Chimera重叠网格技术处理多体相对移动问题,分析了振动频率和振动幅度两个因素对时均气动力的影响.     

方向舵蒙皮的壁板颤振分析

针对某飞机出现的方向舵蒙皮壁板颤振故障引发的振动疲劳裂纹问题，进行了相关的分析计算和修改方案。利用MSC．NASTRAN软件，建立了计算壁板颤振的有限元模型，取前15阶固有振动模态构建模态坐标系下的颤振方程。选用ZONA51超声速气动力计算程序分析壁板的气动力，采用空气密度折半的方法来模拟壁板单面承受气动力的情况。计算结果表明，方向舵蒙皮发生了壁板颤振。通过对方向舵蒙皮结构修改方案的进一步分析，表明修改方案能够满足壁板防颤振设计的要求，并且具有相当大的安全裕量。     

复合材料翼面结构的自振和颤振分析及优化设计

 建立复合材料翼面结构的有限元静力分析模型。采用柔度方法作静力模型到动力模型的转变,进行结构固有特性分析,并研究了频率约束下的优化设计。运用亚音速马蹄涡——振荡偶极子格网法计算了谐振翼面的三维非定常气动力,并对两种真实翼面进行了颤振分析,结果令人满意。推导了非定常气动力和颤振响应参数对设计变量的导数,研究了颤振约束下的优化设计。     

失谐叶栅的受迫振动响应特性分析

采用基于计算流体力学(CFD)方法的降阶气动力模型并耦合结构运动方程,实现了存在外激励载荷时失谐叶栅受迫振动响应的快速分析。针对典型的跨声速叶栅,通过求解其位移响应幅值较系统地研究了失谐方式、失谐强度和叶片质量比对失谐叶栅受迫振动响应幅值的影响。研究表明文中刚度失谐形式可以改善叶栅振动的稳定性,同时导致系统受迫振动响应局部化程度的增加,并且受迫响应的最大振幅放大因子随失谐强度增加或者质量比降低存在先增大后减小的一个峰值,不同失谐形式则对这个峰值的大小有着明显的影响。由于该方法可高效地分析失谐叶栅受迫振动各参数对模态局部化的影响,在工程上有一定的应用价值。     

不同迎角的翼型气弹特性风洞实验研究

基于可在不同迎角下作沉浮、俯仰两自由度运动的翼段振动装置,在低速风洞中分别针对普通薄翼型NA-CA0012和风力机翼型NREL S809进行气动弹性测试,得到不同实验状态的气动弹性振动时域响应。分别观察到经典颤振和失速颤振现象,并证明了迎角改变对两种翼型颤振特性的影响。     

小展弦比压气机转子叶片颤振研究

针对小展弦比压气机转子叶片特点,在能量法计算基础上提出了两种颤振评判方法,即首先利用结构动力学分析结果对叶片频率分离度、振型和共振进行分析,判断某压气机小展弦比转子叶片存在高阶耦合颤振的可能;然后将结构动力学得到的振动位移插值到气动网格上,进行三维非定常黏性流场计算,得到非定常气动功,判断叶片是否存在失速颤振的可能.最后综合利用两部分分析结果和试验及断口分析情况,综合给出该小展弦比压气机叶片因失速颤振和耦合颤振共同作用,导致叶片在短时间内出现了裂纹.      

利用CFD计算技术进行机翼模型气动特性分析

计算流体动力学方法（CFD）已经成为空气动力学领域重要的计算方法。利用该方法进行机翼模型Naca0012的数值模拟,包括计算模型的气动力系数以及气动导数。采用了效果很好的新型网格控制算法,通过自有程序处理流固耦合边界问题以及商用软件的求解器模块进行整体区域的数值计算。在证明了方法的有效性之后,又通过多种振幅奈件下的数值模拟结果,对以往广泛采用的理论方法的有效性进行了探讨,研究了机翼模型气动导数随振幅变化而变化的情况。     

半模机翼振动的风洞实验技术研究

飞行器机模型振动是影响其风洞实验数据可靠性的主要因素之一。为此,通过专门研制的半模振动模型,采用对称翼型NACA0012和层流翼型NACA64-210两种半模机翼进行了低速风洞振动实验。选取五种激振方式,用直接测力法得到模型在静态和不同激振方式下升力特性变化曲线,研究了振动频率、自然转捩和固定转捩、数据采集方式等参数对机翼气动特性测量结果的影响。结果表明：振动对半模机翼气动特性的影响程度因不同翼型构型、不同采集方式和翼面不同流动模式等而产生不同效果。结论对提高飞机模型风洞实验大迎角振动情况下的测量精准度有参考价值。     

半模机翼振动对气动性能影响的风洞试验研究

为了验证翼型振动测压试验结果和有关振动数值模拟文献的计算结果,设计完成了对称翼型NA-CA0012和层流翼型NACA64-210两种半模机翼的低速风洞振动试验。通过专门研制的半模振动模型,选取5种激振方式,用直接测力法得到了模型在静态和不同激振方式下升力特性变化曲线,分别研究了振动频率、雷诺数、自然转捩和固定转捩、数据采集方式等参数对机翼气动特性的影响。试验结果表明：半模机翼振动的气动特性与二元翼型和数值计算情况有所不同,其影响因不同翼型构型、不同采集方式和翼面不同流动模式等会产生不同效果。并对产生原因进行初步探讨。     

二维机翼二自由度跨音速颤振分析

 <正> On the assumption that airfoils take small amplitude ha-rmonic oscillating and middle reduced frequency, the two-dimensional, inviscid, small-disturbance nonlinear unsteady transonic flow equation can be split into two parts. One is the nonlinear steady small-disturbance equation, the other is the time-linearized unsteady small-disturbance equation. The steady equation is solved by Carlson code. The unsteady equation is solved by integral method. The standard v-g method is used to solve the flutter eigenvalue equations. A transonic flutter analysis is performed for a NACA 64A006 airfoil with pitching and plunging degrees of freedom.The Mach numbers considered are 0.7, 0.8, and 0.85. The aerodynamic coefficients are obtained by the above method. For each Mach number, the flutter speed and the corresponding flutter reduced frequencies are achieved by varying the airfoil-airmass ratio, plunge-pitch frequency ratio. The results are compared with other computational method and good agreement has been observe     

篦齿封严结构气弹稳定性数值分析

针对航空发动机中篦齿封严装置的气弹稳定性问题,用所发展的气弹稳定性分析的能量法,考查了两种不同的篦齿封严结构模型的颤振特性,对比分析了在不同节径振动下的气动功分布和模态气动阻尼比。结果表明:能量法可以用于预测篦齿封严结构的气弹稳定性,通过数值计算给出了改进前后篦齿封严结构低阶振动模态的气动功和模态气动阻尼比,改进后的篦齿结构在2节径振型下不再发生颤振,说明改进后的结构比改进前的气弹稳定性更好。      

T型尾翼风洞颤振实验保护装置绕流特性分析

采用SST两方程湍流模型,通过求解非定常Navier-Stokes（N-S）方程,对T型尾翼风洞实验流场进行了模拟,分析了保护装置对T型尾翼风洞实验流场的影响,研究了保护装置几何外形和保护装置基座后移距离对流场影响。通过对平尾气动力的分析以及对非定常流场的对比,可以得出：采用NACA系列翼型对基座进行气动整流后,基座两侧局部超声速区显著减小,局部激波减弱甚至消失,流场品质得到改善。且采用NACA0010翼型对基座修形后的结果最理想。随着保护装置基座后移距离的增加,平尾气动力均方根值和波动值先是急剧减小,达到0.85倍平均气动弦长后开始有所增大,在2.45～4.05倍平均气动弦长范围基本不再变化,稳定到单独T型尾翼模型相应系数1倍左右。此结论对T型尾翼风洞颤振实验保护装置设计具有一定的指导意义。     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

A reduced order model for coupled mode cascade flutter analysis

A Reduced Order Model (ROM) based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper. The unsteady aerodynamic model is established by a system identification technique combined with a set of Aerodynamic Influence Coefficients (AIC). Subsequently, the aerodynamic model is encoded into the state space and then coupled with the structural dynamic equations, resulting in a ROM of the cascade aeroelasticity. The cascade flutter can be determined by solving the eigenvalues of the ROM. Bending-torsional coupled mode flutter analysis for the Standard Configuration Eleven (SC11) cascade is used to validate the proposed method.