提取节点最大应力应变方法

zsh1975

　　用ANSYS的APDL方式直接获得最大应力和最大应变的节点编号及其数值的有好多种方法：
　　如下面的这种，先对节点的值进行SORT，在提取最大的值即可。

　　最大应力节点编号及其数值的提取：
　　ALLSEL
　　NSORT,S,EQV,0,0,ALL
　　*GET,MAX_SEQV,SORT,0,IMAX
　　*GET,MAXSEQV,NODE,MAX_EQV,S,EQV

　　依次类推，可以最大应变节点编号及其值数值的提取：
　　ALLSEL
　　NSORT,EPTO,EQV,0,0,ALL
　　*GET,MAX_EPTOEQV,SORT,0,IMAX
　　*GET,MAXEPTOEQV,NODE,MAX_EPTOEQV,S,EQV

　　最大位移节点编号及其数值的提取：
　　ALLSEL
　　NSORT,U,SUM,0,0,ALL
　　*GET,MAX_U,SORT,0,IMAX
　　*GET,MAXU,NODE,MAX_U,U,SUM

　　还有其他方法：如把所有的节点应力应变数值提取然后找最大值
　　*GET,NUMALL,NODE,0,COUNT
　　*GET,NUMSTART,NODE,0,NUM,MIN
　　SMAX=0
　　INUSE=NUMSTART
　　*DO,I,1,NUMALL,1
　　*GET,SSUM,NODE,INUSE,S,EQV
　　*IF,SMAX,GE,SSUM,THEN
　　SMAX=SMAX
　　*ELSEIF,SMAX,LT,SSUM,THEN
　　SMAX=SSUM
　　*ENDIF
　　INUSE=NDNEXT(INUSE)
　　*ENDDO

　　还有一种是首先在窗口上显示应力应变云图，然后利用GET命令得到最大值。
　　如提取最大应力为：
　　PLNSOL,S,EQV,0
　　*GET, PAR, PLNSOL, 0, MAX
　　或将所有的节点定义成一个路径，*GET, PAR, PATH, 0, MAX，
　　另：获取结果各个时间点的最大应力值
　　*DIM,MAXSEQVNODE,TABLE,63,1,1,,, !定义数组
　　*DO,ISET1,1,63,1 !根据结果存储的数目而定,如63
　　/POST1
　　SET,1,,1,,,,ISET1,
　　PLNSOL,S,EQV,0
　　*GET,MAXSEQV_NODEI,PLNSOL,0,MAX !获取结果,这里是最大MISES等效应力
　　*SET,MAXSEQV_NODE(ISET-1,0,1),(ISET-1)*4E-4 !确定对应的时间
　　*SET,MAXSEQV_NODE(ISET-1,1,1),MAXSEQV_NODEI !写入最大等效应力值
　　*ENDDO


转自：http://blog.sina.com.cn/s/blog_4f7a7a1b0100ay1m.html