铁轨-浮置板-梁 系统受简谐力作用下的反应

yjyin

命令流

finish
/clear,nostart
/config,nres,10000
/title,rail under sym force
/filnam,symf response
/PREP7

!定义材料属性

et,1,beam3  
et,2,combin14,1,,2
ET,3,MASS21,,,3  !浮置板集中质量单元

mp,ex,1,2E+011 !rail的材料属性
mp,prxy,1,0.3  
mp,dens,1,7800

mp,ex,2,3.3e10 !梁的材料属性
mp,prxy,2,0.3  
mp,dens,2,2500

mp,damp,4,0

r,1,8.3E8,0,7.5e6  !下弹簧
r,2,1.9,0.388  !梁截面  
r,3,4.8E7,0,1.5E6  !上弹簧
R,4,14400,27648   !浮置板等效质量和转动惯性距长4.8，宽2，厚0.3，密度5000           
r,5,7.708E-3,3.2E-5 !rail截面

!生成节点

n,1
n,18,10.2,0,0
fill
n,19,0,1,0
n,36,10.2,1,0
fill
n,37,0,2,0
n,54,10.2,2,0
fill

!生成梁单元

type,1
mat,2
real,2
*do,k,1,17,1
e,k,k+1
*enddo
!生成钢轨单元
type,1
mat,1
real,5
*do,k,37,53,1
e,k,k+1
*enddo
!上弹簧单元
type,2
real,3
mat,4
*do,k,19,36,1
e,k,k+18
*enddo
!下弹簧单元
type,2
real,1
mat,4
e,1,19
e,18,36
e,9,27
e,10,28

!集中质量单元（等效于浮置板）
type,3
real,4
e,23
e,32

!定义边界条件
d,1,ux,0,,54,1
d,1,uy,0,,18,17

*do,i,1,4,1
CE,next,0,i+32,UY,-1,32,ROTZ,0.6*i,32,uy,1
CE,next,0,-1*i+32,UY,-1,32,ROTZ,-1*0.6*i,32,uy,1

CE,next,0,i+23,UY,-1,23,ROTZ,0.6*i,23,uy,1
CE,next,0,-1*i+23,UY,-1,23,ROTZ,-1*0.6*i,23,uy,1
*enddo

finish

!定义求解选项

/solu
antype,4
trnopt,full

!定义荷载步

*do,n,1,50,1
ft=-2.5e4*sin(10*n*0.03)
f,41,fy,ft
f,50,fy,ft
time,3e-2*n
deltim,3e-2
autots,-1
kbc,0
ALPHAD,0.2
BETAD,0.2
outres,all
outpr,basic,all

lswrite,n
*enddo

!求解

lssolve,1,50,1

yjyin

没人顶，老大也不给加分。郁闷。

zhouqingcam

研究一下...

liuzongfahjx

先顶一下

gaodukongzhifa

做动力学的吗

asdzhy2005

研究一下

dxdcll

好贴要顶

zy0512

谢谢楼主！

彼得大帝

谢谢楼主，好贴

zuchongzhi

谢谢楼主啊

newagan

加了轨道的不平顺了吗？