fish
在运行土体的模量 最小主应力变化 程序中
gen zone b size 3 3 3
model elas
prop bu 3e7 sh 1e7
ini dens 2000
fix x y z ran z -.1 .1
fix x ran x -.1 .1
fix x ran x 2.9 3.1
fix y ran y -.1 .1
fix y ran y 2.9 3.1
set grav 0 0 10
solve
ini xd 0 yd 0 xv 0 yv 0 zv 0
save 8-4.sav
rest 8-4.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_sigl(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
怎么会出现unbalance quot marks 的提示那？
