台阶法隧道开挖

1、隧道建模是在ansys中完成的，计算部分在Flac3D中完成。以下是具体内容。因为时间仓促， 未做注解，有疑问可关注新浪微博：岩土 1112Ansys建模部分命令流finish/clear/filename,3D*afun,deg/prep7et,1,mesh200,6et,2,solid45mp,ex,1,23e9mp,prxy,1,0.2mp,dens,1,2200mp,ex,2,23e9mp,prxy,2,0.2mp,dens,2,2200mp,ex,3,0.18e9mp,prxy,3,0.35mp,dens,3,1950tb,dp,3tbdata,1,0.18e6,23mp,ex,4,

2、0.18e9mp,prxy,4,0.35mp,dens,4,1950tb,dp,4tbdata,1,0.18e6,23mp,ex,5,0.18e9mp,prxy,5,0.35mp,dens,5,1950tb,dp,5tbdata,1,0.18e6,23k,1,0,0,0k,2,0,11.37,0k,3,3.34,0,0k,4,-3.34,0,0k,101,0,0,1k,102,0,11.37,1k,103,3.34,0,1k,104广3.34,0,1circle,1,5.8,101circle,2,14.4,102circle,3,2.67,103circel,4,2.67,104lcsl,a

3、ll/pnum,line,1 lsel,s,17,22,1 lsel,a,46,48,1 lsel,a,27,28,1 lsel,a,43 lsel,inve ldele,all lsel,all lcomb,19,21 lcomb,17,18 lcomb,27,43 lcomb,28,46 cm,l1,line circle,1,6.1,101 circle,2,14.7,102 circle,3,2.97,103 circle,4,2.97,104 cmsel,s,l1,line lsel,inve lcsl,all lsel,s,23,26,1 lsel,a,18,21,3 lsel,a

4、,51,56,5 lsel,a,33,34,1 lsel,a,54,55,1 cm,l2,line cmsel,a,l1,line lsel,inve ldele,all,1 lsel,all lcomb,23,25 lcomb,18,21 lcomb,33,51 lcomb,34,54 circle,1,8.8,101 circle,2,17.4,102 circle,3,5.67,103 circle,4,5.67,104 cmsel,s,l1,line cmsel,a,l2,line lsel,inve lcsl,all lsel,s,line,29,32,1 lsel,a,line21

5、,25,4 lsel,a,line,63,70,7 lsel,a,line,39,40,1 lsel,a,line,66,67,1 cm,l3,line cmsel,a,l1,line cmsel,a,l2,line lsel,inve ldele,all allsel lcomb,29,31 lcomb,21,25 lcomb,39,63 lcomb,40,66 kdele,101,104,1 rectng,-50,50,30,-45 asbl,all,all /pnum,area,1 allsel csys,0 wpcsys wpoff,1.8*8.8 wprota,90 asbw,all

6、wpcsys wprota,90 asbw,allwpcsys wpoff广1.8*8.8 wprota,90 asbw,allwpcsys wpoff,1.7*8.8 wprota,90 asbw,allwpcsys wprota,90 asbw,allwpcsys wpoff-1.4*8.8 wprota,90 asbw,allwpcsys allselwpcsys allsel lesize,68,6,1 lsel,s,line2,4,2 lsel,a,line,7 lesize,all,12 asel,s,area,1 type,1 mshape,0 mshkey,1 amesh,al

7、l asel,s,area,2 type,1 mshape,0 mshkey,1 amesh,all allsel numcmp,all lesize,20,4 lesize,45,4 lesize,18,4 lesize,44,4 lesize,26,14 lesize,25,14 lccat,20,45 lccat,18,44 type,1 mshape,0 mshkey,1 amesh,21,23,2 lsel,s,lccat ldele,all allsel lesize,21,12,1 lesize,16,12,1 lsel,s,line,55,56,1 lsel,a,line,13

8、 lesize,all,1 type,1 mshape,0 mshkey,1 amesh,15,16,1 allsel type,1 mshape,0mshkey,0 amesh,25,26,1 allsellsel,s,line,57,58 lsel,a,line,23 lesize,all,4,1.8 lsel,s,line,19,27,8 lesize,all,12 type,1mshape,0 mshkey,1 amesh,17,18,1 allsellcomb,30,65 lcomb,28,62 lesize,30,11 lesize,28,11 lesize,36,11 lesiz

9、e,37,11 lesize,29,4,1.8 type,1mshape,0 mshkey,1 amap,28,16,20,10,14 amap,27,16,20,8,11 allsel lsel,s,line,66,67,1 lsel,a,line50,76,26 lesize,all,8,1.8 lsel,s,line48,49lsel,a,line60,63,3 lesize,all,6 lsel,s,line70,71 lsel,a,line,73,75,2 lesize,all,11 allsellccat,49,63 lccat,48,60 lccat,71,75 lccat,70

10、,73 lccat,30,36 lccat,28,37 type,1 mshape,0mshkey,1amesh,22,24,2amesh,10,11,1lsel,s,lccatldele,allallselasel,s,loc,x,-100*8.8,-1.8*8.8lsla,slsel,r,tan1,y,0lesize,all,16asel,s,loc,x,1.8*8.8,100*8.8lsla,slsel,r,tan1,y,0lesize,all,16asel,s,loc,y,1.7*8.8,100*8.8lsla,slsel,r,tan1,x,0lesize,all,8asel,s,lo

11、c,y,-100*8.8,-1.4*8.8lsla,slsel,r,tan1,x,0lesize,all,16asel,s,loc,x,-1.8*8.8,1.8*8.8asel,r,loc,y,-1.4*8.8,1.7*8.8asel,invetype,1mshape,0mshkey,1amesh,allallselaplotk,1000,0,0,-201,1000/view,1,1,1,1/replplotaplotextopt,esize,20extopt,aclear,1type,2mat,3asel,s,2vdrag,all62type,2mat,4asel,s,21,23,2vdra

12、g,all62type,2mat,1asel,s,15,16vdrag,all62type,2mat,2asel,s,25,26vdrag,all62type,2mat,5asel,s,17,20,1asel,a,27,28,1asel,a,22,24,2asel,a,8,11,1asel,a,3,7,1asel,a,12,14,1vdrag,all62allseleplot/eshape,1wpcsys,1,0wprota,0,-90,0csys,wpFlac3D命 令流plot block groupgroup 上台阶衬砌 range group 1group 上台阶土体 range gr

13、oup 3group 下台阶衬砌 range group 2group 下台阶土体 range group 4group 围岩 range group 5mo mohrprop bulk 16.67e6 shear 7.69e6 fric 28 coh 22e3 ten 0 dil 0fix x range x -50.01 -49.99fix x range x 49.99 50.01fix y range y -0.01 0.01fix y range y 19.99 20.01fix z range z -45.01 -44.99ini density 1380set grav 0 0

14、-10solve elassave自重应力场.savset log onplot block groupdef gdsc=0endgdsdef cab_parmarray x1_zhu(1,13) z1_zhu(1,13) x2_zhu(1,13), z2_zhu(1,13)loop n(1,13)x1_zhu(1,n)=6.1*cos(15*(n-1)*degrad)z1_zhu(1,n)=6.1*sin(15*(n-1)*degrad)x2_zhu(1,n)=9.1*cos(15*(n-1)*degrad)z2_zhu(1,n)=9.1*sin(15*(n-1)*degrad)end_lo

15、opendcab_parmdef ini_zhx_d=x1_zhu(1,n)z_d=z1_zhu(1,n)x_2d=x2_zhu(1,n)z_2d=z2_zhu(1,n)enddef ini_maoloop m(3,15)n=m-2ini_zhcommandsel cable id m begin x_d c z_d end x_2d c z_2d nseg 3sel cable prop emod 45e9 xcarea 1.57e-3 gr_per 1. yten 25e4 gr_k 17.5e6 gr_c 20e4 ran id mend_commandend_loopenddef ex

16、cavloop cut(1,20)y0=cut-1y1=cutc=cut-0.5yx0=y0-1yx1=y1-1commandmodel null range reg 上台阶土体 y y0 y1model null range reg 上台阶衬砌 y y0 y1end_commandini_maocommandsel shell id 100 group 围岩 range x -6.11 6.11 z -0.01 6.11 cyl endl 0 y0 0 end2 0 y1 0 & radius 6.1sel delete shell range y -0.01 0.01sel shell p

17、rop iso 10.5e9 0.25 thick 0.3 density 2500 ran y y0 y1end_commandif cut1 thencommandmodel null range reg 下台阶土体 y yx0 yx1model null range reg 下台阶衬砌 y yx0 yx1sel shell id 100 group 围岩 range x -6.11 6.11 z -3.34 0.01 cyl end1 0 yx0 0 end2 0 yx1 0 &radius 6.11sel shell prop iso 10.5e9 0.25 thick 0.3 density 2500 ran y yx0 yx1end_commandend_ifif cut=20 thencommandmodel null range reg 下台阶土体 y 19 20model null range reg 下台阶衬砌 y 19 20sel shell id 100 group 围岩 range x