ANSYS实例分析报告连拱隧道开挖命令流.docx

1、ANSYS实例分析报告连拱隧道开挖命令流!III级围岩，二维，全断面开挖/cle/filname,ch-three-2D,1!几何参数!计算外围区域*set,x1,-60 *set,x2,60*set,y1,-30*set,y2,50*set,l1,3.5 !锚杆长度*set,r1,0.0125 !锚杆半径*set,th1,0.15 !初次衬砌厚度*set,th2,0.05 !中导洞喷层厚!材料参数!中隔墙*set,ec25,29.5e9*set,uc25,0.2*set,dens25,2500!喷射混凝土*set,ec20,21e9*set,uc20,0.2*set,dens20,2200!

2、锚杆*set,el,170e9*set,ul,0.3*set,densl,7800!围岩*set,er,6e9*set,ur,0.28*set,densr,2500*set,cr,0.8e6*set,fir,40/prep7!定义单元类型et,1,plane42keyopt,1,3,2 !平面应变问题et,2,beam3keyopt,2,6,1 !打开输出内力选项et,3,link1!定义实常数r,1,th1,th1/12,th1r,2,th2,th2/12,th2r,3,3.1415926*r1*2!几何建模 !建中隔墙k,1,-3,-1.69k,2,-3,-0.89 k,3,-3,-0.3

3、7 k,4,-1.33552,-0.37 k,5,-0.95134,1.57593k,6,-1.25121,3.53662k,7,1.25121,3.53662 k,8,0.95134,1.57593k,9,1.33552,-0.37k,10,3,-0.37k,11,3,-0.89 k,12,3,-1.69l,1,2l,2,3l,3,4larc,4,6,5 larc,6,7,5,3larc,7,9,8l,9,10l,10,11l,11,12l,12,1!左侧围岩及锚杆圈 k,13,-3,0.81k,14,-2.52524,2.42963 k,15,-1.483023,4.11778k,16,-

4、1.77650,4.67036k,17,-2.058364,5.09399k,18,-2.37658,5.49104k,19,-2.734606,5.864085k,20,-3.12567,6.20235k,21,-5.88951,7.72342 k,22,-9,8.25k,23,-12.11049,7.72342 k,24,-14.8743,6.20235 k,25,-16.57763,3.97894 k,26,-17.02812,1.21456 k,27,-16.9,-0.89k,28,-20.4,-0.89k,29,-20.51567,0.91966 k,30,-19.7906,5.36

5、705 k,31,-17.05,8.944 k,32,-13.26253,11.02838 k,33,-9,11.75k,34,-4.146154,10.80594k,35,0,8.11142 k,36,1.143324,6.85081!导洞左侧圈l,3,13larc,13,6,14!左侧开挖圈larc,6,16,15larc,16,18,17 larc,18,20,19 larc,20,22,21larc,22,24,23larc,24,26,25l,26,27!锚杆圈l,28,29larc,29,31,30larc,31,33,32 larc,33,35,34 l,35,36!内圈线合并l

6、comb,15,16lcomb,15,17lcomb,15,18lcomb,15,19!外圈线合并lcomb,20,21lcomb,20,22lcomb,20,23numcmp,allldiv,14,2ldiv,15,36ldiv,16,36ldiv,17,2kdele,19,22,1kdele,25,27,1numcmp,allkdele,17numcmp,all!锚杆线（40根，含建3维时可用的另外20根锚杆）l,16,21l,22,93l,17,20*do,i,23,57,1l,i,115-i*enddol,18,19l,18,2numcmp,all!内层细化区域k,100,-16.9,

7、-1.69k,101,-20.4,-1.69k,102,-24.9,-1.69k,103,-24.9,-0.89k,104,-24.9,16.25625k,105,16.25625k,106,x1,-1.69k,107,x1,16.25625k,108,-24.9,y2k,109,-24.9,y1l,1,100l,100,101l,101,102l,102,103l,103,104l,104,105l,19,103l,70,104l,18,100l,19,101l,102,109l,102,106l,104,107l,104,108!左侧线镜像至右并建立各面allslsymm,x,11,143

8、,1rectng,x1,x2,y1,y2asbl,1,allnummrg,all,lownumcmp,alllcomb,287,288lcomb,289,290lcomb,294,295lcomb,296,297lcomb,130,292numcmp,allalls!定义材料属性!岩体*do,i,1,6,1mp,ex,i,ermp,prxy,i,urmp,dens,i,densrtb,dp,itbdata,1,cr,fir*enddo!中导洞及隧道初支*do,i,7,10,1mp,ex,i,ec20mp,prxy,i,uc20mp,dens,i,dens20*enddo!锚杆左右侧*do,i,

9、11,12,1mp,ex,i,elmp,prxy,i,ulmp,dens,i,densl*enddo!中隔墙mp,ex,13,ec25mp,prxy,13,uc25mp,dens,13,dens25!中导洞初支!左侧type,2mat,7real,2lsel,s,line,1,2,1lsel,a,line,11,12,1lesize,1,1lesize,2,1lesize,11,2lesize,12,6lmesh,all!右侧mat,8real,2lsel,s,line,8,9,1lsel,a,line,259,260,1lsel,a,line,5lesize,9,1lesize,8,1les

10、ize,260,2lesize,259,6lesize,5,4lmesh,all!左侧隧道初支mat,9real,1lsel,s,line,13,15,1lsel,a,line,18,52,1lsel,a,line,288lesize,13,3lesize,all,1lmesh,all!右侧隧道初支mat,10real,1lsel,s,line,139,140,1lsel,a,line,142,176,1lsel,a,line,286lsel,a,line,292lesize,286,3lesize,all,1lmesh,all!左隧道锚杆type,3mat,11real,3lsel,s,li

11、ne,89,123,2lsel,a,line,263,269,3lsel,a,line,271,279,8lsel,a,line,210lsel,a,line,266lesize,263,3lesize,269,3lesize,271,1lesize,279,1lesize,210,1lesize,266,1lesize,all,5,2lmesh,all!右隧道锚杆type,3mat,12real,3lsel,s,line,212,246,2lsel,a,line,277,280,3lsel,a,line,265lsel,a,line,270lsel,a,line,87lsel,a,line,

12、17lesize,277,3lesize,280,3lesize,265,1lesize,270,1lesize,87,1lesize,17,1lesize,all,5,2lmesh,all!中墙type,1mat,1lsel,s,line,1,10,1lesize,3,2lesize,4,8lesize,6,8lesize,7,2lesize,10,8amesh,97!左中导洞mat,2amesh,2!右中导洞mat,3amesh,3!左隧道内岩体mat,4lsel,s,line,124lesize,124,20amesh,101!右隧道内岩体mat,5lsel,s,line,261lesi

13、ze,261,20amesh,102!导洞上方围岩mat,6lsel,s,line,264,267,3lsel,a,line,281lsel,a,line,268,278,10lsel,a,line,272lesize,264,3lesize,278,3lesize,all,1!左侧锚杆加固圈lsel,s,line,88,122,2 !还有264,267,281为扩展到3维的锚杆lesize,all,5,2lsel,s,line,53,86,1lsel,a,line,16lsel,a,line,287lesize,all,1!右侧锚杆加固圈lsel,s,line,211,245,2 !还有27

14、8,87,17为扩展到3维的锚杆lesize,all,5,2lsel,s,line,177,209,1lsel,a,line,291,294,3lsel,a,line,141lesize,all,1!小外围上部lsel,s,line,64,86,1lsel,a,line,188,209,1lsel,a,line,294lccat,alllsel,a,line,130,132,1lsel,a,line,252lesize,130,46lesize,132,8,2lesize,252,8,2amesh,103!小外围左lsel,s,line,53,63,1lsel,a,line,287lsel,a

15、,line,16lccat,alllsel,s,line,129lsel,a,line,131lesize,129,13lesize,131,8,2amesh,99!小外围右lsel,s,line,291lsel,a,line,177,187,1lsel,a,line,141lccat,alllsel,a,line,250lesize,250,13lsel,a,line,251lesize,251,8,2amesh,100!小外围下lsel,s,line,128lesize,128,1lsel,s,line,127lesize,127,8,2lsel,s,line,133,134,1lesiz

16、e,all,1lsel,s,line,126lesize,126,5,2lsel,s,line,125lesize,125,20lsel,s,line,262lesize,262,20lsel,s,line,253,254,1lesize,all,1lsel,s,line,247lesize,247,5,2lsel,s,line,249lesize,249,1lsel,s,line,248lesize,248,8,2amesh,42,43amesh,1,4,3amesh,82,83!大外围网格划分lsel,s,line,274,276,2lsel,a,line,136,137,1lsel,a,

17、line,256,257,1lsel,a,line,282,285,3lesize,all,18,2.5lsel,s,line,273,283,10lsel,a,line,135,255,120lesize,all,15,2.5lsel,s,line,275,284,9lsel,a,line,138,258,120lesize,all,18,2.5lsel,s,line,290,295,5lesize,290,46lesize,295,74lsel,s,line,289,293lesize,all,14lsel,s,line,128,129lccat,alllsel,s,line,249,25

18、0lccat,alllsel,s,line,125,127,1lsel,a,line,10lsel,a,line,262lsel,a,line,247,248,1lccat,allamesh,allallslsel,s,lccatldele,allnumcmp,allallselsave,data!施加约束allsnsel,s,loc,x,x1nsel,a,loc,x,x2d,all,uxnsel,s,loc,y,y1d,all,uyallsacel,10fini/soluantype,staticpred,onlnsrch,onnlgeom,onnropt,fulloutress,all,a

19、lldeltim,0.1,0.05,0.2 !荷载步1time,1esel,s,type,2,3,1ekill,allesel,allesel,s,live nsle,s nsel,inve d,all,all,0 nsel,allesel,all! iswrite,1 !生成初应力文件solvesave,s1,db!荷载步2time,2esel,s,mat,2,3,1ekill,allesel,s,mat,7,8,1ealive,allnsle,sddele,all,allesel,allesel,s,livensle,snsel,inved,all,allesel,s,mat,1mpchg

20、,13,allallssolvesave,s2,db!荷载步3time,3esel,s,mat,4,7,3ekill,allesel,s,mat,9,11,2esel,a,mat,3ealive,allnsle,sddele,all,allesel,allesel,s,livensle,snsel,inved,all,allallssolvesave,s3,db!荷载步4time,4esel,s,mat,5,8,3esel,a,mat,3ekill,allesel,s,mat,10,12,2ealive,allnsle,sddele,all,allesel,allesel,s,livensle

21、,snsel,inved,all,allallssolvesave,s4,dbfini/post1!载荷步1结果resume,s1,dbfile,ch-three-2D,rst/pbc,all,0 !不显示边界约束等条件set,1,lastesel,s,liveplnsol,u,sum,0,1.0plnsol,u,xplnsol,u,yplnsol,s,1plnsol,s,3plnsol,s,eqvplnsol,eppl,eqv,0,1!载荷步2结果resume,s2,db/pbc,all,0 esel,s,liveesel,u,mat,13plnsol,u,sum,0,1 !未减去自重引起的

22、合位移云图plnsol,u,x !未减去自重引起的位移plnsol,u,y !未减去自重引起的位移plnsol,s,1plnsol,s,3plnsol,s,eqvplnsol,eppl,eqv,0,1esel,s,mat,13plnsol,u,yplnsol,s,1plnsol,s,3plnsol,s,eqv!中导洞初支内力esel,s,liveesel,r,ename,3 !选择beam3单元类型etable,fx_i,smisc,1etable,fx_j,smisc,7etable,fy_i,smisc,2etable,fy_j,smisc,8etable,mz_i,smisc,6etab

23、le,mz_j,smisc,12plls,mz_i,mz_j,-0.2 !显示弯矩plls,fx_i,fx_j,0.2 !显示轴力plls,fy_i,fy_j,0.2 !显示剪力!查看位移场file,ch-three-2D,rstset,2,lastesel,s,liveesel,u,mat,13lcdef,1,1 !定义第1个载荷步为工况1lcoper,sub,1 !当前载荷步的结果减去工况1的结果(位移结果有问题）esel,s,liveplnsol,u,sum !显示合位移云图plnsol,u,xplnsol,u,y!载荷步3结果resume,s3,db/pbc,all,0 esel,s,

24、liveesel,u,mat,13plnsol,s,eqvplnsol,s,1plnsol,s,3plnsol,eppl,eqv,0,1!中墙等效应力场esel,s,mat,13plnsol,s,eqvplnsol,s,1plnsol,s,3!中导洞初支及左隧道衬砌内力esel,s,liveesel,r,ename,3 !选择beam3单元类型etable,fx_i,smisc,1etable,fx_j,smisc,7etable,fy_i,smisc,2etable,fy_j,smisc,8etable,mz_i,smisc,6etable,mz_j,smisc,12plls,mz_i,mz

25、_j,-0.2 !显示弯矩（0.2为显示比例)plls,fx_i,fx_j,0.2 !显示轴力plls,fy_i,fy_j,0.2 !显示剪力!查看锚杆轴力esel,s,liveesel,r,ename,1etable,forx,smisc,1plls,forx,forx,0.05!查看位移场file,ch-three-2D,rstset,3,lastlcdef,1,1 !定义第1个载荷步为工况1lcoper,sub,1 !当前载荷步的结果减去工况1的结果esel,s,liveesel,u,mat,3plnsol,u,sum !显示合位移云图plnsol,u,xplnsol,u,y!载荷步4结

26、果resume,s4,db/pbc,all,0 esel,s,liveesel,u,mat,13plnsol,s,eqvplnsol,s,1plnsol,s,3plnsol,eppl,eqv,0,1!中墙等效应力场esel,s,mat,13plnsol,s,eqvplnsol,s,1plnsol,s,3plnsol,u,y!左及右隧道衬砌内力esel,s,liveesel,r,ename,3 !选择beam3单元类型etable,fx_i,smisc,1etable,fx_j,smisc,7etable,fy_i,smisc,2etable,fy_j,smisc,8etable,mz_i,sm

27、isc,6etable,mz_j,smisc,12plls,mz_i,mz_j,-0.2 !显示弯矩（0.2为显示比例)plls,fx_i,fx_j,0.2 !显示轴力plls,fy_i,fy_j,0.4 !显示剪力!查看锚杆轴力esel,allesel,r,ename,1etable,forx,smisc,1plls,forx,forx,0.05!查看位移场file,ch-three-2D,rstset,4,lastlcdef,1,1 !定义第1个载荷步为工况1lcoper,sub,1 !当前载荷步的结果减去工况1的结果esel,s,liveesel,u,mat,3plnsol,u,sum !显示合位移云图plnsol,u,xplnsol,u,y