数值模拟udec软件学习第一周

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;在指定的区域生成一定宽度的单元

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;继续调用下面程序段

下面是一些小模型的代码：

Round 0.1；设置块体与块体之间的圆角半径为0.1
Block -10,-10 -10,10 10,10 10,-10；在（-10，-10）（-10，10）（10，10）（10，-10）范围里定义一个块体
tunnel 0,0 2,16
jset -70,0 40,0 0,0 40,0 -1,-1
jset -50,0 40,0 0,0 3,0 0, 2
gen edge 2.0
change jmat=2 range angle -51,-49
change jmat=5 range angle -71,-69
pro mat=1 d=2500 b=1.5e9 s=0.6e9
pro jmat=1 jkn=2e9 jks=2e9 jcoh=1e10 jten=1e10
pro jmat=2 jkn=2e9 jks=1e9 jfr=45
pro jmat=5 jkn=2e9 jks=1e9 jfr=10
change cons=0 range -1,1 -1,1

下图是运行出来的结果

round 0.1

block 0 0 0 100 150 100 150 0

crack 0 30 150 30 crack 0 40 150 40

crack 0 60 150 60 crack 70 30 70 33 crack 70 33 74 33 crack 74 33 74 30 crack 74 30 70 30 change cons=3 jregion id=1 0 0 0 30 150 30 150 0 jregion id=2 0 30 0 40 150 40 150 30 jregion id=3 0 40 0 60 150 60 150 40 jregion id=4 0 60 0 100 150 100 150 60 jset 0,0 150,0 0,0 6,0 range jreg 1 jset 90,0 6,0 6,0 15,0 (15 0) range jreg 1 jset 90,0 6,0 6,0 15,0 (7.5 6) range jreg 1 block being split jset 0,0 150,0 0,0 2,0 range jreg 2 jset 90,0 2,0 2,0 4,0 (4 30) range jreg 2 jset 90,0 2,0 2,0 4,0 (2 32) range jreg 2 jset 0,0 150,0 0,0 5,0 range jreg 3 jset 90,0 5,0 5,0 10,0 (10 40) range jreg 3 jset 90,0 5,0 5,0 10,0 ( 5 45) range jreg 3 jset 0,0 150,0 0,0 10,0 range jreg 4 jset 90,0 10,0 10,0 15,0 ( 15 60) range jreg 4 jset 90,0 10,0 10,0 15,0 (7.5 70) range jreg 4 change mat=1 range 0 150 0 30 change mat=2 range 0 150 30 40 change mat=3 range 0 150 40 60 change mat=4 range 0 150 60 100 generate edge 5 range 0 150 0 30 generate edge 2 range 0 150 30 40 generate edge 4 range 0 150 40 60 generate edge 10 range 0 150 60 100 delete range 70 74 30 33 bound stress -10 0 -10 insitu stress -10e6 0 -10e6 szz-10e6 ; rock blocks ; joints (no cohesion) prop m=1 dens=2500 k=5e9 g=3e9 prop jm=1 jkn=10e10 jks=10e10 jfric=30 prop m=2 dens=1000 k=5e8 g=3e8 prop jm=2 jkn=10e9 jks=10e9 jfric=20 prop m=3 dens=1500 k=1e9 g=1e9 prop jm=3 jkn=10e10 jks=10e10 jfric=25 prop m=4 dens=2500 k=5e9 g=3e9 prop jm=4 jkn=10e10 jks=10e10 jfric=32 change jcons=5 set jcondf=5 prop jmat jkn=40000 jks=40000 jfric=35 jrfric=10 & jcoh=2.0 jrescoh=0 & jdil=6 zdil=4e-4 bound xvel 0 range -0.01 0.01 -0.01 100.01 bound xvel 0 range 149.99 150.01 -0.01 100.01 bound yvel 0 range -0.01 150.01 -0.01 0.01 ;set gravity ; set in-situ stresses (total) grav 0 -9.8 hist type=2 hist unbal solve_ratio hist xdis=-40,100 ydis=-40,100 damp auto solve

数值模拟udec软件学习第一周相关推荐

迟到的小熊计算机画图画教案,课题(5、下)：《迟到得小熊—认识计算机画图软件》第一周.doc...
课题(5.下):<迟到得小熊-认识计算机画图软件>第一周上课时间:第一周星期本期第1.2课时 (详)案课型:新授课题:<迟到的小熊-认识计算机画图软件> 主编人:廖雄 ...
迟到的小熊计算机画图画教案,课题(5下)《迟到的小熊—认识计算机画图软件》第一周...
课题(5下)<迟到的小熊-认识计算机画图软件>第一周衡阳市船山实验小学衡阳市船山实验小学 (六)年级((六)年级( 下)期(信息技术)学科备课本电子格式下)期(信息技术)学科备课本电子格 ...
201671010140. 2016-2017-2 《Java程序设计》java学习第一周
java学习第一周本周是新学期的开端,也是新的学习进程的开端,第一次接触java这门课程,首先书本的厚度就给我一种无形的压力,这注定了,这门课程不会是轻松的,同时一种全新的学习方 ...
掉进悬崖的小白，捡到python基础秘籍，学习第一周——语法基础，小游戏，无脑循环
掉进悬崖的小白,捡到python基础秘籍,学习第一周--语法基础,小游戏,无脑循环人生苦短,我用python 语言的种类: 语言的发展: 什么是python 搭建 Python开发环境: 集成开发环 ...
Datacamp 学习第一周
Datacamp 学习第一周总结 Intro to Python for Data Science 列表 numpy package 字典 pandas comparision operators l ...
我的JAVA学习第一周知识总结，个人感受小结
JAVA学习第一周 1,了解三种开发模型 2,建包.建类的规范格式 3,了解不同数据类型 4,JAVA环境搭建 5,JVM.JRE.JDK是什么 6,if/else if语句的使用 7,各个语句swi ...
前端学习第一周-计算机基础，html标签的介绍，其中包括列表，表格以及表单的使用
前端学习第一周-HTML基础文章目录前端学习第一周-HTML基础前言 DAY 1 计算机组成计算机网络 DAY 2 vscode快捷方式 HTML骨架 DAY 3 HTML属性 HTML标签 ...
python学习第一周总结
python学习第一周总结一.Markdown语法 Markdown是一种轻量级标记语言,排版语法简洁,它使用易读易写的纯文本格式编写文档,可与HTML混编,可导出 HTML.PDF 以及本身的 . ...
迟到的小熊计算机画图画教案,课题(5、下)：《迟到及小熊—认识计算机画图软件》第一周.doc...
衡阳市船山实验小学 (六)年级( 下)期(信息技术)学科备课本电子格式上课时间:第一周星期本期第1.2课时 (详)案课型:新授课题:<迟到的小熊-认识计算机画图软件> 主编人:廖 ...

数值模拟udec软件学习第一周

数值模拟udec软件学习第一周相关推荐

最新文章

热门文章