亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

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

        2023-11-16 06:42:40金棟平
        力學學報 2023年10期
        關鍵詞:高維算例乘法

        金棟平

        (南京航空航天大學航空航天結構力學及控制全國重點實驗室,南京 210016)

        引言

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

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

        1 多項式向量

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

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

        其中,k次齊次多項式向量函數為

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

        若采用攝動法,將

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

        2 多項式乘法

        定義T階滯后算子

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

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

        因而

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

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

        類似地,有

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

        3 算例

        (1)算例1

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

        根據式(16),有

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

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

        (2)算例2

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

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

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

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

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

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

        圖1 近似解析解與數值結果對比(=[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 近似解析解與數值結果對比(=[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 近似解與數值結果對比(=[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 結論

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

        猜你喜歡
        高維算例乘法
        算乘法
        我們一起來學習“乘法的初步認識”
        《整式的乘法與因式分解》鞏固練習
        把加法變成乘法
        一種改進的GP-CLIQUE自適應高維子空間聚類算法
        測控技術(2018年4期)2018-11-25 09:46:48
        基于加權自學習散列的高維數據最近鄰查詢算法
        電信科學(2017年6期)2017-07-01 15:44:37
        基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
        互補問題算例分析
        一般非齊次非線性擴散方程的等價變換和高維不變子空間
        基于CYMDIST的配電網運行優(yōu)化技術及算例分析
        国产偷国产偷亚洲高清视频| 午夜国产精品久久久久| 热热久久超碰精品中文字幕| 国语淫秽一区二区三区四区| 国产欧美一区二区三区在线看| 极品少妇被黑人白浆直流| 狼人香蕉香蕉在线28 - 百度| 老熟妇乱子伦av| 亚洲成人免费无码| 中文字幕亚洲视频三区| 亚洲精品在线国产精品| 亚洲日韩av无码中文字幕美国| 精品人妻无码中文字幕在线| 一区二区黄色素人黄色| 欧美xxxxx高潮喷水麻豆| 婷婷五月综合丁香在线| 亚洲av成人一区二区三区网址| 亚洲av中文字字幕乱码软件| 国产精品久久国产精麻豆99网站 | av免费观看网站大全| 极品嫩模大尺度av在线播放| 亚洲成人色区| 在线无码国产精品亚洲а∨| 一区二区亚洲熟女偷拍| 在线观看老湿视频福利| 中文字幕人妻偷伦在线视频| 国产精品一区二区午夜久久| 中文字幕影片免费人妻少妇 | 无码Av在线一区二区三区| 久久久精品国产老熟女| 婷婷色国产精品视频二区| 欧美丰满熟妇性xxxx| 国产成人+亚洲欧洲+综合| 国产自产21区激情综合一区| 日韩女同精品av在线观看| 777午夜精品免费观看| 亚洲中文字幕乱码| 久久婷婷色香五月综合激激情| 亚洲男人天堂黄色av| 国产乱人伦av在线a| 亚洲人成网站在线播放小说|