前后模型不一致?
能否解释清楚点?谢谢!
!ex9.16a 梁上移动质量的瞬态分析-位移荷载+耦合法
finish
/clear
/config,nres,2000
/prep7
lb=10 !梁长度
ne=50 !梁划分的单元数
nn=ne+1 !节点数
nn1=nn+1
ne1=ne+1
p=10000 !荷载值
v=36*1e3/3600 !移动速度36km/h
deltl=lb/ne
deltt=deltl/v !移动一个单元所需时间
em=2.1e11
area=0.18
im=0.0054
deng=7800
gra=9.8
f1=acos(-1)/2/lb/lb*sqrt(em*im/(area*deng))
et,1,beam3
et,2,mass21,,,4
mp,ex,1,em
mp,nuxy,1,0.2
mp,dens,1,deng
r,1,area,im,1.0
r,2,p/gra
*do,i,1,nn
n,i,(i-1)*deltl
*enddo
*do,i,1,ne
e,i,i+1
*enddo
nsel,none
n,nn1
type,2
real,2
e,nn1
nsel,all
d,1,ux,,,,,uy
d,nn,uy
d,nn1,all
/solu
antype,4
sstif,on
timint,off
time,1e-5
acel,,gra
nsubst,2
kbc,1
solve
timint,on
outres,all,all
deltim,deltt/10
kbc,1
autots,on
ddele,nn1,uy
*do,i,1,nn
cpdele,all
d,nn1,ux,(i-1)*deltl
cp,next,uy,i,nn1
time,i*deltt
solve
*enddo
cpdele,all
d,nn1,ux,lb+1
d,nn1,uy
time,lb/v+5/f1
solve
finish
/post26
nc=node(lb/2,0,0)
nsol,2,nc,u,y
plvar,2 |