瞬态分析中，如何提取梁单元的轴力
在作瞬态分析的时候，想提取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
发现，每次只算了第一步的结果。
欢迎大家讨论，希望给出参考性的意见，谢谢
页:
[1]