臺(tái)階法隧道開挖

1、隧道建模是在ansys中完成的，計(jì)算部分在Flac3D中完成。以下是具體內(nèi)容。因?yàn)闀r(shí)間倉促， 未做注解，有疑問可關(guān)注新浪微博：巖土 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 上臺(tái)階襯砌 range group 1group 上臺(tái)階土體 range gr

13、oup 3group 下臺(tái)階襯砌 range group 2group 下臺(tái)階土體 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自重應(yīng)力場(chǎng).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 上臺(tái)階土體 y y0 y1model null range reg 上臺(tái)階襯砌 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 下臺(tái)階土體 y yx0 yx1model null range reg 下臺(tái)階襯砌 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 下臺(tái)階土體 y 19 20model null range reg 下臺(tái)階襯砌 y 19 20sel shell id 100 group 圍巖 range x