李貞坤，程起有，朱艷，錢峰，劉晨，代志雄

（中國(guó)直升機(jī)設(shè)計(jì)研究所 直升機(jī)旋翼動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室，江西 景德鎮(zhèn) 333001）

旋轉(zhuǎn)葉片廣泛應(yīng)用于工業(yè)領(lǐng)域，如風(fēng)力發(fā)電機(jī)葉片、飛機(jī)發(fā)動(dòng)機(jī)葉片、螺旋槳葉片以及直升機(jī)槳葉等。由于共振和顫動(dòng)引起的槳葉失效問題是槳葉故障的主要原因之一，因此深入地了解旋轉(zhuǎn)槳葉的動(dòng)力學(xué)行為有助于減少槳葉的振動(dòng)故障。

細(xì)長(zhǎng)的槳葉可視為旋轉(zhuǎn)的懸臂梁，基于線彈性小變形理論，研究者建立了諸多動(dòng)力學(xué)模型來研究旋轉(zhuǎn)梁的線彈性動(dòng)力學(xué)行為，如旋轉(zhuǎn)槳葉的共振頻率與模態(tài)等問題。其后，研究者考慮了更復(fù)雜的情況，如旋轉(zhuǎn)梁處于熱環(huán)境、采用復(fù)合材料以及考慮裂紋等，以研究旋轉(zhuǎn)梁的動(dòng)特性。Qin等研究了旋轉(zhuǎn)復(fù)合材料薄壁梁的自由振動(dòng)特性，考慮了氣動(dòng)力和濕熱環(huán)境的影響。Xie等采用有限元對(duì)含裂紋的旋轉(zhuǎn)梁進(jìn)行了研究，考慮了離心效應(yīng)與裂紋效應(yīng)的耦合。Liang等研究了彎扭耦合旋轉(zhuǎn)梁，討論了控制方程中螺旋項(xiàng)和離心項(xiàng)的作用。Guo等基于鐵木辛柯梁理論，研究了含預(yù)安裝角和扭轉(zhuǎn)角的旋轉(zhuǎn)漸變梁6個(gè)位移分量之間的耦合效應(yīng)。

高速旋轉(zhuǎn)的槳葉葉尖周圍的流場(chǎng)非常不穩(wěn)定，在強(qiáng)烈的周期性氣流擾動(dòng)下可能產(chǎn)生大振幅的非線性振動(dòng)，而當(dāng)前基于線彈性小變形理論獲得的結(jié)果和結(jié)論不足以解釋觀測(cè)到的非線性問題。因此，研究者對(duì)旋轉(zhuǎn)槳葉的非線性動(dòng)力學(xué)行為（如主共振、次共振等）進(jìn)行了研究，同時(shí)對(duì)功能梯度材料、溫度分布、粘彈性等影響開展了討論。Chen等研究了有預(yù)變形旋轉(zhuǎn)梁的3∶1內(nèi)共振，預(yù)變形來自熱梯度分布，考慮了二次和三次非線性，采用多尺度法求解，并用數(shù)值法進(jìn)行了驗(yàn)證。Thomas等考慮了不同非線性來源導(dǎo)致的幾何非線性，分別建立了3種模型，采用解析和有限元的方法求解，研究了大撓度下的跳躍行為，以及不同轉(zhuǎn)速下的軟化/硬化特性。Emam等對(duì)附著有質(zhì)量-彈簧系統(tǒng)的旋轉(zhuǎn)梁亞結(jié)構(gòu)的主共振進(jìn)行了研究，采用攝動(dòng)法求解，討論了系統(tǒng)在自由振動(dòng)和受迫振動(dòng)下的響應(yīng)。Yao等研究了變轉(zhuǎn)速旋轉(zhuǎn)梁的動(dòng)特性，考慮梁的扭轉(zhuǎn)和高溫高速氣流等的影響，采用數(shù)值方法分析了葉片的主共振及1∶1內(nèi)共振。

在上述研究中，研究人員研究了旋轉(zhuǎn)梁的主共振現(xiàn)象，即激勵(lì)頻率接近葉片的固有頻率。非線性系統(tǒng)中存在更復(fù)雜的共振機(jī)制，其他不接近共振頻率的激勵(lì)頻率也可能導(dǎo)致共振，如超諧波共振、亞諧波共振和組合諧波共振。以組合諧波共振為例，當(dāng)槳葉受2個(gè)及以上頻率激勵(lì)時(shí)，如果2個(gè)激勵(lì)頻率滿足特定的關(guān)系，即使2個(gè)頻率均遠(yuǎn)離槳葉共振頻率，槳葉同樣會(huì)發(fā)生共振。Inoue等在實(shí)驗(yàn)中證實(shí)了風(fēng)輪機(jī)葉片存在超諧波共振。Shahgholi等利用諧波平衡法研究了旋轉(zhuǎn)軸的組合共振和次諧波共振。Zhang等探索了預(yù)扭轉(zhuǎn)旋轉(zhuǎn)梁的主共振、超諧波和組合諧波共振等情形，采用多尺度法求解，并使用Runge-Kutta法對(duì)結(jié)果進(jìn)行了驗(yàn)證。熊春等采用多尺度法開展了旋轉(zhuǎn)Rayleigh梁在軸向力下的參數(shù)穩(wěn)定性分析，研究了超諧波共振和組合共振。伊海銘等針對(duì)雙轉(zhuǎn)子系統(tǒng)，采用四階Runge-Kutta法分析了分岔區(qū)域的組合共振特性。

由以上文獻(xiàn)總結(jié)可知，文獻(xiàn)中對(duì)槳葉的非線性主共振研究較多，次共振特別是組合共振研究較少，而針對(duì)直升機(jī)槳葉的組合共振研究更少。本文對(duì)直升機(jī)槳葉的組合共振進(jìn)行了研究，基于哈密頓原理推導(dǎo)出含幾何非線性與非線性慣性的控制方程。采用伽遼金法離散，并使用攝動(dòng)法中L-P法（Lindstedt-Poincaré method）求解，分析了組合共振發(fā)生條件及相應(yīng)的穩(wěn)態(tài)幅頻響應(yīng)，討論了超/次諧波共振與組合共振共存的可能，并對(duì)組合共振進(jìn)行了參數(shù)分析，探究了激勵(lì)力和阻尼等對(duì)共振的影響。

1 建模與求解

1.1 模型建立與無量綱

如圖1所示，以Bo105槳葉為研究對(duì)象，將槳葉視為旋轉(zhuǎn)的細(xì)長(zhǎng)懸臂梁，繞軸以角速度旋轉(zhuǎn)，槳葉長(zhǎng)度為，揮舞方向（方向）上受2個(gè)周 期變化分布的力、作用，激勵(lì)力頻率分別為、。槳葉軸向和橫向位移分別為(,)、(,)，和分別表示橫坐標(biāo)和時(shí)間，、為作用力離固定端的距離。

圖1 槳葉受2個(gè)不同激勵(lì)頻率的力作用示意及坐標(biāo)關(guān)系 Fig.1 Schematic diagram and coordinate relationship of the blades subjected to two different excitation frequencies

槳葉應(yīng)變能可寫為：

動(dòng)能和外力功分別為：

根據(jù)哈密爾頓原理，對(duì)總能量求變分：

將公式（1）—（3）代入方程（4）進(jìn)行相應(yīng)變換，并保留到三次項(xiàng)，可得旋轉(zhuǎn)梁運(yùn)動(dòng)方程為：

式中：為密度；為梁彎曲剛度，其中為彈性模量，為截面慣性矩；為阻尼系數(shù)；為Dirac函數(shù)。等式（5）左邊與時(shí)間相關(guān)的三次項(xiàng)為非線性慣性，產(chǎn)生軟化效應(yīng)（Softening Effect），其他非線性項(xiàng)則為幾何大變形引起，產(chǎn)生硬化效應(yīng)（Hardening Effect），含項(xiàng)為旋轉(zhuǎn)產(chǎn)生。如果令轉(zhuǎn)速為0，方程退化成經(jīng)典的懸臂梁非線性振動(dòng)控制方程，見式（6）。

使用合適的參數(shù)，式（5）可無量綱化為：

采用的無量綱參數(shù)為：

其中是懸臂梁一階振動(dòng)頻率：

偏微分方程（7）可采用伽遼金法離散成二階常微分方程：

式中：φ()為懸臂梁各階模態(tài)函數(shù)；q()為對(duì)應(yīng)的廣義時(shí)間坐標(biāo)。將式（10）代入式（7），并分別乘以對(duì)應(yīng)的各階模態(tài)函數(shù)φ()，對(duì)從0到1積分可得離散后的方程為（取=1）：

式中：上標(biāo)點(diǎn)號(hào)(·)和撇號(hào)(')分別代表對(duì)時(shí)間和坐標(biāo)的導(dǎo)數(shù)。

1.2 攝動(dòng)法求解

在控制方程已離散的基礎(chǔ)上，采用攝動(dòng)法中的L-P法求解上述常微分方程。根據(jù)L-P法，式（11）可以寫為：

取新的時(shí)間變量=，并將激勵(lì)頻率ω和廣義時(shí)間坐標(biāo)q分別展開為的級(jí)數(shù)，見式（13）。

則關(guān)于時(shí)間的導(dǎo)數(shù)改寫為：

2個(gè)外激勵(lì)頻率改寫為：

將式（13）—（15）代入方程（12），并令方程兩端同次冪系數(shù)相等，得：

方程（16）的解為：

其中，解的第1項(xiàng)cos （+）表示自由振動(dòng)，響應(yīng)頻率為旋轉(zhuǎn)梁一階共振頻率，第2、3項(xiàng)表示受迫振動(dòng)，響應(yīng)頻率與激勵(lì)頻率一致。

將方程的解（18）代入方程（17），提取同時(shí)含有,的項(xiàng)可得：

由式（19）可知，當(dāng),滿足一定的關(guān)系，有可能發(fā)生組合共振，即：

以2+=為例，此時(shí)2+=，為當(dāng)前轉(zhuǎn)速下懸臂梁共振頻率，有：

在式（19）中消除永久項(xiàng)，得幅頻響應(yīng)曲線為：

其中：

1.3 同時(shí)存在組合共振與超/次諧波共振

回到公式（19），只保留了同時(shí)包含和的項(xiàng)，其余的項(xiàng)被舍棄，以研究組合共振發(fā)生條件。在滿足組合共振條件下，如果、還滿足額外條件，則可能在組合共振產(chǎn)生的同時(shí)發(fā)生超諧波或次諧波共振。

將方程的解（18）代入方程（17），在所得的方程中提取包含或的項(xiàng)可得：

當(dāng)滿足式（20）任一條件時(shí)，會(huì)發(fā)生組合共振。若此時(shí)、還滿足以下條件之一，則會(huì)同時(shí)發(fā)生超諧波共振：

類似地，組合共振與次諧波共振會(huì)同時(shí)發(fā)生，若、還滿足：

以同時(shí)存在組合共振與超諧波共振為例，使用公式（20）和式（25），即：

代入式（24）中消除永久項(xiàng)可得幅頻響應(yīng)方程為：

1.4 計(jì)算方法驗(yàn)證

為了驗(yàn)證本計(jì)算方法的有效性，采用L-P法對(duì)文獻(xiàn)[24]中懸臂梁非線性振動(dòng)控制方程進(jìn)行求解，并與文獻(xiàn)中結(jié)果進(jìn)行比較。與文獻(xiàn)一樣采用一階模態(tài)離散，離散后控制方程為：

與1.3節(jié)類似，將激勵(lì)頻率ω和廣義時(shí)間坐標(biāo)q展開為的級(jí)數(shù)，并代入式（29），令方程兩端同次冪系數(shù)相等得：

方程（30）的解為：=cos（+），代入方程（31）消除永久項(xiàng)可得幅頻響應(yīng)方程為：

式（32）與文獻(xiàn)中多尺度法結(jié)果一致，證明本文采用計(jì)算方法的正確性。

2 結(jié)果與討論

2.1 組合共振

根據(jù)式（22）可作無量綱幅頻響應(yīng)曲線（如圖2所示），采用的無量綱參數(shù)為：==0.05，=1.5×10。與線彈性振動(dòng)相比，非線性振動(dòng)幅頻響應(yīng)曲線與線彈性最顯著的區(qū)別為在共振點(diǎn)附近曲線左右不對(duì)稱。由圖2可知，幅頻響應(yīng)曲線展現(xiàn)出向左偏的趨勢(shì)，即表現(xiàn)出軟化特性（Softening Effect），且幅頻曲線的頻率與振幅不再是一一對(duì)應(yīng)，而是表現(xiàn)出多值性，以及跳躍現(xiàn)象。當(dāng)激勵(lì)頻率由小增大，曲線沿著左下半支前行，在虛線處跳躍到右上半支；當(dāng)激勵(lì)頻率由大減小，曲線沿著右半支前行，到頂點(diǎn)后跳躍到左下半支。

圖2 無量綱幅頻響應(yīng)曲線 Fig.2 Dimensionless amplitude frequency resonance response curve

由式（18）可知，非線性振動(dòng)與線彈性振動(dòng)的另一區(qū)別在于槳葉的響應(yīng)是由2部分組成：自由振動(dòng)部分cos（+），頻率為旋轉(zhuǎn)懸臂梁共振頻率；受迫振動(dòng)部分cos（+）+cos （+），響應(yīng)頻率分別對(duì)應(yīng)兩激勵(lì)力頻率和，且發(fā)生組合共振時(shí)存在關(guān)系2+=，令：

取共振時(shí)槳葉的時(shí)間歷程進(jìn)行分析，如圖3所示。圖3a、b為受迫振動(dòng)部分2個(gè)激勵(lì)的響應(yīng)，兩者振幅相同，==0.05，相位不同。圖3c為自由振動(dòng)部分時(shí)間歷程，無量綱振幅=0.11，大于兩激勵(lì)力振幅，且略大于兩者的合振幅。圖3d為旋轉(zhuǎn)梁的合響應(yīng)，振幅為三者的疊加，同時(shí)放上組合共振響應(yīng)作對(duì)比，可以看出，組合共振振幅占總振幅的1/2。

圖3 無量綱時(shí)間歷程 Fig.3 Dimensionless time histories at the resonance point: a) excitation force f1; b) excitation force f2; c) the free vibration part; d) the total response (solid line) and combined resonance response (dashed line)

2.2 參數(shù)分析

圖4與圖5分別展示了無量綱自由振動(dòng)頻響曲線和最大響應(yīng)振幅隨激勵(lì)力、的變化。圖4a和圖5a表明，隨著、增大，兩無量綱響應(yīng)振幅均增大，此為兩者共同點(diǎn)。圖4b和圖5b表明，自由振動(dòng)最大振幅隨受迫激勵(lì)力、增加的程度不一。圖4b中，增大，振幅類似于冪指數(shù)增加；增大，振幅成比例增加。

圖4 不同激勵(lì)力f1下的無量綱頻響曲線和最大振幅 Fig.4 Dimensionless frequency response curve (a) and the maximum amplitude (b) under different excitation force f1

圖5 不同激勵(lì)力f2下的無量綱頻響曲線和最大振幅 Fig.5 Dimensionless frequency response curve (a) and the maximum amplitude (b) under different excitation force f2

作最大振幅與激勵(lì)力大小的比值圖，如圖6所示?？梢钥闯?，/呈線性關(guān)系，/為常數(shù)，表明組合共振振幅與激振力呈線性關(guān)系，卻與呈平方關(guān)系，即在當(dāng)前受力形勢(shì)下（激勵(lì)頻率）對(duì)組合共振的振幅影響更大。

圖6 Amax/ fi隨fi變化（i=1,2） Fig.6 Ratio of Amax/ fi with different fi (i=1,2)

激勵(lì)力位置對(duì)響應(yīng)的作用如圖7所示。無量綱參數(shù)和表示激勵(lì)力、與固定端的距離。由圖7 可知，隨著的增加，即激勵(lì)力遠(yuǎn)離固定端，曲線振幅相應(yīng)增加；的影響趨勢(shì)與一致。仔細(xì)對(duì)比圖7a、b可知，響應(yīng)的振幅隨變化更劇烈。這與前文觀察的現(xiàn)象一致，與相關(guān)聯(lián)，響應(yīng)振幅與成正比，而僅與成正比，因此響應(yīng)隨的變化更大。

圖7 響應(yīng)隨激勵(lì)力位置的變化 Fig.7 Changes of position of the excitation force on the response

如前文所述，組合共振的存在極大地提高了旋轉(zhuǎn)梁橫向運(yùn)動(dòng)的振幅。對(duì)比圖3c、d可知，組合共振使得振幅提高了1倍，極大地增加了槳葉的不穩(wěn)定性和破壞的可能。不同阻尼下自由振動(dòng)部分的無量綱頻響曲線如圖8所示。由圖8b可知，增大阻尼能迅速地減小自由振動(dòng)部分振幅，如無量綱阻尼從3×10增大到6×10，振幅降低了50%。因此，增大阻尼能有效抑制組合共振強(qiáng)度，防止旋轉(zhuǎn)梁振幅過大而導(dǎo)致潛在的破壞。

圖8 不同阻尼下的無量綱頻響曲線和最大振幅 Fig.8 Dimensionless frequency response curve (a) and maximum amplitude (b) under different damping;

2.3 超諧波共振與組合共振共存

圖9 有超諧波共振存在的組合共振（CR+SR）與 普通組合共振（CR）對(duì)比 Fig.9 Comparison of combination resonance (CR+SR) with the coexistence of super-harmonic resonance and commom combination resonance (CR)

3 結(jié)論

本文研究了旋轉(zhuǎn)梁橫向運(yùn)動(dòng)下受到2個(gè)激勵(lì)力的組合共振，采用Hamilton原理推導(dǎo)了旋轉(zhuǎn)梁控制方程，同時(shí)考慮了幾何非線性與非線性慣性。采用L-P法對(duì)控制方程進(jìn)行求解，獲得了組合共振發(fā)生條件和共振時(shí)的幅頻響應(yīng)曲線方程，對(duì)組合方程進(jìn)行參數(shù)分析，并進(jìn)一步研究同時(shí)存在超/次諧波共振與組合共振的條件，研究獲得如下結(jié)論：

1）當(dāng)兩激勵(lì)頻率滿足2±=、2±=或1/2(±)=條件之一時(shí)，旋轉(zhuǎn)梁可能發(fā)生組合共振，此時(shí)槳葉的響應(yīng)既有2個(gè)激勵(lì)頻率成分，又有自由振動(dòng)成分。

2）以2±=條件為例進(jìn)行的分析表明，旋轉(zhuǎn)梁表現(xiàn)出軟化特性，組合共振中自由振動(dòng)部分占主導(dǎo)，其振幅與兩受迫振動(dòng)振幅疊加相當(dāng)。

3）參數(shù)分析表明，在2±=條件下，激振力的大小和位置對(duì)組合共振影響更大，振幅與呈平方關(guān)系，與呈線性關(guān)系。

4）在條件1）滿足的情況下，若=1/3或=1/3，會(huì)同時(shí)存在超諧波共振；若=或=，會(huì)同時(shí)存在次諧波共振。=1/3的結(jié)果表明，同時(shí)存在超諧波共振時(shí)，會(huì)進(jìn)一步增加組合共振的振幅。