我用flac3d计算的隧道施工 拱顶围岩位移为什么和拱顶初衬位移不相等啊 应该两者位移一模一样啊 怎么会不相等啊 请问大家 有谁知道?怎么解决啊
代码如下
rest geo.sav
plo blo gro
;自重应力场求解
set gravity 0 0 -9.81
;施加位移边界
fix z range z -30.1 -29.9
fix x range x 39.9 40.1
fix x range x -40.1 -39.9
fix y range y -0.01 0.01
fix y range y 59.9 60.1
;赋予模型参数
def wycs
e_mod=0.1e9
P_ratio=0.4237
cc=0.1342e6
ff=25.4629
end
wycs
def huansuan
s_mod=e_mod/(2*(1+P_ratio))
b_mod=e_mod/(3*(1-2*p_ratio))
end
huansuan
model mohr ;莫尔~库仑模型
ini density 2000 ;围岩的密度
prop bulk=b_mod shear=s_mod fric=ff coh=cc tension=cc
step 2000
;保存初始应力场结果
save gsol.sav
;绘制云图
plot blo sta
initial xdisp=0 ydisp=0 zdisp=0
;开挖1
model null range group 3 y 0 2
model null range group 5 y 0 2 z 0 5.4
step 1200
save step1ks.save
;支护1
sel shell id 1 group 1 range cyl end1 0 0 0 end2 0 2 0 rad 5.44 cyl end1 0 0 0 end2 0 2 0 rad 5.38 not z 0 5.4 x -5.4 5.4
sel shell prop isotropic 21.0e9,0.2 thickness 0.25 density 2500
;fix z range z -0.05 0.05 x -5.41 -5.39 y 0 2
;fix z range z -0.05 0.05 x 5.39 5.41 y 0 2
step 500
save step1zh.sav
;开挖2
model null range group 3 y 2 4
model null range group 5 y 2 4 z 0 5.4
step 1200
save step2ks.save
;支护2
sel shell id 1 group 1 range cyl end1 0 2 0 end2 0 4 0 rad 5.44 cyl end1 0 2 0 end2 0 4 0 rad 5.38 not z 0 5.4 x -5.4 5.4
sel shell prop isotropic 21.0e9,0.2 thickness 0.25 density 2500
;fix z range z -0.05 0.05 x -5.41 -5.39 y 2 4
;fix z range z -0.05 0.05 x 5.39 5.41 y 2 4
step 500
save step2zh.sav
;开挖3
model null range group 3 y 4 6
model null range group 5 y 4 6 z 0 5.4
step 1200
save step3ks.save
;支护3
sel shell id 1 group 1 range cyl end1 0 4 0 end2 0 6 0 rad 5.44 cyl end1 0 4 0 end2 0 6 0 rad 5.38 not z 0 5.4 x -5.4 5.4
sel shell prop isotropic 21.0e9,0.2 thickness 0.25 density 2500
;fix z range z -0.05 0.05 x -5.41 -5.39 y 4 6
;fix z range z -0.05 0.05 x 5.39 5.41 y 4 6
step 500
save step3zh.sav
;开挖4
model null range group 3 y 6 8
model null range group 5 y 6 8 z 0 5.4
step 1200
save step4ks.save
;支护4
sel shell id 1 group 1 range cyl end1 0 6 0 end2 0 8 0 rad 5.44 cyl end1 0 6 0 end2 0 8 0 rad 5.38 not z 0 5.4 x -5.4 5.4
sel shell prop isotropic 21.0e9,0.2 thickness 0.25 density 2500
;fix z range z -0.05 0.05 x -5.41 -5.39 y 6 8
;fix z range z -0.05 0.05 x 5.39 5.41 y 6 8
step 500
save step4zh.sav
;开挖5
model null range group 3 y 8 10
model null range group 5 y 8 10 z 0 5.4
step 1200
save step5ks.save
;支护5
sel shell id 1 group 1 range cyl end1 0 8 0 end2 0 10 0 rad 5.44 cyl end1 0 8 0 end2 0 10 0 rad 5.38 not z 0 5.4 x -5.4 5.4
sel shell prop isotropic 21.0e9,0.2 thickness 0.25 density 2500
;fix z range z -0.05 0.05 x -5.41 -5.39 y 8 10
;fix z range z -0.05 0.05 x 5.39 5.41 y 8 10
step 500
save step5zh.sav