莘智，侯瑾蓉

（內(nèi)蒙古師范大學(xué)數(shù)學(xué)科學(xué)學(xué)院，內(nèi)蒙古 呼和浩特 010022）

1830年，泊松在研究單擺的振動(dòng)時(shí)，提出將非線性系統(tǒng)的解按小參數(shù)ε的冪次展開的近似計(jì)算方法稱為攝動(dòng)法或小參數(shù)法［1］。1892年，林滋泰德為了消除在天文學(xué)中利用正規(guī)攝動(dòng)法求近似解時(shí)出現(xiàn)的久期項(xiàng)問題，最先提出對正規(guī)攝動(dòng)法進(jìn)行改進(jìn)。后來，龐加萊對改進(jìn)后的方法的合理性進(jìn)行了證明，故而稱為林滋泰德-龐加萊法，它的出現(xiàn)很好地解決了系統(tǒng)共振，即久期項(xiàng)問題。林滋泰德-龐加萊法認(rèn)為解和系統(tǒng)的激勵(lì)頻率都是ε的未知函數(shù)，都要展成ε的冪級數(shù)形式，通過確定級數(shù)形式中的待定系數(shù)，從而求得解析解和頻率［2］。本文主要通過林滋泰德-龐加萊法對達(dá)芬系統(tǒng)的幾種振動(dòng)求解，得到了精確度較高的近似解析解。

1 達(dá)芬系統(tǒng)的自由振動(dòng)

達(dá)芬系統(tǒng)的自由振動(dòng)動(dòng)力學(xué)方程為［1］

規(guī)定初始條件為［2］

此時(shí)將該動(dòng)力學(xué)方程的解展開為冪級數(shù)的形式

同時(shí)，將自由振動(dòng)的頻率ω也展成ε的冪級數(shù)形式

將式（4）兩邊平方，得

將式（3）和式（5）代入式（1），引入新的自變量ψ=ωt，將原來的微分改定義為對ψ的微分，轉(zhuǎn)化為

令等式兩邊ε相同次冪項(xiàng)的系數(shù)相等，可得方程

規(guī)定各方程的初始條件為

由零次近似方程（7）和初始條件可以解出

將式（12）的解代入式（8）的右邊，可得

為消除方程中出現(xiàn)的久期項(xiàng)，需要令方程右邊的cosψ項(xiàng)的系數(shù)等于零，于是，可以推導(dǎo)出

這時(shí)，在滿足此條件的基礎(chǔ)上解出x1的值，設(shè)

在方程（13）滿足初始條件（11）的情況下，可得方程的解為

將式（12）和式（15）代入方程（9）中，整理得

為消除方程中出現(xiàn)的久期項(xiàng)，需要令方程右邊cosψ項(xiàng)的系數(shù)等于零，導(dǎo)出

此時(shí)在滿足此條件的基礎(chǔ)上，設(shè)

在方程必須滿足初始條件情況下，可得方程的解為

將式（12），式（15），式（18）代入方程（10）中并加以整理得

為消除方程中出現(xiàn)的久期項(xiàng)，需要令方程右邊的cosψ項(xiàng)的系數(shù)等于零，導(dǎo)出

解x3的過程與前面解x1和x2的過程相同，重復(fù)操作，可得

由此可得方程（1）的三階精度的解

以及自由振動(dòng)頻率ω的表達(dá)式

2 接近共振的受迫振動(dòng)

討論帶微弱阻尼的達(dá)芬系統(tǒng)接近共振的受迫振動(dòng)，其動(dòng)力學(xué)方程為

設(shè)阻尼項(xiàng)，激勵(lì)頻率以及系統(tǒng)的固有頻率均與ε同數(shù)量級，令

同樣將方程（24）的解設(shè)為式（3）的形式，與式（25）一起代入方程（24）中，再令F0=ω20B1，引入新的自變量ψ=ωt，令ε的同次冪項(xiàng)的系數(shù)相等，可以導(dǎo)出以下一系列的方程。

由方程（26）可以解出

將此解代入方程（27）中并整理得

為消去共振對系統(tǒng)的影響，需要消除久期項(xiàng)，得到

消去久期項(xiàng)后，可以設(shè)方程（30）的解為

把所設(shè)方程的解代入方程（30）中，可以解得

將式（32）與式（29）代入式（28）中得

整理得

同樣地，消去久期項(xiàng)得

設(shè)方程（33）的解為

將式（34）代入方程（33），整理后得

由此得

故而方程（24）的二階近似解便可給出

3 亞諧波共振

討論亞諧波共振的動(dòng)力學(xué)方程

同樣設(shè)該方程的解為式（3）的形式，并令

將式（3）與式（39）同時(shí)代入方程（38）中，得到

令上式中兩邊ε的同次冪項(xiàng)的系數(shù)分別相等，由此可以推導(dǎo)出一系列的方程。

設(shè)方程（41）解的形式為

其中，A13由初始條件決定，將式（44）代入方程（41）中得

將式（44）代入方程（42）中并整理得

為了消除方程中的久期項(xiàng)，得到

或

解出A13的不等于零的解

消去久期項(xiàng)后，設(shè)

將式（50）代入式（46）中，可以得到

將式（50）和（44）代入方程（43）中，得到

整理后得到

為了消除久期項(xiàng)，得

綜上，方程（38）的解可以寫成

4 結(jié)論

林滋泰德為了消除天文學(xué)中的久期項(xiàng)改進(jìn)正規(guī)攝動(dòng)法，龐加萊為改進(jìn)的攝動(dòng)法的合理性進(jìn)行了數(shù)學(xué)證明，林滋泰德-龐加萊法和正規(guī)攝動(dòng)法同稱為小參數(shù)法，也叫PL攝動(dòng)法。該方法的基本思想是認(rèn)為非線性系統(tǒng)的固有頻率ω并不等于派生系統(tǒng)的固有頻率ω0，而應(yīng)該是小參數(shù)ε的未知函數(shù)。因此在將基本解展成ε的冪級數(shù)的同時(shí)，因?qū)㈩l率也寫成ε的冪級數(shù)，冪級數(shù)的待定系數(shù)根據(jù)周期運(yùn)動(dòng)的要求依次確定。本文以達(dá)芬系統(tǒng)的自由振動(dòng)為例，利用林滋泰德-龐加萊法給出了三階精度的近似解析解，相比較文獻(xiàn)［1］從二階提高到了三階，進(jìn)而提高了解的精確程度；討論了帶微弱阻尼的達(dá)芬系統(tǒng)接近共振的受迫振動(dòng)，利用該方法給出了二階精度的近似解，比文獻(xiàn)［1］中的解提高了一階精度；針對亞諧波共振的情形，將文獻(xiàn)［1］中的解的精度由一階提高到了二階。利用林滋泰德-龐加萊法對達(dá)芬系統(tǒng)的三種不同振動(dòng)的解求解，實(shí)現(xiàn)了非線性振動(dòng)系統(tǒng)解的精度的提高，進(jìn)而能夠更好地分析非線性系統(tǒng)的運(yùn)動(dòng)規(guī)律，以及對系統(tǒng)參數(shù)和初始條件的依賴關(guān)系。但是，得到非線性振動(dòng)高階近似解，會(huì)增加很大的工作量，因此往往要截?cái)喔唠A項(xiàng)只保留級數(shù)的前面有限項(xiàng)，同時(shí)解的最后結(jié)果還需要試驗(yàn)來驗(yàn)證。

小參數(shù)法是從事理論研究的重要數(shù)學(xué)工具之一，對于弱非線性問題尤其有作用。理論的研究從實(shí)際問題中來，并最終應(yīng)用到實(shí)際問題中，小參數(shù)法在基礎(chǔ)和應(yīng)用研究中已被廣泛應(yīng)用于微分方程、軌道力學(xué)、非線性振動(dòng)、固體力學(xué)、流體力學(xué)等領(lǐng)域［3-6］，并且隨著科學(xué)的不斷發(fā)展，一定會(huì)研究出更多有關(guān)小參數(shù)的理論，對它的應(yīng)用也會(huì)更加廣泛。