在对结构进行时程分析后,我们经常提取的是全时程最大von Mises stress。
命令流如下:
/POST1
ALLSEL
*DIM,MAX_EQV,,4001
SET,1
PLNSOL,S,EQV,0
*GET,MAX_EQV(1),PLNSOL,0,MAX
MAXEQV=MAX_EQV(1)
NUMSTEP=1
*DO,I,2,4001,1
SET,I
PLNSOL,S,EQV,0
*GET,MAX_EQV(I),PLNSOL,0,MAX
*IF,MAX_EQV(I),GT,MAXEQV,THEN
MAXEQV=MAX_EQV(I)
NUMSTEP=I
*ELSE
MAXEQV=MAXEQV
*ENDIF
*ENDDO
!显示应力云图
SET,NUMSTEP
PLNSOL,S,EQV,0
如需提取某个时间点的最大von Mises stress:
*DIM,MAXSEQVNODE,TABLE,NUMBER,1,1,,, ! 定义数组
*DO,ISET1,1,NUMBER,1 !
/POST1
SET,1,,1,,,,ISET1,
PLNSOL,S,EQV,0
*GET,MAXSEQV_NODEI,PLNSOL,0,MAX !获取von Mises stress
*SET,MAXSEQV_NODE(ISET-1,0,1),(ISET-1)*1E-1! 对应时间
*SET,MAXSEQV_NODE(ISET-1,1,1),MAXSEQV_NODEI
*ENDDO
那么如何提取某一个节点的von Mises stress呢?
我们可以通过选择单元上节点的方法提取节点应力。
1 确定节点所在单元,显示节点编号,例单元号8560,节点号8678。
![9efe44cf623b8a634d324385341378b3.png](https://img-blog.csdnimg.cn/img_convert/9efe44cf623b8a634d324385341378b3.png)
2 进入TimeHist Postpro, 定义变量。
![dc6e5d07ed8454133d2e674fb6c5c67c.png](https://img-blog.csdnimg.cn/img_convert/dc6e5d07ed8454133d2e674fb6c5c67c.png)
![471e874ac9d024acd73b9c6abc2ae12e.png](https://img-blog.csdnimg.cn/img_convert/471e874ac9d024acd73b9c6abc2ae12e.png)
![b3783f9e9e75517497c5a177c49b8b52.png](https://img-blog.csdnimg.cn/img_convert/b3783f9e9e75517497c5a177c49b8b52.png)
3 变量显示。
![9f12d7cf4ab85a18d67bfd64ccd3d0dd.png](https://img-blog.csdnimg.cn/img_convert/9f12d7cf4ab85a18d67bfd64ccd3d0dd.png)
![695761846fa4d424e5bd09d67e7c9c70.png](https://img-blog.csdnimg.cn/img_convert/695761846fa4d424e5bd09d67e7c9c70.png)
命令流如下:
FORCE,TOTAL
SHELL,TOP
LAYERP26,0
ESOL,5,8560,8678,S,EQV,
NUMVAR,200
XVAR,1
PLVAR,4,