ANSYS转子动力学--转子系统不平衡激励的谐响应分析

/PREP7
L1=0.1
L2=0.2
L3=0.3
L4=0.2
L5=0.1
R=0.01
M1=15
R1=0.1
JZ1=1/2*M1*R1**2
JY1=1/2*JZ1
JX1=1/2*JZ1
M2=20
R2=0.2
JZ2=1/2*M2*R2**2
JY2=1/2*JZ2
JX2=1/2*JZ2
KX1=2.63e+7
KY1=2.63e+7
KX2=1.75e+7
KY2=1.75e+7
F1=40E-2
F2=65E-2
ET,1,MASS21
ET,2,BEAM188
ET,3,COMBI214
!*
R,1,M1,M1,M1,JX1,JY1,JZ1,
R,2,M2,M2,M2,JX2,JY2,JZ2,
R,3,KX1,KY1, , , , ,
R,4,KX2,KY2, , , , ,

MP,EX,1, 2.1E11
MP,PRXY,1,0.3
MP,DENS,1,7800

SECTYPE,1, BEAM, CSOLID, , 0
SECOFFSET, ORIG
SECDATA,R,10,10,0,0,0,0,0,0,0,0,0
K,1,,,,
K,2,,,L1,
K,3,,,L1+L2,
K,4,,,L1+L2+L3,
K,5,,,L1+L2+L3+L4,
K,6,,,L1+L2+L3+L4+L5,
L,1,2
L,2,3
L,3,4
L,4,5
L,5,6
ESIZE,0.02,0,
TYPE, 2
MAT, 1
LMESH,ALL

NODEP1=NODE(0,0,L1+L2)
NODEP2=NODE(0,0,L1+L2+L3)

TYPE, 1
REAL, 1
E,NODEP1
REAL, 2
E,NODEP2

NODEZ1=NODE(0,0,L1)
NODEZ2=NODE(0,0,L1+L2+L3+L4)
N,5000,0,0.07,L1
N,5001,0,0.07,L1+L2+L3+L4
TYPE, 3
REAL, 3
E, 5000,NODEZ1
REAL, 4
E, 5001,NODEZ2
ESEL,ALL
ESEL,S,TYPE,,1,2
CM,ZD,ELEM
ALLSEL,ALL
/SOL
ANTYPE,HARMIC
D,ALL,UZ,,,,,ROTZ
D,5000,ALL
D,5001,ALL
CMOMEGA,ZD,0,0,18000/60*2*3.1415
CORIOLIS,1, , ,1,1
F, NODEP1,FX,F1,0
F, NODEP1,FY,0,F1
F, NODEP2 ,FX,F2,0
F, NODEP2, FY,0,F2
OUTRES,ALL,ALL,
HARFRQ,3.8,25.1,
NSUBST,300,
KBC,1
DMPRAT,0.002,
SOLVE
/POST1
SET,1,20
/VIEW,,1,1,1
PLORB !查看转子系统的轨迹图