梁 峰，包日東，金 瑩，蘇 勇

（沈陽化工大學(xué)能源與動力工程學(xué)院，沈陽 110142）

外激勵作用下輸流管道伴隨內(nèi)共振的非線性振動分析

梁 峰，包日東，金 瑩，蘇 勇

（沈陽化工大學(xué)能源與動力工程學(xué)院，沈陽 110142）

利用多元L－P法研究外部周期激勵下兩端固定輸流管道伴隨內(nèi)共振的非線性受迫振動問題。外激勵流固耦合系統(tǒng)固有頻率第二階約為第一階3倍且激勵頻率接近前兩階固有頻率中間值時會發(fā)生伴隨強(qiáng)烈內(nèi)部共振的組合共振，并用多元L－P法求解振動響應(yīng)，分析前兩模態(tài)運動及外激勵幅值對內(nèi)共振的影響。數(shù)值算例揭示出系統(tǒng)因內(nèi)共振發(fā)生的更豐富、復(fù)雜的動力學(xué)行為，隨激勵幅值增大內(nèi)共振發(fā)生趨勢降低，響應(yīng)形式亦發(fā)生變化。用多元L－P法研究非線性動力學(xué)便捷、高效。

輸流管道；非線性振動；內(nèi)共振；外周期激勵；多元L－P法

輸流管道廣泛應(yīng)用于航空航天、水利工程、石油工業(yè)、農(nóng)業(yè)及日常生活中，其振動穩(wěn)定性問題一直被關(guān)注。（作為典型流固耦合系統(tǒng)，輸流管道振動失穩(wěn)通常由與內(nèi)部流體耦合所致。對此，利用理論、實驗方法進(jìn)行的研究發(fā)現(xiàn)包括顫振、屈曲、參數(shù)共振、分岔及混沌等多種動力學(xué)現(xiàn)象［1－3］。諸多工程結(jié)構(gòu)中輸流管道不僅受內(nèi)流作用亦常受外部激勵作用，如橫向風(fēng)載、地震載荷及支承運動時產(chǎn)生的激勵（該載荷通過適當(dāng)變換可簡化為周期載荷）等，雖激勵幅值較小，但系統(tǒng)頻率關(guān)系滿足某些條件，如第二階固有頻率接近第一階的3倍，此時管道系統(tǒng)則會發(fā)生強(qiáng)烈內(nèi)部共振。關(guān)于動力學(xué)系統(tǒng)內(nèi)共振問題，Wang等［4－7］分別利用懸索、梁、板及殼模型進(jìn)行研究，表明內(nèi)共振不僅存在，且其對系統(tǒng)非線性動力學(xué)行為影響明顯。Xu等［8－9］理論上解釋內(nèi)流引起管道內(nèi)共振機(jī)理。Panda等［10－11］用多尺度法研究脈動內(nèi)流引起具有內(nèi)共振的主參數(shù)共振及組合參數(shù)共振，認(rèn)為出現(xiàn)內(nèi)共振，系統(tǒng)動力學(xué)行為更復(fù)雜。以上研究均未考慮外部激勵。而海洋立管、海底管道、橋梁懸跨管道等，對風(fēng)、流引起的橫向振動尤其渦激振動［12－13］、端部激勵引起的參數(shù)激勵振動、考慮內(nèi)流及外部激勵管道的非線性振動研究均取得較多成果。本文對內(nèi)流、外部周期激勵聯(lián)合作用管道出現(xiàn)內(nèi)共振的非線性動力學(xué)行為進(jìn)行深入研究。

組合共振為工業(yè)管道常發(fā)生但得不到重視的非線性振動。由于發(fā)生條件較苛刻，且其振動形態(tài)較復(fù)雜，對其研究較少。事實上，其振動危害并不亞于主共振及其它類共振。尤其伴隨內(nèi)共振發(fā)生時其振幅較大、振型較復(fù)雜，為管道振動疲勞破裂的重要原因之一。Panda等［11，14－15］曾對輸流管道組合共振進(jìn)行研究，分別利用平均法、數(shù)值仿真技術(shù)研究脈動流管道組合參數(shù)共振穩(wěn)定性及響應(yīng)，認(rèn)為在組合共振區(qū)域內(nèi)亦會發(fā)生準(zhǔn)周期、組合周期等多種運動，并證明組合共振的復(fù)雜性。以上研究只考慮內(nèi)流激勵，本文在此基礎(chǔ)上增加外部激勵，重點研究和式組合共振（激勵頻率約為前兩階固有頻率和的一半）響應(yīng)問題。本文研究所用多元L－P法已廣泛用于非線性系統(tǒng)的動力學(xué)特性分析。如Lau等［16－17］利用該法分析動力學(xué)系統(tǒng)非線性振動問題；Yang等［18］用改進(jìn)的L－P法（MLP）研究參數(shù)激勵下強(qiáng)非線性Van der Pol-Duffing型振蕩器主共振響應(yīng)；Pu?enjak［19］用拓展的L－P法（ELP）分析具有三次非線性的動力學(xué)系統(tǒng)非穩(wěn)態(tài)振動。本文用該方法對外部周期激勵下兩端固定輸流管道伴隨內(nèi)共振的組合共振響應(yīng)進(jìn)行研究，獲得幅－頻響應(yīng)曲線，并利用增量諧波平衡法進(jìn)行驗證，分析該振動的響應(yīng)特性，討論外激勵幅值對內(nèi)共振的影響。

1 力學(xué)模型及控制方程

圖1為外部簡諧激勵下兩端固定輸流管道力學(xué)模型。實際上無論外部周期力或支承簡諧運動，經(jīng)一定變換后均可簡化為簡諧激勵。

圖1 外部激勵下輸流管道模型Fig．1 Model of a fluid-conveying pipe under external excitation

設(shè)管道軸線為x軸，管道只發(fā)生橫向面內(nèi)振動y（x，t），U為管內(nèi)定常流速，F(xiàn)為外激勵單位長度幅值，Φ為頻率。管道單元、流體單元受力可用牛頓法推導(dǎo)出其橫向振動控制方程［15］為

式中：EI，A1，N，L，m分別為管道彎曲剛度、管壁橫截面積、預(yù)緊力、長度及單位長度質(zhì)量；M，A分別為管內(nèi)流體單位長度質(zhì)量及通流截面積。

式中：（ ）′，（·）分別表示?（ ）／?ξ及?（ ）／?τ。

由式（3）及研究結(jié)果［15，20－21］知，用2階Galerkin法離散管道系統(tǒng)精度已足夠。本文取

式中：ψi（ξ），qi（τ）（i＝1，2）為兩端固定梁前兩階特征函數(shù)及廣義坐標(biāo)。

輸流管道為細(xì)長結(jié)構(gòu)，其外部激勵通常非同相位（因激勵幅值較小，忽略非同相位激勵引起的波長效應(yīng)），本文考慮二維激振力、模態(tài)cosφ及sinφ正交，故其相位差為π／2；設(shè)激勵幅值列陣，外部激勵F cos（Φt）通過式（2）無量綱后為［f1cosφ，f2sinφ］T，將式（4）代入式（3），用振型正交性可得二元二階常微分方程組為

2 用多元L－P法求解

為用多元L－P攝動法，令式（5）引入攝動小參數(shù)ε，變?yōu)?/p>

式中：ωn為獨立非線性響應(yīng)頻率；τn為多個不同時間變量。則qn由單個時間變量τ的函數(shù)變?yōu)槎鄠€獨立時間變量τn的函數(shù)（本文n＝2）。將qn（τ1，τ2，ε）與ωn（ε）展開為ε的冪級數(shù)為

式中：ω10，ω20為系統(tǒng)前兩階固有頻率；ωnk（k＝1，2，…）為待定參數(shù)。

qn對無量綱時間變量τ的一、二階導(dǎo)數(shù)為

由式（12）得q10，q20代入式（13），消去長期項獲得可解性條件，從而獲得外激勵頻率與響應(yīng)振幅之關(guān)系。

考慮外激勵頻率Ω在前兩固有頻率中間值附近及內(nèi)共振，有

為消去長期項，應(yīng)使B11，B12，B21，B22全部為零，但由此所得方程組會因不相容無法求解。因此設(shè)q11，q21特解為

將式（20）代入式（19），比較cosτ1，cosτ2，sinτ1，sinτ2各項系數(shù)，可得以E11，E12，E21，E22為未知量的方程組，其非零解條件為

判別式值可出現(xiàn)幾種根，即①若Δ＞0，有1個實根與1對共軛虛根；②Δ＝0，且P＝Q＝0，有3個零根；③Δ＝0，且P≠0，Q≠0，有3個實根，其中兩個相等；④Δ＜0，有3個不等實根。可由實根個數(shù)分析系統(tǒng)的內(nèi)共振。

3 IHB法驗證

為驗證多元L－P法的有效性，本文用半解析、半數(shù)值法－增量諧波平衡法（Incremental Harmonic Balancemethod，IHB法）重新求解管道系統(tǒng)受迫振動響應(yīng)。IHB法由Lau等［22］提出，并在具有三次非線性及陀螺特性等強(qiáng)非線性系統(tǒng)中廣泛應(yīng)用。該法將數(shù)值分析中的增量法與傳統(tǒng)諧波平衡法有機(jī)結(jié)合，且對所研究系統(tǒng)非線性的強(qiáng)、弱無限制，為研究非線性振動尤其強(qiáng)非線性振動的有效方法。梁峰等［3］曾利用IHB法研究輸流管道系統(tǒng)參數(shù)共振響應(yīng)，獲得精確結(jié)果。此處用該方法對多元L－P法結(jié)果進(jìn)行驗證。限于篇幅，略去求解過程。

4 算例分析

結(jié)合數(shù)值算例分析輸流管道含內(nèi)共振組合共振響應(yīng)。算例參數(shù)［14－15］為Mr＝0．447，T＝15，u＝4．75，γ＝5 000。求得管道系統(tǒng)固有頻率及它系統(tǒng)參數(shù)值為

由于ω20≈3ω10，系統(tǒng)可能會發(fā)生內(nèi)共振。由式（28）、（29）可得P，Q均為a10，f1，f2的復(fù)雜函數(shù)。由判別式（30）看出?acr或a10＜－acr時有Δ＞0，式（27）僅1個實根，即對給定的a10值，僅1個外激勵頻率Ω值與其對應(yīng)；當(dāng)時有Δ＜0，式（27）有3個不等實根，對給定的a10值，有3個外激勵頻率Ω值與其對應(yīng)，系統(tǒng)會發(fā)生強(qiáng)烈內(nèi)共振。acr為f1，f2的復(fù)雜函數(shù)，f1，f2越大acr值越小。取f1＝f2＝0．4進(jìn)行分析，求得f10＝－3．967 39×10－4，f20＝9．356 74×10－5。而據(jù)p1，p2，f10，f20值，對式（23）、（24）數(shù)值求解比較式（16）中各諧波項幅值知，系統(tǒng)響應(yīng)主要發(fā)生于q1的cosτ1項及q2的sinτ2項，故只分析a10及p2a20兩項幅值。

f1＝f2＝0．4時a10及p2a20隨激勵頻率Ω變化見圖2。由圖2看出，響應(yīng)曲線中出現(xiàn)環(huán)狀，即在一定幅值范圍內(nèi)1個振幅值對應(yīng)3個頻率值，說明系統(tǒng)已發(fā)生內(nèi)共振；各模態(tài)響應(yīng)圖上、下兩半平面基本對稱，表明均已發(fā)生內(nèi)共振；a10，a20均有兩個解，每個解橫跨上、下兩半平面；振動響應(yīng)只發(fā)生在Ω＞2ω10即外激勵頻率大于2倍的第一階固有頻率區(qū)域；外激勵幅值較小時在Ω≈2ω10（Ω＞2ω10）區(qū)域附近a10及a20呈明顯內(nèi)部共振特征，兩模態(tài)相互激勵，大小相互轉(zhuǎn)換。說明出現(xiàn)內(nèi)共振過程中兩模態(tài)間已進(jìn)行能量交換。圖2中利用IHB法求解響應(yīng)結(jié)果驗證多元L－P法結(jié)果并對比知，二者響應(yīng)振幅較小時較吻合，但隨振幅增大二者差別逐漸增大。原因為多元L－P法為攝動法，僅適合求解弱非線性振動問題，而IHB法為半解析半數(shù)值方法，對強(qiáng)、弱非線性振動均適用，故二者在小振幅時能吻合，振幅增大時多元L－P法結(jié)果誤差亦增大。

圖3、圖4為改變激勵幅值后兩模態(tài)響應(yīng)曲線。由圖3看出，f1＝f2＝2時a10及p2a20曲線中環(huán)狀部分（發(fā)生內(nèi)共振區(qū)域）明顯減小；由圖4看出，f1＝f2＝4時兩模態(tài)內(nèi)共振已完全消失。說明隨激勵幅值增大acr值減小，系統(tǒng)內(nèi)共振發(fā)生趨勢降低。激勵幅值大于某個臨界值后內(nèi)共振不再發(fā)生。響應(yīng)曲線形態(tài)也會隨激勵幅值增大發(fā)生變化，均體現(xiàn)伴隨內(nèi)共振的組合共振響應(yīng)特點。

圖2 f1＝f2＝0．4時幅－頻響應(yīng)曲線Fig．2 Frequency-amplitude curves for f1＝f2＝0．4

圖3 f1＝f2＝2時幅－頻響應(yīng)曲線Fig．3 Frequency-amplitude curves for f1＝f2＝2

圖4 f1＝f2＝4時幅－頻響應(yīng)曲線Fig．4 Frequency-amp litude curves for f1＝f2＝4

由式（23）、（24）知，外激勵幅值f1，f2不都為零，則響應(yīng)振幅a10不為零，a20只在某一激勵頻率Ω＝Ω1處可能為零。Ω1可由式（23）、（24）令a20＝0求得，亦為f1，f2的函數(shù)。第1模態(tài)a10的頻率為Ω／2，其振動為次諧波振動，存在于整個頻率域；第2模態(tài)a20的頻率為3Ω／2，其振動為超次諧波振動，存在于除Ω1的頻率域；在Ω＝Ω1處是否存在尚需額外條件判斷；第3模態(tài)f10的頻率為Ω，其振動為恒定振幅的基諧波振動，存在于整個頻率域。系統(tǒng)實際響應(yīng)為3種振動的疊加。

5 結(jié) 論

本文用多元L－P法并結(jié)合數(shù)值算例對外激勵作用下兩端固定輸流管道伴隨內(nèi)共振組合共振響應(yīng)進(jìn)行研究、分析，結(jié)論如下：

（1）對外部周期激勵的兩端固定輸流管道，當(dāng)系統(tǒng)第二階固有頻率約為第一階的3倍且激勵頻率接近前兩階固有頻率中間值時，系統(tǒng)會發(fā)生含強(qiáng)烈內(nèi)部共振的組合共振。此時前兩階模態(tài)相互激勵，振幅相互交換，致系統(tǒng)非線性動力學(xué)行為更豐富、復(fù)雜。

（2）外激勵幅值為影響內(nèi)共振的關(guān)鍵因素。隨激勵幅值增大內(nèi)共振發(fā)生趨勢降低，響應(yīng)形態(tài)亦發(fā)生變化。激勵幅值超過某臨界值后內(nèi)共振完全消失。管道系統(tǒng)受迫振動響應(yīng)通常由基諧波、次諧波及超次諧波振動疊加產(chǎn)生。

（3）利用增量諧波平衡法對多元L－P法結(jié)果驗證發(fā)現(xiàn)，二者結(jié)果在小響應(yīng)振幅時吻合較好，表明多元L－P法在分析多模態(tài)動力學(xué)系統(tǒng)（弱）非線性振動問題的有效性。

［1］Ni Q，Zhang Z L，Wang L．Application of the differential transformationmethod to vibration analysis of pipes conveying fluid［J］．Applied Mathematics and Computation，2011，217（16）：7028－7038．

［2］Wang L，Dai H L，Qian Q．Dynamics of simply supported fluid-conveying pipes with geometric imperfections［J］．Journal of Fluids and Structures，2012，29：97－106．

［3］梁峰，楊曉東，聞邦椿．脈動流激勵下輸流管道的參數(shù)共振IHB方法研究［J］．振動與沖擊，2008，27（9）：44－46．

LIANG Feng，YANG Xiao-dong，WEN Bang-chun．Analysis of parametric resonances of pipes conveying pulsating fluid by IHB method［J］．Journal of Vibration and Shock，2008，27（9）：44－46．

［4］Wang L H，Zhao Y Y．Nonlinear interactions and chaotic dynamics of suspended cables with three-to-one internal resonances［J］．International Journal of Solids and Structures，2006，43（25／26）：7800－7819．

［5］Dwivedy SK，Kar R C．Simultaneous combination and 1：3：5 internal resonances in a parametrically excited beam-mass system［J］．International Journal of Nonlinear Mechanics，2003，38（4）：585－596．

［6］Ribeiro P，Petyt M．Non-linear free vibration of isotropic plates with internal resonance［J］．International Journal of Nonlinear Mechanics，2000，35（2）：263－278．

［7］Thomas O，TouzéC，Chaigne A．Non-linear vibrations offree-edge thin spherical shells：modal interaction rules and 1：1：2 internal resonance［J］．International Journal of Solids and Structures，2005，42（11／12）：3339－3373．

［8］Xu J，Yang Q B．Flow-induced internal resonances and mode exchange in horizontal cantilevered pipe conveying fluid（I）［J］．Applied Mathematics and Mechanics（English Edition），2006，27（7）：935－941．

［9］Xu J，Yang Q B．Flow-induced internal resonances and mode exchange in horizontal cantilevered pipe conveying fluid（II）［J］．Applied Mathematics and Mechanics（English Edition），2006，27（7）：943－951．

［10］Panda LN，Kar R C．Nonlinear dynamics of a pipe conveying pulsating fluid with parametric and internal resonances［J］．Nonlinear Dynamics，2007，49（1／2）：9－30．

［11］Panda LN，Kar R C．Nonlinear dynamics of a pipe conveying pulsating fluid with combination，principal parametric and internal resonances［J］．Journal of Sound and Vibration，2008，309（3／5）：375－406．

［12］Dai H L，Wang L，Qian Q，et al．Vortex-induced vibrations of pipes conveying fluid in the subcritical and supercritical regimes［J］．Journal of Fluids and Structures，2013，39：322－334．

［13］Dai H L，Wang L．Vortex-induced vibration of pipes conveying fluid using the method of multiple scales［J］．Theoretical and App lied Mechanics Letters，2012，2（2）：022006．

［14］Namachchivaya N S，Tien W M．Non-linear dynamics of supported pipe conveying pulsating fluid．2．combination resonance［J］．International Journal of Nonlinear Mechanics，1989，24（3）：197－208．

［15］Jin JD，Song Z Y．Parametric resonances of supported pipes conveying pulsating fluid［J］．Journal of Fluids and Structures，2005，20（6）：763－783．

［16］Lau SL，Cheung Y K，Chen SH．An alternative perturbation procedure of multiple scale for nonlinear dynamics systems［J］．ASME Journal of Applied Mechanics，1989，56（3）：667－675．

［17］Chen SH，Cheung Y K，Lau SL．On the internal resonance of multi-degree-of-freedom systems with cubic nonlinearity［J］．Journal of Sound and Vibration，1989，128（1）：13－24．

［18］Yang X L，Xu W，Sun Z K，et al．Responses of strongly non-linear oscillator parametrically excited by random narrowband noise［J］．Applied Mathematics and Computation，2005，171（2）：885－899．

［19］Pu?enjak R R．Extended Lindstedt Poincar method for nonstationary resonances of dynamical systems with cubic nonlinearities［J］．Journal of Sound and Vibration，2008，314（1／2）：194－216．

［20］Pa?doussis MP，Issid N T．Dynamic stability of pipes conveying fluid［J］．Journal of Sound and Vibration，1974，33（3）：267－294．

［21］Pa?doussis MP，Issid N T．Experiments on parametric resonance of pipes containing pulsatile flow［J］．ASME Journal of Applied Mechanics，1976，43（2）：198－202．

［22］Lau S L，Cheung Y K．Amplitude incremental variational principle for non-linear vibration of elastic systems［J］．ASME Journal of Applied Mechanics，1981，48（4）：959－964．

Nonlinear vibration of a fluid-conveying pipe under external excitation accom panied with internal resonance

LIANG Feng，BAO Ri-dong，JIN Ying，SU Yong
（School of Energy and Power Engineering，Shenyang University of Chemical Technology，Shenyang 110142，China）

An external periodic load was considered to act on a fluid-conveying pipe clamped at both ends，and the nonlinear forced vibration of such a system was explored by the multidimensional Lindstedt-Poincaré（MDLP）method．According to the analysis，when the second natural frequency of the system is nearly thrice the firstone，and the excitation frequency is nearly at the middle of first two natural frequencies，accompanied internal resonance may occur to form a combination resonance．The characteristics of the response were discussed，and the motions of first two modes were investigated in detail．The influence of excitation amplitude on the internal resonance was analyzed．Some numerical examples reveal the rich and complex dynamic behaviors caused by internal resonance and show that the occurrence tendency of internal resonance will die down and the responsemodeswill vary with the increase of excitation amplitude．The convenience and efficiency of the MDLPmethod in predicting nonlinear dynamics are demonstrated by the results of the study．

fluid-conveying pipe；nonlinear vibration；internal resonance；external periodic excitation；multidimensional Lindstedt-Poincarémethod

O326

：A

10．13465／j．cnki．jvs．2014．22．026

國家自然科學(xué)基金（51275315）；遼寧省教育廳科研項目（L2013160）

2013－08－08 修改稿收到日期：2013－12－04

梁峰男，博士，講師，1979年8月生