Research Article | Open Access
Springback Analysis in Sheet Metal Forming Using Modified Ludwik Stress-Strain Relation
This paper deals with the springback analysis in sheet metal forming using modified Ludwik stress-strain relation. Using the deformation theory of plasticity, formulation of the problem and spring back ratio is derived using modified Ludwik stress strain relationship with Tresca and von Mises yielding criteraia. The results have been representing the effect of different value of or ratio, different values of Strain hardening index (), Poisson’s ratio (), and thickness on spring back ratio (). The main aim of this paper is to study the effects of the thickness, ratio, and Poisson’s ratio in spring back ratio.
In sheet metal working, sheets are deformed to cylindrical and helical shapes by plastic bending with the help of a punch and die set. When such bending is properly done, the inside contour of the section matches the surface contour of the die during the forming operation. However, after the release of the applied loads, the contour assumes a different shape than that of the die because of the release of elastic stresses in the metal. This elastic distortion is commonly called springback. Springback complicates tool design in that the die must be designed to compensate for it. Consequently, it is desirable to have a method of quantitatively predicting the magnitude of springback as a function of the properties of the material and geometry involved in the forming operation.
In bending, the springback is a measure of elastic recovery of radius on removal of the applied bending moment or load after the bending section is beyond elastic limit. While designing the die set, the springback factor should be taken into account to avoid mismatch while assembling formed different sections.
Initially, springback studies are with sheet bending operations. Sachs , Schroeder , Gardiner , Singh and Johnson , and others studied the springback considered bending of sheets of different shapes and depicted springback as a function of material thickness, length, and width of the sheets taken. Their studies were limited to V- and U-shaped dies for applying bending loads and they predicted the springback as a measure of change in the curvature distribution. Huth , Nadai , and Upadhyay  have all done a number of excellent work on the elasto-plastic torsion of bars with rectangular sections, but their work has been limited to monotonically increasing loads only. Dwivedi et al. [8, 9] analytically predicted the residual angle of twist and torque relation, and so forth, for bars of elastic strain-hardening materials with narrow rectangular sections. This work, however, has the limitation that it is valid for thin rectangular strips only. Dwivedi et al. [10, 11] dealt with the torsional springback of square-section bars of linear and nonlinear work-hardening materials. Dwivedi et al.  also dealt with the torsional springback of L-shaped section bars of nonlinear work-hardening materials.
An accurate analysis of springback has been made in the past on sheet bending and tube bending operations through experiment [13–18]. Torsional springback in thin tubes with nonlinear work hardening analysis was made by Choubey et al. [19–21].
In the following, approximation equations are derived in an attempt to provide a quantitative method with practical utility for predicting the springback behavior in bent section of sheet metal as a function of die radius, sheet thickness, and stress-strain characteristics of the material.
2. Basic Equations and Formulation
In this paper, springback prediction approach using the modified Ludwik stress-strain relation was studied. The assumption of narrow beam of a wide sheet can significantly alter the magnitude of the predicted results in sheet metal bending because this ignorance of the transverse stresses are present during forming. In this paper, it is assumed that the bend section is a wide sheet of an elastic plastic strain hardening material in an effort to obtain more accurate expression for predicting the springback behavior. In this method, use of applied moment and curvature relation during the formation of a wide sheet of metal around a portion of cylindrical die was done to derive springback relationship. From the bending equation, we directly relate the bending moment and curvature during the forming as follows:
From this equation, the springback is related to the applied bending moment and curvature during the formation. The modified Ludwik stress-strain relation (Figure 2) is given by where , is the elastic part and is the plastic part, and the equation gives the total stress.
3. Assumption for Derivation of Springback Behaviour
(1)Pure bending of the sheet to cylindrical surface.(2)The cross section dimensions of the sheet are such so that the width to thickness ratio is large. (3)The stress-strain characteristics of material are the same in tension and compression.(4)The natural surface is always in the centre of the sheet, and plane sections remain plane during bending.(5)The cross section dimensions of the sheet do not change significantly in bending.(6)The radius of bending is large compared to the thickness of the sheet so radial stresses are assumed to be negligible. (7)The circumferential strains are sufficiently small so that the conventional strain and the true strain are approximately equivalent. (8)The transverse strain is zero at any point in the sheet. (9)The circumferential strain for any fiber does not vary along the bent section.
4. Derivation of Springback Equation for Tresca Yield Criteria
Figure 1 shows the schematic plot of applied bending moment versus curvature during the formation of a wide sheet of metal around a portion of a cylindrical die at point A, the material yields and plastic deformation continuous until the inside surface of the material conforms to the inside surface of the material conforms to die at point B when the applied moment is released elastic springback occurs from B to C. The change in curvature due to this elastic springback is
From Figure 3,
Figure 3 shows a section of bent sheet and the coordinate system adopted. Then calculation of maximum applied bending moment by balancing the external and internal moments is as follows:
Before proceeding, it is desirable to attempt to characterize the stress-strain behaviour of material in a simple tension. From the modified Ludwik stress-strain relation, we get
This is the equation for springback for Modified Ludwik stress-strain relation.
In elastic region, In plastic region, where , Let Then At yield point, the stress in elastic and plastic region is equal; so we can write From (10), stress in yield point is It will be assumed for stresses less than that the stress-strain behaviour of the material in a simple tension test is elastic (6) and that yielding occurs at and (7) is applicable for stresses greater than .
For combined states of stress, the relationship between principal stresses and strain for the elastic.
By Hooke’s law,
From assumptions that the radius of bending is large compared to thickness of the sheet so radial stresses can assume to be negligible, and the transverse strain is zero at any point in the sheet,
And the circumferential strain for any fiber does not vary along the bent section.
So in yield point, the circumferential strain is
Putting yield point stress value from (11), we get This is the approximate value of circumferential strain at elastic plastic interface. Then in elastic region, the circumferential stress is Here, .
Then, for .
From the deformation theory of plasticity,
From assumptions that the radius of bending is large compared to thickness of the sheet so radial stresses can assume to be negligible and the transverse strain is zero at any point in the sheet, .
And the circumferential strain for any fiber does not vary along the bent section. So, .
Put this value in (23). We get
Putting the value of into (23), we found that
From initial assumption that the circumferential strain are sufficiently small so that the conventional strain and true strain are approximately equivalent,
Putting the value of into (26), we get where is valid for plastic region; that is,
From (4), we found maximum bending moment is where we can split the limit into elastic or plastic part, from 0 to elastic region and to plastic region, where is calculated by Now, Put the value of in above equation.
From elementary plate theory  applied moment per unit width during elastic unloading,
The slope of elastic recovery in elastic unloading is
Putting this into (3), we get
This is the equation for springback ratio for modified Ludwik stress-strain relationship using maximum shear stress theory.
5. Derivation of Springback Equation for Von Mises Yielding
So, in yield point, the circumferential strain is
This is the approximate value of circumferential strain at elastic plastic interface.
Then in elastic region, the circumferential stress is Here, .
Putting value of circumferential strain, we get
From the deformation theory of plasticity ,
From assumptions that the radius of bending is large compared to thickness of the sheet so radial stresses can assume to be negligible and the transverse strain is zero at any point in the sheet, And the circumferential strain for any fiber does not vary along the bent section.
So, Put these values in (44). We get
So, Putting the value of into (44), we get
From initial assumption that the circumferential strain are sufficiently small so that the conventional strain and true strain are approximately equivalent, Putting the value of into (49), we get where is valid for plastic region; that is, From (4), maximum bending moment is We can split the limit into elastic or plastic part from 0 to to 0 to elastic region and to plastic region where is calculated by Putting the value of , we get Then
From elementary plate theory  applied moment per unit width during elastic unloading, Now, the slope of elastic recovery in elastic unloading is Putting this into (3), we get This is the equation for springback ratio for modified Ludwik stress-strain relationship using von Mises yield criteria.
6. Results and Discussion
For the sheet metal bending, the springback is calculated from derived equation (using modified Ludwik stress-strain relationship with Tresca and von Mises yield criteria). Here, from (8),
We can see that both equations have a last term that is very negligible because the or and its three times will be ; therefore, we have to neglect the last term of both equations.
Then we write, Here, we can see that both (62) and (63) are the same if we can apply Tresca or von Mises, getting same result but last term is different. The values of last term are much less; therefore, we neglect this term.
Therefore, we can say that if we apply Tresca or von Mises yield criteria, then we are getting same result. This equation is dependent on the , , strain hardening coefficient , and is Possion’s ratio.
Figure 4 (), Figure 5 (), Figure 6 (), and Figure 7 () show the variation in springback ratio along with change in with different ratio. The higher value of (strain hardening coefficient) material to closer to elasto-ideal plastic material, and thickness is one mm, and value is , and (for different material c10100 copper, 1100al, and 1065steel). From the figure, we come to know that as the with increasing ratio of is decrease the springback ratio and with decreasing springback ratio is increases.
Figures 8, 9, and 10 show the variation in springback ratio along with with the different value of with fixed value of , , and , here we take the values of , and 0.4. From the figure, we come to know that with increasing value of , springback ratio is decreasing.
Figures 11, 12, and 13 show the variation in springback along with with different value of Possion’s ratio , , and with fixed , , , and . From the figure, we come to know that with increasing value of , springback ratio is increasing.
Figures 14, 15, and 16 show the variation in springback along with with different value of , 0.2, 0.3, and 0.4 and where , , , , and is fixed. From the figure, we come to know that the springback ratio rapidly increases in range of 0 to 2 mm thickness, then the increase in springback ratio is stable.
Based on the results presented in chapter of results and discussion, following conclusions about the springback analysis in sheet metal forming using nonlinear constitutive equation, the following applies.(i)The theoretical analysis for the sheet metal under pure bending has been done, and it was found that the prediction of springback is quite successfull.(ii)As expected, elastic recovery is found to be more with decreasing values of work hardening coefficient. At lower values of n, the material will approach to an elasto-ideally plastic behavior. (iii)Springback ratio increases with increasing thickness.(iv)Springback ratio increases with decreasing ratio of yield point stress to Young’s modulus of elasticity (v)Springback ratio is increasing with increasing Possion’s ratio.(vi)Springback ratio for Tresca and von Mises yield criteria, we say that is less than 20, there will be slightly change in that is negligible and for high value of or more than 20, differences in springback ratio is increases.
|:||Radius of die|
|:||Radius of final product after removing load|
|:||Maximum applied bending moment per unit width|
|:||Distance from midsection|
|:||Modulus of elasticity|
|:||Strain hardening index|
|:||Yield point stress|
|:||Strain at yield point|
|:||Strain in plastic region.|
- G. Sachs, Principles and Methods of Sheet Metal Fabricating, Reinhold, New York, NY, USA, 1951.
- W. Schroeder, “Mechanics of sheet metal bending,” Transactions of the ASME, vol. 36, pp. 138–145, 1943.
- F. J. Gardiner, “The springback of metals,” Transactions of the ASEME, vol. 49, pp. 1–9, 1958.
- A. N. Singh and W. Johnson, “Springback after cylindrically bending metal strips,” in Proceedings of the Dr Karunesh Memorial International Conference, New Delhi, India, December 1979.
- J. H. Huth, “A note on plastic torsion,” Journal of Applied Mechanics, vol. 22, pp. 432–434, 1955.
- A. Nadai, Theory of Flow and Fracture of Solids, vol. 1, McGraw-Hill, 1950.
- P. C. Upadhyay, Elasto-plastic torsion [M.Tech. thesis], Mechanical Engineering Department, I.I.T. Kanpur, 1970.
- J. P. Dwivedi, A. N. Singh, S. Ram, and N. K. Das Talukder, “Springback analysis in torsion of rectangular strips,” International Journal of Mechanical Sciences, vol. 28, no. 8, pp. 505–515, 1986.
- J. P. Dwivedi, P. K. Sarkar, R. S. Ram. S., N. K. D. Talukder, and A. N. Singh, “Experimental aspects of torsional springback in rectangular strips,” Journal of the Institution of Engineers, vol. 67, no. 4, pp. 70–73, 1987.
- J. P. Dwivedi, A. K. Shukla, and P. C. Upadhyay, “Torsional springback of square section bars of linear wor hardening materials,” Computers and Structures, vol. 45, no. 3, pp. 421–429, 1972.
- J. P. Dwivedi, P. C. Upadhyay, and N. K. Das Talukder, “Torsional springback in square section bars of nonlinear work-hardening materials,” International Journal of Mechanical Sciences, vol. 32, no. 10, pp. 863–876, 1990.
- J. P. Dwivedi, P. C. Upadhyay, and N. K. D. Talukder, “Springback analysis of torsion of L-sectioned bars of work-hardening materials,” Computers and Structures, vol. 43, no. 5, pp. 815–822, 1992.
- Z. T. Zhang and S. J. Hu, “Stress and residual stress distributions in plane strain bending,” International Journal of Mechanical Sciences, vol. 40, no. 6, pp. 533–543, 1998.
- T. Kuwabara, “Advances in experiments on metal sheets and tubes in support of constitutive modeling and forming simulations,” International Journal of Plasticity, vol. 23, no. 3, pp. 385–419, 2007.
- H. K. Yi, D. W. Kim, C. J. van Tyne, and Y. H. Moon, “Analytical prediction of springback based on residual differential strain during sheet metal bending,” Journal of Mechanical Engineering Science, vol. 222, no. 2, pp. 117–129, 2008.
- A. El Megharbel, G. A. El Nasser, and A. El Domiaty, “Bending of tube and section made of strain-hardening materials,” Journal of Materials Processing Technology, vol. 203, no. 1–3, pp. 372–380, 2008.
- E. Da-xin, H. Hau-hui, L. Xiao-yi, and N. Ru-xin, “Spring-back deformation in tube bending,” International Journal of Minerals, Metallurgy and Materials, vol. 16, no. 2, pp. 177–183, 2009.
- D. E. Daxin E and Y. Liu, “Springback and time-dependent springback of 1Cr18Ni9Ti stainless steel tubes under bending,” Materials and Design, vol. 31, no. 3, pp. 1256–1261, 2010.
- V. K. Choubey, M. Gangwar, and J. P. Dwivedi, “Torsional springback analysis in thin tubes with non-linear work hardening,” Journal of Mechanical Engineering, vol. 7, no. 1, pp. 15–34, 2011.
- V. K. Choubey, M. Gangwar, and J. P. Dwivedi, “Springback analysis of thin tubes,” Journal of Mechanical Engineering, vol. 7, no. 2, 2011.
- V. K. Choubey, M. Gangwar, J. P. Dwivedi, and N. K. das Talukder, “Springback analysis of thin tubes with arbitrary stress-strain curves,” Journal of Mechanical Engineering, vol. 8, no. 1, 2011.
Copyright © 2013 Sanjay Kumar Patel et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.