This paper presents a numerical model for investigating the structural dynamics response of a rigid rotor supported on deep-groove ball bearings. The numerical model was used to investigate the influence of race waviness on the dynamic characteristics of a rotor ball bearing system, which is very important from a design viewpoint. The forth-order Runge-Kutta numerical integration technique was applied to determine the time displacement response, Poincarè map, and frequency spectra. The analysis demonstrated that the model can be used as a tool for predicting the nonlinear dynamic behavior of a rotor ball bearing system under different operating conditions. The results of this study may help further understanding of the nonlinear dynamics of a rotor bearing system.
Ball bearings are very important in numerous rotating machinery systems. The prediction and analysis of the dynamic behavior of rotor systems are important because their rotating components possess amounts of energy that can be transformed into vibrations. These vibrations not only affect the performance of rotor system, but may also cause serious damage to the rotating machinery systems. Hence this field has become more challenging since nonlinear analysis is far more difficult compared to the analysis of the linear phenomena. Nonlinearities in rotor systems can be caused by many reasons [1] such as Hertzian contact force, internal radial clearance, surface waviness [2], and stiffness coefficient. Waviness is defined as the geometric imperfection of inner or outer race in ball bearing, and it is also considered one of the sources of vibration [3].The aim of this study is to develop a numerical model for investigating the dynamic properties of rotor system supported by ball bearings under the effect of bearing running waviness. In the analytical formulation, the waviness of rolling elements is modeled by sinusoidal function. The Hertzian contact theory is applied to calculate the elastic deflection and non-linear contact force. A two-degree of freedom system is considered with the assumption that there is no friction between the rolling elements and raceways of ball bearings. The rolling elements are positioned symmetrically such that their motion is synchronized. The 4^{th} order Runge-Kutta numerical integration technique has been applied to solve the nonlinear differential equations. The results are presented in form of time displacement response, Poincarè map, and frequency spectra. The analysis demonstrates that the surface waviness is an important parameter affecting on the dynamic characteristics of rotor bearing system [4, 5]. The model can also be used as a tool for predicting nonlinear dynamic behavior of rotor system under different operating conditions.
2. The Problem Formulation
Since a real rotor bearing system is very complicated and difficult to model [6-10], the analysis of the dynamic behavior of the system is done based on the following assumptions while developing the mathematical model as follows:
1. Torsional vibration of rotor and gyroscopic effects may be neglected and only transverse vibration of rotor should be considered.
2. The outer race of the ball bearing is fixed to a rigid support. The inner race is fixed rigidly to the shaft. There is no slipping of balls.
3. Constant vertical radial force acts on the bearing.
4. Elastic deformation between race and ball gives a nonlinear force deformation relation, which is obtained by using Hertzian theory.
Zeillinger and K Köttritsch [11] noted that additional damping usually occurs between the bearing outer race and its housing. It is known that ball bearing have very low inherent damping. This damping would be effective only for small vibration amplitude [12]. Krämer [13] has provided as estimation of the bearing damping and [14] also the bearing coefficient of ball bearing is well within the range of 33.75 Ns/m to 337.5 Ns/m. A value of c = 200 Ns/m was chosen.
- 2-1. Waviness model
PPT Slide
Lager Image
Waviness is defined as the geometric imperfection of inner or outer race in a ball bearing, and it is considered one of the sources of vibration. These are global sinusoidal shaped imperfections on the outer race surface of the bearing components as in Fig. 1. Waves can be described in terms of two parameters: wavelength (λ), and its amplitude (A). The amplitude of wave at the contact angle can be expressed as follows: where A_{i} is initial wave amplitude and A_{0} is maximum amplitude of wave in radial direction of outer race; L is the arc length of wave of inner race,
where N_{w} is number of wave lobes and θ_{j} is the position of j^{th} ball; r is inner race radius.Because the inner race is moving at the speed of the shaft and the ball center at the speed of the cage so the contact angle (θ_{j}) is given by
PPT Slide
Lager Image
- 2-2. Nonlinear ball bearing force
As shown in Fig. 2, the tangent velocity of the contact point between the ball and the inner race V_{A}, outer race V_{B}, respectively, can be given by
Since the outer race is assumed to be stationary, V_{A} = 0. Therefore, the tangent velocity of the cage is
PPT Slide
Lager Image
PPT Slide
Lager Image
The inner race is fixed to the shaft, ω_{B} = ω_{rotor}. Then, the angular velocity of the cage is written by
PPT Slide
Lager Image
The varying compliance frequency or the ball passage frequency can be expressed in terms of the cage speed times the number of balls, N.According to the Hertzian contact theory [Hertz, 1896], the local Hertzian contact force F_{j} between the j^{th} ball and the race is given as follows [15]:
PPT Slide
Lager Image
The contact stiffness coefficient K, can be given by the stiffness coefficient between the ball and each race [16], k_{i} and k_{0} in series as follows:
PPT Slide
Lager Image
where k_{i} and k_{0} can be determined by elastic modulus and Poisson’s ratio and curvature sum of the contact points as from Harris [15]. The geometrical properties of the system is in Table 1, K = 7.055 × 10^{9} (N/mm^{1.5}).
Based on the Hertzian contact force between inner and outer race and the ball, the total restoring force is the sum of restoring force from each of the rolling elements. The “+” sign subscript in the equation (13, 14) signify this step change in the restoring force expression and models the clearance non-linearity. If the expression inside the bracket is greater than zero, then the ball at the angular location θ_{j} is within the angular contact zone and it is loaded giving rise to a restoring force. If the expression within the brackets is negative or zero, then the ball is not in the load zone, and restoring force is set to zero. So the contact force between the j^{th} ball and inner race can be expressed as,
PPT Slide
Lager Image
PPT Slide
Lager Image
The mathematical model takes into account the sources of nonlinearities in rotor bearing system. After assembling the inertia, ball bearing force, damping force, race waviness and constant vertical force acting on the inner race, the dynamic equations of the system are established as follows:
PPT Slide
Lager Image
PPT Slide
Lager Image
3. Results and Discussion
The nonlinear governing equations of motion (15) are solved by the 4^{th} order Runge-Kutta numerical integration technique to obtain and radial displacement of the rolling element and the shaft. The initial conditions are set to x_{0} = 10^{−6} (m), y_{0} = 10^{−6} (m), x’_{0} = 0 and y’_{0} = 0; K = 7.05 × 10^{9} (N/m^{1.5}). In order to observe the effect of race waviness, time displacement response, Poincarè map, and frequency domain vibrations of shaft rotor bearing system are obtained to determine the dynamic behavior of system [4].The amplitude of race waviness is one of the key factors affecting on the dynamic characteristics of a rotor system. Fig. 3(a, b, c) show the response at a speed of 6,000 rpm without waviness in inner and outer race. For horizontal displacement response, the natural frequency coincides with the varying compliance frequency (ω_{vc} = 354 Hz). The peak amplitude is 0.0846 μm on the horizontal direction. The major peaks amplitude of vibration appear at 3ω_{vc}/2 = 531 Hz, 2ω_{vc} = 708 Hz. The amplitude of peaks are 0.003 μm and 0.0058 μm at 3ω_{vc}/2 = 531 Hz and 2ω_{vc} = 708 Hz, respectively. For horizontal displacement response without outer race waviness and with amplitude of inner race waviness, A_{i} = 0.1 μm as in Fig. 3(d, e, f). The peak amplitude of vibration is 0.0867 μm at ω_{vc} = 354 Hz. The major peaks are at 3ω_{vc}/2 = 531 Hz, and 2ω_{vc} = 708 Hz. The Poincarè maps showed in Fig. 3(b) and 3(e), give an indication of a quasi-periodic response because of “net” structure.
a), b), c) Displacement response with time, Poincarè map for displacement response, and FFT for horizontal displacement response at 6,000 rpm, respectively; when A_{i} = 0, A_{0} = 0. sd), e), f) Displacement response with time, Poincarè map for displacement response, and FFT for horizontal displacement response at 6,000 rpm, respectively; when A_{i} = 0.1, A_{0} = 0.
Fig. 4(a, b, c) show the response with amplitude of inner race waviness, A_{i} = 0.1 μm and amplitude of race waviness, A_{0} = 0.3 at 6,000 rpm. The natural frequency coincides with the varying compliance frequency (ω_{vc} = 354 Hz). The peak amplitude is 0.1018 μm on horizontal direction. The major peaks amplitude of vibration appear at ω_{vc}/4 = 90 Hz, 5ω_{vc}/2 = 890 Hz. Fig. 4(d, e, f) show the response with amplitude of inner race waviness, A_{i} = 0.1 μm and amplitude of race waviness, A_{0} = 0.5 at 6,000 rpm. The natural frequency coincides with the varying compliance frequency (ω_{vc} = 354 Hz). The peak amplitude is 0.111 μm on horizontal direction. The major peaks amplitude of vibration appear at 3ω_{vc}/2 = 531 Hz, 2ω_{vc} = 708 Hz, and 5ω_{vc}/2 = 890 Hz. With the orbits of a dense structure in the Poincarè maps, which give indication of chaotic response, as shown in Fig. 4(b) and 4(e).
a), b), c) Displacement response with time, Poincarè map for displacement response, and FFT for horizontal displacement response at 6,000 rpm, respectively; when A_{i} = 0.1, A_{0} = 0.3. d), e), f) Displacement response with time, Poincarè map for displacement response, and FFT for horizontal displacement response at 6,000 rpm, respectively; when A_{i} = 0.1, A_{0} = 0.5.
It is seen from the solutions of the effect of race waviness by certain conditions, the amplitude response are modulation. Without waviness, the peak amplitudes of vibrations at the varying compliance frequency are insignificant. When increasing amplitude of outer race waviness, the peak amplitudes of vibrations at the varying compliance frequency are more significant. This is also demonstrated by Poincarè maps.
4. Conclusions
In this paper, a two-degree of freedom nonlinear model of a rotor-ball bearing system has been developed to obtain the nonlinear vibration response by varying amplitude of race waviness. Nonlinear analysis of this model was performed numerically using 4^{th} order Runge-Kutta integration method.The results show that the dynamic characteristics behavior of the system is sensitive to small variation of the system parameters. The race waviness of ball bearing is an important parameter for vibration analysis of rotor bearing system and should be considered at the design stages of rotating machinery systems. The model can be used as a tool for predicting nonlinear dynamic behavior of rotor ball bearing system under different operating conditions. Moreover, the study may contribute to a further understanding of the nonlinear dynamics of rotor bearing system.
Nomenclature
m: Mass of rotor, kgc: Equivalent viscous damping factor, Ns/mF_{x}/F_{y}: Hertzian contact force on horizontal/vertical direction, NL: Arc length of wave of surface waviness, μmF_{u}: Force due to unbalance rotor, NK: Hertzian contact stiffness, N/m^{3/2}N/N_{w}: Number of balls / wave lobesR_{B}: Inner race radius, mmR_{A}: Outer race radius, mmr: Inner race radius, mmt: Time, sA_{i}: Initial wave amplitude, μmA_{0}: Maximum wave amplitude, μmV_{cage}: Translational velocity of cage center, mm/sV_{a}: Translational velocity of the inner race, mm/sV_{b}: Translational velocity of the outer race, mm/sW: Radial load, Nγ: Internal radial clearance, μmλ: Wavelength of surface waviness, μmω_{rotor}: Angular speed of rotor, rad/sω_{cage}: Angular speed of the cage, rad/sω_{A}: Angular speed of the inner race, rad/sω_{B}/ω_{b}: Angular speed of the outer race / ball, rad/sθ_{j}: Angular location of jth rolling element, rad/s
Submited paper is presented in 2013 fall conference of KSTLE (2013. 10. 16.~18/yeosu)
Acknowledgements
This research was supported by Yeungnam University Research Fund Program.
Cheng Mei
,
Meng Guang
,
Wu Bingyu
2011
“Nonlinear dynamics of a rotor-ball bearing system with Alford force,”
Journal of Vibration and Control
18
(1)
17 -
27
Hwang Pyung
,
Pan Galina
,
Khan Polina
2004
“Natural Frequency Analysis of Sliders and Head/Disk Interaction Detection by Acoustic Emission,”
J. Korean Soc. Tribol. Lubr. Eng. (International)
5
(1)
28 -
31
Hwang Pyung
,
Nguyen Van Trang
2013
“Dynamic Analysis of Effect of Number of Balls on Rotor-Bearing System,”
J. Korean Soc. Tribol. Lubr. Eng.
29
(4)
248 -
254
DOI : 10.9725/kstle-2013.29.4.248
Fukata S.
,
Gad E. H.
,
Kondou T. A.
,
Tamura H.
1985
“On the radial vibrations of ball bearings (computer simulation),”
Bulletin of the. JSME
28
899 -
904
DOI : 10.1299/jsme1958.28.899
El-Saeidy F. M. A.
1998
“Finite element modeling of rotor-shaft-rolling bearing system with consideration of bearing nonlinearities,”
Journal of Vibration and Control
4
514 -
602
Tiwari M.
,
Gupta K.
,
Prakash O.
2002
“Experimental study of a rotor supported by deep groove ball bearing,”
International Journal of Rotating Machinery
8
(4)
243 -
258
DOI : 10.1155/S1023621X02000234
Tiwari M.
,
Gupta K.
,
Prakash O.
2000
“Effect of radial internal clearance of a ball bearing on the dynamics of a balanced horizontal rotor,”
Journal of Sound and Vibration
238
(5)
723 -
756
DOI : 10.1006/jsvi.1999.3109
Waardle F. P.
1988
“Vibration forces produced by waviness of the rolling surfaces of thrust loaded ball bearings. Part I: Theory,”
Proceeding of the IMechE
202
(C5)
305 -
312
Gupta T. C.
,
Gupta K.
,
Sehgal D. K.
2011
“Nonlinear Dynamics and Chaos of an Unbalanced Flexible Rotor Supported by Deep Groove Ball Bearings with Radial Internal Clearance,”
IUTAM Symposium on Emerging Trends in Rotor Dynamics
Jang G. H.
,
Jeong S. W.
2002
“Nonlinear Excitation Model of Ball Bearing Waviness in a Rigid Rotor Supported by Two or More Ball Bearings Considering Five Degrees of Freedom,”
Transactions of the ASME
124
Citing 'Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing
'
@article{ OHHHB9_2014_v30n1_64}
,title={Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing}
,volume={1}
, url={http://dx.doi.org/10.9725/kstle.2014.30.1.64}, DOI={10.9725/kstle.2014.30.1.64}
, number= {1}
, journal={Journal of the Korean Society of Tribologists and Lubrication Engineers}
, publisher={The Korean Society of Tribologists and Lubrication Engineers}
, author={Hwang, Pyung
and
Nguyen, Van Trang}
, year={2014}
, month={Jan}
TY - JOUR
T2 - Journal of the Korean Society of Tribologists and Lubrication Engineers
AU - Hwang, Pyung
AU - Nguyen, Van Trang
SN - 1229-4845
TI - Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing
VL - 30
PB - The Korean Society of Tribologists and Lubrication Engineers
DO - 10.9725/kstle.2014.30.1.64
PY - 2014
UR - http://dx.doi.org/10.9725/kstle.2014.30.1.64
ER -
Hwang, P.
,
&
Nguyen, V. T.
( 2014).
Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing.
Journal of the Korean Society of Tribologists and Lubrication Engineers,
30
(1)
The Korean Society of Tribologists and Lubrication Engineers.
doi:10.9725/kstle.2014.30.1.64
Hwang, P
,
&
Nguyen, VT
2014,
Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing,
Journal of the Korean Society of Tribologists and Lubrication Engineers,
vol. 1,
no. 1,
Retrieved from http://dx.doi.org/10.9725/kstle.2014.30.1.64
[1]
P Hwang
,
and
VT Nguyen
,
“Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing”,
Journal of the Korean Society of Tribologists and Lubrication Engineers,
vol. 1,
no. 1,
Jan
2014.
Hwang, Pyung
and
,
Nguyen, Van Trang
and
,
“Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing”
Journal of the Korean Society of Tribologists and Lubrication Engineers,
1.
1
2014:
Hwang, P
,
Nguyen, VT
Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing.
Journal of the Korean Society of Tribologists and Lubrication Engineers
[Internet].
2014.
Jan ;
1
(1)
Available from http://dx.doi.org/10.9725/kstle.2014.30.1.64
Hwang, Pyung
,
and
Nguyen, Van Trang
,
“Dynamic Model to Predict Effect of Race Waviness on Vibrations Associated with Deep-Groove Ball Bearing.”
Journal of the Korean Society of Tribologists and Lubrication Engineers
1
no.1
()
Jan,
2014):
http://dx.doi.org/10.9725/kstle.2014.30.1.64