黃建亮 肖龍江

摘要：研究兩端固定屈曲梁這種同時含有2，3次非線性項的系統(tǒng)受基礎簡諧激勵作用下的非線性振動響應及分岔演化過程。對屈曲梁的運動微分方程，利用Galerkin方法分離時間和空間變量，應用增量諧波平衡（IHB）法自動追蹤屈曲梁的非線性振動響應的周期解和倍周期解，并用Floquet理論分析其解的穩(wěn)定性。研究表明屈曲梁對稱模態(tài)的固有頻率隨屈曲程度而變化，反對稱模態(tài)的固有頻率保持不變。研究發(fā)現(xiàn)基礎簡諧激勵作用下，不同屈曲程度時存在兩種截然不同的非線性現(xiàn)象：1）在非共振時，反對稱模態(tài)未能被激發(fā)，系統(tǒng)經(jīng)過發(fā)生倍周期分岔通向混沌運動;2）在1∶1內(nèi)共振條件下，反對稱模態(tài)在一定的頻率區(qū)間里會被激發(fā)，系統(tǒng)發(fā)生Hopf分岔導致具有等間距邊頻帶的準周期運動，最后至混沌運動。利用IHB法得到結果與用Runge-Kutta法得到的數(shù)值結果一致。

關鍵詞：非線性振動;屈曲梁;分岔;增量諧波平衡法;準周期運動

中圖分類號：O322;TU375.1文獻標志碼：A文章編號：1004-4523（2020）04-0698-11

DOI：10.16385/j.cnki.issn.1004-4523.2020.04.007

引言

屈曲梁/板/殼等結構的振動十分普遍，廣泛存在于航空航天、土木、機械、電子等工程中[1-5]。飛機、火箭、導彈、加速推進器、航天結構里的柔性纖維、潛水器等典型例子，其基礎的激勵作用都與結構的橫向振動有關。

由于初始彎曲和中面伸長，屈曲梁的振動方程中同時包含2次和3次非線性項，表現(xiàn)出較為復雜的非線性特性[6-7]。兩端固定屈曲梁，由于邊界條件的特殊性，其振動響應更為復雜[8-9]。Afaneh和Ibrahim利用多尺度法分析了1∶1內(nèi)共振條件下兩端固定屈曲梁的前兩階的振動響應[8]，發(fā)現(xiàn)在此條件下，梁的振動響應會出現(xiàn)周期運動、準周期運動和混沌運動三種形態(tài)，并通過實驗進行驗證。Kreider和Nayfeh通過實驗測定[9]，當外激勵的頻率接近于梁的第1階固有頻率時，隨著外激勵的振幅變化會發(fā)生倍周期分岔，并且在倍周期分岔的基礎上會發(fā)生“不可解釋”的準周期運動，然而沒有給出理論上的分析，只是說明用1階模態(tài)計算得到的結果誤差較大。Emam和Nayfeh通過數(shù)值計算的方法驗證了上述實驗結果的可靠性[10]，同時他們還用多時間尺度法和打靶法研究了屈曲梁的次諧波響應[11]。在中國，也有相關學者在這領域作了大量而出色的研究工作。曾志剛和葉敏采用實驗建模和多尺度法研究了黏彈性屈曲梁的參數(shù)激振[12]，并用數(shù)值模擬的方法研究由倍周期分岔通往混沌的道路。陳得良等用實驗研究了脫層屈曲梁在軸向周期激勵作用下的倍周期分岔與混沌運動[13]。

增量諧波平衡法（IHB法）適用于解決強非線性振動問題[14-17]。運用增量諧波平衡法可追蹤振動響應隨著系統(tǒng)參數(shù)的變化，有利于研究系統(tǒng)的總體特性。Huang等[18-19]用增量諧波平衡法研究了兩端簡支彎曲梁的振動。Lau等[20-21]把增量諧波平衡法和多時間尺度法結合起來，用于計算系統(tǒng)的準周期運動。Huang和Zhu[22-23]改進了文獻[20-21]的方法，用多時間尺度的增量諧波平衡法研究了一端固定一端簡支的彎曲梁的準周期運動。

本文采用增量諧波平衡法分析兩端固定屈曲梁的非線性振動響應和分岔過程。先用線性化的方法計算梁的固有頻率和模態(tài)函數(shù)，再通過Galerkin方法對梁的振動偏微分方程進行離散，使其轉化為以時間為變量的常微分方程，然后運用增量諧波平衡法把此常微分方程再次轉化為代數(shù)方程求解，得到系統(tǒng)振動的周期解，并通過Floquet理論[24-25]對解進行穩(wěn)定性分析。研究了彎曲程度不同時，遠離內(nèi)共振條件下和內(nèi)共振條件下系統(tǒng)的非線性振動響應和分岔演化過程。

周期解發(fā)生倍周期分岔的原因和工況1）不同，在工況1）中是因為對稱模態(tài)的響應發(fā)生倍周期分岔，而在工況2）中則是因為被激發(fā)的反對稱模態(tài)與對稱模態(tài)中包含不同的頻率成分。

追蹤period-2周期解可發(fā)現(xiàn)，在圖6中r點外激勵頻率為45.68處，period-2周期解不穩(wěn)定。在這一點，F(xiàn)loquet乘子中有一對復共軛特征值與單位圓相交，因此，在這一點period-2周期解發(fā)生Hopf分岔導致準周期運動。經(jīng)過一段穩(wěn)定的準周期運動，在外激勵頻率為45.6處，準周期運動失穩(wěn)，演變?yōu)榛煦邕\動。與工況1）不同的是，此時產(chǎn)生的混沌運動并不會隨著時間逐漸衰減為周期運動，而是會持續(xù)保持下去。

圖7-9反映了在1∶1內(nèi)共振條件下，圖6中q點附近外激勵頻率接近梁的第1階固有頻率時，梁的振動響應的分岔演化過程。其中圖7（a）為外激勵頻率為45.9處（圖6中q點右側）穩(wěn)定的period-1周期解的頻譜圖。圖7（b）為外激勵頻率為45.7處（圖6中r點右側）穩(wěn)定的period-2周期解的頻譜圖，可以看出，與圖7（a）中period-1周期解相比，q1的運動幾乎沒有發(fā)生變化，而q2的運動卻從無到有被激發(fā)出來，圖8（a）與（d）分別為此時q1與q2的相圖。圖7（c）為外激勵頻率為45.64處準周期運動的頻譜圖，在不穩(wěn)定的period-2周期解各個諧波項的附近出現(xiàn)了間隔較小的等間距的邊頻帶，邊頻帶中相鄰諧波項的間隔約為4.11，存在等間距邊頻帶為準周期運動頻譜圖所特有的現(xiàn)象[22-23]，圖8（b）和（e）分別為此時q1與q2準周期運動的相圖，圖9（a）和（c）為其龐加萊截面圖。圖8（c）與圖8（f）為外激勵頻率為45.6時，q1與q2混沌運動的相圖，圖9（b）和（d）為其龐加萊截面圖。

當外激勵頻率接近梁的第2階固有頻率時，梁的振動響應的變化和外激勵頻率接近梁的第1階固有頻率。當外激勵頻率從44.68減小至43.68時（圖6中s與t點之間），由period-1周期解發(fā)生倍周期分岔產(chǎn)生的period-2周期解穩(wěn)定。此后隨著外激勵頻率的減小，period-2周期解發(fā)生Hopf分岔，產(chǎn)生準周期運動。繼續(xù)減小外激勵頻率，準周期運動將失穩(wěn)變?yōu)榛煦邕\動。

5結論

本文用增量諧波平衡法追蹤了在基礎諧波激勵作用下兩端固定屈曲梁的振動響應隨外激勵頻率的變化，并分析了振動響應的穩(wěn)定性和分岔演化過程。研究發(fā)現(xiàn)，在僅有梁的彎曲程度不同的兩種情況中，運動失穩(wěn)的方式完全不同。

在屈曲程度較小，非共振的條件下，屈曲梁的反對稱模態(tài)的振動不會被激發(fā)，當外激勵頻率接近梁的第1階固有頻率或其2倍時，梁的對稱模態(tài)振動響應不斷發(fā)生倍周期分岔，從而依次得到period-2k（k=1，2…）周期解，直到某個period-2k周期解不再發(fā)生倍周期分岔而發(fā)生鞍結分岔。在鞍結分岔點處period-2k周期解失穩(wěn)而變?yōu)榛煦邕\動，此混沌運動隨著時間的變化逐漸演變?yōu)橐粋€穩(wěn)定的周期運動，這個周期運動為混沌運動中的一個吸引子。

在屈曲程度較大，存在1∶1內(nèi)共振的條件下，隨著外激勵頻率接近梁的第1階或第2階固有頻率，梁的振動響應也會發(fā)生倍周期分岔。但此時振動響應發(fā)生倍周期分岔的原因是反對稱模態(tài)的振動被激發(fā)，而被激發(fā)的反對稱模態(tài)的振動響應中所含有的頻率成分均為外激勵頻率一半的奇數(shù)倍。由于發(fā)生倍周期分岔，梁的振動響應變?yōu)閜eriod-2周期解。隨著外激勵頻率的減小，period-2周期解會發(fā)生Hopf分岔，導致含有等間距邊頻帶的準周期運動。繼續(xù)減小外激勵的頻率，準周期運動也會失穩(wěn)而變?yōu)榛煦邕\動。

增量諧波平衡法可與Floquet理論相結合，用于定量追蹤非線性系統(tǒng)的分岔演化過程。它的分析結果和數(shù)值方法的結果一致。

參考文獻：

[1]田振國，梁金奎.矩形磁流變彈性薄板的動力穩(wěn)定性分析[J].振動工程學報，2018，31（1）：125-131.

TianZhenguo，LiangJinkui.Dynamicstabilityanalysisofathinrectangularmagnetorheologicalelasticplate[J].JournalofVibrationEngineering，2018，31（1）：125-131.

[2]ZineA，TounsiA，DraicheK，etal.Anovelhigher-ordersheardeformationtheoryforbendingandfreevibrationanalysisofisotropicandmultilayeredplatesandshells[J].SteelandCompositeStructures，2018，26（2）：125-137.

[3]熊柳楊，張國策，丁虎，等.黏彈性屈曲梁非線性內(nèi)共振穩(wěn)態(tài)周期響應[J].應用數(shù)學和力學，2014，35（11）：1188-1196.

XiongLiuyang，ZhangGuoce，DingHu，etal.Steady-stateperiodicresponsesofaviscoelasticbuckledbeaninnonlinearinternalresonance[J].AppliedMathematicsandMechanics，2014，35（11）：1188-1196.

[4]ZareiM，F(xiàn)aghaniG，GhalamiM，etal.Bucklingandvibrationanalysisoftaperedcircularnanoplate[J].JournalofAppliedandComputationalMechanics，2018，4（1）：40-54.

[5]張泉，周麗平，金家楣，等.高速柔性并聯(lián)平臺的動力學分析[J].振動工程學報，2015，28（1）：27-37.

ZhangQuan，ZhouLiping，JinJiamei，etal.Dynamicanalysisofahigh-speedflexibleparallelmanipulator[J].JournalofVibrationEngineering，2015，28（1）：27-37.

[6]NayfehAH，KreiderW，AndersonTJ.Investigationofnaturalfrequenciesandmodeshapesofbuckledbeams[J].AmericanInstituteofAeronauticsandAstronauticsJournal，1995，33（6）：1121-1126.

[7]RoncenT，SinouJJ，LambelinJP.Non-linearvibrationsofabeamwithnon-idealboundaryconditionsanduncertainties—Modeling，numericalsimulationsandexperiments[J].MechanicalSystemsandSignalProcessing，2018，110（15）：165-179.

[8]AfanehAA，IbrahimRA.Nonlinearresponseofaninitiallybuckledbeamwith1∶1internalresonancetosinusoidalexcitation[J].NonlinearDynamics，1993，4（6）：547-571.

[9]KreiderW，NayfehAH.Experimentalinvestigationofsingle-moderesponsesinafixed-fixedbuckledbeam[J].NonlinearDynamics，1998，15（2）：155-177.

[10]EmamSA，NayfehAH.Onthenonlineardynamicsofabuckledbeamsubjectedtoaprimary-resonanceexcitation[J]，NonlinearDynamics，2004，35（1）：1-17.

[11]EmamSA，NayfehAH.Nonlinearresponsesofbuckledbeamstosubharmonic-resonanceexcitations[J]，NonlinearDynamics，2004，35（2）：105-122.

[12]曾志剛，葉敏.黏彈性復合材料屈曲梁非線性參數(shù)振動的穩(wěn)定性和分岔分析[J].振動與沖擊，2012，31（2）：129-135.

ZengZhigang，YeMin.Stabilityandbifurcationofaparametricallyexcitedviscoelasticcompositebuckledbeam[J].JournalofVibrationandShock，2012，31（2）：129-135.

[13]陳得良，付欽，陳昌萍.不同脫層位置下脫層屈曲梁振動實驗研究[J].動力學與控制學報，2015，13（6）：468-474.

ChenDeling，F(xiàn)uQin，ChenChangping.Experimentalstudyonthevibrationofdelaminatedbuckledbeamswithdifferentdelaminationlocations[J].JournalofDynamicsandControl，2015，13（6）：468-474.

[14]姚紅良，王重陽，王帆，等.多頻激勵局部非線性系統(tǒng)響應求解的降維增量諧波平衡法[J].振動工程學報，2015，28（5）：741-747.

YaoHongliang，WangChongyang，WangFan，etal.Thedemension-reductiveincrementalharmonicbalancemethodforsolvingtheresponseoflocalnonlineardynamicsystemwithmulti-frequencyexcitation[J].JournalofVibrationEngineering，2015，28（5）：741-747.

[15]晏致濤，張海峰，李正良.基于增量諧波平衡法的覆冰輸電線舞動分析[J].振動工程學報，2012，25（2）：161-166.

YanZhitao，ZhangHaifeng，LiZhengliang.Gallopinganalysisoficedtransmissionlinesbasedonincrementalharmonicbalancemethod[J].JournalofVibrationEngineering，2012，25（2）：161-166.

[16]王威，宋玉玲，李瑰賢.基于增量諧波平衡法的汽車轉向系非線性動力學特性研究[J].振動工程學報，2010，23（3）：355-360.

WangWei，SongYuling，LiGuixian.Nonlineardynamicsofautomobilesteeringusingincrementalharmonicbalancemethod[J].JournalofVibrationEngineering，2010，23（3）：355-360.

[17]DuYZ，WangWH，WangLL，etal.Nonlineardynamicsofheavemotionofthesandglass-typefloatingbodywithpiecewise-nonlinear，time-varyingstiffness[J].MarineStructures，2018，60：136-150.

[18]HuangJL，SuRKL，LeeYY，etal.Nonlinearvibrationofacurvedbeamunderuniformbaseharmonicexcitationwithquadraticandcubicnonlinearities[J].JournalofSoundandVibration.2011，330（21）：5151-5164.

[19]HuangJL，SuKLR，LeeYYR，etal.Variousbifurcationphenomenainanonlinearcurvedbeamsubjectedtobaseharmonicexcitation[J].InternationalJournalofBifurcationandChaos，2018，28（7）：1830023.

[20]LauSL，CheungYK，WuSY.Incrementalharmonicbalancemethodwithmultipletimescalesforaperiodicvibrationofnonlinearsystems[J].JournalofAppliedMechanics，1983，50（4A）：871-876.

[21]LauSL，CheungYK.IncrementalHamilton′sprinciplewithmultipletimescalesfornonlinearaperiodicvibrationsofshells[J].JournalofAppliedMechanics，1986，53（2）：465-466.

[22]HuangJL，ZhuWD.Anincrementalharmonicbalancemethodwithtwotimescalesforquasiperiodicmotionofnonlinearsystemswhosespectrumcontainsuniformlyspacedsidebandfrequencies[J].NonlinearDynamics，2017，90（2）：1015-1033.

[23]HuangJL，ZhuWD.Anewincrementalharmonicbalancemethodwithtwotimescalesforquasi-periodicmotionsofanaxiallymovingbeamwithinternalresonanceundersingle-toneexternalexcitation[J].JournalofVibrationandAcoustics，2017，139（2）：021010.

[24]Karlicˇic'D，Cajic'M，AdhikariS.Dynamicstabilityofanonlinearmultiple-nanobeamsystem[J].NonlinearDynamics，2018，93（3）：1495-1517.

[25]廖海濤，高歌.求解非線性系統(tǒng)共振峰值的限制優(yōu)化打靶法[J].振動工程學報，2014，27（2）：166-171.

LiaoHaitao，GaoGe.Constrainedoptimizationshootingmethodforfindingtheresonantpeakofnonlinearsystems[J].JournalofVibrationEngineering，2014，27（2）：166-171.

Abstract：Thisstudyinvestigatesnonlineardynamicsandbifurcationofafixed-fixedbuckledbeamsubjectedtoauniformbaseharmonicexcitation，whichisgovernedbyacouplednonlinearequationwithbothquadraticandcubicnonlinearities.TheGalerkinmethodisemployedtoseparatetimeandspacevariablesforadifferentialequationofthebuckledbeam，andtheincrementalharmonicbalance（IHB）methodisusedtotrackperiodicandperiodic-doublingsolutionsofnonlinearresponsesofthebuckledbeamautomatically.Inaddition，theFloquettheoryisusedtoanalyzestabilityofthesolutions.Theinvestigationsuggeststhatthenaturalfrequenciesofthesymmetricalmodeofthebuckledbeamvarywiththebucklinglevel，whilethenaturalfrequenciesoftheanti-symmetricmodeareinvariant.Itisfoundthattherearetwodifferentinterestingnonlinearphenomenawithbucklinglevelunderthebaseharmonicexcitation，1）innonresonancecase，theanti-symmetricmodescannotbeexcitedandthesystembifurcatesintoperiod-doublingsolutions;2）in1：1internalresonancecase，theanti-symmetricmodeswillbeexcitedinacertainfrequencyrange，andHopfbifurcationofthesystemwillleadtoquasiperiodicmotionwithuniformsidebands，untilchaosoccurs.Furthermore，theresultsobtainedbytheIHBmethodagreeverywellwiththoseobtainedbythenumericalintegrationusingtheRunge-Kuttamethod.

Keywords：nonlinearvibration;buckledbeam;bifurcation;incrementalharmonicbalancemethod;quasiperiodicmotion