如何提取mises应力
请教各位,如何提取一点的mises应力阿?谢谢 在Flac里面有 如下: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 像三楼一样,编个fish就可以了。 3# ycding
万分感谢!又学到一些东西!真的很谢谢你! 3# ycding
if vmstr2>0.0 then
这个好像不对吧?vmstr2是开根号出来的,肯定不会小于0阿 3# ycding
if vmstr2>0.0 then
这个好像不对吧?vmstr2是开根号出来的,肯定不会小于0阿
我运行了一下,是可以的!
我想具体知道其中一个单元的mises应力的大小,应该怎么提取阿?
初学,所以不会,还望赐教!谢谢您啦!
页:
[1]