Research Article | Open Access

# Springback Analysis in Sheet Metal Forming Using Modified Ludwik Stress-Strain Relation

**Academic Editor:**C. C. Huang

#### Abstract

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.

#### 1. Introduction

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 [1], Schroeder [2], Gardiner [3], Singh and Johnson [4], 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 [5], Nadai [6], and Upadhyay [7] 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. [12] 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,

Then, from putting the value of (12) in (13), From maximum shear stress theory of failure (Tresca yield criteria),

Then,

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

So,

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.

Then

From elementary plate theory [17] 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

From maximum shear stress theory of failure (von Mises yielding criteria), by using (13) and (14).

Substituting (15) in (11), at yield point,

Now,

So, in yield point, the circumferential strain is

Putting yield point stress value from (38) in (40), we get

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

for .

From the deformation theory of plasticity [16],

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 [17] 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),

Then, .

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.

#### 7. Conclusions

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.

#### Nomenclature

: | Radius of die |

: | Radius of final product after removing load |

: | Maximum applied bending moment per unit width |

: | Distance from midsection |

: | Thickness |

: | True stress |

: | True strain |

: | Modulus of elasticity |

: | Poisson’s ratio |

: | Strain hardening index |

: | Empirical constant |

: | Constant |

: | Yield point stress |

: | Strain at yield point |

: | Strain in plastic region. |

#### References

- 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. View at: Google Scholar - F. J. Gardiner, “The springback of metals,”
*Transactions of the ASEME*, vol. 49, pp. 1–9, 1958. View at: Google Scholar - 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. View at: Google Scholar - J. H. Huth, “A note on plastic torsion,”
*Journal of Applied Mechanics*, vol. 22, pp. 432–434, 1955. View at: Google Scholar - 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. View at: Google Scholar - 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. View at: Google Scholar - 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. View at: Google Scholar - 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. View at: Google Scholar - 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. View at: Google Scholar - 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. View at: Google Scholar - 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. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - 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. View at: Google Scholar - V. K. Choubey, M. Gangwar, and J. P. Dwivedi, “Springback analysis of thin tubes,”
*Journal of Mechanical Engineering*, vol. 7, no. 2, 2011. View at: Google Scholar - 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. View at: Google Scholar

#### Copyright

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.