HISTORY里的后处理nodal的X-Displacement 图形

林进军

例子出于尚晓江的多层框架弹性结构的弹性时程分析，我按着书上的命令流并查看了书上的K文件，发现都没错。
结构在后处理的anim中也能振动，就是HISTORY功能图形有些怪。


命令流如下：
/prep7
et,1,beam161
et,2,shell163
r,1
rmore
rmodif,1,1,5/6,0.5,0.5,0.25,0.25
rmodif,1,6,0,0
r,2
rmore
rmodif,2,1,5/6,0.4,0.4,0.4,0.4
rmodif,2,6,0,0
r,3
rmore
rmodif,3,1,5/6,,0.1,0.1,0.1,0.1
mp,dens,1,2500
mp,ex,1,3e10
mp,nuxy,1,0.2
n,1,
n,2,5,
n,3,10
n,4,0,0,5
n,5,5,0,5
n,6,10,0,5
ngen,7,6,1,6,1,0,3,0
nsel,s,loc,x,9.5,10.5
nsel,r,loc,y,9.5,18.5
ndele,all
allsel,all
n,100,0,50,0
n,101,5,50,5
n,102,10,50,5
type,2
mat,1
real,3
*do,i,7,37,6
e,i,i+3,i+4,i+1
*enddo
*do,i,8,20,6
e,i,i+3,i+4,i+1
*enddo
type,1
mat,1
real,1
*do,i,10,40,6
e,i,i+1,101
e,i-3,i-2,100
e,i,i-3,100
e,i+1,i-2,101
*enddo
*do,i,11,23,6
e,i,i+1,101
e,i-3,i-2,100
e,i+1,i-2,102
*enddo
type,1
mat,1
real,2
*do,i,1,31,6
e,i,i+6,2
e,i+1,i+7,100
e,i+4,i+10,2
e,i+3,i+9,1
*enddo
*do,i,3,15,6
e,i,i+6,2
e,i+3,i+9,3
*enddo
cp,1,ux,1,2,3,4,5,6
cp,3,uz,1,2,3,4,5,6
nsel,s,loc,y,-0.1,0.1
d,all,,0,,,,uy,rotx,roty,rotz,,
*dim,time,array,301,1,1,,,
*dim,acc,array,301,1,1,,,
*do,i,1,301,1
time(i,1,1)=(i-1)*0.01
acc(i,1,1)=sin(100*((i-1)*0.01))+sin(110*((i-1)*0.01))+sin(120*((i-1)*0.01))
acc(i,1,1)=acc(i,1,1)*500*sin(10*((i-1)*0.01))*(1-cos(2*((i-1)*0.01)))
*enddo
cm,base,node
edload,add,ax,,base,time,acc
allsel,all
eddamp,all,0.1,
/sol
time,5,
edopt,add,blank,lsdyna
edrs,,0.005
edht,,0.005
eddump,1,
edwrite,lsdyna,buiding,k,


附上K文件：
本人研究生期间用ansys/ls-dyna模拟斜拉桥的连续倒塌， 希望有研究相关的人可以加下我的QQ，我们一起讨论进步。QQ 362287820

zb580504

估计是你后处理哪里弄错了吧，我把你附件里面的k文件算了下，同样画了40号节点的x方向位移时程，图形如下：

林进军

zb580504 发表于 2017-4-24 18:30
估计是你后处理哪里弄错了吧，我把你附件里面的k文件算了下，同样画了40号节点的x方向位移时程，图形如下： ...

我之后把K文件拿到DYNA去算了下，在的d3plot放到后处理去，在次画了40号节点X向的位移时程，图形跟你一样。之前之所以出现那个原因我觉得应该是我直接在ANSYS/LS-DYNA里计算，计算后生成的d3plot放到后处理那，并且经过那种计算在时间步也会出错，本来设终止时间是1秒的，要一两分钟才结束。
谢谢你啦