完整命令流如下:
gen zon bri size 3 3 3
model elastic
prop bu 3e7 sh 1e7
ini dens 2000
fix x y z ran z -.1 .1
fix 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 10
solve
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
def E_new
p_z=zone_head
miu=0.3
loop while p_z # null
bulm=bulk
iE_new=3*bulm*(1-2*miu)
z_prop(p_z,'young')=iE_new
p_z=z_next(p_z)
endloop
end
E_new 想用fish语言调用计算过程中单元的体积模量,进而求出各单元的弹性模量,但是就是调用不了,知道怎么解决的同仁帮帮忙,多谢