ANSYS实例分析连拱隧道开挖命令流

!III级围岩，二维，全断面开挖/cle
/filname,ch-three-2D,1
!几何参数
!计算外围区域
*set,x1,-60
*set,x2,60
*set,y1,-30
*set,y2,50
*set,l1, !锚杆长度
*set,r1, !锚杆半径
*set,th1, !初次衬砌厚度
*set,th2, !中导洞喷层厚
!材料参数
!中隔墙
*set,ec25,
*set,uc25,
*set,dens25,2500
!喷射混凝土
*set,ec20,21e9 *set,uc20,
*set,dens20,2200
!锚杆
*set,el,170e9
*set,ul,
*set,densl,7800
!围岩
*set,er,6e9
*set,ur,
*set,densr,2500 *set,cr,
*set,fir,40
/prep7
!定义单元类型
et,1,plane42
keyopt,1,3,2 !平面应变问题et,2,beam3
keyopt,2,6,1 !打开输出内力选项et,3,link1
!定义实常数
r,1,th1,th1/12,th1
r,2,th2,th2/12,th2
r,3,*r1**2
!几何建模
!建中隔墙
k,1,-3,
k,2,-3,
k,3,-3,
k,4,,
k,5,,
k,6,,
k,7,,
k,8,,
k,10,3,
k,11,3,
k,12,3,
l,1,2
l,2,3
l,3,4
larc,4,6,5
larc,6,7,5,3
larc,7,9,8
l,9,10
l,10,11
l,11,12
l,12,1
!左侧围岩及锚杆圈k,13,-3,
k,14,,
k,15,,
k,16,,
k,17,,
k,19,,
k,20,,
k,21,,
k,22,-9,
k,23,,
k,24,,
k,25,,
k,26,,
k,27,,
k,28,,
k,29,,
k,30,,
k,31,,
k,32,,
k,33,-9,
k,34,,
k,35,0,
k,36,,
!导洞左侧圈
larc,13,6,14 !左侧开挖圈larc,6,16,15 larc,16,18,17 larc,18,20,19 larc,20,22,21 larc,22,24,23 larc,24,26,25 l,26,27
!锚杆圈
l,28,29 larc,29,31,30 larc,31,33,32 larc,33,35,34 l,35,36
!内圈线合并lcomb,15,16 lcomb,15,17 lcomb,15,18 lcomb,15,19
!外圈线合并
lcomb,20,21
lcomb,20,22
lcomb,20,23
numcmp,all
ldiv,14,,,2
ldiv,15,,,36
ldiv,16,,,36
ldiv,17,,,2
kdele,19,22,1
kdele,25,27,1
numcmp,all
kdele,17
numcmp,all
!锚杆线（40根，含建3维时可用的另外20根锚杆）l,16,21
l,22,93
l,17,20
*do,i,23,57,1
*enddo
l,18,19
l,18,2 numcmp,all !内层细化区域k,100,,
k,101,,
k,102,,
k,103,,
k,104,,
k,105,,
k,106,x1,
k,107,x1,
k,108,,y2
k,109,,y1
l,1,100
l,100,101
l,101,102
l,102,103
l,104,105
l,19,103
l,70,104
l,18,100
l,19,101
l,102,109
l,102,106
l,104,107
l,104,108
!左侧线镜像至右并建立各面alls
lsymm,x,11,143,1
rectng,x1,x2,y1,y2
asbl,1,all
nummrg,all,,,low numcmp,all
lcomb,287,288
lcomb,289,290
lcomb,294,295 lcomb,296,297 lcomb,130,292 numcmp,all
alls
!定义材料属性
!岩体
*do,i,1,6,1
mp,ex,i,er
mp,prxy,i,ur mp,dens,i,densr tb,dp,i
tbdata,1,cr,fir
*enddo
!中导洞及隧道初支*do,i,7,10,1
mp,ex,i,ec20
mp,prxy,i,uc20 mp,dens,i,dens20
*enddo
!锚杆左右侧
*do,i,11,12,1 mp,ex,i,el
mp,prxy,i,ul mp,dens,i,densl
*enddo
!中隔墙
mp,ex,13,ec25 mp,prxy,13,uc25 mp,dens,13,dens25
!中导洞初支
!左侧
type,2
mat,7
real,2
lsel,s,line,,1,2,1
lesize,1,,,1 lesize,2,,,1 lesize,11,,,2 lesize,12,,,6 lmesh,all
!右侧
mat,8
real,2
lsel,s,line,,8,9,1 lsel,a,line,,259,260,1 lsel,a,line,,5 lesize,9,,,1 lesize,8,,,1 lesize,260,,,2 lesize,259,,,6 lesize,5,,,4 lmesh,all
!左侧隧道初支mat,9
real,1
lsel,s,line,,13,15,1
lsel,a,line,,288 lesize,13,,,3 lesize,all,,,1 lmesh,all
!右侧隧道初支mat,10
real,1
lsel,s,line,,139,140,1 lsel,a,line,,142,176,1 lsel,a,line,,286 lsel,a,line,,292 lesize,286,,,3 lesize,all,,,1 lmesh,all
!左隧道锚杆type,3
mat,11
real,3
lsel,s,line,,89,123,2 lsel,a,line,,263,269,3
lsel,a,line,,210 lsel,a,line,,266 lesize,263,,,3 lesize,269,,,3 lesize,271,,,1 lesize,279,,,1 lesize,210,,,1 lesize,266,,,1 lesize,all,,,5,2 lmesh,all
!右隧道锚杆type,3
mat,12
real,3
lsel,s,line,,212,246,2 lsel,a,line,,277,280,3 lsel,a,line,,265 lsel,a,line,,270
lsel,a,line,,87 lsel,a,line,,17 lesize,277,,,3 lesize,280,,,3 lesize,265,,,1 lesize,270,,,1 lesize,87,,,1 lesize,17,,,1 lesize,all,,,5,2 lmesh,all
!中墙
type,1
mat,1
lsel,s,line,,1,10,1 lesize,3,,,2 lesize,4,,,8 lesize,6,,,8 lesize,7,,,2 lesize,10,,,8 amesh,97
!左中导洞mat,2 amesh,2
!右中导洞mat,3 amesh,3
!左隧道内岩体mat,4
lsel,s,line,,124 lesize,124,,,20 amesh,101
!右隧道内岩体mat,5
lsel,s,line,,261 lesize,261,,,20 amesh,102
!导洞上方围岩
mat,6
lsel,s,line,,264,267,3
lsel,a,line,,281
lsel,a,line,,268,278,10
lsel,a,line,,272
lesize,264,,,3
lesize,278,,,3
lesize,all,,,1
!左侧锚杆加固圈
lsel,s,line,,88,122,2 !还有264,267,281为扩展到3维的锚杆lesize,all,,,5,2
lsel,s,line,,53,86,1
lsel,a,line,,16
lsel,a,line,,287
lesize,all,,,1
!右侧锚杆加固圈
lsel,s,line,,211,245,2 !还有278,87,17为扩展到3维的锚杆
lesize,all,,,5,2 lsel,s,line,,177,209,1 lsel,a,line,,291,294,3 lsel,a,line,,141 lesize,all,,,1
!小外围上部
lsel,s,line,,64,86,1 lsel,a,line,,188,209,1 lsel,a,line,,294 lccat,all
lsel,a,line,,130,132,1 lsel,a,line,,252 lesize,130,,,46 lesize,132,,,8,2 lesize,252,,,8,2 amesh,103
!小外围左
lsel,s,line,,53,63,1 lsel,a,line,,287
lccat,all
lsel,s,line,,129 lsel,a,line,,131 lesize,129,,,13 lesize,131,,,8,2 amesh,99
!小外围右
lsel,s,line,,291 lsel,a,line,,177,187,1 lsel,a,line,,141 lccat,all
lsel,a,line,,250 lesize,250,,,13 lsel,a,line,,251 lesize,251,,,8,2 amesh,100
!小外围下
lsel,s,line,,128
lsel,s,line,,127 lesize,127,,,8,2
lsel,s,line,,133,134,1 lesize,all,,,1
lsel,s,line,,126 lesize,126,,,5,2
lsel,s,line,,125 lesize,125,,,20
lsel,s,line,,262 lesize,262,,,20
lsel,s,line,,253,254,1 lesize,all,,,1
lsel,s,line,,247 lesize,247,,,5,2
lsel,s,line,,249 lesize,249,,,1
lsel,s,line,,248 lesize,248,,,8,2 amesh,42,43 amesh,1,4,3 amesh,82,83
!大外围网格划分lsel,s,line,,274,276,2 lsel,a,line,,136,137,1 lsel,a,line,,256,257,1 lsel,a,line,,282,285,3 lesize,all,,,18,
lsel,s,line,,273,283,10 lsel,a,line,,135,255,120 lesize,all,,,15,
lsel,s,line,,275,284,9 lsel,a,line,,138,258,120 lesize,all,,,18,
lsel,s,line,,290,295,5
lesize,290,,,46 lesize,295,,,74
lsel,s,line,,289,293 lesize,all,,,14
lsel,s,line,,128,129 lccat,all
lsel,s,line,,249,250 lccat,all
lsel,s,line,,125,127,1 lsel,a,line,,10
lsel,a,line,,262 lsel,a,line,,247,248,1 lccat,all
amesh,all
alls
lsel,s,lccat
ldele,all numcmp,all
allsel
save,data
!施加约束alls
nsel,s,loc,x,x1 nsel,a,loc,x,x2 d,all,ux nsel,s,loc,y,y1 d,all,uy
alls acel,,10
fini
/solu antype,static pred,on lnsrch,on nlgeom,on nropt,full outress,all,all deltim,,,
!荷载步1
time,1
esel,s,type,,2,3,1
ekill,all
esel,all
esel,s,live
nsle,s
nsel,inve
d,all,all,0
nsel,all
esel,all
! iswrite,1 !生成初应力文件solve
save,s1,db
!荷载步2
time,2
esel,s,mat,,2,3,1
ekill,all
esel,s,mat,,7,8,1
ealive,all nsle,s
ddele,all,all esel,all
esel,s,live nsle,s
nsel,inve
d,all,all
esel,s,mat,,1 mpchg,13,all alls
solve
save,s2,db
!荷载步3 time,3
esel,s,mat,,4,7,3 ekill,all
esel,s,mat,,9,11,2 esel,a,mat,,3 ealive,all nsle,s
ddele,all,all esel,all
esel,s,live
nsle,s
nsel,inve
d,all,all
alls
solve
save,s3,db
!荷载步4
time,4
esel,s,mat,,5,8,3 esel,a,mat,,3 ekill,all
esel,s,mat,,10,12,2 ealive,all
nsle,s
ddele,all,all esel,all
esel,s,live
nsle,s
nsel,inve
d,all,all
alls
solve
save,s4,db
fini
/post1
!载荷步1结果
resume,s1,db
file,ch-three-2D,rst
/pbc,all,,0 !不显示边界约束等条件set,1,last
esel,s,live
plnsol,u,sum,0,
plnsol,u,x
plnsol,u,y
plnsol,s,1
plnsol,s,3
plnsol,s,eqv
plnsol,eppl,eqv,0,1
!载荷步2结果
resume,s2,db
/pbc,all,,0
esel,s,live
esel,u,mat,,13
plnsol,u,sum,0,1 !未减去自重引起的合位移云图plnsol,u,x !未减去自重引起的位移
plnsol,u,y !未减去自重引起的位移
plnsol,s,1
plnsol,s,3
plnsol,s,eqv
plnsol,eppl,eqv,0,1
esel,s,mat,,13
plnsol,u,y
plnsol,s,1
plnsol,s,3
plnsol,s,eqv
!中导洞初支内力
esel,s,live
esel,r,ename,,3 !选择beam3单元类型etable,fx_i,smisc,1
etable,fx_j,smisc,7
etable,fy_i,smisc,2
etable,fy_j,smisc,8
etable,mz_i,smisc,6
etable,mz_j,smisc,12
plls,mz_i,mz_j, !显示弯矩
plls,fx_i,fx_j, !显示轴力
plls,fy_i,fy_j, !显示剪力
!查看位移场
file,ch-three-2D,rst
set,2,last
esel,s,live
esel,u,mat,,13
lcdef,1,1 !定义第1个载荷步为工况1
lcoper,sub,1 !当前载荷步的结果减去工况1的结果(位移结果有问题）esel,s,live
plnsol,u,sum !显示合位移云图
plnsol,u,x
plnsol,u,y
!载荷步3结果
resume,s3,db
/pbc,all,,0
esel,s,live
esel,u,mat,,13
plnsol,s,eqv
plnsol,s,1
plnsol,s,3
plnsol,eppl,eqv,0,1
!中墙等效应力场
esel,s,mat,,13
plnsol,s,eqv
plnsol,s,1
plnsol,s,3
!中导洞初支及左隧道衬砌内力
esel,s,live
esel,r,ename,,3 !选择beam3单元类型etable,fx_i,smisc,1
etable,fx_j,smisc,7
etable,fy_i,smisc,2
etable,fy_j,smisc,8
etable,mz_i,smisc,6
etable,mz_j,smisc,12
plls,mz_i,mz_j, !显示弯矩（为显示比例)
plls,fx_i,fx_j, !显示轴力
plls,fy_i,fy_j, !显示剪力
!查看锚杆轴力
esel,s,live
esel,r,ename,,1
etable,forx,smisc,1
plls,forx,forx,
!查看位移场
file,ch-three-2D,rst
set,3,last
lcdef,1,1 !定义第1个载荷步为工况1 lcoper,sub,1 !当前载荷步的结果减去工况1的结果esel,s,live
esel,u,mat,,3
plnsol,u,sum !显示合位移云图
plnsol,u,x
plnsol,u,y
!载荷步4结果
resume,s4,db
/pbc,all,,0
esel,s,live
esel,u,mat,,13
plnsol,s,eqv
plnsol,s,1
plnsol,s,3
plnsol,eppl,eqv,0,1
!中墙等效应力场
esel,s,mat,,13
plnsol,s,eqv
plnsol,s,1
plnsol,s,3
plnsol,u,y
!左及右隧道衬砌内力
esel,s,live
esel,r,ename,,3 !选择beam3单元类型etable,fx_i,smisc,1
etable,fx_j,smisc,7
etable,fy_i,smisc,2
etable,fy_j,smisc,8
etable,mz_i,smisc,6
etable,mz_j,smisc,12
plls,mz_i,mz_j, !显示弯矩（为显示比例)
plls,fx_i,fx_j, !显示轴力
plls,fy_i,fy_j, !显示剪力
!查看锚杆轴力
esel,all
esel,r,ename,,1
etable,forx,smisc,1
plls,forx,forx,
!查看位移场
file,ch-three-2D,rst
set,4,last
lcdef,1,1 !定义第1个载荷步为工况1 lcoper,sub,1 !当前载荷步的结果减去工况1的结果esel,s,live
esel,u,mat,,3
plnsol,u,sum !显示合位移云图plnsol,u,x
plnsol,u,y。