《辛破繭:辛拓展新層次(增訂版)》通過與實際課題銜接時的典型例題以及初步探討的成果,大體上只是給出思路。至于更一般、深入、詳盡的研究則有待能人發(fā)揮了。
0 多維線性動力學的求解
0.1 線性系統(tǒng)的分離變量法與本征問題
0.2 傳遞辛矩陣的本征問題
一 離散系統(tǒng)的保辛-守恒算法
1.1 坐標變換的Jacobi矩陣
1.2 傳遞辛矩陣,Lagrange括號與Poisson括號
1.3 保辛一守恒的參變量算法
1.4 用辛矩陣乘法表述的正則變換
1.4.1 時不變正則變換的辛矩陣乘法表述
1.4.2 時變正則變換的辛矩陣乘法表述
1.4.3 基于線性時不變系統(tǒng)的時變正則變換
1.4.4 包含時間坐標的正則變換
1.5 保辛-守恒的接觸參變量算法
1.6 保辛攝動多層網(wǎng)格法
1.6.1 多層次有限元
1.6.2 多層次的迭代求解
1.6.3 數(shù)值例題
1.7 傳遞辛矩陣群
二 不同時間的有限元離散
2.1 雙曲型偏微分方程的特征線理論概要
2.2 波動方程
2.3 變動邊界問題與混和元
2.4 剛性雙曲型偏微分方程例題
2.5 物理意義,Lorentz變換
三 不同維數(shù)的有限元離散
3.1 結構力學有限元自動保辛
3.2 波動偏微分方程,不同維數(shù)位錯的轉換
3.3 數(shù)值算例
3.4 辛數(shù)學能改革開放嗎?
3.5 接觸問題
3.5.1 拉壓模量不同材料的參變量變分原理和有限元方法〔24〕
3.5.2 拉壓不同剛度桁架的動力參變量保辛方法〔25〕
3.6 本章結束語
四 界帶與時滯
4.1 結構力學的界帶分析〔30,33〕
4.1.1 結構力學的界帶理論與能帶分析
4.1.2 界帶分析的能量變分法
4.1.3 色散關系
4.1.4 子結構界帶分析
4.1.5 不同原子組成周期鏈的數(shù)值分析
4.1.6 無限長多排原子鏈組合的情況
4.2 時滯與界帶
4.2.1 離散-維鏈系統(tǒng)的模擬
4.2.2 逐步前進的算法
4.3 連續(xù)系統(tǒng)的能量形式
4.3.1 連續(xù)系統(tǒng)動力學的能量形式
五 結束語
附錄
附錄1 SiPESC構造的簡單介紹
附錄2 力學具有基礎與應用學科的兩重性
參考文獻
不同維數(shù)的有限元離散
美國的總統(tǒng)信息技術顧問委員會給總統(tǒng)的報告會強調(diào):“計算-科學同理論和實驗并列,已成為科學事業(yè)的第2支柱!
既然要考慮計算科學,有限元法是自然的選擇。回看結構力學的有限元法,有五花八門的網(wǎng)格自動生成,而根本沒有恒定維數(shù)下同時間離散的限制。動力學在有限元離散方面應力求與結構力學有限元法融合,第三種局限性也要破繭。分析動力學常微分方程組有恒定維數(shù)的限制。要發(fā)展到偏微分方程,采用各種離散于段進行求解是自然的。有限元法提供了思路,可考慮不同維數(shù)的離散。
為清楚起見,這里只講空間一維、時間一維,即時一空2維問題。這樣好講些。本書只求破繭,不求做出全面的推進,以打開思路為目標。因此仍以波動偏微分方程為對象進行分析。
然而,對空間坐標與時間坐標分別離散而生成的離散時間空間格點,是規(guī)則的網(wǎng)格。是否能像結構力學有限元一樣,將兩種坐標混合在一起進行有限元離散呢?前面對于不同時間坐標進行了破繭。這里要考慮不同維數(shù)了。
一旦考慮不同維數(shù),群論就發(fā)生困難了。M.F.Atiyah說:“群在自然中產(chǎn)生,它們是使事物運動的東西,它們是變換或置換……。理解這些東西的本性,并且使用它們才是目的!庇终f:“重要的東西常常不是技術上最困難的即最難證明的東西,而常常是較為初等的部分。因為這些部分與其他領域、分支的相互作用最廣泛,即影響面最大!薄霸谌赫撝杏性S多極端重要的,并且在數(shù)學的各個角落到處都出現(xiàn)的東西。這些是較為初等的東西:群及其同態(tài),表示的基本觀點、一般的性質、一般的方法——這些才是真正重要的。辛矩陣群的乘法本來也有恒定維數(shù)的要求。但偏微分方程的離散求解,時一空混合有限元的網(wǎng)格離散不能拘泥于恒定維數(shù)。這樣,就不是單純的傳遞辛矩陣群了,當然也要破繭。在恒定維數(shù)的傳遞辛矩陣群之外,必定要有某些方法處理維數(shù)變化的護展。
為什么要保辛?這是為了保持保守系統(tǒng)的優(yōu)良性質。保守系統(tǒng)有變分原理,可保持其優(yōu)良性質,例如守恒性等。有限元法是從變分原理推導來的,就繼承了這些優(yōu)良性質,大家愿意用。傳遞辛矩陣群不敷應用需要,而要破繭之時,也應遵循變分原理的思路。
……