隧道建模命令流

/TITLE,calculate stress of tunnel'floor

/PREP7  
!定义单元及材料属性
ET,1,plane42

MP,EX,1,25.5e9
MP,NUXY,1,0.25   
MP,DENS,1,2500  

MP,EX,2,150E6
MP,NUXY,2,0.35
MP,DENS,2,1850

TB,DP,2,,,0
TBMODIF,1,1,100000   
TBMODIF,1,2,35
TBMODIF,1,3,35

!土体及衬砌实体模型
CSYS,0         !激活一个预先已知的坐标系统
!外轮廓截面控制关键点
K,1,8.18,-1.98,,
K,2,8.18,2.32,,
K,3,7.51,5.32,,
K,4,5.37,7.97,,
K,5,0,9.8,,
K,6,-5.37,7.97,,
K,7,-7.51,5.32,,
K,8,-8.18,2.32,,   
K,9,-8.18,-1.98,,
K,10,-5.91,-1.98,,
K,11,0,-2.79,,
K,12,5.91,-1.98,,

!内轮廓截面控制关键点
K,13,5.71,-1.26,,
K,14,6.60,-0.25,,
K,15,6.95,3.82,,
K,16,4.81,7.30,,
K,17,0,9.05,,
K,18,-4.81,7.30,,
K,19,-6.95,3.82,,
K,20,-6.60,-0.25,,
K,21,-5.71,-1.26,,
K,22,0,-2.04,,  

!外轮廓线段的生成
lstr,1,2            
larc,2,4,3         
larc,4,6,5         
larc,6,8,7         
lstr,8,9            
lstr,9,10           
larc,10,12,11      
lstr,12,1           

!内轮廓线段的生成
lstr,13,14         
larc,14,16,15      
larc,16,18,17      
larc,18,20,19      
lstr,20,21         
larc,21,13,22      

!将外轮廓线叠合
lcomb,1,2,0
lcomb,4,5,0

!压缩所有实体的编号
numcmp,all

!将内外轮廓对应关键点相连(6条线)
lstr,12,1                 !line 13
lstr,14,3                 !line 14
lstr,16,5                 !line 15
lstr,18,7                 !line 16
lstr,19,8                 !line 17
lstr,11,10                !line 18

!生成周边土体边界关键点
k,21,40,-20
k,22,40,46
k,23,-40,46
k,24,-40,-20
k,25,-30,-20
k,26,30,-20

!将周边土体边界关键点相连成线
*do,i,21,25
  lstr,i,i+1
*enddo
lstr,26,21

!衬砌外轮廓与周边土体关键点相连
lstr,1,21
lstr,3,22
lstr,5,23
lstr,7,24
lstr,8,25
lstr,10,26

!生成衬砌面
*do,i,8,12
  al,i,i+6,i-7,i+5
*enddo
al,7,13,6,18

!生成土体面
*do,i,1,5
al,i,i+25,i+18,i+24
*enddo
al,6,25,24,30

!边墙部位(及对应两侧土体)网格划分时的线段数
lsel,s,,,1
lsel,a,,,3
lsel,a,,,8
lsel,a,,,10
lsel,a,,,19
lsel,a,,,21
lesize,all,,,33,,,,,1

!拱部(及对应两侧土体)网格划分时的线段数
lsel,s,,,2
lsel,a,,,9
lsel,a,,,20
lesize,all,,,40,,,,,1

!墙脚(及对应两侧土体)网格划分时的线段数
lsel,s,,,11                  !第一次用s
lsel,a,,,4                       以后的用a
lsel,a,,,22
lsel,a,,,7
lsel,a,,,6
lsel,a,,,24
lesize,all,,,5,,,,,1

!仰拱(及对应两侧土体)网格划分时的线段数
lsel,s,,,12
lsel,a,,,5
lsel,a,,,23
lesize,all,,,30,,,,,1

!衬砌内外轮廓连线网格划分分段数
lsel,s,,,13,18
lesize,all,,,5,,,,,1

!衬砌与土体连线网格划分时的分段数
lsel,s,,,25,30
lesize,all,,,30,,,,,1

!划分衬砌单元
type,1
mat,1
mshape,0,2D       !指定划分单元的形状
mshkey,1         !指定划分网格方式
asel,s,,,1,6
amesh,all

!划分土体单元
type,1
mat,2
mshape,0,2D
mshkey,1
asel,s,,,7,12
amesh,all

/solu
acel,,9.8,

!施加底边的竖向和两侧水平约束
nsel,s,loc,y,-20
d,all,,,,,,uy,,,,,
nsel,s,loc,x,-40
nsel,a,loc,x,40
d,all,,,,,,ux,,,,,

allsel,all
!solve