这是《FLAC/FLAC3D基础与工程实例》第8.6部分应用实例的的一个修改，原来例子是土体模量随小主应力变化，我自己建的模型：建模----初始应力---开挖隧道----计算----根据小主应力大小改变模量的值，但是会报错：property ‘young’ not found in zone 1，zone 1就是开挖后的null单元，如果直接建有隧道模型就不会报错，请问这个问题出错的原在哪里，该如何修改下，在此先谢谢了。命令：
new
; 圆柱体外环绕放射网格1
ge zo radc size 5 1 15 10 &
p0=0 0 0 p1=15 0 0 p2=0 1 0 p3=0 0 15 &
di 2 2 2 2 fill gro tunnel
;圆柱体外环绕放射网格2
ge zo radc size 5 1 15 10 &
p0=0 0 0 p1=0 0 -15 p2=0 1 0 p3=15 0 0 &
di 2 2 2 2 fill gro tunnel
;镜像
ge zo ref nor -1 0 0 or 0 0 0
mo el
prop bu 3e8 sh 1e8
ini dens 2000
set grav 0 0 -10
fix z ran z -15.01 -14.99
fix x ran x -15.01 -14.99
fix x ran x 14.99 15.01
fix y
solve
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
mo nu ra gro tunnel
solve
save modify.sav
rest modify.sav
def E_modify
p_z=zone_head
d_k=704
d_n=0.38
d_pa=101325.0
loop while p_z # null
sigma_3=-1.0*z_sig1(p_z)
E_new=d_k*d_pa*(sigma_3/d_pa)^d_n
z_prop(p_z,'young')=E_new
p_z=z_next(p_z)
endloop
end
E_modify
