change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围
;参数设置
prop mat=1 dens=2720 ;物质的密度 zone k=12.12e9 g=10.26e9 fric=33 coh=5.77e6 ten=3.6e6 range mat=1 ;定义物质1的体积模量、剪切模量、内摩擦角、粘聚力、抗拉强度
prop jmat=1 jkn=11e8 jks=12e7 jfric=12 jcoh=0 jten=0
;定义节理的属性法向刚度、切向刚度、内摩擦角、粘聚力、抗拉强度
1) D--质量密度 F--摩擦角 B--体积模量 C--内聚力(粘结力) S--剪切模量
2) 应力--正号代表张力，负号代表压力； 应变--正的应变代表伸长，负的应变代表压缩 重力--正号的重力物质往下拉，负号的重力将物质往上提
！！！密度必须设置density for block materical 1 is zero cannot cycle----块体密度未设
2.4 初始应力 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 ;采用位移法固定边界，让其速度=0。
;边界条件
boun stress 0 0 0 range 0,300 99.9 100.1 ;定义边界条件，三向应力为0，上边界，范围，纵向波动范围0 0 0（x、xy、y）
boun stress 0 0 -8.6095397e6 range ;定义边界应力条件，竖直方向的应力为，负号表示方向向下 insitu str -7.5e6 0 -15e6 szz -7.5e6 ygrad 1.1018e4 0 2.203607e4 ;定义应力属性，sxx方向即水平应力，sxy方向即剪切应力，syy方向即竖直应力，szz方向的应力，ygrad应力梯度
bound xvel=0 range -0.1 0.1 0 100 ;X位移边界,位移波动范围，范围 bound xvel=0 range 299.9 300.1 0 100 ;X位移边界，同上 bound yvel=0 range -0.1 300.1 -0.1,0.1 ;Y位移边界
hist solve_ratio type 1 ;设置不平衡力的精度，普适
6
set jmatdf 2 ;定义节理属性，使符合摩尔库伦准则普适。2代表摩尔库伦准则
hist xdisp 10,20 ;测该点的X方向位移 hist sxx 10,20 ;测该点X方向的应力 hist syy 10,20 ;测该点Y方向的应力 reset disp
reset hist ;位移清零，历史记录清零 his unbal
solve ratio 1e-5 ;记录不平衡力，计算直到精度达到1e-5
;执行计算
solve\\step 5000\\cycle 5000
;进行开挖
;删除水平块体，开挖，不可再充填
delete range x1 x2 y1 y2等价于delete range reg x1,y1 x1,y2 x2,y2 x2,y1
;删除不规则四边形
delete range region x1,y1 x2,y2 x3,y3 x4,y4
3.5 支护部分 ;加支柱，支架
support x,y angle segment witdth mat
;support是指加入支柱，x,y表明支柱所处的位置点。angle表明支柱的轴向与x轴正方向的夹角，以逆时针方向为正。segment指一组支柱的个数，默认值为5。witdth指一组支柱支撑的范围。 ;mat是指赋予支柱的材料性质属于几号？下面的table指支柱的增阻特性。 比如：
sup 10,20 angle=90 seg=5 width=2 mat=6 prop mat=6 sup_kn=-1
table 1 0 3e6 0.005 5e6 0.015 25e6 0.5 40e6
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)特性曲线，后面为四个点的坐标。目前，液压支柱的最大支撑力为12*106N 显示支架：pl bl support red
7
;加锚杆
cable x1,y1 x2,y2 npoint mats asteel matg 等
;x1,y1 x2,y2指锚杆的两个端点;npoint指锚杆的粘结点数;mats指锚杆材料性质,材料编号，按着前面给定的材料号继续往下编； asteel指锚杆横截面积；matg指锚固体性质；
cable (hc1bx,hc1by) (hc1ex,hc1ey) 25 9 314e-6 10 cable (hc2bx,hc2by) (hc2ex,hc2ey) 25 9 314e-6 10 ;锚杆参数
prop m=9 cb_dens 7500 cb_ycomp 6.3e8 cb_yield 0.5e7 cb_ymod 0.5e11 ;锚固体的参数
prop m=10 cb_kbond=1.6e9 cb_sbond=2e6
cable (73.4,2.91) (71.15,4.99) 10 12 380e-6 11 ;锚杆，起始点坐标，后面四个数值可固定 prop mat 12 cb_dens 7800 cb_ycomp 430e6 cb_yield 160e3 cb_ymod 1.3e11 ;锚固体的性质
prop mat 11 cb_kbond 6.3e9 cb_sbond 6e5 ;定义锚杆的属性，密度、抗拉强度、屈服强度
；巷道断面加梯子梁的命令
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-界面
struct gen xc=20 yc=40.24 np=100 fang=-40 theta=260 mat=14 thick=0.2 prop mat=14 st_ymod=13.5e9 st_prat=0.14 st_den=7800 prop mat=14 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7 prop mat=14 if_kn=1.35e9 if_kf=1.35e9 if_tens=0 if_fric=18 if_coh=0 梯形梁b(xc,yc)锚杆a(x1,y1)(x2,y2)巷 道 图形显示
pl bl hist 1,2,... 观测点位移 pl bl nu(mber) 显示块体标号
8
pl bl cab red stru red cab--锚杆 stru--梁 pl bl sxx/syy/szz 显示应力图 pl bl dis 位移矢量图 pl bound pl 显示塑性区 pl bound dis 显示位移矢量图
pl hold bl num 在屏幕上将显示建立的块体单元，并在按
plot block xdis ;显示块体的位移(有x和y方向) plot block pl ;显示块体的塑性区(plastic)
;屏幕观看下沉与应力等值线,int是指interval pl bl ydisp int 0.1 pl bl syy int 1e6
udec 颜色
关键词 black blue green cyan red magenta orange brown dgray
图形输出
;抓取图像，在autocad内打开、编辑 ;输出应力分布 pl bl pline 1 syy red set pl dxf 25 set out f:\\syy100.dxf
pl pen bl pline 1 syy red ；=copy f:\\syy100.dxf
1.把图形保存下来，能在AUTOCAD中
plot block cable red supp ye stru bl;显示块体、锚杆（红色）、支架（黄色）、梁（蓝色）的图形
9
颜 色 黑 蓝 绿 蓝绿 红 品红 橙 褐 灰
关键词 lgray lblue lgreen lcyan lred lmagenta lorange yellow white
颜 色 浅灰 浅蓝 浅绿 浅蓝绿 浅红 浅品红 浅橙 黄 白
set plot 256 dxf；设置图形为256色 set out c:\\ss.dxf
copy c:\\ss.dxf；把这图形以ss.dxf文件保存在C磁盘下
注！！！：开始为dxf格式，移动后另存改为dwg格式，可以减少很大的空间！ 动画制作
;打开图象抓拍功能，设置图片间隔时间以及窗口大小 movie on movie step 1000 window x1 x2 y1 y2 ;movie off
【注：只显示改区间内的内容,分别为x1，x2和y1，y2的坐标！window 6 27.5 -1 13 pl bl syy
movie file kaicai size 800 600 on step 1000 del ran 40 70 12.0 15 step 10000 movie off
；电影将被保存至C:\\Program Files\\Itasca\?10下面，
打开方式C:\\Program Files\\Itasca\\shared\\Utility\\movie，必须点击start才开始播放;
【注：1)：电影的制作过程实际上为抓频的过程，因此，在此期间最好将其它的程序先关掉！ 2): 在进行电影之前，输入pl bl syy,则整个过程将显示应力的变化过程，如果输入pl bl pline syy red
则整个过程将显示测线应力变化的过程。输入先输入pl bl syy 后输入pl bl pline syy red，则将在电影中同时出现！输入pl bound，则将出现裂隙带的范围！
10