哈密頓體系中波的傳播、魯棒控制與辛方法探索.pdf

1、本文在對(duì)偶哈密頓體系的統(tǒng)一框架下，探討了波的傳播、魯棒控制與保守動(dòng)力系統(tǒng)中的辛幾何求解理論和高效高精度的數(shù)值求解方法。主要工作如下： (1)本文將各向同性分層半無限空間中表面波問題導(dǎo)向?qū)ε俭w系，基于狀態(tài)方程矩陣是哈密頓矩陣的性質(zhì)，提出了一種求解分層半無限空間中Love和Rayleigh表面波問題的精確、高效方法。此方法聯(lián)合應(yīng)用了精細(xì)積分法和擴(kuò)展W-W算法。精細(xì)積分法是一個(gè)求解一階線性常微分方程組的精確數(shù)值方法，既可以應(yīng)用于空間域

2、的兩點(diǎn)邊值問題，又可以用于時(shí)間域的初值問題。它可以給出計(jì)算機(jī)精度意義下的精確解，并且保持哈密頓體系的辛結(jié)構(gòu)。應(yīng)用擴(kuò)展的W-W算法求解從精細(xì)積分法得到的超越方程時(shí)，對(duì)任意給定的試驗(yàn)頻率ω給出本征值計(jì)數(shù)J(ω)，通過計(jì)算試驗(yàn)頻率ω處的本征值計(jì)數(shù)可以確保不會(huì)遺漏任何本征值。 (2)各向異性介質(zhì)有21個(gè)獨(dú)立的彈性參數(shù)，并且三個(gè)方向的位移分量都耦合在一起，這使得各向異性介質(zhì)中波的傳播性質(zhì)變得異常復(fù)雜。本文將各向異性分層介質(zhì)中的波傳播問題導(dǎo)

3、向?qū)ε俭w系，基于狀態(tài)方程是哈密頓矩陣的性質(zhì)，進(jìn)一步將保辛的精細(xì)積分法和擴(kuò)展W-W算法推廣應(yīng)用于各向異性分層介質(zhì)的波傳播問題，該方法的有效性和高精確性得益于精細(xì)積分法對(duì)哈密頓系統(tǒng)的辛結(jié)構(gòu)的保持和本征值計(jì)數(shù)方法的應(yīng)用。應(yīng)用本文方法求解各向異性分層介質(zhì)的波傳播問題具有與求解各向同性介質(zhì)波傳播問題同樣簡(jiǎn)潔的公式和統(tǒng)一的計(jì)算步驟。 (3)本文將粘彈性分層介質(zhì)中的平穩(wěn)和非平穩(wěn)隨機(jī)波傳播問題導(dǎo)入到對(duì)偶體系，將隨機(jī)性問題轉(zhuǎn)換為確定性問題，通過變

4、換將粘彈性介質(zhì)中的控制方程轉(zhuǎn)換為哈密頓正則方程，提出了平穩(wěn)和非平穩(wěn)隨機(jī)波傳播問題的虛擬激勵(lì)－保辛精細(xì)積分法。首先，虛擬激勵(lì)法將隨機(jī)波傳播問題轉(zhuǎn)換為波數(shù)－頻率域中的確定性問題進(jìn)行求解。然后證明，雖然巖層是粘彈性介質(zhì)，存在阻尼或者能量耗散機(jī)制，但是將控制方程導(dǎo)入到對(duì)偶體系后仍然是哈密頓正則方程。在對(duì)偶體系下，波數(shù)－頻率域中的狀態(tài)方程是線性哈密頓方程的兩端邊值問題，保辛的精細(xì)積分法是求解此類問題的高精度算法，可以給出計(jì)算機(jī)意義下的精確解。

5、 (4)本文提出了隨機(jī)地震波作用下分層半無限空間與上部結(jié)構(gòu)相互作用問題的高效、精確方法。目前，對(duì)于巖層－結(jié)構(gòu)相互作用的研究都采用確定性激勵(lì)，但地震激勵(lì)本質(zhì)上是隨機(jī)的。本文中將虛擬激勵(lì)方法推廣應(yīng)用于隨機(jī)性巖層－結(jié)構(gòu)相互作用問題。應(yīng)用虛擬激勵(lì)法后，所有控制方程都轉(zhuǎn)換為確定性方程，這時(shí)可以應(yīng)用精細(xì)積分法來積分波數(shù)－頻率域中的控制方程，得到計(jì)算機(jī)意義下的精確解。 (5)本文提出了哈密頓體系下大系統(tǒng)分散魯棒控制的廣義模態(tài)綜合法。在大系

6、統(tǒng)控制中，分散控制方法通過將整個(gè)系統(tǒng)劃分成一系列子系統(tǒng)來研究系統(tǒng)的控制問題，這與結(jié)構(gòu)力學(xué)中的子結(jié)構(gòu)分析技術(shù)類似。本文表明分散控制的提法與結(jié)構(gòu)力學(xué)中的子結(jié)構(gòu)分析相互模擬：靜力子結(jié)構(gòu)分析對(duì)應(yīng)于通常意義的分散控制理論，而動(dòng)力子結(jié)構(gòu)分析則對(duì)應(yīng)于分散H∞控制理論。本文由此建立了H∞分散控制與子結(jié)構(gòu)振動(dòng)問題之間的對(duì)應(yīng)關(guān)系，指出大系統(tǒng)分散控制的最優(yōu)參數(shù)γcr-2對(duì)應(yīng)于整體結(jié)構(gòu)的一階本征值，并利用模態(tài)綜合法計(jì)算出H∞分散控制問題中的最優(yōu)參數(shù)γcr-2。

7、給出了大系統(tǒng)分散魯棒控制廣義模態(tài)綜合法的正交性、完備性和展開定理。 (6)已有哈密頓系統(tǒng)的保辛算法多是基于差分近似的Runge-Kutta方法。保辛Runge-Kutta法的困難之處在于構(gòu)造高精度算法。本文基于分析結(jié)構(gòu)力學(xué)，利用時(shí)間有限元的概念，將廣義位移用形函數(shù)離散，得到近似作用量。利用近似作用量作為生成函數(shù)，構(gòu)造無約束和有約束保守系統(tǒng)的保辛算法。利用多個(gè)節(jié)點(diǎn)構(gòu)造高階形函數(shù)，能顯著提高有限元保辛算法的精度。最終通過天體力學(xué)、分