- 在下面这个最简单的模型中,切割块体涉及选择位置和指定裂缝,通过键入这些命令就可以产生10x10个单位的块体,然后劈裂成两个块体
- crack命令产生一个连续、水平的贯通模型的裂缝
- 圆角长度指定为0.1
round 0.1 ;round 指块体与块体之间的圆角半径,默认值是0.5,其值要求小于模型中最小块体的最短那条边长的二分之一。如:round 0.05 set ovtol=0.5 block(0,0) (0,10) (10,10) (10,0) Crack(0,5) (10,5) ;crack x1,y1 x2,y2 两点划一线 ;通过键入如下命令可以产生一个槽口: crack 2.5,10 5.0,7.5 crack 5.0,7.5 7.5,10 ;通过键入以下命令,就显示出包括块体地址号的块体图形 plot block num ;通过采用delete命令,能从模型中删除一个块体,比如删除槽口 delete range block 647
crack命令用于产生块体中的单一直线特诊的裂缝,裂缝由端点坐标(x1,y1)和(x2,y2)所确定 。 jset命令则是自动节理组生成器。根据所给定的特征参数(即倾角、迹长、岩桥长度、间距和空间位置)产生一组裂缝。
- 两组连续节理组的产生
- 第一个jset命令产生一组与x轴方向夹角为45°具有间距为2个单元的连续节理。
- 第二个jest命令产生与x方向夹角为-10°间距为1.5个单位的连续节理。
- 圆角长度的选择可能影响节理组的产状。
- 节理的位置
new round 0.01 block(0,0) (0,20) (20,20) (20,0) jset(45,0) (5,0.5) (0.5,0) (2,0) ;jset 90,0 4,0 4,0 6,0 0,-50 range jreg 3 jset 倾角,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n jset(-10,0)(5,0.5)(0.2,0) (1.5,0)
遇到了错误却找不对原因/(ㄒoㄒ)/~~ 一个基础案列就在这里卡壳了
一开始准备按照教程的基础案例一个一个敲来进行学习,但后来发现这样可能效率不高,掌握的慢,所以就想直接粘贴一个具体案例来搞明白每一块命令行是啥意思,但两三百行的代码从十几行开始就出错,别人在4.0上可以成功运行并且呈现理想的结果,我的7.0就不行T_T。
下面这图可能是先模拟出不同岩层的模型高度,执行很多次crack命令,接下来就开始出错。目前,找错误原因还是一大难题。
所以先总结一下用到的各种命令行是啥意思。
jregion id n x1,y1 x2,y2 x3,y3 x4,y4 delete
jset 90,0 4,0 4,0 6,0 0,-50 range jreg 3
;jset 倾角,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n
;其中xm,ym为起始点坐标,n为设置的区域标号
gen quad 10 range xl xu yl yu ; 在指定的区域生成一定宽度的单元 prop mat=1 dens=2000 prop mat=2 dens=2650 prop mat=3 dens=2700 ;指定各材料的密度,比如1号材料dens=2000,即1立方米重2吨 prop jmat=1 jkn=0.2e8,jks=0.1e7,jcoh=0,jfric=4,jten=0 prop jmat=2 jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0 prop jmat=3 jkn=20e8,jks=16e7,jcoh=0.4e6,jfric=15,jten=0 ;jkn为节理的法向刚度,jks为节理的切向刚度,jfriction为节理的内摩擦角, ;jcohesion为节理的内聚力,jtension为节理的抗拉强度 set gravity 0,-9.81 ;设置重力加速度,x方向为0,y方向为-9.8 bound xvel=0 range -0.1 3.00 -60.1 20.1 bound xvel=0 range 97 100.1 -60.1 20.1 bound yvel=0 range 0.1 100.1 -60.1 -58 ;采用位移法固定边界 set plin 1 x1,y1 x2,y2 n ;输出观测线的位移与应力数据,默认输出到udec.log文件中,可在excel表中处理数据
:在摸索了许久之后我似乎明白了UDEC是基于命令驱动格式这句话。