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

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

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

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

小展弦比带掠风扇转子叶片的耦合颤振

从结构力学出发,对某发动机整体叶盘进行了盘片耦合振动特性分析,通过分析小展弦比带掠风扇转子叶片的共振转速、耦合颤振和失速颤振,发现了该小展弦比带掠风扇转子叶片存在2阶与3阶动频与静频振型发生变化的特殊现象,揭示了该小展弦比带掠风扇转子叶片为另一种耦合颤振的问题.通过对该小展弦比带掠风扇转子叶片的结构调整,找出了叶片前掠是该叶片发生此类特殊现象的重要影响因素,适当压缩该叶片的前掠后,从设计上避免了发生颤振的可能,提出叶片各阶频率的分离度需满足10％的设计要求.      

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

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

基于CFD技术的叶片颤振分析

采用CFD技术,引入等叶片间相位角假设,运用非耦合的能量法,研究了典型叶轮机叶片的颤振特性。采用有限体积法求解NS方程,运用SST湍流模型,通过恒定叶间相位角简化叶轮机流场,计算了NASARotor67和STCF4（the Fourth Standard Configuration）的振动叶片的流场特性。计算表明：转子67采用碳钢材料时,在效率峰值的工况下,叶片不会发生颤振;而STCF4在一阶弯曲振动中会出现不稳定;叶间相位角（IBPA）对振动叶片的气动弹性稳定性有明显的影响。     

基于功能原理的颤振模态参与度分析方法

颤振边界与颤振耦合机理对飞行器设计与颤振试验设计有着重要意义。模态坐标下的颤振计算通过各阶广义坐标的变化特性来确定颤振边界，并对颤振耦合机理进行分析。然而这常常依赖设计人员的工程经验，难以保证判别标准量化统一。基于功能原理提出了一种新的颤振模态参与度分析方法，通过广义坐标下广义气动力做功的能量累积，实现颤振模态参与度分析。分别采用频域方法与CFD/CSD方法对AGARD445.6标模进行颤振计算，验证了所提模态参与度分析方法的正确性。随后针对双体飞机颤振计算中发现的“漂移频率”现象，使用所提方法进行了解释，凸显了方法的优势。综合表明，所提模态参与度分析方法较好地反映了颤振耦合机理，具有指标正确可靠、结果归一化强、物理意义明确、适用于复杂结构复杂模态的特点。     

基于CFD/CSD技术的压气机叶片流固耦合及颤振分析

针对叶片故障的原因颤振问题,基于计算流体动力学/计算结构力学(CFD/CSD)耦合算法,研究了压气机典型叶片的流固耦合(FSI)问题及颤振特性.使用有限体积法求解非定常三维Navier-Stokes方程和有限元法求解三维结构体模型,在两个求解器间由载荷转换、网格变形传递和同步化方法完成数据交换及求解.以高性能风扇NAS...     

上下游叶排对压气机转子叶片气动弹性稳定性的影响

使用自行开发的流固耦合程序集,分别考虑压气机上下游叶排的影响对转子叶片气动弹性稳定性进行分析.使用全环气动弹性模型通过叶片在非定常流场中振幅随时间的变化历程计算转子叶片不同振型下的气动阻尼,分析了不同的叶排轴向间距下尾迹和势干扰对气动阻尼的影响规律.通过将考虑转静干涉效应的气动阻尼与单转子的结果作比较,总结了转子-静子结构和导叶-转子结构的干涉作用对转子叶片颤振特性的影响规律.结果表明:对于1阶弯曲模态,转子叶片气动弹性稳定性不随轴向间距单调变化;上游尾迹干涉作用和下游静子势干扰的增强会加剧1阶扭转模态失稳,但是却促进1阶弯扭耦合模态气动弹性稳定.      

叶片排间耦合对跨声速压气机转子颤振特性影响研究

高推重比、强紧凑性的设计要求迫使压气机相邻叶片排之间的距离进一步减小，在加剧排间非定常扰动的同时，势必对叶片的气动弹性稳定性产生影响。本文针对1.5级跨声速压气机，基于谐波平衡非定常求解技术，结合单向耦合的能量法，开展排间耦合对1.5级跨声速压气机叶片颤振特性影响研究。通过对比分析非定常/定常混合模型的颤振预测结果可以发现：排间压力波反射对振荡叶排压力脉动特性的影响是造成转子颤振特性差异的主要原因。上、下游静叶对转子气动弹性稳定性恶化程度并不相同，且不符合线性叠加原理。对于近失速工况，非定常压力波反射导致转子叶片气动阻尼系数减小75%，定常/非定常混合预测方法可能导致颤振特性预估过于乐观。     

涡轮机械叶片的流固耦合数值计算方法

发展了一种流固弱耦合数值方法来判断涡轮机械叶片气弹稳定性.通过一种数据交换方法将计算结构动力学(CSD)的节点振动位移施加到计算流体力学(CFD)的叶片表面网格点上,CFD分析采用多层动网格技术实时更新可动域的网格点坐标,并通过有限体积法求解了用k-ε湍流模型封闭的雷诺平均Navier-Stokes(RANS)方程.以NASA 67转子叶片为例,在设计状态下,对比了叶片前4阶模态下一个振动周期内的气动功及模态气动阻尼比,初步分析了叶片的气弹稳定性.计算结果表明:所发展的流固弱耦合数值方法用于判断涡轮机械叶片气弹稳定性是可行的.      

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

为准确计算直升机旋翼在复杂的桨-涡干扰(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/CSD的非线性气动弹性分析方法

提出了一种基于计算流体力学/计算结构动力学(CFD/CSD)的非线性气动弹性分析方法,并应用于切尖三角翼的非线性颤振和极限环振荡(LCO)研究。该方法将非线性有限元(FEM)和CFD计算相结合,并辅以高精度的界面插值,能够分析结构和气动非线性共存的气动弹性问题。结构部分以四边形平板壳元为基础,采用更新的拉格朗日(UL)方法分析结构大变形引起的几何非线性问题。气动部分以Navier-Stokes(N-S)方程作为控制方程,采用CFD方法计算跨声速气动力。机翼的非线性颤振计算表明了方法的有效性。最后应用该方法研究了切尖三角翼的LCO现象,其计算精度明显优于已有结果。     

轴流压气机进口小总压径向畸变对失速颤振的影响

对进口流场畸变失速颤振的影响进行了实验研究。实验结果表明:小的叶尖径向畸变能抑制颤振的发作。对此现象用三维变形激盘法进行了计算。计算结果与实验符合。      

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

叶片失速颤振发作机理的探讨

 本文对于现有的叶片失速颤振计算模型进行了一些改进,以期使其更接近于真实物理图画。计算结果表明,就槽道激波以及物面气流分离对于叶片气动弹性稳定性所起的作用来看,在单自由度纯弯振荡和双自由度弯曲主导振荡两种情况之下是具有实质性差别的。在文中对于一个实际转子的失速颤振边界进行了数值预测,并与相应的实验结果进行了对比。     

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

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

基于CFD的静气动弹性优化设计方法

以三维Navier-Stokes方程为控制方程,数值计算弹性机翼跨音速气动力,耦合结构静平衡方程,研究弹性机翼的静气动弹性变形和真实载荷分布,并在此基础上,对机翼的型架外形进行基于静气动弹性的多学科的优化设计(MDO),以及对机翼进行基于总升力不变的飞行姿态确定,以满足弹性机翼在飞行时的品质设计要求.以某后掠机翼为例,设计结果达到预期目标.     

一种CFD/CSD耦合计算方法

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

非线性壁板颤振分析

利用一种改进的计算流体力学与计算结构动力学(CFD/CSD)耦合方法研究了由气动和结构几何非线性引起的壁板颤振问题。在非定常气动力计算中,考虑了通量分裂格式、隐式时间推进方法和几何守恒律;二维和三维壁板的结构几何非线性建模则采用了有限元的协同旋转理论,并利用一种近似能量守恒算法求解结构的非线性响应。流场和结构求解器采用二阶松耦合方法联立求解,并将其应用于壁板在超声速、跨声速和亚声速的颤振计算中。当考虑结构几何非线性和气动非线性时,出现了典型的极限环振荡现象,并对颤振边界和极限环振荡幅度进行了对比分析。