Question

FLAC3D隧道开挖得出的结果不平衡？

为什么我试着建立了个很简单的模型，按顺序开挖，地表沉降位移云图呈现先开挖的地方沉降小，后开挖的地表沉降大，是什么原因？

我的开挖命令都是两个洞同时开挖，一共50个循环。我以为是约束不对，但是我试着改变model null 全部开挖，结果又是对称的。

我以为是参数循环错了，结果没有。恳求大神指导。
这是我的开挖命令 ;;;;;;;;位移为零认为开挖前的位移为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





Answer 1

也许你对平衡的概念理解有错，人家默认不平衡力比率是1e-5是平衡，很多模型达不到的，就去改改就行 set mech ratio 1e-4试试 或者干脆直接用step 5000或者多少来替代这个，反正影响不会很大的。对了 推荐一下百度里面想程小杰的Flac3D视频教程