- 积分
- 0
- 注册时间
- 2012-9-28
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2014-7-18 11:10:42
|
显示全部楼层
来自 湖北武汉
我知道一个计算塑性区体积的命令流,至于计算面积还没有查。另外,无论是面还是体积,我觉得都可以在顶、底、两帮进行不同分组,然后分别调用命令文件输出应该可以实现。
- ;;;;;;;;;;;;;;;监测塑性区体积
- def suxingqu
- shearnow=1
- tensionnow=2
- shearpast=4
- tensionpast=8
- v_shear_now=0
- v_tension_now=0
- v_shear_past=0
- v_tension_past=0
- p_z=zone_head
- loop while p_z # null
- if and(z_state(p_z,0),shearnow)=shearnow then
- v_shear_now=v_tension_now+z_volume(p_z)
- endif
- if and(z_state(p_z,0)),shearpast)=shearpast then
- v_shear_past=v_shear_past+z_volume(p_z)
- endif
- if and(z_state(p_z,0),tensionpast)=tensionpast then
- v_tension_past=v_tension_past=z_volume(p_z)
- endif
- p_z=z_next(p_z)
- endloop
- ii=out('剪切破坏:’+string(v_shear_now))
- ii=out('拉伸破坏:’+string(v_tension_now))
- ii=out('过去曾剪切破坏:’+string(v_shear_past))
- ii=out('过去曾拉伸破坏:’+string(v_tension_past))
- end
复制代码 |
|