Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 14, Article number: 711 (2024 ) Cite this article journal bearing
In this paper, on one hand, the time-varying characteristics of the heat source and thermal boundary conditions of the high-speed spindle system were analyzed considering the thermal-structural coupling effect. And a transient bearing temperature field prediction method combining the thermal network method and finite element method was proposed. Furthermore, the relationships between time step, calculation efficiency and calculation results were analyzed. On the other hand, a online real-time monitoring system of the transient temperature of the cylindrical roller bearing inner ring for maximum speed of 13,000 r/min was designed and implemented using fibre optic sensing technology. Comparing with the conventional static thermal analysis results, it is verified that the simulation method proposed in this paper has higher accuracy. This paper provides a new approach for analysing and testing the thermal characteristics of high-speed spindle system.
The spindle system is the key unit of CNC machines. The thermal performance of the system directly affects the machining accuracy. Therefore, it is important for analyzing the thermal characteristics of the spindle system components, especially the bearings which endure the cutting force, high rotation speed and generate a lot of heat. Currently, the thermal network and finite element methods are the common thermal characteristics analysis methods for bearings. Su et al.1,2 designed a thermal network including the lubricating oil and contact factors, and predicted the temperature rise of motorized spindles. Liu et al.3 proposed a thermal resistance network model of a spindle-bearing-bearing pedestal, and investigated the thermal characteristics of the motorized spindle system. Yan et al.4 studied the transient characteristics of the spindle-bearing system by the thermal network coupling with the thermal-deformation. Liu et al.5 established a thermal resistance network model and predicted the temperature field distribution in the spindle system with reasonable accuracy. Li et al.6 used the thermal network method in the simulation of power output performance of thermoelectric generation on a rotating spindle. In the meantime, a large number of scholars have reported the thermal analysis of the bearing using finite element methods. Liu et al.7 and Li8 predicted thermal distribution of motorized spindles by FEM technique. Gao et al.9 studied the coupled vibrations of motorized spindle using FEM. Zivkovic et al.10 presented 3-D finite element thermal model for studying the machine tool spindle thermal characteristics.
The above research has made a comprehensive and in-depth study on the heat transfer mechanism, heat-solid coupling behavior and analysis methods of high-speed spindle bearing systems using the thermal network and finite element methods. Thermal network method analysis computes thermal solutions at the nodes; however, it cannot accurately compute the temperature field of the actual structure, particularly at the internal portion of complex objects. The finite element method solves this shortcoming by providing solutions at all locations irrespective of the size and shape of the object. However, FEM can only obtain the results under a certain operating condition without considering boundary condition dynamic changes. In view of this, considering the thermal-structural coupling effect, this paper analyzed the time-varying characteristics of the thermal boundary conditions of the high-speed spindle system in operation. A transient bearing temperature field prediction method combining the thermal network method and finite element method was proposed. The relationship between time step and computation accuracy was reported. And a online real-time monitoring system of the transient temperature of the cylindrical roller bearing inner ring for maximum speed of 13,000 r/min was designed and implemented using fibre optic sensing technology. Finally, the proposed model was verified with the test results of the online transient temperature monitoring system.
The bearing is the main heat source of the spindle system, and the heat generation of the bearing is related to its speed, force, structure, lubrication and other conditions. The system's temperature field and thermal deformation interact through the heat source and heat dissipation boundary conditions, and finally the high-speed spindle system reach a thermal equilibrium state. The spindle system studied in this paper is shown in Fig. 1.
Schematic diagram of the spindle system (1. Shaft end caps 2. Sensor fixing plate 3. Coupling 4. Left oil injection nozzle 5. Test bearing).
The temperature rise in bearing is due to the friction between the contact surfaces during bearing rotation, including the sliding friction heat rate between the roller and the ring Q1, the roller churning friction heat rate Q2, the sliding friction heat rate between the cage and the ring surface Q3, the sliding friction heat rate between the roller and the cage pocket Q4, the elastic hysteresis loss heat rate Q5, rolling and the sliding friction heat rate Q6 between the roller end face and the retaining edge. The frictional heat rate between each part is calculated based on the elastic flow lubrication contact, and the detailed process can be referred to the literature11, so the total heat rate of the bearing Q is as follows:
The viscosity-temperature characteristics of the lubricant have a considerable effect on the heat generation. Therefore, it was considered during subsequent iterations of the calculation. And the physical parameters of the No. 46 lubricant used in this study are shown in Fig. 2.
Variation of lubricant physical parameters with respect to temperature.
The circumferential heat transfer was not considered as the system symmetry along the centre line of the spindle. The thermal nodes are divided according to the test spindle structure and the material characteristics of each portion, as shown in Fig. 3.
Schematic representation of test spindle structure and identification of thermal node arrangements.
The effect of thermal radiation can be ignored as the study was performed below 200 °C. The heat conduction and thermal convection was considered. The formulation of considered thermal resistances has been given in Table 111.
The thermal network calculated the transient temperature considering the heat flow at node i equal to the increase in energy within a specific volume13.
where, qi, ρi, Cρi, Vi and Ti represent the heat flow, material density, specific heat, corresponding volume and temperature at the i-th node, respectively, dTi/dt is the rate of temperature rise at the i-th node.
For example, assuming that heat Q0 from node 0 flows into neighbouring nodes 1, 2, 3 and 4, the transient heat balance equation at that location can be expressed as
where, Ti is the transient temperature of node i, and R0−i is the thermal resistance between the neighbouring nodes.
Based on the above analysis, the general procedure for solving the transient temperature field of the spindle system based on thermal network method is presented in Fig. 414. And the convergence criterion of this model is: until the temperature difference between the bearing balls at two adjacent time points is less than 0.1 ℃, the closed-loop iterative modelling is terminated, and then the results are outputted when the model of the thermal-structural interaction properties is converged. Under a radial force of 400 N and speed of 9000 r/min (used for all subsequent analyses), the transient characteristics of the thermal parameters are shown in Fig. 5.
The flowchart for computing the high-speed spindle thermal structure coupling process.
Variation of thermal parameter transient characteristic with respect to time.
Considering the time cost of the calculation, the 3D model needs to be simplified by ignoring secondary structures, such as springs and nuts which have little influence on the heat transfer of the system. In addition, to ensure the calculation accuracy, the mesh grid of the bearing area is to be refined. The results of the mesh division are shown in Fig. 6, with total of 277,328 cells, 1195,256 nodes and an average mesh mass of 0.9. The material properties of the system components are shown in Table 2.
Finite element model of the spindle bearing system.
The thermal boundary conditions are set as follows.
Loading of heat generation in the form of Heat Flow onto the roller, inner and outer raceway surfaces.
Natural air convection was applied to the housing and stationary surfaces. Forced air convection was applied to the rotating surfaces. Forced convection of lubricant was applied to the inner and outer raceway, roller surfaces. Temperature Dependent was selected as the setting method.
Meanwhile, a novel aspect of this paper is that the dynamic variation of the thermal parameters is considered. Therefore, choosing an optimal fitting function for thermal loading is required. Taking the dynamic power loss in Fig. 5, the comparison of the different curve fitting methods was made, as given in Table 3.
The goodness of fit is inversely proportional to the sum of squared errors SSE and root mean square error RMSE, directly associated with the coefficient of determination R2 and the corrected coefficient of determination Adjusted R2. According to the results of the comparison, the 8th-order Gaussian fitting function as Eq. (4) is chosen to simulate the dynamic changes of the heat source. The parameters in the expression are given in Table 4.
Firstly, the effect of iteration step on temperature rise is analyzed. The transient temperature rise of the bearing outer ring surface with respect to iteration step under ambient temperature of 19.5 ℃ is shown in Fig. 7.
Effect of iterative step size on transient temperature rise of the surface of the outer ring.
As can be seen from Fig. 7, the larger the iteration step of the simulation, the faster the temperature rise curve converges (the automatic time step of the program control is almost always greater than 50 s), and the faster the system reaches thermal equilibrium. The iteration step only affects the initial convergence rate of the temperature rise curve and does not affect the steady-state temperature value. This is due to smaller step sizes need more accurate modelling of the system dynamics. However, smaller step sizes lead to an increase in the time cost of the algorithm execution, as more iterative steps need to be performed to reach the desired time point. Conversely, larger step sizes may fail to capture the details of the system's dynamic changes, leading to compromised numerical accuracy. Therefore, the actual analysis process needs to be combined with the actual to set the iterative step size.
Secondly, the comparison of temperature rise under dynamic and static thermal boundary conditions has been done. Figure 8 showed the variation of temperature rise with static and dynamic boundary conditions at different time steps. At the initial stage, the trend of temperature rise under the two boundaries is consistent. However, the final values results were significantly higher. This was due to the static thermal parameters themselves that did not contain the thermal-structural coupling effects of the system. For example, the viscosity of the lubricant decreases significantly decreases with the increase of temperature, and the bearing power loss is variable with time. Therefore, the constant power losses and heat transfer coefficients calculated under the initial operating conditions were substituted. It also shows that the temperature prediction method proposed in this paper will be more accurate than the static boundary. On the basis of the above analysis, the simulation of bearing temperature field has been completed.
Comparison of static and dynamic thermal parameter simulation results.
The temperature rise test platform was required to verify the accuracy of the previously described model and the proposed simulation method. According to the structure of the spindle system used in this study, a bearing inner ring real-time temperature online monitoring system based on fibre optic sensing and rotary connector (FORJ) was designed (Fig. 9).
Schematic diagram of the temperature test system for bearing inner ring based on FBG.
The center wavelength drift of the fibre grating can be expressed as:
where, ΔT and ε denotes the temperature and strain change. Equation (5) is derived without considering the cross-sensitivity. The temperature sensitivity coefficient is KT = 7.22 × 10–6/°C, and the strain sensitivity coefficient is Kε = 0.784.
To verify the accuracy of temperature measurement, the FBG was placed in the temperature control box for calibrating, and the results are shown in Fig. 10. There was a good linear relationship between the centre wavelength change of the FBG and the temperature rise, which proves FBG can perform the measurement of the bearing temperature rise.
Calibration results of fibre Bragg grating sensors.
Fibre Bragg grating is sensitive to strain and temperature at the same time. For obtaining accurate temperature rise of the bearing inner ring, it was necessary to avoid the effect of strain and the environmental influence, so the packaging method shown in Fig. 11 was designed.
The grating was placed on a thin-walled copper substrate with a groove (red part in Fig. 11), and the thermally conductive silicone grease was injected into the groove to fix the substrate during encapsulation. This structure can reduce the effect of the groove gap on heat transfer and use the silicone grease's viscosity to reduce the fibre's micro-strain due to centrifugation. The test result comparison between the packaged and unpackaged FBG sensor has also been done through the thermal model. The contact interface between the connector and the bearing inner ring end face were set as a constant temperature, and then the temperature change of the FBG surface was monitored. The simulation results are shown in Fig. 12.
Test result comparison between the packaged and unpackaged FBG sensor.
It is evident in Fig. 12 that the timeliness and accuracy of the packaged FBG temperature measurement are better than the unpackaged one.
In addition, to prevent the influence of oil mist and high temperature, the substrate is fixed with the C-3 high temperature resistant inorganic adhesive, as shown in Fig. 13. Noting that the end of the fibre should be kept free when fixing the outer cladding, to prevent the influence of stretching caused by the substrate thermal expansion.
End cap connector and FBG sensor.
After the above treatment, the FBG was used almost strain-free, so that the relative wavelength drift under the effect of temperature can be expressed as given in Eq. (6):
where, α and ξ are the thermal expansion coefficient and thermo-optical coefficient of the fibre core, respectively. Pε denotes the effective elasticity coefficient, and αs represents the thermal expansion coefficient of the substrate material.
The test platform and sensor arrangement are shown in Fig. 14. The specifications, test conditions and other parameters of the spindle system and bearings are shown in Table 5. The results of the simulation and experimental tests under different speeds are shown in Fig. 15.
Comparison of results under different conditions.
The results of the simulation and experimental tests at different speeds are shown in Fig. 15.
It can be seen that the method proposed in this paper has sufficient accuracy to simulate the temperature rise characteristics of the spindle bearing system, and the smaller the iteration step, the higher the simulation accuracy. With the decrease of iteration step, the simulation temperature rise is closer to the real temperature rise in the whole process.
In addition, reducing the iteration step length can get better simulation results. In this paper, the step length was reduced from Δτ = 10 s to Δτ = 5 s, but the computation times is increased by three times. In the meantime, it can also be seen that the iteration step has no effect on the final thermal balance temperature of the bearing, so the iteration step length should be chosen reasonably in practical engineering.
This paper comprehensively considers the thermal-structural coupling effect, analyses the time-varying characteristics of the thermal parameter conditions of the high-speed spindle system under the operating condition, combines the thermal network method with the finite element method for the prediction of the transient temperature field of the bearings, and designs a high speed bearing inner and outer ring temperature measurement test system for verification. The main conclusions are shown as follows.
A method for predicting the transient temperature field of bearings by combining the thermal network method with the finite element method is proposed, and it is demonstrated experimentally that the method has higher accuracy than the simulation using static thermal boundary conditions.
Online monitoring of the inner ring transient temperature of cylindrical roller bearings at high speed (maximum test speed 13,000 r/min) is designed using fibre optic rotary connectors and FBG sensors.
Shortening the iteration step can improve the simulation accuracy in the ramp-up phase using dynamic thermal boundary condition.
Through the study of this paper, it is evident that the bearing temperature prediction will be more accurate by using the dynamic thermal boundary, and the changes in boundary conditions cannot be ignored. Simultaneously, the designed online temperature monitoring system of bearing inner ring can be extended to the dynamic strain test of high speed rotating parts, such as gear, aero-engine blades.
The datasets used and analyzed during the current study available from the corresponding author on reasonable request.
Su, C. & Chen, W. F. An optimized thermal network model to evaluate the thermal behavior on motorized spindle considering lubricating oil and contact factors. Proc. Inst. Mech. Eng. Part C-J. Mech. Eng. Sci. 13(236), 7484–7499 (2022).
Su, C. & Chen, W. F. An improved model of motorized spindle for forecasting temperature rise based on thermal network method. Int. J. Adv. Manuf. Technol. 119(9–10), 5969–5991 (2022).
Liu, Z. F. et al. Thermal characteristic analysis of high-speed motorized spindle system based on thermal contact resistance and thermal-conduction resistance. Int. J. Adv. Manuf. Technol. 9–12(76), 1913–1926 (2015).
Yan, K., Hong, J., Zhang, J. H., Mi, W. & Wu, W. W. Thermal-deformation coupling in thermal network for transient analysis of spindle-bearing system. Int. J. Therm. Sci. 104, 1–12 (2016).
Liu, Y., Ma, Y. X., Meng, Q. Y., Xin, X. C. & Ming, S. S. Improved thermal resistance network model of motorized spindle system considering temperature variation of cooling system. Adv. Manufact. 6(4), 384–400 (2018).
Li, S., Yao, X. H. & Fu, J. Z. Power output characterization assessment of thermoelectric generation in machine spindles for wireless sensor driving. Sens. Rev. 34(2), 192–200 (2014).
Liu, M. Y. et al. Study of the temperature distribution of a machine tool spindle bearing based on FBG quasi-distributed sensing. Int. J. Adv. Manuf. Technol. 98(1–4), 263–274 (2018).
Li, S. S., Shen, Y. & He, Q. Study of the thermal influence on the dynamic characteristics of the motorized spindle system. Adv. Manufact. 4(4), 355–362 (2016).
Gao, F., Cheng, M. K. & Li, Y. Analysis of coupled vibration characteristics of PMS grinding motorized spindle. J. Mech. Sci. Technol. 34(9), 3497–3515 (2020).
Zivkovic, A. M. et al. A study of thermal behavior of the machine tool spindle. Thermal Sci. 23(3), 2117–2130 (2019).
Hao, X. Research on Thermal Characteristics of Cylindrical Roller Bearing and Its Effect on Vibration Response of Rotor System (Dalian University of Technology, 2020).
Takabi, J. & Khonsari, M. M. Experimental testing and thermal analysis of ball bearings. Tribol. Int. 60, 93–103 (2013).
Harris, T. A. & Kotzalas, M. N. Advanced Concepts of Bearing Technology: Rolling Bearing Analysis 5th edn. (CRC Press, 2007).
Liu, J. et al. Thermal-structure interaction characteristics of a high-speed spindle-bearing system. Int. J. Mach. Tool. Manu. 137, 42–57 (2019).
This research was funded by the National Natural Science Foundation of China, Grant number 52205096, the Key Research Projects for Higher Education Institutions of Henan Provincial Education Department, Grant number 23A460007.
School of Mechatronics Engineering, Henan University of Science and Technology, Luoyang, 471003, China
Yanfang Dong, Yafei Ma, Ming Qiu, Feifan Chen & Kai He
Henan Provincial Cooperative Innovation Center for Advanced Manufacturing of Mechanical Equipment, Luoyang, 471003, China
Yanfang Dong & Ming Qiu
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Methodology and writing—original draft preparation, data processing and mining, editing, Y.D. and Y.M.; M.Q., supervision; F.C. and K.H., data checking. All authors have read and agreed to the published version of the manuscript.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Dong, Y., Ma, Y., Qiu, M. et al. Analysis and experimental research of transient temperature rise characteristics of high-speed cylindrical roller bearing. Sci Rep 14, 711 (2024). https://doi.org/10.1038/s41598-024-51255-9
DOI: https://doi.org/10.1038/s41598-024-51255-9
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
insert bearings Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.