Research Article
The bearing parallel misalignment running model–based machine tool spindle thermal error analysis
Advances in Mechanical Engineering 2018, Vol. 10(6) 1–7 Ó The Author(s) 2018 DOI: 10.1177/1687814018778629 journals.sagepub.com/home/ade
Yanfang Dong1 , Zude Zhou2 and Lihai Chen1
Abstract As a key component of the machine tool spindle, bearing has critical influences on the spindle thermal error. In particular, the installation errors of bearing have considerable effects upon the spindle thermal error by altering the bearings’ internal contact angles, contact loads, and friction torques for different ball positions, but have yet to be fully elucidated. In this article, the influence of installation errors on the resulting spindle thermal error was evaluated using both empirical methods and simulation method, with the ultimate aim of reducing installation error. Deviations within the bearing support were used to simulate bearing parallel misalignment; bearing parallel misalignment running model was built, and an analysis and comparison of various conditions were used to determine the influence, showing that the parallel misalignment has significant influence on the spindle Z direction thermal error. Keywords Machine tool, spindle, thermal error, bearing, parallel misalignment
Date received: 12 September 2017; accepted: 30 April 2018 Handling Editor: Xichun Luo
Introduction As the key component of the machine tool, the spindle thermal error influences the working accuracy of a machine tool positively.1 Due to high operating speeds, bearing characteristics are significantly affected by the assembly tolerances, generating detrimental effects on the spindle thermal characteristics and leading to thermal deformation of the shaft. Moreover, the lifespan of a bearing is greatly affected by the distribution and variation of the dynamic loads and contact pressures experienced during the running process.2 The bearing installation misalignment is a critical factor that causes concentrated contact loads between the balls and rings, resulting in only a few balls against the cutting force. Under the condition of bearing existing installation misalignment, the radial clearance has a non-uniform distribution, as shown in Figure 1. In addition, when a load is applied to the bearing, the contact angles and loads at different ball positions will
also have a non-uniform distribution. At present, it has been paid extensive attention in determining the influence of installation misalignment on bearing performance. Teutsch and Sauer3 proposed an improved slicing technology to present a more accurate pressure distribution for tilted bearings using the Lundberg logarithmic profile. TJ Park4 discussed the effect of deflection on contact pressure under elastohydrodynamic lubrication (EHL) and reported very high values due to the tilted misalignment. Zamponi et al.5
1
School of Mechatronics Engineering, Henan University of Science and Technology, Luoyang, China 2 School of Mechanical and Electrical Engineering, Wuhan University of Technology, Wuhan, China Corresponding author: Lihai Chen, School of Mechatronics Engineering, Henan University of Science and Technology, Luoyang 471000, China. Email:
[email protected]
Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).
2
Advances in Mechanical Engineering
Figure 1. The radial clearance distribution under parallel misalignment.
presented a mixed methodology employing both a finite element analysis (FEA) method and Hertz contact theory to analyze the contact pressure in bearings and considered the effects of tilted misalignment on the load distribution for bearings with both rigid and flexible rings. Z Ye et al.6 analyzed the distribution of loads and contact stress in high-speed roller bearings with quasi-dynamic method and FEA method considering effects of tilted misalignment between inner and outer rings for the roller cylindrical bearings. Interestingly, the influence of installation misalignment has rarely been studied for angular contact ball bearings, which are widely used in machine tool spindles driven by computer numerical control (CNC). Heat generated by the bearing is the main heat source within spindle systems. Moreover, transfer of heat from the bearing to the shaft has a considerable influence on thermal deformation of the shaft; therefore, the bearing temperature field is the intermediate bond between the bearing generated heat and the shaft thermal deformation. Y Li et al.7 discussed the relationship between the temperature field and thermal deformation of the spindle system by FEA. Xiang et al.8 built preliminary theoretical models of the temperature field and the thermal deformation based on the size of the spindle and the parameters of the bearing. The centrifugal force and thermal expansion occurring on the bearings and motor rotor change the thermal characteristics of the built-in motor, bearings, and assembly joints.9 And also many researchers have done the spindle thermal error analysis on the basis of bearing statuses calculated accurately.10–12 The bearing temperature field is the internal factor for the shaft thermal deformation. Hence, as an influencing factor of the bearing temperature field, parallel misalignment plays an important role in spindle thermal error analysis. In this article, the influence law of parallel installment misalignment on the internal physical relationships of the bearing and spindle thermal error was analyzed in detail. And an
Figure 2. Simplified schematic diagram of angular contact ball bearing.
action mechanism was proposed for bearing installation parallel misalignment.
Bearing running model with parallel misalignment A simplified schematic diagram of the angular contact ball bearing is shown in Figure 2. The radial clearance between the rings was assumed to move the shaft along the radial direction without shaft tilt; therefore, parallel misalignment of the bearing induced an additional radial displacement for each ball drcoscj with respect to the initial position. When the cutting force was applied to the bearing, the ball gained an additional radial displacement d0r cos cj , as shown in Figure 3. In Figure 3, B is the outer raceway groove curvature center, which is fixed in space, A is the inner raceway groove curvature center moving relative to the fixed center, and O is the ball center. The initial contact angles between the ball and rings were assumed nonuniform and the outer ring curvature center assumed fixed for the ball bearing using quasi-static analysis.13,14 Moreover, the inner ring curvature center moved from A to A# under the parallel misalignment dr and then from A# to A$ with the influence of an external load. The auxiliary parameters A1j, A2j, X1j, and X2j and the equilibrium equations for the entire bearing were defined as A1j = BD sin a0 + da A2j = BD cos a0 + dr cos cj + d0 r cos cj
ð1Þ
X2j X1j , sin aoj = ðfo 0:5ÞD + doj ðfo 0:5ÞD + doj A2j X2j A1j X1j cos aij = , sin aij = ðfi 0:5ÞD + dij ðfi 0:5ÞD + dij cos aoj =
ð2Þ
Dong et al.
3
Figure 3. The relative position for the bearing rings’ curvature centers. Figure 4. Ball loads at angle cj.
where f is the ring groove curvature radius coefficient, D is the bearing outer diameter, and cj is the ball bearing position angle. BD= (fo + fi2 1)D, which is the distance between the outer ring curvature center and the inner ring curvature center. The force applied to every ball is depicted in Figure 4. The force balance equations for each ball are presented as follows
Fr
Z X j=1
lij Mgi sin aij cos cj = 0 ð6Þ Qij cos aij + D Fcj = 2:26 3 1011 dmD3 nmj
ð7Þ
Mgj = 4:5 3 1012 nbj nmj sin bj
ð8Þ
where Fa and Fr are the applied loads on the bearing nbj = dm=D n
(1 D cos aij =dm)(1 + D cos aoj =dm) (1 D cos aij =dm) cos (ao bj ) + (1 + D cos aoj =dm) cos (ai bj )
n (1 D cos aij =dm) cos (aoj bj ) (1 + D cos aoj =dm) cos (aij bj ) + (1 D cos aij =dm) cos (aoj bj ) sin aoj D ,r= bj = cos aoj + r dm
nmj =
Qij cos aij Qoj cos aoj +
Mgj (lij sin aij loj sin aoj ) + Fcj = 0 D
ð3Þ Qij sin aij Qoj sin aoj
Mgj lij cos aij loj cos aoj = 0 ð4Þ D
1:5 Kij=oj = 2:15 3 105 where Qij = Kij d1:5 ij , Qoj = Koj doj , P 1=2 ri=j ðnd Þ3=2 . nd is the auxiliary parameter,
which is calculated by the function with principal curvature difference F(r).15 Combining the entire bearing balance equations (5) and (6) and the calculated equation for both centrifugal force and gyroscopic moment equations (7) and (8), the bearing internal contact angles and loads are determined by using the Newton–Raphson method Fa
Z X j=1
lij Mgj cos aij = 0 Qij sin aij D
ð5Þ
The specific parameters of a 7012AC bearing are shown in Table 1. Contact angles and loads were calculated relative to different positions using Fa= 400 N and Fr= 0 N and shown as Figures 5 and 6, respectively. In Figure 5, a huge difference among the contact angles at different ball positions can be observed and Figure 5 suggests that a larger misalignment may lead to a greater difference. From Figure 6, it can be observed that only a few balls support the applied load. This suggests that there were even fewer balls supporting the load and huge differences among the contact loads at different ball positions in the condition of larger misalignments.
Bearing heat generation distribution In machining process of machine tool operations, major heat sources include the heat generated by the cutting process and the heat produced by the bearings.
4
Advances in Mechanical Engineering
Table 1. Parameter values of the bearing 7012AC. Parameter name
Value
Parameter name
Value
Parameter name
Value
d (mm) Z a0
60 19 25
D (mm) Db ðmmÞ fo
95 10.7 0.525
Bearing width (mm) dm (mm) fi
18 77.5 0.515
Figure 5. Contact angles between the ball and ring.
Figure 6. Contact loads between the ball and ring.
It is assumed that the majority of cutting heat is taken away by coolant, and therefore, the heat generated by the bearings is the dominant factor of thermal deformation. The frictional torque M is the sum of three torques: (1) torque due to applied load, M1; (2) torque due to viscous friction, M0; and (3) torque due to spinning motion at the contact area, Msi(o). The torque due to the applied load M1 and the torque due to viscous friction M0 have previously been determined in the literature.13,14,16 The spinning friction torque Msi(o) is generated by ball spinning motion and has previously been defined13 Ms(i=o) =
3m Qi=o ai=o Ei=o 8
ð9Þ
From equation (9), it can be concluded that the spinning friction torque is proportional to the contact loads
and contact ellipse semi-major axis, which is closely connected to the bearing installation offset. Based on the calculated results of the contact loads and angles, the specific heat generated at different contact areas by the spinning friction torque under different axial force is calculated; axial force Fa = 400 N is shown in Figure 7. The generated heat shows the same trend with the contact angles and loads.
Simulation, experiment, and discussion The experimental spindle system proposed in this study is shown in Figure 8. The fine tooth thread alters the space within the oilfilled chamber and moves the rear bearing support backwards to apply a uniform axial load on the
Dong et al.
5 The spindle thermal error measurement system is shown in Figure 10. The system is symmetrical in the Z direction, such that the thermal error is smallest in the Z direction; therefore, the misaligned bearing supports can be installed parallel to the Z direction to verify the influence of bearing outer ring misalignment on the spindle thermal error. The spindle system thermal error, with deviations, was generated by an FEA model, based on a model previously described in our laboratory15 and Jiang and Mao.16 The comparative results for the simulation versus empirically derived measurements in condition of spindle speed of 3000 r/min and an axial load of 800 N are presented in Figure 11, and the axial load 800 N is selected for limiting the ball skidding which is clarified in Dong et al.17 Normally, the Z direction thermal error should have the smallest slope, because the spindle system has absolute symmetry structure. However, it was observed that the parallel bearing installation misalignment has a greater influence on the spindle thermal error in the Z direction than in the X and Y directions. It is indicated that the misalignment has changed the uniform temperature field on the peripheral direction, giving rise to the Z-direction thermal error, which corresponded with the previous analysis.
Figure 7. Specific heat generated at different contact areas.
bearings. A pressure sensor was used to measure the pressure of oil, which can be used for calculating the axial load indirectly. By changing the front bearing support to create different central axis deviations, various bearing parallel misalignments were simulated, as shown in Figure 9.
Figure 8. Experimental spindle system.
Conclusion In this article, the effect of bearing installation misalignment on the bearing characteristics and thermal error was analyzed for the first time. The condition for a parallel bearing installation misalignment was simulated through changing the bearing support center axis position, and the contact angles and loads were calculated based on a quasi-static analysis with reasonable
6
Advances in Mechanical Engineering Acknowledgements The authors would like to acknowledge the contributions from all collaborators within the projects mentioned.
Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Figure 9. Bearing support with different central axis deviations.
Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by the National Natural Science Foundation Committee (NSFC) of China (grant no. 51475343) and International Science & Technology Cooperation Program of China (grant no. 2015DFA70340).
ORCID iD Yanfang Dong
https://orcid.org/0000-0001-9834-7157
References
Figure 10. Spindle thermal error measurement system.
Figure 11. Comparison of results for the simulation versus empirical measurements.
assumptions. From the analysis results, the bearing misalignment was shown to have a considerable influence on the spindle thermal error as determined using both the FEA model and empirical measurements from the representative system. In conclusion, it is shown that if the bearing misalignment is properly optimized, spindle thermal error can be significantly reduced, and this work needs to be carried out later.
1. Abele E, Altintas Y and Brecher C. Machine tool spindle units. CIRP Ann: Manuf Techn 2010; 59: 781–802. 2. Harsha SP. Nonlinear dynamic analysis of a high-speed rotor supported by rolling element bearings. J Sound Vib 2006; 290: 65–100. 3. Teutsch R and Sauer B. An alternative slicing technique to consider pressure concentrations in non-Hertzian line contacts. J Tribol: T ASME 2004; 126: 436–442. 4. Park TJ. Effect of roller profile and misalignment in EHL of finite line contacts. In: Proceedings of the ASME 2010 10th biennial conference on engineering systems design and analysis, Istanbul, 12–14 July 2010, pp.395– 401. New York: IEEE. 5. Zamponi L, Mermoz E, Linares JM, et al. Impact of geometrical defects on bearing assemblies with integrated raceways in aeronautical gearboxes. Mech Mach Theory 2009; 44: 1108–1120. 6. Ye Z, Wang L, Gua L, et al. Effects of tilted misalignment on loading characteristics of cylindrical roller bearings. Mech Mach Theory 2013; 69: 153–167. 7. Li Y and Zhao W. Axial thermal error compensation method for the spindle of a precision horizontal machining center. In: Proceedings of 2012 IEEE international conference on mechatronics and automation, Chengdu, China, 5–8 August 2012. New York: IEEE. 8. Xiang S, Zhu X and Yang J. Modeling method for spindle thermal error based on mechanism analysis and thermal basic characteristics tests. Proc IMechE, Part C: J Mechanical Engineering Science 2014; 50: 144–152. 9. Chen JS and Hsu WY. Characterizations and models for the thermal growth of a motorized high speed spindle. Int J Mach Tool Manu 2003; 43: 1163–1170. 10. Bossmanns B and Tu JF. A thermal model for high speed motorized spindles. Int J Mach Tool Manu 1999; 39: 1345–1366.
Dong et al. 11. Haitao Z, Jianguo Y and Jinhua S. Simulation of thermal behavior of a CNC machine tool spindle. Int J Mach Tool Manu 2007; 47: 1003–1010. 12. Holkup T, Cao H, Kolar P, et al. Thermo-mechanical model of spindles. CIRP Ann: Manuf Tech 2010; 59: 365–368. 13. Jin C, Bo Wu and Hu Y. Heat generation modeling of ball bearing based on internal load distribution. Tribol Int 2012; 45: 8–15. 14. Ma C, Yang J, Zhao L, et al. Simulation and experimental study on the thermally induced deformations of highspeed spindle system. Appl Therm Eng 2015; 86: 251–268.
7 15. Dong Y, Zhou Z and Liu M. A general thermal model of machine tool spindle. Adv Mech Eng 2016; 8: 1–10. 16. Jiang S and Mao H. Investigation of variable optimum preload for a machine tool spindle. Int J Mach Tool Manu 2010; 50: 19–28. 17. Dong Y, Zhou Z and Liu M. Bearing preload optimization for machine tool spindle by the influencing multiple parameters on the bearing performance. Adv Mech Eng 2017; 9: 1–9.