隧道三维结构该如何建模？

!修改网格及参数（Ⅴ级围岩）
finish
/clear
*afun,deg
x_zs=-50
y_zs=30
x_yx=50
y_yx=-40
num=24
th_s=0
th_f=180-th_s
th=(th_f-th_s)/num
yx_1_x=0
yx_1_y=0
yx_2_x=0
yx_2_y=11.37
yx_3_x=3.34
yx_3_y=0
yx_4_x=-3.34
yx_4_y=0
r_1=5.8
r_2=14.4
r_3=2.67
r_4=2.67
r_jg_1=8.8
r_jg_2=17.4
r_jg_3=5.67
r_jg_4=5.67
!beam3单元厚度（考虑钢纤维、二衬提高0.5）
t1=0.25*1.5
!link1单元半径
r1=12.5e-3
!plane42单元
e_r=0.2e9
u_r=0.35
dens_r=1950
c_r=0.18e6
fi_r=24
!beam3单元
e_c=23e9
u_c=0.2
dens_c=2200
!link1单元
e_a=210e9
u_a=0.25
dens_a=7700
/prep7
!单元类型
et,1,beam3
et,2,plane42,,,2
keyopt,2,3,2
et,3,link1
!实常数
r,1,t1,t1*t1*t1/12,t1
r,2,3.14*r1*r1
!材料属性
!左上支护
mp,ex,1,e_c
mp,prxy,1,u_c
mp,dens,1,dens_c
!右上支护
mp,ex,2,e_c
mp,prxy,2,u_c
mp,dens,2,dens_c
!左下支护
mp,ex,3,e_c
mp,prxy,3,u_c
mp,dens,3,dens_c
!右下支护
mp,ex,4,e_c
mp,prxy,4,u_c
mp,dens,4,dens_c
!左上开挖土体
mp,ex,5,e_r
mp,prxy,5,u_r
mp,dens,5,dens_r
tb,dp,5
tbdata,1,c_r,fi_r
!右上开挖土体
mp,ex,6,e_r
mp,prxy,6,u_r
mp,dens,6,dens_r
tb,dp,6
tbdata,1,c_r,fi_r
!左下开挖土体
mp,ex,7,e_r
mp,prxy,7,u_r
mp,dens,7,dens_r
tb,dp,7
tbdata,1,c_r,fi_r
!右下开挖土体
mp,ex,8,e_r
mp,prxy,8,u_r
mp,dens,8,dens_r
tb,dp,8
tbdata,1,c_r,fi_r
!左上锚杆
mp,ex,9,e_a
mp,prxy,9,u_a
mp,dens,9,dens_a
!右上锚杆
mp,ex,10,e_a
mp,prxy,10,u_a
mp,dens,10,dens_a
!左临时支撑
mp,ex,11,e_c
mp,prxy,11,u_c
mp,dens,11,dens_c
!上临时支撑
mp,ex,12,e_c
mp,prxy,12,u_c
mp,dens,12,dens_c
!右临时支撑
mp,ex,13,e_c
mp,prxy,13,u_c
mp,dens,13,dens_c
!下临时支撑
mp,ex,14,e_c
mp,prxy,14,u_c
mp,dens,14,dens_c
!围岩
mp,ex,15,e_r
mp,prxy,15,u_r
mp,dens,15,dens_r
tb,dp,15
tbdata,1,c_r,fi_r
!实体建模
k,1,yx_1_x,yx_1_y
k,2,yx_2_x,yx_2_y
k,3,yx_3_x,yx_3_y
k,4,yx_4_x,yx_4_y
k,101,yx_1_x,yx_1_y,10
k,102,yx_2_x,yx_2_y,10
k,103,yx_3_x,yx_3_y,10
k,104,yx_4_x,yx_4_y,10
circle,1,r_1,101
circle,2,r_2,102
circle,3,r_3,103
circle,4,r_4,104
lcsl,all
lsel,s,,,17,22,1
lsel,a,,,46,48,1
lsel,a,,,27,28,1
lsel,a,,,43
lsel,inve
ldel,all
lsel,all
lcomb,17,18
lcomb,19,21
lcomb,43,27
lcomb,28,46
cm,zh,line
circle,1,r_jg_1,101
circle,2,r_jg_2,102
circle,3,r_jg_3,103
circle,4,r_jg_4,104
cmsel,s,zh,line
lsel,inve
lcsl,all
lsel,s,,,18,21,3
lsel,a,,,23,26,1
lsel,a,,,33,34,1
lsel,a,,,58,59,1
lsel,a,,,55,62,7
cmsel,a,zh,line
lsel,inve
ldele,all,,,1
lsel,all
lcomb,18,21
lcomb,23,25
lcomb,55,33
lcomb,34,58
allsel
kdele,101,104,1
rectng,x_zs,x_yx,y_zs,y_yx
asbl,all,all
!利用工作平面切割围岩
allsel
asel,s,,,3,4,1
wpcsys
wpoff,-1.8*r_jg_1
wprota,,,-90
asbw,all
wpcsys
wprota,,,-90
asbw,all
wpcsys
wpoff,1.8*r_jg_1
wprota,,,-90
asbw,all
wpcsys
wpoff,,1.7*r_jg_1
wprota,,-90
asbw,all
wpcsys
wprota,,-90
asbw,all
wpcsys
wpoff,,-1.4*r_jg_1
wprota,,-90
asbw,all
wpcsys
allsel
/pnum,line,1
/pnum,area,1
!锚杆线
allsel
asel,s,,,5,6,1
wprota,,-90
wprota,,,-th_s
*do,i,1,num,1
wprota,,,-th
asbw,all
*enddo
allsel
!临时支撑
allsel
larc,81,12,7,1.4*r_jg_1
l,7,5
lsel,s,,,19,23,4
asbl,2,all
!围岩分块，利于划分网格
allsel
lsel,none
l,42,92
asbl,19,all
lsel,none
l,44,70
asbl,18,all
lsel,none
l,55,53
asbl,25,all
lsel,none
l,57,54
asbl,24,all
!赋单元属性及划分网格
allsel
aplot
lesize,131,,,10
lesize,132,,,11
lesize,129,,,10
lesize,130,,,6
!左上土体
!aclear,46
allsel
asel,s,,,46
aplot
lsla,s
csys,1
lsel,r,loc,x,r_1
lplot
lesize,all,,,1
cm,zh_zs,line
lsel,r,loc,y,13*th,18*th
lccat,all
cmsel,s,zh_zs,line
lsel,r,loc,y,18*th,24*th
lccat,all
csys,0
mat,5
type,2
mshape,0
mshkey,1
amesh,all
lsel,s,lccat
ldele,all
!右上土体
allsel
!aclear,47
asel,s,,,47
aplot
lsla,s
csys,1
lsel,r,loc,x,r_1
lplot
lesize,all,,,1
cm,zh_ys,line
lsel,r,loc,y,0*th,7*th
lccat,all
cmsel,s,zh_ys,line
lsel,r,loc,y,7*th,13*th
lccat,all
csys,0
mat,6
type,2
mshape,0
mshkey,1
amesh,all
lsel,s,lccat
ldele,all
!左下土体
!alcear,45
allsel
lesize,27,,,8
lesize,48,,,2
lesize,22,,,2
lccat,22,48
mat,7
type,2
mshape,0
mshkey,1
amesh,45
lsel,s,lccat
ldele,all
!右下土体
!alcear,6
allsel
lesize,28,,,8
lesize,47,,,2
lesize,20,,,2
lccat,20,47
mat,8
type,2
mshape,0
mshkey,1
amesh,6
lsel,s,lccat
ldele,all
!临时支撑
allsel
mat,11
real,1
type,1
lmesh,131
mat,12
lmesh,132
mat,13
lmesh,129
mat,14
lmesh,130
!支护
allsel
cmsel,s,zh_zs
latt,1,1,1
lmesh,all
cmsel,s,zh_ys
latt,2,1,1
lmesh,all
asel,s,,,45
lsla,s
lsel,u,,,130,131,1
cm,zh_zx,line
latt,3,1,1
lmesh,all
asel,s,,,6
lsla,s
lsel,u,,,129,130,1
cm,zh_yx,line
latt,4,1,1
lmesh,all
!左上锚杆和锚区围岩
csys,1
asel,s,loc,x,r_1,r_jg_1
asel,r,loc,y,13*th,24*th
lsla,s
lsel,r,loc,x,(r_1+r_jg_1)/2
csys,0
cm,mg_zs,line
mat,9
real,2
type,3
lesize,all,,,4,1.8
lmesh,all
mat,15
type,2
mshape,0
mshkey,1
amesh,all
!右上锚杆和锚区围岩
csys,1
asel,s,loc,x,r_1,r_jg_1
asel,r,loc,y,0*th,13*th
lsla,s
lsel,r,loc,x,(r_1+r_jg_1)/2
csys,0
cm,mg_ys,line
mat,10
real,2
type,3
lesize,all,,,4,1.6
lmesh,all
mat,15
type,2
mshape,0
mshkey,1
amesh,all
!aclear,22,23,1
allsel
lesize,15,,,4,1.6
lcomb,26,62
lesize,26,,,6
lesize,33,,,6
amap,23,12,34,27,7
allsel
lcomb,24,59
lesize,24,,,6
lesize,34,,,6
amap,22,12,34,21,5
!周围围岩
!aclear,all
allsel
asel,s,loc,x,-1.8*r_jg_1,1.8*r_jg_1
asel,r,loc,y,-1.4*r_jg_1,1.7*r_jg_1
lsel,s,ext
lesize,all,,,6
lsel,s,,,42,65,23
lsel,a,,,19,23,4
lsel,a,,,52,53,1
lsel,a,,,133,134,1
lesize,all,,,9,1.2,1
asll,s
csys,1
lsel,s,loc,x,r_jg_1
lsel,r,loc,y,0*th,6*th
lccat,all
lsel,s,loc,x,r_jg_1
lsel,r,loc,y,6*th,12*th
lccat,all
lsel,s,loc,x,r_jg_1
lsel,r,loc,y,12*th,18*th
lccat,all
lsel,s,loc,x,r_jg_1
lsel,r,loc,y,18*th,24*th
lccat,all
csys,0
mat,15
type,2
mshape,0
mshkey,1
amesh,all
lsel,s,lccat
ldele,all
!其他围岩
!aclear,all
asel,s,loc,x,-100*r_jg_1,-1.8*r_jg_1
lsla,s
lsel,r,tan1,y,0
lesize,all,,,16,0.5,1
asel,s,loc,x,1.8*r_jg_1,100*r_jg_1
lsla,s
lsel,r,tan1,y,0
lesize,all,,,16,0.5,1
asel,s,loc,y,1.7*r_jg_1,100*r_jg_1
lsla,all
lsel,r,tan1,x,0
lesize,all,,,8,0.5,1
asel,s,loc,y,-100*r_jg_1,-1.4*r_jg_1
lsla,all
lsel,r,tan1,x,0
lesize,all,,,13,0.5,1
allsel
asel,s,loc,x,-1.8*r_jg_1,1.8*r_jg_1
asel,r,loc,y,-1.4*r_jg_1,1.7*r_jg_1
asel,inve
mat,15
type,2
mshape,0
mshkey,1
amesh,all
allsel
nummrg,node
nummrg,elem
nummrg,kp
numcmp,all
!施加约束
allsel
nsel,s,loc,x,x_zs
nsel,a,loc,x,x_yx
d,all,ux
nsel,s,loc,y,y_yx
d,all,uy
!指定求解选项
/solu
nlgeom,on
nropt,full
pred,on
lnsrch,on
nsubst,10,50,10
!荷载步1,自重应力场
allsel
esel,s,type,,1,3,2
ekill,all
allsel
acel,,9.8
allsel
solve
save,step1,db
esel,s,live
eplot
!荷载步2,左上开挖
!提取节点力
/post1
esel,s,mat,,5
nsle,s
nsel,r,ext
*get,n,node,,count
*dim,n_nf,array,n,3
*do,i,1,n
*get,n_nf(i,1),node,,num,min
nsel,u,,,n_nf(i,1)
*enddo
!*status,n_nf
*do,i,1,n
nsel,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_nf
parsav,all,002,PARM
/solu
antype,,rest
!杀死左上土体单元
allsel
esel,s,mat,,5
ekill,all
!施加等效释放节点力
parres,,002,PARM
esel,s,mat,,5
nsle,s
nsel,r,ext
*do,i,1,n
f,n_nf(i,1),fx,0.65*n_nf(i,2)
f,n_nf(i,1),fy,0.65*n_nf(i,3)
*enddo
!约束死单元自由度
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all
allsel
solve
save,step2,db
esel,s,live
eplot
!荷载步3
!删除节点力，激活支护、锚杆、临时支撑
allsel
esel,s,mat,,5
nsle,s
nsel,r,ext
fdele,all,all
esel,s,mat,,1,9,8
esel,a,mat,,11,12,1
ealive,all
allsel
solve
save,step3,db
!荷载步4,右上开挖，应力释放
!提取节点力
/post1
esel,s,mat,,6
nsle,s
nsel,r,ext
*get,n,node,,count
*dim,n_nf,array,n,3
*do,i,1,n
*get,n_nf(i,1),node,,num,min
nsel,u,,,n_nf(i,1)
*enddo
!*status,n_nf
*do,i,1,n
nsel,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_nf
parsav,all,002,PARM
/solu
antype,,rest
!杀死右上土体单元
allsel
ealive,all
esel,s,mat,,2,5,1
esel,a,mat,,10
esel,a,mat,,13,14,1
esel,a,mat,,6
ekill,all
!施加等效释放节点力
parres,,002,PARM
esel,s,mat,,6
nsle,s
nsel,r,ext
*do,i,1,n
f,n_nf(i,1),fx,0.65*n_nf(i,2)
f,n_nf(i,1),fy,0.65*n_nf(i,3)
*enddo
!约束死单元自由度
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all
allsel
solve
save,step4,db
esel,s,live
eplot
!荷载步5
!删除节点力，激活支护、锚杆、临时支撑
allsel
esel,s,mat,,6
nsle,s
nsel,r,ext
fdele,all,all
esel,s,mat,,2,10,8
esel,a,mat,,13
ealive,all
allsel
solve
save,step5,db
!荷载步6,左下开挖，应力释放
!提取节点力
/post1
esel,s,mat,,7
nsle,s
nsel,r,ext
*get,n,node,,count
*dim,n_nf,array,n,3
*do,i,1,n
*get,n_nf(i,1),node,,num,min
nsel,u,,,n_nf(i,1)
*enddo
!*status,n_nf
*do,i,1,n
nsel,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_nf
parsav,all,002,PARM
/solu
antype,,rest
!杀死左下土体单元
allsel
ealive,all
esel,s,mat,,3,6,1
esel,a,mat,,14
esel,a,mat,,7
ekill,all
!施加等效释放节点力
parres,,002,PARM
esel,s,mat,,7
nsle,s
nsel,r,ext
*do,i,1,n
f,n_nf(i,1),fx,0.65*n_nf(i,2)
f,n_nf(i,1),fy,0.65*n_nf(i,3)
*enddo
!约束死单元自由度
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all
allsel
solve
save,step6,db
esel,s,live
eplot
!荷载步7
!删除节点力，激活支护、临时支撑
esel,s,mat,,7
nsle,s
nsel,r,ext
fdele,all,all
esel,s,mat,,3,14,11
ealive,all
allsel
solve
save,step7,db
!荷载步8,右下开挖，应力释放
!提取节点力
/post1
esel,s,mat,,8
nsle,s
nsel,r,ext
*get,n,node,,count
*dim,n_nf,array,n,3
*do,i,1,n
*get,n_nf(i,1),node,,num,min
nsel,u,,,n_nf(i,1)
*enddo
!*status,n_nf
*do,i,1,n
nsel,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_nf
parsav,all,002,PARM
/solu
antype,,rest
!杀死右下土体单元
allsel
ealive,all
esel,s,mat,,4,7,1
esel,a,mat,,8
ekill,all
!施加等效释放节点力
parres,,002,PARM
esel,s,mat,,8
nsle,s
nsel,r,ext
*do,i,1,n
f,n_nf(i,1),fx,0.65*n_nf(i,2)
f,n_nf(i,1),fy,0.65*n_nf(i,3)
*enddo
!约束死单元自由度
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all
allsel
solve
save,step8,db
esel,s,live
eplot
!荷载步9
!删除节点力，激活支护
esel,s,mat,,8
nsle,s
nsel,r,ext
fdele,all,all
esel,s,mat,,4
ealive,all
allsel
solve
save,step9,db
!荷载步10，拆除临时支撑
esel,s,mat,,11,14,1
ekill,all
allsel
solve
save,step10,db
esel,s,live
eplot
!后处理
!荷载步3
set,3,last
!应力场
esel,s,live
plnsol,s,1
plnsol,s,2
plnsol,s,3
plnsol,s,int
plnsol,eqv
!支护内力
set,10,last
esel,s,type,,1
etable,fx_i,smisc,1
etable,fx_j,smisc,7
plls,fx_i,fx_j,0.2
etable,fy_i,smisc,2
etable,fy_j,smisc,8
plls,fy_i,fy_j,0.2
esel,s,type,,1
etable,mz_i,smisc,6
etable,mz_j,smisc,12
plls,mz_i,mz_j,-0.2
!锚杆内力
set,10,last
esel,s,type,,3
etable,sa,ls,1
plls,sa,sa,0.1
!锚杆应变
!etable,ybae,lepel,1
!plls,ybae,ybae,0.1
!位移场(减去初始位移场)
resume,step10,db
set,10,last
wpoff,100
/triad,off
esel,s,live
/replot
lcdef,1,1,last
lcdef,10,10,last
lcase,10
lcoper,sub,1
plnsol,u,y
/post26
numvar,30
!拱顶位移
nsol,2,115,u,y
!水平收敛
nsol,3,23,u,x
nsol,4,103,u,x
add,5,3,4,,,,,1,-1
!地表沉降
nsol,11,1511,u,y
nsol,12,1902,u,y
nsol,13,1904,u,y
nsol,14,1893,u,y
nsol,15,1942,u,y
nsol,16,1944,u,y
nsol,17,1750,u,y
prvar,2,3,4,5
prvar,11,12,13,14
prvar,15,16,17
这里有个二维的，三维的可以拉伸，也可以按平面计算（按隧道无限长）