马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我用三维实体模型分析一个转子-轴承系统的模态,转子是中空的,内径一致,带一个盘,总算不对,请问我的命令流中哪里错了?谢谢各位
/PREP7
!*
ET,1,MESH200
!*
ET,2,SOLID45
!*
ET,3,COMBIN14
!*
ET,4,COMBIN14
!*
KEYOPT,1,1,6
KEYOPT,1,2,0
!*
KEYOPT,3,1,0
KEYOPT,3,2,2
KEYOPT,3,3,0
!*
KEYOPT,4,1,0
KEYOPT,4,2,3
KEYOPT,4,3,0
!*
R,1,10000000, , , !弹簧系数
!*
*DIM,SPIN,,4 ! SPIN VELOCITY (RPM)
SPIN(1) = 0.
SPIN(2) = 20000.
SPIN(3) = 50000.
SPIN(4) = 90000.
SPIN(5) = 70000.
SPIN(6) = 90000.
SPIN(7) = 110000.
SPIN(8) = 130000.
!*
MP,EX,1,2.06e11 !弹模
MP,DENS,1,7900 !密度
mp,prxy,1,0.275 !泊松比
K,1,0,0.006,,
K,2,0,0.0175,,
K,3,0.0555,0.0175,,
K,4,0.0555,0.0198,,
K,5,0.1557,0.0198,,
K,6,0.1557,0.0509,,
K,7,0.1657,0.0509,,
K,8,0.1657,0.025,,
K,9,0.1989,0.025,,
K,10,0.1989,0.033,,
K,11,0.5039,0.033,,
K,12,0.5039,0.0198,,
K,13,0.5969,0.0198,,
K,14,0.5969,0.0175,,
K,15,0.6369,0.0175,,
K,16,0.6369,0.006,,
K,17,0.0555,0.006,,
K,18,0.1557,0.006,,
K,19,0.1657,0.006,,
K,20,0.1989,0.006,,
K,21,0.5039,0.006,,
K,22,0.5969,0.006,,
!*
LSTR,1,2
LSTR,2,3
LSTR,3,4
LSTR,4,5
LSTR,5,6
LSTR,6,7
LSTR,7,8
LSTR,8,9
LSTR,9,10
LSTR,10,11
LSTR,11,12
LSTR,12,13
LSTR,13,14
LSTR,14,15
LSTR,15,16
LSTR,16,1
!*
AL,ALL !创建面
!*
LSTR,3,17
LSTR,5,18
LSTR,8,19
LSTR,9,20
LSTR,12,21
LSTR,14,22
!*
FLST,3,6,4,ORDE,2
FITEM,3,17
FITEM,3,-22
ASBL,1,P51X !用线分割面
ESIZE,0.006,0,
LCCAT,3,17
LCCAT,5,18
LCCAT,7,19
LCCAT,9,20
LCCAT,11,21
LCCAT,13,22
!*
TYPE,1
MSHAPE,0,2D
MSHKEY,1
AMESH,ALL
!*
LSEL,S,LCCA !选择连接生成的线
LDELE,ALL !删除连接线
ALLSEL,ALL !选择所有图元
!*
K,40,,,,
K,41,100,,,
!*
TYPE,2
VROTAT,ALL, , , , , ,40, 41,360, ,
VGLUE,ALL
VSWEEP,ALL
!*
N,40001,0.090865, -0.0198, 0
N,40002,0.56202, -0.0198, 0
N,40003,0.090865, 0, 0.0198
N,40004,0.56202, 0, 0.0198
TYPE,3 !创建轴承单元
REAL,1
E,3262,40001
E,3442,40002
TYPE,4
REAL,1
E,2002,40003
E,2182,40004
!*
FINI
/SOLU
D,ALL,UX ! NO TRACTION & NO TORSION
D,40001,ALL
D,40002,ALL
D,40003,ALL
D,40004,ALL
RATIO = 4*ATAN(1)/30
ANTYPE,MODAL
CORIOLIS,ON,,,ON ! CORIOLIS ON IN A STATIONARY REFERENCE FRAME
NBF = 10
MODOPT,QRDAMP,NBF,,,ON
/OUT,SCRATCH.txt
*DO,I,1,8
OMEGA,SPIN(I)*RATIO
MXPAND,NBF
SOLVE
*ENDDO
FINI
/POST1
PRCAMP,,1.,RPM ! PRINT CAMPBELL VALUES FOR SLOPE=1, UNIT= RPM
PLCAMP,,1.,RPM |