跨声速工况下流体诱发叶片振动研究

对同一跨声速风扇的非定常流场进行了跨声与亚声两种工况下的数值模拟.对两个工况下流场计算得到的转子叶片表面非定常气动激振力进行频谱分析,并对其进行对比分析,从而得到跨声工况下气动激振力特性.应用ANSYS软件的瞬态分析方法进行了叶片的强迫响应分析,对跨声和亚声两种工况下的叶片强迫响应进行对比分析,从而得到跨声工况下叶片强迫响应特性.      

气流激励下叶片振动响应分析方法

为开展气流激励下叶片振动响应分析方法研究,建立了气动激振力预估方法,采用非线性谐波法对叶排进行三维非定常流动分析,获得叶片表面的脉动压力,编制流固转换程序,计算叶片所受的气动激振力。建立了叶片气动阻尼分析方法,基于能量法和弱耦合分析法,对叶片与流场进行流固弱耦合分析,将气动力对运动的叶片所做的气动负功等效为黏滞阻尼力所做的功,求得转子叶片的模态气动阻尼比。建立了叶片在气流激励下的振动响应分析方法,基于气动激振力和叶片模态气动阻尼比,采用模态叠加法分析叶片振动响应。使用该方法,针对发动机中1.5级压气机转、静子叶排模型,计算了转子叶片在真实流场中的气动激振力、前8阶模态气动阻尼比以及在气动激振力与气动阻尼共同作用下转子叶片的振动响应,振动应力达到100 MPa。      

气动非谐对叶片非定常气动激振力影响

为降低转子叶片在前排静子叶片尾流激振情况下的共振应力, 提高转子叶片疲劳寿命, 对静子叶片的气动非谐设计进行研究.对静子叶片的气动非谐设计考虑了两种方案, 一种为非均匀叶栅设计, 另一种为不对称叶栅设计.通过对静子-转子模型流场进行非定常计算, 比较了两种非谐设计思想下转子叶片表面气动激振力频谱成分和幅值大小.分析结果标明:非均匀叶栅非谐设计并不能降低转子叶片表面激振力, 同时会引起频率成分和幅值都较大的倍频成分;不对称叶栅非谐设计虽然会引起与转速同频的较大幅值激振力成分, 但从整体考虑可有效降低叶片振动应力.      

航空发动机进气畸变诱导叶片振动的研究

介绍了进气畸变导致的叶片振动响应的一般分析方法 ,讨论了引起叶片振动响应的周向总压畸变导致的激振力谐波的计算方法 ,并结合某型发动机叶片试验结果和计算结果进行对比 ,对叶片在周向总压进气畸变条件下的振动响应进行了分析和讨论     

叶顶间隙对非定常激振力及下游静子振动影响

压气机叶片叶顶间隙会产生泄漏流,对非定常流动及下游静子会产生一定影响。为研究不同叶顶间隙下非定常气流激振力以及静子叶片振动响应,对不同叶顶间隙下单级压气机的非定常流动进行计算。随后分析非定常结果的流动特性。再利用快速傅里叶变换对静子表面的非定常激振力进行频谱分析,并将按不同频率分阶的载荷加载到静叶上以求解振动响应。非定常结果表明随叶顶间隙增大,效率和压比与之呈线性递减关系,叶片表面时间平均的压力也随之降低。而激振力的频域结果表明,激振力主要频率与叶片通过频率接近,且一阶激振力效应最大。振动响应结果显示随叶顶间隙增加,振动幅值有所降低,但在20%~40%叶高处有非线性现象,位移波谷有轻微下降。且前三阶响应结果接近静叶第3阶,第8阶和第17阶模态振型,频率相差为10%左右。     

弯曲叶片造型对涡轮叶栅作用力影响的非定常数值研究

为了研究弯曲叶片对涡轮叶片非定常气动力的控制作用,采用非定常数值模拟方法对高压涡轮两级的内部流场进行了计算,研究了第一级静叶采用正弯叶片时,汽轮机涡轮叶片非定常气动力和非定常压力场的变化.结果表明,采用弯曲静叶的设计方法,通过改变静叶出口处径向压力梯度和涡量分布,能够有效抑制动静叶片表面气动力的非定常脉动作用.但弯叶片给非定常压力场带来的改善是局部的,其对叶片中间截面改善最为显著.从整体上来看,正弯曲的静叶设计对于减小非定常流场中的压力脉动是有益的.     

大跨度桥梁抖振响应的频域分析

抖振是大跨度桥梁易发生的一种风振形式。本文根据风洞试验结果和随机振动理论,探讨了处于自然大气湍流中大跨度桥梁抖振的计算方法。在该方法中,自激力采用了更适合抖振特点的完全非定常自激气动力形式,抖振力采用了一种可显示表示水平阵风和竖向阵风贡献的非定常模型。     

非均匀栅距对压气机叶片非定常气动力的影响

研究了轴流压气机叶片采用非均匀栅距对叶片非定常力影响的问题,利用高分辨率格式求解准三维雷诺平均N-S方程。非定常流场计算采用了双时间步法与LU-SGS隐式解法耦合的方法。研究了带进口导叶的单级三排叶栅非定常流动,其中动叶栅距不变,而导叶和静叶采用非均匀栅距。对不同结构的计算分析表明,在导叶和静叶采用一定的非均匀栅距组合型式后,在保证压气机级的性能变化不大的同时,动叶表面所受非定常气动力的状况得到改善,这将有助于提高动叶片的疲劳寿命。     

非均匀栅距对轴流压气机流动影响的数值研究

本文进行了非均匀栅距叶片排的轴流式压气机流场数值研究,进口导叶采用了不同的非均匀栅距结构,探讨了不同工作状态下对压气机性能及动叶非定常气动力的影响。数值研究中采用LU-SGS隐式解法与双时间方法相结合,求解圆柱面二维流动的雷诺平均非定常N-S方程。数值研究的结果表明,进口导叶非均匀栅距结构对压气机气动性能影响不大,但对动叶表面受到的非定常气动力有一定影响。采用不同的非均匀栅距进口导流叶片排结构型式,动叶表面的非定常气动力的分布是不同的。不同工作状态下的计算结果显示,即使采用相同的非均匀栅距结构,非定常气动力的分布特性也是不相同的。     

三维机翼大迎角低速绕流及其涡控制的数值模拟

通过数值求解三维非定常可压缩雷诺平均Navier-Stokes方程，研究了三维机翼大迎角低速绕流及其带表面周期性吹吸气的流动控制，并探讨了非定常质量引射作为外激发手段对机翼气动力特性的影响，计算表明，在机翼前缘附近的小区域内施加周期性吹吸气进行流动控制，这类外激发手段对三维机翼的大迎角分离流动的控制是可行的，通过选择适当的控制参数，尤其是合适的周期性激振频率，可以有效地改善气动力特性。     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

直升机旋翼对尾桨非定常气动载荷的影响

悬停和侧滑状态的直升机主旋翼桨尖涡将穿透尾桨桨尖平面,由此导致尾桨非定常气动载荷发生明显变化。为更准确地模拟由主旋翼/尾桨干扰产生的尾桨非定常气动载荷变化,通过在面元压力项中增加由旋翼桨尖涡诱导的时变项,体现旋翼桨尖涡速度和几何时变对桨叶非定常压力的影响,同时采用涡面镜像法修正涡粒子法的黏性项,确保桨叶附近区域旋翼涡量守恒,建立旋翼尾迹对尾桨叶的非定常气动干扰模型,并耦合面元/黏性涡粒子法,构建直升机主旋翼/尾桨干扰下的尾桨非定常气动载荷分析方法。通过计算AH-1G旋翼桨叶非定常气动载荷特性,并与实验测量值、计算流体力学(CFD)计算结果对比,验证本文非定常气动干扰模型的有效性。随后基于NASA ROBIN(Rotor Body Interaction)模型分析悬停、侧风和60°右侧滑状态主旋翼对尾桨非定常气动载荷的影响,分析表明主旋翼尾迹对尾桨非定常气动载荷影响显著。悬停状态的主旋翼/尾桨干扰导致尾桨拉力平均值下降、非定常气动载荷显著增加;左侧风状态,主旋翼/尾桨干扰削弱尾桨"涡环"程度,显著增加尾桨拉力和非定常气动载荷;60°右侧滑状态,主旋翼/尾桨干扰导致尾桨拉力损失最大,且在低速侧滑状态出现尾桨拉力"迅速恢复"现象,尾桨非定常气动载荷幅值迅速增加。     

轴流压气机转静二维非定常薄层N-S数值模拟

采用非定常薄层Navier-Stokes紊流运动方程来描述一跨声轴流压气机中径处的二维转-静非定常流场，速度脉动的相关分析表明转子尾迹对静子的热流流场几乎没有什么影响，而粘性的影响取决于叶片排之间的轴向距离。     

Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor

A robust unsteady rotor flowfield solver CLORNS code is established to predict the complex unsteady aerodynamic characteristics of rotor flowfield. In order to handle the difficult problem about grid generation around rotor with complex aerodynamic shape in this CFD code,a parameterized grid generated method is established, and the moving-embedded grids are constructed by several proposed universal methods. In this work, the unsteady Reynolds-Averaged Navier-Stokes(RANS) equations with Spalart-Allmaras are selected as the governing equations to predict the unsteady flowfield of helicopter rotor. The discretization of convective fluxes is accomplished by employing the second-order central difference scheme, third-order MUSCL-Roe scheme, and fifth-order WENO-Roe scheme. Aimed at simulating the unsteady aerodynamic characteristics of helicopter rotor, the dual-time scheme with implicit LU-SGS scheme is employed to accomplish the temporal discretization. In order to improve the computational efficiency of holecells and donor elements searching of the moving-embedded grid technology, the ‘‘disturbance diffraction method" and ‘‘minimum distance scheme of donor elements method" are established in this work. To improve the computational efficiency, Message Passing Interface(MPI) parallel method based on subdivision of grid, local preconditioning method and Full Approximation Storage(FAS) multi-grid method are combined in this code. By comparison of the numerical results simulated by CLORNS code with test data, it is illustrated that the present code could simulate the aerodynamic loads and aerodynamic noise characteristics of helicopter rotor accurately.     

多片后缘小翼对直升机旋翼桨叶动态失速及桨毂振动载荷的控制

 减缓直升机后行桨叶动态失速发生、降低直升机桨毂振动载荷是提高直升机飞行速度、改进直升机飞行性能的重要途径。本文研究了直升机在高速高载情况下利用多片受控的桨叶后缘小翼对直升机的后行桨叶动态失速和桨毂振动载荷同时进行控制的有效方法。建立了弹性桨叶和后缘刚性小翼的结构动力学模型。桨叶剖面气动载荷采用Leishman-Beddoes 二维非定常动态失速模型计算,后缘小翼剖面气动载荷采用Hariharan-Leishman二维亚声速非定常气动模型计算。采用伽辽金和数值积分相结合的方法求解旋翼系统的气弹响应。建立了有效的多片后缘小翼控制策略和控制方法,分析了3片后缘小翼的运动规律及对后行桨叶动态失速和桨毂振动载荷的控制效果,结果表明利用多片小翼的运动是控制桨叶动态失速和桨毂振动载荷的有效方法。     

舰载倾转旋翼机着舰域耦合流场数值模拟研究

倾转旋翼机是当前旋翼飞行器研究的热点,但有关舰载倾转旋翼机着舰域耦合流场的研究还很少。以两栖攻击舰(LHA)和V-22"鱼鹰"倾转旋翼机为研究对象,基于雷诺平均Navier-Stokes(RANS)方程和SST k-ω湍流模型对舰载倾转旋翼机着舰域耦合流场进行数值模拟研究,并探讨了不同着舰高度时机/舰耦合流场的相互作用。结果表明:倾转旋翼尾流会与舰船脱落涡、甲板舷涡以及舰岛艉涡发生较强的"涡-涡干扰"现象,加大了耦合流场的湍流强度;舰船流场的低频非稳态特征会导致旋翼桨盘气动载荷发生显著的波动,不利于飞行操纵;垂直降落过程中,舰船甲板会形成"前低后高"的压力分布特征,倾转旋翼RMS气动载荷值也会明显增加,降低了着舰安全性,且右旋翼RMS气动载荷值比左旋翼平均大一倍以上,这也表明右旋翼面临着更加严峻的气动环境。     

基于时间精确自由尾迹方法的旋翼非定常 气动特性分析

 给出了一种显式修正的二阶向后差分格式(PC2B)的时间精确自由尾迹算法,用于分析旋翼非定常气动特性。由于旋翼自由尾迹控制方程具有非线性特征,求解过程中易产生数值不稳定性问题,因此本文推导了修正PC2B算法的修正方程,分析了该算法的求解精度和数值稳定性,并通过数值试验对其稳定性和计算效率加以验证及分析。基于该算法,首先对悬停和总距突增时旋翼非定常气动特性进行了计算,通过与试验结果的对比验证了模型的有效性;其次研究了总距操纵时旋翼非定常气动特性的变化规律。研究结果表明,旋翼自由尾迹几何形状的瞬态变化是造成总距操纵条件下旋翼非定常气动特性波动的主要原因。     

高负荷跨声速压气机动静叶干扰下叶片气动负荷研究

通过数值方法模拟了某型高负荷跨声速压气机在工作点的非定常流场，对级内动静叶的气动干扰进行了深入分析。通过气动力和力矩的变化对动静叶干扰下的叶片非定常气动负荷进行了详细的分析，研究了尾迹和势流干扰下叶片气动负荷的变化规律，结合频谱分析进一步认识了气动负荷各周期分量，研究结果有助于对高负荷压气机内动静叶非定常气动干扰的认识。     

旋转机翼飞机旋翼模式前飞状态干扰气动特性

和传统直升机相比,旋转机翼（CRW）飞机在旋翼模式前飞时各部件之间存在更为严重的气动干扰。为了获得旋转机翼/机身/鸭翼/平尾之间的非定常气动干扰规律,基于运动嵌套网格技术,通过求解三维非定常雷诺平均Navier-Stokes（URANS）方程,建立了旋翼前飞流场数值模拟方法。首先对传统直升机旋翼/机身干扰模型进行了计算,验证了方法的可靠性,然后对某旋转机翼飞机全机在旋翼模式前飞状态下的非定常流场进行了数值模拟,并对各个气动部件上的非定常气动力和力矩的变化进行了分析。结果表明：飞机在旋翼模式前飞时,机身部件对旋转机翼的干扰较弱,在经过机身上方时拉力峰值仅略有增加;旋转机翼对鸭翼和垂尾干扰较弱,对机身和平尾干扰较强,随着前飞速度增大,旋转机翼对平尾的干扰会产生较大的升力损失和抬头力矩,需要引起重视。计算结果为该类飞行器的总体综合设计提供了参考。     

螺旋桨滑流对自转旋翼气动特性影响分析

自转旋翼机在前飞时螺旋桨滑流穿过桨盘平面,为研究螺旋桨滑流对自转旋翼的非定常气动干扰,基于RANS(雷诺平均Navier-Stokes)方程,采用运动嵌套网格方法建立了适用于自转旋翼-螺旋桨气动干扰流场的计算分析方法,并对模型进行模拟。分析低速情况下,孤立状态自转旋翼和组合状态自转旋翼非定常气动特性及流场特性,研究不同速度及螺旋桨位置对自转旋翼气动特性的影响。计算结果表明螺旋桨滑流会影响自转旋翼在各方位角的升阻力特性,并引起自转旋翼尾迹在0°方位角附近发生畸变;相同升力下,来流速度越大旋翼后倾角越小,螺旋桨滑流对自转旋翼影响越小;增大螺旋桨与自转旋翼间距可以减弱螺旋桨滑流对自转旋翼的气动干扰。