直升机桨叶刚柔耦合特性及计算方法分析

为了研究桨叶刚柔耦合特性,采用中等变形梁模型处理桨叶弹性变形,分别采用小转角和有限转角两种方法处理桨叶绕铰的刚性运动。以桨叶扬起下坠碰撞问题为研究对象。计算分析表明:(1)弹性桨叶存在较明显的刚性运动和弹性变形耦合;(2)挥舞角较小时,两种处理方法计算结果一致,均与试验数据吻合较好;(3)挥舞角较大时,刚柔耦合增强,两种处理方法的计算值差别较为明显;(4)桨叶弯曲刚度对桨叶刚柔耦合特性有较明显影响,刚度越小,耦合越强;(5)桨叶如果比较柔软,小转角处理方法计算值相对有限转角处理方法,相位滞后、幅值增加均较为明显,计算桨叶动响应时需引起注意。      

基于时间步进自由尾迹的旋翼总距突增气动响应分析

自由尾迹的桨叶气动模型中,考虑了由于桨叶预扭导致尾随涡不在同一个平面对控制点诱导速度的影响,翼型气动模型采用Leishman-Beddoes非定常气动模型,并推导得到操纵下的桨叶挥舞运动方程。采用二阶精度时间步进格式,耦合桨叶气动模型、挥舞运动方程,构建时间步进自由尾迹算法。而后计算旋翼尾迹在操纵条件下的变化历程,得到旋翼的气动响应。通过计算结果与实验数据的对比,验证方法的可行性,并分析操纵速率、预扭对气动响应的影响。     

基于Broyden法的旋翼多体系统气动弹性分析

建立了旋翼多体系统气动弹性模型并给出了一种适合于该模型响应计算的数值计算方法。采用柔性多体系统动力学方法建立旋翼气动弹性模型,利用驱动约束显著简化约束方程形式,集成大变形桨叶模型,准确考虑变形的非线性,适合于对采用柔性结构的先进旋翼进行气动弹性分析。基于Broyden法改进隐式积分法积分一步中非线性方程的求解,避免求取切线矩阵和矩阵求逆运算,保持隐式积分法具有较好稳定性的同时提高计算效率,解决了旋翼多体系统气动弹性力学方程隐式表达且具有较强非线性和较高刚性比造成的响应计算困难。通过模型旋翼桨叶响应计算验证了结构模型与气动弹性响应求解方法。采用建立的气动弹性模型计算悬停和前飞状态旋翼气动弹性稳定性,与试验结果对比验证了模型的正确性。研究了不同的稳定性计算方法、桨叶结构模型和入流模型等对悬停和前飞稳定性计算的影响,结果表明本文所采用的结构、气动模型及气动弹性稳定性计算方法提高了气动弹性稳定性分析精度。     

基于气弹耦合的旋翼桨涡干扰噪声预测方法

建立了基于气动/弹性耦合的旋翼桨涡干扰（BVI）气动和噪声分析方法。气动模型包括修正Beddoes尾迹模型和CFD模型,噪声计算采用基于声学类比法推导出的FW H（Ffowcs Williams Hawkings）方程,弹性桨叶动力学建模采用有限元方法。应用所建立的方法,对刚性的OLS（operational load survey）旋翼桨涡干扰状态的气动和噪声特性进行了计算,对比了两种气动模型在研究桨涡干扰问题的有效性;以弹性的HART Ⅱ旋翼为研究对象,分析了桨叶弹性、时间步长对桨涡干扰气动载荷和噪声的影响。结果表明:进行桨涡干扰计算时所采用的时间步长不宜超过2°。CFD方法由于固有的数值耗散,计算出的OLS旋翼噪声声压峰值仅为试验值的60%,而修正Beddoes尾迹模型能够避免数值耗散,且具有高效率的优势。考虑桨叶气动弹性能够提高旋翼桨涡干扰噪声的预测精度。      

基于PC2B格式自由尾迹方法的旋翼响应特性分析

旋翼自由尾迹模型的建模和桨叶动态响应计算的问题是直升机空气动力学领域中富有关键意义的课题。本文结合桨叶气动模型、旋翼涡尾迹模型、桨叶挥舞动力学模型、"PC2B"尾迹求解算法和旋翼配平模型,建立了一个时间精确的旋翼自由尾迹和桨叶动态气动响应的分析方法,并针对旋翼总距突增时的旋翼载荷和桨叶动态响应进行了计算和分析,得出了一些有意义的结果。     

倾转状态的旋翼尾迹形状和挥舞响应分析

本文推导了旋翼在倾转时的桨叶挥舞运动方程,给出了桨叶升力面模型和尾迹时间推进预测—校正算法,并结合涡核模型和旋翼配平分析模型,建立了一个时间推进旋翼自由涡计算方法。通过计算的旋翼诱导速度与实验值的对比,验证了本文方法的有效性。然后利用该方法,计算了旋翼在倾转时的尾迹几何形状畸变和桨叶挥舞响应特性。     

基于变步长CFD/RBD方法的旋转弹轨迹仿真

采用变步长方法实现非定常Navier-Stokes方程和刚体六自由度运动方程的耦合求解,并对美国ARL旋转弹试验模型的飞行轨迹进行了仿真。旋转弹自由飞行过程中受到空气阻力作用,转速下降,气动力变化趋缓。因此轨迹仿真时先选取较小的时间步长,随着转速的下降逐渐增大步长,这样可以减少计算时间。刚体六自由度运动方程的数值求解采用变步长Adams预估校正法。将定步长和变步长两种方法应用于旋转弹实验模型轨迹仿真,并将计算得到的旋转弹空间位置、姿态角与试验结果进行了对比,结果显示:方法能较精确地仿真旋转弹飞行轨迹。在相同初始步长条件下,变步长方法较定步长方法计算时间减半;在相同计算时间条件下,变步长方法较定步长方法精度更高。     

大型航天器再入解体气动力热特性模拟的直接模拟蒙特卡洛方法研究

为模拟大型航天器离轨再入近连续过渡流区高超声速气动力/热绕流特征,构建了基于直接模拟蒙特卡洛法碰撞限制器技术的混合方法,发展了基于密度梯度的动态自适应混合网格处理技术与变时间步长计算方案。利用当地流动梯度的克努森数作为判断连续流失效的参数,将流场划分为不同区域,在连续流区采用碰撞限制器以及大网格尺度和大时间步长,在流场的大梯度区域——包括激波和壁面边界层区域——采用基于当地密度梯度的动态自适应碰撞网格和取样网格处理技术。为保证整个流场范围每个碰撞网格内的模拟粒子数分布更加均匀,采用变时间步长计算方案,并固定当地时间步长与粒子权重的比值,避免了因分子穿越网格界面产生的复制或消失。通过计算类天宫飞行器低密度风洞试验状态的气动力系数,并与试验数据对比,验证了上述算法的高精度模拟能力与可靠性。同时模拟分析了带太阳电池帆板的类天宫飞行器再入85 km高超声速复杂气动力热,及头部对接台与板舱非规则物形绕流所致激波/边界层干扰、流动分离与强气动力热致太阳电池帆板毁坏发生首次解体机制。     

一种无人机螺旋桨安装特性预测方法

针对用风洞试验选配无人机推进式螺旋桨费用较高的问题,提出用一根螺旋桨的装机风洞试验数据来推算几何参数相同而桨叶安装角不同的一组螺旋桨的装机特性的方法。方法由风洞试验数据计算出螺旋桨桨叶的当量拉力和扭矩系数,通过调整桨叶攻角推算安装角不同的螺旋桨的推力和扭矩特性。计算结果与试验结果的对比表明方法在桨叶安装角变化较小时是准确的。     

桨叶结冰对尾桨气动性能的影响

基于某直升机尾桨桨叶地面结冰试验数据,建立了结冰桨叶的气动计算模型。对NACA0012翼型的气动特性计算表明,翼型升力系数和阻力系数的计算结果与试验数据吻合良好。然后根据所建模型,利用Fluent软件分别计算了某尾桨桨叶翼型结冰前后剖面的气动特性,发现结冰使桨叶翼型升力系数降低,阻力系数增大。最后采用动量一叶素理论结合...     

铰接式旋翼舰面瞬态气弹响应及参数研究

为了研究舰面铰接式旋翼桨叶桨尖与机体相碰的原因及参数影响,本文采用中等变形梁模型处理桨叶弹性变形,用有限转角模拟桨叶绕铰刚性运动和非线性准定常气动模型处理气动力。通过与国外计算及试验数据的相关性分析,验证了本文建模及计算方法的正确性。参数分析表明:1)挥舞限动块脱开前,桨叶弹性变形提供的桨尖挥舞位移占主要部分,脱开后,桨叶刚性运动提供的桨尖挥舞位移占主要部分;2)舰面流场采用梯形流时,计算的桨尖挥舞位移最大;3)挥舞限动角大小及离心式限动块释放时间对桨尖最大挥舞位移影响不明显;4)直接增加桨叶挥舞刚度可显著地减小桨尖向下最大挥舞位移;5)随总距增加,桨尖挥舞挠度减小不明显,这是通过增加总距来减小桨叶挥舞位移不可行的另一原因。     

Aeroelastic analysis of a shipboard helicopter rotor with ship motions during engagement and disengagement operations

An aeroelastic simulation of a shipboard helicopter rotor with ship motions during engagement and disengagement operations is investigated to explore the coupled dynamic behavior between the rotor and the ship. A finite element analysis based on a moderate deflection beam model is employed to capture the flap, lag and torsion deflections of the rotor blade. The ship is treated as a six-degree-of-freedom rigid body. By using the Hamilton?s principle, system equations of motion are derived based on the generalized force formulation. The responses agree well with the test data of the rotor blade droop stop impact and the transient aeroelastic response of the shipboard teetering model rotor. Parametric investigations illustrate that the ship pitch motion has significant influence on the maximum negative displacement of the blade tip. Additional over 25% increase of the tip deflection can be introduced by the ship pitch motion. The aerodynamic and inertial couplings between the ship motion and the rotor have significantly nonlinear influence on the transient aeroelastic response. Both terms should be taken into account in the coupled helicopter–ship dynamics model.     

航空发动机风扇叶片与机匣刮蹭分析及结构设计

针对异常载荷下,航空发动机宽弦风扇叶片的叶尖与机匣刮蹭变形及损伤特征缺乏数据支持,而传统理论计算方法存在较大的误差问题,建立了宽弦风扇叶片叶尖刮蹭显式动力学分析模型,采用宽弦风扇叶片与机匣刮蹭试验数据,对分析模型的计算精度进行了验证。基于分析模型进行了仿真参数的敏感度分析,得到了叶片与机匣刮蹭后叶片变形及机匣损伤规律。研究结果表明:叶尖伸长量对转子转速非常敏感,叶尖径向伸长量增加速率远大于转速增加值,因此在叶片设计中应考虑到风扇叶片极限转速下叶尖伸长量。同时需要选取合理的扭转角度以满足叶片安全性和气动性能的要求。在风扇机匣包容区设计中应主动考虑异常载荷的影响,增大安全性设计域度;设计合理的耐磨层材料参数,减小风扇叶片对其冲击损伤。采用该方法可以提高叶尖间隙控制精度,减小刮蹭对叶片和机匣造成的损伤。      

基于CFD/CSD耦合的叶轮机叶片失速颤振计算

为了研究叶轮机叶片的失速颤振特性,发展了一种计算流体力学与计算结构力学(CFD/CSD)时域耦合方法。该方法通过每一物理时刻CFD和CSD的循环迭代实现了耦合计算。在CFD分析中,采用鲁棒性较好的空间离散格式AUSM+-UP,并基于延迟脱体涡模型(DDES)模拟了带分离流动。在结构分析中,通过模态法构建了旋转叶片动力学方程并运用杂交多步方法进行求解。以孤立转子Rotor37为例,计算了不同工况下流场总体与细节参数,与实验结果的对比验证了CFD算法的精度。对某转子叶片进行了颤振特性研究,计算所得的广义位移时间响应曲线表明该叶片在近失速工况下会发生失速颤振,其表现形式为一阶弯曲模态发散且各阶模态之间不耦合。分析表明,流场不稳定和非定常效应是引起失速颤振的关键因素,同时折合频率的降低也会导致原本气动弹性稳定的叶片发生失速颤振。     

一种新的单机场地面等待问题优化算法

对单机场地面等待问题提出了一种新的启发式动态优化算法。该方法通过对时间分段，添加飞机是否起飞标志量，分步计算地面延迟和空中延迟，递推求解，从而达到缩小每次求解规模的目的。采用实际数据，仿真验证了所提出的模型和算法的有效性。还给出了一种简便的近似解算法。     

基于粗糙集理论的终端区飞机排序综合评判

针对日益增长的空中交通流量所带来的严重航班延误,提出了一种全新的终端区（TMA）飞机排序方法。在介绍粗糙集综合评判方法的评判过程基础上,建立了TMA飞机排序问题评判指标体系,进行了飞机排序综合评判过程的算例仿真。仿真结果表明,方法简单可行,能够客观得到各评判因素权重,并且与先到先服务（FCFS）排序方法比较,影响范围减少50%,累计延迟时间减少23%,说明方法排序结果优于FCFS排序结果。     

基于响应面法的鸟撞风扇叶片损伤预测

有限元模拟鸟撞风扇叶片损伤成本高,为解决工程问题,采用经典叶栅鸟撞切割模型建立了鸟撞风扇叶片动载荷数学模型,结合鸟撞部件试验结果,以拟合技术明确风扇叶片损伤程度与最大关键动载荷计算值间的函数关系,形成叶片损伤预测响应面,实现对鸟撞风扇叶片损伤的快速预测,并建立基于响应面法的鸟撞风扇叶片损伤预测工作流程。结合涡扇发动机吞鸟试验技术要求、风扇结构设计特征及已开展的鸟撞部件试验结果,建立叶片损伤预测响应面,初步识别2种鸟撞方案的径向弯曲、弦向弯曲,并计算撕裂范围分别不超过0.3867和0.3941,撕裂与弦向弯曲相关性显著,呈抛物线变化趋势。结果表明：预测的损伤在可接受的安全性水平范围内,预测方法能够识别损伤范围及趋势,可为后续鸟撞有限元模拟、试验策划、安全性分析、风扇叶片抗鸟撞设计等工作提供量化的技术支持。     

冲击载荷作用下的锥壁失效仿真及试验验证

根据航空发动机支承锥壁结构受力特点,对风扇叶片飞失冲击载荷作用下的锥壁失效破坏机理进行了研究.利用显式动力学有限元仿真方法,对冲击载荷作用下的锥壁结构动态失效过程进行了瞬态分析.开展了对锥壁的落锤冲击试验,并与分析结果进行了对比验证.试验和分析结果表明:冲击载荷作用下锥壁减薄处破坏为剪切失效破坏.采用的显式动力学有限元...     

直升机旋翼载荷中的非线性强迫振动方程迭代算法

介绍用迭代方法求解直升机桨叶气动载荷中的非线性强迫振动方程。计算时,先假定诱导下洗已确定,求出桨叶的初始载荷;再通过非线性强迫振动方程求出弯－弯－扭弹性变形。以某一直升机为例,求出了桨叶的各阶谐波载荷和弹性变形。计算结果与实验数据吻合。