ansys读取 某点应力和应力集中系数的问题

ansys读取某点应力和应力集中系数的问题

2009-10-13 10:20提问者：我浩然哦|浏览次数：3086次

热应力计算结束后后处理我想查看模型某点的热应力，该怎么操作？

还有应力集中部位有多个但我想查看其中某个位置的应力集中系数可以查看吗？怎么操作？

我来帮他解答

精彩回答

2009-10-13 22:33

GUI操作：在General Postproc——Query Results——Subgrid Solu，选择你想显示的节点。

命令流：

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