ANSYS实例分析连拱隧道开挖命令流.docx
《ANSYS实例分析连拱隧道开挖命令流.docx》由会员分享，可在线阅读，更多相关《ANSYS实例分析连拱隧道开挖命令流.docx（26页珍藏版）》请在冰豆网上搜索。
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,9,,
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,18,,
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,,
!
导洞左侧圈
l,3,13
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
l,i,115-i
*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,103,104
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
lsel,a,line,,11,12,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,,18,52,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,,271,279,8
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
lsel,a,line,,16
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
lesize,128,,,1
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
升级会员 