转自 simwe.com
【原创】自己做的转子临界转速计算的小程序,请高手指点
请指教为什么刚度的变化(即使在轴承的刚度选取范围内)会对临界转速的影响这么大,用matlab计算影响就没这么大.ps:如有错误请高手指出!
谢谢!!
fini
/clear
r1=.020
r2=.030
r3=.225
r4=.030
r5=.025
r6=.020
r7=.0225
r8=.018
pi=acos(-1)
/prep7
mp,ex,,2.1e11
mp,dens,,7850
mp,prxy,,0.3
!定义节点
n,1
n,5,.030
fill,1,5,3
n,13,.070
fill,5,13,7
n,17,.095
fill,13,17,3
n,21,.135
fill,17,21,3
n,31,.240
fill,21,31,9
n,35,.270
fill,31,35,3
n,43,.300
fill,35,43,7
n,47,.330
fill,43,47,3
n,48,.015,-.003
n,49,.255,-.003
!定义单元类型
et,1,4 !梁单元
et,2,14 !弹簧单元
et,3,21,,,4 !集中质量单元
!定义实常数
R,1,pi*R1**2,pi*R1**4/4,pi*R1**4/4,2*R1,2*R1, ,
RMORE, ,pi*R1**4/2, , , , ,
R,2,pi*R2**2,pi*R2**4/4,pi*R2**4/4,2*R2,2*R2, ,
RMORE, ,pi*R2**4/2, , , , ,
R,3,pi*R3**2,pi*R3**4/4,pi*R3**4/4,2*R3,2*R3, ,
RMORE, ,pi*R3**4/2, , , , ,
R,4,pi*R4**2,pi*R4**4/4,pi*R4**4/4,2*R4,2*R4, ,
RMORE, ,pi*R4**4/2, , , , ,
R,5,pi*R5**2,pi*R5**4/4,pi*R5**4/4,2*R5,2*R5, ,
RMORE, ,pi*R5**4/2, , , , ,
R,6,pi*R6**2,pi*R6**4/4,pi*R6**4/4,2*R6,2*R6, ,
RMORE, ,pi*R6**4/2, , , , ,
R,7,pi*R7**2,pi*R7**4/4,pi*R7**4/4,2*R7,2*R7, ,
RMORE, ,pi*R7**4/2, , , , ,
R,8,pi*R8**2,pi*R8**4/4,pi*R8**4/4,2*R8,2*R8, ,
RMORE, ,pi*R8**4/2, , , , ,
r,9,1.5e7,0.1
r,10,25.670
type,1
real,1
e,1,2
e,2,3
e,3,4
e,4,5
type,1
real,2
e,5,6
e,6,7
e,7,8
e,8,9
e,9,10
e,10,11
e,11,12
e,12,13
type,1
real,3
e,13,14
e,14,15
e,15,16
e,16,17
type,1
real,4
e,17,18
e,18,19
e,19,20
e,20,21
type,1
real,5
e,21,22
e,22,23
e,23,24
e,24,25
e,25,26
e,26,27
e,27,28
e,28,29
e,29,30
e,30,31
type,1
real,6
e,31,32
e,32,33
e,33,34
e,34,35
type,1
real,7
e,35,36
e,36,37
e,37,38
e,38,39
e,39,40
e,40,41
e,41,42
e,42,43
type,1
real,8
e,43,44
e,44,45
e,45,46
e,46,47
type,2
real,9
e,3,48
e,33,49
type,3
real,10
e,15
d,3,ux
d,33,ux
d,48,all
d,49,all
save
finish
/solu
antype,2
MODOPT,lanb,40
EQSLV,FRONT
MXPAND,40, , ,1
LUMPM,0
PSTRES,0
MODOPT,lanb,40,10,40000, ,OFF
solve
finish |