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

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

横列式双旋翼直升机旋翼对机翼的干扰分析

针对横列式双旋翼直升机旋翼下洗流对机翼的气动干扰影响,建立了旋翼对机翼的干扰计算模型。该模型首先基于万向铰旋翼建立了挥舞运动方程,以得到桨叶挥舞角,然后对桨叶采用非定常Beddoes翼型模型计算气动力和力矩,以考虑桨叶大负扭转带来的失速影响,接着引入动力人流模型获得旋翼处的诱导速度。最后运用经典方法,以XV-15倾转旋翼机为算例,计算了配平状态下旋翼对机翼的向下载荷,并与GTRS模型数据进行了对比,验证了计算模型的合理性。     

基于非惯性系的悬停状态旋翼CFD/CSD耦合气动分析

旨在提高先进旋翼气动特性的分析精度,在旋翼高精度CFD分析中耦合气动弹性效应,取代传统方法中的刚性桨叶假设,并考虑悬停状态旋翼流场准定常的特性,在非惯性坐标系下建立了一套适合于悬停状态旋翼气动特性计算的CFD/CSD耦合分析方法。旋翼气动载荷通过求解三维Navier-Stokes方程求得,空间离散及通量计算采用Jameson中心格式,时间方向则选用五步Runge-Kutta迭代求解,湍流模型采用B-L模型;基于Hamilton原理建立了描述旋翼弹性运动的非线性微分方程,针对旋翼悬停状态的工作特点,采用Raphson迭代方法求解获得旋翼桨叶的弹性变形量。在CFD/CSD耦合计算中,旋翼桨叶交接面载荷及变形信息通过CFD与CSD模块进行传递,同时为提高桨叶弹性变形后贴体网格生成的效率和质量,采用基于网格点坐标转换的网格变形方法。在CFD和CSD程序分别验证基础上,采用建立的旋翼CFD/CSD耦合分析方法计算了先进的UH-60A直升机旋翼的表面压强及气动载荷。计算结果表明,与刚性旋翼CFD模拟结果比较,本文建立的CFD/CSD耦合分析模型可以更准确地预估旋翼气动载荷和性能。     

旋翼翼型对铰链力矩的影响

本文介绍了旋翼翼型对铰链力矩影响的研究。计算并测量了三种具有不同翼型的桨叶在悬停状态下的铰链力矩。表明旋翼翼型对铰链力矩的影响是大的。此外,还通过计算研究了悬停状态下桨叶剖面弦向质心位置对铰链力矩的影响。表明桨叶剖面弦向质心位置越靠后铰链力矩的负值越小,而且随着旋翼拉力的增大影响也增大。     

倾转旋翼机旋翼对机翼气动干扰的建模及分析

建立了一个耦合旋翼自由尾迹模型和机翼面元模型的悬停状态旋翼/机翼气动干扰迭代计算方法,用于分析旋翼对机翼的气动干扰影响.在该方法中,为较好地模拟大桨盘载荷及大扭转桨叶的气动特征,旋翼桨叶采用Weissinger-L升力面模型;为考虑厚度效应及机翼的升力影响,建立了包含源面元和偶极子面元的厚度机翼模型;为正确模拟旋翼桨尖涡与机翼表面间的贴近干扰,采用了一个"分析数值匹配法"的"贴近涡/面干扰模型".应用上述方法,对单独倾转旋翼下洗流分布以及旋翼对机翼的气动干扰影响进行了计算.结果表明,在旋翼下洗流场的干扰下,机翼各剖面都产生向下载荷,但并非简单地随拉力系数的增大而增大;机翼受到的旋翼干扰影响与旋翼下洗流沿桨叶展向变化密切相关.     

不同飞行状态下旋翼与大气湍流的干扰噪声

本文简述了一种直升机旋翼与大气湍流相干扰产生的噪声预估方法。这一方法考虑了桨叶上非定常载荷的弦向及展向非紧致性，桨叶－桨叶载荷之间的相关，通过结合大气液流能量谱模型与快速畸变湍流收缩模型来确定旋翼平面处非均匀、各向异性湍流场的特性，使得这种方法可用于直升机在真实飞行状态下的旋翼与湍流场干扰噪声计算，利用文献中模型旋翼试验的结果，验证了本方法的可行性，并计算了某型直升机在六种飞行状态下的流场畸变情况     

带后缘小翼的旋翼振动载荷计算

建立了考虑弹性桨叶、刚性小翼的旋翼气动弹性分析模型和旋翼载荷计算方法.以广义质量和广义力的形式描述小翼惯性力和气动力对系统的影响,以非定常/动态失速模型计算剖面气动力,结合基于实验数据修正的组合气动模型计算带小翼部分的剖面气动力,集成大变形桨叶模型考虑弹性变形的非线性,以力积分法计算桨叶剖面振动载荷.通过计算分析与实验结果相比较,验证了建立的气动弹性模型和载荷计算方法.结果表明:建立的桨叶结构模型精度很高,气弹模型能够准确预测旋翼的振动载荷,挥舞弯矩平均误差控制在9.1%,使用修正的小翼气动模型能有效提高小翼运动时桨叶振动载荷的计算精度.      

直升机急拉杆机动飞行仿真建模与验证

 针对直升机大机动飞行仿真,建立了一个非线性的飞行动力学模型,考虑了翼型非定常/动态失速、机动飞行引起的动态尾迹畸变、桨叶弹性变形效应和发动机动态特性。采用基于有限元分析的挥舞-摆振-扭转耦合的弹性桨叶模型,并利用一种新的数值方法将旋翼/机体耦合运动方程表示为显式形式,整个飞行动力学模型表示为状态空间格式。以UH-60A直升机在高速飞行条件下的急拉杆机动飞行为例进行仿真计算,并与飞行试验数据进行对比验证。分析表明,仿真结果与试验结果吻合,高速飞行条件下机体抬头过程中前行桨叶非定常气动载荷的计算误差是引起旋翼和机体运动仿真误差的主要原因。     

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

为准确计算直升机旋翼在复杂的桨-涡干扰(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导致的升力波动幅值和相位的计算结果均与试验值吻合良好.      

直升机旋翼非定常气动力及结构动力响应计算

以二维非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程为模型方程,采用Ｊａｍｅｓｏｎ有限体积并引入变系数隐式残值光顺及固接网格计算了直升机旋翼不同展向位置处翼剖面的非定常气动力,进而采用插值的方法模拟直升机在悬停及前飞状态下,旋翼的桨叶所受的非常气动力载荷,并将之与桨叶的结构振动方程相耦合,求解其时域内的结构响应,非定常下气动力计算结果与实验结果做了对比,两者比较吻合,旋翼结构响应计算结果合理。     

航空发动机风扇叶片冰雹撞击仿真

为预估航空发动机风扇转子叶片受到冰雹撞击后的损伤情况,基于PAM-CRASH软件进行冰雹撞击风扇转子叶片仿真。采用SPH方法和带失效应变的弹塑性材料模型建立冰雹数值模型,模拟冰雹撞击铝合金平板过程,仿真结果与试验数据吻合较好。针对冰雹撞击旋转状态风扇转子叶片试验,建立3维风扇转子有限元模型,使用带失效模型的J-C本构模型定义叶片材料性能,采用该冰雹模型对试验过程进行仿真,获得的冰雹撞击过程和叶片损伤与试验结果基本相同,叶片凹陷深度误差小于10%。仿真与试验结果对比表明:风扇叶片冰雹撞击仿真方法能够预估叶片被冰雹撞击后的损伤情况,可用于叶片抗冰雹撞击设计与评估。     

转静件碰摩状态下的叶片振动载荷和振动特性测试分析

研究转静件碰摩状态下叶片载荷测试方法,测试叶片在碰摩力的冲击与摩擦力作用下的振动特性,给出叶片和机匣受冲击力的时间载荷和频谱,证明叶片在各种碰摩力作用下不发生共振,只受弯矩作用等现象.为旋转机械转静件碰摩作用下叶片的强度分析提供依据.      

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

根据航空发动机支承锥壁结构受力特点,对风扇叶片飞失冲击载荷作用下的锥壁失效破坏机理进行了研究。利用显式动力学有限元仿真方法,对冲击载荷作用下的锥壁结构动态失效过程进行了瞬态分析。开展了对锥壁的落锤冲击试验,并与分析结果进行了对比验证。试验和分析结果表明:冲击载荷作用下锥壁减薄处破坏为剪切失效破坏。采用的显式动力学有限元仿真方法为准确模拟冲击载荷作用下的锥壁结构失效提供了一种可行的仿真手段,分析获得的峰值加速度与试验结果误差小于5%;利用GISSMO(generalized incremental stress state dependent damage model)可以准确预测锥壁减薄处断裂时间、断裂位置。经过试验验证的分析方法及失效模型可运用到风扇叶片飞失冲击载荷下的锥壁失效设计参数的确定,提高锥壁降载结构设计的精度。      

柔性转子动力学特性及支承结构安全性设计

叶片飞失极限状态下,航空发动机转子系统承受冲击载荷和大不平衡载荷,载荷通过转子结构传递到支承结构,对支承结构造成严重损伤。从安全性的角度看,需要对航空发动机轴承-支承结构在极限载荷下的承载能力进行有效的定量评估。提出一种新型缓冲阻尼支承结构,通过对支点刚度突降和高阻尼设计,实现在止推轴承支点处大幅降低横向冲击载荷的影响。建立考虑支承结构刚度突变和阻尼特性的转子系统动力学方程,并计算突加不平衡激励下,转子系统支点动载荷分布的变化规律。结果表明:通过优化设计支点刚度和阻尼参数,缓冲阻尼支承结构,能够有效降低突加不平衡激励下止推滚珠轴承的支点动载荷,提高支承结构的安全性。     

叶片丢失后发动机整机响应模拟试验与仿真

依据大涵道比涡扇发动机结构设计了叶片丢失试验台,开展了一系列模拟叶片丢失试验,并采用显式有限元方法进行了数值仿真,研究了发动机叶片丢失后整机结构响应与载荷传递规律。结果表明:丢失叶片与机匣存在叶尖与叶身两次撞击,对应的加速度曲线存在两个响应峰值;转子转速越高,加速度响应幅值越大。叶片飞断后转子不平衡载荷传递路径为前轴承支承-中轴承支承-中介支板-机匣结构;叶片撞击机匣导致的冲击载荷则由风扇机匣向后传递,最终传给吊装结构;发动机承受的载荷是由不平衡和冲击影响耦合得到,其中冲击载荷为主要部分。该研究为掌握真实发动机叶片丢失下整机响应规律提供了试验模拟方法与数值仿真分析工具。      

直升机主桨叶的鸟撞有限元数值模拟

以某型直升机主桨叶为例,使用瞬态动力学分析软件(MSC.Dytran),采用先进的任意拉格朗日-欧拉(ALE)流固耦合算法,进行鸟撞过程的数值模拟.模拟结果表明:撞击过程中的最大应力达到199MPa,没有超出桨叶材料的屈服极限.研究结果为桨叶的抗鸟撞损伤设计提供了参考.      

DZ125带小孔构件低循环/保载疲劳试验与分析

进行了定向结晶合金DZ125带小孔构件在高温下的低循环/保载疲劳试验.此构件试验用来部分模拟航空发动机气膜冷却的空心涡轮转子叶片, 试验方案根据空心涡轮转子叶片的几何特征和所承受的温度及载荷状况来设计.介绍了试验方案和试验结果分析, 包括寿命分析、试验裂纹扩展规律分析和金相分析.结果表明小孔应力集中和保载是构件破坏的主要原因, 孔边电火花加工特征对裂纹萌生起到决定性影响.      

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

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