udec 数值模拟(入门学)udec 数值模拟(入门学习)udec 数值模拟(入门学习)udec 数值模拟(入门学习).doc

UDEC 入门 ;new 是刷新udec窗口，从新调用一个程序 ;title 与heading代表标题，后面紧跟标题的名称。如： title hang dao mo ni ;round 指块体与块体之间的圆角半径，默认值是0.5，其值要求小于模型中最小块体的最短那条边长的二分之一。如： round 0.05 set ovtol=0.5 ;此命令是指层与层之间的嵌入厚度 block x1,y1 x2,y2 x3,y3 x4,y4 ;建立模型框架， crack x1,y1 x2,y2 ;两点划一线 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 ;在指定的区域生成一定宽度的单元（xu为x方向的取值） zone model mo range xl xu yl yu ;使指定的区域材料采用摩尔--库仑本构关系计算（即弹塑性） change jcons=2 range xl xu yl yu ;使指定的区域节理遵循摩尔--库仑准则计算（即弹塑性） change mat=1 range xl xu yl yu change mat=2 range xl xu yl yu change mat=3 range xl xu yl yu ;指定各岩层的材料标号 change jmat=1 range xl xu yl yu change jmat=2 range xl xu yl yu change jmat=3 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吨 zone k=0.15e9,g=0.1e9,fric=10.00,coh=0.19e6,ten=0.09e6 range mat=1 zone k=2.8e9,g=2.2e9,fric=30.00,coh=1.5e6,ten=0.4e6 range mat=2 zone k=6.9e9,g=6.6e9,fric=38.62,coh=5.63e6,ten=3.20e6 range mat=3 ;k为材料的法向刚度，g为材料的切向刚度，friction为材料的内摩擦角， ;cohesion为材料的内聚力，tension为材料的抗拉强度 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 ;采用位移法固定边界 solve\step 5000\cycle 5000 ;执行计算 save pingheng.sav ;保存文件， 文件的后缀为.sav，文件名可以自己命名。文件保存后，此刻的各种数据均在里面。 ;想要到某个已有的文件中找数据，可用restore，比如： restore pingheng.sav pause ;表示暂停，此时可以察看任何信息 continue ;继续调用下面程序段， ;udec自身仅提供8M or 16M内存空间，对于尺寸大，单元多的模型，需要额外增加内存，否则无法计算 ;改变内存方法有两种，其一，从Dos系统进入；其二，创建快捷图标，然后击右键属性 ;同时按住shif键与Z字母键，屏幕中会出现十字叉，摁住鼠标左键不放，移动光标直到你满意的窗口为止 ;让屏幕中出现十字叉，再