CDEM的计算原理及计算方法-YongliangWang.PDF

目 录 第一章 力学所CDEM的最新进展5 1.1 离散元方法在地质工程中的应用前景 5 1.2 山体稳定性分析面临的问题 7 1.2.1 滑坡体的复杂破坏现象 8 1.2.2 实验表明滑坡体复杂的破坏现象 10 1.2.3 离散元方法的一些不足 12 1.3 连续介质力学模型离散元方法的基本构成 13 1.3.1 CDEM 可变形块体离散元的计算方法 14 1.3.2 可变形块体离散元的计算步骤 18 1.4 连续介质力学模型离散元方法的基本特点 18 1.5 发展连续介质力学模型离散元方法的目的和意义 25 第二章 离散元方法的基本几何知识27 2.1 平面方程 27 2.2 点与平面的相对位置 29 2.3 线段与平面的交点 29 2.4 块体与平面相交 30 2.5 描述岩块体的基本要素 32 第三章 有限元方法简介33 3.1 离散化 33 3.2 单元的位移模式和形函数 34 3.3 应变 38 3.4 应力 41 3.5 单元势能和单元刚度矩阵 42 第四章 弹簧振子与刚体弹簧系统的振动规律48 4.1 弹簧振子的基本方程 48 4.4.1 单自由度系统 48 4.4.2 多自由度系统 49 4.2 有阻尼的弹簧振动方程 49 4.2.1 单自由度系统 49 4.2.2 多自由度系统 50 4.3 弹簧振动的基本物理量 50 4.4 弹簧振动的基本规律 50 4.4.1 单自由度系统无阻尼自由振动 50 4.4.2 单自由度系统有阻尼的自由振动 51 4.4.3 多自由度系统无阻尼振动 54 4.4.4 多自由度系统经典阻尼振动 56 4.5 无量纲量的基本表达方式 57 第五章 刚性块体离散元59 5.1 刚性块体离散元方法的基本方程 59 5.1.1 运动方程 59 5.1.2 物理方程 60 5.2 块体离散元的参数选取 62 5.2.1 Cundall 模型的取法 62 5.2.2 面-面接触模型中等效接触刚度选取 62 5.2.3 等效刚度的尺寸效应 65 5.2.4 关于阻尼系数的物理意义及其选取方法 66 5.3 块体转动惯量的计算方法 68 5.4 计算流程图 71 5.5 计算时步 71 5.6 破坏准则 73 第六章 基于连续介质力学模型的离散元74 6.1 CDEM 中的几何问题 74 6.2 CDEM 中单元的生成 75 6.2.1 平行六面体单元 75 6.2.2 多尺度平行六面体单元 76 6.2.3 随机平行六面体单元 77 6.2.4 凸多面体单元 78 6.3 质点控制方程 79 6.3.1 弹簧破坏前质点平衡方程 79 6.3.2 块体平衡方程 79 6.3.3 边界条件 79 6.4 质点力的计算 80 6.4.1 结构层的变形力 80 6.4.2 单元变形力 81 6.5 数值破坏模型 82 2 6.5.1 六面体块体结构面破坏的计算模型 82 6.5.2 随机节点后六面体块体结构面破坏的计算模型 84 6.5.3 单元界面断裂计算模型 85 6.5.4 块体内部破坏模型 88 第七章 多尺度块体离散元94 7.1 两尺度模型 95 7.1.1 划分块体的几何描述方法 95 7.1.2 划分块体的原则 95 7.1.3 块体分裂的基本步骤 96 7.1.4 各级块体的几何描述 96 7.1.5 块体接触面的描述 100 7.1.6 不同尺度块体之间力和位移的传递方式 104 7.1.7 两尺度多级划分块体计算方法流程图 109 7.1.8 局部网格细化的实现 110 7.2 多尺度模型 112 7.2.1 块体划分的原则 112