找回密码
 注册
Simdroid-非首页
查看: 64|回复: 1

[命令/FISH] FLAC3D中mises应力

[复制链接]
发表于 2016-1-2 11:30:49 | 显示全部楼层 |阅读模式 来自 甘肃
FLAC3D中mises应力的命令流,希望能帮助大家。(注,在网上找到,不是自己的结果。)



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
mstr = (sxx(p_z) + syy(p_z) + szz(p_z)) / 3.
dsxx = sxx(p_z) - mstr
dsyy = syy(p_z) - mstr
dszz = szz(p_z) - mstr
dsxy = sxy(p_z)
dsxz = sxz(p_z)
dsyz = syz(p_z)
vmstr2 = 1.5 * (dsxx*dsxx + dsyy*dsyy + dszz*dszz)
vmstr2 = vmstr2 + 3. * (dsxy*dsxy + dsxz*dsxz + dsyz*dsyz)
if vmstr2 > 0.0 then
z_extra(p_z,1) = sqrt(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
print max_mises


本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2016-3-7 20:20:13 | 显示全部楼层 来自 北京
Simdroid开发平台
mises应力不知道什么
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 19:21 , Processed in 0.031353 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表