flac3d5.0有几个关于单元弹性应变能和塑性耗散能的命令: z_wevol(ptr), z_weshear(ptr),z_wetot(ptr),z_wpvol(ptr),z_wpshear( pz ),我在调用的时候输出结果都是0,不知道是怎么回事,求大神指点。命令流如下
new
gen zone brick size 20 20 20 p0 0 0 0 p1 4 0 0 p2 0 4 0 p3 0 0 4
model mohr
set large
pro den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10
fix x y z range y -.1 .1
fix x y z range y 3.9 4.1
ini yvel 2.5e-5 range y -.1 .1
ini yvel -2.5e-5 range y 3.9 4.1
step 5000
def get_eng
ad4 = z_near(1,1,1)
p_z = ad4
z_v = z_volume(ad4)
z_v_d = z_vol_deformed(ad4)
w_el_tot = z_wetot(ad4)
w_el_vol = z_wevol(ad4)
w_pl_tot = z_wptot(ad4)
end
@get_eng
list fish