Effects of shrink range and preload on dynamics characteristics of antibacklash geared rotorbearing system with composite mesh stiffness
Zheng Yang^{1} , Jianzhong Shang^{2} , Naihui Yu^{3} , Zirong Luo^{4} , Xiaoming Wang^{5}
^{1, 2, 3, 4, 5}College of Mechatronics and Automation, National University of Defense Technology, Changsha, 410073, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 15, Issue 4, 2013, p. 16421655.
Received 20 July 2013; accepted 5 November 2013; published 31 December 2013
JVE Conferences
Antibacklash gear can improve the static transmission precision of system. Besides, the dynamic characteristics of antibacklash gear system have a significant effect on the performance of overall mechanism, and assembly parameter is an important factor to affect the dynamic characteristic. In order to study the effect of assembly parameter, a singlestage antibacklash transmission box was taken as an example. Then, the influence of assembly shrink range and preload force on bearing contact stiffness was considered, and the calculation formula of bearing radial stiffness was derived. The calculation method of antibacklash gear composite mesh stiffness was introduced. Based on this, the transversetorsion geared rotorbearing dynamics model and differential equations of antibacklash gear system with bearing contact stiffness and composite mesh stiffness were proposed. Numerical calculation and the result analysis of frequency sweep show that increasing the shrink range can increase the bearing radial stiffness and further improve the low order resonant frequency (RF) of the system; the effect of inner and outer ring on resonant frequency is equivalent, but neither is obvious. Increasing preload can enhance the bearing radial stiffness and improve the resonant frequency of the system, and the extent of the effect is greater than that of shrink range on the resonant frequency; the effect of preload force on the resonant frequency is obvious.
Keywords: shrink range, preload, antibacklash, geared rotorbearing, dynamics characteristics, composite mesh stiffness.
1. Introduction
The springloaded (or torsionspringloaded) antibacklash gear used for controlling the return difference was invented in 1952 or earlier [1]. At present it is extensively applied in numerous fields, e.g. precision servos [2, 3], radar antennas [4] and precision machine tools. Many precision mechanisms with antibacklash gear require the dynamic characteristics of high speed, high accuracy and high stability. Under the complex work conditions, thus, the internal excitations such as timevarying mesh stiffness, bearing stiffness have effects on dynamic performance, especially vibration characteristics of whole system, which is worth studying.
In the geared rotorbearing system, stiffness including mesh stiffness of gear pair and axial or radial stiffness of bearing is a kind of greatly important parameter. The stability of power transmission and the dynamic performance of system depend on stiffness. As we all know, the factors such as material properties of parts, structural design and processing have an important effect on mechanism stiffness. Besides, the assembly parameters including the shrink range and preload force also have effect on the stiffness [5]. Especially, to the assembly of existing parts, adjusting the assembly parameters is important technical means to improve the dynamic performance of product.
In the recent two decades, the studies on timevary mesh stiffness of gear pair are deep, and the research methods mainly include formula method, such as WeberBanaschck formula and Ishikawa formula, simplified square wave method and FEM.
In the research field of rolling bearing, Harris [5] reviewed and studied the related theory systematically on the rolling bearing. The effect of shrink range on radial internal clearance was analyzed by him. Generally, in order to get a better effect of pretension, additional preload is applied to rolling bearing. In papers [57], pretension types, calculation of preload force and shrink range were analyzed profoundly, and some approximate calculation formulas were derived. Based on these achievements, the calculation formulas of bearing radial, axial and angular stiffness were derived on the basis of rolling bearing quasistatics analysis theory and raceway control theory by Wang and Xia [8], considering the factors of shrink range and preload. In paper [9], the bearing dynamic stiffness/damping matrix was derived in terms of the bearing motions (displacements/rotations) and the bearing analysis considered the bearing nonlinearities, cage rotation and bearing axial preload.
Kahraman and Singh [1011] launched a series of studies on nonlinear dynamics of geared rotor bearing system. In paper [10], nonlinear frequency response characteristics of system were analyzed using RungeKutta numerical integration method and harmonic balance method (HBM) respectively, and factors of gear side clearance and radial internal clearance were both considered in paper [11]. In paper [12], the effects of axial preload on nonlinear dynamic characteristics of a flexible rotor supported by angular contact ball bearings were studied, and the dynamic model of ball bearings was improved for modeling a 5 DOF rotor bearing system. Li and Ren et al. [13] proposed a multiDOF rotor bearing model considering some nonlinear factors. In this paper, the regularity for change of instable rotation speed of geared rotor bearing as amount of bearing radial clearance and aspect ratio was mainly studied. Chen and Kang et al. [14] made some contributions to the analysis of the interactions between gears, shafts and journal bearings in a geared rotorbearing system. In a word, there were many achievements on dynamic analysis of geared rotor bearing system, however, fewer researches considered the effect of assembly factors on bearing stiffness.
At present, however, most of studies about antibacklash gear focus on structure design and calculations, transmission accuracy and applications. There are fewer studies that mentioned mesh stiffness of antibacklash gear transmission and dynamics of antibacklash geared rotor bearing system. In this field, Naoki and Masayoshi [15] proposed an approximate stiffness model of an antibacklash with the characteristic of three segment flexible joint which was used in a manipulator. A dynamic model with friction of antibacklash gear servomechanism by Kwon [2], and the stiffness model of antibacklash gear were similar as three segment flexible joint characteristic. Allan and Levy [3] proposed a method for estimating the minimum preload torque required to obtain a satisfactory step response from a position control system with springloaded antibacklash gear.
In summary, study on system dynamics of antibacklash gear is rare, especially on the antibacklash geared rotorbearing system dynamics, although system dynamics study of ordinary geared rotorbearing is increasingly improved. In this paper, the transversetorsion dynamics mathematical model of a 3 DOF antibacklash geared rotorbearing transmission system will be proposed, considering the composite mesh stiffness of antibacklash gear system and bearing stiffness. The effect of assembly factors, such as shrink range, radial preloaded force of bearing, and stiffness of preloaded spring, on antibacklash geared rotorbearing system dynamics characteristics will be emphatically studied, and some results will provide a reference or basis for further research of dynamics characteristics, optimization design and assembly performance of antibacklash gear system.
2. Antibacklash gear and the composite mesh stiffness
2.1. The structure of antibacklash gear
As shown in Figure 1, gear g1 is named fixed gear; gear g2 is a free gear named loaded gear. The g1 and g2 connected by two or more springs rotate relatively by an angle of some teeth, and springs are preloaded at this point, then g1 and g2 as an integral gear engage with the pinion. Two sides of a tooth of the pinion contact respectively with a tooth of g1 and a tooth of g2 under the action of spring force. Therefore, there is no backlash in transmission no matter it is in positive rotation or in reversal.
Fig. 1. Structure of antibacklash gear
2.2. The composite mesh stiffness of antibacklash gear
The equivalent mesh stiffness of pinion and gear g1 in an antibacklash gear is established as shown in Figure 2. In this figure ${k}_{1}$, ${k}_{2}$, ${k}_{te}$ are respectively the contact stiffness of pinion and gear g1, the contact stiffness of pinion and gear g2, and the stiffness of spring. ${x}_{s1}$, ${x}_{s2}$ are respectively the composite transmission error of pinion and gear g1, pinion and gear g2; ${x}_{te}$ is the deformation of preloaded spring.
Fig. 2. The equivalent model of antibacklash gear mesh stiffness [16]
According to the Yang and Shang et al.’s paper [16], the composite mesh stiffness ${k}_{s}$ is calculated by:
where ${k}_{1}$ and ${k}_{2}$ are:
where ${k}_{A1},{k}_{B1},{k}_{A2},{k}_{B2}$ are respectively the singlepoint mesh stiffness of 4 mesh points [16], and the calculation formula of mesh stiffness ${k}_{{P}_{m}}$ at point ${P}_{m}$ (i. e. point $A1,B1,A2,B2$) was derived by Kuang and Yang [17] as:
where ${k}_{h}$ is the Hertz contact stiffness per unit width [17]; ${B}_{e}$ is the effective tooth width with the unit of mm; $j$ is the number of tooth, $j=1,2$; ${\stackrel{}{k}}_{i}$ is the bending stiffness per unit width in the line of contact [17]; ${r}_{i{P}_{m}}$ is the distance between mesh point ${P}_{m}$ and the center of gear ${g}_{i}$; ${s}_{i}=0.5\left\mathrm{s}\mathrm{g}\mathrm{n}\right(\Delta {\theta}_{i}\left)\right\left[\mathrm{s}\mathrm{g}\mathrm{n}\right(\Delta {\theta}_{i})+1]$, $i=1,2$; $\mathrm{s}\mathrm{g}\mathrm{n}\left(x\right)$ is the sign function, and $\Delta {\theta}_{i}$ is angle difference between the current rotation angle and the critical angle of singletooth mesh area to doubletooth mesh area.
2.3. Contact stiffness of bearing
In this study axial stiffness and angular stiffness of the bearing are not considered.
The problem of the tight fit of bearing and axel, bearing and bearing seat can be treated as a problem of thick cylinder. When inner ring fits tightly with solid axel, the inner ring is expanding, and the inner raceway bottom diameter is increasing. When outer ring fits tightly with bearing seat, the outer ring is shrinking, and the outer raceway bottom diameter is decreasing. Usually, the force balance equations of bearing are written as:
where ${Q}_{rj}$, ${Q}_{aj}$ are separately the radial force and radial force generated by the normal load of ball, and:
where ${\alpha}_{j}$ is the real contact angle of the $j$th ball and raceway; ${\phi}_{j}$ is the position angle of the $j$th ball; ${Q}_{j}$ is the normal contact load of the $j$th ball, and:
where ${\delta}_{nj}$ is the normal contact deformation of the $j$th ball and raceway; ${K}_{n}$_{}is the stiffness coefficient which depends on the material and shape. But for steel bearing ${K}_{n}$ can be calculated by:
${\delta}^{*}=\frac{2F}{\pi}{\left(\frac{\pi}{2{e}^{2}L}\right)}^{1/3}.$
In Eq. (7), $\sum \rho =\text{1/}{R}_{1x}+\text{1/}{R}_{1y}+\text{1/}{R}_{2x}+\text{1/}{R}_{2y}$ is the curvature function, ${R}_{1x}$, ${R}_{1y}$ and ${R}_{2x}$, ${R}_{2y}$ are separately the principal curvature radius of two contact surfaces; the values of $e$, $F$, $L$ are given by Hamrock and Dowson [18] using curve fitting:
$F=1.5277+0.6023\mathrm{l}\mathrm{n}\left({R}_{y}/{R}_{x}\right),$
$L=1.0003+0.5968\left({R}_{x}/{R}_{y}\right),$
$1/{R}_{x}=1/{R}_{1x}+1/{R}_{2x}c,$
$1/{R}_{y}=1/{R}_{1y}+1/{R}_{2y}.$
The normal contact deformation ${\delta}_{nj}$ and the real contact angle ${\alpha}_{j}$ of the $j$th ball and raceway are calculated by:
${\alpha}_{j}={\mathrm{t}\mathrm{a}\mathrm{n}}^{1}\left(\frac{{A}_{0}\mathrm{s}\mathrm{i}\mathrm{n}\alpha +{\delta}_{a}+{\delta}_{0}}{{A}_{0}\mathrm{c}\mathrm{o}\mathrm{s}\alpha +{\delta}_{r}\mathrm{c}\mathrm{o}\mathrm{s}{\varphi}_{j}}\right),$
where ${A}_{0}=({f}_{i}+{f}_{e}1){D}_{b}$; ${\delta}_{a}$ is the radial deformation; ${\delta}_{r}$ is the axial deformation; $\alpha $ is the initial contact angle of bearing after interference fit and applying axial preload, and:
${\delta}_{0}$ is the axial initial deformation after applying axial preload:
here ${\delta}_{i}$, ${\delta}_{e}$ are separately the increment of inner raceway bottom diameter and the decrement of outer raceway bottom diameter [7]; ${u}_{r}$ is the bearing radial clearance; ${D}_{b}$ is the ball diameter; ${f}_{i}$, ${f}_{e}$ are separately the curvature radius coefficient of inner raceway and outer raceway.
Only considering the bearing radial stiffness and under the radial load ${F}_{r}$, the radial stiffness ${k}_{r}$ can be calculated by:
where $\Delta {F}_{r}$, $\Delta {\delta}_{r}$ are separately the increment of load and deformation.
3. The 3DOF transversetorsion dynamics model
As shown in Figure 3a, this is a transmission case with onestage antibacklash gear transmission system. In order to establish a model with enough precision and make the model not too complex, some assumptions and conditions are given as: (1) assume that g1 and g2 work as a gear $g$, and the pinions $p$ and $g$ are connected by a springdamping connector as shown in Figure 3b; (2) assume that the return difference of gear $g$ is equal to zero, and the bearing clearance is considered when modeling; (3) assume that the stiffnesses of shafts are large enough, i. e. the shafts are rigid, and the axial displacements of shafts are neglected; (4) both the axial and radial load are applied to the contact bearing, and ${F}_{a1}^{\text{'}}\left({F}_{a2}^{\text{'}}\right)$, ${F}_{r1}^{\text{'}}\left({F}_{r2}^{\text{'}}\right)$ are separately the axial and radial loads; (5) when modeling, both the low frequency external excitation caused by input fluctuation and the high frequency internal excitation caused by static transmission error are considered, and the output fluctuation is neglected; (6) no friction is considered.
According to the Newton's second law, the transverse – torsion motion differential equations of the antibacklash gear system are written as:
and:
where $\left(\text{\u02d9}\right)$ stands for the full derivative with respect to time $t$; ${T}_{p}$ is the input torque; ${T}_{pm}$ is the mean of ${T}_{p}$; ${T}_{pd}$ is the input torque fluctuation; ${T}_{g}$ is the output torque; ${x}_{p}$, ${x}_{g}$, ${\theta}_{p}$, ${\theta}_{g}$ are separately the transverse and torsion displacements of gears $p$ and $g$; ${m}_{s}$ is the equivalent mass of system; ${F}_{m}^{\text{'}}$ is the average mesh force; ${F}_{d}$ is the force fluctuation relative to the external input torque excitation; ${I}_{p}$, ${I}_{g}$, ${m}_{p}$, ${m}_{g}$ are separately the rotary inertias and masses of gears $p$ and $g$, and ${I}_{g}={I}_{g2}+{I}_{g3}$, ${m}_{g}={m}_{g2}+{m}_{g3}$; ${r}_{p}$, ${r}_{g}$ are separately the base radius of gears $p$ and $g$; ${k}_{b1}$, ${c}_{b1}$, ${k}_{b2}$, ${c}_{b2}$, ${k}_{s}$, ${c}_{s}$ are separately the stiffnesses and dampings of bearing $b1$, bearing $b2$ and gear transmission system; ${f}_{b1}$, ${f}_{b2}$, ${f}_{s}$ are separately the clearance nonlinear functions of bearing $b1$, bearing $b2$ and gear transmission system.
Fig. 3. Example and its dynamics model: a) gear transmission box, b) transversetorsion dynamic model
a)
b)
The composite transmission error $e$ is defined as:
Then Eq. (13) is transformed to the matrix form as:
where $\mathbf{q}={[{x}_{p},{x}_{g},e]}^{\text{T}}$, $\stackrel{}{\mathbf{f}}={\left[{f}_{b1}\right({x}_{p}),{f}_{b2}({x}_{g}),{f}_{s}(e\left)\right]}^{T}$, ${F}_{s}={m}_{s}\ddot{e}$ and:
${f}_{s}\left(e\right)=\left\{\begin{array}{ll}eb,& e>b,\\ 0,& b\le e\le b,\\ e+b,& e<b,\end{array}\right.$
${f}_{b2}\left({x}_{g}\right)=\left\{\begin{array}{ll}{x}_{g}{u}_{r2},& {x}_{g}>{u}_{r2},\\ 0,& {u}_{r2}\le {x}_{g}\le {u}_{r2},\\ {x}_{g}+{u}_{r2},& {x}_{g}<{u}_{r2}.\end{array}\right.$
4. Numerical calculations and analysis
In order to emphasize the effect of timevarying mesh stiffness, shrink range and preload force, etc. on nonlinear dynamics characteristics of system, and considering the fact that when applying the preload force the bearing clearance decreases, in this paper the nonlinear effect caused by bearing radial clearance is neglected. In addition, assume the bearings are of the same size, though the gears are not the same.
An antibacklash gear the parameters of which are shown in Table 1 is taken as a calculation example. The numerical calculation model is established in Matlab/Simulink with the algorithm of fifthorder RungeKutta and fixedstep of 0.0001 s.
Table 1. Parameters of antibacklash gear system
Parameter

Value


$p$

$g$


Number of teeth ${z}_{i}$

20

40

Modulus $m$, mm

0.6


Pressure angle ${\alpha}_{w},$°

20


Addendum circle coefficient ${h}_{a}$

1


Moment of inertia ${J}_{i}$, kg mm^{2}

0.0515

0.258

Mass ${m}_{i}$, g

33.95

48.27

Effective face width $w$, mm

6.5


Tooth thickness ${s}_{i}$, mm

0.8


Stiffness of spring ${k}_{te}$, N mm^{1}

10000


Preload of torsional spring ${T}_{p}$, N mm

1000


Modification coefficient ${X}_{i}$

0

0

Young's modulus $E$, N mm^{2}

2.1E+05


Poisson's ratio $v$

0.3


Datum viscosity coefficient ${\eta}_{s}$, Pa s

0.04


Lubricant viscosity coefficient ${\alpha}_{l}$, m^{2}/N

1.33E08


Initial rotation speed ${\omega}_{0}$, rpm

1000


Drive torque ${T}_{1}$, N mm

100


Load torque ${T}_{2}$, N mm

100

4.1. The effect of shrink range and preload force on bearing radial stiffness
According to the calculation formulas of bearing radial stiffness and the related parameters of bearing referred to paper [9], the bearing radial stiffness ${k}_{r}$ is calculated and the 3D shaded surface figures of the bearing radial stiffness functions are plotted as shown in Figure 4 and Figure 5.
In Figure 4, ${\delta}_{i}$ and ${\delta}_{e}$ both belong to the range of [2e3, 8e3] mm, and the range is limited by the bearing radial clearance ${u}_{r}$, according to Eq. (10). From the Figure 4, the shrink range indeed has effect on bearing radial stiffness, but the effect is not significant in this example. The bearing radial stiffness ${k}_{r}$ increases from 1.0373e+5 N/mm to 1.0381e+5 N/mm, as the axial and radial deformation increases from the minimum to the maximum. With the increment of shrink range both ${\delta}_{i}$ and ${\delta}_{e}$, the bearing radial stiffness increases almost linearly with a small slope as shown in Figure 4b and Figure 4c.
The range of axial and radial preload force is set as [400, 2000] N. As shown in the Figure 5, the axial preload force applying to bearing also has effect on bearing radial stiffness, but the effect is less than that of radial preload force. In the ${k}_{r}\u2013{F}_{a}$ plane, as the increment of ${F}_{r}$, the functional relation of ${k}_{r}$ and ${F}_{a}$ changes from about the parabolic to the linear, and the slope is very small. In the ${k}_{r}\u2013{F}_{r}$ plane, the function relation of ${k}_{r}$ and ${F}_{r}$ is about parabolic, and the effect of ${F}_{r}$ to ${k}_{r}$ is obviously more significant than that of ${F}_{a}$.
Fig. 4. The effect of shrink range on bearing radial stiffness: a) ISO view, b) ${k}_{r}\u2013{\delta}_{i}$ plane and c) ${k}_{r}\u2013{\delta}_{e}$ plane
a)
b)
c)
Fig. 5. The effect of preload force on bearing radial stiffness: a) ISO view, b) ${k}_{r}\u2013{F}_{a}$ plane and c) ${k}_{r}\u2013{F}_{r}$ plane
a)
b)
c)
4.2. Dynamic characteristics analysis of geared rotor bearing system
4.2.1. Single frequency excitation
The frequency response at the single frequency excitation of 3 kHz of the geared rotor bearing system is obtained by numerical calculation as shown in Figure 7, and at this frequency excitation the curve of timevarying mesh stiffness of antibacklash gear is plotted in Figure 8a. Figures 7a, c, e are the time domain responses of the dynamic displacements $x$, ${x}_{p}$ and ${x}_{g}$. From these figures the amplitude of antibacklash gear dynamic transmission error $x$ is the largest, about 0.017 mm, and the amplitudes of dynamic displacements ${x}_{p}$ and ${x}_{g}$ are respectively about 0.8e4 mm and 2.6e4 mm.
The peak arises at the position of 3 kHz, which is caused by external excitation from the frequency spectrum analysis as shown in Figures 6b, d and f. In addition there are also peaks at the positions of frequency division of the external excitation, such as the position of 1.5 kHz, and under the acting of frequency division excitation, the top 3 orders of resonant frequency are exciting.
4.2.2. Frequency sweep excitation
In order to obtain the dynamic response characteristic of the system excited by the excitation of a certain frequency range, the frequency sweep analysis of antibacklash geared rotor bearing system is completed. Frequency sweep means that the excitation frequency increases in a certain range and in some way such as linear, exponential and so on and it is usually used for analyzing the resonance characteristics of the system. So the linear frequency sweep method in which the excitation frequency increases from 0 Hz to 3 kHz with the step size of 60 Hz is used in this study. The main parameters of the system are listed in Table 1.
The Figure 7 shows the curves of antibacklash gear composite mesh stiffness under the actions of single frequency excitation and frequency sweep excitation, and it is obvious that the frequency of mesh stiffness increases with the increment of excitation frequency.
Fig. 6. Responses of the single frequency excitation
a)
b)
c)
d)
e)
f)
Fig. 7. Mesh stiffness at 3 kHz and during frequency sweep: a) response at 3 kHz and b) response of frequency sweep
a)
b)
The Figure 8 shows the responses of the system under the action of frequency sweep excitation, where the Figures 8a, c and e are the figures of time domain response of the dynamic displacements ${x}_{p}$, ${x}_{g}$ and $x$ and the Figures 8b, d and f are the figures of frequency response.
According to the time domain response, the amplitudes of dynamic displacements in the low frequency zone are generally large, which are small in the high frequency zone, and the amplitude of antibacklash gear transmission error $x$ is the largest. From the Figure 9 the same result is obtained. The Figure 9 shows the RMSs (root mean square) of dynamic displacements in the action of different excitation frequencies, where the legends $\Delta $, $*$, $\u2022$ respectively represent the dynamic displacements $x$, ${x}_{p}$ and ${x}_{g}$. It is easy to see that the $x$ is the largest, followed by ${x}_{p}$, ${x}_{g}$ is the smallest.
The main resonant peak of each dynamic displacement is shown in the figure of frequency domain response, and the resonant frequencies for ${x}_{p}$, ${x}_{g}$ and $x$ are respectively about 810 Hz, 503 Hz and 1637 Hz. From the amplitude of main resonant peak, the amplitude of ${x}_{g}$ is the largest and the amplitude of $x$ is the smallest. The reason is that the RMS of $x$ is the largest and the frequency spectrum energy is more dispersed, so the amplitude of the resonant peak becomes small. Likewise, because the RMS of ${x}_{g}$ is the smallest, the frequency spectrum energy is more concentrated.
Fig. 8. Responses of the frequency sweep
a)
b)
c)
d)
e)
f)
In addition to the main resonant peaks, some other resonant peaks can be found in Figures 8b, d and f. Comparing these peaks with each other, it is not difficult to find that the resonant frequency of corresponding peak is basically equal to each other, and just the amplitude has larger difference. It is maybe because that when some part of the system is in resonance, as an “external excitation” for the other parts, the relatively violent vibration causes the forced vibration of other part at this frequency. Of course, this is also worthy of further study.
Fig. 9. The RMSs of dynamic displacements
4.3. The effect of some assembly parameters on the frequency characteristic of the system
4.3.1. The effect of shrink range
From the Figure 5, the effect of inner and outer ring shrink range on the bearing radial stiffness is almost the same, and both are not significant. So, in order to reduce the amount of calculation, the effect of the outer ring shrink range on the frequency response of the system is only considered in the paper. When calculating, set the shrink range of the inner ring equal to 0.005 mm, the radial and the axial preload force both equal to 1 kN, and other parameters refer to Table 1. When frequency sweep is starting, as the other parameters keep the same, the shrink range of outer ring increases linearly from 0.001 mm to 0.01 mm with the step size of 5e4 mm, and the top 3 orders of resonant frequencies of each dynamic displacement in a certain shrink range of outer ring are obtained. Finally, the curves of the relationship between shrink range and resonant frequency of dynamic displacements are plotted in the Figure 10, where the legends ^, * respectively represent the first and second order resonant frequencies.
Fig. 10. Relationship between shrink range and RF: a) the 2 orders RF of dynamic transmission error $x$, b) the 2 orders RF of dynamic displacement ${x}_{p}$ and c) the first order RF of dynamic displacement ${x}_{g}$
a)
b)
c)
According to the calculation result, the third order resonant frequency of each dynamic displacement does not change significantly, for $x$ the third order resonant frequency is about 1609 Hz, and the third order resonant frequency of ${x}_{p}$ and ${x}_{g}$ both are about 1526 Hz. Thus the third order resonant frequency is not shown in Figure 11. Moreover, the second order resonant frequency of ${x}_{g}$ is also not shown, because it is not obvious.
As shown in Figure 10(a), the dynamic transmission error $x$ increases slowly with the increment of shrink range, and the increment of the first order resonant frequency is about 5 Hz, the increment of the second resonant frequency is about 9.8 Hz, both of them are small, and even worse the shrink range is close to the limit. Likewise, in Figure 10b and Figure 10c the increments of the first and second order resonant frequencies of ${x}_{p}$ and ${x}_{g}$ are small too; the increments of the first order resonant frequencies both are about 5 Hz, and the increment of the second order resonant frequency of ${x}_{p}$ is about 8.6 Hz. In a word, it is useful to improve the resonant frequency of the system by increasing the shrink range of the bearing inner and outer ring, but the effect is small.
4.3.2. The effect of preload force
According to the Figure 5, the effect of axial preload force on the bearing radial stiffness is far less than that of radial preload force, thus the effect of axial preload force on the frequency characteristic is not considered.
When frequency sweep is starting, set the shrink range of the inner and outer rings both equal to 0.005 mm, the axial preload forces both equal to 1.2 kN, and other parameters refer to Table 1. As the other parameters remain the same, the radial preload force increases linearly from 200 N to 2 kN with the step size of 100 N, and the top 3 orders of resonant frequencies of each dynamic displacement in the corresponding preload force are obtained. Finally the curves of the relationship between preload force and resonant frequency of dynamic displacement are plotted in the Figure 11, where the legends ^, *, + respectively represent the first, the second and the third order resonant frequencies.
From the Figure 11(a), as the radial preload force increases from 200 N to 2 kN, the first and second order resonant frequencies of dynamic transmission error $x$ increase gradually, the first order resonant frequency increases from 131.8 Hz to 180.7 Hz, the increment is about 49 Hz, the second order resonant frequency increases from 210 Hz to 293 Hz, the increment is about 83 Hz, and the third order resonant frequency is basically unchanged, about 1620 Hz. Thus by the increasing of the bearing radial preload force, the low order resonant frequency of the system can be improved.
The first and second order resonant frequencies of dynamic displacement ${x}_{p}$ increase with the increment of radial preload force as shown in Figure 11(b), where the first order resonant frequency increases from 127 Hz to 178.2 Hz, the second order resonant frequency increases from 212.4 Hz to 295.4 Hz, and the third order resonant frequency is basically unchanged, about 1570 Hz.
Since the changes of the second and third order resonant frequencies of dynamic displacement ${x}_{g}$ are not obvious, only the curve of the first order resonant frequency is shown in Figure 11(c), and the first order resonant frequency and its variation range of ${x}_{g}$ are basically the same as the first order resonant frequencies of $x$ and ${x}_{p}$.
Fig. 11. Relationship between preload force and RF: a) the 3 orders RF of dynamic transmission error $x$, b) the 3 orders RF of dynamic displacement ${x}_{p}$ and c) the first order RF of dynamic displacement ${x}_{g}$
a)
b)
c)
In summary, increasing the bearing preload force has benefit to improve the bearing contact stiffness, and also has positive effect on the improvement of the system resonance frequency. Although the effect of axial preload force on the improvement of the bearing contact stiffness and the system resonance frequency is less than of the radial preload, as shown in Figure 5, they both are positive.
5. Conclusions
(1) The analysis model of the contact bearing considering the assembly shrink range and the preload force was proposed, and the calculation formula of the bearing radial stiffness was given. The results of calculation and analysis show that the bearing radial stiffness can be improved by increasing the shrink range of the inner and outer ring, and by increasing the preload force of the bearing also can improve the bearing radial stiffness, however the effect of axial preload force on the improvement of the bearing contact stiffness is less than of the radial preload.
(2) A single antibacklash gear transmission box was taken as an example, the calculation method of antibacklash gear composite mesh stiffness was used, and the dynamics model with timevarying composite mesh stiffness and bearing stiffness of the antibacklash geared rotor bearing system was proposed, and the dimensionless transversetorsion dynamics differential equations of the system were derived. The proposed model, including some assembly parameters such as shrink range, preload force and antibacklash spring stiffness, provides a basis for the further study on the dynamic characteristic and the control of assembly dynamic performance.
(3) The effect of bearing assembly shrink range and preload force on the system resonant frequency was analyzed. The analysis result shows that the resonant frequencies of the system can be obtained by frequency sweep, and as an excitation, the dynamics behavior of a part can have effect on the other parts; increasing the shrink range of the bearing inner and outer ring is useful to improve the resonant frequency of the system, but the effect is small; increasing the bearing preload force also has positive effect on the improvement of the resonance frequency, and the effect of axial preload force on the improvement of the resonance frequency is less than of the radial preload.
Acknowledgements
We thank the National Natural Science Foundation of China (Grant No. 51175505) for supporting.
References
 Charles W. Cairnes, et al. AntiBacklash Gearing. US Patent: No. 2663198, ID: 22 Dec 1953. [Search CrossRef]
 Kwon Y. S. Rate loop control based on torque compensation in antibacklash geared servo system. Proceedings of the American Control Conference, Boston, Massachusetts, 2004, p. 33273332. [Search CrossRef]
 Allan P. M., Levy N. M. The determination of minimum preload torque for antibacklash gears in a positional servomechanism. IEEE Transactions on Industrial Electronics and Control Instrumentation, Vol. 27, 1980, p. 12321239. [Search CrossRef]
 Li Qing The double gear being applied in construction of radar. Modern Radar, Vol. 22, Issue 5, 2000, p. 8386. [Search CrossRef]
 Harris T. A. Rolling Bearing Analysis. 3rd Edition, John Wiley and Sons, Inc., New York, 1990. [Search CrossRef]
 Qiu Ming, Jiang XingQi, Du YingHui, et al. Rigidity calculation for high speed angular contact ball bearings. Bearing, Issue 11, 2001, p. 58. [Search CrossRef]
 Okamoto Junzo Ball Bearing Design and Calculation. Machinery Industry Press, Beijing, 2003. [Search CrossRef]
 Wang ShuoGui, Xia YuanMing Effect of the interference fit and axial preload in the stiffness of the highspeed angular contact ball bearing. Journal of University of Science and Technology of China, Vol. 36, Issue 12, 2006, p. 13141320. [Search CrossRef]
 EiSaeidy F. M. A. Timevarying total stiffness matrix of a rigid machine spindleangular contact ball bearings assembly: theory and analytical/experimental verifications. Shock and Vibration, Vol. 18, Issue 5, 2011, p. 641670. [Search CrossRef]
 Kahraman A., Singh R. Nonlinear dynamics of a spur geared pair. Journal of Sound and Vibration, Vol. 142, Issue 1, 1990, p. 4975. [Search CrossRef]
 Kahraman A., Singh R. Nonlinear dynamics of a spur geared rotorbearing system with multiple clearances. Journal of Sound and Vibration, Vol. 144, Issue 3, 1991, p. 469506. [Search CrossRef]
 Bai C. Q., Zhang H. Y., Xu Q. Y. Effects of axial preload of ball bearing on the nonlinear dynamic characteristics of a rotorbearing system. Nonlinear Dynamics, Vol. 53, Issue 3, 2008, p. 173190. [Search CrossRef]
 Li ChaoFeng, Ren ZhaoHui, Liu Jie, et al. Bifurcation and stability of the flexible rotorbearing system by a continuum model. Chinese Journal of Applied Mechanics, Vol. 27, Issue 1, 2010, p. 173178. [Search CrossRef]
 Chen Y. C., Kang C. H., Choi S. T. Vibration analysis of geared rotor system under time varying mesh stiffness effects. Journal of Vibroengineering, Vol. 14, Issue 3, 2012, p. 11411150. [Search CrossRef]
 Naoki I., Masayoshi T. Adaptive control of robot manipulators with antibacklash gears. IEEE International Conference on Robotics and Automation, 1995, p. 306311. [Search CrossRef]
 Yang Zheng, Shang Jianzhong, Luo Zirong, Wang Xiaoming, Yu Naihui Research on synthesis meshing stiffness of torsional springloaded doublegear antibacklash mechanism. Journal of Mechanical Engineering, Vol. 49, Issue 1, 2013, p. 2330. [Search CrossRef]
 Kuang J. H., Yang Y. T. An estimate of mesh stiffness and load sharing ratio of a spur gear pair. Proceedings of ASME 12th International Power Transmission and Gearing Conference, Scottsdale, Arizona, Vol. 43, Issue 1, 1992, p. 110. [Search CrossRef]
 Hamrock B. J., Dowson D. Ball Bearing Lubrication. John Wily and Sons, New York, 1981. [Search CrossRef]