论坛 › G0201：FLAC3D/FLAC2D › 下面是对已经破坏的塑性区的参数重新赋值结果不对请指教

azhekill 发表于 2011-11-26 21:44:08

下面是对已经破坏的塑性区的参数重新赋值结果不对请指教

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

azhekill 发表于 2011-11-26 21:46:07

苦闷中请高手指教

azhekill 发表于 2011-11-28 15:00:24

命令流应该是有问题   但不知道怎么改

tollyee 发表于 2011-11-30 13:06:00

有 and 这种用法吗？

azhekill 发表于 2011-12-6 11:17:09

:(:'(:'(:$有

azhekill 发表于 2011-12-6 11:19:56

但 很少:o:(:(:D

cid2012 发表于 2011-12-6 11:21:28

:)太复杂了 先换一个参数 一种状态 先把这个搞定了。

azhekill 发表于 2011-12-6 11:31:32

已经解决 fish调用的问题这个函数没什么问题   谢谢了

haozhenggong 发表于 2011-12-28 17:06:11

我也遇到类似原因，你怎么解决的啊？？

azhekill 发表于 2012-1-5 10:23:20

:o:(:D就行啦   有些特定的程序不一定自己非得编出来

佟彬1103 发表于 2013-4-16 05:32:00

顶一下，很有用。

shouhuwenming 发表于 2013-8-5 15:18:42

这个命令流是对的

lf_1989 发表于 2014-8-1 11:23:50

顶一下，很有用。

查看完整版本: 下面是对已经破坏的塑性区的参数重新赋值结果不对请指教