利用3DEC仿真三维规则节理岩体的文章

个别元素法于三维岩体力学行为之应用
杨长义 陈志民 陈锦清 淡江大学土木工程学研究所 中兴顾问社大地力学研究中心　
摘 要本文利用三维个别元素分析法程序(3DEC)，仿真三维规则节理岩体在单轴与真三轴应力下之变形与强度特性。主要研究结果如下：(1) 3DEC可用以定性分析三维节理岩体力学机制，利用该程序可简易探讨任何节理分布位态下之力学行为，免除物理模模型仿真试验之困难；(2) 多轴应力下岩体之异向性行为亦可透过3DEC仿真分析，定性上均与物理现象相互一致；(3)在节理间距、劲度比较大的岩盘较需要比对二维与三维分析之差异。
一、前 言　
自然界岩体多处于真三轴应力状态下，以往受限于分析工具与实验设备，岩石工程之分析大多局限于二向度分析，对三维岩体行为之仿真则较少[1]。例如目前可用于分析具大变形特性之离散岩体的程序如DDA[2]、UDEC[3]均局限于对二维问题的解析；而3DEC[4]程序系以个别元素法(distinct element method)在UDEC基础下发展而成之数值分析程序，正可用以仿真三维节理岩体之力学行为：3DEC可将岩体视为由许多离散的完整岩块所组成，各完整岩块可以仿真成刚体或可变形体；而各完整岩块间为节理所分隔。(1)在对节理的仿真方面，主要依据位移－作用力法则，计算在节理面上之剪应力及正向应力，以作为个别岩块之边界应力条件，因此可仿真岩块大位移与转动之情况。(2) 3DEC在仿真可变形岩块时，系根据「edge」指令程序自行将三维岩块再细分成许多四面体状次级块体(sub-block)，次级块体可以为任意形状。每个次级块体可配合所指定之材料组成律及外力情况，计算岩块之受力及应力分布情况；每个次级块体的节点有三个自由度，以计算这些次级块体上节点之运动情形，然后配合材料组成律计算次级块体上之应力应变关系，则可得块体间之作用力，接着配合边界所产生的接触力计算得新合力与加速度，以作为下一时阶计算可变形岩块之边界条件。