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

 针对直升机大机动飞行仿真,建立了一个非线性的飞行动力学模型,考虑了翼型非定常/动态失速、机动飞行引起的动态尾迹畸变、桨叶弹性变形效应和发动机动态特性。采用基于有限元分析的挥舞-摆振-扭转耦合的弹性桨叶模型,并利用一种新的数值方法将旋翼/机体耦合运动方程表示为显式形式,整个飞行动力学模型表示为状态空间格式。以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导致的升力波动幅值和相位的计算结果均与试验值吻合良好.      

直升机动态失速研究

本文简要介绍了旋翼和翼型的动态失速特性、动态失速主动和被动控制、旋翼综合法预测使用的一些动态失速模拟方法.UH-60A直升机飞行试验11029飞行是一个典型的旋翼动态失速实例,表演出旋翼后行桨叶上出现高攻角引起的动态失速,而前行桨叶上则出现因为激波诱导前缘分离引起的动态失速.翼型的动态失速特性方面介绍了典型的动态失速迟滞回线(攻角正弦变化和锯齿变化)、缩减频率和雷诺数对翼型动态失速特性的影响等.动态失速主动控制介绍可动前缘翼型的情况,被动控制则是使用前缘涡发生器的情况.     

加装格尼襟翼旋翼的直升机飞行性能

为研究加装格尼襟翼旋翼的直升机飞行性能,建立了加装格尼襟翼旋翼的直升机飞行动力学模型。采用UH-60A直升机试飞数据验证了计算模型的正确性。在此基础上,分析了样例直升机加装格尼襟翼后重量系数、格尼襟翼高度、沿径向位置和加装方式对旋翼需用功率的影响,以及加装格尼襟翼后旋翼桨叶剖面迎角分布、旋翼操纵量和机身姿态角的变化等。研究表明,直升机在重量系数较大的状态下高速前飞时,旋翼加装格尼襟翼能够明显降低直升机的需用功率,且加装转动格尼襟翼的效果优于加装固定格尼襟翼。功率降低幅值随格尼襟翼高度的增加先增加后减小。格尼襟翼在桨叶上布置的位置越靠近桨尖,其对需用功率的影响越大。直升机在重量系数较大的状态下高速前飞时,加装格尼襟翼能够使旋翼后行侧最大迎角显著减小。加装格尼襟翼后旋翼总距和纵横向周期变距减小。     

直升机飞行动力学高阶线性系统建模

针对常规单旋翼带尾浆直升机,从复杂的飞行动力学非线性数学模型入手,采用数值方法在平衡点处求出线性模型;该线性模型包括传统的６自由度刚性机体模型,还包括主桨和尾浆动力入流、主桨和尾桨挥舞自由度,其状态空间形式具有２５个状态变量（５简称２５状态模型）。以某型直升机为例,对线性模型与非线性模型的时域动态响应,及２５状态高阶线性模型和传统６自由度９状态模型的特征值进行了比较,验证了建模方法的可靠性和精度。     

直升机飞行动力学数学建模问题

直升机飞行动力学数学模型是飞行控制系统设计的基础,也是直升机飞行品质设计和评估的主要手段。直升机是一个多体系统,在直升机飞行动力学建模过程中,必须考虑旋翼、机体与升力面等的运动耦合、惯性耦合、结构耦合和气动耦合以及非定常、非线性特性,给出各个运动部件的物理模型及其数学表达形式,是对不同假设、子模型进行分析和综合的一个复杂的过程。鉴于此,简要回顾了单旋翼带尾桨直升机飞行动力学数学模型的研究现状,着重描述了直升机飞行动力学数学建模中的旋翼气动力建模、直升机气动干扰建模、旋翼/发动机建模以及直升机飞行动力学模型的集成与综合的研究现状与研究进展。最后,针对直升机飞行动力学的数学建模提出了今后的研究重点。     

独立桨距控制对直升机飞行性能的影响

为研究旋翼独立桨距控制对直升机飞行性能的影响,在已验证的直升机性能分析模型基础上,以UH-60A直升机为样例,通过输入不同阶次、幅值和相位角的独立桨距,分析直升机的需用功率和升阻比变化,在此基础上给出了旋翼桨盘内迎角和阻力系数的分布,探讨了独立桨距控制提升直升机性能的机理。分析结果表明：高速时,2阶和3阶的独立桨距控制可降低旋翼需用功率、提高直升机的升阻比,但提升效果有限,直升机起飞重量较大时,效果更明显;2阶耦合3阶的独立桨距控制对旋翼性能的提升效果比单独的2阶或3阶桨距输入更好,样例直升机的需用功率最多降低了4.5%,4阶的独立桨距输入则不利于提升旋翼性能;输入独立桨距后,旋翼桨盘迎角分布改善,后行侧迎角减小,有利于推迟失速,桨盘的后行侧的阻力系数减小,可有效降低旋翼的需用功率,提升直升机飞行性能。     

直升机桨叶气动外形多目标优化设计

基于自由尾迹方法建立了直升机桨叶空气动力学分析模型,应用人工神经网络方法建立代理模型,采用改进的多目标遗传算法构建了优化框架,对直升机的悬停和大速度前飞状态进行优化.以悬停效率、旋翼等效升阻比及桨叶叶素的最大阻力系数为约束,对两个飞行状态的需用功率进行优化,得到了Pareto最优解集.并以UH-60A直升机的桨叶为算例,对其外形进行优化设计,优化结果表明,提出的桨叶气动外形多目标优化框架是有效可行的.      

旋翼结冰对直升机飞行动力学特性的影响

基于旋翼冰风洞试验数据,考虑桨叶表面附着冰脱落及桨叶表面局部温度对结冰的影响,提出了旋翼结冰的工程模型,建立了直升机旋翼结冰后的飞行动力学模型,以UH-60A为样机,研究了直升机旋翼结冰后的平衡特性,分析了结冰对直升机稳定性的影响.根据有人驾驶垂直/短距起落飞机军用品质规范(MIL-F-83300)与军用旋翼飞行器驾驶...     

失衡旋翼的直升机自激振动分析模型

直升机使用中很可能出现旋翼各片桨叶特性不一致的情况,为了研究失衡旋翼对直升机自激振动的影响,建立了适用于地面、悬停及前飞状态的旋翼/机体耦合动稳定性分析模型。采用当量铰旋翼模型,计入动力入流的影响,分别在旋转坐标系和固定坐标系中建立了桨叶及机体的动力学方程。以减摆器失效对直升机地面共振的影响为例,对桨叶及机体的时域响应进行了非线性数值仿真,用Floquet传递矩阵法计算了摆振后退型模态频率及阻尼,并用时域分析进行了检验。结果表明,其中一个减摆器失效后,各片桨叶摆振运动特性相差很大,系统的摆振后退型模态阻尼下降幅度高达60%以上。     

直升机机动飞行的逆模拟

本文给出了一种直升机机动飞行的逆模拟方法以计算跟随预定飞行轨迹的驾驶员操纵,根据这一方法可以确定为完成直升机机动飞行所需的驾驶员操纵输入及直长机的飞行速度、角速度和的变化历程。直同飞行动力学模型没有作任何线化假设,其中考虑了旋翼入流的时滞效应、前行桨叶的压缩性物后行桨叶的失速特性及旋翼桨叶的非定常挥舞运动,引入了旋翼尾迹对直升机机身、尾翼和尾桨的气动干扰。最后以黑鹰直升机为例计算了鱼跃越障机动飞行     

基于矩阵运算的直升机旋翼/机体耦合动力学模型

针对直升机旋翼/机体耦合动力学建模中,小角度线性简化带来的计算精度下降以及适用范围缩小等问题,在计算桨叶速度及加速度时,基于矩阵运算,保留了运动的非线性。对不同机体运动幅值和频率情况下的小角度线性简化对桨叶速度的影响进行了仿真分析,结果发现,随着机体运动幅值的增加,线性简化影响越来越大,滚转幅值 15°时,简化带来的误差超过 5%;俯仰幅值 15°时,简化带来的误差超过 6%,且随着机体运动频率的减小,误差越来越大。对不同桨距的情况进行仿真,结果表明,所建模型可以很好地模拟直升机的动态响应,采用该模型可进行直升机旋翼/机体耦合动态响应仿真分析。     

非线性叶间黏弹减摆器对直升机空中共振的影响分析

 建立带非线性叶间黏弹减摆器的直升机旋翼/机体耦合动稳定性分析模型。与全机飞行力学平衡计算相结合，旋翼/机体耦合动稳定性分析模型考虑前飞状态桨叶变距操纵、机体姿态角和桨毂纵向安装角。针对具有非线性特性的叶间黏弹减摆器，采用基于复模量的非线性VKS改进模型、Simulink时域仿真和多桨叶坐标变换等效阻尼识别法分析直升机悬停、前飞状态下旋翼/机体耦合动稳定性及减摆器双频动幅值，并就减摆器布局、全机总重以及前飞速度对桨叶摆振后退型模态阻尼的影响进行分析。结果表明：由悬停到前飞直升机动稳定性一般均下降，一定速度后又上升；加上减摆器能消除前飞不稳定区；叶间黏弹减摆器抬头连接能提高模态阻尼。     

Flight dynamics characteristics of canard rotor/wing aircraft in helicopter flight mode

The aerodynamic layout of the Canard Rotor/Wing(CRW) aircraft in helicopter flight mode differs significantly from that of conventional helicopters. In order to study the flight dynamics characteristics of CRW aircraft in helicopter mode, first, the aerodynamic model of the main rotor system is established based on the blade element theory and wind tunnel test results. The aerodynamic forces and moments of the canard wing, horizontal tail, vertical tail and fuselage are obtained via theoretical analysis and empirical formula. The flight dynamics model of the CRW aircraft in helicopter mode is developed and validated by flight test data. Next, a method of model trimming using an optimization algorithm is proposed. The flight dynamics characteristics of the CRW are investigated by the method of linearized small perturbations via Simulink. The trim results are consistent with the conventional helicopter characteristics, and the results show that with increasing forward flight speed, the canard wing and horizontal tail can provide considerable lift,which reflects the unique characteristics of the CRW aircraft. Finally, mode analysis is implemented for the linearized CRW in helicopter mode. The results demonstrate that the stability of majority modes increases with increasing flight speed. However, one mode that diverges monotonously,and the reason is that the CRW helicopter mode has a large vertical tail compared to the conventional helicopter. The results of the dynamic analysis provide optimization guidance and reference for the overall design of the CRW aircraft in helicopter mode, and the model developed can be used for control system design.     

Hover performance of helicopter main and tail rotors with swirl velocities

It is important to quickly predict the hover performance of main and tail rotors with sufficient precision for helicopter design. To investigate the effects of swirl velocities on the hover performance of main and tail rotors, and give a better prediction for the hover performance, a flight performance model was derived and a swirl velocity model was coupled into it. The test data of the UH-60A helicopter were used for validation. When the blade loading coefficient of the main rotor was higher than 005, the effects of the swirl velocities on the main rotor power became significant. The swirl velocities increased the profile torque of the main rotor. The increased torque required the tail rotor to produce more thrust with more power consumption. At a higher blade loading coefficient of the main rotor of 012, the swirl velocities increased the main rotor power, tail rotor power and total power by 380%, 524% and 508%, respectively. The profile power increase of the main rotor caused by the profile swirl velocity was less than that of the induced swirl velocity, but the power increase was higher at high rotor blade loadings. Considering the swirl velocities in the main rotor can improve the prediction precision of the hover performance, especially at high blade loadings      

Dynamic response analysis under atmospheric disturbances for helicopters based on elastic blades

A flight dynamics model based on elastic blades for helicopters is developed. Modal shape analysis is used to describe the rotating elastic blades for the purpose of reducing the elastic degrees of freedom for blades. The analytical result is employed to predict the rotor forces and moments. The equilibrium equation of the flight dynamics model is then constructed for the elastic motion for blades and the rigid motion for other parts. The nonlinear equation is further simplified, and the gradient descent algorithm is adopted to implement the trim simulation. The trim analysis shows that the effect of blade elasticity on the accuracy of rotor forces and moments is apparent at high speed, and the proposed method presents good accuracy for trim performance. The time-domain response is realized by a combination of the Newmark method and the adaptive Runge-Kutta method. The helicopter control responses of collective pitch show that the response accuracy of the model at a yaw-and-pitch attitude is improved. Finally, the influence of blade elasticity on the helicopter dynamic response in low-altitude wind shear is investigated. An increase in blade elastic-ity reduces the oscillation amplitude of the yaw angle and the vertical speed by more than 70%. Compared with a rigid blade, an elastic blade reduces the vibration frequency of the angular veloc-ity and results in a fast return of the helicopter to its stable flight.     

Investigation of rotor control system loads

This paper concentrates on the aeroelasticity analysis of rotor blade and rotor control systems. A new multi-body dynamics model is established to predict both rotor pitch link loads and swashplate servo loads. Two helicopter rotors of UH-60A and SA349/2, both operating in two critical flight conditions, high-speed flight and high-thrust flight, are studied. The analysis shows good agreements with the flight test data and the calculation results using CAMRAD II. The mechanisms of rotor control loads are then analyzed in details based on the present predictions and the flight test data. In high-speed conditions, the pitch link loads are dominated by the integral of blade pitching moments, which are generated by cyclic pitch control. In high-thrust conditions, the positive pitching loads in the advancing side are caused by high collective pitch angle, and dynamic stall in the retreating side excites high-frequency responses. The swashplate servo loads are predominated by the rotor pitch link loads, and the inertia of the swashplate has significant effects on high-frequency harmonics of the servo loads.     

旋翼变体技术对直升机性能的提升

为研究旋翼变体技术对直升机性能的提升作用,先建立旋翼模型,然后耦合机体模型,从而建立直升机需用功率计算模型.主要对比了几种不同旋翼变体技术,包括旋翼变直径、旋翼变转速、桨叶变弦长和桨叶变负扭转角在不同飞行状态时对直升机需用功率的影响.前飞速度为130km/h时,10%旋翼转速减小、10%旋翼直径减小、10%桨叶弦长减小和桨叶负扭转角由-12°变为-6°时,需用功率分别降低了15.7%,14.6%,5.8%和3.1%;前飞速度为250km/h时,10%旋翼转速减小和10%旋翼直径减小可分别降低14.5%和23.9%需用功率.结果表明,旋翼变转速明显优于桨叶变弦长和桨叶变负扭转角所取得的性能提升,高速前飞时旋翼变直径降低的需用功率大于旋翼变转速技术.      

重型直升机飞行动力学刚弹耦合建模及空中共振稳定性分析

针对重型直升机（HLH）大重量、低转速的固有特性,提出了一种适用于重型直升机的飞行动力学刚弹耦合建模方法。该方法结合传统直升机飞行动力学与旋翼机体耦合动力学,将传统飞行力学的分析频段拓展到了5 Hz,额外考虑了桨叶和机体的弹性变形,基于阻抗匹配法推导出了显式的旋翼/机体耦合动力学方程,模拟了真实飞行状态下的直升机气弹耦合特性,利用该模型计算并分析了算例重型直升机的悬停飞行特性和空中共振稳定性。结果表明：旋翼机体耦合导致摆振前进型和机体弹性模态的阻尼-转速曲线先相互靠近至同一点再分离,可能引起直升机的高频瞬态振动;在摆振等效阻尼不足时,旋翼摆振后退型是不稳定的,但随着等效阻尼增加,摆振二阶周期型模态和机体弹性模态会出现耦合;桨叶弹性变形与机体弯曲模态及挥舞集合型耦合,但不会引起明显的不稳定现象。