zsh1975 发表于 2016-5-18 13:12

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

  用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

页: [1]
查看完整版本: 提取节点最大应力应变方法