new
gene zone bri size 3 3 3
model elastic
prop bu 3e7 sh 1e7
ini dens 2000
fix x range x -1 1
fix x range x 2.9 3.1
fix y range y -1 1
fix y range y 2.9 3.1
fix z range z -1 1
set grav 10
solve
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
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
红体字一行事什么意思?young代表什么?、