/COM,structural !指定结构分析
/TITLE ,Tunnel Construct Analysis !定义工作标题
/FILNAM,SUPPORT,1 !定义工作文件名
/prep7 !前处理
!定义分析参数
*afun,deg ! 设定角度单位
! 定义衬砌、围岩、回填材料的单元类型（四级围岩）
ET,1,BEAM3 ! 衬砌单元，
KEYOPT,1,6,1
ET,2,PLANE42 ! 围岩
KEYOPT,2,3,2 ! 平面应变模型
ET,3,PLANE42 ! 回填层
KEYOPT,2,3,2
!定义材料属性
MP,EX,1,3.565E10
MP,PRXY,1,0.2 ! 衬砌c30的弹性模量和泊松比
!围岩材料
MP,EX,2,1.3E9
MP,PRXY,2,0.35 ! 四级围岩参数
MP,DENS,2,2150
TB,DP,2
TBDATA,1,2E5,27 !粘聚力和内摩擦角
!挖去土体材料
MP,EX,3,1.3E9
MP,PRXY,3,0.35
MP,DENS,3,2150
TB,DP,3
TBDATA,1,2E5,27 !粘聚力和内摩擦角
! 回填层材料参数
MP,EX,4,12E9 !研究 弹性模量 四级围岩最大值到衬砌之间变化 也可以外延
MP,PRXY,4,0.35 ! 区间取值 ，研究什么参数就改变那一个，区间取规范上围岩参
MP,DENS,4,2000 ! 回填层材料密度 规范
TB,DP,4
TBDATA,1,2E5,27
!定义实常数
R,1,0.485,0.485*0.485*0.485/12,0.485 !衬砌实常数
!建立几何模型
k,1,0,4.906,0
k,2,-6.495,1.156,0
k,3,-7.163,-0.753,0
k,4,-7.224,-1.563,0
k,5,-7.087,-2.774,0
k,6,-6.745,-3.438,0
k,7,-6.125,-3.854,0
k,8,0,-4.906,0
k,9,6.125,-3.854,0
k,10,6.745,-3.438,0
k,11,7.087,-2.774,0
k,12,7.224,-1.563,0
k,13,7.163,-0.753,0
k,14,6.495,1.156,0
k,15,0,5,0
k,16,-6.577,1.203,0
k,17,-7.256,-0.739,0
k,18,-7.318,-1.563,0
k,19,-7.179,-2.795,0
k,20,-6.815,-3.5,0
k,21,-6.156,-3.943,0
k,22,0,-5,0
k,23,6.156,-3.943,0
k,24,6.815,-3.5,0
k,25,7.179,-2.795,0
k,26,7.318,-1.563,0
k,27,7.256,-0.739,0
k,28,6.577,1.203,0
larc,1,2,3,7.5
larc,2,3,4,5.438
larc,3,4,5,5.438
larc,4,5,6,5.438
larc,5,6,7,1.5
larc,6,7,8,1.5
larc,7,8,9,18.366
larc,8,9,10,18.366
larc,9,10,11,1.5
larc,10,11,12,1.5
larc,11,12,13,5.438
larc,12,13,14,5.438
larc,13,14,1,5.438
larc,14,1,2,7.5
larc,15,16,17,7.594
larc,16,17,18,5.532
larc,17,18,19,5.532
larc,18,19,20,5.532
larc,19,20,21,1.594
larc,20,21,22,1.594
larc,21,22,23,18.46
larc,22,23,24,18.46
larc,23,24,25,1.594
larc,24,25,26,1.594
larc,25,26,27,5.532
larc,26,27,28,5.532
larc,27,28,15,5.532
larc,28,15,16,7.594
！细分隧道分析线模型
k,29,-50,55,0
k,30,-50,-55,0
k,31,50,-55,0
k,32,50,55,0
K,33,10,8,0
K,34,10,-8,0
K,35,-10,8,0
K,36,-10,-8,0
k,37,-50,-8,0
k,38,-50, 8,0
k,39,50,8,0
k,40,50,-8,0
k,41,-10,55,0
k,42,10,55,0
k,43,10,-55,0
k,44,-10,-55,0
l,30,44 !从下到上横线
l,44,43
l,43,31
l,37,
36
l,36,34
l,34,40
l,38,35
l,35,33
l,33,39
l,29,41
l,41,42
l,42,32 !从左到右竖线
l,29,38
l,38,37
l,37,30
l,41,35
l,35,36
l,36,44
l,42,33
l,33,34
l,34,43
l,32,39
l,39,40
l,40,31
a,29,30,31,32 !通过4个关键点生成一个面积区域
asbl,all,all !通过线分割面生成新面
adele,1 !删除面
NUMCMP,AREA !压缩面编号
!划分网格生成有限元模型
!划分梁单元
lsel,s,,,1,7,6
lsel,a,,,8,14,6
lesize,all,,,8
lsel,s,,,2,13,11
lesize,all,,,3
lsel,s,,,3,6,1
lsel,a,,,9,12,1
lesize,all,,,2
mat,1 !指定梁单元材料特性
type,1
allsel
lsel,s,,,1,14,1
lmesh,all !划分所有线
allsel
!划分开挖掉土体单元网格
mat,3
type,2
mshkey,0 !设定自由网格划分
mshape,0 !设定四边形网格划分
esize,
amesh,9
!划分9号面最里面的面
!划分回填层单元网格
!设置网格份数
lsel,s,,,15,21,6
lsel,a,,,22,28,6
lesize,all,,,8
lsel,s,,,16,27,11
lesize,all,,,3
lsel,s,,,17,20,1
lsel,a,,,23,26,1
lesize,all,,,2
mat,4 !赋予回填层属性
type,2
mshkey,0
mshape,0
allsel
asel,s,,,11 !回填层面积编号11
amesh,11 !划分11面
allsel
!划分围岩网格
!设置网格份数
lsel,s,,,45,48,3 !选择线
lesize,all,,,20 !把所选择线分为20份
lsel,s,,,33,36,3
lesize,all,,,22
mat,3 !赋予围岩单元属性
type,2
mshkey,0 !设定自由网格划分
mshape,0 !设定四边形网格划分
allsel
asel,s,,,10
amesh,all
allsel
lsel,s,,,42,51,9
lesize,all,,,20
lsel,s,,,30,39,9
lesize,all,,,22
lsel,s,,,41,50,3
lsel,a,,,43,52,3
lesize,all,,,30
lsel,s,,,29,38,3
lsel,a,,,31,40,3
lesize,all,,,25
mat,3 !赋予围岩单元属性
type,2
mshkey,0 !设定自由网格划分
mshape,0 !设定四边形网格划分
allsel
asel,s,,,1,8,1
amesh,all
allsel !划分所有面
!保存网格模型
save,Tunnel-grid.db
!施加约束
csys,0
nsel,s,loc,x,-50 !选择X=-50线上所有节点
nsel,a,loc,x,50
!选择X=5线上所有节点
d,all,ux !对所选择节点约束X方向位移
allsel
nsel,s,loc,y,-55 !选择Y=--55线上所有节点
d,all,uy !对所选择节点约束Y方向位移
!施加重力加速度
acel,,9.8
/solu
!初始应力模拟
antype,static !
nsubst,100
pred,on !
nropt,full
nlgeom,on !
lnsrch,on !
outres,all,last !
cnvtol,f,,0.02,2,0.5 !
cnvtol,m,,0.01,2,1
esel,s,mat,,1 !选择材料号为1单元
ekill,all !杀死所选择的单元
allsel
solve !求解
!隧道开挖与支护模拟
Antype,,rest,,
esel,s,mat,,1
ekill,all
allsel
asel,s,,,9
esla,s
ekill,all
esel,s,mat,,1 !
ealive,all
asel,s,,,11
esla,s
mpchg,4,all !
allsel
solve !求解
/post1
lcdef,1,1,,
lcdef,2,2,, !2代表步数
lcase,2,
lcoper,sub,1,,,
esel,s,live
/post1
!初始应力模拟结果分析
!显示位移云图
plnsol,u,sum,
esel,s,live
!绘制总位移矢量云图
plnsol,u,x
!绘制X方向位移云图
plnsol,u,y
!显示应力云图
plnsol,s,x !绘制X方向应力云图
plnsol,s,y !绘制Y方向应力云图
plnsol,s,z
plnsol,s,1
plnsol,s,2 !绘制第2主应力云图
plnsol,s,3
plnsol,s,eqv !绘制等效应力云图
!显示梁支护内力
esel,s,mat,,1 !选择材料号为1的单元
etable,IMOMEMT,SMISC,6 !定义弯矩表
etable,JMOMEMT,SMISC,12
etable,ISHEAR,SMISC,2 !定义剪力表
etable,JSHEAR,SMISC,8
etable,ZHOULI-I,SMISC,1 !定义轴力表
etable,ZHOULI-J,SMISC,7
/TITLE,BENDING MOMENT distribution !定义弯矩分布标题
PLLS,IMOMEMT,JMOMEMT,-0.5 !绘制弯矩分布图
/TITLE,SHEAR force distribution !定义剪力分布标题
PLLS,ISHEAR,JSHEAR,0.5 !绘制剪力分布图
/TITLE,ZHOULI force distribution !定义轴力分布标题
PLLS,ZHOULI-I,ZHOULI-J,0.1
!绘制轴力分布图
/COM,structural !指定结构分析
/TITLE ,Tunnel Construct Analysis !定义工作标题
/FILNAM,SUPPORT,1 !定义工作文件名
/prep7 !前处理
!定义分析参数
*afun,deg ! 设定角度单位
! 定义衬砌、围岩、回填材料的单元类型（四级围岩）
ET,1,BEAM3 ! 衬砌单元，
KEYOPT,1,6,1
ET,2,PLANE42 ! 围岩
KEYOPT,2,3,2 ! 平面应变模型
ET,3,PLANE42 ! 回填层
KEYOPT,2,3,2
!定义材料属性
MP,EX,1,3.565E10
MP,PRXY,1,0.2 ! 衬砌c30的弹性模量和泊松比
!围岩材料
MP,EX,2,1.3E9
MP,PRXY,2,0.35 ! 四级围岩参数
MP,DENS,2,2150
TB,DP,2
TBDATA,1,2E5,27 !粘聚力和内摩擦角
!挖去土体材料
MP,EX,3,1.3E9
MP,PRXY,3,0.35
MP,DENS,3,2150
TB,DP,3
TBDATA,1,2E5,27 !粘聚力和内摩擦角
! 回填层材料参数
MP,EX,4,12E9 !研究 弹性模量 四级围岩最大值到衬砌之间变化 也可以外延
MP,PRXY,4,0.35 ! 区间取值 ，研究什么参数就改变那一个，区间取规范上围岩参
MP,DENS,4,2000 ! 回填层材料密度 规范
TB,DP,4
TBDATA,1,2E5,27
!定义实常数
R,1,0.485,0.485*0.485*0.485/12,0.485 !衬砌实常数
!建立几何模型
k,1,0,4.906,0
k,2,-6.495,1.156,0
k,3,-7.163,-0.753,0
k,4,-7.224,-1.563,0
k,5,-7.087,-2.774,0
k,6,-6.745,-3.438,0
k,7,-6.125,-3.854,0
k,8,0,-4.906,0
k,9,6.125,-3.854,0
k,10,6.745,-3.438,0
k,11,7.087,-2.774,0
k,12,7.224,-1.563,0
k,13,7.163,-0.753,0
k,14,6.495,1.156,0
k,15,0,5,0
k,16,-6.577,1.203,0
k,17,-7.256,-0.739,0
k,18,-7.318,-1.563,0
k,19,-7.179,-2.795,0
k,20,-6.815,-3.5,0
k,21,-6.156,-3.943,0
k,22,0,-5,0
k,23,6.156,-3.943,0
k,24,6.815,-3.5,0
k,25,7.179,-2.795,0
k,26,7.318,-1.563,0
k,27,7.256,-0.739,0
k,28,6.577,1.203,0
larc,1,2,3,7.5
larc,2,3,4,5.438
larc,3,4,5,5.438
larc,4,5,6,5.438
larc,5,6,7,1.5
larc,6,7,8,1.5
larc,7,8,9,18.366
larc,8,9,10,18.366
larc,9,10,11,1.5
larc,10,11,12,1.5
larc,11,12,13,5.438
larc,12,13,14,5.438
larc,13,14,1,5.438
larc,14,1,2,7.5
larc,15,16,17,7.594
larc,16,17,18,5.532
larc,17,18,19,5.532
larc,18,19,20,5.532
larc,19,20,21,1.594
larc,20,21,22,1.594
larc,21,22,23,18.46
larc,22,23,24,18.46
larc,23,24,25,1.594
larc,24,25,26,1.594
larc,25,26,27,5.532
larc,26,27,28,5.532
larc,27,28,15,5.532
larc,28,15,16,7.594
！细分隧道分析线模型
k,29,-50,55,0
k,30,-50,-55,0
k,31,50,-55,0
k,32,50,55,0
K,33,10,8,0
K,34,10,-8,0
K,35,-10,8,0
K,36,-10,-8,0
k,37,-50,-8,0
k,38,-50, 8,0
k,39,50,8,0
k,40,50,-8,0
k,41,-10,55,0
k,42,10,55,0
k,43,10,-55,0
k,44,-10,-55,0
l,30,44 !从下到上横线
l,44,43
l,43,31
l,37,
36
l,36,34
l,34,40
l,38,35
l,35,33
l,33,39
l,29,41
l,41,42
l,42,32 !从左到右竖线
l,29,38
l,38,37
l,37,30
l,41,35
l,35,36
l,36,44
l,42,33
l,33,34
l,34,43
l,32,39
l,39,40
l,40,31
a,29,30,31,32 !通过4个关键点生成一个面积区域
asbl,all,all !通过线分割面生成新面
adele,1 !删除面
NUMCMP,AREA !压缩面编号
!划分网格生成有限元模型
!划分梁单元
lsel,s,,,1,7,6
lsel,a,,,8,14,6
lesize,all,,,8
lsel,s,,,2,13,11
lesize,all,,,3
lsel,s,,,3,6,1
lsel,a,,,9,12,1
lesize,all,,,2
mat,1 !指定梁单元材料特性
type,1
allsel
lsel,s,,,1,14,1
lmesh,all !划分所有线
allsel
!划分开挖掉土体单元网格
mat,3
type,2
mshkey,0 !设定自由网格划分
mshape,0 !设定四边形网格划分
esize,
amesh,9
!划分9号面最里面的面
!划分回填层单元网格
!设置网格份数
lsel,s,,,15,21,6
lsel,a,,,22,28,6
lesize,all,,,8
lsel,s,,,16,27,11
lesize,all,,,3
lsel,s,,,17,20,1
lsel,a,,,23,26,1
lesize,all,,,2
mat,4 !赋予回填层属性
type,2
mshkey,0
mshape,0
allsel
asel,s,,,11 !回填层面积编号11
amesh,11 !划分11面
allsel
!划分围岩网格
!设置网格份数
lsel,s,,,45,48,3 !选择线
lesize,all,,,20 !把所选择线分为20份
lsel,s,,,33,36,3
lesize,all,,,22
mat,3 !赋予围岩单元属性
type,2
mshkey,0 !设定自由网格划分
mshape,0 !设定四边形网格划分
allsel
asel,s,,,10
amesh,all
allsel
lsel,s,,,42,51,9
lesize,all,,,20
lsel,s,,,30,39,9
lesize,all,,,22
lsel,s,,,41,50,3
lsel,a,,,43,52,3
lesize,all,,,30
lsel,s,,,29,38,3
lsel,a,,,31,40,3
lesize,all,,,25
mat,3 !赋予围岩单元属性
type,2
mshkey,0 !设定自由网格划分
mshape,0 !设定四边形网格划分
allsel
asel,s,,,1,8,1
amesh,all
allsel !划分所有面
!保存网格模型
save,Tunnel-grid.db
!施加约束
csys,0
nsel,s,loc,x,-50 !选择X=-50线上所有节点
nsel,a,loc,x,50
!选择X=5线上所有节点
d,all,ux !对所选择节点约束X方向位移
allsel
nsel,s,loc,y,-55 !选择Y=--55线上所有节点
d,all,uy !对所选择节点约束Y方向位移
!施加重力加速度
acel,,9.8
/solu
!初始应力模拟
antype,static !
nsubst,100
pred,on !
nropt,full
nlgeom,on !
lnsrch,on !
outres,all,last !
cnvtol,f,,0.02,2,0.5 !
cnvtol,m,,0.01,2,1
esel,s,mat,,1 !选择材料号为1单元
ekill,all !杀死所选择的单元
allsel
solve !求解
!隧道开挖与支护模拟
Antype,,rest,,
esel,s,mat,,1
ekill,all
allsel
asel,s,,,9
esla,s
ekill,all
esel,s,mat,,1 !
ealive,all
asel,s,,,11
esla,s
mpchg,4,all !
allsel
solve !求解
/post1
lcdef,1,1,,
lcdef,2,2,, !2代表步数
lcase,2,
lcoper,sub,1,,,
esel,s,live
/post1
!初始应力模拟结果分析
!显示位移云图
plnsol,u,sum,
esel,s,live
!绘制总位移矢量云图
plnsol,u,x
!绘制X方向位移云图
plnsol,u,y
!显示应力云图
plnsol,s,x !绘制X方向应力云图
plnsol,s,y !绘制Y方向应力云图
plnsol,s,z
plnsol,s,1
plnsol,s,2 !绘制第2主应力云图
plnsol,s,3
plnsol,s,eqv !绘制等效应力云图
!显示梁支护内力
esel,s,mat,,1 !选择材料号为1的单元
etable,IMOMEMT,SMISC,6 !定义弯矩表
etable,JMOMEMT,SMISC,12
etable,ISHEAR,SMISC,2 !定义剪力表
etable,JSHEAR,SMISC,8
etable,ZHOULI-I,SMISC,1 !定义轴力表
etable,ZHOULI-J,SMISC,7
/TITLE,BENDING MOMENT distribution !定义弯矩分布标题
PLLS,IMOMEMT,JMOMEMT,-0.5 !绘制弯矩分布图
/TITLE,SHEAR force distribution !定义剪力分布标题
PLLS,ISHEAR,JSHEAR,0.5 !绘制剪力分布图
/TITLE,ZHOULI force distribution !定义轴力分布标题
PLLS,ZHOULI-I,ZHOULI-J,0.1
!绘制轴力分布图