ansys 怎样提取某个节点的应力和应变？
时间:2010-11-14 来源:网络 浏览次数: 次
1. 最简单的办法是使用NSORT，打印出结果，可以通过控制使其输出到文件
2. 使用apdl能复杂一点，下面是以前经常用的一段命令流，参考着修改一下吧 *CREATE,GET_node_inf,mac,
*GET,Nnod,NODE,0,COUNT !获取所选择的节点总数 *DIM,S_Xyz,ARRAY,NNOD,5 !定义1个数组存放数据 *GET,Nd,NODE,0,NUM,MIN !获取最小的节点编号 *DO,I,1,Nnod,1
S_Xyz(I,1)=Nd !将节点列表放数组第1列 S_Xyz(I,2)=NX(Nd) !节点的X坐标放数组第2列 S_Xyz(I,3)=NY(Nd) !节点的Y坐标放数组第3列 S_Xyz(I,4)=NZ(Nd) !节点的Z坐标放数组第4列
!*GET,S_Xyz(I,5),NODE,ND,S,EQV !节点的 von mises值放数组第5列 *GET,S_Xyz(I,5),NODE,ND,U,SUM !节点的总变形值值放数组第5列 Nd=NDNEXT(Nd) !读出下一个节点编号 *ENDDO *END
*CREATE,OUT_node_inf,mac, *CFOPEN,node_info,txt,,
*VWRITE,S_Xyz(1,1),S_Xyz(1,2),S_Xyz(1,3),S_Xyz(1,4),S_Xyz(1,5) (F10.0,3F15.4,E15.5) *CFCLOS *END
GET_node_inf OUT_node_inf
/delete,GET_node_inf,mac /delete,OUT_node_inf,mac
另附
1.先对节点的值进行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
2.如把所有的节点应力应变数值提取然后找最大值
*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
3.首先在窗口上显示应力应变云图，然后利用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