Design and manufacturing of roller bearing polymeric cages and development of a theoretical model for predicting the roller push-out force | Scientific Reports

News

HomeHome / News / Design and manufacturing of roller bearing polymeric cages and development of a theoretical model for predicting the roller push-out force | Scientific Reports

Oct 14, 2024

Design and manufacturing of roller bearing polymeric cages and development of a theoretical model for predicting the roller push-out force | Scientific Reports

Scientific Reports volume 12, Article number: 1017 (2022) Cite this article 2577 Accesses 4 Citations 1 Altmetric Metrics details In the present work, polymeric cages with 18 different pocket

Scientific Reports volume 12, Article number: 1017 (2022) Cite this article

2577 Accesses

4 Citations

1 Altmetric

Metrics details

In the present work, polymeric cages with 18 different pocket geometries are developed to investigate the effects of geometrical parameters and material properties on the amount of roller push-out force. An experimental setup including a specialized injection molding tool is designed and fabricated and three sets of polymeric cages are manufactured using the selected materials (PA46, PA66, PPA). Force measurements are carried out five times on each pocket and three cages for each material are tested. Considering three different materials, a total of 810 force measurements are performed. A theoretical model is developed to predict the roller push-out forces in polymeric cages with different materials and pocket geometries. The model is developed by estimating the deformed region of the cage as a cantilever beam with a parabolic profile. An empirical coefficient is reposed in the model to compensate for the assumptions applied to the model. Experimental results showed that a fixed coefficient gives accurate results for all the geometries and materials, which confirms the validity of the approach adopted in this paper for modeling such problems. Considering the geometrical and material tolerances, force limits predicted by the model cover all the forces measured for a specific pocket with excellent accuracy and consistency.

Roller bearings are prominent components in rotating machinery, and their cost and performance directly affect the entire system. Typical roller bearings comprise an outer ring, an inner ring, balls or rollers as rolling elements, and a cage. The cage (also known as a bearing retainer or roller separator) is a component in roller bearings, and the primary functions of the cage are:

separating the rolling elements to reduce friction in the bearing

optimizing the load distribution by evenly spacing the rolling elements

guiding the rolling elements

keeping the rolling elements of separable bearings in position when a bearing ring is removed for mounting or dismounting

providing lubrication by a solid film (coating or cage material) or functioning as an oil reservoir

All cages produce frictional forces that increase the starting and running torque of the bearings. Therefore, cages are designed in various configurations and materials to improve the performance of the bearings in different applications such as high speed, torque-sensitive, and extreme environment applications1.

There are four types of cages, namely: stamped metal, machined metal, pin-type, and polymeric cages. The polymeric cages, among others, produce less friction due to their good sliding properties and enable high speeds. These cages lower the risk of seizure and secondary damage in poor lubrication conditions because they can survive with a low amount of lubrication for some time. Furthermore, polymeric cages can be complex in shape while being inexpensive. These cages are usually made of polyamide 46 (PA46), polyamide 66 (PA66), polyphthalamide (PPA), polyether ether ketone (PEEK), polyethylene (PE), and polytetrafluoroethylene (PTFE). These cages are usually made in very large quantities, so the costs of molds and process automation are very high. Therefore, various research studies have focused on prototype samples to understand or verify the effects of manufacturing processes, geometrical dimensions, and material properties on different functions of the cages. For instance, Gaydamaka2 determined the relationships between structural elements and their functions in roller bearings and studied the effects of increasing the number of rollers in the cage. Williams3 investigated a new retainer for cylindrical roller bearings made of fiberglass reinforced nylon. The cage characteristics set up as design goals in this study are the applicability of one retainer for various bearings, maximizing the number of rollers, being race or roller guided, minimizing the roller/cage drop, and low cost. Ghaisas et al.4 developed a six Degrees of Freedom (DOF) model to study cage instabilities in roller bearings considering the roller-race and roller-cage pocket clearances for light-load and high-speed conditions. The results show that the cage exhibits stable motion for small values of roller-race and roller-cage pocket clearances. Kohar et al.5 considered large, tapered roller bearing with flexible body cages to provide information about kinematic and dynamic relationships of steel and plastic cages under various operating conditions. Cui et al.6 presented the nonlinear dynamic differential equations of high-speed roller bearing considering the impact of roller dynamic unbalance. Ripanu et al.7 considered different manufacturing parameters on the precision of the windows main dimensions that exist in the bearing cages in the case of the double row tapered roller bearings. Empirical mathematical models able to offer information concerning the impact of manufacturing parameters on three main dimensions of the bearing cage windows were determined in this research. Cui et al.8 presented a method for measuring the imbalance in a small-sized cylindrical roller and considered the impact of the roller imbalance on the dynamic characteristics of a cage. Moreover, there are some studies on the static9 and dynamic10 numerical analysis of fiber-reinforced polymeric materials. The study presented here, focusing on the effects of material properties and geometrical parameters on the roller push-out forces in polymeric cages, addresses the gaps in current literature. This exploration is essential because the amount of push-out force should be kept in a certain range. This range is defined by compromising between the highest allowable force required to eject the mold inserts and the lowest acceptable force needed to keep the rollers in position. To answer this question, a polymeric cage with 18 different pocket geometries is designed to investigate the effect of four main geometrical parameters. An innovative tooling system is designed and fabricated, enabling low-volume production of polymeric cages. Next, cages with different materials (PA46, PA66, PPA) are manufactured through injection molding, and the push-out forces are measured through a set of measurements. Finally, a mathematical model considering different cage geometries and materials is developed, and the results are validated by the experiments.

To experimentally investigate the effect of the selected materials and geometrical parameters on the push-out force variation and later validate the theoretical model, a special polymeric cage was designed in this work. Another purpose of this design was to explore how the selected design parameters and their interactions can affect the production of the polymeric cages, i.e., the injection molding process and the required tooling system. Among all the possible design parameters, four geometrical parameters were selected based on their influence on the push-out force and considering mold design limitations. The selected design parameters, their variation levels, and the graphical representation of each parameter are indicated in Fig. 1.

Selected geometrical parameters, their variation levels, and graphical representation.

Lube grooves are semi-circular grooves considered on the outer surface of the cage mainly for allowing additional lubrication into the bearing. The existence of this groove and its depth (A1) can affect the push-out force by allowing the retention tab to deform while the roller is moving outward. In a similar mechanism, another function of this feature is to assist mold inserts being removed from the pockets during the molding process—this action is here called mold insert extraction. The angle of the retention tab (F1) is another important geometrical parameter that directly affects the push-out force. This angle differentiates the retention tab from the rest of the pocket and creates the undercut feature that limits the mold insert extraction. Parts of the retention tab are usually cut off to reduce the undercut area, thus easing the mold insert extraction. The length and location of this cut-off area (C5) also impacts the push-out force in the same manner. The radius of the fillet on the outward edge of the retention tab, which is here called outward pocket rounding (E1), is another important geometrical parameter that influences both the insertion and extraction of the rollers in and out of the pocket. The larger fillets also reduce the severity of the undercut feature and prevent the retention tab edge from being distorted during the molding process. Based on the variation levels considered for these geometrical parameters, a special polymeric cage was designed with 18 different pockets, as listed in Table 1.

The designed polymeric cage and its overall dimensions are shown in Fig. 2 from several points of view. A pocket number was marked on the side of the cage to clearly distinguish each pocket. This numbering was also useful during experimentally measuring push-out forces to make sure proper values are reported for each pocket. As shown in Fig. 2, pocket 18 has no lube groove on one side and a deep lube groove on the other side. This is not a usual case in actual design of polymeric cages and causes unbalanced conditions in this pocket. Also, minor imbalance occurs in pockets 4, 8, and 12 due to the need for change in lube groove designs. Hence, pockets 4, 8, 12, and 18 were later excluded from the push-out force investigation to prevent erroneous feeding data into the modeling process.

The design of the polymeric cage with 18 different pockets.

An injection molding tool was specially designed and fabricated to manufacture the designed cage. This specialized tool and its main parts are demonstrated in Fig. 3. A modular approach was employed in this design of the tool by considering several inserts that can be modified and replaced to perform experimental iterations and develop new polymeric cages (beyond the scope of this work). Moreover, to reduce the complexity in this prototyping tool and cost-efficiently fabricate it, the mold was designed without considering any side actuator to extract the pocket inserts during the molding process. Hence, the injected part is ejected out of the cavity together with the mold inserts inside the pockets (see pocket inserts in Fig. 3). The pocket inserts were later manually removed from the cage, as discussed below.

The specialized injection mold and its main parts (the models were not presented on the same scale).

After the mold was fabricated and tested, three series of polymeric cages were injection molded using the selected materials (PA46, PA66, and PPA). Prior to molding, pellets of the three selected materials were dried in an oven at 80 or 120 °C for 2–4 h to remove any residual moisture. The polymeric cages were then manufactured using an Engel Victory 30 injection molding machine. A minimum of ten cages were manufactured for each of the three selected materials with relevant molding parameters listed in Table 2.

The overall procedure of manufacturing and assembling these cages is shown in Fig. 4, along with the final number of molded parts for each material. In the first step, the mold is in the closed position and ready for injection. After all the sequences of the injection molding are completed, the injected and solidified part (the cage) is ejected out of the cavity while the pocket inserts are still in the pockets. In the next step, the pocket inserts are removed from the cage using a manual tool designed for this purpose. Finally, all the pocket inserts are placed back inside their specific slots by matching their numbers. Now, the mold is ready for the next cycle and following the same procedure. The overall time of producing each cage was about 7 ± 1 mi, which is reasonable for prototyping the polymeric cages for research and development.

The procedure of manufacturing the polymeric cages and the final number of the molded parts for each material.

The test setup is illustrated in Fig. 5. A Shimpo (FGE-10X) with a capacity of 50 N, resolution of 0.01 N, and accuracy of ± 0.2% is used to measure the push-out forces. This gauge can display real-time and peak forces. In this study, the peak value is recorded as a measure of the push-out force. A special component is designed and attached to the gauge hook. This component applies uniform pressure on the roller radially from inside the cage towards out of it without any contact with the cage. The force is applied to the roller gradually until the roller pops out, and the peak force is recorded.

Experimental setup of the push-out force measurement.

The cages in this study are made of three different materials, and the same experiments are repeated for each material. The forces are measured five times for each pocket to find the push-out force variations in a specific cage. Moreover, three different samples for each material are tested to determine the manufacturing variations. The force measurements are performed on 18 pockets, out of which 14 completely symmetric pockets are reported here. Therefore, 210 measurements for each material and a total of 630 force measurements are performed on nine cages with different materials (three cages for each material). The results of the push-out force measurements are shown in Tables 3, 4 and 5.

A predictive model is needed to calculate the push-out force in the roller bearings and characterize the effect of each design parameter on the push-out force. Being analytical, this model explicitly reveals the effect of each design parameter on the push-out force, which will be helpful in design optimization. Furthermore, the critical dimensions and their manufacturing tolerances can be introduced based on the constraints on the push-out forces.

The first step in modeling this complex problem is to simplify and estimate the problem. Since the push-out force is uniformly applied to the roller and there is no change in parameters along the roller axis, the modeling is carried out based on the cage profile in the plane normal to the roller axis. By considering the deformation mechanism of the cages when the roller is being pushed out, it is realized that when the force is applied to the roller, the cage retention is deformed, and the roller comes out (Fig. 6). This deformation is confined to some region near the retention tip, as shown in Fig. 7.

Push-out force and deflection path directions.

Region of significant deformation.

Since the roller and the cage surfaces are very smooth, and there is not much friction between them, the frictional forces are neglected. The part of the material that undergoes significant deformation while the roller is pushed out is considered as a cantilever beam to model the cage deformation and relate it to the force required to push out the roller. The Region of Significant Deformation (ROSD) is considered as the region above the line connecting the bottom of the lube groove to the bottom of the retention, as shown in Fig. 7.

To model the ROSD as a cantilever beam, a reference height and length are needed. The ROSD height is termed as effective height (the length of the line connecting the bottom of the lube groove to the bottom of the retention). The ROSD length is termed as effective length (the length of the line connecting the middle of the effective height line to the tip of the retention), as shown in Fig. 8. The cage geometry is not regular and changes from cage to cage. Therefore, a profile that adapts itself to the cage geometry and provides an accurate estimation is needed. Here, a parabolic profile with the effective height as its height and the effective length as its length is introduced to estimate the ROSD and capture various cage profiles (see Fig. 8).

Modeling the deformed areas of the cages as a parabolic curve.

Therefore, the problem is considered as a cantilever beam with a parabolic profile and a concentrated load at the tip, as depicted in Fig. 9.

Cantilever beam with a parabolic profile.

The constant \(a\) in the equation \({y}^{2}=ax\) shown in Fig. 9 can be found by substituting effective length (\({l}_{eff}\)) for \(x\), and half of the effective height (\({h}_{eff}\)) as \(y\) and solving for the constant \(a\), which gives \(a=\frac{{h}_{eff}^{2}}{4{l}_{eff}}\).

Next, the equation for the deflection of the parabolic cantilever beam is derived based on the strain energy and Castigliano’s theorem. Since the width of the beam is large compared to other dimensions, the effects of shear forces are also considered. The strain energy of a beam due to bending is defined as:

And, the strain energy due to shear for a rectangular cross-section is defined as:

where \(M\) is the moment, \(F\) is the shear force, \(E\) is Young’s modulus, \(G\) is the shear modulus, and \(I\) and \(A\) are the second moment of area and the cross-sectional area, respectively, and are defined as follows:

where \({B}_{5}\) is the width of the beam. The total strain energy of the beam would be:

Performing the integrations gives:

Next, Castigliano’s theorem can be used to calculate the deflection of the beam. Castigliano’s theorem states: “the first partial derivative of the total strain energy in a structure with respect to the force applied at any point is equal to the deflection at the point of application of that force in the direction of its line of action”. Thus, the deflection of the beam can be calculated as follows:

Solving for the force \((F)\) and introducing a correction factor \((C)\) to compensate for the simplifications and estimations gives:

The deflection \(\delta\) here is the difference between the roller diameter and the outer pocket width (\({D}_{1}\)). The notations used here are shown in Fig. 10.

The important geometrical features and dimensions of polymeric cages notations.

The names of the features in Fig. 10 are listed in Table 6.

Furthermore, some of the cage pockets have a discontinuity (\({A}_{55}\) in Fig. 11) in their retention that must be considered in the modeling.

Pocket width and retention discontinuity.

Now, some extra terms are needed to account for the end rim width (\({D}_{5}\)) and retention discontinuity (\({A}_{55}\)). Also, since the fraction \(\frac{{l}_{eff}}{{h}_{eff}}\) is sometimes less than one, and other times very close to one, and the results showed more stability with the 0.5 power for the ratio of the effective length over effective height, the power three is replaced by 0.5 in Eq. (8), and the final push-out force model is cast as follows:

where the effective height and length are shown in Fig. 8, and are defined as:

In Eqs. (9)–(11), \(E\) and \(G\) are material properties, and all other parameters are standard geometrical parameters that are readily available. The correction factor \(C\) is determined by substituting the geometrical parameters and material properties in Eq. (9) and setting the magnitude of the push-out force \(F\) to the mean of the experimental values. There are 14 different pocket geometries and three materials, which means there are 42 different cases. Interestingly the correction factors for all 42 cases were very close to each other that we could set one single number for \(C,\) which works for all materials and geometries. This is a significant achievement proving the ability of the model to consider all the parameters and their effects correctly. The correction factor turned out to be \(C=0.005,\) and the final model can be written as:

In this section, the results of the measurements and calculations based on the proposed model are presented. It is noteworthy that the most critical dimension affecting the push-out force is the pocket width (\({D}_{1}\)), so the tolerances of this dimension must be tight to keep the push-out forces in the allowable region. The shrinkage of the materials has a significant effect on the pocket width, so it should be considered very carefully. The pocket widths of the cages in this study vary slightly from material to material because one mold is used for all the materials, deviating from nominal dimensions. Hence, the actual pocket width is measured in different cages, and the average values are used for each material. The diagram of the measured forces (experimental) and calculated forces (theoretical) for different materials are depicted in Figs. 12, 13 and 14.

Measured and calculated push-out forces in PA46 cages.

Measured and calculated push-out forces in PA66 cages.

Measured and calculated push-out forces in PPA cages.

The average values for material properties and pocket width are considered in all the plots. However, the effects of manufacturing variations on push-out forces must be investigated to define the tolerances of different parameters. Five different coupons for each material are made, and tensile testing is performed to measure the material properties. The standard deviation of the material properties is calculated, and six-sigma is considered as the range of the material properties. As the pocket width is the critical dimension, its tolerance is defined as 0.02 mm, and all other parameters in Eq. (12) are considered to have 0.05 mm tolerance. Based on the tolerances mentioned above, the upper and lower limits of calculated forces along with minimum and maximum measured push-out forces are plotted in Figs. 15, 16 and 17.

The range of measured push-out forces and calculated upper and lower bonds with six sigma variations in material properties, 0.02 mm in pocket width, and 0.05 mm in other dimensions for PA46 cages.

The range of measured push-out forces and calculated upper and lower bonds with six sigma variations in material properties, 0.02 mm in pocket width, and 0.05 mm in other dimensions for PA66 cages.

The range of measured push-out forces and calculated upper and lower bonds with six sigma variations in material properties, 0.02 mm in pocket width, and 0.05 mm in other dimensions for PPA cages.

As can be seen from the plots, all the measured forces are in the range predicted by the model. Therefore, this model can predict the minimum and maximum push-out forces for a specific design with certain material properties and cage dimensions tolerances. This is very helpful in the design stage because the designers must compromise between the minimum push-out force (for the roller not to pop out when delivered to the customer) and the maximum push-out force (the amount of force required to eject the mold inserts without damaging the cage). Moreover, the manufacturing tolerances of the cages can be defined based on the given range of push-out forces.

The fluctuations of measured forces in different samples root in the variations in material properties and, more importantly, in the dimensions of the cages. PPA has the highest consistency in terms of material properties followed by PA66, and PA46 showed the highest deviation in the tensile tests. However, the higher shrinkage and dimensional consistency come together in which PA46 is at the top of the list, followed by PA66 and PPA. This might be because the higher shrinkage makes ejecting the mold inserts easier. The plots reveal that the choice of material does not affect the push-out forces as the dimensional parameters do. For instance, increasing the outward retention angle (F1), which results in smaller pocket width, substantially raises the push-out force. A deeper lube groove reduces the effect of retention angle on the push-out force. Also, a deeper lube groove decreases the stiffness of the cage in the retention tab region, which leads to a decrease in the push-out force. Increasing the retention tab axial gap (C5) reduces the push-out force, while adding a local gap at the middle of the retention tab increases this force.

In this paper, a special polymeric cage is designed and manufactured through injection molding to study the effects of geometrical dimensions and material properties on the amount of push-out force. A specialized injection molding tool was designed and fabricated to realize manufacturing the cages made of three selected materials (PA46, PA66, PPA). The push-out forces were experimentally measured several times on different cages to determine manufacturing and measurement variations. An analytical model for predicting the push-out forces in roller bearings was presented. An empirical coefficient was introduced to the model, which is determined by the experimental push-out force measurements. The correction factor for all the geometries and materials were very close to each other such that one single number could be used as the correction factor in the model. This acknowledges the procedure adopted here to model such problems. Finally, the calculated forces were compared to measured forces, and a very good agreement was achieved. The proposed model is analytical, which shows the effects of each parameter on the amount of push-out force and helps engineers in the design stage where the dimensions of the cages are being determined, and size optimization is performed. Furthermore, the manufacturing tolerances can be defined using this model by introducing the push-out force limits.

The results of this study show that it is better to select the cage material based on the operational conditions of the roller bearing as the amount of push-out force can be adjusted by the dimensional parameters more effectively. The most effective dimension on the push-out force is the pocket width which is controlled by the outward retention angel. Being heavily dependent on the retention angle, the push-out force fluctuates with small variations in this parameter. A deep lube groove is an effective way to reduce the sensitivity of the push-out force to the retention angel and decrease the deviation of the measured forces. Therefore, it is advisable to maximize the lube groove depth and minimize the outward retention angle to achieve the best consistency in the push-out forces.

AST Bearings LLC. Ball Bearing Cages, Retainers, and Ball Separators. Document no: ENB-04-0573 REV. B. (2010).

Gaydamaka, A. Improving the design of cages for bearings of high technical level. Eastern-Eur. J. Enterprise Technol. 84, 4–8 (2016).

Article Google Scholar

Williams, C.W. Development of a new nylon-fiberglass composite cylindrical bearing retainer. SAE Technical Paper 770747. https://doi.org/10.4271/770747 (1977).

Ghaisas, N., Wassgren, C. R. & Sadeghi, F. Cage instabilities in cylindrical roller bearings. J. Tribol. 126, 681–689 (2004).

Article Google Scholar

Kohár, R., Brumerčík, F., Lukáč, M. & Nieoczym, A. Numerical analysis of roller bearing. Appl. Comput. Sci. 12(1), 5–16 (2016).

Google Scholar

Cui, Y., Deng, S., Zhang, W. & Chen, G. The impact of roller dynamic unbalance of high-speed cylindrical roller bearing on the cage nonlinear dynamic characteristics. Mech. Mach. Theory 118, 65–83 (2017).

Article Google Scholar

Rîpanu, M., Nagî, G., Slatineanu, L. & Dodun, O. The dimensional precision of forming windows in bearing cages. Machines. 6(1), 9. https://doi.org/10.3390/machines6010009 (2018).

Article Google Scholar

Cui, Y., Deng, S., Deng, K., Liao, H. & Zhang, W. Experimental study on impact of roller imbalance on cage stability. Chin. J. Aeronaut. https://doi.org/10.1016/j.cja.2020.05.034 (2020).

Article Google Scholar

Zarei, A. & Khosravifard, A. A meshfree method for static and buckling analysis of shear deformable composite laminates considering continuity of interlaminar transverse shearing stresses. Compos. Struct. 209, 206–218 (2019).

Article Google Scholar

Zarei, A. & Khosravifard, A. Meshfree investigation of the vibrational behavior of pre-stressed laminated composite plates based on a variationally consistent plate model. Eng. Anal. Bound. Elem. 111, 118–133 (2020).

Article MathSciNet Google Scholar

Download references

JTEKT North America Corporation/Koyo Bearings North America and Robert Patrick Jenkins Endowed Professorship.

Department of Automotive Engineering, Clemson University, Greenville, SC, USA

Alireza Zarei, Saeed Farahani, Sai Aditya Pradeep & Srikanth Pilla

Clemson Composites Center, Clemson University, Greenville, SC, USA

Alireza Zarei, Saeed Farahani, Sai Aditya Pradeep & Srikanth Pilla

JTEKT North America Corporation, Greenville, SC, USA

John Driscoll & Rob Lukasiewicz

Department of Mechanical Engineering, Clemson University, Clemson, SC, USA

Srikanth Pilla

Department of Materials Science and Engineering, Clemson University, Clemson, SC, USA

Srikanth Pilla

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

A.Z. conducted the modeling study. S.F. and S.A.P. designed and completed the experiments. J.D. and R.L. conducted force measurements. A.Z. and S.F. drafted the manuscript. S.P. edited the manuscript, initiated the concept and directed the research, and advised each step of the study. Figures 1, 2 and 3 were prepared by S.F., Fig. 4 by S.F. and S.A.P., Figs. 5 and 10 by J.D. and R.L., and Figs. 6, 7, 8, 9, 11, 12, 13, 14, 15, 16 and 17 by A.Z. Table 1 was prepared by S.F., Tables 2, 3 and 4 by J.D. and R.L., and Tables 5 and 6 by A.Z.

Correspondence to Srikanth Pilla.

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/.

Reprints and permissions

Zarei, A., Farahani, S., Pradeep, S.A. et al. Design and manufacturing of roller bearing polymeric cages and development of a theoretical model for predicting the roller push-out force. Sci Rep 12, 1017 (2022). https://doi.org/10.1038/s41598-022-04959-9

Download citation

Received: 23 August 2021

Accepted: 03 January 2022

Published: 19 January 2022

DOI: https://doi.org/10.1038/s41598-022-04959-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