辛體系算法在波的傳播與振動問題中的應(yīng)用.pdf

1、隨著科技的不斷發(fā)展，各種工程問題逐漸趨向復(fù)雜化，非線性領(lǐng)域也引起越來越多的重視。因此得益于計算機技術(shù)發(fā)展的科學(xué)計算得到了飛速的發(fā)展?？茖W(xué)計算已經(jīng)同理論與實驗共同構(gòu)成當代科學(xué)研究的三大支柱。本博士論文在Hamilton體系框架下，研究了辛幾何算法在部分波傳播問題和振動問題中的應(yīng)用。辛幾何算法相對傳統(tǒng)算法有其獨特的優(yōu)越性，因為保守體系可用Hamilton體系的方法描述，其特點是保辛。保辛給出保守體系結(jié)構(gòu)最重要的特性。而對于某些非保守系統(tǒng)則也

2、可通過轉(zhuǎn)化為保守系統(tǒng)進行分析。針對某些波傳播問題和振動問題，在辛體系下應(yīng)用精細積分法、擴展的Wittrick-Williams算法、子結(jié)構(gòu)分析、界帶分析等理論方法提出了一系列算法。大量的數(shù)值算例表明，本文給出的辛幾何算法在效率上和精度上有較大的優(yōu)勢。本論文的主要研究工作如下： (1)將精細積分法應(yīng)用于二階橢圓函數(shù)的求解，給出了二階橢圓函數(shù)的精細積分算法。大量的算例驗證了算法的可行性和正確性，并通過奇點處理使算法進一步完善。通過與

3、其他的現(xiàn)有算法和軟件對比，發(fā)現(xiàn)橢圓函數(shù)的精細積分算法無論是在精度上、速度上還是在適用范圍上都優(yōu)于其他算法。隨后在分層介質(zhì)光波導(dǎo)問題的應(yīng)用中再次證明該算法在處理工程問題時存在很大的優(yōu)勢。 (2)利用電磁波導(dǎo)的辛體系理論分析并求解了分層介質(zhì)中光波導(dǎo)的問題。首先在辛體系下，對線性分層介質(zhì)波導(dǎo)問題進行了分析，指出其本征值問題的復(fù)雜性。進而將精細積分法與擴展的Wittrick-Williams算法相結(jié)合給出了一種高精度的求解方法。利用本征

4、值計數(shù)可以做到指定范圍準確求解而且決不丟根，這是本文算法的一大特點。在此基礎(chǔ)上還對非線性分層介質(zhì)光波導(dǎo)問題作了簡單的分析，并利用精細積分方法求解了非線性Kerr材料的分層介質(zhì)波導(dǎo)問題。 (3)在Lagrange坐標下，利用位移法提出了一個描述淺水波問題的新方程。該方程不同于傳統(tǒng)的Euler坐標系下的淺水波方程，位移法淺水波方程可以用其相應(yīng)的Lagrange函數(shù)、變分原理進行分析，沿分析結(jié)構(gòu)力學(xué)的路走，從而達到保辛。孤波現(xiàn)象在淺水

5、波問題中受到廣泛重視，位移法淺水波方程也可以得出孤波解。大量的數(shù)值算例發(fā)現(xiàn)：在淺水波的基本假設(shè)下，位移法孤波解與傳統(tǒng)的KdV方程的孤波解非常接近。本文還分析了兩種方法的特點與區(qū)別，指出位移法淺水波方程為利用分析結(jié)構(gòu)力學(xué)研究淺水波問題開辟了一條新路。 (4)本文研究了非線性Duffing方程的保辛數(shù)值算法，并指出利用時間有限元等數(shù)值方法求解時須注意的保辛性問題?？臻g有限元是自動保辛的，因為有限元是基于變分原理的；而從變分原理導(dǎo)出的

6、時間有限元矩陣也有對稱性，從而達到保辛。對于非線性振動問題的數(shù)值求解常用到攝動法，傳統(tǒng)攝動法采用加法攝動，無法實現(xiàn)保辛。本文在處理該問題時采用乘法攝動，從而保證了傳遞矩陣保持為辛型，實現(xiàn)了保辛。通過和傳統(tǒng)Runge-Kutt法的大量數(shù)值對比，驗證了本文算法的正確性和穩(wěn)定性。 (5)對于工程問題中廣泛存在的跨域影響問題，將子結(jié)構(gòu)間的分界面擴展為有一定寬度的分界帶從而建立了界帶結(jié)構(gòu)，并給出了相應(yīng)的子結(jié)構(gòu)界帶分析和色散分析的理論和算法