适于时变幅值分析的直升机黏弹减摆器模型

吴靖， 胡国才， 刘湘一

(海军航空大学 航空基础学院, 烟台 264001)

结构简单、重量较轻、易于维护及可靠性强的黏弹减摆器现广泛应用于无铰式、无轴承式及星型柔性旋翼直升机桨叶减摆。黏弹减摆器的非线性特性对直升机旋翼/机体耦合动稳定性影响较大，而耦合系统的非线性动稳定性分析一般采用时域分析法，因此建立一种能用于直升机旋翼/机体耦合非线性动稳定性时域分析的黏弹减摆器模型尤为重要。

直升机动力学研究中常用的黏弹减摆器模型有VKS(Voigt-Klevin)模型[1-3]、ADF(Anelastic Displacement Fields)模型[4-7]和迟滞非线性模型[8-11]等。在直升机旋翼/机体耦合动稳定性分析中，一般对黏弹减摆器进行线性化处理后采用特征分析法：Gandhi和Chopra[12]对黏弹减摆器线性化处理后，分析了其对直升机悬停空中共振的影响；笔者课题组[13]采用线性化的VKS模型分析了黏弹减摆器对直升机前飞空中共振的影响；王波等[14]采用线性化的VKS模型分析了黏弹减摆器的不同连接方式对直升机地面共振的影响。而在保留黏弹减摆器非线性情况下进行直升机旋翼/机体耦合动稳定性分析的研究较少，王波等[15]将基于复模量的非线性VKS模型用于直升机空中共振时域分析；卫丽君等[16]则将该方法用于分析黏弹减摆器不同连接方式对直升机空中共振的影响。但基于复模量识别的非线性VKS模型不能很好地预估时域响应，而基于时域响应识别的VKS模型却难以描述其幅变特性[2]。

为了适用较宽的动幅值范围，黏弹减摆器一般引入动幅值变量[11]，但很难将其用于幅值变化的直升机旋翼/机体耦合动稳定性时域分析。为解决上述问题，对动幅值参量进行改进，并将改进后的黏弹减摆器模型用于直升机地面共振非线性时域分析，验证上述方法的可行性。

1 黏弹减摆器改进模型

现有的黏弹减摆器模型为了适用于宽的动幅值范围，一般在模型中引入动幅值参量δ，如式(1)所示一种双线性迟滞模型，在周期激振x=x0+δcos(ωt)作用下(x0为静位移,ω为响应频率，t为时间)，黏弹减摆器恢复力表示为

(1)

其中：xs为滑移极限位移；ks=zs/xs为滑移线性刚度，zs为滑移极限恢复力。

黏弹减摆器的复模量经推导可得

(2)

式中：G′和G″分别为储能模量和耗能模量。

在直升机旋翼/机体耦合动稳定性时域分析中，减摆器响应的动幅值一般是变化的，而减摆器刚度会随动幅值变化而变化，这就导致频率也是变化的，采用类似上述包含动幅值及频率参量的模型难以进行时域分析，需要进行改进。

1.1 单频激振条件

减摆器在收敛、中性稳定和发散3 种情况下受到扰动后的时域响应如图1所示。

减摆器等效线性化后，其响应可表示为

x=x0+δ0e-ξωntsin(ωdt+φd)

(3)

(4)

其中：Ib和Sb分别为桨叶对摆振铰的惯性矩和静矩；kb为除减摆器外的摆振刚度；e为摆振铰外伸量；Ω为旋翼转速；Rd为摆振铰至减摆器距离。

由式(3)可知，响应振动幅值可表示为

(5)

摆振阻尼较大时，响应很快衰减，没有振荡特性，故暂针对小阻尼情况进行研究。此时，响应振动幅值可近似表示为

(6)

图1中所示的点线即为采用式(6)计算的振动幅值，结果说明了式(6)计算的准确性。

将式(6)代入式(1)即可消除黏弹减摆器模型中的动幅值及频率参量。

1.2 双频激振条件

直升机前飞或滑跑时，桨叶摆振面由于存在1Ω的哥氏力矩作用，将产生1Ω的强迫振动响应作用于减摆器。因此，黏弹减摆器受到扰动后实际上是处于双频激振条件下，其动态响应可表示为

(7)

式中：xΩ为强迫振动位移；δΩ为强迫振动幅值；φΩ为强迫振动初始相位角。

1Ω强迫振动在摆振阻尼较小的情况下，强迫振动幅值可表示为

(8)

式中：Mc为哥氏力矩幅值。

此时，扰动响应振动幅值δd可表示为

(9)

完整响应振动幅值可表示为

(10)

式(10)中频率是作用于速度项，因此将其表示成关于强迫振动和扰动运动速度幅值比的指数函数。当δΩ=0时，即没有强迫振动的情况下，式(10)与式(6)相同，表示单频条件下的振动幅值。

在1.1节单频激振条件的基础上叠加频率为6 Hz、幅值为1 mm的强迫振动，3种情况下的响应如图2所示，图中实线表示位移，采用式(10)计算的幅值如点线所示。

由图2可以看出，式(10)所计算的幅值曲线较好地反映了双频条件下黏弹减摆器响应幅值在时域上的变化趋势。

因为幅值计算需要频率参数，而频率又受幅值的影响，因此初始幅值用初始位移与静位移差值的绝对值|xt=0-x0|代替。

2 计入黏弹减摆器的地面共振分析

2.1 计入黏弹减摆器的地面共振模型

黏弹减摆器安装在旋翼上主要用于提高桨叶摆振阻尼，计入黏弹减摆器作用后的直升机地面共振当量平面模型为

(11)

2.2 仿真分析主要参数

某型黏弹减摆器经识别后所得的参数为：a1=-0.346 2，a3=-5.087，a5=-57.83，p=0.798 2 mm，k1=307.6 N/mm，k3=11.76 N/mm3，k5=-0.046 8 N/mm5，b=3.103，q=1.367 mm，ce=541.8 N/mm，zs=505.6δN，xs=0.234 8δmm。

采用式(2)对静位移为6 mm时，不同动幅值下黏弹减摆器的复模量进行计算，并与实验值对比，结果如图4所示。

由图4可知，采用本文黏弹减摆器双线性迟滞模型计算的复模量与实验值吻合较好，证明了该模型的准确性。后续也将采用静位移为6 mm状态下的黏弹减摆器特性进行计算。

某型直升机旋翼及机体的主要参数如表1所示，除减摆器外的摆振刚度及阻尼较小，可以忽略不计，即kb=cb=0。

参 数数 值桨叶片数4桨叶质量/kg42.3桨叶对摆振铰静矩/(kg·m)123.7桨叶对摆振铰惯性矩/(kg·m2)457摆振铰外伸量/m0.23减摆器到摆振铰距离/m0.35机体纵向当量质量/kg2000机体纵向固有频率/Hz1机体纵向当量阻尼/(N·s·m-1)1000机体横向当量质量/kg800机体横向固有频率/Hz1.5机体横向当量阻尼/(N·s·m-1)500

2.3 地面共振仿真分析

旋翼摆振后退型模态在直升机地面共振分析中较为重要，为分析旋翼摆振后退型模态，首先需要通过对桨叶摆振方向以摆振固有频率施加激振力矩的方式激出旋翼的摆振后退型模态。由图4可知，黏弹减摆器在不同动幅值下具有不同的复模量，因此采用的激振力矩需要使得激出的桨叶摆振幅值(减摆器振动幅值)与激振频率(减摆器储能模量)符合图 4中的复模量曲线。根据强迫振动理论[17]可估计需要的激振力矩幅值Mr为

(12)

黏弹减摆器取不同动幅值下的复模量确定激振频率，采用式(12)确定激振力矩幅值，激出不同转速下旋翼的摆振后退型模态后，确定模拟出的减摆器响应幅值如图5所示。

由图5可知，模拟出的减摆器响应幅值δr与最初选择的振动幅值(理论值)基本吻合，虽然个别点由于旋翼摆振后退型模态与机体运动模态耦合引起的误差稍大，但最大误差仍不超过6%，证实了采用式(12)确定激振力矩幅值的可行性。

选取动幅值分别为0.5、2和6 mm的3种情况下对应的复模量状态，针对黏弹减摆器线性化(复模量取常值)的情况，在激出旋翼摆振后退型模态后，采用基于多桨叶坐标变换的阻尼识别法[18]对瞬态响应进行模态阻尼识别，并与采用特征值法计算的模态阻尼进行对比，具体结果如图6所示，图中BX和BY分别表示机体的纵横向运动模态，RL表示旋翼摆振后退型模态，fm和σ分别为模态频率和模态阻尼。

由图6可知，在黏弹减摆器取不同复模量的情况下，通过阻尼识别所得的模态阻尼与采用特征值法计算所得的模态阻尼基本吻合，说明所采用的基于多桨叶坐标变换的阻尼识别法的有效性。由图中还可以看出，在旋翼摆振后退型模态与机体模态耦合时，激出旋翼摆振后退型响应的同时也激出了机体响应，使得所识别的摆振后退型模态阻尼误差稍大。

以旋翼转速Ω=300 r/min为例，取黏弹减摆器响应幅值δ=6 mm时的储能模量确定激振频率，以此激振频率激出系统的摆振后退型模态后，计算黏弹减摆器线性化(复模量取常值)和非线性2种情况下的瞬态响应，由于各桨叶减摆器响应仅存在相位差的区别，因此只给出第1个减摆器的位移时域响应和功量图，如图7、图8所示。

由图7可以看出，与黏弹减摆器线性化的响应相比，在保留减摆器非线性的情况下，其响应衰减得更快，且振动频率更大。据分析可知，是衰减过程中减摆器阻尼及刚度随响应幅值的减小而增加引起的。

由图8可见，保留黏弹减摆器非线性的功量图与线性化后差异较大，说明在进行黏弹减摆器载荷计算时不能忽略其非线性的影响。

旋翼摆振后退型模态响应ζRL如图9所示，采用阻尼识别法对旋翼摆振后退型模态阻尼进行识别，其随时间的变化如图10所示。

由图10可知，由于减摆器非线性的影响，在响应衰减的过程中，随着黏弹减摆器振动幅值的减小，其耗能模量增加，使得桨叶的摆振阻尼增加，旋翼摆振后退型模态阻尼在时域上呈增加趋势，减摆器线性化后，摆振后退型模态阻尼保持不变，这也就是图9中保留减摆器非线性后，旋翼摆振后退型响应衰减的更快的原因。

据分析，在图4所示的黏弹减摆器非线性条件下，系统初始稳定能确保整个过程的稳定。结合图6的结果可知，减摆器受到的扰动越大，旋翼摆振后退型模态阻尼越小，情况越危险。因此，在地面共振试验或理论分析中应该重点关注大扰动下系统初始阶段的响应及其模态阻尼。

3 结 论

1) 针对小摆振阻尼比的情况，给出了单频及双频条件下黏弹减摆器模型中动幅值参量的计算方法，采用该方法计算系统在收敛、中性稳定及发散3种情况下的幅值曲线能较好地反映响应幅值在时域上的变化趋势，说明该方法可使带动幅值参量的黏弹减摆器模型用于直升机旋翼/机体耦合动稳定性时域分析。

2) 给出了计入非线性黏弹减摆器后，激出旋翼摆振后退型模态所需激振力矩的计算方法，在不同转速不同复模量状态下，采用该方法确定的激振力矩对桨叶进行激振激出的响应幅值与预期值误差不超过6%，说明了所述方法的准确性。

3) 将改进后的黏弹减摆器模型用于直升机地面共振时域分析，对摆振后退型响应进行分析可知，系统稳定时，与线性化结果相比，计入黏弹减摆器非线性后，旋翼摆振后退型响应衰减更快，其模态阻尼在时域上呈增加趋势。