• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      非線性振動系統(tǒng)的多項式向量方法1)

      2023-11-16 06:42:40金棟平
      力學(xué)學(xué)報 2023年10期
      關(guān)鍵詞:高維算例乘法

      金棟平

      (南京航空航天大學(xué)航空航天結(jié)構(gòu)力學(xué)及控制全國重點(diǎn)實驗室,南京 210016)

      引言

      非線性振動系統(tǒng)近似解法主要有基于派生系統(tǒng)解特征而衍生出的方法,如平均法等,以及基于攝動展開的方法,如Lindstedt-Poincaré方法、KBM 法、多尺度法和源于濾波的方法,如諧波平衡法等[1-5].學(xué)者們采用上述方法成功揭示了工程非線性振動中的復(fù)雜動力學(xué)行為,如滯后非線性中的分岔[6]、相軌線跳躍現(xiàn)象[7],并用于非線性隔振器設(shè)計[8]等.現(xiàn)代工程結(jié)構(gòu)和系統(tǒng)大多帶有柔性附件、振動控制器[9],使得受控動力學(xué)模型表現(xiàn)為很高維數(shù)之特征.借助動力學(xué)等效等模型進(jìn)行降階,可以大大降低動力學(xué)模型的維數(shù)[10-11],使得由周期單元組成的復(fù)雜結(jié)構(gòu)非線性振動能夠進(jìn)行半解析和半數(shù)值分析[12-13].然而,當(dāng)非線性振動系統(tǒng)的維數(shù)較高時,上述攝動方法、多尺度方法等需逐步、依次求解微分方程組,加上計算機(jī)符號軟件可以不支持海量指數(shù)的提取[14-15]等運(yùn)算,需要人工參與,以致工作量巨大、計算效率低.

      本文針對高維非線性振動系統(tǒng)的近似求解問題,通過常數(shù)矩陣和多項式向量函數(shù)的乘法,將非線性部分表示成為矩陣和向量積形式,以致冪次近似解均滿足一組非齊次狀態(tài)方程,進(jìn)而根據(jù)一階非齊次方程組的求解方法,一次性地獲得全部的解析解.由于帶有小參數(shù)的多項式向量函數(shù)之間的乘法運(yùn)算可以借助Toeplitz 矩陣表示成矩陣之間的乘法運(yùn)算,借助計算機(jī)符號軟件等工具可以實現(xiàn)計算過程的程序化.最后,通過一個四維彈簧擺振動的算例,展示了多項式向量方法在解決高維非線性振動時的可行性.

      1 多項式向量

      研究n維非線性振動系統(tǒng),狀態(tài)方程為

      一般地,可以將方程(1)寫成多項式向量形式

      其中,k次齊次多項式向量函數(shù)為

      式中,ei=[0 0···1···0 0]T表示第i列元素為1 的單位向量,mi為非負(fù)整數(shù).

      若采用攝動法,將

      代入方程(4)并比較方程兩邊 ε同次冪系數(shù),得

      2 多項式乘法

      定義T階滯后算子

      該算子由單位矩陣IT=[e1e2···eT]通過消去第一列并在最后一列添加零列向量形成,在次對角線上的元素為1、其余均為0.此外,左下角元素為1;對r≥0,.例如

      利用上述滯后算子,可以將多項式相乘轉(zhuǎn)換成Toeplitz 型下三角矩陣的乘法形式.將小參數(shù)視為變量,令

      因而

      即x1(L4)x2(L4)第1 列元素χ(q)正是兩個多項式相乘后的 εq次冪的系數(shù).根據(jù)方程式(5),第k次齊次多項式函數(shù)的 εq次冪的系數(shù)為

      式中,T=q+1.依據(jù)式(16),即可確定方程式(7)~式(10)中任意次冪系數(shù)的多項式函數(shù),繼而依次求出各階近似解.

      類似地,有

      可見,非線性振動系統(tǒng)中齊次多項式乘法運(yùn)算可以轉(zhuǎn)換成一系列Toeplitz 矩陣的乘法運(yùn)算,以致原本繁瑣的逐個攝動求解轉(zhuǎn)變?yōu)橐幌盗芯仃囆问降姆驱R次方程求解問題.這樣做的好處就是,對于高維非線性振動系統(tǒng),借助計算機(jī)符號軟件對于多項式向量函數(shù)和矩陣運(yùn)算的優(yōu)勢,通過多項式矩陣乘法的運(yùn)算,無需對二階微分方程逐個求解,從而一次性地獲得形式上結(jié)構(gòu)簡潔的全部近似解.值得注意的是,直接攝動法得到的解在t∈[0,1/ε]內(nèi)有很好的近似.

      3 算例

      (1)算例1

      首先,考慮二維平方非線性振動系統(tǒng)

      根據(jù)式(16),有

      因而,一次近似解為x=x(0)+εx(1),該結(jié)果與文獻(xiàn)[4]直接攝動法給出的結(jié)果相同.可以繼續(xù)求滿足初始條件的二次冪方程的解,即

      可見,借助Maple 等計算機(jī)符號軟件工具,上述建立在多項式向量求解方法上的矩陣運(yùn)算程序很容易編寫和執(zhí)行.注意到,式(28)中出現(xiàn)了九期項tsint.下面研究更高維數(shù)的非線性振動問題.

      (2)算例2

      考慮靜平衡位置附近彈簧擺的面內(nèi)微振動,此時近似到二階的非線性振動方程為

      寫成狀態(tài)方程形式,令x=x1,=x2,y=x3,=x4,有

      根據(jù)式(16),對于二次近似解,有

      為簡潔起見,取 ω=1,有

      將式(30)和式(32)代入式(36),得到二次冪方程的解析表達(dá)式.為簡潔起見,取 α=1,β=2,則

      因而,二次近似解為x=x(0)+εx(1)+ε2x(2).注意到,式(37)中出現(xiàn)了九期項tsint和tcost.

      圖1 近似解析解與數(shù)值結(jié)果對比(=[0.01 0 0.01 0]T)Fig.1 Comparison of the proposed method and simulation results for the differential small parameters (=[0.01 0 0.01 0]T)

      圖1 近似解析解與數(shù)值結(jié)果對比(=[0.01 0 0.01 0]T)(續(xù))Fig.1 Comparison of the proposed method and simulation results for the differential small parameters (=[0.01 0 0.01 0]T)(continued)

      圖2 近似解與數(shù)值結(jié)果對比(=[0.1 0 0.1 0]T)Fig.2 Comparison of the proposed method and simulation results for the differential small parameters (=[0.1 0 0.1 0]T)

      4 結(jié)論

      對于齊次多項式作為非線性的高維振動系統(tǒng),可以通過多項式向量將非線性部分表示成常數(shù)矩陣和多項式向量之積,進(jìn)而將非線性振動系統(tǒng)表示成一系列矩陣和向量的乘積形式.在采用基于攝動展開的方法時,小參數(shù)作為冪級數(shù)出現(xiàn),因而可以采用Toeplitz 矩陣實現(xiàn)冪次近似解作為元素的多項式之間的乘法運(yùn)算.因此,非線性振動系統(tǒng)的解完全由矩陣和多項式向量之間的乘法運(yùn)算獲得,從而有望解決高維非線性振動系統(tǒng)的求解難題.

      猜你喜歡
      高維算例乘法
      算乘法
      我們一起來學(xué)習(xí)“乘法的初步認(rèn)識”
      《整式的乘法與因式分解》鞏固練習(xí)
      把加法變成乘法
      一種改進(jìn)的GP-CLIQUE自適應(yīng)高維子空間聚類算法
      基于加權(quán)自學(xué)習(xí)散列的高維數(shù)據(jù)最近鄰查詢算法
      基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
      互補(bǔ)問題算例分析
      一般非齊次非線性擴(kuò)散方程的等價變換和高維不變子空間
      基于CYMDIST的配電網(wǎng)運(yùn)行優(yōu)化技術(shù)及算例分析
      海安县| 宁明县| 司法| 黄龙县| 荥阳市| 黔南| 霍城县| 五台县| 贵溪市| 四会市| 永顺县| 武穴市| 长沙市| 江达县| 西林县| 洞口县| 桂东县| 漯河市| 乾安县| 徐水县| 高雄县| 镇安县| 宿州市| 河东区| 昭觉县| 浙江省| 正蓝旗| 正安县| 绍兴县| 湟中县| 宜春市| 清苑县| 洪泽县| 仙游县| 固原市| 明溪县| 泾源县| 外汇| 柘荣县| 长武县| 武邑县|