求助！！汽车悬架仿真分析

又个小例子，不知道对你又没有用，共享给你
/clear
/units,si
ks=2535000 !悬架刚度
cs=191000 !悬架阻尼
kt=4280000 !轮胎刚度
ct=98000 !轮胎阻尼
m1=38500 !车体质量
aa=2446000 !车体点头刚度
m2=4330 !悬架及车轮质量
/prep7
et,1,combin14,,2 !1-D纵向弹簧－阻尼单元（UY自由度)
et,2,mass21,,,3 !带转动惯量的2-D质量
et,3,mass21,,,4 !2-D质量不带转动惯量
et,4,mpc184,1 !刚臂梁单元
r,1,ks,cs
r,2,kt,ct
r,3,m1,aa
r,4,m2
r,5
n,1,0
n,2,8.4,0
ngen,3,10,1,2,1,,0.5
n,100,4.2,1
type,1
real,1
e,11,21
e,12,22
type,1
real,2
e,1,11
e,2,12
type,2
real,3
e,100
type,3
real,4
e,11
e,12
type,4
real,5
e,21,100
e,22,100
d,1,uy,,,,,ux
d,2,uy
d,11,ux
d,12,ux
d,100,ux
finish
/solu
antype,modal
modopt,lanb,4
mxpand,4
solve
finish
/post1
set,list
set,first
pldisp,0
anmode,10,0.05