求助：临界转速在不同转速下一致是什么原因？

qianhuihua

我计算了一个转子的临界转速，压气机-轴-涡轮三部分，两个轴承，计算文件见附件。问题：这样计算合理吗？为什么给不同转速时临界转速未变？希望指导。q_huihua121@sohu.com

ldhltw

直接把命令流或是图形贴出来

qianhuihua

!考虑陀螺力矩，单位制（mm，kpa，kg，s，1e-3N，rad/s,kgm**2）
fini
/clear
FINISH  
/CLEAR,NOSTART  
/COM,ANSYS RELEASE 11.0    UP20070125       17:03:50    08/25/2010
/PREP7  

K,1,1.2,,,  
K,100,-94,,,
K,101,-15.8,,,  
K,200,-70,,,
K,201,-44,,,
K,2,-120,,,
K,3,-113,,,
K,4,-86,,,  
K,5,-81,,,  
K,6,-80,,,  
K,7,-76.4,,,
K,8,-75,,,  
K,9,-36,,,  
SAVE
!*SET,spin,1
/prep7  

!*   轴材料属性
mp,ex,1,2.06e8  
mp,prxy,1,0.3   
mp,dens,1,7.82e-6   

!*  单元类型

ET,1,BEAM188
!*  
ET,2,MASS21
!*  
ET,3,COMBIN14   

!*   集中质量  转动惯量
R,1,0.03,0.03,0.03,4.87e6,3.62e6,3.62e6,
!*  
!*  
R,2,0.07,0.07,0.07,7.73e6,7.48e6,7.48e6,
!*  轴承刚度
R,3,10000000, , ,   

!*  轴的截面尺寸和分网格
SAVE
SECTYPE,   1, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,4.75,0,0,0,0,0,0,0,0,0  
SECTYPE,   2, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,3,0,0,0,0,0,0,0,0,0
SECTYPE,   3, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,4.8,0,0,0,0,0,0,0,0,0   
SECTYPE,   4, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,9,0,0,0,0,0,0,0,0,0
SECTYPE,   5, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,4.5,0,0,0,0,0,0,0,0,0   
SECTYPE,   6, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,7.5,0,0,0,0,0,0,0,0,0   
SECTYPE,   7, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,4,0,0,0,0,0,0,0,0,0
SECTYPE,   8, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,6.7,0,0,0,0,0,0,0,0,0   
LSTR,       2,       3  
LSTR,       3,     100  
LSTR,     100,       4  
LSTR,       4,       5  
LSTR,       5,       6  
LSTR,       6,       7  
LSTR,       8,       7  
LSTR,       8,     200  
LSTR,     200,     201  
LSTR,     201,       9  
LSTR,       9,     101  
LSTR,     101,       1  
ESIZE,0.5,0,


CMSEL,S,_Y1
LATT,1,1,1, , , ,1  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,2   
FITEM,5,-3  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,1,1,1, , , ,2  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
  
CM,_Y,LINE  
LSEL, , , ,       4
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,1,1,1, , , ,3  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,LINE  
LSEL, , , ,       5
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,1,1,1, , , ,4  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,LINE  
LSEL, , , ,       6
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,1,1,1, , , ,5  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,LINE  
LSEL, , , ,       7
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,1,1,1, , , ,6  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,3,4,ORDE,2   
FITEM,5,8   
FITEM,5,-10
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,1,1,1, , , ,7  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  

FLST,5,2,4,ORDE,2   
FITEM,5,11  
FITEM,5,-12
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,1,1,1, , , ,8  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,2,12,4,ORDE,2  
FITEM,2,1   
FITEM,2,-12
LMESH,P51X  
elist,all,,,0,0

!*                                                 压气机集中质量划分网格

!*  
CM,_Y,KP
KSEL, , , ,     100
CM,_Y1,KP   
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
KATT,       1,       1,   2,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*                                                      涡轮
CM,_Y,KP
KSEL, , , ,     101
CM,_Y1,KP   
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
KATT,       1,       2,   2,       0
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
KMESH,     100  
KMESH,     101  
!*                                                           Y方向
KEYOPT,3,1,0
KEYOPT,3,2,2
KEYOPT,3,3,0
!*  
ET,4,COMBIN14   
!*                                                          Z方向
KEYOPT,4,1,0
KEYOPT,4,2,3
KEYOPT,4,3,0
!*  
N,2001,-70,20,,,,,,
N,2002,-70,,20,,,,,
N,2011,-44,20,,,,,,
N,2012,-44,,20,,,,,
TYPE,   3   
REAL,       3   
e,node(-70,0,0),2001
TYPE,   4   
REAL,       3   
e,node(-70,0,0),2002
TYPE,   3   
REAL,       3   
e,node(-44,0,0),2011
TYPE,   4   
REAL,       3   
e,node(-44,0,0),2012

  
ALLSEL,ALL  
FINISH  
/SOL
FLST,2,247,1,ORDE,6
FITEM,2,1   
FITEM,2,-243
FITEM,2,2001
FITEM,2,-2002   
FITEM,2,2011
FITEM,2,-2012   
!*  
/GO
D,P51X, , , , , ,UX,ROTX, , , ,
FLST,2,4,1,ORDE,4   
FITEM,2,2001
FITEM,2,-2002   
FITEM,2,2011
FITEM,2,-2012   
!*  
/GO
D,P51X, , , , , ,ALL, , , , ,   
ESEL,S,TYPE,,1  
CM,shaft,ELEM   
SAVE


ALLSEL,ALL  
ANTYPE,2
MODOPT,QRDAMP,10
EQSLV,FRONT
MXPAND,10, , ,0
LUMPM,0
PSTRES,0
!*  
MODOPT,QRDAMP,10,0,0,1,OFF  
coriolis,on,,,on
!   
mxpand,10,,,yes
cmomega,shaft,0
solve   
mxpand,10,,,yes
cmomega,shaft,200000000*2*3.14/60  
solve   
fini
/post1  
set,list
set,last,last   
plcamp,on,1,rds,0,shaft
prcamp,on,1,rds,0,shaft
plorb   
FINISH

qianhuihua

已解决，可能是单位制问题出错导致的。