CC BY
43
УДК 534.1
NONLINEAR DYNAMICS OF A ROTOR WITH CANTILEVERED DISK RESTING ON ANGULAR CONTACT BALL BEARINGS
S. Filipkovskyi, Assoc. Prof., Ph. D. (Eng.), Kharkov National Automobile and Highway University
Abstract. The mathematical model of nonlinear oscillations of the rotor resting on angular contact ball bearings is developed. The disc is fixed on the console end of the shaft. The deflection of the shaft, and the elastic deformation of the bearings have the same order. Analysis of free oscillations is carried out, using nonlinear normal modes. The modes and backbone curves of rotor nonlinear oscillations are calculated. The system has soft characteristics.
Key words: rotor, ball bearing, oscillations, nonlinear normal modes, backbone curves.
НЕЛИНЕЙНАЯ ДИНАМИКА РОТОРА С КОНСОЛЬНО ЗАКРЕПЛЁННЫМ ДИСКОМ НА РАДИАЛЬНО-УПОРНЫХ ШАРИКОПОДШИПНИКАХ
С.В. Филипковский, доц., к.т.н., Харьковский национальный автомобильно-дорожный университет
Аннотация. Получена математическая модель нелинейных колебаний ротора на радиально-упорных шарикоподшипниках. Диск закреплён на консольном конце вала. Прогибы вала одного порядка с упругими деформациями подшипников. Анализ свободных колебаний выполнен методом нелинейных нормальных форм. Рассчитаны формы и скелетные кривые колебаний ротора.
Ключевые слова: ротор, шарикоподшипник, колебания, нелинейные нормальные формы, скелетные кривые.
НЕЛ1Н1ЙНА ДИНАМ1КА РОТОРА З КОНСОЛЬНО ЗАКР1ПЛЕНИМ ДИСКОМ НА РАД1АЛЬНО-УПОРНИХ ШАРИКОП1ДШИПНИКАХ
С.В. Фшшковський, доц., к.т.н., Харкчвський нащональний автомобшьно-дорожнш ушверситет
Анотаця. Отримано математичну модель нелттних коливань ротора на радiально-упорних шарикоподшипниках. Диск закртлено на консольному ктщ вала. Прогини вала одного порядку i3 пружними деформащями тдшипниюв. Аналiз вшьних коливань виконано методом нелiнiйних нормальних форм. Розраховано форми i скелетт ^rni коливань.
Ключов1 слова: ротор, шарикотдшипник, коливання, нелiнiйнi нормальш форми, скелетн кривi.
Introduction
Analysis of nonlinear dynamics of machines allows predicting destructive oscillations in conditions, which are safe from the point of view of the linear model of the system, and due to more precise definition to reduce their materials consumption and terms of design. Practically, all vehicles contain rotors supported by nonlinear bearings. Nonlinearity of ball-bearings is
caused by clearances between the balls and races, and the nonlinear dependence of deformations on contact forces. Nonlinear analysis of rotors on ball-bearings with clearances are investigated in articles [1,2,3]. Closing of the radial internal clearance in ball-bearings causes shock loads and excessive vibrations. In order to reduce them, axial preload of ball bearing are used. Nonlinear dynamics of such rotors is investigated in papers [4, 5].
In the majority of works, the oscillations of rotors in which one disk is in the middle between the supports are considered, and oscillations are caused by unbalance or bearings defects. In such models, the deformed shaft center line is approximated by the harmonic functions as a rule. In many machines the rotor resting on axial preloaded angular contact ball-bearings has a disk on the console end, and the bearings are mounted on the vibrating basis. Using such a model it is possible to provide the rotors of helicopter turbines which are used as engines of dump trucks, or electric motors with impellers of pumps which are used in other vehicles.
Problem formulation
Influence of vibration of the basis on nonlinear dynamics of a rotor, which has a disk on console part of a shaft, supported by ball-bearings is insufficient investigated. Therefore the problems of creation of mathematical model and development of technique to research oscillations of a rotor on axial preloaded angular contact ballbearings, and also the analysis of dynamics of a rotor when the frequency of its rotation is in frequency band of vibration of the basis are solved.
Equations of rotor oscillations
For the rotor of such a structure it is difficult to apply harmonic functions, therefore we use the finite element method for approximation of a deformed shaft. The design model of a rotor is presented in fig. 1. Finite elements approximate the sites of a shaft of a constant section. Disks and supports are placed in nodes. Numbers of nodal sections are denoted in fig. 1 by digits 1-5. We consider the forces and the moments of forces of inertia of a disk, and also the contact forces arising in bearings as boundary conditions in the corresponding nodes.
Fig. 1. Finite elements of a rotor and nodes on the ends of elements (1-5)
Coordinate axes are directed, as it is shown in Fig. 1. Generalized coordinates which are the components of a vector of nodal displacements of a node i are the following sequence ui 1 = uix,
Ui,2 = 6i,y :
u, 3 = u,
Ui,4 =6i,x : Ui,5 = Ui,z • In-
terpolation polynoms of finite element are the functions of a bending line of a beam at single movements of nodal sections [7]
C 2 C 3
N,1 (C)=1 -3 7T+21r>
Ne,2 (C) = l
i i2 i3
e,3 (
NeA (C) = I
с2 с3 2 V,
l2 l3 '
fc3 с21
113 l2 J
(2)
where I is the element length, Z is the coordinate along an element axis. Deflections of finite element between nodes i, i +1 are determined by polynoms
uy = Ne,1Ui,1 + Ne,2Ui,2 + Ne,3Ui+U + Ne,4Uг+J,2,
Free oscillations of a shaft of constant section are described by the following equations [6]
EI
El-
aU
ÔÇ 4
54u
+ pF
ас4
- + pF-
5t2
а V.
5t1
= 0,
= 0,
(1)
where I and F - are the second moment of area and the area of the shaft, respectively, E and p -are Young's modulus and the mass density of the shaft, respectively.
uy = Ne,1Ui,3 + Ne,2u i,4 + Ne,3Ui+1,3 + Ne,4Ui+1,4. (3)
Values uiz depend only on time, because the shaft is not deformed along a rotation axis.
The equations of oscillations of a shaft are obtained by Galerkin's method at simultaneous approximation of both the equations and boundary conditions [8]
j WeR0 dO + J WeRT dT = 0,
(4)
where R0 is a residual of the solution of the equation; Rr is a residual in boundary conditions; We and We are the weight functions in the area and on border, respectively; e is a number of finite element. As weight functions in this method we take interpolation polynoms We - Ne .
If expressions (1), (2) and (3) substituted into the first integral (4) we receive the following integrals longwise of an element
J[n г
0
J N ]T
( ш 0 U F д U Л
EI-X + pF x
El-
dC
д 4u
4
+ pF
dt2 j д ч ^
dt2
dC, dC.
(5)
The first boundary condition on the end of a shaft with a disk is equality of the bending moment to the moment of forces of inertia of a disk
EI dX++1 oQd uy
dC2
1 dCdt2 0
d 2uy EI-У +11
d 3u
dC
2
dCdt
-1oO,
2 ^
dCdt
d 4
dCdt
= 0,
C=o
= 0, (8)
C=o
where I1 and I0 are the diametrical and polar moments of inertia of a disk, respectively, Q is an angular speed of rotor. If expressions (2), (3), (8) and Ç = 0 substituted into the second integral (4) we receive an additive to a matrix of masses [MI1 ] and a gyroscopic matrix [Gj ] for degrees of freedom of the corresponding node
Carrying out in (5) integration by parts for terms, which containing derivatives on coordinate Z, we receive
d4 ^
f[N ]T I EI 4 JL eJ { dC4
= [Ke ][ui,1 ui,2
J [N1
0
= K I
dC =
d V
EI y
\
dC4
dC =
j
(6)
where \Ke ] is a stiffness matrix of finite element. Carrying out in (5) integration for terms with derivatives on time we receive
d2u ^
J[N î{pF^ dC =
= Me k,
!]T
J[N ]
0
= [М I
pF
d4 ^ dt2
dC =
v j
ui,3 ui,4 ui+1,3 u i+1,4
(7)
],
N ]T ¡ydu-
dC dt2
+ Ne ]T I/4
C=0
dC dt2
C=0
(9)
[MI1 ][u1,1 u 1,2 «1,3 u 1,4]T
[Ne ]T 10^d2uy
C t [G1 ][u1,1 u1,2
+ [N ]T 10adu
C=0
C t
C=0
. (10)
'l 1,4 F .
The second boundary condition on the end of a shaft with a disk is equality of lateral force to force of inertia of a disk
EI
d ux
EI
dC
3u
3 + m0
C
d\ dt2 _
d \ 3- + m0~dtT
= 0,
C=0
= 0,
(11)
C=0
where m0 is the mass of a disk. If expressions (2), (3), (11) and Ç = 0 substituted into the second integral (4), we receive an additive to a matrix of masses [Mm1 ]
where \Me ] is a matrix of mass of finite element. The components of lines and columns of matrixes corresponding to movement uz will be zero because the shaft is not deformed along a rotation axis, except a diagonal component of a matrix of masses which is equal to the mass of finite element.
[Ne ]
it d 2ux
el
dt
+ [Ne ]T m 0^
C=0
dt2
C=0
(12)
\Mm1 ][u1,1 u 1,2 u1,3 u 1,4 f.
If the disk is fixed on a shaft with eccentricity a, then equations (11) must change so
]
T
]
T
Ч
EI x
El-
ас
а 3u
3 + m0
ас
3- + m0
et2 а 2u
y
~дё
- m0 aQ2 cos Qt = 0,
с=о
- m0 aQ2 sinQt = 0 .
с=0
(13)
If expressions (2), (3), (13) and Ç = 0 substituted into the second integral (4) we receive a right-hand vector of the equations of oscillations {HD (Q, t )} which is caused by a disk unbalance
besides the matrix [Mm1 ]
{HD(Q,t)}= m0aQ2[cosQt 0 sinQt 0]T . (14)
For the node which is fixed in the bearing, boundary conditions on axes x, y are
' ri аuл r
El x
v
ас
3
+
Ji
El
а 4 ^
v
ас3
Px (ux : uy , uz )= 0:
Ji+1
3.. Л
El-
а 3u
y
ас
3
3
El-
а 3u
v J i v
-Py (ux , uy, uz )= 0:
ас3
J i+1
(15)
where Px . and Py . are functions of bearing restoring forces which are received in paper [9]. If expressions (2), (3), (15) and Ç = 0, in case of the left node of an element, or Ç = l, in case of the right node of an element, substituted into the second integral (4) we receive a vector function of the bearing restoring forces
[N ]T Px. (Ux , Uy , U, )Ç=l +
+ [N ]t P . (u , u , u ) =
1/ e jî y,i v*' y ' ^ ç=i
= [Px,i (uг,1, Uг,3, Uz ) 0 Py,i (U^ Ui,3, U, ) 0f =
= K n(U )}.
(16)
If the rotor is mounted on the vibrating base, the vector of kinematic excitation of oscillations is added on the right side of equation [10]
H п(ш, t )}=-[M К(ш, t )},
(17)
where [M ] is a matrix of masses, {/4n(ra, t)} is a vector of vibration accelerations of support,
œ is the angular frequency of vibration of support. Damping forces are concentrated in bearings. Therefore the vector of damping forces has the same structure as the vector {Kn(U)}. In this paper we accept model of viscous damping, then coefficients of damping matrix [C] placed on its diagonal in the same lines as similar component of the vector {Kn (U)}. Assembling the matrixes which received on formulas (6), (7), (9), (10), (12) and vectors (14), (16), (17) we receive the equation of oscillations.
M ]M+G]{U}+[C ]U }+[K ]{UZ }+ + {K n(U )} = {H d (Q, t )}+{H п(ш, t )}•
(18)
Analysis of free oscillations of a rotor
The equation of free oscillations without damping has a form
[M]{t/}+ [G]{U}+ [K]{U}+ {Kn(U)}= 0 . (19)
To analyse free oscillations, we use the method of nonlinear normal modes, which allows to transform the analysis of finite degree of freedoms system to the analysis of single degree of freedom oscillator [9, 11].
Multiply (19) by [M ]- and write the vector of generalized velocities {V }= {U} we receive the first order system of equations
V }+G']{V }+[K ']{U }+{K^(u)}= 0, (20)
where [M j-1 [G] = [G'], [M j-1 [K ] = [K'], [M ]-1{K n(U )}={K n(u)}.
We write all phase coordinates as the functions of one pair of phase coordinates which can be chosen arbitrary [11]
P( p, q) 1 Q( p, q)i
(21)
where p is displacement and q = p is velocity. We can express the components of vector functions {p(p, q)} = [p1,..., pN ]T and
Q(p,q)}=[q1,...,qN]T in the form of a Taylor series:
pn (p, q)=з„,1 p+^„,2 q+s„,3 p2 +
+з„,4 pq2 +^„,6 p 3 +
+^„,7 p 2 q+^n,8 pq2 +^„,9 q3, q„ q )=9n+„,1 p+9n+„,2 q+^N+„,3 p2 +
+ 9N+n,4 pq + ^N+n,sq 2 +9N+n,6 p 3 + (22)
+ 9N+n,7 p 2 q + ^N+n,8 pq 2 +9N+n,9 q
where « = 1,...,J-1,J +1,...,N are numbers of degree of freedoms, J is number of the chosen basic generalized coordinate, p = pJ, q = qJ .
For determination of coefficients of power series (22) we take derivatives of components of vector functions (21):
p„(p,qh^nOMp + q ,
op oq
qn (p, q) = 0q4p^ p +0q«(M)q, (23)
dp
dq
We can express a component of any line of a vector V }= U} in a formula (20) as a function of the chosen pair of phase coordinates. Taking into account series (22) we can write this function in the form:
IN
(pj , qj )=-X g„ ,
n J N
m=1 N
(24)
K'n,mpm K'n,m (PPXP|Jm •
\,mr m m=1 m=1
At the same time this function is the left side of the second formula (23). If expression (24) substituted into (23) the following equations are obtained
p„ ^ q)=qn q )= _ dpn (p, q) ^ dPn q X
-q + -
'Cp' q)'
dp dq
qn (p, q )=u n (p, q )=
= q + uj (p, ,)■
dp dq
(25)
Gathering terms of like powers in p and q we receive system of algebraic equations for definition of coefficients S nm [11].
We substitute the calculated coefficients &nm
into that equation of system (20) which phase coordinates are chosen as p and q. Executing transformations we receive one differential equation of the movement on the chosen vibration mode:
p + B1 p + B2 p + B3 p2 + B4 pp + B5 p2 +
+ Вб p3 + B7 p2 p + B8 pp2 + B9 p3 = 0.
(26)
The equation (26) is solved by method of harmonious balance.
Results of numerical researches
Rotor parameters are as follows: L = 0.34 m is a shaft length; l = 0.06 m is a length of the console end; di = 0.025 m is a diameter of the console end of a shaft; d2 = 0.025 m is a diameter of a shaft between supports; E = 2.11011 Pa and p = 0.3; m = 5.0 kg, Ix = 0.1 kg-m2,
Iz = 0.2 kg-m2; f = Q/2n = 50 Hz is a rotation frequency of a rotor. The standard angular contact ball-bearing parameters are as follows: a = 15° is contact angle; R2 = 27.525 mm is the radius of outer race; R1 = 16.000 mm is the radius of inner race; RK = 5.930 mm is the race radius of curvature; dB = 11.510 mm is diameter of a ball; Nb = 7 is number of balls; E = 2,M0n Pa; M = 0,3.
Frequencies of transverse oscillations of the linearized system are 70,54 Hz, 110,63 Hz, 196,22 Hz and 202,11 Hz. Frequency of longitudinal oscillations of the linearized system is 101,70 Hz. Backbone curves of a rotor are shown in fig. 2.
Fig. 2. Backbone curves of a rotor
The system has soft characteristics. Curves 1 and 2 correspond to oscillations in the funda-
mental mode. With a bigger frequency the curved shaft axle rotates in the direction of shaft rotation, and with a smaller frequency in opposite direction. Curves 3 and 4 correspond to similar oscillations in the second mode.
The fundamental mode of an elastic shaft axle at transverse oscillations when frequency is near 110 Hz is shown in fig. 3 (corresponds to curve 2 in fig. 2).
Д
; V
: i \ 0,2 z
: o,i 4. I о'.з _[
Fig. 3. Fundamental mode of shaft oscillations
At oscillations in the fundamental mode the shaft spindles are from the opposite sides from an axis of bearings. The mode of an elastic shaft axle when frequency is near 202 Hz is shown in fig. 4 (corresponds to curve 4 in fig. 2). In this case the shaft spindles are on the same side from an axis of bearings.
Fig. 4. Second mode of shaft oscillations
Conclusions
Oscillations of a rotor supported by the preloaded angular contact ball bearings are investigated. The disk is fixed on the console end of a shaft. Backbone curves and non-linear normal modes by Shaw and Pierre are obtained. The system has
soft characteristics. Resonant oscillations can occur throughout the whole frequency range below the principal resonance frequency.
References
1. Villa C. Stability and vibration analysis of a
complex flexible rotor bearing system /
C. Villa, J.-J. Sinou, F. Thouverez // Communications in Nonlinear Science and Numerical Simulation. - 2008. - Vol. 13(4). -P. 804-821.
2. Yadav H.K. Study of Effect of Unbalanced
Forces for High Speed Rotor / H.K. Yadav, S.H. Upadhyay, S.P. Harsha // Procedia Engineering. - 2013. - Vol. 64. - P. 593602.
3. Babu C.K. Nonlinear vibration analysis of
an elastic rotor supported on angular contact ball bearings considering six degrees of freedom and waviness on balls and races / C.K. Babu, N. Tandon, R. K. Pandey // Journal of Vibration and Acoustics. - 2014.
- Vol. 136(4). - P. 044503-1-5.
4. Panda K.C. Optimum support characteristics
for rotor-shaft system with preloaded rolling element bearings / K.C. Panda, J.K. Dutt // Journal of Sound and Vibration. -2003. -Vol. 260(4). - P. 731-755.
5. Bai C., Zhang H., Xu Q. Effects of axial
preload of ball bearing on the nonlinear dynamic characteristics of a rotor-bearing system // Nonlinear Dynamics. - 2008. -Vol. 53. - P. 173-190.
6. Dimentberg, F.M. Flexural vibrations of
rotating shaft / F.M. Dimentberg - London: Butterworths, 1961. - 243 p.
7. Бидерман В.Л. Теория механических ко-
лебаний / В.Л. Бидерман. - М.: Высшая школа, 1980. - 408 с.
8. Zienkiewicz, O.C., Finite elements and ap-
proximation / O.C. Zienkiewicz, K. Morgan. - N.Y.: John Wiley & Sons, 1983. -350 p.
9. Filipkovskii, S.V. Nonlinear Free Vibrations
of Multi-Disk Rotors on Ball Bearings / S.V. Filipkovskii, K.V. Avramov // Strength of Materials. - 2013. - Vol. 45(3).
- P. 316-323.
10. Timoshenko, S.P. Vibration Problems in
Engineering / S.P. Timoshenko,
D.H. Young, W. Weaver. - N.Y.: John Wiley & Sons, 1974. - 521 p.
11. Shaw S.W. Normal modes of vibration for
non-linear continuous systems / S.W. Shaw, C. Pierre // Journal of Sound
and Vibration. - 1994. - Vol. 169(3). -P.319-347.
12. Avramov K.V. Nonlinear dynamics of one disk asymmetrical rotor supported by two journal bearings / K.V. Avramov, O.V. Bo-rysiuk // Nonlinear Dynamics. - 2012. -Vol. 67. - P. 1201-1219.
References
1. Villa C., Sinou J.-J., Thouverez F. Stability and vibration analysis of a complex flexible rotor bearing system. Communications in Nonlinear Science and Numerical Simulation, 2008, vol. 13(4). pp. 804-821.
2. Yadav H. K., Upadhyay H. Study of effect of unbalanced forces for high speed rotor. Procedia Engineering, 2013, vol. 64, pp.593-602.
3. Babu C. K., Tandon N.R., Pandey K. Nonlinear vibration analysis of an elastic rotor supported on angular contact ball bearings considering six degrees of freedom and waviness on balls and races. Journal of Vibration and Acoustics, 2014, vol. 136(4), pp. 044503-1-5.
4. Panda K.C., Dutt J.K. Optimum support characteristics for rotor-shaft system with preloaded rolling element bearings. Journal of Sound and Vibration, 2003, vol. 260(4). pp. 731-755.
5. Bai C., Zhang H., Xu Q. Effects of axial preload of ball bearing on the nonlinear dynamic characteristics of a rotor-bearing
system. Nonlinear Dynamics, 2008, vol. 53, pp. 173-190.
6. Dimentberg F.M. Flexural vibrations of rotating shaft. London, Butterworths, 1961. 243 p.
7. Biderman V.L. Teorija mehanicheskih kole-banij [Theory of mechanical oscillations], Moscow, Vysshaja shkola Publ., 1980. 408 p.
8. Zienkiewicz O.C., Morgan K. Finite elements and approximation. John Wiley & Sons, 1983. 350 p.
9. Filipkovskii S.V., Avramov K.V. Nonlinear Free Vibrations of Multi-Disk Rotors on Ball Bearings. Strength of Materials, 2013, vol. 45, no. 3, pp. 316-323.
10. Timoshenko S.P., Young D.H., Weaver W. Vibration Problems in Engineering. John Wiley & Sons, 1974. 521 p.
11. Shaw S.W., Pierre C. Normal modes of vibration for non-linear continuous systems. Journal of Sound and Vibration, 1994, vol. 169(3). pp. 319-347.
12. Avramov K.V., Borysiuk O.V. Nonlinear dynamics of one disk asymmetrical rotor supported by two journal bearings. Nonlinear Dynamics, 2012, vol. 67, pp. 1201-1219.
Рецензент: В.Д. Мигаль, профессор, д.т.н., ХНАДУ.
Статья поступила в редакцию 9 июня 2016 г.
