粘弹减摆器连接形式对旋翼机体耦合动稳定性的影响.docxVIP

粘弹减摆器连接形式对旋翼机体耦合动稳定性的影响 直升机耦合动稳定性分析模型 飞机的“地面共振”是飞机在地面上行驶时，自身的倾斜力和身体之间的耦合运动的现象。这是飞机动力学的基本问题之一。目前采取的主要措施是在旋翼桨毂中增加粘弹减摆器。 文献采用二维当量模型进行了直升机“地面共振”分析,其机体当量模态参数由试验得到。胡国才建立了计入桨叶挥舞/摆振结构耦合、非定常空气动力影响以及机体瞬时转动模态的旋翼/机体耦合动稳定性分析模型,采用基于Coleman变换的特征分析法得到模态阻尼。考虑桨毂结构尺寸及直升机飞行性能的要求,粘弹减摆器在安装时存在多种连接形式,如普通型、几何耦合型和叶间型。粘弹减摆器几何耦合可能会减小减摆器的静态位移,提高直升机悬停时摆振后退型模态阻尼。N.M.Sela和A.Rosen通过等效刚度系数的对比,发现当桨叶片数较多时改进的叶间减摆器连接形式更有效,减摆器的桨叶-桨叶连接比桨叶-桨毂连接有优势。韩景龙等在桨叶的运动方程中未考虑机身运动,旋翼系统中包含有粘弹减摆器(每桨叶一个)。 考虑到特征分析法概念清晰的特点,本文将在文献的基础上,建立带叶间粘弹减摆器和普通连接粘弹减摆器的直升机旋翼/机体耦合动稳定性分析模型,考虑了机体绕瞬时转动中心的滚转和俯仰运动、桨叶的挥舞自由度和非定常的气动力。将常用的二维当量模态数据转换成绕机体瞬时转动中心的模态数据。根据试验提供的旋翼升力和粘弹减摆器的复模量,计算了两种型号直升机地面运转时的稳定性,相对于时域内系统自由度的仿真和模态识别法,本文的处理方式比较简便和快捷。 1 减摆器的力学模型 对直升机地面共振来说,粘弹减摆器的复模量(储能模量及耗能模量)对其动稳定性有很大的影响。目前工程应用时,往往将粘弹减摆器等效成一线性刚度和黏性阻尼,其力学模型为 f=G′x0+G′?x+G″˙?x/ω(1)f=G′x0+G′x?+G′′x?˙/ω(1) 式中G′,G″分别为减摆器的储能模量和耗能模量,?xx?为减摆器的动态位移,ω为振动圆频率(rad/s)。 将减摆器对摆振铰的力矩、对桨叶作用力在当量铰处的分量对机体的力矩分别加到文献的旋翼摆振运动方程、机体运动方程中,即可组成带粘弹减摆器的旋翼/机体耦合运动方程。粘弹减摆器与旋翼连接形式不同,单片桨叶所受的减摆器作用力、减摆器对摆振铰的作用力矩也不同。 1.1 带减摆器的比合方程 以某带星型柔性桨毂的直升机为例,粘弹减摆器与旋翼的连接形式如图1所示。因粘弹减摆器不存在几何耦合,它只承受由减摆器运动引起的剪切作用。设减摆器与摆振铰的距离为Rd,则减摆器对摆振铰的力矩Md为 Μd=Rdf(x0,?x)(2)Md=Rdf(x0,x?)(2) 摆振角ζ与减摆器x之间的关系为 x=Rdζ0+Rd?ζ(3)x=Rdζ0+Rdζ?(3) 直升机地面运转平衡计算时,减摆器只考虑静态位移,即将Md加入桨叶摆振运动方程中,采用牛顿迭代法求解。其中 Μd=R2dG′ζ0(4)Md=R2dG′ζ0(4) 粘弹减摆器的线化小扰动力矩为 δΜd=Rd(G′δx+G″δ˙x/ω)=R2d(G′δζ+G″δ˙ζ/ω)(5)δMd=Rd(G′δx+G′′δx˙/ω)=R2d(G′δζ+G′′δζ˙/ω)(5) 式中δζ=δζ0+δζccosψk+δζssinψk,δ˙ζk=Ω[δ˙ζ0+(δ˙ζc+δζs)cosψ+(δ˙ζs-δζc)sinψk]δζ=δζ0+δζccosψk+δζssinψk,δζ˙k=Ω[δζ˙0+(δζ˙c+δζs)cosψ+(δζ˙s?δζc)sinψk]。 ω可以用旋转桨叶一阶摆振固有圆频率ωζ0来近似代替,并计入减摆器储能模量G′的影响。将减摆器的储能模量G′转换为桨叶摆振面的等效刚度K0=R2d2dG′。 ωζ0=√[Κζ+Κ0Κh/(Κ0+Κh)+Ω2eSb]/Ιb(6)ωζ0=[Kζ+K0Kh/(K0+Kh)+Ω2eSb]/Ib?????????????????????????????√(6) 式中 Kζ为弹性轴承摆振面的约束刚度,Kh为星型柔性支臂摆振面的刚度。 展开式(5),得到带减摆器的旋翼/机体耦合扰动运动方程,各系数不再详细列出。 δΜk=xk1δζ0+xk2δζc+xk3δζs+xc1δ˙ζ0+xc2δ˙ζc+xc3δ˙ζs(7)δMk=xk1δζ0+xk2δζc+xk3δζs+xc1δζ˙0+xc2δζ˙c+xc3δζ˙s(7) 1.2 小扰动方程推导 减摆器采用叶间连接形式时,一片桨叶受到减摆器的作用力矩与三片桨叶的挥/摆/扭运动有关,采用矢量法推导小扰动方程非常困难。在地面运转状态,减摆器的动幅值主要受桨叶摆振运动影响,因此在实际工程应用中,一般不考虑桨叶变距、挥舞运动对减摆器的影响。 1.2.1 减摆器长度值 采用三角形