临时支撑CD法

1、!修改网格及参数（级围岩）finish/clear*afun,degx_zs=-50y_zs=30x_yx=50y_yx=-40num=24th_s=0th_f=180-th_sth=(th_f-th_s)/numyx_1_x=0yx_1_y=0yx_2_x=0yx_2_y=11.37yx_3_x=3.34yx_3_y=0yx_4_x=-3.34yx_4_y=0r_1=5.8r_2=14.4r_3=2.67r_4=2.67r_jg_1=8.8r_jg_2=17.4r_jg_3=5.67r_jg_4=5.67!beam3单元厚度（考虑钢纤维、二衬提高0.5）t1=0.25*1.5!link1单元

2、半径r1=12.5e-3!plane42单元e_r=0.2e9u_r=0.35dens_r=1950c_r=0.18e6fi_r=24!beam3单元e_c=23e9u_c=0.2dens_c=2200!link1单元e_a=210e9u_a=0.25dens_a=7700/prep7!单元类型et,1,beam3et,2,plane42,2keyopt,2,3,2et,3,link1!实常数r,1,t1,t1*t1*t1/12,t1r,2,3.14*r1*r1!材料属性!左上支护mp,ex,1,e_cmp,prxy,1,u_cmp,dens,1,dens_c!右上支护mp,ex,2,e_cm

3、p,prxy,2,u_cmp,dens,2,dens_c!左下支护mp,ex,3,e_cmp,prxy,3,u_cmp,dens,3,dens_c!右下支护mp,ex,4,e_cmp,prxy,4,u_cmp,dens,4,dens_c!左上开挖土体mp,ex,5,e_rmp,prxy,5,u_rmp,dens,5,dens_rtb,dp,5tbdata,1,c_r,fi_r!右上开挖土体mp,ex,6,e_rmp,prxy,6,u_rmp,dens,6,dens_rtb,dp,6tbdata,1,c_r,fi_r!左下开挖土体mp,ex,7,e_rmp,prxy,7,u_rmp,dens,7

4、,dens_rtb,dp,7tbdata,1,c_r,fi_r!右下开挖土体mp,ex,8,e_rmp,prxy,8,u_rmp,dens,8,dens_rtb,dp,8tbdata,1,c_r,fi_r!左上锚杆mp,ex,9,e_amp,prxy,9,u_amp,dens,9,dens_a!右上锚杆mp,ex,10,e_amp,prxy,10,u_amp,dens,10,dens_a!左临时支撑mp,ex,11,e_cmp,prxy,11,u_cmp,dens,11,dens_c!上临时支撑mp,ex,12,e_cmp,prxy,12,u_cmp,dens,12,dens_c!右临时支撑m

5、p,ex,13,e_cmp,prxy,13,u_cmp,dens,13,dens_c!下临时支撑mp,ex,14,e_cmp,prxy,14,u_cmp,dens,14,dens_c!围岩mp,ex,15,e_rmp,prxy,15,u_rmp,dens,15,dens_rtb,dp,15tbdata,1,c_r,fi_r!实体建模k,1,yx_1_x,yx_1_yk,2,yx_2_x,yx_2_yk,3,yx_3_x,yx_3_yk,4,yx_4_x,yx_4_yk,101,yx_1_x,yx_1_y,10k,102,yx_2_x,yx_2_y,10k,103,yx_3_x,yx_3_y,1

6、0k,104,yx_4_x,yx_4_y,10circle,1,r_1,101circle,2,r_2,102circle,3,r_3,103circle,4,r_4,104lcsl,alllsel,s,17,22,1lsel,a,46,48,1lsel,a,27,28,1lsel,a,43lsel,inveldel,alllsel,alllcomb,17,18lcomb,19,21lcomb,43,27lcomb,28,46cm,zh,linecircle,1,r_jg_1,101circle,2,r_jg_2,102circle,3,r_jg_3,103circle,4,r_jg_4,10

7、4cmsel,s,zh,linelsel,invelcsl,alllsel,s,18,21,3lsel,a,23,26,1lsel,a,33,34,1lsel,a,58,59,1lsel,a,55,62,7cmsel,a,zh,linelsel,inveldele,all,1lsel,alllcomb,18,21lcomb,23,25lcomb,55,33lcomb,34,58allselkdele,101,104,1rectng,x_zs,x_yx,y_zs,y_yxasbl,all,all!利用工作平面切割围岩allselasel,s,3,4,1wpcsyswpoff,-1.8*r_jg_

8、1wprota,-90asbw,allwpcsyswprota,-90asbw,allwpcsyswpoff,1.8*r_jg_1wprota,-90asbw,allwpcsyswpoff,1.7*r_jg_1wprota,-90asbw,allwpcsyswprota,-90asbw,allwpcsyswpoff,-1.4*r_jg_1wprota,-90asbw,allwpcsysallsel/pnum,line,1/pnum,area,1!锚杆线allselasel,s,5,6,1wprota,-90wprota,-th_s*do,i,1,num,1 wprota,-th asbw,al

9、l*enddoallsel!临时支撑allsellarc,81,12,7,1.4*r_jg_1l,7,5lsel,s,19,23,4asbl,2,all!围岩分块，利于划分网格allsellsel,nonel,42,92asbl,19,alllsel,nonel,44,70asbl,18,alllsel,nonel,55,53asbl,25,alllsel,nonel,57,54asbl,24,all!赋单元属性及划分网格allselaplotlesize,131,10lesize,132,11lesize,129,10lesize,130,6!左上土体!aclear,46allselasel

10、,s,46aplotlsla,scsys,1lsel,r,loc,x,r_1lplotlesize,all,1cm,zh_zs,linelsel,r,loc,y,13*th,18*thlccat,allcmsel,s,zh_zs,linelsel,r,loc,y,18*th,24*thlccat,allcsys,0mat,5type,2mshape,0mshkey,1amesh,alllsel,s,lccatldele,all!右上土体allsel!aclear,47asel,s,47aplotlsla,scsys,1lsel,r,loc,x,r_1lplotlesize,all,1cm,zh

11、_ys,linelsel,r,loc,y,0*th,7*thlccat,allcmsel,s,zh_ys,linelsel,r,loc,y,7*th,13*thlccat,allcsys,0mat,6type,2mshape,0mshkey,1amesh,alllsel,s,lccatldele,all!左下土体!alcear,45allsellesize,27,8lesize,48,2lesize,22,2lccat,22,48mat,7type,2mshape,0mshkey,1amesh,45lsel,s,lccatldele,all!右下土体!alcear,6allsellesize,

12、28,8lesize,47,2lesize,20,2lccat,20,47mat,8type,2mshape,0mshkey,1amesh,6lsel,s,lccatldele,all!临时支撑allselmat,11real,1type,1lmesh,131mat,12lmesh,132mat,13lmesh,129mat,14lmesh,130!支护allselcmsel,s,zh_zslatt,1,1,1lmesh,allcmsel,s,zh_yslatt,2,1,1lmesh,allasel,s,45lsla,slsel,u,130,131,1cm,zh_zx,linelatt,3,1

13、,1lmesh,allasel,s,6lsla,slsel,u,129,130,1cm,zh_yx,linelatt,4,1,1lmesh,all!左上锚杆和锚区围岩csys,1asel,s,loc,x,r_1,r_jg_1asel,r,loc,y,13*th,24*thlsla,slsel,r,loc,x,(r_1+r_jg_1)/2csys,0cm,mg_zs,linemat,9real,2type,3lesize,all,4,1.8lmesh,allmat,15type,2mshape,0mshkey,1amesh,all!右上锚杆和锚区围岩csys,1asel,s,loc,x,r_1,

14、r_jg_1asel,r,loc,y,0*th,13*thlsla,slsel,r,loc,x,(r_1+r_jg_1)/2csys,0cm,mg_ys,linemat,10real,2type,3lesize,all,4,1.6lmesh,allmat,15type,2mshape,0mshkey,1amesh,all!aclear,22,23,1allsellesize,15,4,1.6lcomb,26,62lesize,26,6lesize,33,6amap,23,12,34,27,7allsellcomb,24,59lesize,24,6lesize,34,6amap,22,12,34

15、,21,5!周围围岩!aclear,allallselasel,s,loc,x,-1.8*r_jg_1,1.8*r_jg_1asel,r,loc,y,-1.4*r_jg_1,1.7*r_jg_1lsel,s,extlesize,all,6lsel,s,42,65,23lsel,a,19,23,4lsel,a,52,53,1lsel,a,133,134,1lesize,all,9,1.2,1asll,scsys,1lsel,s,loc,x,r_jg_1lsel,r,loc,y,0*th,6*thlccat,alllsel,s,loc,x,r_jg_1lsel,r,loc,y,6*th,12*th

16、lccat,alllsel,s,loc,x,r_jg_1lsel,r,loc,y,12*th,18*thlccat,alllsel,s,loc,x,r_jg_1lsel,r,loc,y,18*th,24*thlccat,allcsys,0mat,15type,2mshape,0mshkey,1amesh,alllsel,s,lccatldele,all!其他围岩!aclear,allasel,s,loc,x,-100*r_jg_1,-1.8*r_jg_1lsla,slsel,r,tan1,y,0lesize,all,16,0.5,1asel,s,loc,x,1.8*r_jg_1,100*r_j

17、g_1lsla,slsel,r,tan1,y,0lesize,all,16,0.5,1asel,s,loc,y,1.7*r_jg_1,100*r_jg_1lsla,alllsel,r,tan1,x,0lesize,all,8,0.5,1asel,s,loc,y,-100*r_jg_1,-1.4*r_jg_1lsla,alllsel,r,tan1,x,0lesize,all,13,0.5,1allselasel,s,loc,x,-1.8*r_jg_1,1.8*r_jg_1asel,r,loc,y,-1.4*r_jg_1,1.7*r_jg_1asel,invemat,15type,2mshape,

18、0mshkey,1amesh,allallselnummrg,nodenummrg,elemnummrg,kpnumcmp,all!施加约束allselnsel,s,loc,x,x_zsnsel,a,loc,x,x_yxd,all,uxnsel,s,loc,y,y_yxd,all,uy!指定求解选项/solunlgeom,onnropt,fullpred,onlnsrch,onnsubst,10,50,10!荷载步1,自重应力场allselesel,s,type,1,3,2ekill,allallselacel,9.8allselsolvesave,step1,dbesel,s,liveepl

19、ot!荷载步2,左上开挖!提取节点力/post1esel,s,mat,5nsle,snsel,r,ext*get,n,node,count*dim,n_nf,array,n,3*do,i,1,n*get,n_nf(i,1),node,num,minnsel,u,n_nf(i,1)*enddo!*status,n_nf*do,i,1,nnsel,s,n_nf(i,1)fsum*get,n_nf(i,2),fsum,item,fx*get,n_nf(i,3),fsum,item,fy*enddo!*status,n_nfparsav,all,002,PARM/soluantype,rest!杀死左

20、上土体单元allselesel,s,mat,5ekill,all!施加等效释放节点力parres,002,PARMesel,s,mat,5nsle,snsel,r,ext*do,i,1,nf,n_nf(i,1),fx,0.65*n_nf(i,2)f,n_nf(i,1),fy,0.65*n_nf(i,3)*enddo!约束死单元自由度allselesel,s,livensle,snsel,inved,all,allallselsolvesave,step2,dbesel,s,liveeplot!荷载步3!删除节点力，激活支护、锚杆、临时支撑allselesel,s,mat,5nsle,snsel

21、,r,extfdele,all,allesel,s,mat,1,9,8esel,a,mat,11,12,1ealive,allallselsolvesave,step3,db!荷载步4,左下开挖，应力释放!提取节点力/post1esel,s,mat,7nsle,snsel,r,ext*get,n,node,count*dim,n_nf,array,n,3*do,i,1,n*get,n_nf(i,1),node,num,minnsel,u,n_nf(i,1)*enddo!*status,n_nf*do,i,1,nnsel,s,n_nf(i,1)fsum*get,n_nf(i,2),fsum,it

22、em,fx*get,n_nf(i,3),fsum,item,fy*enddo!*status,n_nfparsav,all,002,PARM/soluantype,rest!杀死左下土体单元allselealive,allesel,s,mat,2,5,1esel,a,mat,10esel,a,mat,13,14,1esel,a,mat,7ekill,all!施加等效释放节点力parres,002,PARMesel,s,mat,7nsle,snsel,r,ext*do,i,1,nf,n_nf(i,1),fx,0.65*n_nf(i,2)f,n_nf(i,1),fy,0.65*n_nf(i,3)*

23、enddo!约束死单元自由度allselesel,s,livensle,snsel,inved,all,allallselsolvesave,step4,dbesel,s,liveeplot!荷载步5!删除节点力，激活支护、临时支撑esel,s,mat,7nsle,snsel,r,extfdele,all,allesel,s,mat,3,14,11ealive,allallselsolvesave,step5,db!荷载步6,右上开挖，应力释放!提取节点力/post1esel,s,mat,6nsle,snsel,r,ext*get,n,node,count*dim,n_nf,array,n,3

24、*do,i,1,n*get,n_nf(i,1),node,num,minnsel,u,n_nf(i,1)*enddo!*status,n_nf*do,i,1,nnsel,s,n_nf(i,1)fsum*get,n_nf(i,2),fsum,item,fx*get,n_nf(i,3),fsum,item,fy*enddo!*status,n_nfparsav,all,002,PARM/soluantype,rest!杀死右上土体单元allselealive,allesel,s,mat,2esel,a,mat,4,7,1esel,a,mat,10,13,3ekill,all!施加等效释放节点力pa

25、rres,002,PARMesel,s,mat,6nsle,snsel,r,ext*do,i,1,nf,n_nf(i,1),fx,0.65*n_nf(i,2)f,n_nf(i,1),fy,0.65*n_nf(i,3)*enddo!约束死单元自由度allselesel,s,livensle,snsel,inved,all,allallselsolvesave,step6,dbesel,s,liveeplot!荷载步7!删除节点力，激活支护、锚杆、临时支撑allselesel,s,mat,6nsle,snsel,r,extfdele,all,allesel,s,mat,2,10,8esel,a,m

26、at,13ealive,allallselsolvesave,step7,db!荷载步8,右下开挖，应力释放!提取节点力/post1esel,s,mat,8nsle,snsel,r,ext*get,n,node,count*dim,n_nf,array,n,3*do,i,1,n*get,n_nf(i,1),node,num,minnsel,u,n_nf(i,1)*enddo!*status,n_nf*do,i,1,nnsel,s,n_nf(i,1)fsum*get,n_nf(i,2),fsum,item,fx*get,n_nf(i,3),fsum,item,fy*enddo!*status,n

27、_nfparsav,all,002,PARM/soluantype,rest!杀死右下土体单元allselealive,allesel,s,mat,4,7,1esel,a,mat,8ekill,all!施加等效释放节点力parres,002,PARMesel,s,mat,8nsle,snsel,r,ext*do,i,1,nf,n_nf(i,1),fx,0.65*n_nf(i,2)f,n_nf(i,1),fy,0.65*n_nf(i,3)*enddo!约束死单元自由度allselesel,s,livensle,snsel,inved,all,allallselsolvesave,step8,dbesel,s,liveeplot!荷载步9!删除节点力，激活支护esel,s,mat,8nsle,snsel,r,e