- 积分
- 237
- 注册时间
- 2003-6-27
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2019-5-26 10:13:09
|
显示全部楼层
来自 上海
我不熟悉 workbench,但是你这个问题可以在 workbench 中插入 APDL 命令来实现,大致思路如下:
在 /POST1 中,读取所需的结果组后:
ALLSEL, all ! 选择所有实体
*get, emax, ELEM, 0, NUM, MAX ! 获取最大单元编号
*dim, e_seqv, array, max, 1, 1 ! 定义一个数组,春芳个定义的 misses 应力值
*do, i, 1, emax ! 对定义进行循环
*if, esel(i),gt,0,then ! 如果单元 i 存在
*get, seqv_i, ELEM, i, s, EQV ! 获取单元 i ,重心处的 misses 应力
e_seqv(i) = seqv_i
*else
e_seqv(i) = 0
*endif
*enddo
e_n = 0 ! 准备存放满足应力条件的单元个数
! 假设 misses 应力范围为: s_min <= misses 应力 <= s_max
*do, i, 1, emax
seqvi = e_seqv(i)
*if, seqvi, ge, s_min, and, seqvi, le, s_max, then
e_n = e_n + 1
*endif
*enddo
至此, e_n 为mauzu条件的单元个数;数组 e_seqv 中为各个单元重心出的 Mises 应力。
没有验证过,仅供参考。请自己进行验证和修改完善。
|
|