地鐵車站FLAC命令流

1、new title 樁錨支護基坑初始應力模擬 ;調入模型重新指定分組 impgrid 10mmodel.flac3D generate merge 0.02 group 粉質粘土 ran group 3 group 雜填土 ran group 1 group 粉土 ran group 2 group 中砂 ran group 4 group 粉粘 ran group 5 group 開挖1 ran group 11 group 開挖2 ran group 10 group 開挖3 ran group 9 group 開挖4 ran group 8 group 開挖5 ran group 7 g

2、roup 開挖6 ran group 6 ;指定 材料模型 model elas define mqb1 ;雜填土 me1=7.8e7 mv1=0.35 ; mbulk1=me1/(3.0*(1-2.0*mv1 mshear1=me1/(2.0*(1+mv1 ;粉土 me2=6.36e7 mv2=0.4 ; mbulk2=me2/(3.0*(1-2.0*mv2 mshear2=me2/(2.0*(1+mv2 ;粉粘上 me3=8e7 mv3=0.35 ; mbulk3=me3/(3.0*(1-2.0*mv3 mshear3=me3/(2.0*(1+mv3 ;中砂 me4=20e7 mv4=0.

3、4 ; mbulk4=me4/(3.0*(1-2.0*mv4 mshear4=me4/(2.0*(1+mv4 ;粉粘下 me5=9.74e7 mv5=0.3 ; mbulk5=me5/(3.0*(1-2.0*mv5 mshear5=me5/(2.0*(1+mv5 end mqb1 ;定義材料參數(shù) property density 1850 bulk mbulk1 shear mshear1 ran group 雜填土 any group 開挖1 any property density 1950 bulk mbulk2 shear mshear2 ran group 粉土 any group

4、開挖2 any property density 1980 bulk mbulk3 shear mshear3 ran group 粉質粘土 any group 開挖3 any property density 1980 bulk mbulk3 shear mshear3 ran group 開挖4 any group 開挖5 any property density 1980 bulk mbulk3 shear mshear3 ran group 開挖6 property density 2000 bulk mbulk4 shear mshear4 ran group 中砂 property

5、 density 1970 bulk mbulk5 shear mshear5 ran group 粉粘 ;定義邊界條件 fix x y z ran z 0 fix x ran x 0 fix x ran x 140 fix y ;初始應力狀態(tài) set gravity 0 0 -10 history unbal solve ;塑性計算 model mohr property density 1850 bulk mbulk1 shear mshear1 friction 15 cohesion 1e4 tension 0 ran group 雜填土 any group 開挖1 any prope

6、rty density 1950 bulk mbulk2 shear mshear2 friction 28 cohesion 2.46e4 tension 0 ran group 粉土 any group 開挖2 any property density 2000 bulk mbulk4 shear mshear4 friction 35 cohesion 0 tension 0 ran group 中砂 property density 1970 bulk mbulk5 shear mshear5 friction 19.2 cohesion 2.41e4 tension 0 ran gr

7、oup 粉粘 property density 1980 bulk mbulk3 shear mshear3 friction 16.3 cohesion 2.31e4 tension 0 ran group 粉質粘土 any group 開挖3 any property density 1980 bulk mbulk3 shear mshear3 friction 16.3 cohesion 2.31e4 tension 0 ran group 開挖4 any group 開挖5 any property density 1980 bulk mbulk3 shear mshear3 fric

8、tion 16.3 cohesion 2.31e4 tension 0 ran group 開挖6 solve save initial.sav ;調入初始應力計算結果 rest initial.sav title 基坑分步開挖支護效果 set log on set logfile jisuanguocheng.log initial xdisp 0 ydisp 0 zdisp 0 initial xvel 0 yvel 0 zvel 0 initial sta 0 ;開挖; ; ;設置監(jiān)控變量 ;豎向位移 hist id=2 gp zdisp 93.1 5 50 hist id=3 gp z

9、disp 89.1 5 50 hist id=4 gp zdisp 84.1 5 50 hist id=5 gp zdisp 79.1 5 50 hist id=6 gp zdisp 74.1 5 50 hist id=7 gp zdisp 69.1 5 50 hist id=8 gp zdisp 59.1 5 50 hist id=9 gp zdisp 49.1 5 50 hist id=10 gp zdisp 39.1 5 50 hist id=11 gp zdisp 29.1 5 50 ;樁的水平位移 hist id=12 gp xdisp 93.1 5 50 hist id=13 gp

10、 xdisp 100 0 43 hist id=14 gp xdisp 100 0 41 hist id=15 gp xdisp 100 0 39 hist id=16 gp xdisp 100 0 37 hist id=17 gp xdisp 100 0 35 hist id=18 gp xdisp 100 0 33 hist id=19 gp xdisp 100 0 31 hist id=20 gp xdisp 100 0 29 hist id=21 gp xdisp 100 0 27 hist id=22 gp xdisp 100 0 25 hist id=23 gp xdisp 100

11、 0 23 hist id=24 gp xdisp 100 0 21 ;第一步開挖; model null ran group 開挖1 ;安裝第一道土釘 ;定義土釘相關參數(shù) def canshu ;-可以修改的量- _mLength=9.0 ;錨桿1,2長度 _mLength3=6.0 ;錨桿3長度 _mAngle=5.0 ;錨桿與水平方向的夾角 _msInterval=2.0 ;錨桿水平方向的間隔 _mzInterval=1.5 ;錨桿縱向的間隔 _mnseg=9 ;12錨桿劃分單元數(shù) _mnseg3=6 ;3錨桿劃分單元數(shù) _zm1t=48 ;第一排錨桿頂點高度，每排錨桿是在同一高度的。

12、_zm2t=46.5 ;第二排錨桿頂點高度 _zm3t=45 ;第三排錨桿頂點高度 _x1mt=95.1 ;第一排錨桿的頂點x坐標 _x2mt=95.847 ;第2排錨桿的頂點x坐標 _x3mt=96.6 ;第3排錨桿的頂點x坐標 ;-源于上面參數(shù)- _xVariation=_mLength*cos(_mAngle*pi/180 ;1,2錨桿x方向的增量,FLCA內置三角函數(shù)是以弧度制計算的 _yVariation=_mLength*cos(_mAngle*degrad ;1,2錨桿y方向的增量,degrad=pi/180,內置常規(guī)變量，弧度與度的轉換 _zVariation=_mLength

13、*sin(_mAngle*degrad ;1,2錨桿z方向的增量 _xVariation3=_mLength3*cos(_mAngle*pi/180 ;3錨桿x方向的增量,FLCA內置三角函數(shù)是以弧度制計算的 _yVariation3=_mLength3*cos(_mAngle*degrad ;3錨桿y方向的增量,degrad=pi/180,內置常規(guī)變量，弧度與度的轉換 _zVariation3=_mLength3*sin(_mAngle*degrad ;3錨桿z方向的增量 _x1mb=_x1mt- _xVariation ;第一排錨桿的底點x坐標 _zm1b=_zm1t-_zVariatio

14、n ;第一排錨桿的底點高度 _x2mb=_x2mt- _xVariation ;第2排錨桿的底點x坐標 _zm2b=_zm2t-_zVariation ;第2排錨桿的底點高度 _x3mb=_x3mt- _xVariation3 ;第3排錨桿的底點x坐標 _zm3b=_zm3t-_zVariation3 ;第3排錨桿的底點高度 end canshu sel cable id=1 begin _x1mt 0 _zm1t end _x1mb 0 _zm1b nseg=_mnseg sel cable id=1 begin _x1mt 2 _zm1t end _x1mb 2 _zm1b nseg=_m

15、nseg sel cable id=1 begin _x1mt 4 _zm1t end _x1mb 4 _zm1b nseg=_mnseg sel cable id=1 begin _x1mt 6 _zm1t end _x1mb 6 _zm1b nseg=_mnseg sel cable id=1 begin _x1mt 8 _zm1t end _x1mb 8 _zm1b nseg=_mnseg sel cable id=1 begin _x1mt 10 _zm1t end _x1mb 10 _zm1b nseg=_mnseg ;定義參數(shù) sel cable id 1 prop xcarea

16、6.2e-4 emod 200e9 yte ns 5.4e5 ycomp 1.0e5 sel cable id 1 prop gr_k 2.0e7 gr_coh 1.0e6 gr_per 0.377 gr_fric 25.0 plo con zdisp solve save kaiwa1.sav rest kaiwa1.sav ;第二步開挖 model null ran group 開挖2 sel cable id=2 begin _x2mt 0 _zm2t end _x2mb 0 _zm2b nseg=_mnseg sel cable id=2 begin _x2mt 2 _zm2t end

17、 _x2mb 2 _zm2b nseg=_mnseg sel cable id=2 begin _x2mt 4 _zm2t end _x2mb 4 _zm2b nseg=_mnseg sel cable id=2 begin _x2mt 6 _zm2t end _x2mb 6 _zm2b nseg=_mnseg sel cable id=2 begin _x2mt 8 _zm2t end _x2mb 8 _zm2b nseg=_mnseg sel cable id=2 begin _x2mt 10 _zm2t end _x2mb 10 _zm2b nseg=_mnseg ;定義參數(shù) sel c

18、able id 2 prop xcarea 6.2e-4 emod 200e9 ytens 5.4e5 ycomp 1.0e5 sel cable id 2 prop gr_k 2.0e7 gr_coh 1.0e6 gr_per 0.377 gr_fric 25.0 plo con zdisp solve save kaiwa2.sav rest kaiwa2.sav ;第三步開挖 model null ran group 開挖3 sel cable id=3 begin _x3mt 0 _zm3t end _x3mb 0 _zm3b nseg=_mnseg3 sel cable id=3 b

19、egin _x3mt 2 _zm3t end _x3mb 2 _zm3b nseg=_mnseg3 sel cable id=3 begin _x3mt 4 _zm3t end _x3mb 4 _zm3b nseg=_mnseg3 sel cable id=3 begin _x3mt 6 _zm3t end _x3mb 6 _zm3b nseg=_mnseg3 sel cable id=3 begin _x3mt 8 _zm3t end _x3mb 8 _zm3b nseg=_mnseg3 sel cable id=3 begin _x3mt 10 _zm3t end _x3mb 10 _zm

20、3b nseg=_mnseg3 ;定義參數(shù) sel cable id 3 prop xcarea 6.2e-4 emod 200e9 ytens 5.4e5 ycomp 1.0e5 sel cable id 3 prop gr_k 2.0e7 gr_coh 1.0e6 gr_per 0.377 gr_fric 25.0 plo con zdisp solve save kaiwa3.sav ;直墻段開挖支護 rest kaiwa3.sav ;第四步開挖 model null ran group 開挖4 ;定義錨索參數(shù) def inip x11= 100 x12= 100-9*cos(15*de

21、grad x13= 100-25*cos(15*degrad x21= 100 x22= 100-7*cos(15*degrad x23= 100-28*cos(15*degrad x31= 100 x32= 100-5*cos(15*degrad x33= 100-25*cos(15*degrad z11= 43.5 z21= 39 z31= 35 z12= 43.5-9*sin(15*degrad z13= 43.5-25*sin(15*degrad z22= 39-7*sin(15*degrad z23= 39-28*sin(15*degrad z32= 35-5*sin(15*degr

22、ad z33= 35-25*sin(15*degrad xt=100 ;樁的參數(shù) zt=43.8 ;樁的參數(shù) xb=100 ;樁的參數(shù) zb=21.3 ;樁的參數(shù) pnseg=20 ;樁的參數(shù) end inip ;加第一道錨索 def ins_cab loop j(1,7 yj=1.5*(j-1 i=7+j j1=j command sel cable id=j1 begin x11 yj z11 end x12 yj z12 nseg 9 ;錨索自由段 sel cable id=j1 begin x12 yj z12 end x13 yj z13 nseg 16 ;錨索錨固段 end_com

23、mand p_se=findselbycen(x11,yj,z11 mybeg1=sel_cid(p_se myend1=mybeg1+8 mybeg2=mybeg1+9 myend2=mybeg1+24 p_sn1=nd_near(x11,yj,z11 ;findsnbypos(x11,yj,z11 myid=nd_id(p_sn1 p_sl1=sn_link(p_sn1 mylink1=sl_cid(p_sl1 command sel cable id=j1 prop emod=200e9 yTens=190e7 xcarea=5.44e-4 & gr_k=1 gr_coh=1 g

24、r_per=0.471 gr_fric 1 ycomp 95e7 range cid mybeg1 myend1 ;自由段參數(shù) sel cable id=j1 prop emod=200e9 yTens=190e7 xcarea= 5.44e-4 & gr_k=5.35e9 gr_coh=4.2e5 gr_per=0.471 gr_fric 35 ycomp 95e7 range cid mybeg2 myend2 ;錨固段參數(shù) ; sel delete link range id mylink1 sel link id=mylink1 myid target zone sel lin

25、k attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id mylink1 sel cable id=j1 pretension 150e3 range cid mybeg1 myend1 end_command command sel pile id=i begin xt yj zt end xb yj zb nseg=22 ;灌注樁 sel pile prop emod=36.6e9 nu=0.2 xcarea=0.503 xciz=2.01e-2 xciy=2.01e-2 x

26、cj=4.02e-2 density=2.5e3 perim=2.513 & cs_scoh=5.027e5 cs_sfric=0.0 cs_sk=4e8 cs_ncoh=0.0 cs_nfric=0.0 cs_ngap=off cs_nk=4e8 slide=on end_command end_loop end ins_cab ;建立樁錨之間的連接 def merge_link1 node_num=0 node_pnt1 = nd_head loop while node_pnt1 # null ;尋找總節(jié)點數(shù)，注：不能自己任生成node，程序缺省的方式為連續(xù)生成無不連續(xù) node

27、_num = node_num+1 node_pnt1 = nd_next(node_pnt1 endloop node_num_minus1 = node_num-1 link_id=30000 loop ii (1,node_num_minus1 node_pnt1 = nd_find(ii xx1 = nd_pos(node_pnt1,2,1 yy1 = nd_pos(node_pnt1,2,2 zz1 = nd_pos(node_pnt1,2,3 ii_plus1 = ii+1 loop jj (ii_plus1,node_num node_pnt2 = nd_find(jj xx2

28、= nd_pos(node_pnt2,2,1 yy2 = nd_pos(node_pnt2,2,2 zz2 = nd_pos(node_pnt2,2,3 node_dist = sqrt(xx1-xx22+(yy1-yy22+(zz1-zz22 dist_tol = 1e-2 if node_dist < dist_tol then link_pnt1 = nd_link(node_pnt1 link_pnt2 = nd_link(node_pnt2 if link_pnt1 # null then temp1 = lk_delete(link_pnt1 endif if link_pn

29、t2 # null then temp2 = lk_delete(link_pnt2 endif link_id = link_id+1 command ;生成新link(6自由度全固結 ，大的node的id作為target node，小的node的id作為source node，需注意不同情況下的靈活調整 sel set link node_tol=dist_tol sel link id=link_id ii target = node tgt_num = jj ;指定link的ID ;sel link ii target = node tgt_num = jj ;不指定link的id，自

30、動生成 sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id=link_id endcommand endif endloop endloop end merge_link1 ;plo sk ;plo sel geom black black node off shrinkfac 0 plo con zdisp solve save kaiwa4.sav ;第五步開挖 rest kaiwa4.sav model null ran group 開挖5 ;加第二道錨

31、索 def ins_cab2 loop j2(1,7 yj2=1.5*(j2-1 jj2=j2+100 command sel cable id=jj2 begin x21 yj2 z21 end x22 yj2 z22 nseg 7 ;錨索自由段 sel cable id=jj2 begin x22 yj2 z22 end x23 yj2 z23 nseg 21 ;錨索錨固段 end_command p_se=findselbycen(x21,yj2,z21 mybeg21=sel_cid(p_se myend21=mybeg21+6 mybeg22=mybeg21+7 myend22=my

32、beg21+26 p_sn21=nd_near(x21,yj2,z21 ;findsnbypos(x21,yj2,z21 myid 2=nd_id(p_sn21 p_s21=sn_link(p_sn21 mylink21=sl_cid(p_s21 command sel cable id=jj2 prop emod=200e9 yTens=190e7 xcarea=5.44e-4 & gr_k=1 gr_coh=1 gr_per=0.471 gr_fric 1 ycomp 95e7 range cid mybeg21 myend21 ;自由段參數(shù) sel cable id=jj2 pr

33、op emod=200e9 yTens=190e7 xcarea=5.44e-4 & gr_k=5.35e9 gr_coh=4.2e5 gr_per=0.471 gr_fric 35 ycomp 95e7 range cid mybeg22 myend22 ;錨固段參數(shù) ; sel delete link range id mylink21 sel link id=mylink21 myid2 target zone sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid

34、range id mylink21 sel cable id=jj2 pretension 200e3 range cid mybeg21 myend21 end_command end_loop end ins_cab2 ;建立樁錨之間的連接 def merge_link12 node_num2=0 node_pnt21 = nd_head loop while node_pnt21 # null ;尋找總節(jié)點數(shù)，注：不能自己任生成node，程序缺省的方式為連續(xù)生成無不連續(xù) node_num2 = node_num2+1 node_pnt21 = nd_next(node_pnt21 end

35、loop node_num_minus21 = node_num2-1 link_id2=40000 loop ii2 (1,node_num_minus21 node_pnt21 = nd_find(ii2 xx21 = nd_pos(node_pnt21,2,1 yy21 = nd_pos(node_pnt21,2,2 zz21 = nd_pos(node_pnt21,2,3 ii_plus21 = ii2+1 loop jj2 (ii_plus21,node_num2 node_pnt22 = nd_find(jj2 xx22 = nd_pos(node_pnt22,2,1 yy22 =

36、 nd_pos(node_pnt22,2,2 zz22 = nd_pos(node_pnt22,2,3 node_dist2 = sqrt(xx21-xx222+(yy21-yy222+(zz21-zz222 dist_tol2 = 1e-2 if node_dist2 < dist_tol2 then link_pnt21 = nd_link(node_pnt21 link_pnt22 = nd_link(node_pnt22 if link_pnt21 # null then temp21 = lk_delete(link_pnt21 endif if link_pnt22 # nu

37、ll then temp22 = lk_delete(link_pnt22 endif link_id2 = link_id2+1 command ;生成新link(6自由度全固結 ，大的node的id作為target node，小的node的id作為source node，需注意不同情況下的靈活調整 sel set link node_tol2=dist_tol2 sel link id=link_id2 ii2 target = node tgt_num = jj2 ;指定link的ID ;sel link ii target = node tgt_num = jj ;不指定link的id

38、，自動生成 sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id=link_id2 endcommand endif endloop endloop end merge_link12 plo con zdisp solve save kaiwa5.sav ;第六步開挖 rest kaiwa5.sav model null ran group 開挖6 ;加第三道錨索 def ins_cab3 loop j3(1,7 yj3=1.5*(j3-1 jj3=j3+200

39、 command sel cable id=jj3 begin x31 yj3 z31 end x32 yj3 z32 nseg 5 ;錨索自由段 sel cable id=jj3 begin x32 yj3 z32 end x33 yj3 z33 nseg 20 ;錨索錨固段 end_command p_se=findselbycen(x31,yj3,z31 mybeg31=sel_cid(p_se myend31=mybeg31+4 mybeg32=mybeg31+5 myend32=mybeg31+24 p_sn31=nd_near(x31,yj3,z31 ;findsnbypos(x3

40、1,yj3,z31 myid3=nd_id(p_sn31 p_s31=sn_link(p_sn31 mylink31=sl_cid(p_s31 command sel cable id=jj3 prop emod=200e9 yTens=190e7 xcarea=5.44e-4 & gr_k=1 gr_coh=1 gr_per=0.471 gr_fric 1 ycomp 95e7 range cid mybeg31 myend31 ;自由段參數(shù) sel cable id=jj3 prop emod=200e9 yTens=190e7 xcarea=5.44e-4 & gr_k=5.35e9 gr_coh=4.2e5 gr_per=0.471 gr_fric 35 ycomp 95e7 range cid mybeg32 myend32 ;錨固段參數(shù) ; sel delete link range id mylink31 sel link id=mylink31 myid3 target zone sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id mylink31 sel cable id=