interface面的破裂状态使用哪个命令啊？急急急。。。

guo1yan2

请问：z_state（p_z, 0）可以获得区的破裂状态，但interface面的破裂状态使用哪个命令啊？急急急。。。
怎样用f内部函数来表示 iinterface shear slip的破裂状态啊，

cloud328

plot interface slip shear/normal
shear 表示剪切滑动
normal表示法向是否接触脱开

qq511639960

亲，您好，请问您的问题解决了吗？最近我也遇到了这个问题，想用FISH语言控制interface的破裂状态。还请您不吝赐教！！！

qq511639960

想了好久终于自己想明白了，我是用下面这个函数来对模型中单个interface进入塑型区之后进行变参的，4.92Mpa是我初始的c，当剪应力达到11.5MPa时，进行变参
现在贴出来给大家做个参考
有错误之处还请给位大神点出来！
发个牢骚，感觉现在逛论坛的人明显比05到10年之间少了很多呀
def iget_plast
    while_stepping
    p_i=i_head
    p_in=i_node_head(p_i)
    loop while p_in # null
      if in_prop(p_in,'coh') = 4920000 then
            local ipd = in_sstr(p_in,1)         
            if ipd > 11500000 then
              in_prop(p_in,'fric') = 13
              in_prop(p_in,'coh') = 0
            endif
      endif
      p_in= in_next(p_in)
     end_loop
end

qq511639960

想了好久终于自己想明白了，我是用下面这个函数来对模型中单个interface进入塑型区之后进行变参的，4.92Mpa是我初始的c，当剪应力达到11.5MPa时，进行变参
现在贴出来给大家做个参考
有错误之处还请给位大神点出来！
发个牢骚，感觉现在逛论坛的人明显比05到10年之间少了很多呀
def iget_plast
    while_stepping
    p_i=i_head
    p_in=i_node_head(p_i)
    loop while p_in # null
      if in_prop(p_in,'coh') = 4920000 then
            local ipd = in_sstr(p_in,1)         
            if ipd > 11500000 then
              in_prop(p_in,'fric') = 13
              in_prop(p_in,'coh') = 0
            endif
      endif
      p_in= in_next(p_in)
     end_loop
end

qq511639960

对上面的进行了一下改进，1.41422是我的节理角度为45°，ipd = in_sstr(p_in,1)里的1，2，3分别指的是整体坐标系的xyz轴。
希望能对大家有所帮助！！！！！！
def iget_plast
    while_stepping
    p_i=i_head
    p_in=i_node_head(p_i)
    loop while p_in # null
      if in_prop(p_in,'coh') = 4820000 then
               ipd = in_sstr(p_in,1)
                        ipdd=ipd*1.41422
                        innstr=in_nstr(p_in)
                        incoh=in_prop(p_in,'coh')
                        infric=in_prop(p_in,'fric')
                        infric1=infric*pi/180
                        infric2=tan(infric1)
                pohuai=innstr*infric2+incoh
            if ipdd > pohuai then
              in_prop(p_in,'fric') = 1
              in_prop(p_in,'coh') = 0
              in_prop(p_in,'ten') = 0
            endif
      endif
      p_in= in_next(p_in)
     end_loop
end






不好意思，网络卡了发重了，版主看到可以帮忙删一下

cugyi

qq511639960 发表于 2019-10-8 18:23
对上面的进行了一下改进，1.41422是我的节理角度为45°，ipd = in_sstr(p_in,1)里的1，2，3分别指的是整体 ...

层主，你有没有不规则结构面直剪的例子