!ex9.17有预应力时梁上移动荷载的瞬态分析finish/clear/config,nres,4000/prep7lb=20!梁长度ne=100!梁划分的单元数nn=ne+1!节点数p=10000!荷载值v=18*1e3/3600!移动速度18km/hdeltl=lb/nedeltt=deltl/v!移动一个单元所需时间em=2.1e11area=0.2im=0.006deng=7800gra=9.8et,1,beam3mp,ex,1,emmp,nuxy,1,0.2mp,dens,1,dengr,1,area,im,1.0*do,i,1,nnn,i,(i-1)*deltl*enddo*do,i,1,nee,i,i+1*enddod,1,ux,,,,,uyd,nn,uyfinish/soluantype,trans!静力分析sstif,ontimint,offtime,1e-5acel,,graf,node(5,0,0),fy,-20e3f,node(10,0,0),fy,-60e3f,node(15,0,0),fy,-30e3nsubst,2kbc,1lswrite,1!瞬态分析timint,onoutres,all,alldeltim,deltt/10kbc,1autots,ontime,1*delttf,1,fy,-plswrite,2*do,i,2,nntime,i*deltt*get,pi1,node,i-1,f,fy*get,pi,node,i,f,fyf,i-1,fy,pi1+pf,i,fy,pi-plswrite,i+1*enddo!荷载移出后的瞬态分析fdele,nn,alltime,lb/v+1.0lswrite,nn+2lssolve,1,nn+2finish/post26nc=node(lb/2,0,0)nsol,2,nc,u,yplvar,2esol,3,50,51,m,zplvar,3