上图是地表的沉降图
我定义里平衡，工法都是对称的，约束我可以保证假的也对。我不知道为什么模型的数值位移出现这种情况。这是我的地应力平衡和开挖命令。
;;;;;;;;;;;位移为零认为开挖前的位移为0 ，便于分析开挖后的位移
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
ini state 0
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;四号线隧道台阶法开挖
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;变量初值
def bianliangchuzhi
bushu=1
saveanme='开挖'
end
bianliangchuzhi
def jisuan
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;开挖计数
shibeimishu=16*bushu-0
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;开挖
zx_kaiwa0=-1.6+1.6*bushu
zx_kaiwa1=0+1.6*bushu
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;盾构管片
zx_guanpian0=-1.6+1.6*bushu
zx_guanpian1=0+1.6*bushu
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;注浆
zx_zhujiang0=-4.8+1.6*bushu
zx_zhujiang1=-3.2+1.6*bushu
;;;;;;;;;;;;;;;;;;;;;;;;;;;文件名
savename=string(shibeimishu)+'.sav'
end
jisuan
;开挖
def kaiwa
loop bushu(1,55)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;开挖施支
command
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;计算数据
jisuan
;左洞开挖
model null range group zx_kaiwa y zx_kaiwa0 zx_kaiwa1
model null range group zx_guanpian y zx_kaiwa0 zx_kaiwa1
;右洞开挖
model null range group yx_kaiwa y zx_kaiwa0 zx_kaiwa1
model null range group yx_guanpian y zx_kaiwa0 zx_kaiwa1
step 10
;左洞支护管棚1.1
model mohr range gr zx_zhujiang y zx_guanpian0 zx_guanpian1
pro bulk 5.6e7 shear 4.05e7 fric=33 coh=5.67e5 tension=3e4 dens 2340 dil=0 range gr zx_zhujiang y zx_guanpian0 zx_guanpian1
model elas range gr zx_guanpian y zx_guanpian0 zx_guanpian1
pro bu 1.86e10 sh 1.39e10 dens 2500 range gr zx_guanpian y zx_guanpian0 zx_guanpian1
;右洞支护管棚1.1
model mohr range gr yx_zhujiang y zx_guanpian0 zx_guanpian1
pro bulk 5.6e7 shear 4.05e7 fric=33 coh=5.67e5 tension=3e4 dens 2340 dil=0 range gr yx_zhujiang y zx_guanpian0 zx_guanpian1
model elas range gr yx_guanpian y zx_guanpian0 zx_guanpian1
pro bu 1.86e10 sh 1.39e10 dens 2500 range gr yx_guanpian y zx_guanpian0 zx_guanpian1
step 30
;大循环结束
pl cont zd out on
hist unbal
pl add hist 1
set grav 0 0 -10.0
;set mech ratio 5e-5
step 200
save savename
endcommand
endloop
end
kaiwa
我自己试着一次性model null 隧道开挖部分，结果出现的结果是对称的
恳求大神回复谢谢，为什么分布开挖就得到那个结果
