前飞状态直升机旋翼／机身非定常气动干扰的分析

一个计算旋翼／机身／尾迹间非定常气功干扰的分析方法是建立在二阶升力线／全展自由涡模型和机射面元模型基础之上的。通过迭代机身在桨盘平面、阡迹定位点的诱导速度和旋翼／尾迹在机身表面的诱导速度,形成一个耦合的分析模型。在分析中计入了非定常项。作为算例,对两种孤立机身表面的平均压强系数分布以及旋翼机身组合体中机身表面的非定常压强系数分布进行了计算,其结果与实验值相吻合。     

直升机旋翼/机身非定常气动干扰数值分析

基于非定常面元/时间步进全展自由尾迹建立了旋翼/机身非定常气动干扰分析方法。方法中耦合了非定常面元法和时间步进自由尾迹,以准确模拟旋翼非定常气动力、旋翼尾迹及桨叶对机身的非定常干扰效应。为模拟前飞状态下具有升力的机身,将机身离散为汇/偶极子面元,并采用涡线镜像法模拟旋翼尾迹靠近机身表面产生的加速效应。通过计算前飞状态的Maryland、ROBIN(Rotor Body INteraction)旋翼/机身干扰下的非定常压力分布,并与可得到的实验值、CFD计算结果对比,验证方法的准确性。随后分析前飞速度、旋翼与机身距离对旋翼/机身非定气动干扰的影响。计算结果表明机身头部和中部非定常压力主要受桨叶的通过性影响,而机身尾梁主要受尾迹/机身干扰影响,机身非定常气动力频率为桨叶片数的倍频。随前飞速度的增加,机身非定常压力幅值增加,尾梁压力幅值先增加后减小;增加旋翼与机身距离将减小机身和尾梁非定常压力幅值。     

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

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

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

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

基于非定常面元/黏性涡粒子混合法的旋翼/平尾非定常气动干扰

旋翼/平尾非定常气动干扰是导致直升机纵向“抬头（Pitch-up）”现象的主要原因。为在直升机设计阶段准确分析旋翼/平尾非定常气动干扰以及由此引起的低速纵向操纵特性变化,通过涡量等效原则和Neumann物面边界条件建立了适用于旋翼/平尾气动干扰分析的非定常面元/黏性涡粒子混合法。该方法耦合了考虑尾迹时变效应的非定常面元法、黏性涡粒子法及涡量镜面法,以准确模拟旋翼和平尾的非定常气动载荷、旋翼尾迹的非定常特性以及旋翼尾迹对平尾的气动干扰效应。首先通过计算NASA ROBIN（Rotor Body Interaction）旋翼尾迹几何和诱导速度分布,并与实验测量值、时间精确自由尾迹及CFD计算结果对比验证方法的准确性。相比于时间精确自由尾迹,本文方法计算精度更高。随后分析了旋翼/平尾非定常气动干扰对平尾向下气动载荷和气动导数的影响,并分析了平尾构型对旋翼/平尾非定常气动干扰的影响规律。分析表明：旋翼尾迹与平尾干扰导致低速状态的平尾载荷突增,气动导数反号;低平尾气动载荷突增较大,高平尾较小,但高速气动导数反号;前置平尾载荷突增量减小,但对应速度范围较宽;右旋直升机右平尾载荷突增量较小,但气动导数特性基本不变。     

旋翼地面效应的气动建模与特性

地面对旋翼气动特性影响明显,且导致旋翼流场更加复杂。为分析地面效应下的旋翼桨尖涡和流场变化特性,基于涡面和无滑移边界条件,求解第2类Fredholm方程获得地面涡面矢量分布,且将涡面矢量按涡扩散方程扩散到流体中,建立考虑黏性效应的地面气动模型,并耦合非定常面元/黏性涡粒子混合法以体现旋翼桨叶气动特性和旋翼尾迹的非定常效应,构建旋翼地面效应气动分析方法。通过计算Lynx尾桨地面效应下的性能和桨尖涡轨迹,并计算Maryland大学模型旋翼和NASA缩比旋翼地面效应下的垂向、径向速度分布,且与试验和CFD计算结果对比,验证了本文方法能较好捕捉地面效应下的旋翼尾迹变化特性和复杂速度场特性,且结果表明本文方法能较好模拟地面效应下旋翼桨尖涡的收缩、扩散、井喷、地面射流等物理现象。     

神经网络在旋翼／机身气动干扰模型中的应用

由于直升机自身的特点，旋翼／机身气动干扰呈现非线性，且受多种因素的影响，用神经网络来解决这一非线性问题是一个很好的办法。将旋翼／机身气动干扰试验数据构造的学习样本，对网络参数进行学习，可以得到旋翼／机身气动干扰神经网络模型，进而可以用该神经网络模型研究直升机机身受旋翼气流干扰时的空气动力特性。参数训练好的旋翼／机身气动干扰神经网络模型可直接用于直升机气动设计和实时仿真。笔者在对旋翼／机身气动干扰神     

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

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

直升机机身干扰对旋翼气动与噪声特性影响研究

本文采用基于CFD/FW-H的旋翼气动/噪声计算模型,针对悬停及中小速度前飞时,考虑机身对旋翼涡流场、非定常气动载荷和噪声特性的干扰影响,开展了直升机机身干扰状态下的旋翼气动和噪声特性数值计算分析研究,获得了不同飞行状态下,机身对旋翼涡流场、气动力和噪声特性的影响规律.计算结果表明,机身对旋翼气动力的影响主要存在于0°...     

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

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

基于Fw-H方程的直升机尾桨与涡线干扰噪声分析

本文通过对直升机尾桨噪声的机理进行分析,建立了尾桨与涡线干扰噪声的计算模型.计算模型中包括气动计算和声学计算,气动计算的结果作为声学计算的输入量.气动计算中分别采用三维非定常面元法计算桨叶表面压力和有扰动薄翼理论计算涡线干扰下桨叶的压力增量;声学计算中采用推导自FW-H方程的Farassat 1a公式,获得声压的时间历程.通过算例分析得到了一些有意义的结论.     

Numerical simulation of aerodynamic interaction for a tilt rotor aircraft in helicopter mode

A rotor CFD solver is developed for simulating the aerodynamic interaction phenomenon among rotor, wing and fuselage of a tilt rotor aircraft in its helicopter mode. The unsteady Navier–Stokes equations are discretized in inertial frame and embedded grid system is adopted for describing the relative motion among blades and nacelle/wing/fuselage. A combination of multi-layer embedded grid and ‘‘extended hole fringe" technique is complemented in original grid system to tackle grid assembly difficulties arising from the narrow space among different aerodynamic components, and to improve the interpolation precision by decreasing the cell volume discrepancy among different grid blocks. An overall donor cell searching and automatic hole cutting technique is used for grid assembly, and the solution processes are speeded up by introduction of Open MP parallel method. Based on this solver, flow fields and aerodynamics of a tilt rotor aircraft in hover are simulated with several rotor collective angles, and the corresponding states of an isolated rotor and rotor/wing/fuselage model are also computed to obtain reference solution.Aerodynamic interference influences among the rotor and wing/fuselage/nacelle are analyzed,and some meaningful conclusions are drawn.     

Non-synchronous blade vibration analysis of a transonic fan

This study numerically investigates the aeromechanic behavior of a transonic fan model with a flat tip-leading-edge on the NASA rotor 67 test case. Single-passage unsteady calculations at a near stall operating point of 82% design speed show that the dominant frequencies of mass flow were not the harmonics of the rotor rotational frequency. A full-annulus fluid–structure interaction analysis was subsequently carried out to examine the unsteady flows and their interactions with blade vibrations. The results show that the modal displacement of the backward traveling seventh nodal diameter of the second torsion mode grew exponentially, which reveals that the blade vibration was non-synchronous. The vibration pattern indicates that the aerodynamic mode was resonant with the structural vibration mode. Around the rotor tip, the circumferential vortical propagation induced by interactions among the main flow, tip leakage flow, and tip clearance vortex was the source of aerodynamic excitation. To clarify the mechanism of the non-synchronous vibration, the coupling between aerodynamic disturbance and structural response, i.e., aliasing, was summarized. The frequency spectra of the fluctuating pressure show that an aerodynamic Backward Traveling Wave (BTW) was co-aliased to a structural BTW due to the propagation of the circumferential vortex. The correlation between the frequency and free convective speed of the aerodynamic disturbance determined the directions of aliasing.     

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。      

倾转旋翼机多部件对机翼气动干扰的分析及优化

为研究倾转旋翼机各部件对机翼的气动干扰效应,建立了一套适用于倾转旋翼机流场CFD求解方法,提出并生成了一套由三棱柱/四面体组成的适用于倾转旋翼机多种飞行状态的非结构混合网格系统.以三维Navier-Stokes方程为主控方程,采用动量源方法进行倾转旋翼的模拟,并引入了高效的OpenMP并行加速技术等,提出了一套新型的动量源项添加及搜索方法.计算分析了倾转旋翼机旋翼/机身/短舱对于机翼气动特性的干扰影响,得出了一些对设计有指导意义的结论.采用基于径向基函数（RBF）的代理模型方法,根据机翼不同展向位置的干扰程度对各段翼型配置进行气动优化.优化结果表明:考虑气动干扰作用,在巡航速度下优化后的全机升阻比增长达到了36.78%.