APDL如下：
finish
/clear
keyw,pr_struc,1 !选择结构分析
/prep7
ET,1,PLANE42
!第一步开挖
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,3
KEYOPT,1,5,0
KEYOPT,1,6,0
!第二步开挖
ET,2,PLANE42
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,3
KEYOPT,2,5,0
KEYOPT,2,6,0
!第三步开挖
ET,3,PLANE42
KEYOPT,3,1,0
KEYOPT,3,2,0
KEYOPT,3,3,3
KEYOPT,3,5,0
KEYOPT,3,6,0
!外部围岩
ET,4,PLANE42
KEYOPT,4,1,0
KEYOPT,4,2,0
KEYOPT,4,3,3
KEYOPT,4,5,0
KEYOPT,4,6,0
R,1,4
R,2,4
R,3,4
R,4,4
ET,5,LINK10 !第一道横撑
keyopt,5,3,1
!*
ET,6,LINK10 !第二道横撑
keyopt,6,3,1
ET,7,LINK10 !第三道横撑
keyopt,7,3,1
ET,8,BEAM3 !围护
MP,EX,1,5e10
MP,PRXY,1,0.3
MP,DENS,1,18e2
MP,EX,2,30e11
MP,PRXY,2,0.25
MP,DENS,2,25e2
mp,ex,3,26e17
mp,prxy,3,0.3
mp,dens,3,78e2
R,8,4*1,4/12,1,,,,
R,5,3.14*(0.6*0.6-0.5*0.5)/4,-200e3/(26e11*(3.14*(0.6*0.6-0.5*0.5)/4)),
R,6,3.14*(0.6*0.6-0.5*0.5)/4,-200e3/(26e11*(3.14*(0.6*0.6-0.5*0.5)/4)),
R,7,3.14*(0.6*0.6-0.5*0.5)/4,-200e3/(26e11*(3.14*(0.6*0.6-0.5*0.5)/4)),
SECTYPE,5,BEAM,CTUBE,,0
SECOFFSET,CENT
SECDATA,0.5,0.6,0,0,0,0,0,0,0,0
SECTYPE,6,BEAM,CTUBE,,0
SECOFFSET,CENT
SECDATA,0.5,0.6,0,0,0,0,0,0,0,0
SECTYPE,7,BEAM,CTUBE,,0
SECOFFSET,CENT
SECDATA,0.5,0.6,0,0,0,0,0,0,0,0
SECTYPE,8,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,4,1,0,0,0,0,0,0,0,0
k,1,-40,-40
k,2,-40,-20
k,3,-40,-16
k,4,-40,-12
k,5,-40,-6
k,6,-40,0
allsel,all
KSEL,s,,,1,6,,
KGEN,2,all,,,30,,,,0
allsel,all
KSEL,s,,,1,6,,
KGEN,2,all,,,50,,,,0
allsel,all
KSEL,s,,,1,6,,
KGEN,2,all,,,80,,,,0
allsel,all
*do,i,1,5,1
a,i,i+1,i+7,i+6
*enddo
*do,i,7,11,1
a,i,i+1,i+7,i+6
*enddo
*do,i,13,17,1
a,i,i+1,i+7,i+6
*enddo
lplot
/PSYMB,LDIR,1
LREVERSE,16,1
LREVERSE,13,1
LREVERSE,10,1
LREVERSE,7,1
/PSYMB,LDIR,0
allsel,all
lsel,s,loc,x,-40
lsel,a,loc,x,-10
lsel,a,loc,x,10
lsel,a,loc,x,40
lsel,r,loc,y,-20,0
lesize,all,1
allsel,all
lsel,s,loc,x,-40
lsel,a,loc,x,-10
lsel,a,loc,x,10
lsel,a,loc,x,40
lsel,r,loc,y,-40,-20
lesize,all,5
allsel,all
lsel,s,loc,y,0
lsel,a,loc,y,-6
lsel,a,loc,y,-12
lsel,a,loc,y,-16
lsel,a,loc,y,-20
lsel,a,loc,y,-40
lesize,all,5
allsel,all
asel,s,,,1,5
asel,a,,,6,7
asel,a,,,11,15
AATT,1,4,4,
amesh,all
allsel,all
asel,s,,,8
AATT,1,3,3,
amesh,all
allsel,all
asel,s,,,9
AATT,1,2,2,
amesh,all
allsel,all
asel,s,,,10
AATT,1,1,1,
amesh,all
allsel,all
lsel,s,loc,x,-10
lsel,a,loc,x,10
lsel,r,loc,y,-20,0
LATT,2,8,8,,,,8
lmesh,all
allsel,all
lsel,s,loc,y,0
lsel,r,loc,x,-10,10
LATT,3,5,5,,,,5
lmesh,all
allsel,all
lsel,s,loc,y,-6
lsel,r,loc,x,-10,10
LATT,3,6,6,,,,6
lmesh,all
allsel,all
lsel,s,loc,y,-12
lsel,r,loc,x,-10,10
LATT,3,7,7,,,,7
lmesh,all
ALLSEL,ALL
NSEL,S,LOC,x,-40
NSEL,A,LOC,x,40
d,all,ux,0
ALLSEL,ALL
NSEL,S,LOC,y,-40
d,all,uy,0
ALLSEL,ALL
acel,,10
/SOLU
NLGEOM,1
NROPT,FULL,,ON
LUMPM,1
EQSLV,,,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
SSTIF
PSTRES
TOFFST,0,
NSUBST,5,5,1
time,1 !初始应力状态
allsel,all !杀死围护、横撑
esel,s,type,,8
esel,a,type,,7
esel,a,type,,6
esel,a,type,,5
ekill,all
esel,all
esel,s,live
allsel,all
solve
save
time,2 !第一步开挖
allsel,all
esel,s,type,,8
esel,a,type,,5
ealive,all
allsel,all
esel,s,type,,1
ekill,all
esel,all
esel,s,live
allsel,all
solve
save
time,3 !第二步开挖
allsel,all
esel,s,type,,6
ealive,all
allsel,all
esel,s,type,,2
ekill,all
esel,all
esel,s,live
allsel,all
solve
save
time,4 !第三步开挖
allsel,all
esel,s,type,,7
ealive,all
allsel,all
esel,s,type,,3
ekill,all
esel,all
esel,s,live
allsel,all
solve
save
算到time3时就不熟收敛了，有什么需要修正的么？
