(1) Elastic Joint Model
(2) Mohr-Coulomb Joint Model
(3) Softening-Healing Mohr-Coulomb Joint Model
(4) Bilinear Mohr-Coulomb Joint Model
(5) Power Law Creeping Joint Model
(6) Continuously Yielding Joint Model
(7) Nonlinear Joint Model
2.1 弹性节理模型
弹性节理模型(Elastic joint model)用来模拟构建几何模型时所需的虚拟节理 (construction joints),但这些节理并不对应实际的滑动节理或断层,其行为由节理的法向和切向刚度控制,不会发生屈服现象【弹性节理模型(Elastic Joint Model)与块体连接方法】。
2.2 莫尔-库仑节理模型
block contact jmodel assign mohr
block contact property stiffness-normal=100000 ...
100000 friction=30.0 ... =
dilation 15 dilation-zero 6e-4
block contact material-table default property ...
100000 stiffness-shear=100000 ... =
friction=30.0 dilation 15 dilation-zero 6e-4
(1) 当节理开始滑移时,在用户指定的临界滑移距离 Dc上,剪切强度会从峰值演变为残余值,这意味着并非所有滑移事件都是地震。
(2) 当滑动停止时,剪切强度会瞬间恢复到峰值。这意味着可能会出现节理粘滑行为。
剪切强度的演变有线性和非线性两种形式,对于非线性公式,一个指数(α>1)决定了强度下降的严重程度。
block contact jmodel assign softening-mohr
block contact property stiffness-normal 1e11 ...
stiffness-shear 1e11 friction 25 ...
friction-residual 20 tension 1e12 dc 1e-4
block contact jmodel assign bilinear-mohr
block contact property stiffness-normal 1e5 ...
stiffness-shear 1e5 tension 1e6
block contact property friction 30 friction-2 20 ...
cohesion 1 cohesion-2 2 ...
range position-x -0.1 0.1 position-y -0.1 0.1
2.5 幂律蠕变节理模型
block contact jmodel power
block contact property stiffness-normal [_jkn] ...
stiffness-shear [_jks] ...
cohesion [_jcoh] friction [_jfri] tension [_jten] ...
friction-residual [_jfri]...
constant-1 [_jpow_a1] exponent-1 [_jpow_n1]
2.6 连续屈服模型
(1) 剪应力/剪位移曲线总是趋向于节理的“目标”剪切强度,即曲线的瞬时梯度直接取决于强度与应力之间的差值。
(2) 目标剪切强度随着累积塑性位移 (损伤的量度)的增加而持续减小。
(3) 膨胀角被视为表观摩擦角 (由当前剪应力和法向应力确定)与残余摩擦角之间的差值。
基于这些假设,该模型自动展现出岩石节理常见的峰值/残余强度行为。此外,无论应变水平多么小,在所有应变水平的卸载和重新加载循环中都会表现出滞回特性【连续屈服节理模型(continuously yielding joint model)】。
block contact jmodel assign cyjm
block contact property kn-initial=100000 ...
100000 friction-residual=30.0 =
block contact property exponent-normal=0.0 ...
0.0 =
block contact property roughness=1.0e-4
block contact property friction-initial=59.3
block contact jmodel assign nonlinear
block contact property stiffness-normal=1e4 ...
1e4 friction=30.0 ... =
exponent-normal 1.1 kn-maximum 1e6
block contact material-table default jmodel nonlinear
block contact material-table default property ...
1e4 stiffness-shear=1e4 ... =
friction=30.0 exponent-normal 1.1 ...
kn-maximum 1e6 dilation 15