Next Article in Journal
Comparative Analysis of Deep Learning Models for Optimal EEG-Based Real-Time Servo Motor Control
Next Article in Special Issue
Assessing CNN Architectures for Estimating Correct Posture in Cruise Machinists
Previous Article in Journal
On the Importance of Solar Radiation and Shading for the Alkali–Aggregate Reaction Prediction of Concrete Arch Dams
Previous Article in Special Issue
Development, Designing and Testing of a New Test Rig for Studying Innovative Polycrystalline Diamond Bearings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Prediction of Spring-Back in Cylindrical Helical Springs Using Analytical and Numerical Models

Department of Mechanical and Industrial Engineering, University of Brescia, Via Branze, 38, 25123 Brescia, Italy
*
Author to whom correspondence should be addressed.
Eng 2024, 5(3), 1696-1707; https://doi.org/10.3390/eng5030089
Submission received: 5 July 2024 / Revised: 29 July 2024 / Accepted: 30 July 2024 / Published: 2 August 2024
(This article belongs to the Special Issue Feature Papers in Eng 2024)

Abstract

:
This research focuses on cylindrical helical springs with circular cross-sections made from carbon steel (SH 0.82% C) and stainless steel (AISI 302). The transformation from a linear bar to a circular spiral involves numerous factors such as material mechanical behavior, stress–strain relationships and residual stresses. This research investigates the spring-back phenomenon, which affects the final diameter of helical springs post-manufacture, using analytical, experimental and numerical methods. An analytical model, derived from the mechanical bending process, was proposed to predict spring-back, and its accuracy was validated against experimental data. This study also employed finite element simulations to analyze elastic recovery, confirming the analytical predictions. Results indicated that the spring-back ratio k could be expressed as an exponential function of the spring index C (the ratio between the final diameter of the spring D2 and the diameter of the wire DW), with a maximum error of 4.80% for stainless steel and 3.62% for carbon steel. This study’s findings provide valuable insights into optimizing the spring manufacturing process, enhancing the precision of spring diameter predictions, and potentially reducing production errors and material waste.

1. Introduction

Springs are one of the most important components used in various industrial applications. Broadly speaking, it can be said that there are no machines that do not utilize this component. The primary applications of springs include addressing vibration-related problems to reduce dynamic effects, preventing shock loads, applying force in a specific pattern, and more.
For the purposes mentioned above, numerous types of springs are used in the automotive sector. For example, Salah [1] discussed the design of coil springs for automotive suspension systems, while Solazzi [2] presented the design of a torsion spring bar for the cab suspension system of an industrial vehicle.
One of the most important parameters in the design of springs is the evaluation of their dynamic behavior, which is related to both the natural frequencies of vibration and the displacement of the spring under a specific load. This aspect is crucial because it can either increase or decrease the displacement or transmission force between the two elements connected by the spring [3,4,5].
The static axial load also influences the magnitude of natural frequencies [6]. The design problem of coil springs, particularly the estimation of deflections or the stress state of the wire, is highly complex due to its nonlinear nature [7]. This can be evaluated using analytical or numerical methods [8,9,10]. Recently, to reduce component weight, innovative materials have been introduced for the construction of coil springs [11,12,13].
Both the design and manufacturing processes of coil springs are very complicated. This complexity arises because the process begins with a linear straight bar that is transformed into a circular spiral, forming the coil spring. This process involves many factors, such as the mechanical behavior of the material, including the stress–strain curve (particularly Young’s modulus and yield stress [14], plasticity, and residual stress [15,16]. Additionally, whether the material exhibits kinematic or isotropic hardening [14] significantly affects the mechanical behavior.
Helical compression springs for automotive suspensions are made using either cold coiling or hot coiling processes. In hot coiling, steel wire is heated above the austenitizing temperature, coiled, quenched, and tempered, which can cause decarburization and surface scaling, detrimental to fatigue performance. In cold coiling, wire is coiled at room temperature, preset with torsion residual stress, and shot peened, avoiding decarburization and scaling since stress relief occurs below the austenitizing temperature. Optimizing stress management in springs involves balancing stress relief temperature, wire strength, shot peening intensity, and presetting conditions. Understanding the effects of each processing step on residual coiling stresses is essential for achieving optimal residual stress in helical springs. Kobelev [16] presented a model to simulate residual stresses in helical springs due to pre-setting, using the deformational theory of plasticity. It provided analytical formulas to assess residual stresses, changes in spring shape, and pre-setting force over spring travel. The method allowed expressing torque as a function of spring travel and residual stresses in a closed form. The analysis included elastic–plastic bending and torsion, offering closed-form solutions for preset moments and elastic spring-back.
There is little literature on experimental measurements of spring-back in helical springs and its analytical prediction based on experimental values. In this paper, we focused on cylindrical helical springs with circular cross-sections made of two different steels. The bending theory of plates was applied to verify its applicability in spring production. In plate production, elastic recovery can simply be measured as the variation in the angle between the bent part before and after recovery, in relation to the angle produced by the tool, or by the ratio of these angles. Numerous diagrams allow manufacturers to account for recovery in sheet production, typically reporting the recovery with respect to the ratio between the angle after recovery and plate thickness. Marciniak [17] and Gardiner [18] first proposed an analytical solution to predict elastic recovery, but they only hypothesized a perfectly plastic elastic behavior. Gonzalez-Coleo et al. [15] later introduced a model that considered hardening after yielding. The main drawback of these studies is that they are based on the elastic curve formulation, which is consistent with small thickness geometries; indeed, in these studies, the error in the estimation of recovery increased with increasing sheet thickness. Furthermore, larger thicknesses are affected by different hardening areas that are not considered in the models.
In this article, the bending theory of plates is used to predict the elastic recovery of cylindrical springs with circular cross-sections. First, the machine used to produce these springs will be described, and data will be collected for different sizes and two materials. Then, a model for estimating recovery will be proposed, and the equation proposed in [15] will be validated. Finally, as an alternative to analytical models, a numerical finite element model with a nonlinear solver will be proposed.

2. Experimental Details

Modern machines used to produce cylindrical helical springs are numerically controlled with 4 to 12 axes. These machines consist of a straightening and feeding area for the wire and a winding-cutting area for the spring. The wire, unwound from a coil, passes through a series of straightening rollers and then through two or four pairs of rollers that press the wire from above and below, rotating and pushing it to the next stage. Figure 1 shows a machine that produces springs from wire with diameters ranging from 2.5 to 8 mm.
Figure 2 shows a detailed view of the feeding and winding section with the following components: (1) feed rollers, which pull the wire from the coil (fed from the left in this photo) and push it to the right; (2) vertical pitch, which moves up and down to deform the wire and create the spring’s pitch; (3) cutter, which cuts the wire against the central mandrel; (4) horizontal pitch, which has the same function as the vertical pitch but is used for springs with larger diameters; (5) coiling fingers, which guide and bend the wire to form the helix; and (6) central mandrel. Smaller machines have six axes, while larger ones have ten axes, with four motors controlling the diameter (radial and tangential movement for the two winding tips) and an additional one for the rotation of the lower tip.
To regulate the coiling diameter on a machine and control the diameter axes via the system interface, manual adjustments are required for the axial position and vertical plane rotation at the first and second winding points. These tasks are crucial during winding to center the spring relative to the axis of the central pin, ensuring that the wire is cut at the highest horizontal point of the coil and that the contact point of the wire with the tips is centered within the channel. A crucial aspect of spring production is the diameter variation due to spring-back phenomena.
Given the manual nature of these adjustments, achieving accurate and repeatable positioning of the winding tips to form an exact circle is challenging. Indeed, with identical diameter axis positions displayed on the control, incorrect adjustments can result in springs with diameters differing by up to several millimeters. To quantify this effect, for each wire size, the position of the diameter axis relative to the winding diameter was documented alongside the winding and Δ D graphs (where Δ D is the diameter variation due to spring-back phenomena).
The relationship between the electronically set diameter interval and the measured diameter difference exhibits a consistent and repeatable trend.
Accurate measurements of the cut spring can be obtained using a centesimal caliper. However, measuring the winding diameter is more challenging; only half a coil can be measured on the machine, and the space for the caliper is very limited for small spring diameters. Additionally, small wire diameters are very flexible, and even slightly excessive pressure from the caliper can result in inaccurate measurements. The maximum estimated error for these measurements was approximately 0.2 mm, an acceptable margin for this work.
Springs with wire diameters ranging from 0.05 to 12 mm for this study were selected. The investigated springs were made of carbon steel type SH (0.82% C) and stainless steel AISI 302, which are among the most common materials for producing small and medium-sized springs. The chemical compositions of these steels are reported in Table 1 and Table 2. Table 3 reports the ultimate tensile stress and the yield stress.
Various springs were manufactured for each wire size and material type, with a spring index ranging from 4 to 25. The spring index C is defined as the ratio between the final diameter of the spring D 2 and the diameter of the wire DW.
For each spring size, we measured the diameter across the two winding points to calculate the diameter increase. Figure 3 illustrates the diameter variation ( Δ D = D 2 D 1 ) of the spring before ( D 1 ) and after spring-back ( D 2 ) for carbon steel springs (a) and stainless-steel springs (b) according to the spring index C . The higher elastic return in stainless steel can be ascribed to the lower Young modulus (185 GPa for the stainless steel and 206 GPa for the carbon steel) combined with the lower yield strength.

3. Analytical Approaches

This section reviews various approaches for predicting the spring-back.
The first method we investigate is derived from the mechanical bending process and based on the spring-back ratio k , defined as follows [22]:
k = D 1 + D W D 2 + D W = α 2 α 1
where α 1 and α 2 are for the bending angles (see Figure 4). This approach has been widely adopted for several materials, but curves for the steels investigated in this paper are not available.
Figure 5 shows the spring-back parameter k as a function of the spring index C for the two materials. An exponential function was chosen to fit the points:
k = a 1   e a 2 C
where Table 4 reports the estimated variables for carbon steel and stainless steel, with each coefficient of determination R2.
Table 5 and Table 6 report the spring geometrical parameters, the k parameter derived by applying Equation (1) and using the diameters D 1 and D 2 estimated during the spring production process ( k e x p ), and the k parameter calculated from Equation (2). The maximum error is about 3.62% and 4.80% for carbon and stainless steels, respectively.
Gonzalez-Coleo et al. [15] recovery model was based on the plasticity constitutive law of the steel:
1 r 1 r 2 = 8 1 n n + 2 r 1 σ y D W E 3 + 3 2 1 n n + 2 r 1 σ y D W E 1 n
where n is the nondimensional hardening exponent describing the plastic trend of the tensile test, E the Young’s modulus, σ y the yield stress, r 1 the radius curvature before the spring-back, r 2 the radius curvature after the spring-back. The yield stresses are reported in Table 4, while the strain hardening exponent is 0.1. Table 7 shows the relative error between the nominal diameter experimentally found and that calculated from Equation (3); the error increases with the spring index C . The error can be ascribed to the fact that Equation (3) focuses solely on the bending phenomenon. As the spring diameter increases, it is overly simplistic to consider only bending; torque should be included as well, for instance.

4. Numerical Simulations

We conducted finite element simulations to investigate the effect of elastic return in spring production using Ansys. This analysis requires a nonlinear solver. The model contains six components, as labeled in Figure 6: spring steel wire (1), coiling mandrel (2), and two rollers that bend the wire on the inner disc (3), which are joined to the support braces (4). Except for the wire, all components are designated as rigid, allowing the mesh to be generated solely on the surfaces in contact with the wire, preventing deformations of the rollers or arms that could result in inaccurate calculations and extended solving times.
SH class steel was characterized using the tensile test chart of the same type and batch of material used in the coiling tests previously reported. Table 8 lists the mechanical properties of the steel. To account for the plastic behavior, we implemented an elastic–plastic constitutive model; in detail, a bilinear law approximated the stress–strain tensile curve.
The two support braces are attached to the central disc using revolute joints, allowing them to rotate. Similarly, the two rollers are attached to the upper pivots of the arms using revolute joints, enabling their rotation. The central disc is fixed to the ground. The wire is constrained from moving vertically, keeping it centered on the disc and rollers throughout the process. The bending operation is performed by rotating the arms: the right arm rotates clockwise up to 90° and then back to 75°, while the left arm rotates counterclockwise up to 90° and then back to 75°. This motion retracts the rollers, releasing the wire.
The interaction between the wire and the rollers is set through a contact constraint with an Augmented Lagrange method. This was chosen because it provides more accurate results than those obtained through a Pure Penalty method and shorter calculation times, much less than what is required by Normal Lagrange. Additionally, the chosen method does not present any risk of contact vibration (which is possible with Normal Lagrange) and has low permissible contact penetration values. The pure penalty method allows high overlapping values, which result in inaccurate diameter measurements.
The contact behavior is set to asymmetric, and the friction coefficient is set at 0.4. The coiling tools of the actual machine are made of lapped tungsten carbide, while the wire is covered with a thin layer of stearate, which reduces friction during the forming process.
The components are discretized using tetrahedral elements. For the simulation with 6 mm diameter wire and 65 mm coiling diameter, the number of nodes is 26,894, and the number of elements is 5871, as shown in Figure 7. These belong mostly to the steel wire since all other elements are rigid, and the mesh is generated only on the surfaces of the components that come into contact with the wire. Figure 8 and Figure 9 illustrate the progression of the wire throughout the simulations.
The simulations results were calculated based on the parameter k . Firstly, the angles α 1 and α 2 were estimated from the simulations as shown in Figure 9. Subsequently, k F E M using Equation (1) and, finally, the final diameter D 2 ,   F E M according to the following equations:
D 2 , F E M = 1 k F E M D 1 + s s  
where D 2 , F E M is the diameter after the spring-back estimated by the finite element mode, D 1 is the diameter before the spring-back and D w the wire diameter.
Table 9 shows the simulation results. The extent of spring-back increases nonlinearly as the spring index increases. Of note, the relative errors found in the numerical models are consistent with those calculated with Equation (3). The relative errors are not more than 8%. However, the computation time (in general more than 5 h) makes simulations a less practical method compared to an analytical prediction.

5. Conclusions

This research provided an in-depth analysis of the spring-back phenomenon in cylindrical helical springs made of carbon steel and stainless steel. By employing both experimental measurements and analytical models, the exponential relationship between the spring-back ratio k (Equation (1)) and the spring index C (the ratio between the final diameter of the spring D 2 and the diameter of the wire 2s) is determined. The maximum errors for the predicted spring-back were 4.80% for stainless steel and 3.62% for carbon steel, demonstrating the model’s accuracy. Furthermore, a model requiring the initial diameter, thickness, and plastic properties of the material as input data were employed. The outcomes of this model were dependent on the spring index C , with the error increasing as C increased.
Additionally, finite element simulations corroborated the analytical predictions, though they required significant computation time. These findings offer valuable insights for optimizing spring manufacturing processes, enhancing precision, and minimizing material waste.
Future work could explore more geometries and materials to further refine predictive models and improve manufacturing efficiency.

Author Contributions

Conceptualization, L.S.; methodology, N.Z. and L.S.; validation, N.Z. and L.S.; investigation, N.Z. and L.S.; resources, L.S.; data curation, N.Z.; writing—original draft preparation, N.Z.; writing—review and editing, N.Z. and L.S.; visualization, N.Z.; supervision, L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

The authors would like to thank you very much to Marco Biggiogero for the experimental support for this research.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Salah, H.R.A.; Mostafa, A.; Eslam, E.; Youssef, T.; Ibrahim, S.; Mina, S. 3D Design of helical spring for automotive independent suspension system under fatigue test conditions. In Applied Mathematics, Modelling and Computer Simulation; IOS Press: Amsterdam, The Netherlands, 2023. [Google Scholar] [CrossRef]
  2. Solazzi, L. Design of a truck cab torsion bar. Int. J. Heavy Veh. Syst. 2019, 26, 692–706. [Google Scholar] [CrossRef]
  3. Yildirim, Y. Expressions for predicting fundamental natural frequencies of non-cylindrical helical springs. J. Sound Vib. 2002, 252, 479–491. [Google Scholar] [CrossRef]
  4. Yildirim, V. Free vibration of uniaxial composite cylindrical helical springs with circular section. J. Sound Vib. 2001, 239, 321–333. [Google Scholar] [CrossRef]
  5. Yu, A.M.; Hao, Y. Free vibration analysis of cylindrical helical springs with noncircular cross-sections. J. Sound Vib. 2011, 330, 2628–2639. [Google Scholar] [CrossRef]
  6. Kobelev, V. Effect of static axial compression on the natural frequencies of helical springs. Multidiscip. Model. Mater. Struct. 2014, 10, 379–398. [Google Scholar] [CrossRef]
  7. Zhuo, Y.; Qi, Z.; Zhang, J.; Wang, G. A geometrically nonlinear spring element for structural analysis of helical springs. Arch. Appl. Mech. 2022, 92, 1789–1821. [Google Scholar] [CrossRef]
  8. Yu, A.; Yang, C. Formulation and evaluation of an analytical study for cylindrical helical springs. Acta Mech. Solida Sin. 2010, 23, 85–94. [Google Scholar] [CrossRef]
  9. Čakmak, D.; Tomičević, Z.; Senjanović, I.; Wolf, H.; Božić, Ž.; Semenski, D. A Study on Cylindrical Coil Spring Deflection and Stress Done Using Analytical and Numerical Methods. Int. J. Eng. Model. 2022, 35, 57–81. [Google Scholar] [CrossRef]
  10. Kamil, F.; Abdulshaheed, A.G.; Kadhom, M.A. Development of a computational system to design a helical spring. Int. Rev. Appl. Sci. Eng. 2014. [Google Scholar] [CrossRef]
  11. Chen, L.; Wu, L.; Fu, H.; Tang, Y. Design and Performance Evaluation of Polymer Matrix Composite Helical Springs. Polymers 2022, 14, 3900. [Google Scholar] [CrossRef] [PubMed]
  12. Wu, L.; Chen, L.; Fu, H.; Jiang, Q.; Wu, X.; Tang, Y. Carbon fiber composite multistrand helical springs with adjustable spring constant: Design and mechanism studies. J. Mater. Res. Technol. 2020, 9, 5067–5076. [Google Scholar] [CrossRef]
  13. Nirala, A.; Kumar, N.; Bandhu Singh, D.; Kumar Singh, A.; Kumar Sharma, S.; Kumar Yadav, J.; Prasad, H.; Chandan, S.; Kumar Shrivastava, A. Simulation analysis of composite helical spring for compression, torsional and transverse mode. Mater. Today Proc. 2020, 28, 2263–2267. [Google Scholar] [CrossRef]
  14. Kobelev, V. Elastic–plastic deformation and residual stresses in helical springs. Multidiscip. Model. Mater. Struct. 2020, 16, 448–475. [Google Scholar] [CrossRef]
  15. González-Coneo, J.; Gonzalez-Olier, C.; Pedraza-Yepes, C.; Higuera-Cobos, O.; Troncoso-Palacio, A. Analytical method for elastic recovery prediction of air bending sheet. J. Appl. Eng. Sci. 2023, 21, 973–981. [Google Scholar] [CrossRef]
  16. Kobelev, V. Elastic-plastic work-hardening deformation under combined bending and torsion and residual stresses in helical springs. Int. J. Mater. Form. 2010, 3 (Suppl. S1), 869–888. [Google Scholar] [CrossRef]
  17. Marciniak, Z.; Duncan, J.L.; Hu, S.J. Mechanical of Sheet Metal Forming; Elsevier: Oxford, UK, 2002. [Google Scholar] [CrossRef]
  18. Gardiner, F.J. The springback of metals. Trans. ASME 1957, 79, 1–7. [Google Scholar] [CrossRef]
  19. EN 10270-1; Steel wire for mechanical springs. Part 1: Patented cold drawn unalloyed spring steel wire. iTeh Standards: Newark, DE, USA, 2017.
  20. EN 10270/3; Steel wire for mechanical springs. Part 3: Stainless spring steel wire. iTeh Standards: Newark, DE, USA, 2011.
  21. DIN EN ISO 6931-1; Stainless steels for springs—Part 1: Wire. ISO: Geneva, Switzerland, 2016.
  22. Kalpakjian, S.; Schmid, S. Manufacturing Engineering & Technology, 7th ed.; Pearson Publishing Company: Upper Saddle River, NJ, USA, 2014. [Google Scholar]
Figure 1. Bench machine for 2.5–8 mm springs.
Figure 1. Bench machine for 2.5–8 mm springs.
Eng 05 00089 g001
Figure 2. Detail of machine coiling systems: (1) feed rollers, (2) vertical pitch, (3) cutter, (4) horizontal pitch, (5) coiling fingers, and (6) central mandrel.
Figure 2. Detail of machine coiling systems: (1) feed rollers, (2) vertical pitch, (3) cutter, (4) horizontal pitch, (5) coiling fingers, and (6) central mandrel.
Eng 05 00089 g002
Figure 3. Experimental correlation between index factor C and spring-back.
Figure 3. Experimental correlation between index factor C and spring-back.
Eng 05 00089 g003
Figure 4. Definition of angles and radii r1 and r2 before recovery (1) and after recovery (2).
Figure 4. Definition of angles and radii r1 and r2 before recovery (1) and after recovery (2).
Eng 05 00089 g004
Figure 5. Spring-back parameter k of carbon steel and stainless steel as a function of spring index C (dotted lines represent the interpolation curves).
Figure 5. Spring-back parameter k of carbon steel and stainless steel as a function of spring index C (dotted lines represent the interpolation curves).
Eng 05 00089 g005
Figure 6. Finite element model: (1) spring steel wire, (2) coiling mandrel, (3) rollers, and (4) support braces.
Figure 6. Finite element model: (1) spring steel wire, (2) coiling mandrel, (3) rollers, and (4) support braces.
Eng 05 00089 g006
Figure 7. Detail of model mesh: (a) overall model; (b) wire detail.
Figure 7. Detail of model mesh: (a) overall model; (b) wire detail.
Eng 05 00089 g007
Figure 8. Deformation steps during the simulation.
Figure 8. Deformation steps during the simulation.
Eng 05 00089 g008
Figure 9. Elastic return detail.
Figure 9. Elastic return detail.
Eng 05 00089 g009
Table 1. Chemical composition of the SH steel (EN 10270-1) [19].
Table 1. Chemical composition of the SH steel (EN 10270-1) [19].
CSiMnPSCu
0.35–1.000.10–0.300.50–1.200.035 max0.035 max0.20 max
Table 2. Chemical composition of the AISI 302 steel (EN 10270/3 [20]; DIN EN ISO 6931-1 [21]).
Table 2. Chemical composition of the AISI 302 steel (EN 10270/3 [20]; DIN EN ISO 6931-1 [21]).
CSiMnPSNCrMoNi
0.05–0.15≤2.00≤2.00≤0.045≤0.015≤0.1016.0–19.0≤0.806.0–9.5
Table 3. Mechanical properties according to the wire diameters.
Table 3. Mechanical properties according to the wire diameters.
Carbon SteelStainless Steel
D W   [ m m ] σ R   [ M P a ] σ y   [ M P a ] D W   [ m m ] σ R   [ M P a ] σ y   [ M P a ]
0.25269018830.2523101617
0.6242316960.620601442
0.9210016061.519001330
1.5209814692.517841249
2.518861320
616481154
Table 4. Coefficient of interpolation for Equation (2).
Table 4. Coefficient of interpolation for Equation (2).
Material a 1 a 2 R 2
Carbon steel1.0309−0.0090.96
Stainless steel1.0257−0.010.98
Table 5. Spring geometry details and results from spring-back recovery k for carbon steel ( k e x p is the parameter experimentally found; k i n t is the one calculated from Equation (2)).
Table 5. Spring geometry details and results from spring-back recovery k for carbon steel ( k e x p is the parameter experimentally found; k i n t is the one calculated from Equation (2)).
DW [mm]D2 [mm]D1 [mm]kexpkintError kDW [mm]D2 [mm]D1 [mm]kexpkintError k
0.2543.400.8590.8863.13%2.51211.700.9790.9850.60%
0.2554.070.8230.8533.62%2.51514.530.9730.9740.09%
0.2564.720.7950.8213.23%2.51615.530.9750.9700.44%
0.2586.000.7580.7610.43%2.51817.280.9650.9630.20%
0.25107.070.7140.7051.26%2.52019.000.9560.9560.01%
0.654.640.9360.9531.81%2.52220.700.9470.9480.15%
0.665.510.9260.9381.29%2.52422.360.9380.9410.33%
0.687.120.8980.9081.20%2.52825.700.9250.9270.26%
0.6108.600.8680.8801.40%2.53228.980.9120.9130.06%
0.61210.000.8410.8531.36%2.53632.050.8970.8990.21%
0.9109.100.9170.9281.13%2.54035.100.8850.8860.11%
0.91210.660.8960.9081.37%2.54236.600.8790.8790.04%
0.91513.000.8740.8800.67%2.54438.150.8740.8720.21%
0.91613.750.8670.8710.46%2.54841.000.8610.8590.26%
0.91815.200.8520.8530.10%2.55042.350.8540.8530.19%
1.5109.640.9690.9680.09%2.55243.940.8520.8460.69%
1.51211.360.9530.9560.32%64240.930.9780.9651.32%
1.51514.030.9410.9380.38%64442.660.9730.9621.18%
1.51614.850.9340.9320.27%64846.110.9650.9560.97%
1.51816.500.9230.9200.33%65047.830.9610.9530.90%
1.52018.200.9160.9080.86%65249.540.9580.9500.83%
1.52219.780.9060.8970.94%66056.320.9440.9380.70%
1.52421.380.8970.8861.29%67064.770.9310.9230.89%
1.52824.420.8790.8641.72%68073.420.9230.9081.63%
1.53227.320.8600.8422.13%610089.000.8960.8801.80%
1.53630.060.8420.8212.46%
Table 6. Spring geometry details and results from spring-back recovery k for stainless steel ( k e x p is the parameter experimentally found; k i n t is the one calculated from Equation (2)).
Table 6. Spring geometry details and results from spring-back recovery k for stainless steel ( k e x p is the parameter experimentally found; k i n t is the one calculated from Equation (2)).
DW [mm]D2 [mm]D1 [mm]kexpkintError kDW [mm]D2 [mm]D1 [mm]kexpkintError k
0.2543.410.8620.8741.91%2.51211.770.9840.9781.13%
0.2554.100.8290.8401.89%2.51514.500.9710.9660.84%
0.2564.760.8010.8071.23%2.51615.390.9670.9620.77%
0.2585.960.7530.7451.26%2.51817.160.9590.9540.63%
0.25107.030.7100.6884.80%2.52018.900.9510.9470.51%
0.654.650.9380.9440.59%2.52220.610.9430.9390.41%
0.665.460.9170.9281.29%2.52422.290.9360.9320.32%
0.687.000.8840.8981.91%2.52825.580.9210.9170.15%
0.6108.480.8560.8681.94%2.53228.750.9060.9020.01%
0.6129.870.8310.8401.64%2.53631.810.8910.8880.11%
1.5109.890.9900.9603.32%2.54034.760.8770.8740.23%
1.51211.520.9640.9471.87%2.54236.190.8690.8670.28%
1.51513.910.9340.9280.49%2.54437.590.8620.8600.34%
1.51614.690.9250.9220.18%2.54840.310.8480.8470.44%
1.51816.240.9100.9100.28%2.55041.630.8410.8400.50%
1.52017.770.8960.8980.57%2.55242.920.8330.8330.55%
1.52219.260.8830.8860.72%
1.52420.730.8720.8740.78%
1.52823.590.8510.8510.65%
1.53226.350.8310.8290.28%
1.53629.000.8130.8070.27%
Table 7. Spring recovery according to Equation (3) for carbon steel and relative error with respect to experimental results.
Table 7. Spring recovery according to Equation (3) for carbon steel and relative error with respect to experimental results.
DW [mm] σ y [MPa]C [-]D2 exp. [mm]D2 calc. [mm]Error D2DW [mm] σ y [MPa]C [-]D2 exp. [mm]D2 calc. [mm]Error D2
0.2518838.044.39%2.514692.41212.54%
0.25188310.055.59%2.514693.01515.85%
0.25188312.066.711%2.514693.21617.06%
0.25188316.089.417%2.514693.61819.16%
0.25188320.01012.121%2.514694.02021.26%
0.616964.255.35%2.514694.42223.36%
0.616965.066.47%2.514694.82425.46%
0.616966.788.68%2.514695.62829.76%
0.616968.31010.99%2.514696.43234.17%
0.6169610.01213.19%2.514697.23638.47%
0.916065.61010.66%2.514698.04042.87%
0.916066.71212.87%2.514698.44245.07%
0.916068.31516.28%2.514698.84447.38%
0.916068.91617.49%2.514699.64851.78%
0.9160610.01819.710%2.5146910.05053.88%
1.513203.31010.55%2.5146910.45256.38%
1.513204.01212.54%611543.54244.35%
1.513205.01515.85%611543.74446.35%
1.513205.31616.95%611544.04850.35%
1.513206.01819.06%611544.25052.45%
1.513206.72021.26%611544.35254.45%
1.513207.32223.46%611545.06062.64%
1.513208.02425.67%611545.87073.14%
1.513209.32830.07%611546.78084.25%
1.5132010.73234.48%611548.3100105.05%
1.5132012.03638.88%6115410.0120127.46%
Table 8. Mechanical properties implemented in finite element simulations.
Table 8. Mechanical properties implemented in finite element simulations.
PropertiesValue
Density [kg/m3]7850
Young’s modulus [GPa]197
Poisson’s ratio [-]0.3
Yield strength [MPa]1154
Tangent modulus [MPa]117
Table 9. Detail of finite elements model results.
Table 9. Detail of finite elements model results.
D 1
[mm]
D 2 , e x p
[mm]
α 1
[°]
α 2
[°]
k F E M
[-]
D 2 , F E M
[mm]
%   Error   D 2
5052.690.6289.360.98650.793%
6570.391.5389.70.98066.455%
7582.691.6900.98376.447%
8088.891.38900.98581.328%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zani, N.; Solazzi, L. Analysis and Prediction of Spring-Back in Cylindrical Helical Springs Using Analytical and Numerical Models. Eng 2024, 5, 1696-1707. https://doi.org/10.3390/eng5030089

AMA Style

Zani N, Solazzi L. Analysis and Prediction of Spring-Back in Cylindrical Helical Springs Using Analytical and Numerical Models. Eng. 2024; 5(3):1696-1707. https://doi.org/10.3390/eng5030089

Chicago/Turabian Style

Zani, Nicola, and Luigi Solazzi. 2024. "Analysis and Prediction of Spring-Back in Cylindrical Helical Springs Using Analytical and Numerical Models" Eng 5, no. 3: 1696-1707. https://doi.org/10.3390/eng5030089

APA Style

Zani, N., & Solazzi, L. (2024). Analysis and Prediction of Spring-Back in Cylindrical Helical Springs Using Analytical and Numerical Models. Eng, 5(3), 1696-1707. https://doi.org/10.3390/eng5030089

Article Metrics

Back to TopTop