物質(zhì)點(diǎn)法改進(jìn)算法及其工程應(yīng)用
定 價(jià)：65 元
- 作者：孫政，周曉敏著
- 出版時(shí)間：2021/11/1
- ISBN：9787548746614
- 出 版 社：中南大學(xué)出版社
- 中圖法分類：O302
- 頁(yè)碼：193頁(yè)
- 紙張：膠版紙
- 版次：1
- 開本：16開
本書主要研究物質(zhì)點(diǎn)法的算法和擴(kuò)展其工程應(yīng)用。在算法方面：首先，基于粒子子域積分代替物質(zhì)點(diǎn)積分，提出粒子子域積分物質(zhì)點(diǎn)法，消除物質(zhì)點(diǎn)法的積分誤差，提高計(jì)算精度和計(jì)算收斂性；其次，基于B樣條基函數(shù)代替物質(zhì)點(diǎn)法的線性插值形函數(shù)，歸結(jié)了B樣條物質(zhì)點(diǎn)法的基本概念，對(duì)B樣條物質(zhì)點(diǎn)的計(jì)算精度、計(jì)算收斂性和計(jì)算效率等方行了系統(tǒng)研究而，分別基于截?cái)鄬哟蜝樣條、局部加密B樣條和橋域法，提出了B樣條物質(zhì)點(diǎn)法的背景網(wǎng)格局部細(xì)化算法，通過具體算例驗(yàn)證算法有效性。在工程應(yīng)用方面：一是將蒙特卡洛模擬和物質(zhì)點(diǎn)法相耦合，發(fā)展和提出隨機(jī)物質(zhì)點(diǎn)法，并對(duì)土質(zhì)滑坡問行定量風(fēng)險(xiǎn)評(píng)估；二是基于人工狀態(tài)方程，發(fā)展和提出弱可壓物質(zhì)點(diǎn)法，并應(yīng)用于求解牛頓、非牛頓流體的流動(dòng)問題以及流固耦合問題。
數(shù)值模擬方法在工程應(yīng)用和科學(xué)研究等領(lǐng)域起著施來趨重要的作用。然而，傳統(tǒng)的有網(wǎng)格類方法，如有限無法、有限差分法、有限體積法等，在處理大變形、材料破壞、穿孔、多相介質(zhì)界面?zhèn)鞑サ葐栴}時(shí)，往往需要借助網(wǎng)格重構(gòu)技術(shù)以避免網(wǎng)格畸變，相應(yīng)地將遺成算法復(fù)雜、精度下降、效率低甚至數(shù)值不穩(wěn)定等問題。不依賴網(wǎng)格的無網(wǎng)格類方法可很好地求解這類有網(wǎng)格方法難以求解的復(fù)雜問題。物質(zhì)點(diǎn)法是一種無網(wǎng)格粒子類方法，由流體流動(dòng)隱式粒子法演變而來，其來用拉格朗日粒子來離散求解域，同時(shí)采用歐拉背景網(wǎng)格求解系統(tǒng)的控制方程，既發(fā)揮了拉格朗日算法和歐拉算法的優(yōu)點(diǎn)，又避免了各自的不足。目前，物質(zhì)點(diǎn)法已在各類復(fù)雜問題的求解中得到了廣泛應(yīng)用，并展現(xiàn)了強(qiáng)大的求解優(yōu)勢(shì)。
傳統(tǒng)物質(zhì)點(diǎn)法采用線性插值函數(shù)來實(shí)現(xiàn)物質(zhì)點(diǎn)粒子與背景網(wǎng)格節(jié)點(diǎn)之間的信息映射，即物質(zhì)點(diǎn)法僅具有C連續(xù)性，因此當(dāng)物質(zhì)點(diǎn)粒子穿越背景網(wǎng)格邊界時(shí)將會(huì)產(chǎn)生所謂的“網(wǎng)格穿越誤差”，從而造成應(yīng)力的非物理振蕩和求解精度不足等問題。同時(shí)，在傳統(tǒng)物質(zhì)點(diǎn)法中，物質(zhì)點(diǎn)粒子既被用于離散求解域，又被作為積分點(diǎn)求積節(jié)點(diǎn)內(nèi)力，其中加權(quán)函數(shù)為物質(zhì)點(diǎn)所占的體積域。因此，在求解大變形問題時(shí)，隨著物體的運(yùn)動(dòng)和變形，物質(zhì)點(diǎn)的位置和體積變化較大，從而帶來積分誤差。如何消除傳統(tǒng)物質(zhì)點(diǎn)法的網(wǎng)格穿越誤差和積分誤差，從而提高其計(jì)算精度以一步擴(kuò)展其工程應(yīng)用，是當(dāng)前物質(zhì)點(diǎn)法研究的熱點(diǎn)問題。
本書圍繞物質(zhì)點(diǎn)法的算法及其工程應(yīng)用的擴(kuò)展兩個(gè)方面展開論述。在算法方面，首先，基于粒子子域積分代替物質(zhì)點(diǎn)積分，提出粒子子域積分物質(zhì)點(diǎn)法，消除物質(zhì)點(diǎn)法的積分誤差，提高計(jì)算精度和計(jì)算收斂性；其次，基于B樣條基函數(shù)代替物質(zhì)點(diǎn)法的線性插值函數(shù)，發(fā)展和提出了B樣條物質(zhì)點(diǎn)法，并對(duì)B樣條物質(zhì)點(diǎn)的計(jì)算精度、計(jì)算收斂性和計(jì)算效率等方行了系統(tǒng)研究；后，分別基于藏?cái)鄬哟蜝樣條、局部加密B樣條和橋城法，提出了B樣條物質(zhì)點(diǎn)法的背景網(wǎng)格局部細(xì)化算法，通過具休算例驗(yàn)證算法的有效性。在工程應(yīng)用方面，一是將蒙特卡洛模擬和物質(zhì)點(diǎn)法相耦合，發(fā)展和提出了隨機(jī)物質(zhì)點(diǎn)法，并對(duì)土質(zhì)滑坡問行了定量風(fēng)險(xiǎn)評(píng)估；二是基于人工狀態(tài)方程，發(fā)展和提出了弱可壓物質(zhì)點(diǎn)法，并應(yīng)用于求解牛頓、非牛頓流體的流動(dòng)問題以及流固耦合問題。
本書由江西理工大學(xué)孫政和周曉敏共同撰寫完成，其中孫政撰寫本書的、三、四、五章節(jié)，共約12萬(wàn)字，周曉敏撰寫本書的第二、六、七章共約10萬(wàn)字。本書涉及的研究得到了國(guó)家自然科學(xué)基金青年項(xiàng)目（11902127）、江西省自然科學(xué)基金青年項(xiàng)目（20192BAB212010）和江西省教育廳科技項(xiàng)目（GJJ20083pan>）等的支持，江西理工大學(xué)博士啟動(dòng)基金（JXXJBS18042）和江西理工大學(xué)清江學(xué)術(shù)文庫(kù)為本書的出版提供了經(jīng)費(fèi)支持，在此一并表示感謝。向在寫作過程中參考的國(guó)內(nèi)外文獻(xiàn)作者一并表示感謝。感謝筆者的家人和朋友在本人工作、學(xué)習(xí)中給予的關(guān)心和支持。限于筆者的學(xué)識(shí)，書中難免存在疏漏之處，懇請(qǐng)讀者和學(xué)界同行不吝指出。
1緒論 /001
1.1研完背景和意義/ 001
1.2物質(zhì)點(diǎn)法研究現(xiàn)狀 1002
1.3數(shù)值模擬在巖土工程問題中的應(yīng)用1009
1.4數(shù)值模擬在流圖耦合問題中的應(yīng)用1012
1.5本書主要研究?jī)?nèi)容1015
1.6本章小結(jié)1017
2物質(zhì)點(diǎn)法基本理論及現(xiàn)有算法 /018
2.1控制方程/019
2.2物質(zhì)點(diǎn)離散/023
2.3算法實(shí)現(xiàn)1027
2.4現(xiàn)有物質(zhì)點(diǎn)算法概述1030
2.5本章小結(jié)/037
3基于粒子子域積分的物質(zhì)點(diǎn)法 /038
3.1引言/038
3.2物質(zhì)點(diǎn)法的內(nèi)力積分誤差1038
3.3高斯積分簡(jiǎn)介/040
3.4粒子子域積分物質(zhì)點(diǎn)法/041
3.5數(shù)值算例/045
3.6本章小結(jié)/054
4基于B樣條基函數(shù)的物質(zhì)點(diǎn)法/055
4.1 引言 / 055
4.2B樣條基函敲/055
4.3 B樣條物質(zhì)點(diǎn)法 /059
4.4數(shù)值算例/062
4.5本章小結(jié)/080
5B樣條物質(zhì)點(diǎn)法背景網(wǎng)格局部細(xì)化算法 /082
5.1引言/082
5.2基于THB和LRB的BSMPM背景網(wǎng)格局部細(xì)化算法 / O83
5.3基于橋城法的BSMPM背景網(wǎng)格局部細(xì)化算法 /O87
5.4 數(shù)值算例/090
5.5本章小結(jié)/104
6蒙特卡洛物質(zhì)點(diǎn)法在巖土工程問題中的應(yīng)用/106
6.1 引言 /106
6.2蒙特卡洛模擬/107
6.3 隨機(jī)場(chǎng)理論 / 109
6.4蒙特卡洛物質(zhì)點(diǎn)法/113
6.5土質(zhì)滑坡及其定量風(fēng)險(xiǎn)評(píng)估/114
6.6本章小結(jié)/125
7弱可壓物質(zhì)點(diǎn)法在流固耦合問題中的應(yīng)用 /126
7.1 引言 / 126
7.2弱可壓物質(zhì)點(diǎn)法/127
7.3牛頓流體流動(dòng)問題/129
7.4非牛頓流體流動(dòng)問題/141
7.5牛頓流體與彈性體的相互耦合/149
7.6 非牛頓泥石流沖擊剛性障礙物 /161
7.7本章小結(jié)/166
參考文獻(xiàn) /167
目前，數(shù)值計(jì)算已廣泛應(yīng)用于求解各類科學(xué)問題和工程技術(shù)問題，與理論分析和試驗(yàn)研究構(gòu)成了現(xiàn)代科學(xué)技術(shù)的三大支柱，并具有快捷、和低成本的優(yōu)勢(shì)，同時(shí)可以模擬復(fù)雜材料、復(fù)雜結(jié)構(gòu)及復(fù)雜的變形和運(yùn)動(dòng)。
數(shù)值計(jì)算方法中，描述運(yùn)動(dòng)的主要形式有拉格朗日法和歐拉法。在拉格朗日法中，求解域與計(jì)算網(wǎng)格相固連，材料和網(wǎng)格之間不存在相對(duì)運(yùn)動(dòng)，控制方程中不存在對(duì)流項(xiàng)簡(jiǎn)化了控制方程的建立和求解，便于追蹤物質(zhì)信息間和空間的變化，同時(shí)易于處理與變形歷史相關(guān)的材料本構(gòu)模型，計(jì)算固體力學(xué)多采用拉格朗日法；但當(dāng)涉及材料特大變形，如破碎、卷曲等強(qiáng)非線性現(xiàn)象時(shí)，由于網(wǎng)格畸變的影響，傳統(tǒng)的拉格朗日網(wǎng)格法將難以適用。在歐拉法中，計(jì)算網(wǎng)格固定于計(jì)算空間內(nèi)，不隨物體的運(yùn)動(dòng)而運(yùn)動(dòng)，通過計(jì)算網(wǎng)格邊界上質(zhì)量、動(dòng)量和能量的通量得到所求解問題的各物理量的空間分布，計(jì)算網(wǎng)格在計(jì)算過程中保持不變，因此不存在網(wǎng)格畸變問題，在計(jì)算流體力學(xué)中多采用歐拉法，但歐拉法追蹤的是網(wǎng)格邊界上質(zhì)量、動(dòng)量和能量的通量流動(dòng)，需求解非線性對(duì)流項(xiàng)，增加了求解難度，且難以追蹤自由表面和材料交界面，同時(shí)不易于追蹤各質(zhì)點(diǎn)的運(yùn)動(dòng)時(shí)間歷程。
將拉格朗日法和歐拉法有機(jī)結(jié)合，發(fā)揮各自的優(yōu)點(diǎn)、克服各自的缺點(diǎn)，則可求解一系列單一拉格朗日法或歐拉法難以求解的工程問題，物質(zhì)點(diǎn)法正是一種結(jié)合了拉格朗日法和歐拉法優(yōu)勢(shì)，并避免了各自不足的無網(wǎng)格粒子型方法。物質(zhì)點(diǎn)法，首先，將問題域離散成一系列拉格朗日質(zhì)點(diǎn)，質(zhì)點(diǎn)攜帶材料區(qū)域的物質(zhì)信息，并跟隨物體運(yùn)動(dòng)而運(yùn)動(dòng)，便于追蹤物質(zhì)信息間和空間的變化，同時(shí)易于處理與變形歷史相關(guān)的材料本構(gòu)模型，且不存在網(wǎng)格畸變問題；其次，采用歐拉背景網(wǎng)行空間導(dǎo)數(shù)和動(dòng)量方程的求解，以及實(shí)現(xiàn)各質(zhì)點(diǎn)間的相互作用和聯(lián)系，減少了粒子搜索算法的耗時(shí)，同時(shí)易于本質(zhì)邊界條件的施加。物質(zhì)點(diǎn)法自提出之目起，已在機(jī)高速碰撞，沖擊侵徹，牌堆沖擊，由體滑坡、流團(tuán)耦合等一系列復(fù)雜同題中得到了廣泛應(yīng)用，并展現(xiàn)了強(qiáng)大的求解能力和求解優(yōu)勢(shì)。
但受限于線性插值形雨數(shù)影啊，物質(zhì)點(diǎn)法僅具有C”連續(xù)性，當(dāng)質(zhì)點(diǎn)算越背景網(wǎng)格邊界時(shí)將會(huì)產(chǎn)生所謂的“網(wǎng)格穿穩(wěn)誤差”，造成求解精度的降低。在物質(zhì)點(diǎn)法中，物質(zhì)點(diǎn)同時(shí)也被作為積分點(diǎn)求解節(jié)點(diǎn)內(nèi)力，其加權(quán)函數(shù)為物質(zhì)點(diǎn)所古的體積城。在計(jì)算過程中，尤其涉及大變形同題時(shí)，隨著物質(zhì)點(diǎn)的運(yùn)動(dòng)和變形，物質(zhì)點(diǎn)體積城變化較大，從麗帶來積分誤差。因此如何消除和網(wǎng)格穿越誤差，減小內(nèi)力積分誤差一步發(fā)展高精度物質(zhì)點(diǎn)算法和擴(kuò)展其工程應(yīng)用一直是計(jì)算力學(xué)研究領(lǐng)域的熱點(diǎn)問題之一。
1.2物質(zhì)點(diǎn)法研究現(xiàn)狀
物質(zhì)點(diǎn)法是Sulsky（蘇爾斯凱）等提出的一種粒子型無網(wǎng)格算法，其可追溯到流體力學(xué)領(lǐng)城中的質(zhì)點(diǎn)網(wǎng)格法（partical in cell，PIC）和流體隱式粒子法（luid implieit paricle method,FLIP)。Sulsky等通過流體隱式粒子法，將其運(yùn)用到求解固體力學(xué)問題上，并命名為物質(zhì)點(diǎn)法（material pointmethod，MPM），其主要有以下幾個(gè)方面：所有的物理信息及本構(gòu)方程分別由質(zhì)點(diǎn)攜帶和計(jì)算，因此可方便處理與歷史相關(guān)的本構(gòu)材料模型；基于質(zhì)點(diǎn)離散，采用等效積分弱形式，重新推導(dǎo)了動(dòng)量方程的離散格式；采用顯式時(shí)間積分對(duì)控制方行求解。
物質(zhì)點(diǎn)法隸屬于伽遼金型無網(wǎng)格法，其采用有限元法中的線性函數(shù)作為插值函數(shù)，采用拉格朗日和歐拉雙重描述，由一系列的拉格朗日質(zhì)點(diǎn)離散求解域，并攜帶材料的所有物質(zhì)信息括位置、速度、動(dòng)量、應(yīng)力應(yīng)變等歷史變量，跟隨物體運(yùn)動(dòng)而運(yùn)動(dòng)；采用歐拉背景網(wǎng)格覆蓋整個(gè)求解域，在背景網(wǎng)格行空間導(dǎo)數(shù)和控制方程的求解，并實(shí)現(xiàn)相鄰物質(zhì)點(diǎn)之間的相互作用與聯(lián)系。
……
書單推薦
