UDEC 1

   

  • 在下面这个最简单的模型中,切割块体涉及选择位置和指定裂缝,通过键入这些命令就可以产生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是基于命令驱动格式这句话。

原文地址:https://www.cnblogs.com/raindawn/p/11923333.html