Elastic-plastic solution to stamping thin strip on elastic foundation
LI Yu-ping(李玉平)1, ZHOU Li-qun(周里群)1, ZHOU Yi-chun(周益春)2
1. School of Mechanical Engineering, Xiangtan University, Xiangtan 411105, China;
2. Key Laboratory of Advanced Materials and Rheological Properties of Ministry of Education,
Xiangtan University, Xiangtan 411105, China
Received 2 June 2005; accepted 16 September 2005
Abstract: An analytical method was proposed to solve the mechanical problems of stamping a thin strip on an elastic foundation. The thin strip was divided into four parts according to its deformation and contact with the punch and the elastic foundation, especially an elastic-plastic part was considered in the deflection of the thin strip. Analytical solutions were derived individually for each part and two models were established with the help of elastic and plastic large deflection theories. Compatibility conditions between the neighboring parts of the thin strip constructed the non-linear equation group. Solutions were carried out by programming with a software. The deformation shape, the membrane force, and the moment and shear force of the deformed thin strip were obtained. The results of the two models were compared. The study shows that the method is effective.
Key words: elastic-plastic solution; beam bending; thin strip stamping; contact mechanics
1 Introduction
In the recent years, there has been growing need for improving quality and decreasing costs of products of the metal forming industry[1-3]. The thin strip metal forming with elastic foundation has been widely used in engineering. A rigid punch presses the thin strip into a pliable polyurethane pad which is displaced, and applies compressive force onto the work-piece surface. This problem is a mechanical problem relating to three-body contact. Contact mechanics has been studied extensively [4-6], but this problem is still challenging because the deformation of the system involves strong non-linearity of geometry and materials.
ZHANG et al[7] and LIN et al[8] studied the problem of stamping a thin strip on elastic foundation, but their studies were only limited to elastic deformation. ZHANG[9] also studied the elastic-plastic deformation of thin strip metals stamped by rigid punch, but the research was only suitable for small deflection. Geiger et al[10] investigated the bending of U-shaped parts with a rigid punch and an elastic die by use of finite element method, but it required great user-expertise and computational cost. Till present, any analytical solution to the problem has not been very satisfactory, mainly because:
1) The forming process is of large deflection with non-linear material, thus any linear elastic model or small strain elastic-plastic model is not consistent with practice.
2) In analytical approach, the distributions of the interface contact stresses are assumed a priori, which limits the applicability of the solution in terms of the indenter profiles, the ratio of the thickness of the thin strip to punch radius, and the variation of the material properties of the thin strip and elastic foundation.
On the basis of former researches, the authors proposes a new mechanical model which considers large elastic-plastic deflection and non-linearity of materials for the process of thin strip stamped on elastic foundation, and corresponding computation method is performed.
2 Mechanical model
Consider the stamping of a thin strip (thickness: h) on an elastic foundation by a half-circular rigid punch (radius: a). The deformation of the thin strip is in plane stress and is symmetrical to its central line, as shown in Fig.1. The contact zones and interface stresses between the punch and thin strip and the thin strip and foundation are unknown in advance, which are functions of the punch stroke in the stamping process. Because of the symmetry of the deformation, however, we can study half of the thin strip only. In order to analyze the problem properly, we divide the thin strip into four parts as shown in Fig.2:
1) The contact part AB, where the thin strip is in perfect contact with both the punch and the foundation surface;
2) The elastic-plastic part BC, where the thin strip is in perfect contact with the foundation but has no contact with the punch, its deformation is elastic-plastic;
3) The elastic part CD, just like part BC but the deformation is only elastic;
4) The free part DF, where no contact takes place with either the punch or the foundation.
Fig.1 Stamping thin strip on elastic foundation
Fig.2 Deformation process and four parts of thin strip
The curvature of the contact part is a known function which is identical to that of the punch outer surface. The contact stresses on the thin strip of this part are the normal contact stresses between the thin strip and punch, Qnp, and the normal one between the thin strip and elastic foundation, Qne. The frictional force is neglected as compared with Qnp, Qne, as shown in Fig.3.
The elastic-plastic part and the elastic part can be modelled as a cantilever beam subjected to normal stress Qne, due to the interaction between the thin strip and foundation as shown in Fig.4. The end B of this part is the contact-off point between the punch and the thin strip, and the end D is the contact-off point between the thin strip and the elastic foundation. The boundary conditions of the cantilever must be specified to guarantee the continuity of stresses and deformation across these two ends, i.e.
1) The bending moment, membrane and shear force are zero at point D.
2) The deflection and its slope, bending moment and membrane force are equal to those of the contact part AB at B.
Fig.3 Mechanical model for contact part AB
Fig.4 Mechanical model for parts BC, CD and DF
The free part of the thin strip, DF, does not deform during the stamping process. Its displacement relies on the deflection and deflection slope at end point D, and can be calculated easily when the solution to CD is obtained. Hence, we will ignore this part in the following analysis.
To simplify the computation of the contact stresses, we assume that the normal reaction of the elastic foundation follows the Winkler’s hypotheses.
2.1 General mechanical description
It is postulated that the thin strip metals can be characterized by the ideal elastic-plastic model[6]. The elastica of the thin strip bending of large deflection can be expressed as[11]
M/Me≤1 (1)
The plastica of the thin strip bending of large deformation can be expressed as[12]
1≤M/Me≤ (2)
where X and Y are Cartesian coordinates, θ is the tangent angle of the deformed thin strip surface with the positive direction of X-axis(Fig.4), dS is the length of an infinitesimal thin strip element(Fig.4), E is the elastic modulus of the thin strip material, I is the second moment of area per unit width of the thin strip cross section, M is the bending moment of the thin strip, R is the curvature radius of the thin strip, and Me is the bending moment of elastic limit. The equilibrium equation of the thin strip with large deflection can be written as[13]
(3)
where T is the membrane force of the thin strip, N is the shear force of the thin strip, Qn is the normal stress, Qn=Qnp–Qne, and is the included angle of the external normal of the deformed thin strip surface at θ with the Y-axis.
Above Eqns.(1)-(3) can also be written as non-dimensional forms to simply the solutions as follows:
(4)
(5)
(6)
(7)
where
,
where σy is the yield stress of the thin strip, and Ed is the elastic modulus of the elastic foundation.
2.2 Solution to part AB
The deformation of the thin strip in this zone exactly follows the geometrical profile of the punch, and its deformation is plastic. Using the geometrical equation of a circle, we have
x=sin θ, y=–cos θ, r=1 (8)
w=wr+cos θ (9)
where w is the perpendicular displacement at any point of part AB, wr is the punch displacement at the punch center(Fig.2). The force in the thin strip can easily be determined when Eqns.(8) and (9) are substituted into Eqns.(6) and (7) and noticing .
n=0, t=t0 (10)
where t0 is a dimensionless membrane force at the central point A of the thin strip. If t0 is known, the normal contact stress between the thin strip and the punch is
(11)
2.3 Solution to part BC
When a force (or several forces) acts on a beam, the produced bending moment can always be written as the linear function of the coordinates x and y[12], we use this relation here.
(12)
where A, B and C are different constants, which are related to the quantity and position of the forces.
Using the third equation of Eqn.(5), we have
Thus
(13)
Substituting of Eqn.(12) into Eqn.(13), and then using the first two equations of Eqn,(5), yield
(14)
where C1 and C2 are integration constants.
Substituting Eqn.(14) into the first two equations of Eqn.(5), thus we obtain the thin strip equation:
(15)
(16)
where C3 and C4 are also integration constants.
According to the third equation of Eqn.(6), we have the shear force of the thin strip in part BC,
(17)
Using the second equation of Eqn.(6), the membrane force of the thin strip in part BC is,
(18)
2.4 Solution 1 to part CD
Only elastic deformation appears in this part, the bending deflection of this part is far smaller than its length. Even though the thin strip has deformed, it can approximately be seen as a straight line. Hence, the coordinate x of the deformed thin strip is proportional to its coordinate y. Thus according to Eqn.(12), the produced moment of the thin strip in this part can be written as
(19)
where g and p are constants. We use this expression for distributed load.
Substitution of Eqn.(19) into the third equation of Eqn.(4), yields
Notice
Thus
Integrating this equation, and noticing C5=2η1g? cos θd, we have
(20)
Assuming
then yields
(21)
In point C, .In point D, . Thus integration of Eqn.(21) may obtain the arc length of part CD,
(22)
In Eqn.(22),
is called the complete elliptical integration of the first kind; and
is called the incomplete elliptical integration of the first kind. Substitution of Eqn.(21) into the first two equations of Eqn.(4), yields
(23)
(24)
where xc and yc are coordinates of point C in terms of Eqns.(15) and (16), and
is called the incomplete elliptical integration of second kind.
In point D, =π/2, and its longitudinal coordinate yd=wr. According to Eqn.(24), we have
(25)
Because D is a free point, no drawing force acts on it, and the pressure of the elastic foundation is perpendicular to the thin strip, thus the membrane force of the thin strip in part CD in fact is zero. That is
t=0 (26)
The moment in point D is zero, and the longitudinal coordinate yd=wr. Substitution of these conditions into Eqn.(19) may attain the moment of the thin strip in part CD,
(27)
According to the third equation of Eqn.(5), we have the shear force of the thin strip in part CD,
(28)
2.5 Solution 2 to part CD
The numerical results of the solution 1 to part CD show that the elastic part CD is very straight. Experience tells us that the elastic part is also straight. So as to simplify the computation, and to avoid special functions, we further assume the elastic part CD is a straight line. That is to say, the tangent angle of the thin strip in part CD θ=θc=θd. Thus according to the first two equations of the Elastica 4, we may obtain the coordinates of the part CD,
(29)
(30)
where s is the arc length of the thin strip in part CD from the beginning of point C.
The moment of the thin strip in part CD is still the Eqn.(27). Because the point C is a connective point from elastic state into plastic state, thus the moment of the thin strip in this point is the elastic limited moment. So when y=yc, we have . Substituting this condition into Eqn.(19), we obtain the moment of the thin strip in part CD
(31)
And in terms of Eqn.(28), the shear force of the thin strip in the part CD is
(32)
The membrane force of the thin strip in part CD is zero, too.
3 Numerical method
Above analyses may be carried out by programming, but the constants which appear in the above equations should be determined at first. They may be decided by the compatibility conditions of the parts of the thin strip.
3.1 Numerical model 1
This model uses the solution 1 of the part CD. In the point B of part AB, the tangent angle θ=θb. According to Eqns.(9) and (11), we have the expressions of
In the point B of part BC, also θ=θb. According to Eqns.(12), (15)-(18), we have the expressions of .
In the point C of part BC, the tangent angle θ=θc. According to Eqns.(12), (16)-(18), we have the expressions of .
In the point C of part CD, also θ=θc, .
Using the solution 1 of part CD, that is in terms of Eqns.(26)-(28), we have the expressions of .
According to Eqns.(23) and (24), the compatibility of coordinates in the point C of part BC and part CD is satisfied, but according to Eqn.(25), we can get another equation of .
From above expressions, we may construct 10 non-linear equations as listed in Table 1. In the 10 non-linear equations, there are 11 unknown parameters, so it is necessary to add another equation. From Eqn.(11) we know that at θ=0, the normal contact stress qnp between the thin strip and the punch is maximal. It will affect the plastic forming of the thin strip, and results in the contact part AB being in plastic state. So we assume that when θ=0, qnp=1.Thus there are 10 unknowns in the 10 non-linear equations in Table 1. The problem may be solved.
Table 1 Equations and unknowns of numerical model 1
3.2 Numerical model 2
This model uses the solution 2 of the part CD. Part equations in the above numerical model 1 are still valid here. According to the solution 2 to part CD, we can obtain the additional boundary conditions in the point C of part CD are . Thus we can attain 8 non-linear equations as listed in Table 2. Unknowns in the equations are just 8, that is . The model can be solved.
Table 2 Equations and unknowns of numerical model 2
4 Numerical results
The proposed models may be accomplished in Matlab software with programming. The non-linear equations may be solved by a least squares method, and the numerical integral is calculated by the adaptive Simpson quadrature[14]. As an numerical example, the thin strip is a low carbon steel of No.15, its yield stress σy=226 MPa, and its elastic modulus E=200 GPa; the elastic foundation is a kind of polyurethane 80 A, its elastic modulus Ed=30 MPa. The punch radius a=50 mm. The thin strip thickness h=3.0 mm. The computed results are shown in Figs.5-10 at the punch stroke wr=1.
Fig.5 Deflection process and shape of thin strip of model 1
Figs.5 and 6 show the deflective shapes of the thin strip of model 1 and model 2 when the punch stroke is 1. It can be seen that the deflective shape of the thin strip just likes a V word. They obviously have three parts, including the contact part, the elastic-plastic part and the elastic part. Furthermore, the length of the elastic-plastic part is not very short, it can not be neglected. It is clear that the elastic part is very straight, which matches the case of practice. The span of the V-shape is related to the elastic modulus of the elastic foundation. The difference of model 1 and model 2 reflected to the deflection of the thin strip is that the deflection of the thin strip is smooth in model 2, but not smooth and has a inflexion in model 1. That is to say, model 2 seems to be more reasonable.
Fig.6 Deflection process and shape of thin strip of model 2
Fig.7 Moment and shear force of deformed thin strip of model 1
Fig.8 Moment and shear force of deformed thin strip of model 2
Figs.7 and 8 show the moments and shear forces of the deformed thin strip. The curve for shear force is ideal. Theoretically, the shear force for contact part AB should exactly be zero from Eqn.(11), and at the end point D of the elastic part it is also zero since DF part is free. These are all satisfied in the computed curve. This curve also shows that the shear force in the elastic-plastic part and the elastic part is very small.
The moment of the contact part AB is a constant, and the moment of the elastic part CD is a linear function of the coordinates. These are all met. The moment of the elastic-plastic part BC should satisfy the boundary conditions at both points B and C, but computation shows they are very difficult to meet. To speak exactly, the convergent precision of the equations y5 and y9 in model 1 and y5 and y8 in model 2 are not very good. In the optimal processes of solving the non-linear equations, we already have magnified these equations 10 times. If they are magnified more, it will deteriorate the convergent precision of other equation. The discontinuity error of the moment in both points B and C is about 20%. The cause resulting in this instance seems to be that Eqn.(12) is not strict under the distributed load.
The membrane forces of the thin strip are shown in Figs.9 and 10. The curves are continuous. The membrane force of the contact part AB is a constant, and it is zero in the elastic part CD. The curves for the elastic-plastic part are smooth. The numerical results are good. The difference between model 1 and 2 is that the arc length of corresponding part is different. In model 1, the arc length of the contact part is about 0.6, the elastic-plastic part is about 0.63. In model 2, the arc length of the contact part is about 0.6, the elastic-plastic part is about 0.35. The cause resulting in this instance is because the convergent precision of the two models is different. Since the convergent precision of model 2 is 10-4, and better than model 1, the results of model 2 are believable and more accurate.
Fig.9 Membrane force of deformed thin strip of model 1
Fig.10 Membrane force of deformed thin strip of model 2
5 Conclusions
An analytical method is proposed to solve the mechanics problem of stamping a thin strip on an elastic foundation. Two models are established with the help of elastic-plastic theory of large deformation. Especially the elastic-plastic part is considered in the deflection of the thin strip. The results of the two models are compared. The result of the simplified model is more satisfactory. The method adopted in this paper may be extended to consider clad thin strip’s stamping.
References
[1] LI R D, ZHOU P. Numerical modeling of microstructure forming process for Al-Al3Fe eutectic alloy[J]. Trans Nonferrous Met Soc China, 2003, 13(4): 849-854.
[2] Wang X Y, Wang Z J, Wang Z R. Numerical simulation of aluminum alloy ladder parts with viscous pressure forming[J]. Trans Nonferrous Met Soc China, 2003, 13(2): 391-397.
[3] Zhang X M, Chen M A. Influence of process parameters on hybrid forming of aluminum sheet[J]. Trans Nonferrous Met Soc China, 2001,11(6): 879-883.
[4] Gladwell G M L. Contact Problems in Elasticity[M]. Beijing: University of Science and Technology of Beijing, 1991.
[5] Pempsey J P, Zhao Z G, Li H. Axisymmetric indentation of an elastic layer supported by a winkler foundation[J]. Int J Solids Structures, 1991, 27(1): 73-87.
[6] Girija C V, Vallabhan, Das Y C. A refined model for beams on elastic foundations[J]. Int J Solids Structures, 1991, 27(5): 629-637.
[7] Zhang L C, Lin Z Q. A new mechanics model of stamping a thin strip on an elastic foundation[J]. Int J Solids Structures, 1997, 34(3): 327-339.
[8] Lin Z Q, Zhang J W, Yu X B. Mechanical modeling of cold metal sheet forming process on elastic deformable die[J]. Chinese Journal of Mechanical Engineering, 1999, 35(1): 61-67.
[9] Zhang L C. A mechanics model for sheet-metal stamping using deformable dies[J]. Journal of Materials Processing Technology, 1995, 53: 798-810.
[10] Geiger M., Engel V, Ende A V. Investigations on the sheet bending process with elastic tools[J]. Journal of Materials Processing Technology, 1991, 27: 265-277.
[11] YU T X, ZHANG L C. Plastic Bending: Theory and Applications[M]. Beijing: Science Press, 1992.(in Chinese)
[12] YU T X, Johnson W. The plastica: the large elastic-plastic deflection of a strut[J]. Int J Non-linear Mechanics, 1982, 17: 195-209.
[13] ZHOU L Q. Impact Performance and Stamping Forming Processes Analysis for Electrodeposited Nickel Coating[D]. Xiangtan: Xiangtan University, 2004. (in Chinese)
[14] WANG M R. MATLAB 5.x and Scientific Calculation[M]. Beijing: Tsinghua University Press, 2000. 145-149, 242-252.(in Chinese)
Foundation item: Project(04JJ40035) supported by the Natural Science Foundation of Hunan Province, China; Project(2003AA331090) supported by the High-tech Research and Development Program of China
Corresponding author: ZHOU Li-qun; Tel: +86-732-2371301; Fax: +86-732-8292210; E-mail: Johnzlq@163.com
(Edited by YUAN Sai-qian)