1、yuanyan.log print pline 1 syy 【记录测线1应力】 print pline 2 ydisset log off ;设置观测线set pline 1 x1,y1 x2,y2 n (n-观测线分的段数)set pline 25,20 300,20 10 ;定义测线,起始点坐标,测线分段2.2 图形划分块体命令 block x1,y1 x2,y2 x3,y3 x4,y4 【建立模型框架,四个坐标角点必须按“顺时针”排列,也可以为五个坐标点。】 图2 5个坐标点 pl bl显示块体;pl hold bl =pl bl +pause屏幕上直接显示划分的块体,关闭显示后将继续运
2、行; del range x1 x2 y1 y2删除块体;fix range x1 x2 y1 y2 固定块体不动画圆弧arc xc,yc xb,yb theta ncracks【xc,yc圆弧中心,xb,yb指画圆弧的起始点,theta沿逆时针方向圆弧的角度;ncracks指构成指圆弧的折线段数,数目越多,圆弧越圆; 圆弧命令必须配合Crack命令,但tunnel 命令可以在输入后直接显示为圆形;先划分巷道,然后划分其它岩层节理】 ;开挖巷道节理,形成巷道空间,巷道必须先开挖出来arc 52.17 85.59 53.97 85.59 180 8cra 50.37 83.59 50.37 85
3、.59 cra 50.37 83.59 53.97 83.59 cra 53.97 83.59 53.97 85.59建圆形隧道、圆形巷道tunnel xc,yc r n【xc,yc圆的中心,r半径,n构成圆的线段数。tunnel 4,5 2 45 先输tunnel,然后输crack,裂隙不贯穿隧道;隧道开挖:方法一:删除编号pl bl num delete block 1920 方法二:delete range -1 1 -1 1方法三:指定窗口 windows x1 x2 y1 y2 del range window,不一定能全部删除,可以继续删除编号方法四:定义区域,删除区域,jreg
4、n x1,y1 x2,y2 x3,y3 x4,y4 delete jreg n】节理划分、定义crack x1,y1 x2,y2【划分节理后,最好输入pause暂停,pl bl显示划分的块体是否正确或文本中输入pl hold bl】jregion id n x1,y1 x2,y2 x3,y3 x4,y4 【节理区域定义,赋予节理区域地址号,可以删除任何区域delete jreg n 】jset 90,0 4,0 4,0 6,0 0,-50 range jreg n【在定义的节理区域内进行节理的再划分;jset 、t、g、s (xo,yo) n :倾角,节理长度, 岩桥长度,垂距,起始点坐标,n
5、区域标号】 图3 节理组定义;不规则块体划分 jreg id 4 14.65 0 73.0 160 78.16 160 19.96 0 ;定义节理区域4change mat 4 range 0 100 -65 -50 ;定义vor edg 5 iter 20 rou 0.1 ran jreg 4【edge平均边长,最少20倍rounding长度;iterations 迭代次数,默认5;round,产生块体的圆角长度,默认0.5;rtol 容许偏差,默认0.02*round;seed n】2.3 参数定义gen quad 10 range xl xu yl yu在指定的区域生成一定宽度的单元zo
6、ne 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
7、xl xu yl yu 指定各岩层的节理标号prop mat=1 dens=2000 prop mat=2 dens=2650prop mat=3 dens=2700 指定各材料的密度,比如1号材料dens=2000(kg/m3),即1立方米重2吨zone k=0.15e9,g=0.1e9,fric=10.00,coh=0.19e6,ten=0.09e6 range mat=1zone k=2.8e9,g=2.2e9,fric=30.00,coh=1.5e6,ten=0.4e6 range mat=2zone k=6.9e9,g=6.6e9,fric=38.62,coh=5.63e6,ten=
8、3.20e6 range mat=3k为材料的法向刚度,g为材料的切向刚度,friction为材料的内摩擦角,cohesion为材料的内聚力,tension为材料的抗拉强度prop jmat=1 jkn=0.2e8,jks=0.1e7,jcoh=0,jfric=4,jten=0prop jmat=2 jkn=8e8,jks=5e7,jcoh=0.1e6,jfric=8,jten=0prop jmat=3 jkn=20e8,jks=16e7,jcoh=0.4e6,jfric=15,jten=0 jkn为节理的法向刚度,jks为节理的切向刚度,jfriction为节理的内摩擦角,jcohesion
9、为节理的内聚力,jtension为节理的抗拉强度块体力学参数k-体积模量;g-剪切模量;fric-摩擦角;coh-粘结力(内聚力);ten-抗拉强度;d-体积力接触面力学参数jkn-法向刚度;jks-切向强度;jc-粘结力;jf-摩擦角;jt-抗拉强度jreg id 1 0,0 0,15 300,15 300,0 delete ;定义节理,命名1,定义范围,四点坐标,delete为常规语言 jset 90,0 5,0 5,0 5,0 0,0 rang jreg 1 ;定义节理角度90节理的长度;节理的距离,即纵向间隔;横向距离,即横向隔5一个;起始点坐标为0,0jset 90,0 5,0 5,
10、0 5,0 2.5,5 rang jreg 1 ;定义节理的另外一项表1 岩石块体参数岩石名称 弹性模量E泊松比体积模量K剪切模量G砂岩19.30.3826.87.0粉砂岩26.30.2215.610.8石灰岩28.50.2922.611.1页岩、泥岩8.84.3大理岩55.80.2537.222.3花岗岩73.843.930.2内摩擦角粘聚力抗拉强度27.827.21.1732.134.714.438.4石英岩48.070.642.06.721.5851.055.1玄武岩31.066.213.1体积模量K和剪切模量G之间的关系: gen quad 10 range 0 300 0 100 ;
11、当0 300 0 100范围内的块体的对角长度小于10时就是可变形体。zone model mo range 0 300 0 100 ;定义摩尔库伦模型的范围,X的范围,Y的范围change jcons=5 range 0 300 0 100 ;定义节理的某个属性change mat=1 range reg 0,0 0,15 300,15 300,0 ;定义物质1的范围,即赋予这个范围为物质1,只赋名,没有定义属性change mat=2 range reg 0,15 0,20 300,20 300,15change mat=3 range reg 0,20 0,23 300,23 300,2
12、0change mat=4 range reg 0,23 0,28 300,28 300,23change jmat=1 range reg 0,0 0,15 300,15 300,0 ;定义节理1的范围set jmatdf 2 ;定义节理属性,使符合摩尔库伦准则普适。2代表摩尔库伦准则 参数设置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
13、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.8bound xvel=0 range -0.1 3.00
14、-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.1
15、018e4 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 ;设置不平衡力的精度,普适hist xdisp 10,20 ;测该点的X方向位移hist sxx 10,20 ;测
16、该点X方向的应力hist syy 10,20 ;测该点Y方向的应力reset dispreset hist ;位移清零,历史记录清零his unbalsolve ratio 1e-5 ;记录不平衡力,计算直到精度达到1e-5执行计算solvestep 5000cycle 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,y43.5 支护部分加支柱,支架support x
17、,y angle segment witdth matsupport是指加入支柱,x,y表明支柱所处的位置点。angle表明支柱的轴向与x轴正方向的夹角,以逆时针方向为正。segment指一组支柱的个数,默认值为5。witdth指一组支柱支撑的范围。mat是指赋予支柱的材料性质属于几号?下面的table指支柱的增阻特性。比如:sup 10,20 angle=90 seg=5 width=2 mat=6prop mat=6 sup_kn=-1 table 1 0 3e6 0.005 5e6 0.015 25e6 0.5 40e64.液压支架的命令supp xc yc wid l seg n ma
18、t j; (xc,yc)是指中心点坐标,l指支架的宽度,n指分段数,mat j指支柱材料性质为jprop 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加锚杆cable x1,y1 x2,y2 npoint
19、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_k
20、bond=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
21、段,a b 指角度prop mat=16 st_ymod=13.5e9 st_prat=0.14 st_den=7800prop mat=16 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7prop 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.2prop mat=14 st_ymod=13.5e9 st_prat=0.14
22、 st_den=7800prop mat=14 st_yield=6e7 st_yresid=6e7 st_ycomp=2.5e7prop mat=14 if_kn=1.35e9 if_kf=1.35e9 if_tens=0 if_fric=18 if_coh=0图形显示pl bl hist 1,2,. 观测点位移pl bl nu(mber) 显示块体标号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
23、 num 在屏幕上将显示建立的块体单元,并在按键后退出后,或直接继续执行命令 。相当于pau+pl bl num;删除块体 del bl +地址号plot block stress ;显示块体的应力plot block xdis ;显示块体的位移(有x和y方向)plot block pl ;显示块体的塑性区(plastic)屏幕观看下沉与应力等值线,int是指intervalpl bl ydisp int 0.1pl bl syy int 1e6udec 颜色关键词颜 色black黑lgray浅灰blue蓝lblue浅蓝green绿lgreen浅绿cyan蓝绿lcyan浅蓝绿red红lred浅
24、红magenta品红lmagenta浅品红orange橙lorange浅橙brown褐yellow黄dgray灰white白图形输出抓取图像,在autocad内打开、编辑输出应力分布pl bl pline 1 syy redset pl dxf 25set out f:syy100.dxfpl pen bl pline 1 syy red ;=copy f:1.把图形保存下来,能在AUTOCAD中 plot block cable red supp ye stru bl;显示块体、锚杆(红色)、支架(黄色)、梁(蓝色)的图形 set plot 256 dxf;设置图形为256色 set out
25、 c:ss.dxf copy c:ss.dxf;把这图形以ss.dxf文件保存在C磁盘下注:开始为dxf格式,移动后另存改为dwg格式,可以减少很大的空间!动画制作打开图象抓拍功能,设置图片间隔时间以及窗口大小movie on movie step 1000window x1 x2 y1 y2movie off【注:只显示改区间内的内容,分别为x1,x2和y1,y2的坐标!window 6 27.5 -1 13pl bl syy movie file kaicai size 800 600 on step 1000del ran 40 70 12.0 15step 10000电影将被保存至C:
26、Program FilesItascaudec310下面, 打开方式C:Program FilesItascasharedUtilitymovie,必须点击start才开始播放;1):电影的制作过程实际上为抓频的过程,因此,在此期间最好将其它的程序先关掉! 2): 在进行电影之前,输入pl bl syy,则整个过程将显示应力的变化过程,如果输入pl bl pline syy red则整个过程将显示测线应力变化的过程。输入先输入pl bl syy 后输入pl bl pline syy red,则将在电影中同时出现!输入pl bound,则将出现裂隙带的范围!4 实例分析: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
copyright@ 2008-2023 冰点文库 网站版权所有
经营许可证编号:鄂ICP备19020893号-2