;grid xzones yzones
gen zone cylinder p0 0 0 0 p1 1 0 0 p2 0 4 0 p3 0 0 1 size 2 2 3
gen zone cylinder p0 0 0 0 p1 0 0 1 p2 0 4 0 p3 -1 0 0 size 2 2 3
gen zone cylinder p0 0 0 0 p1 -1 0 0 p2 0 4 0 p3 0 0 -1 size 2 2 3
gen zone cylinder p0 0 0 0 p1 0 0 -1 p2 0 4 0 p3 1 0 0 size 2 2 3
model mech cam-clay
;prop shear @m_g0 bulk_bound @m_kc
prop mm @mcsF lambda @lamcsF kappa @kapcsF;mm:frictional constant,M
prop mpc @m_pc mp1 1. mv_l 3.32
; --- boundary and initial conditions ---
fix x y
;ini xvel -0.5e-4 range x 0.9 1.1
;ini yvel -0.5e-4 range y 0.9 1.1
;ini zvel -0.5e-4 range z 0.9 1.1
def camclay_inip
local p_z = zone_head
loop while p_z # null
m_p=-(z_sxx(p_z)+z_syy(p_z)+z_szz(p_z))/3
z_prop(p_z,'cam_cp') = m_p
m_q1=1.5*((z_sxx(p_z)-m_p)^2+(z_syy(p_z)-m_p)^2+(z_szz(p_z)-m_p)^2)
m_q=sqrt(m_q1)
z_prop(p_z,'cq')=m_q
ev=z_volume(p_z)
z_prop(p_z,'cam_ev')=ev
p_z=z_next(p_z)
end_loop
end
为什么我的程序输出以后，m_q和m_p均为0呢？不知道有没有大神知道的，能不能告诉我一下，小弟拜谢~
