时程分析后多个BEAM188单元的应力提取后处理问题

bambooli

BEAM188单元的应力提取比较麻烦，结合论坛上一位仁兄的帖子以及王新敏的书，本人写了一段程序来提取时程分析后BEAM188单元的应力。共有16个单元的应力需要提取，一个思路就是每一个荷载步作用下，结构有一个应力响应，提取出最大值和最小值，然后在1000个荷载步里面找出最大最小值。
命令流如下：
/POST1
  set,last !读取最后一个荷载步的解
  *GET,NUM_SUB,ACTIVE,0,SOLU,NCMLS !获取荷载步的总步数，Cumulative number of load steps，提取后值为1000
  *DIM,DMAX,ARRAY,NUM_SUB  ！存储1000个荷载步每步的最大应力
  *DIM,S_MAX,ARRAY,16             ！存储所有1000个荷载步中中某单元的最大应力，共16个单元
  *DIM,DMIN,ARRAY,NUM_SUB  ！存储1000个荷载步每步的最小应力
  *DIM,S_MIN,ARRAY,16             ！存储所有1000个荷载步中某单元的最大应力

  *DO,i,1,16
   S_MAX(i)=0     !数组初始化为0
   S_MIN(i)=0       !数组初始化为0
  *ENDDO
  
*DO,i,1,16         
*DO,j,1,NUM_SUB
  SET,j  !读取第j个荷载步的解
  *GET,DMAX(j),SECR,DunTop(i),S,X,MAX   ! SX最大应力值
                                     ！Duntop里面保存的是需要提取的单元之单元号
  *GET,DMIN(j),SECR,DunTop(i),S,X,MIN     ! SX最小应力值
  *IF,DMAX(j),GE,S_MAX(i),THEN
      S_MAX(i)=DMAX(j)
    *ELSE
      S_MAX(i)=S_MAX(i)
  *ENDIF
  *IF,DMIN(j),LE,S_MIN(i),THEN
      S_MIN(i)=DMIN(j)
    *ELSE
      S_MIN(i)=S_MIN(i)
  *ENDIF
*ENDDO
*ENDDO
运行完毕后，发现只有第一个单元即i=1时对应的单元有值提取出来，但是j只运行了一次。output窗口见附件。
附件：

请大神帮忙看看，若问题解决了，本人定当将修改好的命令流发到论坛与大家共享。
谢谢了！

bambooli

运行完毕后，output窗口正常。感觉是两层DO循环那里出问题了。请大神指教。

zhangxw023

两个循环。。。不简单呀
我最近正在做一个beam161单元的动力学分析，也需要提取所有载荷步中的最大塑性变形。正在学习中，嘿嘿

Timhhu

能否发个命令流学习一下

laohan0916

没看懂。。。

baiyin511

*do,i,1,16,1
是吧?

fengying52

领教高手了