如何提取mises应力

gnik77

请教各位，如何提取一点的mises应力阿？谢谢

范金泉

在Flac里面有

ycding

如下：
config zextra 1
def mises
; --- calculate and store Von Mises stress in zone extension 1 ---
   p_z = zone_head
   max_mises = 0.0
   loop while p_z # null
           temp = (z_sig1(p_z)-z_sig2(p_z))^2+(z_sig2(p_z)-z_sig3(p_z))^2+(z_sig3(p_z)-z_sig1(p_z))^2
           vmstr2 = sqrt(0.5*temp)
      if vmstr2 > 0.0 then
         z_extra(p_z,1) = vmstr2
      else
         z_extra(p_z,1) = 0.0
      endif
      max_mises = max(max_mises,z_extra(p_z,1))
   p_z = z_next(p_z)
   end_loop
end
mises
plot cont zextra 1 alias 'Von Mises Stress' average shade on outl on
print max_mises

cloud328

像三楼一样，编个fish就可以了。

gnik77

3# ycding
万分感谢！又学到一些东西！真的很谢谢你！

gnik77

3# ycding
if vmstr2>0.0 then
这个好像不对吧？vmstr2是开根号出来的，肯定不会小于0阿

gnik77

3# ycding
if vmstr2>0.0 then
这个好像不对吧？vmstr2是开根号出来的，肯定不会小于0阿
我运行了一下，是可以的！
我想具体知道其中一个单元的mises应力的大小，应该怎么提取阿？
初学，所以不会，还望赐教！谢谢您啦！