FINISN
/CLEAR
/PREP7
!创建几何模型
K,1,0,2.7388794
K,2,-2.0711036,1.7033277 K,3,-2.6388794,0
K,4,-2.3359835,-1.9455903 K,5,-1.6736199,-2.6205221 K,6,0,-2.8972106
K,1000,,,
LARC,1,2,1000,2.589 LARC,2,3,1000,2.839 LARC,3,4,1000,6.4 LARC,4,5,1000,1.05 LARC,5,6,1000,5.2 LSYMM,X,ALL NUMMRG,ALL NUMCMP,ALL LCOMB,1,2
LCOMB,6,7 NUMMRG,ALL NUMCMP,ALL LCOMB,3,4
LCOMB,7,8 NUMMRG,ALL NUMCMP,ALL LCOMB,2,3
LCOMB,5,6 NUMMRG,ALL NUMCMP,ALL
AL,ALL
BLC4,-36,-69,72,84 ASBA,2,1,,,KEEP NUMCMP,ALL
!切分
ASEL,S,,,2 WPROTA,,,90
ASBW,ALL WPOFFS,,,6
ASBW,ALL WPOFFS,,,-12
ASBW,ALL WPCSYS,-1
ASBW,ALL
WPOFFS,,,6
ASBW,ALL
WPOFFS,,,-12
ASBW,ALL
WPCSYS,-1
ALLS
NUMMRG,ALL
NUMCMP,ALL
!定义材料属性
*SET,TH1,0.3
ET,1,BEAM3
ET,2,PLANE42
KEYOPT,2,3,2
R,1,TH1,TH1*TH1*TH1/12,TH1,,
MP,EX,1,31E9
MP,PRXY,1,0.3
MP,DENS,1,2500
MP,EX,2,1.40E8!IV级围岩物理力学参数MP,PRXY,2,0.3
MP,DENS,2,2049
TB,DP,2
TBDATA,1,3.4E4,30
!赋予属性
LSEL,S,,,1,4
LATT,1,1,1
ASEL,ALL
AATT,2,,2
!隧道周围划分网格
LSEL,S,,,1,3,2
LESIZE,ALL,,,8
LSEL,S,,,2,4,2
LESIZE,ALL,,,8
LSEL,S,,,1,4
LMESH,ALL
LSEL,S,,,24,25
LSEL,A,,,35
LSEL,A,,,44
LESIZE,ALL,,,8,3
LSEL,S,,,29,30
LESIZE,ALL,,,4
LSEL,A,,,32,34,2
LSEL,S,,,38,39 LSEL,A,,,41,43,2 LESIZE,ALL,,,4 ASEL,S,,,1 AMESH,ALL ALLS
AMAP,16,1,29,19,2 AMAP,17,4,17,29,1 AMAP,12,3,24,19,2 AMAP,13,4,17,24,3 LSEL,S,,,15 LSEL,A,,,18 LESIZE,ALL,,,4 LSEL,S,,,16 LSEL,A,,,19 LESIZE,ALL,,,4 LSEL,S,,,12 LSEL,A,,,13 LSEL,A,,,23 LSEL,A,,,27 LSEL,A,,,36 LESIZE,ALL,,,20 LSEL,S,,,5 LSEL,A,,,7 LSEL,A,,,11 LSEL,A,,,14 LSEL,A,,,20 LESIZE,ALL,,,20 LSEL,S,,,9,10 LSEL,A,,,37 LSEL,A,,,21,22 LESIZE,ALL,,,6 LSEL,S,,,42 LSEL,A,,,40 LESIZE,ALL,,,4 LSEL,S,,,33 LSEL,A,,,31 LESIZE,ALL,,,4 LSEL,S,,,17 LSEL,A,,,26 LSEL,A,,,28 LSEL,A,,,8 LSEL,A,,,6
ALLS
ASEL,U,,,1
ASEL,U,,,12,13
ASEL,U,,,16,17
AMESH,ALL
ALLS
ACEL,,9.8
ASEL,S,,,1
ESLA,S
EKILL,ALL
ALLS
!2,施加边界条件
NSEL,S,LOC,X,-36
D,ALL,UY,0
NSEL,S,LOC,X,36
D,ALL,UY,0
NSEL,S,LOC,Y,-69
D,ALL,ALL
ALLS
!3,读入加速度时程
!施加水平地震荷载，命令流如下：
!先把SHUIPING.TXT的EI-CENTRO波水平方向加速度时程文件放在E:\盘中即工作目录下ALLS
*DIM,SHUIPING,ARRAY,2,1000,0,,,
*CREATE,ANSUITMP
*VREAD,SHUIPING(1,1),'SHUIPING','txt','E:\',1000,,,,
(E9.3,E11.3)
*END
/INPUT,ANSUITMP
!6.7.5加载与求解
!模型建立的施加水平加速度荷载完毕后，进行计算求解分析，命令流如下：
/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0
*DO,T,1,75,1
TIME,0.02*T
KBC,1
ACEL,SHUIPING(2,T)
ALPHAD,0.069
BETAD,0.035
SOLVE
*ENDDO
*DO,T,76,150,1
TIME,0.02*T
KBC,1
NSUB,50
ACEL,SHUIPING(2,T)
ALPHAD,0.069
BETAD,0.035
SOLVE
*ENDDO
*DO,T,151,750,1
TIME,0.02*T
KBC,1
NSUB,600
ACEL,SHUIPING(2,T)
ALPHAD,0.069
BETAD,0.035
SOLVE
*ENDDO
FINISH
!6.7.6结果分析
!地表中点竖向位移，速度及加速度分析/POST26
/AXLAB,Y,UY
NSOL,2,2811,U,Y,UY_2
PRV AR,2
PLV AR,2
/AXLAB,Y,VY
DERIV,3,2,1,,,,,1
PLV AR,3
PRV AR,3
/AXLAB,Y,AY
DERIV,4,3,1,,,,,1
PLV AR,4
PRV AR,4
!地表中点竖向位移，速度及加速度分析/POST26
/AXLAB,Y,UX
NSOL,2,2811,U,X,UX_2
PLV AR,2
/AXLAB,Y,VX DERIV,3,2,1,,,,,1 PLV AR,3
PRV AR,3
/AXLAB,Y,AX DERIV,4,3,1,,,,,1 PLV AR,4
PRV AR,4。