ansys怎样提取某个节点的应力和应变

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)

,,

*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