- 积分
- 0
- 注册时间
- 2011-8-12
- 仿真币
-
- 最后登录
- 1970-1-1
|
def get_plast
shearnow = 1
tensionnow = 2
shearpast = 4
tensionpast = 8
p_z = zone_head
loop while p_z # null
if and(z_state(p_z,0),shearnow) = shearnow then
z_prop(p_z,'bulk')=1.7*15.98148e9
z_prop(p_z,'shear')=1.7*6.537878e9
z_prop(p_z,'cohesion')=1.7*21.82e6
z_prop(p_z,'tension')=0.7*8e6
endif
if and(z_state(p_z,0),tensionnow) = tensionnow then
z_prop(p_z,'bulk')=1.3*15.98148e9
z_prop(p_z,'shear')=1.3*6.537878e9
z_prop(p_z,'cohesion')=1.3*21.82e6
z_prop(p_z,'tension')=0.3*8e6
endif
if and(z_state(p_z,0),shearpast) = shearpast then
z_prop(p_z,'bulk')=1.5*15.98148e9
z_prop(p_z,'shear')=1.5*6.537878e9
z_prop(p_z,'cohesion')=1.5*21.82e6
z_prop(p_z,'tension')=0.5*8e6
endif
if and(z_state(p_z,0),tensionpast) = tensionpast then
z_prop(p_z,'bulk')=1.3*15.98148e9
z_prop(p_z,'shear')=1.3*6.537878e9
z_prop(p_z,'cohesion')=1.3*21.82e6
z_prop(p_z,'tension')=0.3*8e6
endif
p_z = z_next(p_z)
endloop
end
get_plast
|
|