Skip to content
BY-NC-ND 4.0 license Open Access Published by De Gruyter Open Access November 10, 2018

Failure prediction of axi-symmetric cup in deep drawing and expansion processes

  • Ján Slota EMAIL logo , Marek Šiser and Ivan Gajdoš
From the journal Open Engineering

Abstract

Presented work deals with the prediction of the forming limit of food can obtained with deep drawing, reverse drawing and expansion operations. Two commonly used materials of TH330 steel and AA5352 aluminium alloy for packaging can production were studied. Finite element simulation (FEM) is an essential tool in the packaging industry to prevent different sheet metal forming difficulties such as failure under complex nonlinear strain paths and plastic anisotropic earing. In order to characterise the material plastic properties and to specify failure criteria, the uniaxial tensile, hydraulic bulge test, as well as the routines for obtaining forming limit curves were carried out. The input material data required for various material models are also described. Utilisation of advanced material models in numerical simulation require a large number of input data. Prediction of failure location in drawing and expansion of axi-symmetric cups were estimated for each material. The FEM results were verified by real experiments.

1 Introduction

The deep drawing process is one of the most used progressive forming methods with a wide range of applications in different fields such as packaging, automotive and aerospace industries. In this process, a sheet metal is transformed into a hollow part using forming tools whose geometry is consistent with the shape of the part. Several process variables influence the success of a deep drawing process. These include the tool geometry, the blank-holding force, the friction conditions, mechanical properties of the blank and others.

In recent years, many research works have been conducted to study the influence of process variables on forming process using finite element simulations.

Numerical simulation of deep drawing axi- symmetric cups provides significant information for processors in the packaging industry. Thanks to FEM simulation, we can predict many different errors that greatly affect the cups production. These errors include extreme thinning, earing, which is affected by plastic anisotropy, wrinkling, binding, cosmetic defects on cups and cracks or fractures started from various combinations of stress states on the cup, such as, for example, during deep-drawing and cup expansion.

Based on the background knowledge and experiences [1, 2, 3, 4], it is important to select a suitable material model to predict the material behaviour during plastic deformation to achieve adequate simulation results. One of the most significant factors that are included in the material model is the hardening curve, the yield function, and the forming-limit curve. In order to achieve a better correlation of prediction with experimental results, several works have focused on improving the quality of constitutive models. Authors [5, 6, 7] focused on the effect of the yield function, in the works [8, 9] the influence of the hardening curve was examined. A large number of research works [10, 11, 12, 13] deal with predicting failures, where forming limit diagram (FLD) is used. In many cases, the authors deal only with deep drawing processes. However, there is relatively little research on combined processes such as deep drawing, reverse drawing and expansion. Our study aims at predicting sheet behaviour after deep drawing, reverse drawing and expansion using FEM simulation. The main objective of this study was focused on predicting the failure point after the expansion operation.

2 Methodology

Experiments with two different materials were carried out. Medium strength steel grade TH330 is double reduced and tin coated sheet, specifically intended for round or rectangular cans. It is the most common material used in can industry [14]. The alternative is aluminium alloy AA5352, used for packaging production processes, but also in automotive and electrotechnical industry. This aluminium alloy has good formability and high hardening rate. Due to the presence of magnesium, it is susceptible to oxidation and has to be heat treated [15]. In order to determine input data required for FEM simulation, basic mechanical and technological tests were carried out. In the present study, the material properties were characterized by both tensile and hydraulic bulge tests. The tensile test was carried out according to standard ISO 6892-1:2016. From the test, basic mechanical properties as yield point, ultimate tensile stress, the strain-hardening exponent, uniform elongation, Lankford’s coefficients and others material parameters were obtained. The mechanical properties obtained from the tensile test are shown in Table 1.

Table 1

Basic mechanical properties of studied materials

MaterialDirectionYieldTensileElongation [%]Normal
Strength [MPa]Strength [MPa]UniformTotalanisotropy
AA53520197.59239.298.21911.2300.535
45193.94231.549.58013.0901.105
90198.74239.0411.03216.1102.270
TH3300258.87366.8715.88129.0421.449
45255.25371.0215.66529.5051.266
90249.42367.0815.77129.3471.510

Thin sheet metal usually splits at a strain of approximately 10% for aluminium and 20% for steel in a uniaxial tensile test. In order to obtain the material data at higher strains, hydraulic bulge test is commonly used to extend the tensile test data. The use of this test contributes to reducing discrepancies on data extrapolation of hardening curves too. Scheme of the hydraulic bulge test is shown in Figure 1. Following equations were used to calculate the equivalent stress:

Figure 1 Scheme of the hydraulic bulge test
Figure 1

Scheme of the hydraulic bulge test

(1)σb=pR2t
(2)t=t0expεt

where σb is effective stress [MPa], p is hydraulic pressure [MPa], R is the radius of curvature [mm], εt is current thickness strain, t and t0 are actual and initial thicknesses [mm].

Using values from hydraulic bulge test the strain-hardening curve can be established. The parameters for determination of hardening curve and results of hydraulic bulge test are shown in Table 2.

Table 2

Material properties obtained from hydraulic bulge test and Krupkowski hardening curve

MaterialThickness (t) [mm]Hydraulic bulge testKrupkowski model σ=K(ε0p)n
σbεtKε0n
[MPa][-][MPa][-][-]
AA53520.279304.240.1200362.70.01610.1535
TH3300.270579.860.4642575.80.00450.1578

In sheet metal forming, the FLD is commonly used to predict material failure in forming operations. The failure of sheet material can be induced by strain paths ranging from uniaxial tension to plane strain and biaxial tensile loading (different ratios between major and minor strains). Forming limit curve (FLC) represents the major and the minor strains in the plane of a sheet and corresponds to the maximum admissible local strains achievable just before the occurrence of visible defects in the sheet metal like fracture or necking.

FLC’s were determined according to the ISO 120042:2008 for both sheet materials. The universal sheet metal testing machine Erichsen 145-60 was used for performing the tests. Experimental samples with different radii were cut and Nakajima test with punch diameter of 100 mmwas applied. The machine can provide controllable slide velocity, and punch speed was set to 1.5 mm/s. The test was measured by 3D optical system Aramis M5. The experimental device is shown in Figure 2.

Figure 2 Device for estimating of FLC by Nakajima test
Figure 2

Device for estimating of FLC by Nakajima test

The experiment consisted of three operations. Cup was formed in a two-stage process of drawing and reverse redrawing using a combination tool (Figure 3), which occur sequentially during a single stroke of the drawing operation punch. The formed cup was subsequently subjected to die expansion operation. The diameter of the blank was Ø162.96 for both materials. Die extension operation was realized in the separate tool (Figure 4). Inner die diameter in deep drawing operation was Ø116.84 mm. Cup without flange with an average height of 27.5 mm was drawn. After cup was fully drawn followed redrawing operation. As it is shown in Figure 3, in the combination tool the punch from the first operation was used as die in the second operation. Inner die diameter was Ø87.85 mm and cup was drawn without failure with approximately 52.5 mm of height.

Figure 3 Combination tool for deep drawing and redrawing operations
Figure 3

Combination tool for deep drawing and redrawing operations

Figure 4 Experimental tool for expansion operation
Figure 4

Experimental tool for expansion operation

After drawing operations, the die expansion operation was performed. From practical knowledge we know, that force required to expand cup up to failure will cause the material flow to the bottom of the cup. This cause wrinkling in the cup bottom and reduce cup height. To avoid this, the bottom of the cup was clamped by inner clamp plate attached by a bolt to cup support tool. The preformed cup was placed between inner clamp plate and cup support tool as shown in Figure 4. In the centre of the bottom of cup is placed hole with Ø9 mm diameter. Clamping plate was screwed to cup support tool (Figure 4). Bleed at punch prevents accumulation of air. The same clamp arrangement is used for both materials - steel and aluminium alloy. Some process parameters defined in FEM code are shown in Table 3.

Table 3

Parameters defined in Pam-Stamp 2G software

UnitsDeep drawingRedrawingExpansion
Drawing speed[mm/min]400400200
Friction coefflcient[-]0.050.050.05
Holding force[kN]21.116.6-

FEM software with the explicit time-integration scheme was used. Parameters defined in numerical simulation are shown in Table 4.

Table 4

Parameters defined in FEM software

ParameterValueParameterValue
Mesh typeQuadrilateralLevel of refinement5
Element typeShellHardening curveKrupkowski
Friction coefflcient0.05Hardening RuleIsotropic
Yield function (TH330)Hill48Yield function (AA5352)Barlat2000
Number of integration points5Failure criteriaStrain-based FLD

The tools geometry together with material input data were imported into FEM software. Krupkowski hardening curve was used for both materials because it shows the best correlation with experimental curves. For each material, different yield function was used. Hill48 yield function best describes the behaviour of the conventional steel materials during plastic deformation. This material law is most frequently used for steel materials and therefore been used in this case for steel tinplate TH330. However, this yield function is not suitable for aluminium alloys [16]. Therefore, more sophisticated yield function e.g. Barlat2000 is needed [16]. This law requires determining both, the exponent m and eight parameters used in the model. Exponent m is defined on the basis crystalline structure of the material. For metals with so-called BCC unit cell like steel is m = 6, while for FCC metals like aluminium alloys is exponent m = 8. These eight coefficients are required for this function and they have been gained numerically, based on the mechanical tests referred to Table 5 [17]. For anisotropic yield models, the biaxial r value (rb) is required and obtained from disk compression tests.

Table 5

Summary of experimental data used in calibration of the Yld2000-2d model

TestFlow stress ratioNormal anisotropy – r
Uniaxial tensile test (7x@0, 15, 30, 45, 60, 75, 90)
Disk compression-
Hydraulic bulge test-

Krupkowski hardening curve is defined as [17]:

(3)σ=Kεp+εon

where εp is plastic strain [-], εo is offset strain [-], n is strain hardening exponent [-] a K is material constant [MPa].

For plane stress conditions Hill48 yield function may be written as [16]:

(4)σK2=σ122r01+r0σ1σ2+r0r901+r901+r0σ22

where σK is yield stress [MPa], σ1 is major principal stress [MPa], σ2 is minor principal stress [MPa] and r0, r45, r90 are Lankford’s coefficients [-].

Barlat2000 (Yield2000 2nd) yield function can be written as [18]:

(5)ϕ=S1S2M;=2S2+S1M+2S1+S2M
(6)ϕ=S1S2M+2S2+S1M+2S1+S2M=2σyM

The Siand Siare the principal values of stress tensors S′ and S′′, which can be expressed by:

(7)S~1,2=12s~xx+s~yy±s~xxs~yy2+4s~xy2

The Sijand Sijare the linear functions of the stress deviators.

(8)sxxsyysxy=L11L120L21L22000L66sxxsyysxy
(9)sxxsyysxy=L11L120L21L22000L66sxxsyysxy
(10)L11L12L21L22L66=2/3001/30001/3002/30001α1α2α7
(11)L11L12L21L22L66=2282014440444402822000009α3α4α5α6α8

The coefficients Lij and Lij are described by relationships of a set of eight coefficients αk as shown in (10) and (11). The coefficients α1α8 are reduced to one coefficient if an isotropic case is considered.

Even if the aim of this section is to failure prediction during expansion operation, other important values were analysed too. Overall, there are four parameters evaluated in this work:

  1. Ear profile after deep drawing,

  2. Wall thickness after deep drawing in 25 mm and 45 mm distance from cup bottom,

  3. Punch force,

  4. The position of crack.

3 Results and discussion

Since we used anisotropic models (anisotropy in 0, 45 and 90 according to rolling direction) results will be same in every quadrant (every 90) of the specimen. For this reason, all compared results will be only in the first quadrant (from 0 to 90). As the first, the earing profile after redrawing cycle was examined. Cup height was measured around the whole specimen as shown in Figure 5. The biggest impact in ear creation has material anisotropy. In the directions where the coefficient of normal anisotropy was higher, ear occurs on cups.

Figure 5 Ear measurement [19]
Figure 5

Ear measurement [19]

The following data was measured in experimental cups:

  1. Steel cups:

    • – Creation of 4 ears in 0, 90, 180 and 270 direction

    • – Average cup height was 52.61 mm

    • – Amplitude (max–min) was 1.14 mm.

  2. Aluminium alloy cups:

    • – Creation of 8 ears in 0,45, 90, 135, 180, 225, 270 and 315 direction

    • – Average cup height was 52.05 mm

    • – Amplitude (max–min) was 4.13 mm

    • – As a result of blankholder load acting in the directions 0 and 180 was material highly thinned, and secondary ears were created (Figure 6).

Figure 6 Ear creation on aluminium cup [19]
Figure 6

Ear creation on aluminium cup [19]

Results of the experiment and FEM simulation are compared in the Figure 79. The cup heights of studied materials are shown in Figure 7 At the left side is the comparison of aluminium alloy experimental cup height with FEM simulation. The main difference is in 0 direction. That is due to the above-mentioned blankholder force applied in the rolling direction. In the directions 0 and 180 regarding to rolling direction was material highly thinned and secondary ears were created. This phenomenon is closely linked to the so-called ironing effect, which is an

Figure 7 Ear profile of the cup from aluminium alloy AA5352 (left) and steel sheet TH330 (right)
Figure 7

Ear profile of the cup from aluminium alloy AA5352 (left) and steel sheet TH330 (right)

Figure 8 Wall thickness of AA 5352 (left) and TH330 steel (right) in 25 mm (top) and 45 mm (bottom) from cup base
Figure 8

Wall thickness of AA 5352 (left) and TH330 steel (right) in 25 mm (top) and 45 mm (bottom) from cup base

Figure 9 Punch force during expansion operation for AA5352 (left) and TH330 sheets (right)
Figure 9

Punch force during expansion operation for AA5352 (left) and TH330 sheets (right)

increase in wall thickness (Figure 8). For deep drawing was used the same tool, with the gap between the punch and the die of 0.605 mm for drawing and 0.63 mm for redrawing. For material with a higher thickness, which was an aluminium alloy, the effect is more significant. The average cup height of aluminium was 51.81 mmand steel 52.56 mm. The amplitude for steel was 1.29 mm and for aluminium 2.22 mm.

Another compared parameter was wall thickness after redrawing cycle in 25 mmand 45 mmdistance from the cup base. Measurement of experimental samples was carried out on the whole circumference at a distance of 30. Results of numerical simulation were compared with experimental results and are shown in Figure 8 Wall thickness 25 mm from cup base is same for both materials and shown a little deviations.

Next, the punch force during the expansion operation was evaluated. Experimental results are compared with numerical simulation in Figure 9. On the right are results for steel TH330 and on the left for aluminium alloy AA5352.

The failure occurred during the last operation. The distance from the base to the expansion punch at the moment of crack formation for each material was measured (Figure 10). In aluminium alloy, the crack appeared when the punch stroke reached 25.55 mm, corresponding to distance from expansion punch to the bottom of cup 24.27 mm. As regards the steel TH330, the crack appeared when punch stroke reached 35.33 mm, which is 14.46 mmfrom the base. When the failure occurs, cup height reached 49.79 mm for TH330 and 49.82 mm for AA5352 material. Figure 11 shows the position of failures.

Figure 10 Expansion punch stroke at the moment of failure occurred (a.), drawn cups with failure of TH330 (b.) [19] and AA5352 (c.) [19]
Figure 10

Expansion punch stroke at the moment of failure occurred (a.), drawn cups with failure of TH330 (b.) [19] and AA5352 (c.) [19]

Figure 11 Failure position for AA5352 (a) and TH330 (b) sheets
Figure 11

Failure position for AA5352 (a) and TH330 (b) sheets

As shown in Figure 11 numerical results more accurately predicted the position of the crack at the AA5352 aluminium alloy as for TH330 steel. This may result from the yield function. Aluminium alloys require a more sophisticated description of yield function such as Barlat2000. This anisotropic model gives a more precise prediction of material behaviour. For most conventional steels, the Hill48 yield law is sufficient, but for high-strength steels or aluminium alloys, this model is not appropriate. To describe material behaviour with satisfactory accuracy, Hill90, BBC, Vegter or other more sophisticated models should be used.

4 Conclusions

Simulations of the sheet metal parts forming processes using FEM are performed in the early design phases. Many commercial FE programs are today capable of reliably predicting possible crack and tear formation, wrinkle, earing from plastic anisotropy and many different sheet metal defects and instabilities that significantly affect the efficient production of these parts. Approximation of the simulation model is satisfactory to achieve good forming simulation results. However, for a more accurate calculation, the material model used and the hardening rule for describing sheet metal behaviour have great influence. This study has the main objective of predicting the failure point after deep-drawing, reverse redrawing and expansion operations. Two commonly used materials of TH330 steel and AA5352 aluminium alloy for packaging can production were studied. For each sheet metal, the material model for FEM analysis was created. Results of numerical simulation were compared with experimental results.

Prediction of earing after the reverse redrawing operation due to the plastic anisotropy of the material showed good correlation with experimental results. For TH330, the experimentally measured average cup height was 52.61 mm and the difference in max and min height was 1.14 mm, with numerical simulation predicting a difference in max and min height of 1.3 mm and an average cup height of 52.27 mm. For the aluminium alloy, an average cup height was 52.05 mmand a difference in max and min height was 4.13 mm, which in numerical simulation corresponds to an average cup height was 50.76 mm and a difference in max and min height was 2.22 mm. As a result, the selected material model for TH330 steel was more accurate in this case.

Experimentally measured values of the wall thickness correlated with the FEM simulation results more at a distance of 25 mm from the cup bottom than 45 mm from the bottom. This applied to both investigated materials. However, in the case of aluminium alloys, deviations were not so pronounced. The more sophisticated material model thus predicted more accurately the thickness of the cup wall. The failure position after the expansion operation for AA5352 material was more accurately predicted. Our findings can be useful for achieving a more accurate prediction of the forming processes.


Ján Slota: Tel.: +421 55 602 3545

Marek Šiser: Tel.: +421 55 602 3519

Ivan Gajdoš: Tel.: +421 55 602 3518


Acknowledgement

This paper is the result of the implementation of projects VEGA 1/0441/17 and APVV-14-0834.

References

[1] Šiser M., Slota J., Material Model of AW 5754 H11 Al Alloy for Numerical Simulation of Deep Drawing Process, Acta Mechanica Slovaca, 2016, 20 no.2, 32-4010.21496/ams.2016.012Search in Google Scholar

[2] Slota J., Šiser M., Advancedmaterial models for stamping of AW 5754 aluminum alloy, Strength of Materials, 2016, 48 no. 4, 487-49410.1007/s11223-016-9790-zSearch in Google Scholar

[3] Sumikawa S., Ishiwatari A., Hiramoto J., Urabe T., Improvement of springback prediction accuracy using material model considering elastoplastic anisotropy and Bauschinger effect, Journal of Materials Processing Technology, 2016, 230, 1–710.1016/j.jmatprotec.2015.11.004Search in Google Scholar

[4] Mulder J., Vegter H., et. al., Accurate determination of flow curves using the bulge test with optical measuring systems, Journal of Materials Processing Technology, 2015, 226, 169-18710.1016/j.jmatprotec.2015.06.034Search in Google Scholar

[5] Williams B. W., Boyle K.P., Characterization of anisotropic yield surfaces for titanium sheet using hydrostatic bulging with elliptical dies, International Journal of Mechanical Sciences, 2016, 114, 315-32910.1016/j.ijmecsci.2016.05.022Search in Google Scholar

[6] Reis L.C., Oliveira M.C., Santos A.D., Fernandes J.V., On the determination of the work hardening curve using the bulge test, International Journal of Mechanical Sciences, 2016, 105, 158-18110.1016/j.ijmecsci.2015.11.009Search in Google Scholar

[7] Chen K., Scales M., Kyriakides S., Corona E., Effects of anisotropy on material hardening and burst in the bulge test, International Journal of Solids and Structures, 2016, 82, 70-8410.1016/j.ijsolstr.2015.12.012Search in Google Scholar

[8] Liu K., Lang L., Cai G., Yang X., Guo C., Liu B., A novel approach to determine plastic hardening curves of AA7075 sheet utilizing hydraulic bulging test at elevated temperature, International Journal of Mechanical Sciences, 2015, 100, 328-33810.1016/j.ijmecsci.2015.07.002Search in Google Scholar

[9] Tian H., Kang D., A study on determining hardening curve for sheet metal, International Journal of Machine Tools & Manufacture, 2003, 43, 1253-125710.1016/S0890-6955(03)00132-9Search in Google Scholar

[10] Esmaeilpour R., et. al., Comparison of 3D yield functions for finite element simulation of single point incremental forming (SPIF) of aluminum 7075, International Journal of Mechanical Sciences, 2017, 133, 544-55410.1016/j.ijmecsci.2017.09.019Search in Google Scholar

[11] Cardoso R.P.R., Adetoro O.B., A generalisation of the Hill’s quadratic yield function for planar plastic anisotropy to consider loading direction, International Journal of Mechanical Sciences, 2017, 128-129, 253-26810.1016/j.ijmecsci.2017.04.024Search in Google Scholar

[12] Banabic, D., et al., FLD theoretical model using a new anisotropic yield criterion, Materials Processing Technology, 2004, 273-281.10.1016/j.jmatprotec.2004.09.015Search in Google Scholar

[13] Cao J., Li F., Ma X., Sun Z., A modified elliptical fracture criterion to predict fracture forming limit diagrams for sheet metals, Journal of Materials Processing Tech, 2017, 252, 116-12710.1016/j.jmatprotec.2017.09.018Search in Google Scholar

[14] Steel for packing: Product catalogue. Arcelor Mittal, Flat Carbon Europe, http://packaging.arcelormittal.com/repository2/Unassigned/ArcelorMittal%20Packaging%20%20product%20catalogue.pdfSearch in Google Scholar

[15] Jin, H., Wagstaff B., Gallerneault, M., Lloyd, D. J., Microstructure and Deformation Properties of AA5352/Al-6%Mg Sheet Produced by the Direct Chill Casting Fusion Process, Proceedings of the 12th International Conference on Aluminium Alloys, 2010, 2186-2191Search in Google Scholar

[16] Mulidrán, P.; Šiser, M.; Slota, J.; Spišák, E.; Sleziak, T., Numerical Prediction of Forming Car Body Parts with Emphasis on Springback, Metals, Vol. 8 (6), 435, 2018, 1-15, ISSN 2075-4701.10.3390/met8060435Search in Google Scholar

[17] Slota J., Numerical simulation of sheet metal forming processes, (in Slovak) 1st ed., Technical University, Košice, 2016Search in Google Scholar

[18] Panich, S., Barlat, F., Uthaisangsuk, V. et. al.: Experimental and theoretical formability analysis using strain and stress based forming limit diagram for advanced high strength steels, Materials & Design, Vol 51, 2013, 756-76610.1016/j.matdes.2013.04.080Search in Google Scholar

[19] Watson M., Dick R., Huang Y. H., Lockley A., Cardoso R., Santos A., Failure Prediction after Cup Drawing, Reverse Redrawing and Expansion Part, Journal of Physics: Conference Series 734, 2016, 1-8510.1088/1742-6596/734/2/022001Search in Google Scholar

Received: 2017-11-15
Accepted: 2018-03-08
Published Online: 2018-11-10

© 2018 J. Slota et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Downloaded on 29.3.2024 from https://www.degruyter.com/document/doi/10.1515/eng-2018-0045/html
Scroll to top button