|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
在作瞬态分析的时候,想提取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
发现,每次只算了第一步的结果。
欢迎大家讨论,希望给出参考性的意见,谢谢 |
|