UDEC模拟
udec命令
1.把图形保存下来，能在AUTOCAD中打开、编辑 plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形 set plot 256 dxf；设置图形为256色 set out c:\ss.dxf copy c:\ss.dxf；把这图形以ss.dxf文件保存在C磁盘下 2.把数据导出 set log on print pline 1 ydis;把pline 1的y方向位移的数据导出到UDEC目录下udec.txt文件中 set log off 3.plot block stress ;显示块体的应力 plot block dis ;显示块体的位移(有x和y方向) plot block pl ;显示块体的塑性区(plastic) save xx.save ;保存计算结果 res xx.save ;调用 4.液压支架的命令 supp xc yc wid l seg n mat j;(xc,yc)是指中心点坐标,l指支架的宽度,n指分段数,mat j指支柱材料性质为j prop mat j sup_kn -1 ;sup_kn指支架的刚度,-1与表1相对应 del range x1 x2 y1 y2 ;挖掘范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体,由支架支撑 table 1 0 42.e6 0.05 5.0e6 0.1 6.0e6;表1表示的是液压支架的(P-DS)特性曲线 5 巷道施工中断面加梯子梁的命令(见图1) stru gen xc yc np 100 fa a thetra b mat=16 thick=0.2 ;100 指分100段,a b 指角度 prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800 prop mat=16 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7 prop mat=16 if_kn=1.35e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;interface-界面 例子讲解: 建模(以米为单位,;后为解释部分) round=0.1 ;方块的圆角块半径为0.1米 set ovtol=1.0 ;块体与块体之间相互嵌入量最大值为1米 bl 0,0 0,26 50,26 50,0 ;在(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)生成块体范围 crack 0,25 50,25 ;在(x1, y1)(x2 ,y2)两点间画直线 jregion id 1 0,0 0,10 50,10 50,0 delete ;删除此区域的块体 jset 90,0 2.5,0 2.5,0 3,0 0,0 range jregion 1 ;jset A,0 a,0 b,0 c,0 x0,y0 range jregion 1 jset 90,0 2.5,0 2.5,0 3,0 1.5,2.5 range jregion 1 可画成列的线段(见图2) jset 0,0 26,0 0,0 2.5,0 0,0 range jregion 1 pa ;当程度运行至此时暂停,可以看看你所建的部分模型,用continue继续运行下面部分 save t.save ;建模保存在t.save中 岩层赋属性及原岩力平衡计算 res t.save ;调用已建好的模型 gen quad 12 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体 zone model mo range 0,180 0,45.2 ;第一个x，第二个y在范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体符合库仑准则 change jcons=2 range 0,180 0,45.2 ;节理面间接触-coulmb滑移 ;下面主要是讲岩层赋属性
UDEC学习整理资料
1、角点必须按顺时针方向排列； 2、Crack 命令用于产生块体中单一直线特征的裂缝。裂缝由端点坐标（x1，y1）和（x2，y2）所确定。 3、Jset 命令则是自动节理组生成器。根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。 4、round d---d是圆角距离，建议在block命令前指定圆角长度。 5、DELETE 命令，能从模型中删除一个块体。例如，为了删除槽口块体，delete range 4.5，5.5 8，10。 6、GEN命令激活三角形网格有限单元自动生成器。命令GEN edge v 将作用于任意形状的块体。其v值定义三角形单元的最大边长，即v值越小，块体中的单元越小。应当注意的是：具有高的边长比值的块体并不能产生单元，其极限的比重近似为1:10。 7、采用命令GEN quad v，指定模型为塑性材料模型的单元。该类型的单元提供了对于塑性问题的精确解。然而，GEN quad 命令可能对某些形状的块体不起作用。在此情况下，应当采用GEN edge 8、Change 命令改变块体为指定的变形块体。Cons=0意味着模型块体材料被移出或开挖。Cons=1 改变块体为各向同性弹性特性；而Cons=3则改变块体为摩尔－库仑模型，考虑塑性特性。缺省值为所有变形体则自动改变为Cons=1。P21 9、cha nge jcons=2，所以不连续结构面的缺省模型是Jcons=2。 10、可用以下命令检查材料号Plot block mat
12、INSITU命令用来初始化应力。采用该命令，可以赋值初始应力。 13、hist xvel 5, 5 hist ydisp 0, 11 第一个是记录位移坐标（x=5，y=5）附近结点x方向的速度，而第二个是记录接近坐标（x=0，y=11）位置处y方向的位移。 14、set grav 0.0 , -9.81第一个是x方向的加速度，第二个值为y方向的加速度为9.81m/sec2（向下作用）。 15、UDEC计算可在任何时间通过按
FLAC3D常见命令与使用技巧
FLAC3D常见命令与使用技巧 1、FLAC3D常见命令： 是有限元程序吗答：不是！是有限差分法。 2.最先需要掌握的命令有哪些 答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 3.怎样看模型的样子答：plo blo gro可以看到不同的group的颜色分布 4.怎样看模型的边界情况答：plo gpfix red 5.怎样看模型的体力分布答：plo fap red 6.怎样看模型的云图答：位移：plo con dis (xdis, ydis, zdis)应力：plo con sz (sy, sx,sxy, syz, sxz) 7.怎样看模型的矢量图答：plo dis (xdis, ydis, zdis) 8.怎样看模型有多少单元、节点答：pri info 9.怎样输出模型的后处理图 答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件 10.怎样调用一个文件答：File/call或者call命令 10.如何施加面力答：app nstress 11.如何调整视图的大小、角度答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键 12.如何进行边界约束答：fix x ran（约束的是速度，在初始情况下约束等效于位移约束） 13.如何知道每个单元的ID答：用鼠标双击单元的表面，可以知道单元的ID和坐标 14.如何进行切片 答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 15.如何保存计算结果答：save +文件名. 16.如何调用已保存的结果答：rest +文件名；或者File / Restore 17.如何暂停计算答：Esc 18.如何在程序中进行暂停，并可恢复计算答：在命令中加入pause命令，用continue进行继续 19.如何跳过某个计算步答：在计算中按空格键跳过本次计算，自动进入下一步 20. Fish是什么东西 答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能
UDEC经典学习总结 张科学修改 完成
1.把图形保存下来，能在AUTOCAD(图片)中打开、编辑 plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形 set plot dxf(jpg) 256；设置图形为256色 set out c:\ss.dxf(jpg) (set out 1.dxf 这个文件保存在当前带数值模拟的文件夹里) copy c:\ss.dxf(jpg) ( copy 1.dxf 这个文件保存在当前带数值模拟的文件夹里) ；把这图形以ss.dxf文件保存在C磁盘下 2.把数据导出 set log on print pline 1 ydis;把pline 1的y方向位移的数据导出到UDEC目录下udec.txt文件中 set log off 3.plot block stress ;显示块体的应力 plot block dis ;显示块体的位移(有x和y方向) plot block pl ;显示块体的塑性区(plastic) save xx.save ;保存计算结果 res xx.save ;调用 4.液压支架的命令 supp xc yc wid l seg n mat j;(xc,yc)是指中心点坐标,l指支架的宽度,n指分段数,mat j指支柱材料性质为j prop mat j sup_kn -1 ;sup_kn指支架的刚度,-1与表1相对应 del range x1 x2 y1 y2 ;挖掘范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体,由支架支撑 table 1 0 42.e6 0.05 5.0e6 0.1 6.0e6;表1表示的是液压支架的(P-DS)特性曲线 5 巷道施工中断面加梯子梁的命令(见图1) stru gen xc yc np 100 fa a thetra b mat=16 thick=0.2 ;100 指分100段,a b 指角度 prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800 prop mat=16 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7 prop mat=16 if_kn=1.35e9 if_ks=1.35e9 if_tens=0 if_fric=18 if_coh=0;interface-界面 例子讲解: 建模(以米为单位,;后为解释部分) round=0.1 ;方块的圆角块半径为0.1米 set ovtol=1.0 ;块体与块体之间相互嵌入量最大值为1米 bl 0,0 0,26 50,26 50,0 ;在(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)生成块体范围 crack 0,25 50,25 ;在(x1, y1)(x2 ,y2)两点间画直线 jregion id 1 0,0 0,10 50,10 50,0 delete ;删除此区域的块体 jset 90,0 2.5,0 2.5,0 3,0 0,0 range jregion 1 ;jset A,0 a,0 b,0 c,0 x0,y0 range jregion 1 jset 90,0 2.5,0 2.5,0 3,0 1.5,2.5 range jregion 1 可画成列的线段(见图2) jset 0,0 26,0 0,0 2.5,0 0,0 range jregion 1 pa ;当程度运行至此时暂停,可以看看你所建的部分模型,用continue继续运行下面部分 save t.save ;建模保存在t.save中 岩层赋属性及原岩力平衡计算 res t.save ;调用已建好的模型 gen quad 12 ;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体 zone model mo range 0,180 0,45.2 ;在范围(x1, y1)(x1 ,y2)(x2, y2)(x2, y1)的块体符合库仑准则 change jcons=2 range 0,180 0,45.2 ;节理面间接触-coulmb滑移
udec命令总结精华-正宗
《Udec 命令总结》精华版 1. 安装、打开、保存 (1) 安装：①执行Udec 3.1→将Crack文件中的内容替换； ②复制Udec.exe,粘贴为快捷方式→属性，目标，加入空格256【开始内存8M，将内存 改为256】，可以复制快捷方式至桌面或硬盘。 (2) 打开：①直接在udec> 命令行输入； ②写好程序，udec>命令输入call，然后将*.txt文件拖入命令行，执行 (3) 保存：输入save d:\kaicai.sav, 调用命令rest，将kaicai.sav 拖入命令行 【！！！文件名最好不用汉字，有时候不识别】 注释：如果保存为save d:\111\kaicai.sav, 注意其中111文件必须提前建好，否则无法保存或者保存错误 (4) 操作：Ctrl+Z选中图像可以放大，Ctrl+Z 双击复原，屏幕中会出现十字叉，按住鼠标左键不放， 移动光标直到你满意的窗口为止；pause暂停，此时可以察看任何信息；continue继续调 用下面程序段。ESC可以随时进行停止，但不能继续；英文分号; 表示注释不运行命令。 2. 基本命令 2.1 基本设置 Udec>n【new刷新窗口，从新调用一个程序，修改后的*.txt文件必须输入n，重新运行文本文件】 Udec>title 【或heading代表标题，后面紧跟标题的名称。如：hang dao mo ni 或济宁三号井围岩变形破坏规律研究】 Udec> round d 【“圆角”命令，Udec中所有的块体都有圆角，目的是为防止块体悬挂在有棱角的节点上，由于块体悬挂将产生应力集中。d指块体与块体之间的圆角半径，默认值是0.5，其值要求小于模型中最小块体的最短那条边长的二分之一，最大圆角长度不能超过块体平均棱长的1%。在block 命令前指定圆角长度。如：round 0.05】 圆角 图1 10×10块体圆角 Udec>set ovtol 0.5 【此命令是指层与层之间的嵌入厚度，当提示为“overlap too large”时就需要修改此值更大一些，可以显示设计的块体，plot overlap！嵌入太大的原因可能为块体强度太小】
UDEC4.0使用说明
菜单驱动模式运行离散元 1、菜单驱动模式运行离散元 对于Itasca加码图形界面是一个菜单驱动的图形界面开发，以协助助用户掌握Itasca代码。在UDEC中，UDEC—GIIC很容易与点和点击式操作，以访问所有的命令和设施。该GIIC结构是专门用来模拟预期的Windows功能，并允许显示的项目相对应的离散元操作的一般性鼠标性操作。你可以能够立即使用UDEC解决问题，无须通过命令来选择你需要的分析。本节提供了一个GIIC的介绍，并包括一个简单的教程，以帮助您开始。你会注意到在GIIC主菜单栏中一个帮助菜单。帮助按钮还包括在GIIC中的每一个工具，并且帮助窗格可以通过在模型工具标签上右击打开。咨询帮助意见可以得到具体的GIIC功能的详细信息。 图1-1 UDEC—GIIC主窗口 在利用UDEC进行全面的项目分析之前，我们强烈建议你阅读离散元用户指南和核查问题和示例应用程序，从而对离散元模型、分配材料特性、模型的初始条件和计算程序获得一个一般性的认识和理解。 1.1进入GIIC并选择分析选项 在开始/程序/Itasca/udec菜单中，当UDEC加载时，你可以选择“UDEC 4.01 with GIIC”，自动启动GIIC。或者，您可以在打开UDEC时，如果你在文本模式下，你可以在“udec>”提示下键入“giic”命令。GIIC主窗口如图1.1所示。 该代码名称和当前版本号印在标题栏中的窗口顶部，主菜单栏位于标题栏下方的位置。在主菜单栏下方有两个窗口：一个资源窗口和一个模型视图窗口。资源窗口包括四个以文本为基本信息的标签。“console”（控制台）标签显示文本输出和允许命令行输入（在窗口的底部）。“record”窗口显示生成当前模型项目状态的命令的记录。该记录以“项目树”的形式，显示保存文件之间的变化。保存状态显示在树状结构。数据可以作为离散元组命令形式导出到数据文件，命令组代表所分析的问题。“FISH”窗口可以打开FISH编辑器，能方便执行FISH的功能。项目说明显示在“Notes”窗口。 “Model-vies pane”（模型视图窗口）显示模型的图形化视图，附加标签查看可以添加到这个窗口，它显示用户自定义的图表。在模型视图窗口上部是包含建模阶段标签的标签栏。当你点击模型阶段标签的时候，就会打开一个工具栏，这包括访问模型工具的按钮。在图1.1中，生成的工具栏显示在模型选项对话框后面，当在您单击工具栏按钮，这将打开一个建模阶段的
UDEC模拟实例与解析
UDEC 实例翻译与命令解析中铁隧道集团科研所——珠穆浪玛UDEC 实例翻译与命令解析 翻译：珠穆朗玛
1 地震诱发地层坍塌 Seismic-Induced Groundfall 1.1 问题描述 本例展示使用 UDEC 模拟分析地震诱发地层坍塌的一类的问题,模型见图 1.1,该模型基 于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿 34-1-554 切割断面的一个剖面图的结构和 尺寸. 用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高 5m,宽 10m. 假定两个连续节理交叉平面分析:一个角度为 45 度,另一个为-9 度,两者节理间距均为 5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳 定性。围岩参数来自试验室平均测试数值，假定岩石块体参数如下： 假定块体仅具有弹性行为，节理假定符合库伦滑动准则，选择典型的教课书数值作为节 理参数，如下：
初始应力状态按各向同性估计为24Mpa（假定垂直荷载由覆盖深度大约800m 的岩层产生）。 1.2 UDEC 分析 UDEC 模拟顺序分三个阶段,首先,模型在初始应力状态下进行无超采固结.其次,进行开挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图1.2.超采正上方和下方的块体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题域的外周边界引入粘滞边界用以消除波的反射.从而模拟有限的岩体,地震事件用施加到模型顶部y 方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中,施加1.25Mpa 的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部, 施加的有效影响应力应该是1.25 MPa/2, or 0.625 MPa.0.02 秒后的开挖面拱顶的应力分布见图1.3,两点的位移被监测,1 点位于开挖面的左角,点2 位于拱顶块体的右角, 图1.4 的位移时间曲线显示两点本质上是弹性反应 .
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 ;在指定的区域生成一定宽度的单元 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
UDEC学习整理资料
t i m e a n d A l 1、角点必须按顺时针方向排列；2、Crack 命令用于产生块体中单一直线特征的裂缝。裂缝由端点坐标 （x1，y1）和（x2，y2）所确定。 3、Jset 命令则是自动节理组生成器。根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。 4、round d---d 是圆角距离，建议在block 命令前指定圆角长度。 5、DELETE 命令，能从模型中删除一个块体。例如，为了删除槽口块体， delete range 4.5，5.5 8，10。 6、GEN 命令激活三角形网格有限单元自动生成器。命令GEN edge v 将作用 于任意形状的块体。其v 值定义三角形单元的最大边长，即v 值越小，块体中的单元越小。应当注意的是：具有高的边长比值的块体并不能产生单元，其极 限的比重近似为1:10。 7、采用命令GEN quad v ，指定模型为塑性材料模型的单元。该类型的单元提 供了对于塑性问题的精确解。然而，GEN quad 命令可能对某些形状的块体不 起作用。在此情况下，应当采用GEN edge 8、Change 命令改变块体为指定的变形块体。Cons=0意味着模型块体材料被 移出或开挖。Cons=1 改变块体为各向同性弹性特性；而Cons=3则改变块体为 摩尔－库仑模型，考虑塑性特性。缺省值为所有变形体则自动改变为 Cons=1。P21 9、cha nge jcons=2，所以不连续结构面的缺省模型是Jcons=2。10、可用以下命令检查材料号Plot block mat 11、边界条件命令总结 命 令 效 果 Stress 施加总应力到刚体或变形体块体的边界上Xload 施加刚体或变形体边界的x 方向的荷载Yload 施加刚体或变形体边界的y 方向的荷载Xvel 施加变形体边界的x 方向的速度（位移）BOUNDARY Yvel 施加变形体边界的y 方向的速度（位移） FIX 固定刚体边界的速度（位移）FREE 释放刚体的速度（位移） LOAD Xload 施加x 方向的荷载到刚体的边界
UDEC滑坡实例步骤
1、加载UDEC进入DOC环境后输入giic或者gui命令，然后进入主菜单 2、Model option 选择合适条件通常情况下，你可以使用默认域联系（domain-logic）检 测模式。如果你想监测任何块体的位移，这些块体可能从隧道顶部分离或掉落，你应该使用“cell-space detection”模式跟踪位移和下落块体的潜在接触。 3、命名并且保存文件 4、New block 建模，根据需要设置模型的长30 宽15 415 410 405 400 395 390 385 380 375 370 415 410 405 400 395 390 385 380 375 370 420 425 430 435 440 445 450 455 m 5、Bound 调节边界，与实际相符 6、Crack 添加节理，（层状岩体是否按节理处理？）岩层20°∠34°，J1产状60°∠15°J2产状为35°∠47°，J3产状为95°∠89°（怎么将不同产状节理进行转换？） 路线设计好，为后来开挖做好准备。 7、execute 执行文件 8、zone 执行长度为0.5的最大区域边界，划分网格 9、Zone material 创建一个或者几个块体材料属性，选择一种本构模型，本次选择的是 Mohr-Coulomb模型 prop mat 1 den--2143 bu=30e9 sh=18e9 c=1.2e5 f=21 t--2e5 prop mat 2 den=2260 bu=40e9 sh=24e9 c=1.5e5 f=28 t--2.5e5 prop mat 3 den--2300 bu=50e9 sh=28e9 c=3.5e5 f=32 t--3.5e5
UDEC60常用命令集锦
rest fall1.sav delete range -2,2 -2,2 reset disp ;位移归0 reset hist hist unbal hist ydis 0,2 step 2000 plot hold block stress disp edge Joints jregion
Systems of Units Hist hist solve type 1 solve rat 1e-5 save tun1.sav hist ydis 42 -67 sxx 42 -67 hist ydis 42 -73 sxx 42 -73 reset disp jdisp solve rat 1e-5 save tun2.sav
Plot pl syy fill projected pl syy fill bfill pl ydisp fill pl bl id black plot pen stress block;显示应力 pen作用是udec4.0窗口显示pl se ;看裂隙 pl
;抓取图像，在autocad内处理set pl dxf 256 copy xxx.dxf
Zone joint 力学属性 zone model elastic density 2E3 bulk 8E9 shear 5E9 range group ’block’group joint ’joint’ joint model area jks 5E10 jkn 2.5E10 jfriction 30 r ange group ’joint’ ; new contact default set jcondf joint model area jks=5E10 jkn=2.5E10 jfriction=30 insitu stress -2500000.0,0.0,-1.0E7
UDEC实例详解[1]
实例1 荷载作用下边坡稳定问题 朔准线黄河大桥，边坡高105m，层理间距取3m，岩层视倾角6°；主节理间距取2m，视倾角85°，计算模型如下： DEM计算模型 相关计算参数取值如下：桥基长15m，宽12m，荷载180000kN；灰岩弹性模量E=4×104MPa，泊松比μ=0.23，岩体密度γ=2.7×103kg/m3。 new round 0.05 set delc off block 0 0 0 50 50 50 62.6 115.2 68.9 123.9 85.5 132.5 & 105.5 140.1 125.4 156.4 180 156.4 180 0 jset 85,0 300,0 0 0 3,0 (50 50) jset 6 0 300 0 0 0 4 0 (50 50) gen auto 5 change jmat=1 range angle 84 86 change jmat=2 range angle 5 7 prop mat=1 den=2700 b=1.5e9 s=0.6e9 prop jmat=1 jkn=2e9 jks=2e9 jcoh=1e9 jfr=25 jten=1e5 节理参数不同prop jmat=2 jkn=2e9 jks=2e9 jcoh=1e10 jfr=35 grav 0 -10 hist n=100 xdis 68.9 123.9 ydis 68.9 123.9 hist n=100 xdis 125.4 156.4 ydis 125.4 156.4 hist unbal bound xvel=0 range -1 0.2 -1 50.1
UDEC 数值模拟(入门学习)
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 ;指定各岩层的节理标号
UDEC模拟实例与解析
^` UDEC 实例翻译与命令解析 翻译：珠穆朗玛
1 地震诱发地层坍塌 Seismic-Induced Groundfall 1.1 问题描述 本例展示使用 UDEC 模拟分析地震诱发地层坍塌的一类的问题,模型见图 1.1,该模型基 于加拿大安大略省萨德伯里市鹰桥公司弗雷则矿 34-1-554 切割断面的一个剖面图的结构和 尺寸. 用二维平面应变模型代表垂直于超采轴向方向的平面效应,超采面高 5m,宽 10m. 假定两个连续节理交叉平面分析:一个角度为 45 度,另一个为-9 度,两者节理间距均为 5m,为了演示的目的,一个近似垂直的“虚拟节理”也被添加到块体内开挖面顶部以增强不稳 定性。围岩参数来自试验室平均测试数值，假定岩石块体参数如下： 假定块体仅具有弹性行为，节理假定符合库伦滑动准则，选择典型的教课书数值作为节 理参数，如下：
初始应力状态按各向同性估计为24Mpa（假定垂直荷载由覆盖深度大约800m 的岩层产生）。 1.2 UDEC 分析 UDEC 模拟顺序分三个阶段,首先,模型在初始应力状态下进行无超采固结.其次,进行开挖并且模型循环至平衡状态.本阶段超采面周围的应力分布见图1.2.超采正上方和下方的块体滑动后稳定.在第三阶段.估计了两个不同的峰值速度的地震事件.对所有地震模拟,在问题域的外周边界引入粘滞边界用以消除波的反射.从而模拟有限的岩体,地震事件用施加到模型顶部y 方向的正弦应力波表现.应力波被叠加到已存在的初始地应力上.在第一个模拟中,施加1.25Mpa 的峰值应力,应当注意的是,由于粘滞边界条件实际是在模型顶部, 施加的有效影响应力应该是1.25 MPa/2, or 0.625 MPa.0.02 秒后的开挖面拱顶的应力分布见图1.3,两点的位移被监测,1 点位于开挖面的左角,点2 位于拱顶块体的右角, 图1.4 的位移时间曲线显示两点本质上是弹性反应 .
(完整word版)Udec计算算例
1.3 UDEC算例 1.3.1工程概况 某隧道位于一包含高角度连续节理岩体内，节理倾角为50度，平均间距为7m，隧道为一半径为9m的圆形隧道。贯穿于开挖面内的一垂直断层，在隧道拱顶形成了一个三角楔形体。 本算例使用UDEC的结构单元逻辑来模拟喷射混凝土和锚杆联合支护的圆形隧道开挖问题。 1.3.2构建模型 隧道埋深451m，为半径9m圆形隧道，本次计算模型左右边界取41m，隧道至上下边界也取41m。总的来说，模型长100m，宽100m。计算模型如图1-1所示。 图1-1 UDEC计算模型 1.3.3计算参数 在包含高倾角节理和垂直断层的岩体内进行圆形开挖的UDEC模型岩体、节理和断层参数如下所示： 表1-1 完整岩石物理力学参数 密度Dens （kg/m3）体积模量K （Gpa） 剪切模量G （Gpa） 2500 1.5 0.6
表1-2 节理、断层物理力学参数 表1-3 喷射混凝土物理力学参数 表1-4 岩体和喷射混凝土接触面物理力学参数 表1-5 锚杆物理力学参数 作为演示的目的，隧道开挖和支护是瞬时发生的。本算例共两种支护分析被计算：第一，只施加喷射混凝土衬砌；第二，喷射混凝土和锚杆联合提供支护。 为了在第二种支护情况分析中可以更清晰的看到锚杆提供的支护，算例采用喷射混凝土的抗压强度被设置成一个很低的值，且厚度仅取为10cm。 1.3.4模拟步骤 1.建立模型 在UDEC中输入以下命令可建立隧道结构模型及边界。如图1-2所示。 new round 0.1 block -50,-50 -50,50 50,50 50,-50 jset -50,0 100,0 0,0 7,0 ；刷新UDEC窗口，重新调用一个新程序 ；块与块之间的圆角半径，必须小于块体最小边的1/2 ；建立模型框架 ；设置节理
UDEC指令 基本命令
;new titlemoni;定义名称moni 块体都有“圆角”，其目的在于避免块体悬挂在有棱角的节点上。由于块体悬挂引起应力集中。然而，圆角值存在与模型有关的上限值。对于变形块体，最大圆角长度应当不超过块体平均棱长的1％。圆角长度可以如下命令加以改变：roundd round=0.025；d是圆角距离（缺省值是d=0.5）。模型中的所有圆角长度都是相同的。 Forexample,ifthecommands SETedge0.4andROUND0.1arespecified,thenblockedgelengthss mallerthan0.4willnotbecreated,andtheroundinglengthforblockswillbe0.1.Thesecommandsmustbegi venbeforetheBLOCKcommand. setovtol=1;块体与块体之间相互嵌入量最大值为1米 bl0,00,100300,100300,0;定义范围，四点坐标，顺时针方向 cr0,15300,15;Crack命令用于产生块体中单一直线特征的裂缝，裂缝由端点坐标（x1，y1）和（x2，y2）所确定。 jregid10,00,15300,15300,0delete;定义节理，命名1，定义范围，四点坐标，delete为常规语言 jset90,05,05,05,00,0rangjreg1;Jset命令则是自动节理组生成器。根据所给定的特征参数（即倾角、迹长、岩桥长度、间距和空间位置）产生一组裂缝。定义节理角度90°；节理的长度；节理的距离，即纵向间隔；横向距离，即横向隔5一个；起始点坐标为0,0 jset90,05,05,05,02.5,5rangjreg1;定义节理的另外一项 genquad106range03000100;定义块体最大变形,若没有此语句,刚所有块体均为刚性块体；automaticgenerationofdiagonallyopposedtriangularzonestoimproveplasticflowcalculation.(对角三角形区域的自动生成改善塑性流动计 算)Parametersxwandywarezonewidthsinthex-andy-directions zonemodelmorange03000100;定义摩尔库伦模型的范围，X的范围，Y的范围 changejcons=5range03000100;定义节理的某个属性 changemat=1rangereg0,00,15300,15300,0;定义物质1的范围，即赋予这个范围为物质1,只赋名，没有定义属性 changemat=2rangereg0,150,20300,20300,15 changejmat=1rangereg0,00,15300,15300,0;定义节理1的范围 setjmatdf2;定义节理属性，使符合摩尔库伦准则 savech.sav;保存文件 calprop-25.txt;调出文件 resch.sav;调出前面保存的文件接着计算 ;参数设置
udec注意事项
3.9 提示和建议 当用UDEC求解问题时，重要的事情是为有效地进行问题求解而进行模型优化设计。本节为改进运行模型提出一些建议。而且，还给出了进行UDEC计算时应避免可能出现的一些陷井。 3.9.1 节理几何形状的选择 对于选择输入到UDEC模型中的节理几何形状是分析中的关键一步。对具有10至数百条节理的多裂隙岩体，考虑所有的节理结构是不现实的。通常，为实际分析产生具有合理大小的模型和运行速度，仅很小比例的节理能够被输入模型中。因此，为模拟力学响应，模拟者必须对节理的几何数据进行过滤，选择关键节理。通常从钻孔和节理图开始，需要一个反复过程。选择关键节理的困难集中在确定性的UDEC模型中如何表征这些统计数据。对于多组节理组，应当假设给定的一组节理，其产状是相同的，且变化较小（比如说小于10o到15o）。节理岩体范围可由特定的节理被划分成子区域，且在子域内的节理被认为是连续的。这是对有限的长度和不连续性特性影响的一个上限估计。 下一步是应用一个滤波准则来识别那些在给定的荷载作用下，容易滑移或张开的节理。该准则可能涉及范围从：（1）识别是否提供充分的自由度允许滑移，到（2）利用块体理论来检查节理位移的可动性，到（3）与现场观察和记录进行比较－例如，从微震监测识别关键节理。 3.9. 2 设计模型 模拟者总是试图建立尽可能考虑详细地质结构的UDEC模型。反对这种研究的主要理由在于：（1）期望采用充分的数据模拟节理岩体的每一细节是无用的；（2）对于工程问题，计算机硬件需要一个快速模型，接近典型的工程实际；（3）最重要的是，考虑得较为详细节理对模拟结果的理解与控制的效率会更低。 产生UDEC模型时应考虑两个方面：首先是不连续分析是否真正的需要。这在大部分情况下，取决于研究实际模型尺寸与节理的平均间距之比。例如，如果包含一组平均间距不大于1m的岩体中开挖，开挖的最小尺寸是10m，然而，用堆砌节理材料模型进行连续分析可能是合理的。在这种情况下，连续分析产生的响应与模拟节理的计算结果在总体上是等效的。用UDEC模型分析给出破坏机理较为详细分析，但比连续分析花费更多的计算时间。在实际模型尺寸与节理间距之比值大于10:1的情况下，
UDEC常用命令
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 jset 90,0 4,0 4,0 6,0 0,-50 range jreg 3 ;jset 倾角,0 线段长,0 线段与线段轴向间隔长,0 垂向间距,0 xm,ym range jregion n 每一数据对中的第一个值是均值，而第二个是对应于均值的最大均方差。第一组数据对是节理迹线与x-坐标轴的正方向的夹角。第二对数据节理迹线长度；第三组数据是不连续节理的岩桥长度；第四组数据是节理间距。还有一些选择参数，可以用于产生一组比较复杂的节理模式。 ;其中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