目錄
序
前言
第一部分 基礎(chǔ)篇
第1章 離散元法的原理和實(shí)現(xiàn) 3
1.1 離散元法應(yīng)用領(lǐng)域與數(shù)值計(jì)算軟件 3
1.2 離散元法的基本原理 4
1.2.1 單元的接觸模型 5
1.2.2 兩個(gè)不同單元的連接 6
1.2.3 時(shí)間步迭代算法 6
1.3 單元的阻尼簡(jiǎn)諧振動(dòng) 7
1.3.1 單元的阻尼力 7
1.3.2 阻尼簡(jiǎn)諧振動(dòng)方程 8
1.3.3 單元振動(dòng)周期與半徑的關(guān)系 8
1.3.4 單元運(yùn)動(dòng)過(guò)程的微分 9
1.3.5 阻尼作用和最優(yōu)阻尼系數(shù) 9
1.4 離散元系統(tǒng)的能量轉(zhuǎn)換與能量守恒 11
1.4.1 系統(tǒng)的機(jī)械能 11
1.4.2 熱量的計(jì)算 12
1.4.3 能量守恒與外力做功 14
1.5 離散元法的計(jì)算尺度和工程應(yīng)用可行性 14
1.6 離散元法三大問(wèn)題的認(rèn)識(shí)和解決 15
1.6.1 矩陣離散元和高性能矩陣計(jì)算 16
1.6.2 離散元法的定量化建模問(wèn)題 17
1.6.3 基于離散元法的多場(chǎng)耦合方法 19
1.7 矩陣離散元的發(fā)展和展望 20
第2章 MatDEM的基本結(jié)構(gòu) 22
2.1 軟件基礎(chǔ) 22
2.1.1 運(yùn)行環(huán)境和軟件安裝 22
2.1.2 MatDEM窗口程序 23
2.1.3 軟件運(yùn)行和中止 27
2.1.4 軟件應(yīng)用示例簡(jiǎn)介 29
2.1.5 幫助文件簡(jiǎn)介 31
2.2 MatDEM程序結(jié)構(gòu) 32
2.2.1 MatDEM軟件文件夾 32
2.2.2 MatDEM的層次結(jié)構(gòu) 32
2.2.3 MatDEM中主要的類 34
2.3 MatDEM中的單元類型 35
2.3.1 活動(dòng)單元、固定單元、虛單元 35
2.3.2 鎖定單元自由度 36
2.3.3 單元編號(hào)的規(guī)則 37
2.4 MatDEM的數(shù)據(jù)結(jié)構(gòu) 37
2.4.1 單元屬性數(shù)組 37
2.4.2 鄰居矩陣和連接信息矩陣 38
2.4.3 組的數(shù)據(jù)結(jié)構(gòu)和操作 40
2.5 單元的接觸模型 42
2.5.1 線彈性模型 42
2.5.2 赫茲接觸模型 43
2.6 彈性clump團(tuán)簇 44
2.6.1 彈性clump的原理 44
2.6.2 彈性clump的使用 45
第3章 幾何建模和材料設(shè)置 48
3.1 建立數(shù)值模擬箱 48
3.1.1 建立堆積模型 48
3.1.2 生成空箱子 51
3.2 模塊化建�；�礎(chǔ) 52
3.2.1 利用結(jié)構(gòu)體建模 52
3.2.2 單元過(guò)濾器和單元篩選 57
3.2.3 Tool_Cut和數(shù)字高程建模 59
3.3 基于數(shù)字圖像建模 60
3.3.1 導(dǎo)入黑白圖像切割模型 60
3.3.2 導(dǎo)入彩色圖像分塊建模 62
3.4 材料設(shè)置 66
3.4.1 材料設(shè)置概述 66
3.4.2 直接輸入材料性質(zhì) 66
3.4.3 自動(dòng)訓(xùn)練材料 67
3.5 平衡模型 70
3.5.1 迭代計(jì)算函數(shù)和標(biāo)準(zhǔn)平衡 70
3.5.2 強(qiáng)膠結(jié)平衡 71
3.5.3 鄰居單元檢索和零時(shí)平衡 72
3.5.4 模型平衡狀態(tài)判斷標(biāo)準(zhǔn) 73
3.6 裂隙和節(jié)理的設(shè)置 74
3.6.1 通過(guò)組來(lái)設(shè)置軟弱層和裂隙 75
3.6.2 利用三角面來(lái)設(shè)置裂隙和節(jié)理 76
3.6.3 使用Tool_Cut來(lái)設(shè)置裂隙和節(jié)理 76
3.6.4 連接過(guò)濾器的定義和使用 76
第4章 荷載設(shè)置和數(shù)值計(jì)算 78
4.1 數(shù)值計(jì)算初始化和設(shè)置 78
4.1.1 模型參數(shù)初始化 78
4.1.2 計(jì)算相關(guān)參數(shù)設(shè)置 78
4.2 邊界和荷載 79
4.2.1 邊界條件 79
4.2.2 應(yīng)力荷載 80
4.2.3 位移荷載 82
4.2.4 振動(dòng)荷載 82
4.2.5 其他荷載 84
4.3 迭代計(jì)算時(shí)間和設(shè)置 85
4.3.1 迭代計(jì)算的運(yùn)行時(shí)間 85
4.3.2 單元半徑和計(jì)算耗時(shí) 86
4.3.3 動(dòng)態(tài)問(wèn)題的設(shè)置 88
4.3.4 準(zhǔn)靜態(tài)問(wèn)題的設(shè)置和模擬精度 88
4.4 模擬參數(shù)的定義和修改 90
4.4.1 自定義參數(shù)的創(chuàng)建 90
4.4.2 組單元屬性的修改 91
4.4.3 時(shí)間步、單元半徑和剛度的設(shè)定 93
4.5 文件的壓縮、保存和讀取 94
4.5.1 文件的壓縮 94
4.5.2 文件的保存和讀取 95
第5章 后處理和系統(tǒng)函數(shù) 96
5.1 后處理窗口界面 96
5.1.1 后處理的主窗口 96
5.1.2 后處理圖件的保存 98
5.1.3 GIF動(dòng)畫制作窗口 99
5.2 后處理繪圖函數(shù) 101
5.2.1 d.show通用繪圖函數(shù) 101
5.2.2 切片顯示和過(guò)濾顯示 106
5.2.3 單元位置和受力顯示 108
5.3 數(shù)據(jù)處理和曲線繪制 108
5.3.1 利用保存的數(shù)據(jù)繪制曲線 108
5.3.2 利用自動(dòng)記錄數(shù)據(jù)繪制曲線 109
5.3.3 在MATLAB里處理數(shù)據(jù)和繪圖 111
5.4 系統(tǒng)函數(shù)與功能 112
5.4.1 GPU計(jì)算設(shè)置和狀態(tài)查看 112
5.4.2 函數(shù)的定義和運(yùn)行 114
5.4.3 代碼文件的批處理 115
5.4.4 隨機(jī)種子和隨機(jī)模型 116
5.4.5 計(jì)時(shí)函數(shù) 116
5.5 利用系統(tǒng)底層函數(shù)來(lái)建模 116
5.5.1 兩球碰撞過(guò)程的底層建模 117
5.5.2 利用曲線圖來(lái)分析碰撞過(guò)程 119
第二部分 實(shí) 踐 篇
第6章 巖土工程基礎(chǔ)應(yīng)用 123
6.1 樁土作用 123
6.1.1 堆積地層模型 123
6.1.2 建立樁土作用模型 125
6.1.3 拔樁過(guò)程的數(shù)值模擬 129
6.2 隧道建模 132
6.2.1 堆積地層模型 132
6.2.2 建立隧道模型 133
6.2.3 施加荷載和數(shù)值模擬 136
6.3 盾構(gòu)滾刀破巖 138
6.3.1 堆積地層模型 139
6.3.2 建立滾刀破巖模型 139
6.3.3 滾刀破巖數(shù)值模擬過(guò)程 143
6.3.4 提高滾刀破巖計(jì)算的速度 146
第7章 巖土體離散元試驗(yàn) 147
7.1 直剪和扭剪試驗(yàn) 147
7.1.1 定義試驗(yàn)參數(shù)和堆積試樣 147
7.1.2 制作剪切盒和切割試樣 149
7.1.3 將試樣放入剪切盒 152
7.1.4 材料設(shè)置和數(shù)值模擬 155
7.2 真三軸試驗(yàn)和節(jié)理建模 158
7.2.1 構(gòu)建真三軸試驗(yàn)箱 158
7.2.2 利用三角面和多邊形定義裂隙面 159
7.2.3 利用Tool_Cut定義復(fù)雜的節(jié)理面 161
7.2.4 施加真三軸應(yīng)力 165
第8章 復(fù)雜三維模型的建模 167
8.1 用數(shù)字高程定義三維層面 167
8.1.1 利用離散點(diǎn)獲得數(shù)字高程 167
8.1.2 層面的數(shù)字高程數(shù)據(jù)處理 169
8.2 創(chuàng)建薄殼模型 174
8.2.1 建立幾何模型 174
8.2.2 切割和堆積模型 176
8.3 建立三維邊坡模型 179
8.3.1 導(dǎo)入材料 180
8.3.2 設(shè)置地層的材料 180
8.4 滑坡運(yùn)動(dòng)過(guò)程數(shù)值模擬 182
8.4.1 數(shù)值模擬參數(shù)設(shè)置 183
8.4.2 迭代計(jì)算和模擬結(jié)果 183
第9章 動(dòng)力作用數(shù)值模擬 185
9.1 隕石撞擊地面 185
9.1.1 堆積地層模型 185
9.1.2 建立隕石模型 186
9.1.3 隕石撞擊過(guò)程模擬 189
9.2 礦山斜坡爆破 191
9.2.1 建立斜坡模型 191
9.2.2 設(shè)置爆破點(diǎn)和爆破能量 193
9.2.3 迭代計(jì)算和模擬結(jié)果 195
9.3 地震動(dòng)力作用 196
9.3.1 建立地形和分層 197
9.3.2 地震波產(chǎn)生和傳播過(guò)程 199
9.2.3 迭代計(jì)算和模擬結(jié)果 195
9.3 地震動(dòng)力作用 196
9.3.1 建立地形和分層 197
9.3.2 地震波產(chǎn)生和傳播過(guò)程 199
第10章 多場(chǎng)耦合數(shù)值模擬 203
10.1 滑坡滑帶摩擦生熱 203
10.2 微波輔助破巖 205
10.2.1 建立clump團(tuán)簇堆積模型 205
10.2.2 輝石和長(zhǎng)石的分組和材料設(shè)置 206
10.2.3 輝石受熱膨脹數(shù)值模擬 209
10.3 能源樁熱力耦合 211
10.3.1 建立能源樁-地層模型 212
10.3.2 熱力耦合過(guò)程數(shù)值模擬 214
10.4 地面沉降和地裂縫 217
10.4.1 切割地層模型 217
10.4.2 水位下降過(guò)程數(shù)值模擬 219
參考文獻(xiàn) 221
附錄 屬性、函數(shù)和常見(jiàn)問(wèn)題 224
附錄A 類的屬性 224
附錄B 主要函數(shù) 230
附錄C 常見(jiàn)問(wèn)題解答 243