Effect of Die Design Parameters on Thinning of Sheet Metal in the Deep Drawing Process
1Mechanical Design and Production Department, Faculty of Engineering, Cairo University, Giza, Egypt
Prediction of the forming results, determination of the thickness distribution and the thinning of the sheet metal blank will decrease the production cost through saving material and production time. In this paper, A Finite Element (FE) model is developed for the 3-D numerical simulation of deep drawing process (Parametric Analysis) by using ABAQUS/EXPLICIT FEA program with the proper material properties (anisotropic material) and simplified boundary conditions. The FE results are compared with experimental results for validation. The developed model can predict the thickness distribution and thinning of the blank with the die design parameters (geometrical and physical parameters). Furthermore, with numerical simulation, working parameters can be optimized without expensive shop trials.
At a glance: Figures
Keywords: finite element, parametric analysis, geometrical parameters, thickness distribution, thinning
American Journal of Mechanical Engineering, 2013 1 (2),
Received December 30, 2012; Revised April 14, 2013; Accepted April 15, 2013Copyright © 2014 Science and Education Publishing. All Rights Reserved.
Cite this article:
- Zein, H., et al. "Effect of Die Design Parameters on Thinning of Sheet Metal in the Deep Drawing Process." American Journal of Mechanical Engineering 1.2 (2013): 20-29.
- Zein, H. , El-Sherbiny, M. , Abd-Rabou, M. , & Shazly, M. E. (2013). Effect of Die Design Parameters on Thinning of Sheet Metal in the Deep Drawing Process. American Journal of Mechanical Engineering, 1(2), 20-29.
- Zein, H., M. El-Sherbiny, M. Abd-Rabou, and M. El Shazly. "Effect of Die Design Parameters on Thinning of Sheet Metal in the Deep Drawing Process." American Journal of Mechanical Engineering 1, no. 2 (2013): 20-29.
|Import into BibTeX||Import into EndNote||Import into RefMan||Import into RefWorks|
Currently, the demand for metal containers in the world is in the order of 410 billion units per annum. Out of this, drink cans account for 320 billion and processed food cans account for 75 billion. It is also known, that for the majority of food and drink containers, the cost of the processed metal accounts for 50-70% of the total cost . From the manufacturer’s point of view, a reduction in material usage is of massive importance to cost reductions.
Reduction of material can be done by decreasing the wall thickness of the containers, whilst also retaining adequate strength to allow the container to serve its purpose, without fear of failure. The designer must meet these expectations, but also, must determine the material requirements and properties which are suitable for the food or drink being packaged. Traditionally, this process, as well as most metal forming techniques have been tested experimentally using trial-and-error or empirical methods, which are expensive and time consuming approaches, as dies, blank holders and punches need to be manufactured.
The objective of the work is to successfully simulate the deep drawing process and validate the results, by making a comparison with previously published literature. The literature contains both experimental and numerical results which are in agreement with each other and therefore is being set as the benchmark to for comparison purposes.
The incentive for doing this research is that deep drawing has come to a stage in the current industrialized world that requires the most efficient, low-cost, manufacturing route to be taken at all times.
The successful simulation of the deep drawing process is very much dependent on successful modeling of the deep drawing process for numerical analysis.
By making use of finite element analysis and statistical methods, the prediction of results such as the punch force, the blank holder force, the thickness distribution through sections of the metal and the lubrication requirements can be determined. This can significantly reduce the production costs, for higher quality containers by reducing the lead time to production and provides engineers the ability to promptly respond to market changes. In doing this, the level of knowledge in how various materials interact at the contact surfaces are enhanced, and the data for dealing with specific materials are also increased, which is another positive outcome.
The first of these works is published by Colgan and Monaghan  which has taken a statistical approach, based on experimental design using orthogonal arrays to ascertain which factors most influence the deep drawing process. The geometries for the deep drawing process have been reproduced, and a critical comparison was made for the products. 400 tests were carried out on blanks cut from drawing quality mild steel EN10130 FeP01 of 1mm thickness. The authors explain that for the experimental component of the research, hydraulic cylinders were used to deliver the blank holder force. They also reported that a pressure gauge was used to measure the pressure in the cylinders in conjunction with a needle valve. Although providing adequate results, they suggested, it would be desirable to employ the use of strain gauges in this instance, so as to gain a higher level of accuracy and confidence in the values obtained.
This work established that the punch/die radii have the greatest effect on the thickness of the formed mild steel cups. It was made apparent that the smaller the die radii, the greater the force on the blank, resulting in thinner wall thicknesses. It was also shown that the type of lubrication is very much affecting the force on the punch. Although Colgan and Monahan’s work , concluded that the die radii were noted as a prominent factor, they did not provide enough substance to account for other factors affecting the deep drawing process.
The problem is also addressed, in another approach, by Demirci et al,  for the deep drawing process of AL1050 by performing experiments and finite element analysis using ANSYS/LS-DYNA, investigating the effects of the blank holder forces on differences in the cup wall thickness. The paper shows that under constant pressure the base of drawn cups remains at a constant value until thickness drops dramatically near the cup corners and follows exponential variation towards the outer edge on the cup flange. Taking into account the anisotropic properties of the material they concluded that forces exceeding 10MPa may result in tearing in the cup walls.
A paper published by Waleed Khalid Jawad  confirms the previously stated point that increasing the punch radii can slightly decrease the punch load and vice versa. He investigated the effects of punch radii on maintaining the interfacial contact between the punch and the blank, and punch load on thickness through a section of the drawn cup, and finally predicted the resulting localized strain and stress distribution. He concluded that frictional forces act mainly through the edge of the punch in a shearing manner, with little effect on the flat section.
Vladimirov et al  presented the derivation of a finite strain material model for plastic anisotropy and nonlinear kinematic and isotropic hardening. The work is applied numerically to the drawing of cylindrical and square shaped cups. They show that the numerical simulations, can be suitably extracted from the complex mathematical material models and that the phenomena of anisotropic properties can still be accounted for even with large deformations. They further go on to show the applicability of the mathematical work to account for the occurrence of earring during the drawing process.
Fereshteh-Saniee and Montazeran  predicted forming load, by using various finite element types, and the strain thickness distributions are produced, as reported in . The authors also show that the use of shell 51 elements of the finite element package give a much higher agreement with the experimental results, as compared with the results using visco solid 106 elements. In making these comparisons, Fereshteh-Saniee and Montazeran use a formula derived by Siebel in  to produce an analytical solution to which the level of agreement between the various methods can generally be accepted.
Meguid and Refaat  used a method of variation inequalities to develop an analytical method to model the frictional contact in elastoplastic models which undergo large deformations. The effect of nonlinearities arising due to either the geometry or the materials used in the model is handled by using the updated Lagrange formulation.
Having arrived at a point where the material can be modeled to a very high degree of accuracy and the forming loads can be predicted to complete the drawing process, one may consider various frictional combinations at the contacting surfaces aiming at minimizing such forming loads, and hence cost analysis can be carried out for the manufacturing process of the container.
This paper intends to approach the deep drawing process of thin walled, mild steel, cylindrical containers, by means of a finite element analysis.
2. Simulation of Deep Drawing Process2.1. Finite Element Model
Figure 1 shows a sketch for drawing a circular cup. The process is described in the following steps: The circular blank is inserted between die and blank holder; the punch move downwards; the diameter of the blank decreases with the blank being drawn into the die cavity; the flange of the blank form the perpendicular wall of the cup and the center part of the blank form the bottom of the cup. The important dimensions, Figure 2, of the blank, die, punch, and blank holder are shown in Table 1.
Due to the symmetry, the numerical analysis of the deep-drawing process was performed by using only one quarter of 3D numerical model to reduce the computational time. The model is shown in Figure 3. Discrete rigid form was used to model the punch, die and holder, whose motion was governed by the motion of a single node, known as the rigid body reference node. Die, punch, and holder were meshed with R3D4 elements. Therefore, only the blank sheet metal (224mm diameter x 1mm thickness) was considered deformable with a planar shell base and meshed with reduced integration S4R shell type element .
The blank is made of mild steel . The material is modeled as an elastic-plastic material with isotropic elasticity, using the Hill anisotropic yield criterion for the plasticity to describe the anisotropic characteristics of the sheet metal within the simulation program (ABAQUS/Explicit). Figure 4 shows the true stress, true strain curve of the material behavior, whilst Table 2 shows the material properties.
The die design parameters are classified into two categories.
1st. category is the tool geometry parameters in drawing operations
(1). Die shoulder radius (rd)
(2). Punch nose radius (rp)
(3). Sheet metal thickness (t)
(4). Radial clearance (wc)
2nd category is the process operating parameters
(5). Blank holder force (BHF)
• Coefficient of friction between punch/blank (μp),
• Coefficient of friction between holder/blank (μh),
• And coefficient of friction between die/blank (μd)
The study of the influence of the tool geometry parameters and the physical parameters on thinning of sheet metal in the deep drawing process have been made, according to recommended geometrical values of these parameters from the code, Table 3.
3. Validation of The Model
Colgan and Monaghan , performed an experimental and numerical analysis of the process and therefore a comparison between the present work against the values which have been recorded for the blank thickness of the deep drawn cup may be used to validate the model. Figure 5 shows the typical shape characteristics and thickness variations of a partially drawn blank. Thickening occurs in the flange area and some stretching and necking occurs in the sidewall, just above the punch nose radius .
They measured the sheet metal blank thickness at 8 locations. Therefore, 8 values of thickness were taken from the generated data, at distances corresponding to these lengths as shown in Figure 6 and Figure 7.
The values obtained from the present results correspond to the thickness at the closest node to the corresponding locations given by Colgan et al.  along the cup cross section. Obviously these locations cannot be exactly the same due to automatic mesh generation, and therefore small differences in thickness variations are expected. Figure 8 depicts the location of reference nodes given by Colgan et al. , whilst both locations and values from the present model are demonstrated in Figure 9.
It is shown in Table 6 that the average thickness distribution in the blank of the present FEA model is more close to the average obtained from the experimental results . Furthermore the present results are very close to their numerical results with differences which do not exceed 2.9 %, on average. It is even more closer to their experimental results with differences less than 0.6 % on average.
4. Results and Discussion4.1. Tool Geometry Parameters
4.1.1. Die Shoulder Radius (rd)
The geometry of die influences the thickness distribution and thinning of sheet metal blank in the deep drawing processes. Table 7 and Figure 10 show thickness distribution for different values of the die shoulder radius (rd). Figure 11 shows thinning of sheet metal with diverse values of the die shoulder radius (rd).
4.1.2. Punch Nose Radius (rp)
The geometry of punch influnces the thickness distribution and thinning of sheet metal blank in the deep drawing processes. Table 8 and Figure 12 show thickness distribution with different values of the punch nose radius (rp). Figure 13 shows thinning of the sheet metal with diverse values of the punch nose radius (rp).
It is shown that for a punch nose radius (rp) that is less than three times the thickness of the blank (t), the cup fails due to increased thinning, whilst for (rp) greater than 3t, thinning is slightly stable.
4.1.3. Blank Thickness (t)
The original blank thickness has some effect on the thickness distribution and thinning of sheet metal blank in the deep drawing processes. Table 9 and Figure 14 show thickness distribution with several values of the blank thickness (t). Figure 15 shows the percentage of thinning of sheet metal with different blank thickness (t).
It has been shown that the average distribution of the blank thickness is increasing with increasing of the blank thickness. Also, the % of thinning is increasing with increasing of the blank thickness. Taking into account, the blank thickness and the punch diameter affect the Limiting drawing ratio (LDR) decreases as the relative punch diameter increases; see the Figure 16 . Slightly thicker materials can be gripped better during the deep drawing process. Also, thicker sheets have more volume and hence can be stretched to a greater extent with increasing in thinning.
4.1.4. Radial Clearance (wc)
Radial clearance is an important parameter, formulated as the difference between die radius and punch radius (wc = DR – PR). Table 10 and Figure 17 show thickness distribution with different values of the radial clearance (wc). Figure 18 shows thinning of the blank with diverse values of the radial clearance (wc).
It is shown that the distribution in sheet metal thickness is increasing with reducing the radial clearance (wc). In addition, for the radial clearance (wc) that is less than the blank thickness (t), the cup fails due to increased thinning. Whilst for the radial clearance (wc) greater than the blank thickness (t), thinning is stable. The radial clearance which is less than (0.5t) is not acceptable because the percentage of reduction in thickness is more than 45%, while the maximum allowable percentage of reduction in thickness is 45% .4.2. Physical Parameters
4.2.1. Blank Holder Force (BHF)
The blank holder force (BHF) required to hold a blank flat for a cylindrical draw varies from very little to a maximum of one third of the drawing pressure . Table 11 and Figure 19 show thickness distribution with different values of the blank holder force (BHF). Figure 20 shows thinning of the blank with the variation of the blank holder force (BHF).
The higher the blank-holder force, the greater will be the strain over the punch face, however the process is limited by the strain in the side-wall. If the tension reaches its maximum value, the side-wall will fail by splitting . It has been shown that the cup collapse due to thinning with the increase of the blank holder force (BHF) over 0.5 ton (recommended value from code).
4.2.2. Friction Parameters
In Sheet Metal Forming processes, such as deep drawing, friction plays an important role. Together with the deformation of the sheet, the friction determines the required punch force and the blank holder force. Consequently, the friction influences the energy which is needed to deform a sheet material. Friction also influences the stresses and strains in the work piece material and, hence, the quality of the product. Therefore, it is important to control the friction between the tool and the work piece.
During deep drawing, hollow bodies are produced from metal blanks using punches, dies and blank holders. In no other forming operation are the friction and, as a result, the lubricating conditions so complex. In one drawing operation, a particularly low coefficient of friction is required in one area and a particularly high coefficient of friction in another.
The friction conditions when drawing can be best explained by means of the example of drawing a cylindrical cup. Figure 21 shows the different friction areas and generally illustrates the deep drawing operation .
184.108.40.206. Coefficient of Friction Between Punch/Blank (μp)
The drawing force in the flange, necessary to form the sheet metal, is applied by the punch on the base of the cup and transferred through the wall into the flange. This transmission of force calls for the highest possible coefficient of friction on the punch edge. This demonstrates the first significant rule for lubrication. In this case: neither the punch nor the sheet metal blank should be lubricated in this area. Even if this consideration was optimal for an isolated forming operation and force transfer, one still must consider punch wear. This friction area requires lubrication in the boundary mode with a high coefficient of friction and anti-wear behavior .
Table 12 and Figure 22 show thickness distribution with different values of the coefficient of friction between punch/blank (μp). Figure 23 shows thinning of the blank with the variation of the coefficient of friction between punch/blank (μp).
For fluid lubricant, it has been shown that the distribution in sheet metal thickness is decreasing with increasing of the coefficient of friction between punch/blank (μp). For solid lubricant and dry lubricant, the distribution in sheet metal thickness is stable with increasing the coefficient of friction between punch/blank (μp).
For fluid lubricant, it has been shown that the thinning of sheet metal is increasing by small values with increasing the coefficient of friction between punch/blank (μp). On the other hand, for solid lubricant and dry lubricant, the thinning growing with rising of the coefficient of friction between punch/blank (μp).
220.127.116.11. Coefficient of Friction Between Holder/Blank (μh)
The main aim of deep drawing lubrication is to achieve minimum friction in the blank holder area. This includes the lubrication on the blank holder side, which, as much as possible, should be a lubricant ring in order to reduce the friction in the region of the punch edge as little as possible. However, in the case of high viscosity oils and drawing pastes, if the applied lubricant film is excessive, there is the risk that hydrostatic effects on the blank holder side of the flange can cause reduction of sheet metal contact with the blank holder resulting in wrinkling .
Table 13 and Figure 24 show thickness distribution with different values of the coefficient of friction between holder/blank (μh). Figure 25 shows thinning of the blank with the variation of the coefficient of friction between holder/blank (μh).
It is noted that the suitable lubricant for the holder/blank friction is the fluid lubricant. When the coefficient of friction between holder/blank (μh) is rising, the average value of the thickness distribution is slightly stable. But the thinning of the cup is growing with rising of the coefficient of friction between holder/blank (μh).
18.104.22.168. Coefficient of Friction Between Die/Blank (μd)
In this case the aim of lubrication is to achieve a minimum coefficient of friction with minimum wear. Even excessive lubrication of the sheet metal surface on the die side normally does not create a problem. Only on parts with large surface areas and a high percentage of stretch drawing, such as car body parts can excessive lubricant quantities lead to unacceptable deviations in form .
Table 14 and Figure 26 show thickness distribution with different values of the coefficient of friction between die/blank (μd). Figure 27 shows thinning of the blank with the variation of the coefficient of friction between die/blank (μd).
These results show that the proper lubricant for the die/blank friction is the fluid lubricant. In addition, the average value of the thickness distribution is reduced with increasing the coefficient of friction between die/blank (μd). Finally, the thinning in the drawn cup is increasing with increasing the coefficient of friction between die/blank (μd).
Successful deep drawing depends on many parameters (geometrical parameters and physical parameters). The finite element analysis simulation (FEAS) is a powerful tool in the sheet metal die design, reducing time and trial and error efforts. A Finite Element Model is developed to reach the optimum solution without many costly trials of production which is a common practice in the traditional production approaches.
The results show that:
1. The die shoulder radius is recommended to be about 8 times sheet thickness.
2. The punch nose radius is recommended to be greater than 3 times sheet thickness.
3. The thicker sheet metal is softer due to its increased volume which increases the thinning in sheet metal.
4. The radial clearance is recommended to be greater than the value of the sheet thickness [If the clearance is not large enough, ironing or thinning will occur].
5. Blank holder force (BHF) is recommended to be less than 3 tons to avoid increased thinning.
6. The fluid lubricant with (μ p= 0 to 0.3), is more suitable to apply between punch and blank to reduce the thinning of the cup.
7. The fluid lubricant with (μ h= 0 to 0.2), is more suitable to apply between holder and blank. (μ h) should be about 0.1 to reduce the thinning of the cup.
8. The fluid lubricant with (μ d= 0 to 0.2), is more suitable to apply between die and blank. (μ d) is recommended to be about 0.1 to reduce the thinning of the cup.
|||Mark, Richard Coles, Kirwan, J. , Food and Beverage Packaging Technology, Blackwell Publishing, London, 2011, 344.|
|||Colgan, M., and Monaghan, J., "Deep Drawing Process: Analysis and Experiment," Journal of Materials Processing Technology, 132, 35-41, 2003.|
|||Ibrahim Demirci, H., Yaşar, M., Demiray, K., "The Theoretical and Experimental Investigation of Blank Holder Forces Plate Effect in Deep Drawing Process of AL 1050 Material," Journal of Materials & Design, 29(2), 526-532, 2008.|
|||Jawad, Waleed Khalid "Investigation of Contact Interface between the Punch and Blank in Deep Drawing Process," Journal of Engineering and Technology, 25(3), 370-382, 2007.|
|||Vladimirov, Ivaylo N., Pietryga, Michael P. Reese, Stefanie, "Anisotropic Finite Elastoplasticity with Nonlinear Kinematic and Isotropic Hardening and Application to Sheet Metal Forming", International Journal of Plasticity, 26,659-687, 2010.|
|||Fereshteh-Saniee, F., Montazeran, M.H., "A Comparative Estimation of the Forming Load in the Deep Drawing Process," Journal of Materials Processing Technology, 140, 555-561, 2003.|
|||Lange, Kurt, Handbook of Metal Forming, Society of Manufacturing Engineers, 1985.|
|||Meguid, S.A., Refaat, M.H. "Finite Element Analysis of the Deep Drawing Process using Variational Inequalities," Journal of Finite Elements in Analysis and Design, 28, 51-67, 1997.|
|||Abaqus Analysis User’s Manual, Volume IV: Elements, 2011.|
|||Taylor, L. M., Cao, J., Karafillis, A. P., Boyce, M. C. "Numerical Simulations of Sheet Metal Forming," Proceedings of 2nd International Conference, NUMISHEET 93, Isehara, Japan, Ed. A. Makinovchi, 1993.|
|||Suchy, Ivana Handbook of Die Design, 2nd Edition, McGraw-Hill Publishing, 2006, 383- 400.|
|||Dieter, George E. Mechanical Metallurgy, SI Metric Edition, McGraw-Hill Series in Material Science and Engineering, 1988, 669.|
|||Grote K.-H., Antonsson E.K., Springer Handbook of Mechanical Engineering, Springer Science and Business Media, 2008, 617.|
|||David S., Jeff L., John G., Fundamentals of tool design, American Society of Tool and Manufacturing Engineers. New Delhi: Prentice Hall, 2003.|
|||Marciniak, Z. Duncan, J.L. Hu, S.J. Mechanics of Sheet Metal Forming, 2nd Edition, published by Butterworth-Heinemann, 2002, 53.|
|||Mang, Theo, Dresel, Wilfried Lubricants and Lubrication, 2nd Edition, WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, 2007, 527-530.|