没考虑你说的土体损失率，土体损失率是啥？我看了相关文章，他们使用数值模拟得到的结果是埋深越大，沉降越小。不知道我的为什么不行。代码如下
new 开始新的分析
title;
围岩稳定性参数敏感性分析001
expgrid xiaozhu.flac3d
plot zone
group tunnel range group 1
group surock range group 2
plot zone wireframe width 2 color blue
define para_001
young_e=240e6
p_ratio=0.3
end
@para_001
define para_002
bulk_e=young_e/(3*(1-2*p_ratio))
shear_e=young_e/(2*(1+p_ratio))
cohesion_001=30e3
friction_001=32
tension_001=20e3
density_001=2.1e3
end
@para_002
;指定岩土体本构模型
model mohr
;指定岩土体参数
prop bulk @bulk_e shear @shear_e density @density_001 &
cohesion @cohesion_001 friction @friction_001 &
tension @tension_001 dilation 0 range group tunnel
prop bulk @bulk_e shear @shear_e density @density_001 &
cohesion @cohesion_001 friction @friction_001 &
tension @tension_001 dilation 0 range group surock
;设定位移边界条件
fix x range x -22.51 -22.49
fix x range x 22.49 22.51
fix y range y -0.01 0.01
fix y range y 9.99 10.01
fix z range z -22.51 -22.49
set gravity 0 0 -10
solve elas
save 自重应力场001.sav
plot zcon szz
plot con dis
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
hist id 1 gp zdis 0 5 2.5
hist id 2 gp xdis -2.5 5 0
hist id 3 gp xdis 2.5 5 0
hist id 4 gp zdis 0 5 -2.5
efine excav
loop n(1,5)
y1=2*(n-1)
y2=2*n
filename='循环开挖'+string(n)+'.sav'
command
model null range group tunnel y @y1 @y2
;model null range reg tunnel y y1 y2 (3.0)
sel shell id 1 group surock range cyl end1 0 @y1 0 end2 0 @y2 0 radius 2.5
sel shell id 1 prop isotropic 10.5e9 0.25 thickness 0.25 density 2500 range y @y1 @y2
solve
endcommand
command
save @filename
endcommand
endloop
end
@excav
;写入监测数据
hist write 1 2 3 4 file 监测数据0510hist.txt
table 1 name "沉降值"
def cross_curve
loop ii(1,21)
xx=2.25*(ii-11)
yy=5
zz=12.5
p_gp=gp_near(xx,yy,zz)
xtable(1,ii)=xx ;xtable(表编号，行号)=xx x代表x列
ytable(1,ii)=gp_zdisp(p_gp)
;ytable(表编号，行号)=gp_zdisp(p_gp) y代表y列
endloop
end
@cross_curve
set log on
set logfile 地面不同点沉降位移监测.log
print table 1
set log off
