# Analyzing Vibration Mechanism of Angular Contact Ball Bearing with Compound Faults on Inner and Outer Rings

Authors: Lihai Chen; Ma Fang; Ming Qiu; Yanfang Dong; Xiaoxu Pang; Junxing Li; Chuanmeng Yang
Journal: Shock and Vibration
 (2021)

---

## Abstract

This paper investigates a method to dynamically model compound faults on the inner and outer rings of an angular contact ball bearing as well as their effects on its dynamic behavior. Gupta’s dynamic modeling method is used to consider changes in the deformation and direction of the contact load when the ball passes through the damaged area and to develop a dynamic model of compound faults in the angular contact ball bearing. The step-changing fourth-order Runge–Kutta method is used to solve the dynamic compound fault model. The time-domain signal of vibration responses in the case of a single fault in the inner and outer rings exhibited a certain periodicity, and the frequency of faults in the envelope spectrum was clear. By comparison, the periodicity of compound faults was not clear. The signals of compound faults were decomposed by the dual-tree complex wavelet transform to identify their characteristic frequency. Errors occurred between the characteristic frequency of the theoretical fault and its simulated value. They increased with the rotational speed and decreased with an increase in axial load, whereas the influence of radial load on them was minor. For compound faults on the inner and outer rings of an angular contact ball bearing, this study provides a modeling method that can describe changes in the deformation and direction of the contact load when the ball passes through the damaged area of the inner and outer rings. The work here can provide an important foundation for fault identification in angular contact ball bearings.


---

## Body

## 1. Introduction

Angular contact ball bearings (ACBBs) play an important role in the safety and reliability of mechanical systems. However, localized defects due to fatigue, abrasion, and skidding scratches easily occur on raceways and can significantly affect the performance of such systems. To reduce the possibility of damage to the system owing to localized defects in the bearings, their precise and early recognition is important. Many signal processing techniques have been used to identify the characteristic frequency of the localized defects based on the impulse characteristics of the vibration signals of defective bearings [1], including wavelet and correlation filtering [2, 3], time-frequency ridge enhancement [4], intrinsic time-scale decomposition [5], integrated spectral coherence [6], spectral kurtosis [7], and the sparsity-based algorithm [8].A dynamic model of ACBBs with localized defects can be used to investigate and predict their vibrational characteristics. According to whether the rotational characteristics of the bearing elements are considered, such dynamic models can be divided into two types: the simplified spring-mass model and the complete dynamics model. The simplified spring-mass model uses a nonlinear spring and damping to handle the contact characteristics between the balls and the raceways. The dynamic model has been proposed to investigate the effects of multipoint faults in the races on vibrations in deep groove ball bearings [9]. A dynamic model has been proposed for a ball bearing with combined defects, including surface waviness and localized defects [10, 11]. The statistical features of ball bearings based on a simplified spring-mass model have been examined to assess spalling propagation [12]. A complete dynamic model was first presented by Gupta [13], in which the motions of each element were described by differential equations. Many dynamic fault models have been built based on Gupta’s model. To investigate the dynamics and vibrational responses of ACBBs with ball defects, a dynamic model with five degrees of freedom was proposed [14]. An improved time-dependent displacement excitation model was established to consider coupling errors, including roundness and waviness in the inner and outer races of an ACBB [15]. A high-speed fault dynamic model of ACBBs with defects in the races was built by using a B-spline fitting displacement excitation method to represent fault excitation [16]. The effects of cage pocket wear on the dynamic characteristics of ball bearings have been investigated by using a complete dynamics model [17, 18].Some limitations persist in the literature: (1) Simplified models cannot describe complex dynamic behavior and usually focus on the characteristics of planar motion. They are thus not suitable to deal with ACBBs. (2) In the processing for complete dynamic modeling, the fault excitation is characterized by some impulse function that cannot describe changes in the deformation and direction of the contact load when the ball passes through the damaged area. (3) Methods to analyze the vibration mechanism of ACBBs with compound faults on inner and outer rings have rarely been studied.In this paper, a dynamic model of compound faults in ball bearings based on Gupta’s model is established by considering changes in the deformation and direction of the contact load when the ball passes through the damaged area. Simulations are conducted to analyze the mechanism of vibration of ACBBs with compound faults on the inner and outer rings. The correctness of the proposed model was verified through experiments, and errors between the theoretical and the simulated characteristic fault frequencies are explained. This can provide an important foundation for fault diagnosis in ball bearings.


## 2. Dynamic Model of Ball Bearing

We used Gupta’s method [13] to formulate a dynamic model of angular contact ball bearings with compound faults on the inner and outer rings. We assumed that the centers of mass of elements of the bearing coincide with their geometric centers and that the cage pocket is cylindrical. Hertz contact theory and Coulomb’s friction law are used to calculate the normal and tangential forces of all elements.


### 2.1. Coordinate Systems

Some coordinate systems are defined to describe the relative position and motion of elements of the bearing, as shown in Figure1, as follows:(1)
O-X-Y-Z is the inertial coordinate system.(2)
ob-xb-yb-zb is the body-fixed coordinate system of the ball.(3)
or-xr-yr-zr is the body-fixed coordinate system of the ring.(4)
oc-xc-yc-zc is the body-fixed coordinate system of the cage.(5)
op-xp-yp-zp is the coordinate system of the cage pockets.Figure 1
Coordinate systems.Other local coordinate systems are constructed as needed.The transformation matrix A is obtained by using the Cardin angles[13] as follows:(1)Aη,ξ,λ=cosξcosλcosηsinλ+sinηsinξcosλsinηcosλ−cosηsinξcosλ−cosξsinλcosηcosλ−sinηsinξsinλsinηcosλ+cosηsinξsinλsinξ−sinηcosξcosηcosξ.


### 2.2. Interactions between Ball and Rings

According to the geometric relationship shown in Figure2, the position vector from the ball to the inner ring in the body-fixed coordinate system of the ring can be expressed as follows:(2)rbrr=Airrbi−rri,Air=Aφr1,φr2,φr3,where Air is the transformation matrix from the O-X-Y-Z coordinate system to the or-xr-yr-zr coordinate system, and φr1, φr1, and φr3 are the attitude angles of the inner ring.Figure 2
Positional relationships of the ball and rings.Then, the angular positionθbr of the ball in the inner race coordinate system is obtained as follows:(3)θbr=arctan−rbr2rrbr3r.The vector from the center of curvature of the inner race to the center of the ball in the azimuth coordinate system is(4)rbca=AiaArirbrr−rcrr,Aia=Aθbr,0,0,rcrr=0,−rfsinθbr,rfcosθbr,rf=0.5dr+fD,where dr is the diameter of the bottom of the groove of the raceway, f is the curvature factor of the bottom of the groove of the inner raceway, and D is the diameter of the ball.The contact angle between the inner race and the ball can be calculated as follows:(5)α=arctanrbc1arbc3a.The vector from the center of curvature of the inner race to the center of the ball in the contact coordinate system is(6)rbccon=Aaconrbca,Aacon=A0,α,0.Contact deformation between the inner ring and the ball can be given as follows:(7)δ=rbc3con−f−0.5D.Although there is a film of oil between the ball and the raceway, it does not affect the distribution of the contact load. The contact force can thus be obtained according to Hertz’s theory as follows:(8)Fn=0δ≤0Kδ3/2δ>0,K=πkE′Rε4.5ζ31/2,where K is the load deformation coefficient and can be calculated as proposed by Harmrock [19].Because the size of contact ellipse in theycon direction is generally very small, the change in the sliding velocity along this direction can be ignored, as shown in Figure 2. Then, the region of contact ellipse between the ball and the inner ring can be divided into narrow strips in the xcon direction. The vector from the narrow strips P to the center of the ball in the contact coordinate system is(9)rpbcon=x,y,R′2−x2−R′2−a2+1/4D2−a2,R′=2fD2f+1,where a is the semimajor of the contact ellipse. The vector from the narrow strips P to the center of the inner ring in the inertial coordinate system is(10)rpri=AriAarAconarpbcon+rbri.The sliding velocity vector of the inner race relative to the ball at strips Pvrbcon can be calculated as follows:(11)vrbcon=vprcon−vpbcon,where vpbcon=Aiconx˙b,0,r˙bT+Aaconωba×rpbcon(12)vprcon=Aiconvri+Ariωrr−θ˙b,0,0T×rpri,vpbcon=Aiconx˙b,0,r˙bT+Aaconωba×rpbcon,Aicon=AaconAraAir.The tangential forceTP on the strips of P can be obtained by(13)TP=μpFP,FP=3Fn4a1−xP2a2am,where xp is the distance between the center of P and that of the ellipse, m is the number of strips in the semimajor axis, and μp is the equivalent friction coefficient obtained by(14)μp=μbdqbd+μhd1−qbd.In case of a fully flooded EHL,qbd = 0, and then μp can then be written as follows:(15)μp=μhd=A+Bvrbcone−Cvrbcon+D,where A, B, C, and D are the coefficients of a particular lubricant that can be obtained as in Ref. [20]. The components of the tangential force along the semimajor and semiminor axes can be written as follows:(16)TPx=−TPsinθpTPy=TPcosθp,θp=arcsin−vrb1convrb2con.The moment vectors of the contact strip acting on the ball and the ring in the contact coordinate system are(17)Mbrpcon=rpbcon×Tpx,Tpy,0T,Mrbpcon=rpbcon+Aiconrbri×Tpx,Tpy,0T.The total forces and moments on the ball and ring are(18)Fbra=Acona∑p=1mTpx,∑p=1mTpy,FnTFrbi=−AaiFbra,Mbra=Acona∑p=1mMbrpconMrbr=Aconr∑p=1mMrbpcon.


### 2.3. Interactions between Ball and Cage

As shown in Figure3, rb and rc are vectors from the ball and the cage to the center of the inertial coordinate system, respectively, and rpc is the vector from the center of the cage pocket to that of the cage.Figure 3
Positional relationships of the ball and the cage.Then, the vector from the center of the ball to that of the cage pocket can be obtained as follows:(19)rbpp=AcpAicrbi−rci−rpcc.The position of contact between the ball and the pocket can be obtain as follows:(20)ϕ=arctan−rbp1prbp2p.The clearance between the ball and the pocket is(21)h0=Δbp−rbp1p2+rbp2p2,Δbp=dp−D2,where dp is the diameter of the pocket. The critical thickness of the film of oil in the contact state transition Δr is known; if h0>Δr, there is no contact between the ball and pocket; if h0≤Δr, there is a contact, and its force and moment can be calculated as described in Section 2.2.(22)Fbca=AgaFbcgFcbi=Agi−Fbcg,Mbca=Agargbg×FbcgMcbc=Agcrgpg+rpcc×Agc−Fbcg.


### 2.4. Interactions between Cage and Ring

As shown in Figure4, rjc is the vector from the contact point j to the center of the cage and can be written as(23)rjcc=−Lc2,Rcgsinϕ,RcgcosϕT,where Lc is the width of the cage, Rcg is its guide surface, and ϕ is the angle between rjc and the zc-axis. The vector from the contact point j to the center of the ring can be written as(24)rjrr=rcrr+Acrrjcc.where Rrg is the radius of the lip of the raceway.Figure 4
Positional relationships of the cage and the ring. The minimum clearance between the cage and the ring is(25)h0=Rrg−rjr2r2+rjr3r2,Ifh0>Δr, the force vector on the cage can be calculated by(26)Fcrxe=2πη0vcr1eRcg2LcRcg−RrgFcrye=−πη0ucrε2Lc34Rcg−Rrg21−ε23/2Fcrze=−η0ucrε2Lc3Rcg−Rrg21−ε22,vcre=vce−vre,vce=Aievci+AceωccRcg,vre=Aievri+AreωrrRrg,where ucr is(27)ucr=Rcgωc1c+Rrgωr1r.Ifh0≤Δr, the force vector on the cage can be calculated by Hertz contact as described in Section 2.2.The force and moment on the cage can be calculated as follows:(28)Fcri=AeiFcreMcrc=rgcc×AecFcre.


### 2.5. Dynamic Equations and Solutions

To study the law of fault vibration response, two translational degrees of freedom were added to the outer ring of the bearing to simulate the measurement of the acceleration sensor in an empirical test, as shown in Figure5. The differential equation of motion of the outer ring can be written as follows:(29)moz¨o+czz˙o+kzzo=Foz,moz¨o+cyy˙o+kyyo=Foy,where Foy and Foz are the loads of the outer ring in the Y and the Z directions, respectively, mo is the mass of the outer ring, yo and zo are the displacements of the outer ring in the Y and the Z directions, respectively, cy and cz are the respective damping coefficients in the Y and the Z directions, ky and kz are the respective stiffness coefficients in the Y and the Z directions, and the coefficients of stiffness and damping were obtained by the modal test of bearing the pedestal of the test rig for bearing.Figure 5
Model of the outer ring.After all the forces and moments acting on each element of the ball bearing have been calculated, the dynamic equations of each element can be established in different coordinate systems.Newton’s law is used to express the translational motion of any element in the bearing, and the differential equation of motion of the ball can be written as follows:(30)x¨=Fbro1a+Fbri1a+Fbc1amb,r¨b=Fbro3a+Fbri3a+Fbc3a+mrbθ˙bmb,θ¨b=−Fbro2a−Fbri2a+Fbc2a−2mr˙bθ˙bmrb..The differential equation of motion of the cage can be written as follows:(31)x¨c=Fcr1i+∑k=1nFcb1kimc,y¨b=Fcr2i+∑k=1nFcb2kimc,z¨b=Fcr3i+∑k=1nFcb3kimc..Moreover, the equation of the rotational motion is given by Euler’s kinetic equations. Therefore, the equation of rotational motion for each element of the bearing can be written as follows:(32)Ixω˙x−Iy−Izωyωz=∑Mx,Iyω˙y−Iz−Ixωxωz=∑My,Izω˙z−Ix−Iyωxωy=∑Mz..The dynamic equations of the ball bearing system can be numerically solved using the step-changing fourth-order Runge–Kutta method. A flowchart of the dynamic simulation program is shown in Figure6.Figure 6
Flowchart of the dynamic simulation.


### 2.6. Dynamic Modeling of Compound Faults on Inner and Outer Rings

As shown in Figure7, the width of the damaged area, and the depths of the inner and outer rings are wdi, wdo, and hid, respectively, and hod. θie, and θoe are half the center angles corresponding to the damaged areas of the inner and outer rings. θid and θod are the angles of rotation of the center of the damaged area of the inner and outer rings, respectively, and θib and θi+1b are the angles of revolution of the ith and i+1th balls, respectively.Figure 7
Relationship between ball and fault.To determine whether the ball enters the damaged area, it is necessary to calculate the center angle corresponding to the damage and the center of the ball:(33)θbdi=modθi+1b,2π−modθid,2π,θbdo=modθib,2π−modθod,2π..Ifθbdi˂θie, the ball enters the damaged area of the inner ring and otherwise does not. And if θbdo˂θoe, the ball enters the damaged area of the outer ring and otherwise does not.In case the ball does not enter the damaged area, the contact deformations of the ball and the ring can be calculated according to Section2.2. If it does enter the damaged area, the contact deformations need to be calculated by considering the additional deformation caused by the damage. The total deformations are calculated as(34)δdi=δi−hid,θbdi<θie,δi,else,δdo=δo−hod,θbdo<θoe,δo,else,.The contact loads are still calculated according to the Hertz contact load formula in Section2.2.Moreover, as shown in Figure8, the direction of the contact load changes when the ball enters the damaged area. The contact points in the inner and outer rings are Gi and Go, respectively. According to the normal bearing model, the direction of the contact load refers to the center of curvature of the ring and is described in the contact coordinate system. However, once damage has been introduced, the contact load refers to the center of mass of the ball through Gi or Go.Figure 8
Change in the direction of the contact force.For the convenience of description, the damage contact coordinate systems ofGi-xidc-yidc-zidc and Go-xodc-yodc-zodc are established at contact points Gi and Go, respectively.The contact loads can then be obtained as follows:(35)Fidc=Kδdi,Fodc=Kδdo.To introduce the mathematical expression of damage to the dynamic bearing model, it is necessary to transform the damage contact coordinate system into the normal bearing contact coordinate system. For this, it is necessary to obtain the transformation matrix between the coordinate systems. The radius of the ball and the radius of curvature of the rings are known. The angle Ψi(o) between the zi(o)c-axis and the zi(o)dc-axis can be easily determined according to the geometric relationship. The transformation matrix can be expressed as(36)Aiodcc=Aψio,0,0.The flowchart of the simulation program of the dynamic model of compound faults on the inner and outer rings is shown in Figure9.Figure 9
Flowchart of the dynamic fault simulation.


## 2.1. Coordinate Systems

Some coordinate systems are defined to describe the relative position and motion of elements of the bearing, as shown in Figure1, as follows:(1)
O-X-Y-Z is the inertial coordinate system.(2)
ob-xb-yb-zb is the body-fixed coordinate system of the ball.(3)
or-xr-yr-zr is the body-fixed coordinate system of the ring.(4)
oc-xc-yc-zc is the body-fixed coordinate system of the cage.(5)
op-xp-yp-zp is the coordinate system of the cage pockets.Figure 1
Coordinate systems.Other local coordinate systems are constructed as needed.The transformation matrix A is obtained by using the Cardin angles[13] as follows:(1)Aη,ξ,λ=cosξcosλcosηsinλ+sinηsinξcosλsinηcosλ−cosηsinξcosλ−cosξsinλcosηcosλ−sinηsinξsinλsinηcosλ+cosηsinξsinλsinξ−sinηcosξcosηcosξ.


## 2.2. Interactions between Ball and Rings

According to the geometric relationship shown in Figure2, the position vector from the ball to the inner ring in the body-fixed coordinate system of the ring can be expressed as follows:(2)rbrr=Airrbi−rri,Air=Aφr1,φr2,φr3,where Air is the transformation matrix from the O-X-Y-Z coordinate system to the or-xr-yr-zr coordinate system, and φr1, φr1, and φr3 are the attitude angles of the inner ring.Figure 2
Positional relationships of the ball and rings.Then, the angular positionθbr of the ball in the inner race coordinate system is obtained as follows:(3)θbr=arctan−rbr2rrbr3r.The vector from the center of curvature of the inner race to the center of the ball in the azimuth coordinate system is(4)rbca=AiaArirbrr−rcrr,Aia=Aθbr,0,0,rcrr=0,−rfsinθbr,rfcosθbr,rf=0.5dr+fD,where dr is the diameter of the bottom of the groove of the raceway, f is the curvature factor of the bottom of the groove of the inner raceway, and D is the diameter of the ball.The contact angle between the inner race and the ball can be calculated as follows:(5)α=arctanrbc1arbc3a.The vector from the center of curvature of the inner race to the center of the ball in the contact coordinate system is(6)rbccon=Aaconrbca,Aacon=A0,α,0.Contact deformation between the inner ring and the ball can be given as follows:(7)δ=rbc3con−f−0.5D.Although there is a film of oil between the ball and the raceway, it does not affect the distribution of the contact load. The contact force can thus be obtained according to Hertz’s theory as follows:(8)Fn=0δ≤0Kδ3/2δ>0,K=πkE′Rε4.5ζ31/2,where K is the load deformation coefficient and can be calculated as proposed by Harmrock [19].Because the size of contact ellipse in theycon direction is generally very small, the change in the sliding velocity along this direction can be ignored, as shown in Figure 2. Then, the region of contact ellipse between the ball and the inner ring can be divided into narrow strips in the xcon direction. The vector from the narrow strips P to the center of the ball in the contact coordinate system is(9)rpbcon=x,y,R′2−x2−R′2−a2+1/4D2−a2,R′=2fD2f+1,where a is the semimajor of the contact ellipse. The vector from the narrow strips P to the center of the inner ring in the inertial coordinate system is(10)rpri=AriAarAconarpbcon+rbri.The sliding velocity vector of the inner race relative to the ball at strips Pvrbcon can be calculated as follows:(11)vrbcon=vprcon−vpbcon,where vpbcon=Aiconx˙b,0,r˙bT+Aaconωba×rpbcon(12)vprcon=Aiconvri+Ariωrr−θ˙b,0,0T×rpri,vpbcon=Aiconx˙b,0,r˙bT+Aaconωba×rpbcon,Aicon=AaconAraAir.The tangential forceTP on the strips of P can be obtained by(13)TP=μpFP,FP=3Fn4a1−xP2a2am,where xp is the distance between the center of P and that of the ellipse, m is the number of strips in the semimajor axis, and μp is the equivalent friction coefficient obtained by(14)μp=μbdqbd+μhd1−qbd.In case of a fully flooded EHL,qbd = 0, and then μp can then be written as follows:(15)μp=μhd=A+Bvrbcone−Cvrbcon+D,where A, B, C, and D are the coefficients of a particular lubricant that can be obtained as in Ref. [20]. The components of the tangential force along the semimajor and semiminor axes can be written as follows:(16)TPx=−TPsinθpTPy=TPcosθp,θp=arcsin−vrb1convrb2con.The moment vectors of the contact strip acting on the ball and the ring in the contact coordinate system are(17)Mbrpcon=rpbcon×Tpx,Tpy,0T,Mrbpcon=rpbcon+Aiconrbri×Tpx,Tpy,0T.The total forces and moments on the ball and ring are(18)Fbra=Acona∑p=1mTpx,∑p=1mTpy,FnTFrbi=−AaiFbra,Mbra=Acona∑p=1mMbrpconMrbr=Aconr∑p=1mMrbpcon.


## 2.3. Interactions between Ball and Cage

As shown in Figure3, rb and rc are vectors from the ball and the cage to the center of the inertial coordinate system, respectively, and rpc is the vector from the center of the cage pocket to that of the cage.Figure 3
Positional relationships of the ball and the cage.Then, the vector from the center of the ball to that of the cage pocket can be obtained as follows:(19)rbpp=AcpAicrbi−rci−rpcc.The position of contact between the ball and the pocket can be obtain as follows:(20)ϕ=arctan−rbp1prbp2p.The clearance between the ball and the pocket is(21)h0=Δbp−rbp1p2+rbp2p2,Δbp=dp−D2,where dp is the diameter of the pocket. The critical thickness of the film of oil in the contact state transition Δr is known; if h0>Δr, there is no contact between the ball and pocket; if h0≤Δr, there is a contact, and its force and moment can be calculated as described in Section 2.2.(22)Fbca=AgaFbcgFcbi=Agi−Fbcg,Mbca=Agargbg×FbcgMcbc=Agcrgpg+rpcc×Agc−Fbcg.


## 2.4. Interactions between Cage and Ring

As shown in Figure4, rjc is the vector from the contact point j to the center of the cage and can be written as(23)rjcc=−Lc2,Rcgsinϕ,RcgcosϕT,where Lc is the width of the cage, Rcg is its guide surface, and ϕ is the angle between rjc and the zc-axis. The vector from the contact point j to the center of the ring can be written as(24)rjrr=rcrr+Acrrjcc.where Rrg is the radius of the lip of the raceway.Figure 4
Positional relationships of the cage and the ring. The minimum clearance between the cage and the ring is(25)h0=Rrg−rjr2r2+rjr3r2,Ifh0>Δr, the force vector on the cage can be calculated by(26)Fcrxe=2πη0vcr1eRcg2LcRcg−RrgFcrye=−πη0ucrε2Lc34Rcg−Rrg21−ε23/2Fcrze=−η0ucrε2Lc3Rcg−Rrg21−ε22,vcre=vce−vre,vce=Aievci+AceωccRcg,vre=Aievri+AreωrrRrg,where ucr is(27)ucr=Rcgωc1c+Rrgωr1r.Ifh0≤Δr, the force vector on the cage can be calculated by Hertz contact as described in Section 2.2.The force and moment on the cage can be calculated as follows:(28)Fcri=AeiFcreMcrc=rgcc×AecFcre.


## 2.5. Dynamic Equations and Solutions

To study the law of fault vibration response, two translational degrees of freedom were added to the outer ring of the bearing to simulate the measurement of the acceleration sensor in an empirical test, as shown in Figure5. The differential equation of motion of the outer ring can be written as follows:(29)moz¨o+czz˙o+kzzo=Foz,moz¨o+cyy˙o+kyyo=Foy,where Foy and Foz are the loads of the outer ring in the Y and the Z directions, respectively, mo is the mass of the outer ring, yo and zo are the displacements of the outer ring in the Y and the Z directions, respectively, cy and cz are the respective damping coefficients in the Y and the Z directions, ky and kz are the respective stiffness coefficients in the Y and the Z directions, and the coefficients of stiffness and damping were obtained by the modal test of bearing the pedestal of the test rig for bearing.Figure 5
Model of the outer ring.After all the forces and moments acting on each element of the ball bearing have been calculated, the dynamic equations of each element can be established in different coordinate systems.Newton’s law is used to express the translational motion of any element in the bearing, and the differential equation of motion of the ball can be written as follows:(30)x¨=Fbro1a+Fbri1a+Fbc1amb,r¨b=Fbro3a+Fbri3a+Fbc3a+mrbθ˙bmb,θ¨b=−Fbro2a−Fbri2a+Fbc2a−2mr˙bθ˙bmrb..The differential equation of motion of the cage can be written as follows:(31)x¨c=Fcr1i+∑k=1nFcb1kimc,y¨b=Fcr2i+∑k=1nFcb2kimc,z¨b=Fcr3i+∑k=1nFcb3kimc..Moreover, the equation of the rotational motion is given by Euler’s kinetic equations. Therefore, the equation of rotational motion for each element of the bearing can be written as follows:(32)Ixω˙x−Iy−Izωyωz=∑Mx,Iyω˙y−Iz−Ixωxωz=∑My,Izω˙z−Ix−Iyωxωy=∑Mz..The dynamic equations of the ball bearing system can be numerically solved using the step-changing fourth-order Runge–Kutta method. A flowchart of the dynamic simulation program is shown in Figure6.Figure 6
Flowchart of the dynamic simulation.


## 2.6. Dynamic Modeling of Compound Faults on Inner and Outer Rings

As shown in Figure7, the width of the damaged area, and the depths of the inner and outer rings are wdi, wdo, and hid, respectively, and hod. θie, and θoe are half the center angles corresponding to the damaged areas of the inner and outer rings. θid and θod are the angles of rotation of the center of the damaged area of the inner and outer rings, respectively, and θib and θi+1b are the angles of revolution of the ith and i+1th balls, respectively.Figure 7
Relationship between ball and fault.To determine whether the ball enters the damaged area, it is necessary to calculate the center angle corresponding to the damage and the center of the ball:(33)θbdi=modθi+1b,2π−modθid,2π,θbdo=modθib,2π−modθod,2π..Ifθbdi˂θie, the ball enters the damaged area of the inner ring and otherwise does not. And if θbdo˂θoe, the ball enters the damaged area of the outer ring and otherwise does not.In case the ball does not enter the damaged area, the contact deformations of the ball and the ring can be calculated according to Section2.2. If it does enter the damaged area, the contact deformations need to be calculated by considering the additional deformation caused by the damage. The total deformations are calculated as(34)δdi=δi−hid,θbdi<θie,δi,else,δdo=δo−hod,θbdo<θoe,δo,else,.The contact loads are still calculated according to the Hertz contact load formula in Section2.2.Moreover, as shown in Figure8, the direction of the contact load changes when the ball enters the damaged area. The contact points in the inner and outer rings are Gi and Go, respectively. According to the normal bearing model, the direction of the contact load refers to the center of curvature of the ring and is described in the contact coordinate system. However, once damage has been introduced, the contact load refers to the center of mass of the ball through Gi or Go.Figure 8
Change in the direction of the contact force.For the convenience of description, the damage contact coordinate systems ofGi-xidc-yidc-zidc and Go-xodc-yodc-zodc are established at contact points Gi and Go, respectively.The contact loads can then be obtained as follows:(35)Fidc=Kδdi,Fodc=Kδdo.To introduce the mathematical expression of damage to the dynamic bearing model, it is necessary to transform the damage contact coordinate system into the normal bearing contact coordinate system. For this, it is necessary to obtain the transformation matrix between the coordinate systems. The radius of the ball and the radius of curvature of the rings are known. The angle Ψi(o) between the zi(o)c-axis and the zi(o)dc-axis can be easily determined according to the geometric relationship. The transformation matrix can be expressed as(36)Aiodcc=Aψio,0,0.The flowchart of the simulation program of the dynamic model of compound faults on the inner and outer rings is shown in Figure9.Figure 9
Flowchart of the dynamic fault simulation.


## 3. Results and Discussion

The angular contact ball bearing H7005 C was taken as a research object. H7005 C is often used as the spindle bearing in high-speed machine tools; its limit speed is 33000r/min; the range of speed of the tool spindle is generally within 12000r/min. The parameters of its geometry and the operating conditions are summarized in Table1.Table 1
Parameters of bearing and operating conditions.
ParametersValueDiameter of inner ring/mm25Diameter of bottom of groove of inner ring/mm30.425Density of inner ring material/g/cm37.83Diameter of outer ring/mm47Diameter of bottom of groove of outer ring/mm41.562Density of outer ring material/g/cm37.83Diameter of balls/mm5.556Number of balls16Width of bearing/mm12Contact angle/°15Inner diameter of cage/mm36Outer diameter of cage/mm39.27Diameter of cage pocket/mm5.89Density of cage material/g/cm31.3Dynamic viscosity of grease/mm2/sec5.1Radial load/N300Axial load/N100Initial rotational speed of inner ring/r/min1000Stable rotational speed of inner ring/r/min (linear increase to stable speed within 0.1s)10000Initial step size0.0001sSimulation duration0.5sDiameter of damaged area in inner ring/mm0.2Depth of damaged area in inner ring/mm0.1Diameter of damaged area in outer ring/mm0.2Depth of damaged area in outer ring/mm0.1Figures10 and 11 show the vibration signals in case of a single fault in the inner and outer rings in the time domain, and Figure 12 shows the signals in case of a compound fault. The periodicity of the former set of signals was clear.Figure 10
Vibration signals of single inner fault in the time domain.Figure 11
Vibration signals of single outer fault in the time domain.Figure 12
Vibration signals of compound fault in the inner and outer rings in the time domain.A comparison shows that the components of noise of vibration signals in case of a compound fault increased and thus did not have significant periodicity.The theoretical characteristic frequency of single faults in the inner and outer rings can be calculated according to the following formula:(37)f0=Z2×n60×1−Dwdmcosα,fi=Z2×n60×1+Dwdmcosα,where fi = 1532.1 Hz and fo = 1134.6 Hz.Figures13 and 14 show the spectral analysis of vibration signals in case of a single fault in the inner and outer rings. The frequency components in the figures are clean, and the characteristic fault frequency is prominent. These results were consistent with the theoretical calculations, verifying the correctness of the model.Figure 13
Envelope spectra of single fault in the inner ring.Figure 14
Envelope spectra of single fault in the outer ring.To overcome the influence of noise as shown in Figure12, the compound fault signals were decomposed by dual-tree complex wavelet transform, and a series of decomposed components were obtained, as shown in Figure 15. According to the kurtosis criterion, the maximum component of kurtosis was selected for spectral analysis.Figure 15
Signal decomposition based on dual-tree complex wavelet transform.A spectral analysis of the second decomposed component in Figure15 was then conducted. Figure 16 shows the characteristic frequency of compound faults of the inner and outer rings and the doubling of their frequencies.Figure 16
Envelope spectra of compound fault in the inner and outer rings.Different sizes of damage to the inner and outer rings were considered to verify the accuracy of the proposed method. A three-factor four-level orthogonal experiment was arranged, as shown in Table2. As shown in Figure 17, the nine groups of simulations obtained a clear characteristic frequency of the compound fault by using the method proposed in this paper.Table 2
Three-factor four-level experimental design.
Serial numberDiameter of damaged area in inner ring/mmDepth of damaged area in inner ring/mmDiameter of damaged area in outer ring/mmDepth of damaged area in outer ring/mm10.10.050.10.0520.10.10.20.130.10.150.30.1540.20.050.30.150.20.10.10.1560.20.150.20.0570.30.050.20.1580.30.10.30.0590.30.150.10.1Figure 17
envelope spectra of compound fault in the inner and outer rings of different sizes of damage obtained in the simulation.To verify the correctness of the proposed model and the results of analysis, experiments were conducted on a bearing test rig, as shown in Figure18. The bearing test rig consisted of a mechanical structure, driving system, and loading device. The mechanical structure was a cantilever structure, two sets of bearing were used to support the spindle, and the test bearing was located at one end of it. The drive system used the direct drive of the motor to rotate the spindle. The loading device was loaded by a hydraulic cylinder that could apply axial and radial loads on the bearing. The load was reflected through the pressure of the hydraulic cylinder, as shown in Figure 19.Figure 18
Bearing test rig.Figure 19
Loading diagram of bearing test rig.As shown in Figure20, the tested bearing was an angular contact ball bearing with a compound faults in the inner and outer rings constructed by using EDM technology. And the parameters of the bearing are listed in Table 1.Figure 20
Morphology of compound faults. (a)Inner ring (b)Outer ring.
(a)(b)During data processing, the dual-tree complex wavelet transform was used to decompose the vibration signals, and spectral analysis was carried out. As shown in Figure21, in comparison with Figure 16, the frequency component was more complex, and many unknown frequencies appeared. Although some errors occurred, the characteristic frequency of the test was consistent with the simulation value. The correctness of the proposed was thus verified.Figure 21
Envelope spectra of test acceleration responses of the bearing with a compound fault in the inner and outer raceways.To further verify the accuracy of the proposed method, we first tested the bearing according to the parameters listed in Table2 and then examined it on the bearing test rig. Finally, as shown in Figure 22, the nine groups of tests yielded the characteristic frequency of the compound fault by using the proposed method.Figure 22
Envelope spectra of compound fault in the inner and outer rings of different sizes of damage obtained through experiments.Errors occurred between the theoretical and the simulated values. The influence of the operating parameters on the errors was thus examined. As shown in Figure23, the simulation value was greater than the theoretical value in case of a single fault in the inner ring, and the error increased with the rotating speed.Figure 23
Influence of speed of rotation on the characteristic frequency of single fault in the inner ring (Fr = 300 (N) and Fa = 100 N).Similarly, as shown in Figure24, the error increased with rotating speed in case of a single fault in the outer ring but the simulation value was less than the theoretical value.Figure 24
Influence of speed of rotation on the characteristic frequency of single fault in the outer ring (Fr = 300 (N) Fa = 100 N).As shown in Figures25 and 26, regardless of whether the inner or outer ring was involved, the error decreased with an increase in the axial load. In case of a fault in the inner ring, the simulation value was greater than the theoretical value, while for a fault in the outer ring, the simulation value was less than the theoretical value.Figure 25
Influence of axial load on the characteristic frequency of single fault in the inner ring (Fr = 300 (N) and (n) = 10000 r/min).Figure 26
Influence of axial load on the characteristic frequency of single fault in the outer ring (Fr = 300 (N) and (n) = 10000 r/min).From Figures27 and 28, it is clear that the radial load had little effect on the error.Figure 27
Influence of radial load on the characteristic frequency of single fault in the inner ring (Fa = 100 (N) and (n) = 10000 r/min).Figure 28
Influence of radial load on the characteristic frequency of single fault in the outer ring (Fa = 100 (N) and (n) = 10000 r/min).As for the theoretical value and simulation value of characteristic frequency, the law is consistent with the previous conclusion.


## 4. Conclusions

This paper proposed and verified a dynamic model of an angular contact ball bearing, based on Gupta’s dynamic modeling method, by considering changes in the deformation and the direction of contact load when the ball passed through a damaged area. The vibration response signals in case of single and compound faults in the inner and outer rings were studied by using this model.(1)
For a single fault in the inner and outer rings, the signal of the vibration response of the bearing in the time domain exhibited a certain periodicity, and the frequency of fault in the envelope spectrum was clear.(2)
For compound faults in the inner and outer rings, the noise component of the signal increased in the time domain, and a period was not prominent. Once the vibration signal had been processed by the dual-tree complex wavelet method, the characteristic frequency of the complex fault was identified.(3)
Errors occurred between the theoretical and the simulated values and increased with an increase in the rotational speed and decreased with an increase in the axial load. The influence of radial load was minor, however. In case of a fault in the inner ring, the simulation value was greater than the theoretical value, whereas it was smaller than the theoretical value in case of a fault in the outer ring.



---

*Source: 9951110-2021-12-02.xml*