瞬态分析中，如何提取梁单元的轴力

single2581

在作瞬态分析的时候，想提取188梁轴力，在post26中，不知道有没有办法提取梁单元的轴力呢？
   我采用的是post1中的etable命令，由于是瞬态分析，荷载步比较多，我采用了*do循环，发现etable命令好象对循环无效，提取命令如下：

*dim,time,array,3,1
time(1)=0.02,0.05,0.1     !假设有3步
*do,j,1,3,1
SET,j
!------------------建立单元表
ETABLE,NI,SMISC,1     !单元I点轴力
ETABLE,NJ,SMISC,14    !单元J点轴力－梁单元的轴力

fx_res(j+1,1)=time(j)                !定义数组的第一列为时间
!------------------将轴力赋值给数组
*do,i,1,num(1)
    fx_res(1,2*i)=elm(i,1)             !定义数组的第一行为节点号
  fx_res(1,2*i+1)=elm(i,1)
  *GET,fx_res(j+1,2*i),ELEM,i,ETAB,ni
  *GET,fx_res(j+1,2*i+1),ELEM,i,ETAB,nj
*enddo
发现，每次只算了第一步的结果。
欢迎大家讨论，希望给出参考性的意见，谢谢