第1章 引言 1
1.1 拉格朗日法 1
1.2 欧拉法 3
1.3 混合方法 5
1.3.1 任意拉格朗日-欧拉法 5
1.3.2 质点网格法(PIC) 6
1.3.3 物质点法(MPM) 7
1.4 无网格法 13
第2章 控制方程 15
2.1 物体运动和变形的描述 15
2.2 变形梯度 17
2.3 变形率 19
2.4 柯西应力 20
2.5 焦曼应力率 21
2.6 更新拉格朗日格式 23
2.6.1 积分的物质导数 23
2.6.2 质量守恒 24
2.6.3 动量方程 24
2.6.4 能量方程 25
2.6.5 控制方程 26
2.7 冲击波 27
2.7.1 Hugoniot(雨贡纽)方程 27
2.7.2 人工体积粘性 29
2.8 爆轰波 31
2.8.1 CJ模型 31
2.8.2 ZND模型 34
2.9 温度计算 36
第3章 物质点法 38
3.1 更新拉格朗日格式的弱形式 38
3.2 物质点离散 39
3.3 运动方程的解法 42
3.3.1 显式求解 42
3.3.2 隐式求解 50
3.4 广义插值物质点法(GIMP) 52
3.5 接触算法 57
3.5.1 接触界面条件 58
3.5.2 接触判据 60
3.5.3 接触力 61
3.5.4 刚柔接触 63
3.5.5 算法实现 65
3.6 自适应 66
3.6.1 质点自适应 66
3.6.2 网格自适应 67
3.7 无反射边界 76
3.8 三维显式物质点法程序MPM3D 77
3.8.1 自由格式读入模块 78
3.8.2 基本数据的封装 80
3.8.3 程序流程 87
第4章 物质点有限元法 103
4.1 有限元法 103
4.1.1 更新拉格朗日格式的弱形式 104
4.1.2 有限元法离散 104
4.1.3 运动方程的解法 107
4.1.4 算法实现 115
4.1.5 与物质点法的比较 116
4.2 物质点有限元法 117
4.2.1 背景网格区域的求解 118
4.2.2 有限元网格区域的求解 120
4.2.3 时间积分 121
4.2.4 算法实现 121
4.2.5 算法应用 121
4.3 耦合物质点有限元法 121
4.3.1 接触探测 122
4.3.2 接触法线计算 123
4.3.3 接触力 124
4.3.4 时间积分 125
4.3.5 算法实现 125
4.3.6 算法应用 126
4.4 自适应物质点有限元法 127
4.4.1 转化算法 128
4.4.2 耦合算法 129
4.4.3 算法实现 131
4.4.4 算法应用 131
4.5 杂交物质点有限元法 131
第5章 材料模型 135
5.1 应力更新 135
5.2 本构模型 137
5.2.1 弹性模型 137
5.2.2 单塑性模型 138
5.2.3 返回映射法 141
5.2.4 J2流动理论 147
5.2.5 Johnson-Cook模型 152
5.2.6 Drucker-Prager模型 154
5.2.7 Gurson模型 161
5.2.8 Deshpande-Fleck模型 164
5.2.9 HJC混凝土材料模型 168
5.2.10 RHT混凝土材料模型 172
5.2.11 JH-2陶瓷材料模型 177
5.2.12 牛顿流体材料模型 180
5.2.13 炸药材料模型 181
5.3 状态方程 182
5.3.1 多方过程 183
5.3.2 不可压缩流体 183
5.3.3 线性多项式状态方程 184
5.3.4 JWL状态方程 185
5.3.5 Mie-Grüneisen状态方程 186
5.3.6 P-α状态方程 188
5.4 失效模型 189
5.5 材料模型程序实现 190
5.5.1 MaterialData模块 191
5.5.2 MaterialModel模块 192
第6章 物质点法的应用 216
6.1 动态断裂问题 216
6.2 碰撞问题 219
6.3 爆炸问题 225
6.4 流固耦合问题 230
6.5 多尺度模拟 234
6.6 生物力学问题 239
6.7 其他超大变形问题 242
附录A MPM3D-F90程序 245
A.1 程序使用说明 245
A.1.1 命令行执行方式 245
A.1.2 文件类型 245
A.2 输入文件 246
A.2.1 单位 246
A.2.2 关键字 246
A.2.3 输入文件范例 251
A.3 典型算例 253
A.3.1 一维板条TNT爆轰 253
A.3.2 泰勒杆碰撞 253
A.3.3 卵形弹斜侵彻 255
A.3.4 砂土边坡失效 256
索引 259
参考文献 262