基于能量法的跨声速风扇叶片气弹稳定性研究

基于能量法原理,采用弱耦合的方法对跨声速风扇叶片进行气弹稳定性分析;采用3维线性插值算法编程实现CSD/CFD数据交换。利用有限元法计算叶片的模态振型和固有振动频率;应用FLUENT动网格技术,对NASA R 67风扇叶片在不同模态振动下的非定常流场进行数值模拟,给出在前3阶模态振动下叶片表面的周期累积气动功和气动阻尼,并探讨了激波对跨声速风扇叶片气弹稳定性的影响。     

转子叶片气弹稳定性与强迫响应分析

针对由叶轮机械转子叶片流致振动引起的叶片高循环疲劳失效的突出问题,采用数值模拟方法对叶片流致振动特性进 行预测分析。常见的流致振动问题包括气弹稳定性和强迫响应,分别以压气机和高压涡轮转子叶片为例给出了气弹稳定性和强 迫振动响应的分析方法及其流程。对于压气机叶片,采用能量法和特征值法进行了气弹稳定性计算。结果表明：气弹稳定性结果 相互吻合,在高阶模态下的气动阻尼比在低阶模态下的更小,根据振型特征确定危险点位于叶尖近前缘和近尾缘部位,在该部位 可能导致叶片掉角;对于高压涡轮叶片,分析了其强迫振动响应特征。结果表明：叶片非定常气动激励主要来源于上游导叶,在设 计状态下其频率距离转子叶片第5阶模态频率较近,并对危险模态阶次的共振应力进行分析,在0.14%阻尼下振动应力最大可达 134 MPa。     

前后掠风扇叶片颤振特性对比分析

采用基于能量法的流固耦合数值预测方法对比研究了模态、叶间相位角对前后掠风扇转子叶片气动弹性稳定性的影响.结果表明:对于两种掠形叶片,在所研究的前3阶振动模态下,前掠叶片对应于弯扭耦合振型的气动模态阻尼比最大,其数值为0.801%,而后掠叶片对应于弯扭耦合振型的气动模态阻尼比最小,其数值为0.248%;叶间相位角对两种掠形叶片气动模态阻尼比都有显著影响,在1阶振动模态相同叶间相位角(节径)下,前掠叶片前行波对应的气动模态阻尼比小于逆行波对应的气动模态阻尼比,而后掠叶片与此相反;1阶振动模态所有叶间相位角下,前掠叶片比后掠叶片气弹稳定性更好.      

压气机转子叶片颤振边界的预测方法

通过发展的基于能量法的颤振数值预测方法得到了压气机转子叶片的颤振边界.将计算结构动力学(CSD)分析得到的叶片表面节点位移插值到耦合面的流体网格点上,并将设计的多层动网格技术应用于计算流体动力学(CFD)方法,实现叶栅振荡作用下的非定常分析,得到叶片表面的非定常气动功以及模态气动阻尼比.以某第一级压气机转子叶片为例,对...     

一种基于N—S方程的CFD／CSD耦合计算方法

对于气动弹性问题,在跨音速和大变形情况下,模态叠加法失效,为此需要探讨全场求解耦合方程的新方法。利用耦合计算流体动力学和计算结构动力学（CFD／CSD）方法,数值模拟了机翼的气弹响应问题。气动力的求解采用三维Navier—Stokes方程的有限体积法,实现了一种多块结构网格的动网格方法生成气动网格,采用有限元方法计算应...     

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

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

基于振动时滞法的非零叶间相位角叶片颤振

发展了一种利用叶片延迟振动设置叶间相位角的振动时滞法和多通道叶片非同相振动的流固耦合颤振分析模型。模型通道数选取相邻节径线之间通道数的两倍,在循环扇区的不同通道中,令叶片的各阶振动模态位移滞后于前一叶片,结合基于虚拟弹性体的快速动网格算法实现流场及叶片网格的高效更新。针对Rotor 37多通道模型,研究了不同叶间相位角对叶片气弹稳定性及通道流场特性的影响。结果表明:多通道方法与全环叶片颤振分析的计算结果基本一致,而18节径振动下多通道方法的计算时间是全环分析的1486%;节径振动形式对气动阻尼有显著影响,且在2节径时发生气弹失稳;叶间相位角引起流道内激波位置和强度变化和非定常激波脉动异相冲击,是影响颤振的主要原因。     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。     

叶栅流场中的气动阻尼研究

从气动阻尼产生的机理出发,建立合理的气动阻尼模型;利用动网格技术,结合等效模态气动阻尼比的计算式,以Rotor37跨声叶栅为例,用弱耦合方法数值计算得到一弯、一扭与弯扭耦合模态的气动阻尼比,并考虑流场压力水平及叶片振幅对气动阻尼的影响.计算结果表明,在保证进出口压比不变的情况下,气动阻尼与压力成线性关系;在小振幅下,气动阻尼基本与叶片振幅无关.      

压气机叶片流固耦合数值计算

采用流固耦合方法对压气机单排跨声叶片进行气动弹性数值模拟.首先对压气机流场进行三维定常求解,并与试验结果进行比对.在定常流场分析的基础上,进行流固耦合计算,采用非定常方法分别求解流体域Navier-Stokes(N-S)方程和固体域结构动力学方程.耦合过程中,流体域和固体域在物理边界上进行能量传递.计算结果表明:所采用的流固耦合数值方法用于判断叶轮机械叶片气弹稳定性是可行的,可以得到叶片的瞬态响应,从而判断叶片是否发生颤振.      

基于高拟真度模型的高超声速飞行器静气动弹性优化

为解决基于计算流体动力学(CFD)和计算结构动力学(CSD)高拟真度模型的静气动弹性优化过程中模型更新自动化困难、求解速度慢的难点,提出了基于几何主模型技术的CFD/CSD一体化参数化建模方法和基于网格单元修正的常体积四面体(CVT)数据交换技术,并引入响应面模型来降低气动弹性优化求解难度.通过AGARD 445.6机翼静气动弹性分析对上述方法的可行性和有效性进行了验证;并以某高超声速飞行器为例,采用基于二次响应面的多目标优化算法进行了CFD/CSD气动弹性分析与优化,优化后飞行器升阻比增加16%,结构质量减少9%,且响应面模型精度拟合误差不超过1.5%.计算结果表明:所发展的CFD/CSD一体化参数化建模与优化方法能够有效地解决高拟真度模型的静气动弹性优化问题.      

直升机旋翼桨-涡干扰状态非定常气弹载荷高精度预估

为准确计算直升机旋翼在复杂的桨-涡干扰(BVI)状态下的气弹载荷,在刚性旋翼计算流体力学方法中引入桨叶弹性变形的影响,建立了一套适合于弹性旋翼BVI状态气动特性分析的计算流体力学/计算结构力学(CFD/CSD)耦合方法.CFD模块对Reynolds averaged Navier-Stokes(RANS)/Euler方程进行求解,并采用双时间法推进和Baldwin-Lomax(B-L)湍流模型.CSD模块采用中等变形梁假设的有限元模型,通过Newmark-Beta方法求解桨叶运动方程.通过代数变换方法进行桨叶网格变形,并建立一个适于流场/结构信息交换的CFD/CSD耦合方法.在分别验证CFD和CSD模块的有效性的基础上,开展UH-60A直升机旋翼的BVI 状态载荷分析,并与飞行测试数据进行了对比.计算结果表明:相比于旋翼综合分析中的升力线理论和刚性旋翼CFD方法,耦合的CFD/CSD方法可以更准确地预测BVI状态气弹载荷,并有效地模拟桨叶前行侧方位角和后行侧方位角附近的BVI现象,对BVI导致的升力波动幅值和相位的计算结果均与试验值吻合良好.      

主动流动控制技术的初步数值研究

利用亚跨超CFD软件平台2.0版本(trip2.0)对几种常见的流动主动控制技术进行了数值模拟,对包括进气道和襟翼等几个算例的流场进行了定常和非定常求解,获得了合理的计算结果,为今后的进一步深入研究打下了良好基础.     

考虑喷流效应的机载导弹发射及气动干扰数值模拟

研究比较了机载导弹发射时是否考虑导弹尾喷流情况下导弹运动特性及机-弹气动干扰影响.通过动态结构嵌套网格技术实现了导弹运动的模拟,结合刚体六自由度运动方程求解三维非定常Euler方程得到流场信息.导弹推力由发动机燃烧室喷流计算得到;当不考虑喷流效应时,推力通过直接在导弹尾部给定非定常作用力实现.应用此方法模拟了基于类"全...     

悬停状态共轴刚性双旋翼非定常流动干扰机理

基于运动嵌套网格方法,建立了一套适合于悬停状态下共轴刚性双旋翼非定常干扰流场分析的计算流体力学(CFD)方法。首先,基于高效的运动嵌套网格技术,采用积分形式的可压雷诺平均Navier-Stokes(RANS)方程作为双旋翼非定常流场求解控制方程,湍流模型选用Baldwin-Lomax模型,时间推进采用双时间方法。在CFD方法的验证基础之上,对干扰过程中的桨尖涡涡核位置及强度演变规律进行了细致分析,揭示了共轴双旋翼非定常干扰流场中上、下旋翼桨尖涡与双旋翼桨叶之间的贴近干扰、碰撞现象,以及上、下旋翼桨尖涡之间的相互干扰机理。然后,进一步研究了不同总距角下的共轴旋翼系统中上、下旋翼的非定常气动特性以及影响规律。计算结果表明:上旋翼桨叶的桨尖涡会直接与下旋翼桨叶发生碰撞,导致下旋翼桨叶拉力损失;上旋翼桨叶的桨尖涡和下旋翼桨叶的桨尖涡相互干扰,改变了桨尖涡的强度和轨迹;上、下旋翼桨叶相互靠近时,上、下旋翼桨叶的拉力均会上升,之后相互远离时上、下旋翼桨叶拉力均会先下降再上升。     

基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

基于动网格技术的飞机舱弹分离过程 CFD/ RBD耦合数值仿真

飞机舱弹分离瞬间特性严重影响导弹/炸弹的打击精度。采用计算流体动力学(CFD)和刚体动力学(RBD)模型耦合求解的数值仿真方法,应用网格弹性光顺和局部重构相结合的动网格技术,有效地避免了刚体运动而引起的网格畸变,仿真计算了炸弹舱弹分离后俯仰角度变化过程 .仿真结果表明,舱弹分离时向炸弹施加 1个合适的抛 射力矩,有助于炸弹在下降过程中快速地调整其打击姿态。     

Active aeroelastic control over a four control surface wing model

This paper describes the procedure implemented to design, develop and test an aeroelastic control system installed on a forward swept wing of the aeroelastic demonstrator X-DIA. A control method directly based on Nissim aerodynamic energy concept has been chosen. Two different modeling techniques have been adopted for the calculation of generalized aerodynamic forces, such as doublet lattice method and computational fluid dynamics and the obtained results are finally compared. The latter approach, applied to better estimate the control surfaces effectiveness, requires the capability to correctly model the control surface rotation and the grid deformation, usually addressed as non-trivial problems in CFD based aeroelastic analysis. A genetic algorithm optimization technique has been adopted to state and refine all the control gains.     

多块多重网格法及其跨声速转子内流并行模拟

为提高压气机内部复杂流动计算的收敛速度,发展了一种多块多重网格上的可压缩流动并行计算流体动力学(CFD)方法,数值模拟了跨声速转子NASA Rotor 35内部流动及气动性能,重点分析了对接分区方式和内界面通信模式对并行性能的影响。结果表明:对接分区方法和多块多重网格内界面数据处理方法能够保证串、并行CFD方法结果基本一致,且与实验数据吻合良好;并行精度基本不受分区数目和通信模式的影响,并行效率随着相对通信负载的增大而降低。所发展的并行计算方法对跨声速转子气动分析和设计具有一定的可靠性,为开发大型流体机械复杂流动的基础并行算法及软件提供了参考价值。