-
引言
-
一直以来,常规构型直升机都存在着航程短、巡航速度低的缺点[1],世界各大直升机研发机构从未停止过探索高速直升机,以突破直升机速度限制,其中采用前行桨叶概念的[2]共轴高速直升机是未来直升机发展的重要方向之一,该构型高速直升机历经几十年发展,目前已突破370 km/h巡航速度,共轴高速直升机也历经XH-59A,X-2,S-97和SB>1等多型迭代(SB>1于2019年初实现首飞[3]),正朝着400 km/h速度发展.共轴高速直升机是未来直升机发展的重要方向之一.
-
但恶劣的旋翼振动载荷问题也始终给共轴高速直升机研发带来挑战.因其非对称性的流场环境,使得旋翼一直是直升机主要振源,而其中气流分离、动态失速以及桨涡干扰等非定常气动问题又引起更为复杂的高阶结构载荷问题[4].降低旋翼转速,推迟前行侧激波阻力发生是高速直升机实现高速前飞的策略之一,但低转速下的旋翼整体处在大前进比状态,且桨盘后行侧存在大面积反流区现象,桨叶内段经过该区域时会产生高负升力和高阻力,引起周期性脉冲式交变气动载荷,进而引起以2/rev(基频为旋翼转速频率)为主要频率的挥舞、摆振结构载荷[5,6],也使得桨尖挥舞形变也具有2/rev主要变化特征[7].升力偏置是高速直升机实现高速前飞的另一大策略,但在高升力偏置下会造成显著的以1/rev为主要特征的挥舞结构载荷以及桨尖挥舞形变问题[8,9],同时减小双旋翼桨尖间距,存在潜在的安全问题[10].
-
XH-59A飞行试验[11]显示驾驶舱3/rev(3片桨叶)垂向振动在461.7 km/h状态下达到0.5g而横向振动高达0.75g.X-2载荷实测结果也表明桨毂4/rev(4片桨叶)振动载荷水平在444 km/h飞行状态下急剧增加[12].
-
在美国德克萨斯州立大学的共轴双旋翼振动载荷试验[13]中发现上下旋翼桨叶交叉作用最大能引起相当于10%旋翼拉力静态值和30%旋翼拉杆载荷静态值的冲击效应.当前共轴刚性旋翼振动载荷分析仍然是动力学研究中的难点之一.
-
针对共轴刚性旋翼结构载荷问题,本文建立共轴刚性旋翼气弹综合分析模型,针对共轴刚性旋翼挥舞、摆振的低阶动载荷展开详细的参数影响分析,包括旋翼转速、前进比以及升力偏置重要参数,通过参数影响分析总结共轴刚性旋翼挥舞、摆振低阶动载荷特性及变化规律,并得出相应结论.
-
1 共轴刚性旋翼气弹模型
-
1.1 共轴双旋翼干扰模型
-
基于Pitt-Peters线性入流模型[14]并引入双旋翼气动干扰因子建立共轴双旋翼入流模型,其干扰形式下的旋翼桨盘诱导速度表达式如下:
-
式中诱导速度v下标1和2代表上、下旋翼(1为下),vu,t为上、下旋翼最终诱导速度,依据动量守恒定律引入气动干扰因子[7]12(上旋翼对下旋翼干扰),共轴双旋翼几何结构和气动干扰示意分别如图1和图2所示.
-
图1 双旋翼干扰
-
Fig.1 Coaxial rotor interaction
-
图2 共轴双旋翼直升机
-
Fig.2 Coaxial rotor helicopter
-
1.2 共轴刚性旋翼气弹模型
-
本文基于中等变形梁理论[15]建立桨叶结构动力学模型,通过哈密尔顿原理建立旋翼系统动力学方程:
-
其中δU为系统应变能变分项,其变分形式如下,其中ε为应变矩阵,σ为材料应力矩阵.
-
δT为动能变分项,包括桨叶旋转和桨叶变形运动产生的动能,其变分形式如下∶qb为广义自由度矩阵,Mb为桨叶惯性运动引起的质量矩阵,FM动能的非线性广义力.
-
δW为外载荷做功项,包括气动力和气动力矩两部分.其中Fα为气动力,为Mα气动力矩,为Gb广义力系数矩阵.
-
采用5节点15自由度梁单元对桨叶进行离散(图3).其中单元自由度中u,v,w,φ表示拉伸、摆振、挥舞和扭转位移,方向定义为桨叶向上挥舞运动w、逆旋转方向摆振运动v以及桨叶抬头扭转φ为正.根据哈密顿原理推导并得到每个单元的质量阵M、阻尼阵C、刚度阵K以及外载荷项F,最后通过有限元总体阵组集得到最终旋翼动力学方程,其中q为系统所有自由度的列向量.
-
图3 有限元单元自由度
-
Fig.3 Degree of Finite element
-
气动力模型采用非线性准定常模型,并针对共轴高速直升机的特殊工作环境,进行大前进比反流区气动力计算修正以及建立适用于大前进比状态的提升收敛性的配平迭代方法[5].
-
与常规单旋翼直升机操纵方式不同,共轴高速直升机操纵策略中通过设置差动变距来调节上下旋翼操纵量的差异,也通过差动变距的方式设置升力偏置,进而实现高速前飞,其上下旋翼操纵表达式如下:
-
上式中,下标U,L分别表示上、下旋翼,其中θ0为总距,用于控制旋翼总拉力; Δθ0为差动总距,控制上下旋翼拉力差,调节双旋翼扭矩,产生偏航力矩; A1为纵向周期变距,A′1为差动纵向周期变距,二者用来调节旋翼总俯仰力矩与上下旋翼俯仰力矩差; B1为横向周期变距,B′1为差动纵向周期变距,用来调节旋翼总滚转力矩与上下旋翼滚转力矩差; Γ为提前操纵角,反映驾驶员操纵与旋翼响应的滞后关系.
-
在本文中设置4个配平目标,分别是双旋翼总拉力Fz,总俯仰力矩My和总滚转力矩Mx以及升力偏置(Lift Offset)[表达式如式(4)],对应的配平操纵量依次为θ0,A1,B1,B′1.
-
上式中,T表示旋翼拉力,上标U,L分别代表上、下旋翼,R表示旋翼半径.
-
2 共轴刚性旋翼低阶动载荷计算分析
-
本文以XH-59A共轴刚性旋翼为分析对象[16],旋翼总体参数见表1.着重分析旋翼转速、前飞速度以及升力偏置对共轴刚性旋翼挥舞、摆振弯矩的低阶动载荷影响.以下算例分析均按照旋翼拉力T=3500kg进行配平计算.
-
首先针对本文建立的共轴刚性旋翼气弹分析方法进行验证.计算XH-59A的旋翼前5阶模态频率并绘制共振图(图4),图中黑点为试验值[16],黑线为计算值.为旋翼转速,额定转速为345rpm,E表示摆振频率,F表示挥舞频率,T表示扭转频率.前5阶模态频率包括挥舞前2阶、摆振前2阶、扭转1阶频率,图中计算值基本贴合试验值,具有较高的准确性,尤其是挥舞1阶、摆振1阶频率结果对比,表明本文建立的旋翼气弹分析方法准确可靠.
-
图4 旋翼共振图
-
Fig.4 Resonance diagram of rotor
-
2.1 旋翼转速参数影响分析
-
图5 上旋翼桨根挥舞弯矩低阶谐波幅值
-
Fig.5 The lower harmonic amplitude of root flap bending moment in upper rotor
-
在前飞速度v=60m/s(前进比μ=0.302)状态下,由于旋翼尾迹被快速吹离桨盘,此时上下旋翼可忽略气动干扰[7],此时分析转速对旋翼动载荷影响.图5和图6是不同转速下的上旋翼挥舞、摆振弯矩的低阶谐波幅值计算结果.结果表明,在低转速下挥舞、摆振弯矩的2/rev谐波幅值变化显著,尤其在挥舞方向.在低转速大前进比状态下,桨叶在前行侧桨尖处以及后行侧桨根处存在较大的负升力区域,且桨叶根部经历周期性交变脉冲气动载荷[5],为满足桨盘整体升力的横向平衡,后行侧桨叶升力向桨尖集中,而前行侧整理升力降低以平衡后行侧桨叶升力,这一现象造成周期内的挥舞弯矩以2/rev谐波为主.图5中相比于100%转速状态,60%转速状态下的挥舞2Ω谐波幅值增加约82.2%,呈指数上升特点.另外,挥舞、摆振弯矩的1/rev、3/rev谐波幅值受转速影响不大.
-
图6结果也有2/rev谐波幅值增加的特点.相比于100%转速状态,60%转速状态下的摆振2/rev谐波幅值增加高达824.4%,尤其是低转速下幅值急剧增加.一方面是因为摆振运动受到挥舞运动影响,2/rev的挥舞运动引起2/rev哥氏力,进而引起2/rev摆振载荷; 另一方面因为降低转速后,摆振1阶频率接近2/rev(见图4),存在共振现象,进而摆振动载荷急剧增加.
-
图6 上旋翼桨根摆振弯矩低阶谐波幅值
-
Fig.6 The lower harmonic amplitude of root lag bending moment in upper rotor
-
2.2 前飞速度参数影响分析
-
图7和图8是前飞速度影响下的上旋翼桨根挥舞、摆振弯矩计算结果.在低速飞行时由于上下旋翼气动干扰效应导致上旋翼拉力大于下旋翼拉力,即上旋翼拉力Fz>T/2,前飞速度增加,后旋翼尾迹加快远离桨盘区域,气动干扰效应减弱,上下旋翼拉力差(ΔFz)逐渐降低,即每幅旋翼拉力都趋于T/2.这就使得上旋翼桨根挥舞、摆振弯矩1/rev谐波幅值先增大后降低.而前飞速度进一步增加,使得旋翼处在大前进比状态,与降低转速原因类似,导致挥舞弯矩2/rev谐波幅值急剧增加.但不同的是,前飞速度增加后挥舞弯矩各低阶谐波幅值都增加,尤其是2、3阶谐波幅值.
-
图7中相比v=15m/s(前进比μ=0.075)时刻,v=120m/s(前进比μ=0.604)时刻的挥舞弯矩2/rev谐波幅值增长约13.7倍,也暴露出大前进比下共轴刚性旋翼振动载荷恶劣的严峻问题.摆振弯矩仍然以1/rev谐波为主,且前飞速度对摆振1/rev谐波载荷影响较大,都存在大前进比状态下各低阶谐波幅值增加的现象.
-
图7 上旋翼桨根挥舞弯矩低阶谐波幅值
-
Fig.7 The lower harmonic amplitude of root flap bending moment in upper rotor
-
图8 上旋翼桨根摆振弯矩低阶谐波幅值
-
Fig.8 The lower harmonic amplitude of root lag bending moment in upper rotor
-
2.3 升力偏置下的旋翼低阶动载荷分析
-
共轴刚性旋翼通过升力偏置策略实现提升前行侧桨叶升力,同时卸载后行侧升力,以提高旋翼气动效率,同时实现高速前飞.但与此同时造成了更大的旋翼滚转力矩以及桨盘横向不平衡性,带来新的动力学问题.
-
图9和图10是两种前飞速度(30m/s和90m/s)下的上旋翼桨根挥舞、摆振弯矩结果图.计算结果表明,在升力偏置影响下,桨根挥舞弯矩峰峰值增加,其弯矩最大值在90°和270°两个方位角处,挥舞弯矩随升力偏置增加而逐渐呈现1/rev谐波为主的变化趋势,且该现象愈发明显.这主要因为升力偏置是通过调节差动横向周期变距来施加,通过1/rev的周期变距变化直接影响旋翼滚转力矩、挥舞响应以及挥舞弯矩,使得挥舞弯矩整体变为1/rev为主的谐波特性.
-
图8 上旋翼桨根摆振弯矩低阶谐波幅值
-
Fig.8 The lower harmonic amplitude of root lag bending moment in upper rotor
-
图9 上旋翼桨根弯矩
-
Fig.9 The root bending moment in upper rotor
-
摆振方面载荷主要受挥舞运动以及挥舞载荷影响,图9(b)和图10(b)显示,在升力偏置影响下,摆振弯矩逐渐降低,尤其是在大速度状态.值得注意的是:大速度状态下由于桨盘后行侧区域的负升力及脉冲式气动效应,直接造成摆振弯矩在270°陡增,而通过施加升力偏置,卸载后行侧升力同时降低了后行侧的阻力,进而摆振弯矩降低,这种陡增现象得到有效缓解,侧面表明升力偏置具有降低摆振载荷的潜在能力.
-
图10 上旋翼桨根弯矩
-
Fig.10 The root bending moment in upper rotor
-
图11 低阶谐波幅值
-
Fig.11 The lower harmonic amplitude
-
图12 上旋翼桨根弯矩
-
Fig.12 The root bending moment in upper rotor
-
图13 上旋翼桨根弯矩
-
Fig.13 The root bending moment in upper rotor
-
为更深入研究升力偏置的影响,选取v=90m/s(前进比μ=0.45)下的挥舞、摆振弯矩进行谐波分解并绘制图11.图中显示,在大速度状态下,升力偏置几乎不改变挥舞弯矩静态值,但随着升力偏置增加,挥舞2阶载荷逐渐降低,而1/rev载荷急剧增加.而摆振弯矩中所有成分均降低,再次说明了升力偏置具有降低摆振载荷的潜在能力.
-
针对不同转速下的升力偏置影响也展开相应的计算分析,图12和图13为前飞速度60m/s下80%和60%两档转速下的上旋翼桨根挥舞、摆振弯矩计算结果.降低旋翼转速和增加前飞速度都会使得旋翼处于大前进比状态下,在80%转速结果中,升力偏置的增加导致挥舞弯矩的1/rev谐波载荷明显增加,2/rev谐波载荷降低,这一现象在60%的挥舞弯矩结果[图13(a)]更为突出.摆振弯矩随升力偏置的增加而整体降低,尤其是2/rev谐波载荷[图13(b)],表明在大前进比状态下,采用升力偏置策略能有效降低摆振弯矩.
-
3 结论
-
本文建立了共轴双旋翼气动干扰模型,集成共轴刚性旋翼综合气弹分析模型,针对共轴刚性旋翼低阶动载荷特性,开展了包括前飞速度、转速以及升力偏置在内的详细的参数影响分析,并得出以下结论:
-
(1)降低旋翼转速和提高飞行速度使得旋翼处在大前进比状态,此时的桨根挥舞、摆振弯矩大大增加,且均以2/rev谐波载荷为主要载荷.
-
(2)在大前进比状态下,采用升力偏置会使得挥舞弯矩1/rev谐波载荷增加,而2/rev谐波载荷降低; 并且升力偏置能整理抑制摆振弯矩.
-
参考文献
-
[1] 吴希明.高速直升机发展现状、趋势与对策 [J].南京航空航天大学学报,2015,47(2):173-179.WU X M.Current Status,development trend and countermeasure for high-speed rotorcraft [J].Journal of Nanjing University of Aeronautics & Astronatics,2015,47(2):173-179.
-
[2] DENG J K,FENG F,LIU P A,et al.Aerodynamic characteristics of rigid coaxial rotor by wind tunnel test and numerical calculation [J].Chinese Journal of Aeronautics,2019,32(2):568-576.
-
[3] SCOTT B,ERIC K.Design and development of the main rotor gearbox for the sikorsky boeing SB>1 DEFIANT JMR technology demonstrator aircraft [C]//The Proceedings of 75th Vertical Flight Society Annual Forum.Philadelphia:Vertical Flight Society,2019:2949-2957.
-
[4] 余智豪,周云,宋彬,等.基于流固耦合的旋翼结构振动载荷计算分析 [J].振动工程学报,2020,32(02):285-294.YU Z H,ZHOU Y,SONG B,et al.Computational analysis on structural vibration loads of rotor based on liquid-solid coupling method [J].Journal of Vibration Engineering,2020,32(02):285-294.
-
[5] 余智豪,周云,宋彬.大前进比变转速旋翼气弹动力学建模与载荷特性分析 [J].振动与冲击,2021,40(04):17-22.YU Z H,ZHOU Y,SONG B.Aeroelastic modeling and load analysis of a variable speed rotor in high advance ratio [J].Journal of Vibration and Shock,2021,40(04):17-22.
-
[6] XING W,LAUREN T,INDERJIT C.Refined performance results on a slowed mach-scaled rotor at high advance ratios [J].Journal of the American Helicopter Society,2020,65(1):1203-12016.
-
[7] 余智豪,周云.含升力偏置的共轴刚性旋翼配平与响应分析 [J].科学技术与工程,2020,20(32):13466-13471.YU Z H,ZHOU Y.Trim and response analysis of coaxial rigid rotor with lift offset [J].Science Technology and Engineering,2020,20(32):13466-13471.
-
[8] JOSEPH H.Aeromechanics of a high speed coaxial helicopter rotor [D].Maryland:University of Maryland,2017.
-
[9] SCHMAUS J,CHOPRA I.Aeromechanics for a high advance ratio coaxial helicopter [C]//The Proceedings of 71st American Helicopter Society International Annual Forum.Fairfax:American Helicopter Society,2015:1139-1153.
-
[10] 周云,胡和平,孟微.共轴刚性旋翼桨尖间距建模与参数影响研究 [J].动力学与控制学报,2020,18(4):52-60.ZHOU Y,HU H P,MENG W.Modeling and parametric study of tip clearance of coaxial rigid rotor [J].Journal of dynamics and control,2020,18(4):52-60.
-
[11] ROBERT K,BURGESS.The ABCTM rotor a historical perspective [C]//The Proceedings of 60th American Helicopter Society International Annual Forum Vertical Flight Transformation,Baltimore:American Helicopter Society,2004:2306-2352.
-
[12] BLACKWELL R,MILLOTT T.Dynamics design characteristics of the sikorsky X2 technologyTM demonstrator aircraft [C]//The Proceeding of 64th American Helicopter Society International Annual Forum.Montreal:American Helicopter Society,2008:886-898.
-
[13] DAIJU U,JAYANT S,ROLAND F,et al.Blade passage loads and deformation of a coaxial rotor system in hover [J].Journal of Aircraft,2019,56(6):2144-2157.
-
[14] PITT D M,PETERS D A.Theoretical prediction of dynamic inflow derivatives [J].Vertica,1981,5(1):21-34.
-
[15] 翁晨.直升机稳态平飞状态旋翼桨叶振动载荷计算的气弹模态法 [D].南京:南京航空航天大学,2014.WENG C.The prediction and analysis of rotor vibration loads based on aeroelastic characteristics in the helicopter steady level flight condition [D].Nanjing:Nanjing University of Aeronautics and Astronautics,2014.
-
[16] FELKER F F.Performance and loads data from a wind tunnel test of a full-scale,coaxial,hingeless rotor helicopter:NASA TM8129 [R].Washington:NASA,1981.
-
摘要
针对共轴刚性旋翼低阶动载荷问题开展分析研究.建立考虑升力偏置的共轴刚性旋翼气弹配平计算方法,建立含干扰因子的共轴双旋翼气动干扰模型与考虑双旋翼构型的旋翼动力学模型,最终集成共轴双旋翼考虑升力偏置的共轴刚性旋翼气弹配平模型.在分析前进比、转速以及升力偏置参数影响下的旋翼低阶挥舞、摆振弯矩变化特性中发现:在大前进比状态下,桨根挥舞、摆振弯矩大大增加,且均以2/rev谐波载荷为主;而施加升力偏置使得挥舞1/rev谐波载荷增加,而2/rev谐波载荷降低,能整体抑制摆振弯矩.
Abstract
An aeroelastic trim of a coaxial rigid rotor with lift-offset was developed and applied to analyzing the lower harmonic load of coaxial rigid rotor. The structural dynamics model and the double-rotor aerodynamics model were jointly used to develop the coaxial rigid rotor’s vibration load computing method. The lower harmonic load of blade root flap bending moment and blade root lag bending moment were analyzed in detail under the influence of advance ratio, rotor speed and lift offset. The results show that the predominant 2/rev trend in blade root flap bending moment and edge bending moment at high advance ratio, and the predominant 1/rev trend in blade root flap bending moment at high lift offset, the edge bending moment would be suppressed at high lift offset.