本帖最后由 赤血冰霜 于 2010-11-19 21:16 编辑
impulse 发表于 2010-11-19 21:03
回复 1 # 赤血冰霜 的帖子
首先,我不是高手,只是觉得lz至少要贴出代码、命令流之类的吧,不如人家怎么给 ...
好吧,下边是我的命令流,大家可以运行一下,先谢谢了
!ex9.16a 梁上移动质量的瞬态分析-位移荷载+耦合法
finish
/clear
/config,nres,2000
/prep7
lb=10 !梁长度
ne=10 !梁划分的单元数
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.015
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
d,1,ux,,,,,uy
d,nn,uy
/solu
antype,2
modopt,lanb,3,,,,1
mxpand,3,,,yes
lumpm,oFF
outres,all,all
!acel,,gra
solve
|