使用隧道开挖顺序对比的问题
浏览：159756
一个简单的模型，使用重启动分析模拟两种开挖顺序，一个先左后右，一个先右后左，（都是一下挖通，不考虑支护）两种开挖顺序在第二步开挖时应力变形都是一样的，请问我是哪错了？
第一种：第一步挖左，第二步挖左右
第二种：第一步挖右，第二步挖左右
fini$/cle$/prep7
et,1,45
et,2,63
et,3,200,6
r,1,0.33 !车站站初喷
r,2,0.27 !行车道初喷
r,3,0.24 !匝道初喷
!r,4,0.02,6.7e-5,1.7e-5,0.2,0.1 !20a钢拱架(纵1m)
!r,5,3.14*0.011*0.011 !锚杆(横纵均1m)
mp,ex,1,2.9e9 !围岩
mp,prxy,1,0.25
mp,dens,1,2500
tb,dp,1
tbdata,1,1.3e6,37
mp,ex,2,3.94e8 !填土
mp,prxy,2,0.35
mp,dens,2,2000
tb,dp,2
tbdata,1,2e4,30
mp,ex,3,27.5e9 !初支
mp,prxy,3,0.28
mp,dens,3,2800
mp,ex,4,28e9 !二衬
mp,prxy,4,0.25
mp,dens,4,2500
mp,ex,5,2e11 !锚杆等效
mp,prxy,5,0.2
mp,dens,5,7850
!tb,dp,5
!tbdata,1,2e7,40
mp,ex,6,2.8e10
mp,prxy,6,0.26
mp,dens,6,2600
!tb,dp,6
!tbdata,1,1.8e7,28 !超前支护
mp,ex,7,4e9 !土体加固
mp,prxy,7,0.30
mp,dens,7,2200
tb,dp,7
tbdata,1,2.4e4,35
blc4,-50,-50,100,80
cyl4,,,5,,8
cyl4,15,4,4,,6
alls
aovlap,all
mshkey,2$mshape
aatt,1,,3$esize,5
amesh,all
wpstyl,defa
*do,i,1,4
extopt,esize,1
extopt,aclear,1
mat,1
asel,s,loc,z,(1-i)*5
vext,all,,,,,-5
*enddo
alls
vsel,s,,,1,16,5
vsel,a,,,3,18,5
eslv,s
cm,d1,elem
alls
vsel,s,,,1,16,5
eslv,s
cm,z1,elem
alls
vsel,s,,,2,17,5
vsel,a,,,4,19,5
eslv,s
cm,d2,elem
alls
vsel,s,,,2,17,5
eslv,s
cm,z2,elem
!alls
!vsel,s,loc,y,,100
!eslv,s
!mpchg,2,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/sol
alls
csys
wpstyl,defa
acel,,9.8
asel,s,loc,x,-50
asel,a,loc,x,50
da,all,ux
asel,s,loc,y,-50
da,all,uy
asel,s,loc,z
asel,a,loc,z,-20
da,all,uz
antype,static
time,1
nsubst,10
autots,1 !使用自动时间步
nlgeom,on !打开大位移效果
nropt,full !采用完全NR法
pred,on !打开时间步长预测器
lnsrch,on !打开线性搜索
outres,all,last
alls
solve
fini
/sol
antype,,rest,1
cmsel,s,d1
ekill,all
alls
solve
fini
/sol
antype,,rest
cmsel,s,d1
cmsel,a,d2
ekill,all
alls
solve
/post1
esel,s,live
plnsol,s,1 !查看大主应力
fini
/sol
antype,,rest,1
cmsel,s,d2
ekill,all
alls
solve
fini
/sol
antype,,rest
cmsel,s,d1
cmsel,a,d2
ekill,all
alls
solve
/post1
esel,s,live
plnsol,s,1 !和上面结果一样的
第一种：第一步挖左，第二步挖左右
第二种：第一步挖右，第二步挖左右
fini$/cle$/prep7
et,1,45
et,2,63
et,3,200,6
r,1,0.33 !车站站初喷
r,2,0.27 !行车道初喷
r,3,0.24 !匝道初喷
!r,4,0.02,6.7e-5,1.7e-5,0.2,0.1 !20a钢拱架(纵1m)
!r,5,3.14*0.011*0.011 !锚杆(横纵均1m)
mp,ex,1,2.9e9 !围岩
mp,prxy,1,0.25
mp,dens,1,2500
tb,dp,1
tbdata,1,1.3e6,37
mp,ex,2,3.94e8 !填土
mp,prxy,2,0.35
mp,dens,2,2000
tb,dp,2
tbdata,1,2e4,30
mp,ex,3,27.5e9 !初支
mp,prxy,3,0.28
mp,dens,3,2800
mp,ex,4,28e9 !二衬
mp,prxy,4,0.25
mp,dens,4,2500
mp,ex,5,2e11 !锚杆等效
mp,prxy,5,0.2
mp,dens,5,7850
!tb,dp,5
!tbdata,1,2e7,40
mp,ex,6,2.8e10
mp,prxy,6,0.26
mp,dens,6,2600
!tb,dp,6
!tbdata,1,1.8e7,28 !超前支护
mp,ex,7,4e9 !土体加固
mp,prxy,7,0.30
mp,dens,7,2200
tb,dp,7
tbdata,1,2.4e4,35
blc4,-50,-50,100,80
cyl4,,,5,,8
cyl4,15,4,4,,6
alls
aovlap,all
mshkey,2$mshape
aatt,1,,3$esize,5
amesh,all
wpstyl,defa
*do,i,1,4
extopt,esize,1
extopt,aclear,1
mat,1
asel,s,loc,z,(1-i)*5
vext,all,,,,,-5
*enddo
alls
vsel,s,,,1,16,5
vsel,a,,,3,18,5
eslv,s
cm,d1,elem
alls
vsel,s,,,1,16,5
eslv,s
cm,z1,elem
alls
vsel,s,,,2,17,5
vsel,a,,,4,19,5
eslv,s
cm,d2,elem
alls
vsel,s,,,2,17,5
eslv,s
cm,z2,elem
!alls
!vsel,s,loc,y,,100
!eslv,s
!mpchg,2,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/sol
alls
csys
wpstyl,defa
acel,,9.8
asel,s,loc,x,-50
asel,a,loc,x,50
da,all,ux
asel,s,loc,y,-50
da,all,uy
asel,s,loc,z
asel,a,loc,z,-20
da,all,uz
antype,static
time,1
nsubst,10
autots,1 !使用自动时间步
nlgeom,on !打开大位移效果
nropt,full !采用完全NR法
pred,on !打开时间步长预测器
lnsrch,on !打开线性搜索
outres,all,last
alls
solve
fini
/sol
antype,,rest,1
cmsel,s,d1
ekill,all
alls
solve
fini
/sol
antype,,rest
cmsel,s,d1
cmsel,a,d2
ekill,all
alls
solve
/post1
esel,s,live
plnsol,s,1 !查看大主应力
fini
/sol
antype,,rest,1
cmsel,s,d2
ekill,all
alls
solve
fini
/sol
antype,,rest
cmsel,s,d1
cmsel,a,d2
ekill,all
alls
solve
/post1
esel,s,live
plnsol,s,1 !和上面结果一样的
