幾何力學(xué)作為基于動力學(xué)系統(tǒng)的對稱性(或?qū)ΨQ性破缺)與保守(或非保守)性質(zhì)之間映射關(guān)系的新興學(xué)科,主要任務(wù)是發(fā)展盡可能多地保持動力學(xué)系統(tǒng)的固有(保守或非保守)性質(zhì)的數(shù)值方法。
1687年,隨著《自然哲學(xué)的數(shù)學(xué)原理》(Philosophiae Naturalis Principia Mathematica)一書的出版,牛頓力學(xué)宣告誕生,其最重要的貢獻(xiàn)之一是建立了力與運動之間微分形式的解析關(guān)系。之后,拉格朗日力學(xué)和哈密頓力學(xué)兩個基本框架相繼為動力學(xué)問題的描述提供了新的描述范式。在拉格朗日力學(xué)框架下,動力學(xué)系統(tǒng)優(yōu)美的數(shù)學(xué)對稱性被忽略了,因此,以拉格朗日力學(xué)為基礎(chǔ)發(fā)展的數(shù)值方法不具備保結(jié)構(gòu)特性。在動力學(xué)系統(tǒng)中引入對偶變量,可將動力學(xué)系統(tǒng)在哈密頓規(guī)范形式下進(jìn)行重新表述。從根本上說,牛頓力學(xué)在拉格朗日框架和哈密頓框架中的數(shù)學(xué)表達(dá)是完全等價的。但正如馮康先生所理解的,對于同一個動力學(xué)問題,由等價的數(shù)學(xué)表述形式得到的分析結(jié)果并不一定完全等價。
長期以來,經(jīng)典的數(shù)值方法,無論是著名數(shù)學(xué)家歐拉(Leonhard Euler)提出的歐拉差分格式,還是計算數(shù)學(xué)中廣泛使用的龍格庫塔(RungeKutta)法,共同的目標(biāo)都是提高數(shù)值方法的精度。然而在這個過程中,數(shù)學(xué)模型所描述的力學(xué)系統(tǒng)的幾何特性卻被忽略了。從根本上講,描述動力學(xué)系統(tǒng)的微分方程是連續(xù)的,而其數(shù)值方法對應(yīng)的系統(tǒng)是離散的;诖擞^點,數(shù)值模擬理應(yīng)在力學(xué)系統(tǒng)的同一幾何框架下進(jìn)行,并盡可能地保持原力學(xué)系統(tǒng)的定性性質(zhì),以提高數(shù)值解的長期數(shù)值穩(wěn)定性。
在忽略所有耗散效應(yīng)的情況下,一切物理過程都可以表述為能量守恒的哈密頓形式。這一結(jié)論的提出,一方面哈密頓力學(xué)的重要地位得到了提升; 另一方面也對哈密頓系統(tǒng)的數(shù)值分析提出了更高的要求,即對哈密頓系統(tǒng)的數(shù)值分析結(jié)果應(yīng)能夠再現(xiàn)其幾何性質(zhì),包括第一積分、辛結(jié)構(gòu)以及能量守恒定律。20世紀(jì)80年代,馮康提出了有限維哈密頓系統(tǒng)的辛方法,開創(chuàng)了計算幾何力學(xué)(又稱為保結(jié)構(gòu)方法)。在過去的半個世紀(jì)里,辛方法的保結(jié)構(gòu)特性已經(jīng)被許多研究者所報道。為了研究連續(xù)力學(xué)系統(tǒng)的局部動力學(xué)行為,Thomas J.Bridges和Jerrold E.Marsden將有限維哈密頓系統(tǒng)的辛方法推廣至無限維哈密頓系統(tǒng)的多辛方法。其重要貢獻(xiàn)在于,對無限維哈密頓系統(tǒng)構(gòu)造了時空聯(lián)合辛結(jié)構(gòu),稱為多辛結(jié)構(gòu); 并證明了在保持多辛結(jié)構(gòu)的情況下,數(shù)值分析可以很好地再現(xiàn)無限維哈密頓系統(tǒng)的局部動力學(xué)行為。
回顧哈密頓系統(tǒng)的定義不難發(fā)現(xiàn),哈密頓力學(xué)的局限性在于忽視了力學(xué)系統(tǒng)的耗散效應(yīng),這意味著基于哈密頓系統(tǒng)的分析方法不能用于解決實際工程中存在各種耗散效應(yīng)的力學(xué)問題;诖,本書作者提出了廣義多辛框架來解決這一問題,這正是本書的主要理論貢獻(xiàn)。廣義多辛理論框架在幾何力學(xué)和工程問題之間架起了一座橋梁,基于這一紐帶,本書介紹了大量廣義多辛分析方法在實際工程問題中的應(yīng)用實例。
本書章節(jié)結(jié)構(gòu)如下。第1章簡要介紹了幾何力學(xué)的起源和發(fā)展。第2章介紹了辛方法的數(shù)學(xué)基礎(chǔ)和幾個辛方法的例子。第3章中回顧了多辛理論,并介紹了多辛方法在無限維哈密頓系統(tǒng)中的若干應(yīng)用。第4章將保守系統(tǒng)保結(jié)構(gòu)思想推廣到非保守系統(tǒng),介紹了廣義多辛積分方法及其在波傳播問題中的應(yīng)用。第5~7章介紹了保結(jié)構(gòu)方法,包括辛方法、多辛方法、廣義多辛方法和復(fù)合保結(jié)構(gòu)方法在沖擊動力學(xué)問題、微/納米動力學(xué)系統(tǒng)和航天動力學(xué)系統(tǒng)中的應(yīng)用。當(dāng)然,限于篇幅,本書并未包含幾何力學(xué)的全部理論和應(yīng)用。
本書的主要內(nèi)容是基于作者團(tuán)隊近20年的研究成果整理而成。為使本書的內(nèi)容更加系統(tǒng),本書對其他著作和論文闡述的幾何力學(xué)相關(guān)基礎(chǔ)知識也進(jìn)行了綜述。需要說明的是,雖然在與本書相關(guān)的論文中,中國兵器科學(xué)研究院肖川研究員未被列為合著者,但為了感謝他在將保結(jié)構(gòu)方法應(yīng)用于沖擊動力學(xué)問題方面的重要貢獻(xiàn),將他列為本書的合著者。
感謝大連理工大學(xué)鐘萬勰院士、浙江大學(xué)朱位秋院士、德國錫根大學(xué)張傳增院士、中國科學(xué)院崔俊芝院士、利物浦大學(xué)歐陽華江教授、加拿大英屬哥倫比亞大學(xué)James J.Feng教授、英國薩里大學(xué)Thomas J.Bridges教授、美國圣澤維爾大學(xué)AbdulMajid Wazwaz教授、清華大學(xué)馮西橋教授、香港城市大學(xué)林志華教授、美國得克薩斯大學(xué)里約熱內(nèi)盧格蘭德谷分校喬志軍教授、大連理工大學(xué)彭海軍教授、河北經(jīng)貿(mào)大學(xué)王崗偉博士、日本慶應(yīng)大學(xué)彭林玉博士和美國中佛羅里達(dá)大學(xué)Brian Moore博士在本書的編寫工作中給予的無私幫助。
本書的出版受到西安理工大學(xué)西北旱區(qū)生態(tài)水利國家重點實驗室資助。同時感謝國家自然科學(xué)基金(12172281、11972284、11672241、11432010、11872303、11372253、11002115)、陜西省杰出青年科學(xué)基金(2019JC29)、軍科委基礎(chǔ)加強(qiáng)計劃173基金(2021JCJQJJ0565)、陜西省科技創(chuàng)新團(tuán)隊(2022TD61)和陜西高校青年教師創(chuàng)新團(tuán)隊的資助。
特別感謝陜西省科技創(chuàng)新團(tuán)隊(先進(jìn)設(shè)備關(guān)鍵動力學(xué)與控制)和陜西高校青年教師創(chuàng)新團(tuán)隊(空間太陽能電站動力學(xué)與控制)成員在本書的出版中做出的貢獻(xiàn)。謹(jǐn)向師俊平教授、曹小杉教授、胡義鋒副教授、張帆副教授、王震博士、徐萌波博士、淮雨露博士、惠小健博士、章培軍博士和團(tuán)隊所有同學(xué)表示衷心的感謝,謝謝大家無私的幫助和鼓勵。
限于作者的知識,書中的錯漏之處難以避免,其錯漏之處將在以后的工作中得以訂正。
西安理工大學(xué)胡偉鵬
中國兵器科學(xué)研究院肖川
西北工業(yè)大學(xué)鄧子辰
中國西安
2022年1月
第1章緒論
1.1幾何力學(xué)的生命力
1.1.1從線性諧振子的歐拉法開始
1.1.2數(shù)學(xué)擺模型St-rmer-Verlet格式的探討與改進(jìn)
1.2從拉格朗日力學(xué)到哈密頓力學(xué)
1.2.1拉格朗日力學(xué)
1.2.2哈密頓力學(xué)
1.3幾何力學(xué)的靈魂幾何積分
參考文獻(xiàn)
第2章有限維系統(tǒng)的辛算法
2.1辛方法的數(shù)學(xué)基礎(chǔ)
2.2典型的辛離散化方法
2.2.1辛龍格庫塔法
2.2.2分裂離散方法
2.3辛方法在力學(xué)問題中的應(yīng)用
2.3.1起落架折疊和展開過程辛精細(xì)積分方法研究
2.3.2航天動力學(xué)問題的辛龍格-庫塔法
參考文獻(xiàn)
第3章無限維哈密頓系統(tǒng)的多辛方法
3.1波動方程的多辛描述
3.2多辛理論的數(shù)學(xué)基礎(chǔ)
3.2.1辛和逆辛的對合與可逆性
3.2.2動量與能量守恒性
3.2.3多辛結(jié)構(gòu)與多辛守恒律
3.2.4哈密頓泛函
3.2.5多辛理論的一個更普遍的描述
3.3典型的多辛離散方法
3.3.1顯式中點格式
3.3.2歐拉Box格式
3.4多辛方法在波傳播問題中的應(yīng)用
3.4.1膜自由振動方程的多辛分析方法
3.4.2廣義五階KdV方程的多辛方法
3.4.3廣義(2 1)維KdV-mKdV方程的多辛方法
3.4.4朗道-金茲堡-希格斯方程的多辛龍格-庫塔法
3.4.5廣義波希尼斯克方程的多辛方法
3.4.6(2 1)維波希尼斯克方程孤立波共振的多辛模擬方法
3.4.7準(zhǔn)Degasperis-Procesi方程peakon-antipeakon碰撞的多辛模擬方法
3.4.8對數(shù)KdV方程高斯孤立波解的多辛分析
參考文獻(xiàn)
第4章非保守系統(tǒng)的動力學(xué)對稱破缺和廣義多辛方法
4.1動力學(xué)對稱破缺簡介
4.2從多辛積分到廣義多辛積分
4.3無限維動力學(xué)系統(tǒng)的對稱破缺
4.4廣義多辛分析方法在波傳播中的保結(jié)構(gòu)性質(zhì)初探
4.4.1關(guān)注伯格斯方程局部守恒性質(zhì)的隱式差分格式
4.4.2KdV-伯格斯方程中的幾何色散與黏性耗散的競爭關(guān)系
4.4.3復(fù)合KdV-伯格斯方程的廣義多辛離散化
4.4.4周期擾動下具有弱線性阻尼的非線性薛定諤方程近似保結(jié)構(gòu)分析
參考文獻(xiàn)
第5章沖擊動力學(xué)系統(tǒng)的保結(jié)構(gòu)分析方法
5.1沖擊動力學(xué)研究進(jìn)展介紹
5.1.1受軸向沖擊的柱和殼
5.1.2橫向沖擊載荷作用下的梁和板
5.1.3沖擊或爆炸載荷作用下的夾層結(jié)構(gòu)
5.1.4沖擊載荷下的多孔材料
5.2脈沖爆震發(fā)動機(jī)中燃料黏度引起的能量損失
5.3沖擊作用下非均勻中心對稱阻尼板內(nèi)的波傳播問題
5.4沖擊作用下非均勻非對稱圓板內(nèi)的波傳播問題
參考文獻(xiàn)
第6章微納米動力學(xué)系統(tǒng)的保結(jié)構(gòu)分析
6.1嵌入式單壁碳納米管中的混沌現(xiàn)象
6.2阻尼懸臂單壁碳納米管振蕩器的能量耗散
6.3嵌入式載流單壁碳納米管的混沌特性
6.4彈性約束的單壁碳納米管的混沌特性
6.5嵌入式單壁碳納米管軸向動力學(xué)屈曲的復(fù)合保結(jié)構(gòu)分析方法
參考文獻(xiàn)
第7章航天動力學(xué)系統(tǒng)的保結(jié)構(gòu)分析
7.1空間柔性阻尼梁的耦合動力學(xué)行為研究
7.2非球攝動下空間柔性阻尼梁動力學(xué)行為
7.3空間柔性梁所需的最小振動控制能量問題
7.4空間在軌繩系系統(tǒng)的能量耗散/轉(zhuǎn)移與穩(wěn)定姿態(tài)
7.5空間繩系系統(tǒng)中柔性梁的內(nèi)共振現(xiàn)象
7.6中心剛體-主動伸長柔性梁系統(tǒng)的耦合動力學(xué)行為
7.7由四根彈簧單邊約束的空間柔性阻尼板內(nèi)的彈性波傳播特性研究
參考文獻(xiàn)