Material parameter identification for the numerical simulation of deep-drawing drawing of aluminium alloys

Transcription

1 Material parameter identification for the numerical simulation of deep-drawing drawing of aluminium alloys B.M. Chaparro, J.M. Antunes, A.J. Baptista, D.M. Rodrigues,, J.V. Fernandes, L.F. Menezes M2D'26 - PORTO, JULY University of Coimbra CENTRO DE ENGENHARIA MECÂNICA DA UNIVERSIDADE DE COIMBRA I. Introduction II. III. IV. Constitutive models Inverse Analysis Demeri Test V. Conclusions OUTLOOK 2

2 Multiplicity of constitutive models F ( σ X L ) σ ( σ X M ) Y ( ) ',, = ',, = k k1 k 2 Equivalent stress associated with the yield criteria Equivalent stress associated with the work hardening law In this study, an automatic procedure for the identification of material parameters according to several phenomenological elastoplastic constitutive models was developed. Introduction 3 Hill, 1948 (Hill48) ( ) ( ) ( ) yy zz + zz xx + xx yy + 2 yz + 2 zx + 2 xy = F σ σ G σ σ H σ σ Lτ Mτ Nτ K Barlat, 1991 (YLD91) m m m = 2 S S S S S S K m Plastic surface s = L :σ (Isotropic Equivalent Stress) ( ) c2 + c3 / 3 c3 / 3 c2 / 3 c3 / 3 ( c3 + c1 ) / 3 c1 / 3 c2 / 3 c1 / 3 ( c1 + c2 ) / 3 L = c4 c5 c6 m=6 for BCC materials m=8 for FCC materials Constitutive Models - Yield Criteria 4

3 Isotropic hardening - Voce Law σ = Y + R sat (1 e p n v ε ) Y - Initial yield stress R sat - saturation stress n - constant isotropic hardening Kinematic hardening - Lemaitre & Chaboche model X& = Cx X sat ( σ X ) p σ X & ε C x and X sat are material parameters kinematic hardening Constitutive Models Work hardening law 5 The calculation methodology is based on the best fit between each model and the full set of available experimental results. 2 2 σ r erroryield _ criteria = w1 1 w exp exp σ r σ r exp error exp σ is the yield stress in tension for a specific orientation () with the rolling direction (RD) σ and r denotes the corresponding values obtained from the yield criteria wi exp is the r-ratio obtained in tension performed at a specific angle with the RD are weighting factors work hardening σ = 1 exp σ is the experimental equivalent stress 2 σ denotes the corresponding values obtained from the hardening law. Inverse Analysis 6

4 Materials: EN AW-5182-H111 EN AW-612-T4 Inverse analysis - Results 7 Experimental curves obtained in the tensile, shear and Baushingers shear tests for the EN AW-5182-H111 aluminium alloy Stress [MPa] σ11 / τ Strain ε 11 / γ _T_ 5182_T_ _T_3 5182_T_ _T_6 5182_T_ _T_9 5182_S_ 5182_S_ _S_3 5182_S_ _S_6 5182_S_ _S_9 5182_B_1 5182_B_2 5182_B_3 8

5 Experimental yield stresses and r values for the EN AW-5182-H111 aluminium alloy Angle with RD [º] r values Yield stress [Mpa] σ r erroryield _ criteria = w1 1 w exp exp σ r Hill48 F G N YLD91 c 1 c 2 c 3 c 6 m Material parameters identified for the Hill48 and YLD91 yield criteria 9 Comparison between experimental and numerical results Yield stresses Txy=. Txy=.25 Txy=.433 Txy=.5 EXP Txy=. Txy=.25 Txy=.433 Txy=.5 EXP Sy.6 Sy Sx Sx Hill48 YLD91 1

6 Comparison between experimental and numerical results r-values 1..9 r EXP HILL48 YLD Angle with RD 11 Isotropic work hardening parameters identified using the von Mises, Hill48 and YLD91 yield models error work hardening σ = 1 exp σ 2 Yield model R sat [MPa] nv Von Mises Hill YLD

7 Experimental and numerical equivalent stress-strain curves considering the von Mises yield criterion and pure isotropic work hardening _T_ 5182_T_ _T_3 5182_T_ _T_6 5182_T_ _T_9 5182_S_ 5182_S_ _S_3 5182_S_ _S_6 5182_S_ _S_9 DD3MAT Voce Iso 13 Experimental and numerical equivalent stress-strain curves considering the Hill 48 yield criterion and pure isotropic work hardening _T_ 5182_T_ _T_3 5182_T_ _T_6 5182_T_ _T_9 5182_S_ 5182_S_ _S_3 5182_S_ _S_6 5182_S_ _S_9 DD3MAT Voce Iso 14

8 Experimental and numerical equivalent stress-strain curves considering the YLD91 yield criterion and pure isotropic work hardening _T_ 5182_T_ _T_3 5182_T_ _T_6 5182_T_ _T_9 5182_S_ 5182_S_ _S_3 5182_S_ _S_6 5182_S_ _S_9 DD3MAT Voce Iso 15 Isotropic and Kinematic work hardening parameters identified using the von Mises, Hill48 and YLD91 yield models error work hardening σ = 1 exp σ 2 Yield model R sat [MPa] nv C x X sat [MPa] Mises Hill YLD

9 Experimental and numerical equivalent stress-strain curves considering the von Mises yield criterion and isotropic and kinematic work hardening _T_ 5182_T_ _T_3 5182_T_ _T_6 5182_T_ _T_9 5182_S_ 5182_S_ _S_3 5182_S_ _S_6 5182_S_ _S_9 5182_B_1 5182_B_2 5182_B_3 DD3MAT_B_1 DD3MAT_B_2 DD3MAT_B_3 17 Experimental and numerical equivalent stress-strain curves considering the Hill 48 yield criterion and isotropic and kinematic work hardening _T_ 5182_T_ _T_3 5182_T_ _T_6 5182_T_ _T_9 5182_S_ 5182_S_ _S_3 5182_S_ _S_6 5182_S_ _S_9 5182_B_1 5182_B_2 5182_B_3 DD3MAT_B_1 DD3MAT_B_2 DD3MAT_B_3 18

10 Experimental and numerical equivalent stress-strain curves considering the YLD91 yield criterion and isotropic and kinematic work hardening _T_ 5182_T_ _T_3 5182_T_ _T_6 5182_T_ _T_9 5182_S_ 5182_S_ _S_3 5182_S_ _S_6 5182_S_ _S_9 5182_B_1 5182_B_2 5182_B_3 DD3MAT_B_1 DD3MAT_B_2 DD3MAT_B_3 19 Experimental curves obtained in the tensile, shear and Baushingers shear tests for the EN AW-612-T4 aluminium alloy Results - EN AW- 612 T4 2

11 Experimental yield stresses and r values for the EN AW-612 T4 aluminium alloy Angle with RD [º] r value Yield stress [MPa] σ r erroryield _ criteria = w1 1 w exp exp σ r Hill48 F G N YLD91 c1 c2 c3 c6 m Material parameters identified for the Hill48 and YLD91 yield criteria Results - EN AW-612 T4 21 Comparison between experimental and numerical results Yield stresses Txy=. Txy=.25 Txy=.433 Txy=.5 EXP Txy=. Txy=.25 Txy=.433 Txy=.5 EXP Sy.6 Sy Sx Sx Hill48 YLD91 Results - EN AW-612 T4 22

12 Comparison between experimental and numerical results r-values 1..9 r EXP HILL48 YLD Angle with RD Results - EN AW-612 T4 23 Isotropic work hardening parameters identified using the von Mises, Hill48 and YLD91 yield models error work hardening σ = 1 exp σ 2 Yield model R sat [MPa] nv Mises Hill YLD Results - EN AW-612 T4 24

13 Experimental and numerical equivalent stress-strain curves considering the von Mises yield criterion and pure isotropic work hardening _T_ 612_T_15 612_T_3 612_T_45 612_T_6 612_T_75 612_T_9 612_S_ 612_S_15 612_S_3 612_S_45 612_S_6 612_S_75 612_S_9 DD3MAT Voce Iso Results - EN AW-612 T4 25 Experimental and numerical equivalent stress-strain curves considering the Hill 48 yield criterion and pure isotropic work hardening _T_ 612_T_15 612_T_3 612_T_45 612_T_6 612_T_75 612_T_9 612_S_ 612_S_15 612_S_3 612_S_45 612_S_6 612_S_75 612_S_9 DD3MAT Voce Iso Results - EN AW-612 T4 26

14 Experimental and numerical equivalent stress-strain curves considering the YLD91 yield criterion and pure isotropic work hardening _T_ 612_T_15 612_T_3 612_T_45 612_T_6 612_T_75 612_T_9 612_S_ 612_S_15 612_S_3 612_S_45 612_S_6 612_S_75 612_S_9 DD3MAT Voce Iso Results - EN AW-612 T4 27 Isotropic and Kinematic work hardening parameters identified using the von Mises, Hill48 and YLD91 yield models error work hardening σ = 1 exp σ 2 Yield model R sat [MPa] nv von Mises Hill YLD C x X sat [MPa] Results - EN AW-612 T4 28

15 Experimental and numerical equivalent stress-strain curves considering the von Mises yield criterion and isotropic and kinematic work hardening _T_ 612_T_15 612_T_3 612_T_45 612_T_6 612_T_75 612_T_9 612_S_ 612_S_15 612_S_3 612_S_45 612_S_6 612_S_75 612_S_9 612_B_1 612_B_2 612_B_3 DD3MAT_B_1 DD3MAT_B_2 DD3MAT_B_3 Results - EN AW-612 T4 29 Experimental and numerical equivalent stress-strain curves considering the Hill 48 yield criterion and isotropic and kinematic work hardening _T_ 612_T_15 612_T_3 612_T_45 612_T_6 612_T_75 612_T_9 612_S_ 612_S_15 612_S_3 612_S_45 612_S_6 612_S_75 612_S_9 612_B 1 612_B 2 612_B 3 DD3MAT_B_1 DD3MAT_B_2 DD3MAT_B_3 Results - EN AW-612 T4 3

16 Experimental and numerical equivalent stress-strain curves considering the YLD 91 yield criterion and isotropic and kinematic work hardening _T_ 612_T_15 612_T_3 612_T_45 612_T_6 612_T_75 612_T_9 612_S_ 612_S_15 612_S_3 612_S_45 612_S_6 612_S_75 612_S_9 612_B 1 612_B 2 612_B 3 DD3MAT_B_1 DD3MAT_B_2 DD3MAT_B_3 Results - EN AW-612 T4 31 Springback evaluation Demeri Test 1 25 a) Drawn cup b) Ring cutting c) Ring specimen d) Ring springback Difference between the ring diameters, before and after splitting, gives a direct measure of the springback effect Formability Analysis 32

17 Benchmark experimental apparatus Final stage of the stamping operation After the stamping operation Formability Analysis 33 NUMERICAL SIMULATIONS PROCEDURE DD3IMP DD3TRIM Blank discretization Full model used in the simulations Formability Analysis 34

18 Punch force during the stamping operations Force [kn] EXP EN AW-5182-H111 NUM EN AW-5182-H111 2 EXP EN AW-612-T4 1 NUM EN AW-612-T Displacement [mm] Formability Analysis 35 Springback results Difference between the ring diameters, before and after splitting Material Experimental ring opening Numerical ring opening EN AW-5182-H mm 87.4 mm EN AW-612-T mm 47.4 mm Formability Analysis 36

19 Independently of the constitutive model considered, the material parameter identification error reported was always less than 5%. For the EN AW-5182 alloy, the YLD91 criteria enable an accurate description of the plastic behaviour of this alloy for all the testing directions. For the EN AW-612 alloy, none of the plasticity criteria enables the correct description of the entire range of stress-strain curves. This can be related with the strong variation of the plastic anisotropy coefficients in the sheet plane. A more flexible criterion is needed to improve the description of the plastic behaviour of this material. The numerical program developed for the material parameter identification tasks proved to be an efficient tool that can be used to select the best phenomenological model. Conclusions 37