ANSYS提取层间位移命令流

ANSYS提取层间位移命令流
!提取层间位移时程曲线
!注意：最大变量号按需输入，必须大于：(层数*4)+(层数*5)+层数
! 末尾两行按需输入
! 本命令用于提取x向，如需提取y向，请参考第37行注释
! cg——层高；cs——层数；xxkj——x向跨距；xxks——x向跨数；
! yxkj——y向跨距；yxks——y向跨数
! a,b,c,d用于循环和参数控制，不可更改
! 本文仅做学习使用，如造成其余不良后果本人概不负责
/POST26
NUMVAR,200 !定义最大变量号为200
cg=3
cs=2
xxkj=6
xxks=1
yxkj=6
yxks=1
*DIM,jdh,ARRAY,4,cs !得到每层4个角柱顶点节点号
*DO,i,1,cs,1
ALLSEL
NSLE,S,ALL
NSEL,R,LOC,X,0,xxks*xxkj
NSEL,U,LOC,X,0.01,xxks*xxkj-0.01
NSEL,U,LOC,Y,0.01,yxks*yxkj-0.01
NSEL,R,LOC,Z,i*cg
*GET,NODE,NODE,,COUNT
*GET,MINND,NODE,,NUM,MIN
*DO,ii,1,NODE,1
jdh(ii,i)=MINND
MINND=NDNEXT(MINND)
*ENDDO
NSEL,U,,,ALL
*ENDDO
ALLSEL,ALL
a=2 !得到每个节点的x向位移
*DO,i,1,cs,1
*DO,ii,2,5,1
NSOL,a,jdh(ii-1,i),U,X !如需提取y向，把x改为y即可
a=a+1
*ENDDO
*ENDDO
c=a !将每一层的4个节点相加再除以4
b=2
*DO,i,1,cs,1
ADD,a+1,b,b+1
ADD,a+2,b+2,b+3
ADD,a+3,a+1,a+2
ADD,a+4,a+3,,,,,,1/4
a=a+4
b=b+4
*ENDDO
d=a !用上一层的位移减去下一层的位移得到层间位移
*DO,i,1,cs-1,1
ADD,a+1,c+i*4,c+i*8,,,,,-1,1
a=a+1
*ENDDO
PLVAR,c+4,d+1,d+2,d+3,d+4,d+5 !显示每一层的层间位移，如
有10层，就一直写到d+9
第 1 页
PRVAR,c+4,d+1,d+2,d+3,d+4,d+5 !c+4不能变第 2 页。