命令流:
finish /clear
/COM,Structural ! 指定结构分析 /TITLE,Tunnel Construct Modeling Analysis ! 定义工作标题
/FILNAM,support,1 !锚杆
MP,EX,4,200E9 MP,PRXY,4,0.3 MP,DENS,4,7960 !定义实常数
R,1,0.3,0.3*0.3*0.3/12,0.3, !衬砌支护实常数
! 定义工作文件名 ！进入前处理器 /PREP7
!定义单元类型
ET,1,BEAM3 !定义衬砌支护单元
KEYOPT,1,6,1 局部坐标系输出力、弯矩
ET,2,PLANE42 !定义围岩单元
KEYOPT,2,3,2 ! Plane strain
ET,3,LINK1 !定义LINK1单元 !定义材料属性 !衬砌支护
MP,EX,1,2.75E10 MP,PRXY,1,0.2 MP,DENS,1,2500
TB,DP,1 !采用D-P模型 TBDATA,1,2.42E6,54 !围岩材料
MP,EX,2,1.3E9 MP,PRXY,2,0.3 MP,DENS,2,2450 TB,DP,2
TBDATA,1,0.15E6,39 ！内粘聚力 内摩擦角 !挖去土体材料 MP,EX,3,4.5E8 MP,PRXY,3,0.3 MP,DENS,3,2450 TB,DP,3
TBDATA,1,0.1E6,24.3
R,2,3.14*0.025*0.025/4,, !锚杆实常数 NUMSTR,KP,1 NUMSTR,LINE,1 NUMSTR,AREA,1 NUMSTR,VOLU,1 !CSYS,0 k,1,0,14.475 k,2,0,-3.525
k,3,5.6098,-2.6285 k,4,1.4247,-3.4685 k,5,5.1423,-1.2032 k,6,6.6179,-1.4726 k,7,5.9709,-2.4536 k,8,-1.4488,0 k,9,6.7512,0
k,10,6.6761,-1.1071 k,11,0.9012,0
k,12,4.8631,4.3042 k,13,6.6267,1.2006 k,14,0,-0.9791 k,15,0,6.2016
k,16,3.8021,5.1124 k,17,0,-7.525
k,18,6.8564,-6.4293 k,19,1.7413,-7.456 k,20,10.5529,-2.191 k,21,8.1805,-5.7879 k,22,10.7512,0
k,23,10.6395,-1.6472 k,24,7.572,7.2473 k,25,10.5415,2.0215 k,26,0,10.2016 k,27,5.9201,8.5057 k,28,-4.8631,4.3042
k,29,-1.3279,6.0778 k,30,-0.9012,0 k,31,-6.7512,0
k,32,-5.6621,3.3995 k,33,1.4488,0
k,34,-6.6179,-1.4726 k,35,-6.7429,-0.37 k,36,-5.1423,-1.2032 k,37,-5.6098,-2.6285 k,38,-6.4937,-1.8542 k,39,-4.2384,-3.0189 k,40,-6.8564,-6.4293 k,41,-5.1803,-6.9064 k,42,-10.5529,-2.191 k,43,-10.0974,-3.5901 k,44,-10.7512,0
k,45,-10.7388,-0.5505 k,46,-7.572,7.2473 k,47,-8.9174,5.7239 k,48,-2.0677,10.0088 k,49,12.8586,50 k,50,29.6587,50 k,51,45.6987,50 k,52,0,0
k,53,-46.4587,50 k,54,-12.8512,50
k,55,-12.8512,12.3016 k,56,-46.4587,12.3016 k,57,-12.8512,-9.625 k,58,-46.4587,-9.625 k,59,-12.8512,-25.525 k,60,-46.4587,-25.525 k,61,12.8586,-25.525 k,62,12.8586,-9.625 k,63,12.8586,12.3016 k,64,45.6987,-25.525 k,65,45.6987,-9.625 k,66,45.6987,12.3016 k,67,0,12.3016 k,68,0,-9.625 k,69,12.8586,0 k,70,-12.8512,0 l,49,50, l,50,51, l,31,52, l,52,9, l,53,54, l,55,56, l,57,58, l,59,60, l,59,61, l,57,62, l,55,63, l,54,49, l,61,64, l,62,65, l,63,66, l,51,66, l,66,65, l,65,64, l,60,58, l,58,56, l,56,53, l,59,57, l,57,55, l,55,54, l,61,62, l,62,63, l,63,49, l,26,67, l,17,68, l,22,69, l,44,70, LARC,2,3,4 LARC,3,6,7 LARC,6,9,10 LARC,9,12,13 LARC,12,15,16 LARC,17,18,19 LARC,18,20,21 LARC,20,22,23 LARC,22,24,25 LARC,24,26,27 LARC,15,28,29 LARC,28,31,32 LARC,31,34,35 LARC,34,37,38 LARC,37,2,39
LARC,40,17,41 LARC,42,40,43 LARC,44,42,45 LARC,46,44,47 LARC,26,46,48 AL,5,24,6,21 AL,12,27,11,24 AL,1,2,16,15,27
a,56,60,64,66 !通过4个关键点生成一个面积区域
asbl,all,all !通过线分割面生成新面
adele,4 !删除面 NUMCMP,AREA !压缩面编号 !生成锚杆
asel,s,,,15 !选择面15
csys,1 !激活柱坐标系
wprota,,-90 !工作平面绕X轴旋转-90度
*do,i,1,30,1 !循环控制
wprota,,,-6 !工作平面绕Y轴旋转-6度asbw,all !用工作平面切割所选择的所有面 *enddo
Lcomb,89,92,0 Lcomb,186,190,0 Lcomb,64,70,0 Lcomb,82,77,0 Lcomb,204,199,0 Lcomb,211,216,0 Lcomb,99,95,0 Lcomb,181,185,0 Lcomb,49,65,0 Lcomb,36,101,0 Lcomb,178,187,0 Lcomb,223,217,0
l,55,157 l,57,96 l,62,155 l,63,98 allsel !删除面
FLST,2,4,5,ORDE,2 FITEM,2,10 FITEM,2,-13 ADELE,P51X !生成面
FLST,2,11,4 FITEM,2,60 FITEM,2,28 FITEM,2,146 FITEM,2,152 FITEM,2,158 FITEM,2,164 FITEM,2,170 FITEM,2,176 FITEM,2,182 FITEM,2,181 FITEM,2,65 AL,P51X FLST,2,10,4 FITEM,2,56 FITEM,2,65 FITEM,2,194 FITEM,2,200 FITEM,2,197 FITEM,2,212 FITEM,2,209 FITEM,2,214 FITEM,2,227 FITEM,2,31 AL,P51X FLST,2,10,4 FITEM,2,31 FITEM,2,26 FITEM,2,49 FITEM,2,72 FITEM,2,78 FITEM,2,84 FITEM,2,90 FITEM,2,87 FITEM,2,70 FITEM,2,57 AL,P51X
FLST,2,11,4 FITEM,2,36 FITEM,2,41 FITEM,2,104 FITEM,2,110 FITEM,2,116 FITEM,2,122 FITEM,2,128 FITEM,2,143 FITEM,2,29 FITEM,2,58 FITEM,2,70 AL,P51X FLST,2,11,4 FITEM,2,137 FITEM,2,37 FITEM,2,149 FITEM,2,155 FITEM,2,161 FITEM,2,167 FITEM,2,173 FITEM,2,178 FITEM,2,82 FITEM,2,59 FITEM,2,29 AL,P51X FLST,2,10,4 FITEM,2,82 FITEM,2,62 FITEM,2,30 FITEM,2,229 FITEM,2,217 FITEM,2,218 FITEM,2,202 FITEM,2,206 FITEM,2,43 FITEM,2,38 AL,P51X FLST,2,10,4 FITEM,2,92 FITEM,2,63 FITEM,2,30 FITEM,2,55 FITEM,2,40 FITEM,2,11 FITEM,2,68 FITEM,2,75 FITEM,2,80 FITEM,2,96 AL,P51X FLST,2,11,4 FITEM,2,28 FITEM,2,61 FITEM,2,92 FITEM,2,95 FITEM,2,100 FITEM,2,107 FITEM,2,113 FITEM,2,119 FITEM,2,125 FITEM,2,131 FITEM,2,141 AL,P51X
!划分网格生成有限元模型!划分梁单元
mat,1 !指定梁单元材料特性 type,1 Lsel,s,,,54 lsel,a,,,35 lsel,a,,,45,46 lsel,a,,,10 lsel,a,,,74 lsel,a,,,89
lsel,a,,,102,138,6 lsel,a,,,140 lsel,a,,,145 lsel,a,,,51
lsel,a,,,151,175,6 lsel,a,,,188 lsel,a,,,186 lsel,a,,,193 lsel,a,,,196 lsel,a,,,205 lsel,a,,,208 lsel,a,,,34 lsel,a,,,226 lsel,a,,,23
lsel,a,,,66 lsel,a,,,64 lsel,a,,,71 lsel,a,,,77 lsel,a,,,83 lsel,a,,,86 lsel,a,,,93 lsel,a,,,47 lsel,a,,,105 lsel,a,,,111 lsel,a,,,117 lsel,a,,,123 lsel,a,,,129 lsel,a,,,142 lsel,a,,,136 lsel,a,,,32 lsel,a,,,148 lsel,a,,,154 lsel,a,,,160 lsel,a,,,166 lsel,a,,,172 lsel,a,,,179 lsel,a,,,50 lsel,a,,,191 lsel,a,,,199 lsel,a,,,33 lsel,a,,,211 lsel,a,,,224 lsel,a,,,228 lesize,all,,,1 lmesh,all
mat,4 !指定锚杆单元材料特性 type,3 real,2
lsel,s,,,203 lsel,a,,,219 lsel,a,,,225 lsel,a,,,220 lsel,a,,,53 lsel,a,,,44 lsel,a,,,48 lsel,a,,,69 lsel,a,,,76 lsel,a,,,81 lsel,a,,,97 lsel,a,,,103 lsel,a,,,109 lsel,a,,,115 lsel,a,,,121 lsel,a,,,127 lsel,a,,,133 lsel,a,,,139 lsel,a,,,135 lsel,a,,,147 lsel,a,,,153 lsel,a,,,159 lsel,a,,,165 lsel,a,,,171 lsel,a,,,177 lsel,a,,,183 lsel,a,,,189 lsel,a,,,195 lsel,a,,,201 lsel,a,,,198 lsel,a,,,213 lsel,a,,,210 lsel,a,,,215 lsel,a,,,39 lsel,a,,,52 lsel,a,,,67 lsel,a,,,73 lesize,all,,,3 lmesh,all allsel
!划分开挖区 mat,3 type,2
FLST,5,2,4,ORDE,2 FITEM,5,3 FITEM,5,-4 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !*
LESIZE,_Y1, , ,10, , , , ,1 !*
