Experimental and numerical analyses of the displacements and stiffnesses of a needle roller bearing under different conditions | Scientific Reports

Feb 23, 2025

Scientific Reports volume 15, Article number: 6364 (2025) Cite this article

Metrics details

The displacement and stiffness of a roller bearing affect the reliability and remaining life of the roller bearing. In order to establish a method for predicting the displacement and stiffness of a needle roller bearing, experimental and numerical analyses were performed on an NA4848 needle roller bearing under different conditions. The tested displacement varies from 0.0029 mm to 0.0125 mm for the radial load in a range of 10-50kN. Test results also showed that the effect of the tested rotational speeds on the radial stiffness is not significant. The minimum stiffness based on the tested displacement occurs when the radial loading is 20kN. Finite element analysis was performed, from which the displacements and von Mises/contact stress distributions of the NA4848 needle roller bearing under different conditions can be obtained. The FE predicted relative displacements and stiffnesses are compared with the analytical calculations and experimental results, which showed that the predicted values of the relative displacements were higher than the analytical calculations and experimental results by 0.009–0.018 mm, and the predicted values of the stiffnesses were higher than the analytical calculations and experimental results by 36 × 106-164 × 106 N/mm. In addition, it was found that the predicted values of the finite element method for the contact stresses were lower than the analytical calculations by about 400–600 MPa.

Roller bearings, which possess the advantage of high load carrying capacity and the ability to accommodate high rotational speeds, are widely used in rotating machinery. The stiffness characteristics of roller bearings are essential to reduce noise, heat, energy consumption, and wear. In practice, the stiffness of roller bearings can be influenced by many factors, such as rotational speed, clearance, and temperature. It is therefore necessary to ensure that the accurate stiffness of the roller bearing is estimated under various influencing factors to improve the reliability of the rotating machinery and determine the remaining life of the needle roller bearing.

Early studies on the stiffness of rolling element bearings, which were carried out by many researchers1,2,3, established different nonlinear relationships between the deformation and the external load for roller element bearings. For instance, Harris and Kotzalas3 proposed a power-law relationship between the deformation and the external load for a roller bearing. Recently, Demirhan and Kanber4 numerically investigated the stress and displacement distributions on a cylindrical roller bearing by using a two-dimensional finite element (FE) model with the consideration of steel shaft-bearing system.

Wei et al.5 carried out FE analyses on the pre-loaded hollow cylindrical roller bearing to investigate the effect of load, interference magnitude, and hollow ratio on the radial stiffness of the bearing. Guo and Parker6 developed a three-dimensional FE/contact mechanics model to determine the fully-populated stiffness matrix corresponding to two radial, one axial, and two angular coordinates of rolling element bearings. Chen and Wang7 proposed a method to study the radial stiffness of a cylindrical roller bearing with a corrected roller generator, which was derived by coupling this load-approach function into a complete mechanical model of the bearing. Hao et al.8 developed a finite element model for obtaining the displacement, raceway stress, and bearing stiffness characteristics of the inner and outer rings of the N1013M cylindrical roller bearing mounted in a bearing housing. Hao and his colleagues9 also experimentally investigated the time-varying stiffness of the cylindrical roller bearing due to thermal behavior, i.e. the warm-up process and different lubricant oil temperatures. Han et al.10 developed an analytical method for calculating the radial stiffness, load distribution, and displacement of the rotor center for the radially loaded ball and roller bearings, with consideration of two boundary positions of the inner ring supported on an even or odd number of rolling elements. Ri et al.11 investigated a numerical method for creating a comprehensive stiffness model of a cylindrical roller bearing, taking into account the radial clearance. The model was based on the Hertz contact and elastohydrodynamic lubrication theories. Hu et al.12 established a four-degree-of-freedom quasistatic model based on the basic analysis theory of rolling bearings. They analyzed the deformation, load distribution, stiffness, and fatigue life of an NU1014 cylindrical roller bearing with different roller profiles under various load and rotating speed conditions.

A failure analysis was carried out for a NA4848 needle roller bearing in a previous study13. In this study, experimental tests were carried out on a NA4848 needle roller bearing under different clearances, rotational speeds, and applied loads. The characteristics of the tested radial stiffness were analyzed. A finite element (FE) model of the NA4848 was established, and a mesh sensitivity study was performed to obtain mesh-independent FE results, such as relative displacements and contact stresses. The analytical and FE predictions were compared to the test results. The analytical predictions show good agreement with the test results for both relative displacements and stiffnesses. The finite element predicted values of relative displacements are 0.009–0.018 mm higher than analytically calculated and experimental results, and the finite element predicted values of stiffness are 36 × 106-164 × 106 N/mm higher than analytically calculated and experimental results. In addition, it is found that the predicted values of the finite element method for the contact stresses are about 400–600 MPa lower than the analytically calculated results. However, when taking the scatter into account, the FE analyses are likely to be used for predicting the remaining life of the NA4848 needle roller bearing in future studies.

The test rig consisted primarily of a servo-electric motor, a universal joint coupling, a reducer, a fixture, the loading system, and a data acquisition system, as can be seen in Fig. 1. The needle roller bearing being tested was installed on the test rig using the fixture. The tested needle roller bearing was installed on the test rig by using the fixture. A radial load, in the range of 10-50kN, and a rotational speed, in the range of 0-15 rpm, were applied to the tested needle roller bearing. Two different assembly clearances, i.e. 0.05 mm and 0.1 mm, were tested for the tested needle roller bearing, respectively. Two Bojke BLG-30 N-485 laser displacement sensors were used to monitor the displacements of the inner and outer rings of the tested bearing, i.e., 1 and 2. All tests were carried out at room temperature, and lubricating oil was used for lubrication and cooling.

Hao et al.8 experimentally identified the nonlinear stiffness characteristics of the bearing. They used the function obtained from curve-fitting the data of radial load versus relative displacement of the inner ring. The relative displacement of the inner ring is r = 2 - 1. The stiffness can be calculated as the reciprocal of the differentiated curve-fitting function, as shown in Eq. 1.

where, Kr is the radial stiffness of the bearing, with a dimension of N/mm, and \(\:\frac{d{\delta\:}_{r}}{d{F}_{r}}\) is the reciprocal of the relative radial displacement of the bearing to the radial load, with a dimension of 1.

The test rig.

Figure 2 shows the experimental relationship between the radial relative displacement, r, and the radial loading, Fr, at different rotational speeds and assembly clearances. It can be seen in Fig. 2 that the effects of the tested rotational speeds and clearances on the radial relative displacements are not significant. However, the radial relative displacements for the assembly clearance of 1.0 mm are slightly higher than those for the assembly clearance of 0.05 mm. Figure 3 shows the experimentally determined stiffnesses based on Eq. 1 versus radial loads at different rotational speeds and radial clearances. As shown in Fig. 3, the stiffness increases with the increase in radial load at different rotational speeds and assembly clearances.

Experimental relationship between the radial relative displacement, r, and the radial loading, Fr.

Experimentally determined stiffnesses versus radial loads at different rotational speeds and radial clearances.

A finite element model for the NA4848 was established using the commercial code ABAQUS 202214. The penalty function method was adopted for each contact pair, and the friction coefficient for each contact pair was taken as 0.1. Mesh sensitivity analyses were carried out using three types of mesh designs for the needle roller, as shown in Fig. 4. A general purpose linear brick element, with reduced integration (i.e. C3D8R element) was used to mesh the FE models15,16,17. The total number of C3D8R elements for the three designs were 80,640, 127,232, and 166,656, respectively. A radial load was applied to a central reference point on the inner ring, which was coupled to the inner surface of the ring. An encastre boundary condition was applied to the outer surface of the outer ring, constraining all degrees of freedom of the nodes on the bottom surface. A radial load of 50kN was applied to the FE model for the mesh sensitivity study. According to the test results, the rotational speed (in the range of 0-15 rpm) had no significant effect on displacement and stiffness. Therefore, a maximum experimental rotational speed of 15 rpm was used for all FE analyses. Figure 5 shows the contact pressure values corresponding to the number of elements. It can be seen in Fig. 5 that 127,232 C3D8R elements were found to be capable of giving sufficiently accurate results using relatively less calculation time.

Three types of mesh designs for the needle roller for total elements: (a) 80,640; (b)127,232; (c) 166,656.

Contact pressure values corresponding to the number of elements.

The FE-predicted von Mises stress distributions of each component of the NA4848 needle roller bearing with two radial clearances of 0.05 mm and 0.1 mm, a radial load of 50kN, and a rotational speed of 15 r/min are shown in Figs. 6 and 7, respectively. From Figs. 6 and 7, it can be seen that the von Mises stresses in the loaded needle roller are the highest with 428.5 MPa and 451.7 MPa for both clearances, followed by the outer ring with 194.2 MPa and 157.6 MPa for both clearances, the inner ring with 131.6 MPa and 155.4 MPa for both clearances, and the cage with the lowest stresses with 66.78 MPa and 53.69 MPa for both clearances, respectively. There is both rolling and sliding between the needle roller and the inner/outer rings, and the contact area is similar to an elongated ellipse.

The FE-predicted von Mises stress distributions of each component of the NA4848 needle roller bearing (clearances: 0.05 mm ): (a) needle roller; (b) outer ring; (c) inner ring; (d) cage.

The FE-predicted von Mises stress distributions of each component of the NA4848 needle roller bearing (clearances: 0.1 mm): (a) needle roller; (b) outer ring; (c) inner ring; (d) cage.

Figure 8 (a) and (b) show the von Mises stress versus rotational angle results with two assembly clearances for the nodes of each component of the NA4848 needle roller bearing that have the maximum von Mises stress. It can be seen in Fig. 8 that the levels of stress for the needle rollers are much higher than those for the other components of the NA4848 needle roller bearing for two radial clearances. The von Mises stress values of the rolling needle under two types of radial clearances are 530.453 MPa and 555.078 MPa, respectively. The von Mises stress values of the outer ring under two types of radial clearances are 198.687 MPa and 171.786 MPa, respectively. The von Mises stress values of the inner ring under two types of radial clearances are 138.316 MPa and 161.574 MPa, respectively. The von Mises stress values of the cage under two types of radial clearances are 71.4292 MPa and 70.8446 MPa, respectively.

Von Mises stress versus rotational angle results: (a) 0.05 mm radial clearance; (b) 0.1 mm radial clearance.

The contact stress distributions of the needle roller with two radial clearances of 0.05 mm and 0.1 mm, a radial load of 50kN, and a rotational speed of 15 r/min are shown in Fig. 9. It can be seen that the contact stress of the needle roller with a radial clearance of 0.05 mm is 59.76 MPa, which is lower than the 65.69 MPa with a radial clearance of 0.1 mm.

The contact stress distributions of the needle roller with two radial clearances: (a) 0.05 mm and (b) 0.1 mm.

Figure 10 show the von Mises stress and contact stress distributions along the needle roller total length under different radial clearances. It can be seen that the maximum stress occurs at the ends of the needle rollers at about 12 mm. The von Mises stress and contact stress profiles for the radial clearances of 0.05 mm are almost consistent with those for the radial clearances of 0.01 mm, except for the radial load is in a value of 30kN. Figure 11 shows the FE-predicted maximum von Mises/contact stresses versus radial loads under two radial clearances. As shown in Fig. 11, the relationship between the FE predicted maximum stress and radial load for the radial clearance of 0.1 mm is linear while that for radial clearance of 0.05 mm is nonlinear.

(a) von Mises stress and (b) contact stress distributions along the needle roller total length under different radial clearances.

FE-predicted maximum von Mises/contact stresses versus radial loads under two radial clearances.

Harris and Kotzalas3 suggested a correlation between the applied radial load, Fr, and the maximum radial deflection, δmax, as shown in Eqs. 2and 3.

Where Fr is the applied radial load on the needle roller bearing, Z is the number of the rollers in the needle roller bearing, Kn is the load-deflection factor, δr is the ring radial shift, Pd is the clearance of the bearing, Jr(ε) is the radial integral which is a function of load distribution factor, ε.

For the needle roller bearing with a given clearance and load, Eq. 2 can be solved by trial and error method, as demonstrated in Ref3. in detail. The value of parameter Kn in Eq. 2 is shown in Eqs. 4 and 5.

For steel roller-steel raceway contact,

where l is the effective length of the roller. The total roller lengths and logarithmic roller profiles for NA4848 were determined in our previous study13. Therefore, l can be taken as the straight portion of the crowned roller, with a value of 30 mm.

Figures 12 and 13 show the comparison between the tested and FE results of relative displacement and stiffness for different rotational speeds and radial clearances and those obtained from Harris and Kotzalas’s method, respectively. Multiply the finite element prediction of the maximum strain of the needle roller with the radius of the needle roller to obtain the finite element prediction of the relative displacement. It can be seen that the FE predicted the relative displacements are almost keep constant and larger than those obtained from test data and analytical results. The analytical calculated relative displacements show good agreement with test results for two different radial clearances. It can be seen that the predicted values of the FE method of stiffness are 36 × 106-164 × 106 N/mm higher than the analytical calculations and experimental results.

Comparison between the tested and FE results of relative displacement and stiffness for different rotational speeds and radial clearances and those obtained from Harris and Kotzalas’s method.

Comparison between the tested and FE results of stiffness for different rotational speeds and radial clearances and those obtained from Harris and Kotzalas’s method.

Based on Hertz’s contact theory18,19,20,21, the maximum contact stress between the needle roller and the bearing ring can be given by Eq. 6 and Eq. 7, respectively.

where Q is the normal force between rolling body and raceway; b is the short half-axis of projection touching ellipse; Σρ is the sum of the curvatures, \(\:\varvec{\varSigma\:}\varvec{\rho\:}=\frac{1}{{\varvec{R}}_{\varvec{i}\varvec{n}\varvec{n}\varvec{e}\varvec{r}}}+\frac{1}{{\varvec{R}}_{\varvec{r}\varvec{o}\varvec{l}\varvec{l}\varvec{e}\varvec{r}}}\)

Figure 14 shows the comparison between the analytical and FE results of contact stress for different radial clearances. Contact stresses predicted by the FE method are 400–600 MPa lower than those calculated analytically.

Comparison between the analytical and FE results of contact stress for different radial clearances.

This study proposed a numerical method for predicting the displacement and stiffness of a needle roller bearing under different conditions. Experimental and analytical validations were carried out to validate the numerical predictions. The concluding remarks are as follows:

(1) A testing system was established to obtain the displacement of the inner and outer rings of the tested bearing under different conditions. The experimental results show that the effect of the tested rotational speeds on the radial stiffness is not significant, and the minimum radial stiffness occurs when the radial loading is 20 kN.

(2) Finite element analyses were performed to predict the displacement and contact stress of the needle roller bearing under a load of 50kN and a rotational speed of 15 rpm. The FE results show that the maximum von Mises stress occurs on the needle rollers for both radial clearances. The von Mises stress and contact stress distributions along the total length of the needle roller show that the effect of radial clearance on the stress profile is insignificant, except for a radial load of 30kN. The relationships between the FE predicted maximum stress and radial load for radial clearances of 0.05 mm and 0.1 mm are linear and nonlinear, respectively.

(3) The FE and Harris and Kotzalas’s analytical predictions of the relative displacement and stiffness for different rotational speeds and radial clearances were compared with test results. Comparing with the test results, the FE results show significant overpredictions, while the analytical predictions show good agreement. The FE predicted contact stresses are lower than those obtained from analytical predictions.

(4) Due to the insufficient testing apparatus, the experimental tests were only carried out under relatively low rotational speeds, i.e., ≤ 15 rpm. Thus, we suggest considering an optimized testing method for high rotational speed conditions in future studies. The power loss and load distribution, which are essential for evaluating the remaining life of bearings, should also be investigated in the future.

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Palmgren Ball and Roller Bearing Engineering (3rd Edition) 34–41Burbank, Philadelphia, (1959).

Harris, T. A. Rolling bearing analysis (2nd Edition) (John Wiley & Sons, (1984).

Harris, T. A. & Kotzalas, M. N. Advanced Concepts of Bearing Technology: Rolling Bearing Analysis (5th Edition)CRC press, (2007).

Necdet, D. & Bahattin, K. Stress and displacement distributions on cylindrical Roller Bearing rings using FEM. Mechanics based design of structures and machines. Int. J. 36:1, 86–102. https://doi.org/10.1080/15397730701842537

Wei, Y., Liu, Y. & Zhang, X. Theoretical Research on the Radial Stiffness of Pre-Loaded Hollow Cylindrical Roller Bearings. Key Engineering Materials (Trans Tech Publications Ltd, 2011). https://doi.org/10.4028/www.scientific.net/kem.464.362

MATH Google Scholar

Parker, G. R. G. Stiffness matrix calculation of rolling element bearings using a finite element/contact mechanics model. Mechanism and Machine Theory: Dynamics of Machine Systems Gears and Power Trandmissions Robots and Manipulator Systems Computer-Aided Design Methods. 51, (2012).

Chen, G. & Wang, H. Contact stress and radial stiffness of a cylindrical roller bearing with corrected roller generator. Trans. Can. Soc. Mech. Eng. 40, 725–738. https://doi.org/10.1139/tcsme-2016-0059 (2016).

Article MATH Google Scholar

Hao, X., Gu, X., Zhou, X., Liao, X. & Han, Q. Distribution characteristics of stress and displacement of rings of cylindrical roller bearing. Proc. Institution Mech. Eng. Part. C: J. Mech. Eng. Sci. 233 (12), 4348–4358. https://doi.org/10.1177/0954406218820551 (2019).

Article MATH Google Scholar

Hao, Y. & H, Q. Time-varying stiffness characteristics of roller bearing influenced by thermal behavior due to surface frictions and different lubricant oil temperatures. Tribol. Int. 144 (2020).

Han, Y., Yang, L. & Xu, T. Analysis of static stiffness fluctuation in radially loaded ball and roller bearings. Arch. Appl. Mech. 91 (4), 1–16. https://doi.org/10.1007/s00419-020-01853-6 (2021).

Article MATH Google Scholar

Ri, C. U. et al. Comprehensive stiffness analysis of the cylindrical roller bearing for aircraft engines considering the radial clearance. Aircraft engineering and aerospace technology. (2023). https://doi.org/10.1108/AEAT-01-2022-0026

Hu, X. et al. Mechanical Properties Analysis of Cylindrical Roller Bearings with different profiles under combined load. ASTM Int. J. Test. Eval. 52 (1), 417–445. https://doi.org/10.1520/JTE20230230 (2024).

Article MATH Google Scholar

Li, R. et al. Failure analysis of a needle roller bearing in a Megawatt Reciprocating Pump. J. Fail. Anal. Preven. 24, 108–115. https://doi.org/10.1007/s11668-023-01821-6 (2024).

Article MATH Google Scholar

Abaqus Online documentation. Dassault Systèmes;2022. (2022).

Hyde, T. H., Saber, M. & Sun, W. Testing and modelling of creep crack growth in compact tension specimens from a P91 weld at 650°C. Eng. Fract. Mech. 77 (15), 2946–2957. https://doi.org/10.1016/j.engfracmech.2010.03.043 (2010).

Article MATH Google Scholar

Hyde, T. H., Li, R., Sun, W. & Saber, M. A simplified method for predicting the creep crack growth in P91 welds at 650℃. Proceedings of the Institution of Mechanical Engineers, Part L, Journal of Materials: Design and Application. (4), 224. (2010). https://doi.org/10.1243/14644207JMDA330

Li, R., Hyde, T. H. & Sun, W. Finite element prediction of fatigue crack growth in super CMV hollow shafts with transverse holes under combined torsional and axial loading. J. Strain Anal. Eng. Des. 48 (7), 457–469. https://doi.org/10.1177/0309324713492318 (2013).

Article MATH Google Scholar

Gawande, S. H., Balashowry, K., Raykar, K. A. & Munde, K. H. Estimation of contact stresses in EN31 rolling contact bearings for screw compressor using Gauss quadrature & statistical analysis. Australian J. Mech. Eng. 1–17. https://doi.org/10.1080/14484846.2023.2195098 (2023).

Mitul, T. S. & Dipak, V. A finite element analysis of an elastic contact between a layered cylindrical hollow roller and flat contact. Industrial Lubrication Tribology. 69 (1). https://doi.org/10.1108/ILT-11-2015-0189 (2017).

Ji, P., Gao, Y., Ma, F. & An, Q. Influences of roller diameter error on contact stress for cylindrical roller bearing. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology. 229, 689–697. (2014). https://doi.org/10.1177/1350650114559617

Han, C., Yu, C., Li, Y., Yan, J. & Zhang, J. Mechanical performance analysis of hollow cylindrical roller bearing of cone bit by FEM. Petroleum 1 (4), 388–396. https://doi.org/10.1016/j.petlm.2015.10.012 (2015).

Article MATH Google Scholar

Download references

This research was supported by CCTEG Project (2023-TD-QN004, 2023-TD-QN006) and TMIC Project (2022TM-167 M) funding, China.

Graduate School, Coal Research Institute, 100013, Beijing, China

Ran Li, Fengcai Liu, Tianze Cao & Wenbo Li

State Key Laboratory of Intelligent Coal Mining and Strata Control, 100013, Beijing, China

Ran Li, Wenshu Wei, Hao Liu, Jian Ye, Yuehua Lai, Jianfeng Li & Haicheng Lu

Beijing Tianma Intelligent Control Technology Co., Ltd, 101399, Beijing, China

Ran Li, Wenshu Wei, Hao Liu, Jian Ye, Yuehua Lai, Jianfeng Li & Haicheng Lu

College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, 266590, Qingdao, China

Ran Li

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

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

Conceptualization, R. Li; methodology, R. Li; formal analysis, F. C. Liu; data curation, F. C. Liu, H. Liu, J. Ye, Y. H. Lai, J. F. Li, T. Z. Cao and W. B. Li; writing—original draft preparation, R. Li; writing—review and editing, H. Liu, supervision, W. S. Wei.; project administration, H. C. Lu; funding acquisition, R. Li, Y. H. Lai. All authors have read and agreed to the published version of the manuscript.

Correspondence to Yuehua Lai.

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-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.

Reprints and permissions

Li, R., Wei, W., Liu, F. et al. Experimental and numerical analyses of the displacements and stiffnesses of a needle roller bearing under different conditions. Sci Rep 15, 6364 (2025). https://doi.org/10.1038/s41598-025-87276-1

Download citation

Received: 06 November 2024

Accepted: 17 January 2025

Published: 21 February 2025

DOI: https://doi.org/10.1038/s41598-025-87276-1

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