Next Article in Journal
Selective Laser Spectroscopy of the Bixbyite-Type Yttrium Scandate Doped by Rare-Earth Ions
Next Article in Special Issue
Relative Entropy Application to Study the Elastoplastic Behavior of S235JR Structural Steel
Previous Article in Journal
The Influence of Different Curing Environments on the Mechanical Properties and Reinforcement Mechanism of Dredger Fill Stabilized with Cement and Polypropylene Fibers
Previous Article in Special Issue
Effect of Austempering Processes on the Tensile Properties and the Work-Hardening Behavior of Austempered Bainitic Steels Below the Martensite Start Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plastic Evolution Characterization for 304 Stainless Steel by CQN_Chen Model under the Proportional Loading

1
College of Mechanical and Vehicle Engineering, Chongqing University, Chongqing 400044, China
2
China Automotive Engineering Research Institute Co., Ltd., Chongqing 401122, China
3
School of Mechanical Engineering, Xi’an Jiaotong University, 28 Xianning West Rd, Xi’an 710049, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(21), 6828; https://doi.org/10.3390/ma16216828
Submission received: 16 September 2023 / Revised: 16 October 2023 / Accepted: 18 October 2023 / Published: 24 October 2023

Abstract

:
In this paper, the CQN_Chen function is used to characterize the plastic anisotropic evolution of 304 stainless steel (SS304). The uniaxial tensile tests along different loading directions are conducted to experimentally investigate the anisotropic hardening behavior for SS304. The experimental data indicates that the anisotropy of SS304 is weak. The convexity analysis is carried out by the geometry-inspired numerical convex analysis method for the CQN_Chen yield locus during plastic deformation. The Hill48, SY2009 and CQN functions are used as the comparison to evaluate the accuracy of the CQN_Chen function in characterizing plastic evolution. The predicted values are compared with the experimental data. The comparison demonstrates that the CQN_Chen function can accurately characterize anisotropic hardening behavior under uniaxial tension along distinct loading directions and equibiaxial tension. Simultaneously, the CQN_Chen model has the capacity to adjust the yield surface shape between uniaxial tension and equibiaxial tension. The CQN_Chen model is recommended to characterize plastic evolving behavior under uniaxial tension along different directions and equibiaxial tension.

1. Introduction

Due to carbon neutrality, lightweight materials are widely used to reduce weight, such as aluminum alloy, magnesium alloy and stainless steel [1]. During the rolling forming process, the textures of these sheet metals form preferred orientation, which leads to anisotropy. It is vital to accurately characterize the plastic anisotropic evolution in sheet-metal forming numerical simulation.
The Hill48 function was the pioneer of anisotropic yield functions [2]. However, the Hill48 function cannot distinguish the difference of yield surface of metals with different crystal structures because of its quadratic form. Hosford [3,4] combined the exponent with the yield surface shapes for body-centered cubic (BCC) and face-centered cubic (FCC) metals. Barlat et al. [5] extended the isotropic function to anisotropy by introducing anisotropic coefficients through a linear transformed tensor. More anisotropic coefficients were introduced through two linear transformed tensors to characterize strongly anisotropic metals and consider r-values, such as Yld2000-2d [6] and Yld2004-18p [7]. Cazacu et al. [8] introduced an orthotropic yield criterion in the form of the principal values of the stress deviator to capture the anisotropy and the asymmetry between tension and compression. He et al. [9] proposed an enhanced constitutive model based on the Yld2000-2d model. Recently, many yield functions were developed based on stress invariants. Cazacu and Barlat [10] extended the isotropic Drucker [11] function to orthotropy. Cazacu and Barlat [12] modelled the asymmetry and anisotropy for pressure insensitive metals. Gao et al. [13] developed a plasticity model with the form of the hydrostatic stress as well as the second and third invariants of the stress deviator. Aretz and Barlat [14] developed a very flexible anisotropic yield function with 27 anisotropic parameters. Yoshida et al. [15] developed an anisotropic yield function in a form of the stress invariants. Hu et al. [16] proposed a normalized stress invariant-based yield criterion to describe both asymmetry and Lankford coefficients. Lou et al. [17] extended the Drucker yield function to anisotropy by linear transformed tensor, and calibrated the yield surface shape control parameters for BCC and FCC materials. When a linear transformed tensor is used to introduce anisotropic coefficients, it is generally necessary to use an optimization algorithm to optimize error functions to calibrate anisotropic coefficients under specific plastic work. Xu et al. [18] reduced the Yld2011-27p function to decrease the number of experiments for the calibration of anisotropic parameters. Therefore, the anisotropic yield functions commonly cannot accurately describe the change of the yield surface with increasing plastic strain under the isotropic hardening.
The anisotropic hardening model captures the evolution of the yield surface by analytically describing the anisotropic coefficient. Namely, the hardening curves under different stress states are directly used to calculate the anisotropic coefficients with the increase of the plastic deformation. The anisotropic hardening was experimentally observed by Choi et al. [19], Khan et al. [20,21,22], Lou et al. [23], etc. Stoughton and Yoon [24] replaced the anisotropic coefficients of the Hill48 function with hardening curves under uniaxial tension along the rolling direction (RD), diagonal direction (DD) and transverse direction (TD) as well as equibiaxial tension, which achieved the accurate description of anisotropic hardening behavior (SY2009). The SY2009 function inherited the characteristic of the Hill48 function. Namely, the shape of the yield surface cannot be adjusted by the SY2009 model. Lee et al. [25] multiplicatively coupled SY2009 with the Hosford72 [26] yield criteria, in which the plastic evolution was characterized by SY2009 function, and the shape of the yield surface was adjusted by the Hosford72 function (CQN). Chen et al. [27] coupled SY2009 with Drucker function to characterize anisotropic hardening for BCC and FCC metals (CQN_Chen). Hou et al. [28] coupled a quadratic asymmetric yield function with a non-quadratic function by multiplication to characterize the anisotropic–asymmetric hardening behavior and adjust the yield surface shape. Hou et al. [29] proposed an anisotropic hardening model by coupling an asymmetric Hill48 function with an isotropic stress-invariant-based yield function. Hu et al. [30] coupled a fourth order polynomial yield criterion with a non-quadratic yield function under associated flow rule to describe the evolution of anisotropic yielding behavior analytically. Zhang and Lou [31] characterized the evolving plastic behavior by coupling an enhanced pDrucker function and SY2009 model. Zhou et al. [32] proposed a new analytical asymmetric yield criterion base on the SY2009. In addition to constructing anisotropic hardening model through multiplicative coupling method, a direct analytical description of the yield function is developed by the hardening curves or r-value evolution of different stress states to describe the anisotropic hardening. Hu et al. [33] analytically described the Yoon2014 [34] function by considering anisotropic hardening. Hu et al. [35] analytically described polynomial yield criterion by considering both plane strain and pure shear states. Lou et al. [36] converted the stress-invariant-based function into stress triaxiality, Lode parameter and von Mises stress forms, and analytically described the parameters through four of the five stress states. Hou et al. [37] proposed a constitutive model to accurately describe the anisotropic behaviors of sheet metals in terms of the yield stress and plastic flow under plane strain loading. Lou and Yoon [38] constructed an anisotropic-asymmetric hardening model by the additive coupling of two Hill48 functions with Lode-dependent weight functions. The more parameters in the analytical description of the yield function, the more anisotropic hardening behaviors of the stress states can be characterized.
The anisotropic models are widely used to predict the formation of sheet metals. Clausmeyer and Svendsen [39] compared two of the models for anisotropic hardening and yield surface evolution in BCC sheet steels. Min et al. [40] used the Yld89 [41] function to consider anisotropic hardening and sheet anisotropy for a multi-phase 980 MPa steel. Li et al. [42] characterized the anisotropic hardening behavior of metals using the interpolation method. Hao and Dong [43] proposed an interpolation-based anisotropic yield and hardening models. Fu et al. [44] proposed a method to simultaneously identify the anisotropic yield and hardening parameters from a single test. Li et al. [45] developed an improved yield criterion to characterize the anisotropic and tension-compression asymmetric behavior of magnesium alloy. Lou et al. [46] developed a reduced form to the Yld2004-18p function to decrease the experimental cost for the calibration of anisotropic parameters. Wang et al. [47] investigated the earing characteristics of the 6K21 aluminum alloy in a circular-cup deep drawing by virtue of experiments and multi-scale simulations. Hou et al. [48] described the evolving yield surfaces of dual-phase steels. Mu et al. [49] characterized the anisotropic hardening and evolution of r-values for sheet metal based on the evolving non-associated Hill48 model. Du et al. [50] characterized the asymmetric evolving yield and flow of the 6016-T4 aluminum alloy and DP490 steel by the several existing asymmetric yield criteria under the associated and non-associated flow rules. Mamros et al. [51] captured the plastic anisotropy evolution of stainless steel 316L under proportional loadings. Yang et al. [52] investigated the anisotropic plastic flow of low/medium carbon steel plates in different loading conditions. Lee et al. [53] utilized the evolving Hill48 function to describe the evolution of the distortional yield surface of the Al6014-T4 alloy. Gawad et al. [54] developed an evolving plane stress yield criterion based on crystal plasticity virtual experiments.
This research investigates the plastic anisotropic evolution of 304 stainless steel (SS304) by using the CQN_Chen function. The uniaxial tensile tests along RD, DD and TD are conducted. The convex domain of CQN_Chen function is determined by the geometry-inspired numerical convex analysis (GINCA) method with the increasing plastic evolution. The Hill48, SY2009 and CQN models are selected as a comparison to reveal the capability of CQN_Chen model in characterizing anisotropic hardening and adjusting yield surface shape. The experimental yield stresses are fully compared to the predicted values by the four models with the increase of the plastic deformation, including 3D plane stress yield surfaces σ x x , σ y y , σ x y , 2D yield locus under biaxial loading, predicted uniaxial tensile yield stresses, and hardening curves under uniaxial tension along RD, DD and TD as well as equibiaxial tension, which is used to evaluate the accuracy of CQN_Chen function in describing plastic evolution.

2. CQN_Chen Anisotropic Hardening Model

2.1. CQN_Chen Anisotropic Hardening Model

The CQN_Chen function is proposed by Chen et al. [27], as shown in Equation (1). The CQN_Chen is formulated by coupling SY2009 and Drucker function. Hence, the function can accurately capture anisotropic hardening and adjust the yield surface shape for BCC and FCC metals under the proportional loading conditions. Meanwhile, the function is convenient for numerical simulation and has high simulation efficiency.
f C Q N _ C h e n ( σ , λ ¯ ) = [ f H i l l 48 ( σ , λ ¯ ) f D r u c ker ( σ ) ] 1 8
with
f H i l l 48 ( σ , λ ¯ ) = F ( λ ¯ ) ( σ 22 σ 33 ) 2 + G ( λ ¯ ) ( σ 33 σ 11 ) 2 + H ( λ ¯ ) ( σ 11 σ 22 ) 2 + 2 L ( λ ¯ ) σ 23 2 + 2 M ( λ ¯ ) σ 31 2 + 2 N ( λ ¯ ) σ 12 2
f D r u c ker ( σ ) = a ( J 2 3 c J 3 2 )
where λ ¯ denotes the plastic compliance factor; ( σ 11 , σ 22 , σ 33 , σ 12 , σ 23 , σ 13 ) present Cauchy stress components; J 2 and J 3 indicate the second and third invariants of stress deviator, respectively; c controls the yield surface shape, and c = 1.5776 for BCC and c = 2.5116 for FCC; a = 729 27 4 c ; and F ( λ ¯ ) , G ( λ ¯ ) , H ( λ ¯ ) , L ( λ ¯ ) , M ( λ ¯ ) and N ( λ ¯ ) are anisotropic parameters of the Hill48 function. The six parameters are functions of the hardening curves along different loading directions, which is depicted as follows:
F ( λ ¯ ) = 1 2 ( 1 [ σ 90 ( λ ¯ ) ] 8 + 1 [ σ b ( λ ¯ ) ] 8 1 [ σ 0 ( λ ¯ ) ] 8 )
G ( λ ¯ ) = 1 2 ( 1 [ σ 0 ( λ ¯ ) ] 8 + 1 [ σ b ( λ ¯ ) ] 8 1 [ σ 90 ( λ ¯ ) ] 8 )
H ( λ ¯ ) = 1 2 ( 1 [ σ 0 ( λ ¯ ) ] 8 + 1 [ σ 90 ( λ ¯ ) ] 8 1 [ σ b ( λ ¯ ) ] 8 )
L ( λ ¯ ) = M ( λ ¯ ) = N ( λ ¯ ) = 1 2 ( 4 [ σ 45 ( λ ¯ ) ] 8 1 [ σ b ( λ ¯ ) ] 8 )
where σ 0 , σ 45 , σ 90 and σ b are the yield stresses under uniaxial tension along the RD, DD and TD as well as equibiaxial tension, respectively.

2.2. Convexity Analysis

Equations (2) and (3) observe that the c and λ ¯ values are related to the convexity of the CQN_Chen yield surface. Namely, the c adjusts the yield surface curvature, and the flow stresses under distinct plastic strain are obtained by the λ ¯ to compute the anisotropic parameters of Hill48 function. The convexity of the yield locus is usually validated by the Hessian matrix with the positive semi-definite. The Hessian matrix is obtained by calculating the second order partial derivative of the yield function, and the convexity of the yield surface is judged according to the eigenvalue of the Hessian matrix. It is inconvenient and time-consuming to verify the convexity of the yield locus during the plastic deformation process. Lou et al. [36] proposed a geometry-inspired numerical convex analysis (GINCA) method to determine the convexity of the yield surface by calculating the equivalent stress, as depicted in Figure 1. The points A, B and C are on the 3D plane stress yield surface. Point D is located at the center of the infinitesimal curved surface A B C ˜ . The coordinates of points A, B and C are presented as ( x A , y A , z A ) , ( x B , y B , z B ) and ( x C , y C , z C ) , respectively. Therefore, the coordinate of point D is denoted as ( ( x A + x B + x C ) / 3 , ( y A + y B + y C ) / 3 , ( z A + z B + z C ) / 3 ) . Due to points A, B and C being located on the yield surface, the equivalent stresses at these points are f A = f B = f C = 1 . The equivalent stress at point D needs to be calculated to judge the convexity of the A B C ˜ . If point D is located inside the curved surface A B C ˜ , then A B C ˜ is convex. Namely, the curved surface A B C ˜ is convex if f D 1 . On the contrary, the curved surface A B C ˜ is concave if f D > 1 . It is convenient to use the GINCA method to analyze the convexity of the yield surface during the plastic evolving process. It is only necessary to calculate the equivalent stress of the yield function without the computation of the second–order partial derivative compared to the complicated Hessian matrix.

3. Plastic Evolution Characterization for SS304

3.1. Uniaxial Tensile Experiment

The CQN_Chen function is analytically described by the hardening curves of uniaxial tension along the RD, DD and TD as well as equibiaxial tension for SS304 in Equations (4)–(7). Hence, the dog–bone specimens with a thickness of 0.4 mm are cut by a laser cutting machine along the RD, DD and TD for uniaxial tension. The dimensions of the dog–bone specimen are shown in Figure 2. The uniaxial tensile tests are under quasi–static loading with a tensile speed of 3.6 mm/min. The uniaxial tensile test is repeated five times in each direction to ensure the reliability of the experimental results. The surface of the specimen is sprayed with speckle uniformly, and the deformation process is recorded by the XTOP digital image correction system. The virtual extensometer is set to 30 mm to obtain the load–stroke curve in the longitudinal direction, as presented in Figure 2a.
Figure 2a shows that the load–stroke curves in the RD and DD are almost identical, and the load–stroke curve in the TD is slightly higher than that in the RD and DD. The engineering stress–engineering strain curve is depicted in Figure 2b. According to the principle of plastic work equivalence, the uniaxial tensile hardening curves along different directions are converted into functions of λ ¯ in order to calculate the parameters of the CQN_Chen function, as shown in Figure 3. The hardening curve of TD is slightly higher than that of RD and DD at 0 . 15 λ ¯ 0 . 35 . The anisotropy is not obvious for SS304. The equivalent strain field in Figure 3 shows that the equivalent plastic strain of SS304 at fracture is about 0.6, which indicates that SS304 has good plastic deformation capacity. The uniaxial tensile hardening curves along the RD, DD and TD are fitted by the Hockett–Sherby/Hollomon hardening model. The parameters of the hardening law are listed in Table 1. The yield stress under equibiaxial tension is assumed as σ b = ( σ 0 + 2 σ 45 + σ 90 ) / 4 .

3.2. Anisotropic Evolution Characterization of SS304

The convex domain of the CQN_Chen yield surface is investigated by using the GINCA method under 3D–plane stress space with plastic evolution. Under the assumption of isotropy, the CQN_Chen function is reduced to a Cazacu2018 function [55]. Hence, the convex domain of the CQN_Chen function should be within the convex domain of the Cazacu2018 function. The comparison between the convex domain of CQN_Chen and the convex domain −5.4 c 3 of the Cazacu2018 function is shown in Figure 4. The blue solid line and the red solid line in Figure 4 represent the convex domains of the Cazacu2018 and CQN_Chen functions, respectively. The comparison indicates that the upper and lower limits of parameter c for CQN_Chen are located in the convex domain of Cazacu2018 function. The parameter c of the CQN_Chen function is recommended to be 1.5776 for SS304 of BCC metal. c = 1.5776 is always in the convex domain of the parameter c at 0 λ ¯ 10 . Therefore, the CQN_Chen yield surface remains convex during plastic deformation for SS304.
Figure 5 presents the 3D plane stress yield surfaces predicted by the CQN_Chen function at λ ¯ = 0 , 0.1, 0.15 and 0.2. The yield locus evolves with an increasing λ ¯ value. The experimental yield stresses are precisely located on the yield surfaces for uniaxial tension along the RD, DD and TD as well as equibiaxial tension at different λ ¯ . This is because these yield stresses are used to analytically describe the parameters of the CQN_Chen function. The anisotropic hardening behavior for SS304 is accurately characterized by the CQN_Chen model.
In order to describe the comparison more intuitively between the predicted yield stresses and the tested values, the 2D yield surface under biaxial loading space is extracted, which is denoted by the solid red line highlighted by <Biaxial loading along RD and TD> in Figure 5. The Hill48 anisotropic function, SY2009 and CQN anisotropic hardening models are selected as a comparison to evaluate the ability of the CQN_Chen model to characterize anisotropic hardening and adjust yield surface curvature. Figure 6 shows the comparison of 2D yield surfaces predicted by different models with experimental data under biaxial loading space. For the Hill48 function with the assumption of isotropic hardening, only the 2D yield surface at λ ¯ = 0 accurately captures the experimental yield stresses of equibiaxial tension and uniaxial tension along the RD, DD and TD. The uniaxial tensile yield stress along the RD is accurately predicted as the plastic strain increases. However, the predicted yield stresses are underestimated compared to the experimental values under unaixial tension along the TD, and to equibiaxial tension with the increase in plastic deformation. For the three anisotropic hardening models, the predicted yield surfaces at different plastic strain all precisely pass the experimental yield stresses for uniaxial tension along RD and TD as well as equibiaxial tension. The predicted yield surfaces of the CQN and CQN_Chen functions are almost identical. Compared to the predicted yield surface of CQN and CQN_Chen, the predicted yield stresses of SY2009 are significantly smaller than that of CQN and CQN_Chen under plane strain tension along the RD and TD as well as shear along DD. This is because the SY2009 function with the quadratic form cannot distinguish the yield surface differences of metals with different crystal structures. Figure 7 presents 2D yield surface evolution predicted by the CQN_Chen function at 0 λ ¯ 0 . 6 with an interval of 0.02. The predicted yield surfaces by CQN_Chen are in good agreement with the experimental data of uniaxial tension along the RD and TD as well as the equibiaxial tension at different λ ¯ values. The CQN_Chen function accurately modeled anisotropic hardening behavior and adjusted the shape of the yield surface for SS304.
The solid pink lines in Figure 5 are extracted and represents the predicted uniaxial tensile yield stresses along the different directions. Figure 8 indicates the comparison of experimental data with the uniaxial tensile yield stresses along distinct loading directions predicted by the Hill48, SY2009, CQN and CQN_Chen models. The uniaxial tensile yield stresses predicted by the four models accurately pass the experimental values of uniaxial tension along the RD, DD and TD at λ ¯ = 0 . With the increase in plastic strain, the uniaxial tensile yield stresses along distinct loading directions predicted by the Hill48 function are significantly lower than the tested values. However, the uniaxial tensile yield stress along the RD can be accurately predicted by the Hill48 function. The predicted uniaxial tensile yield stresses along different directions nearly overlap for the three anisotropic hardening models. These three models accurately capture the experimental uniaxial tensile yield stresses along RD, DD and TD. The uniaxial tensile yield stresses along RD, DD and TD are precisely characterized at the different plastic strain by CQN_Chen for SS304.
Figure 9 illustrates the comparison of predicted and tested values for the hardening curves under uniaxial tension along the RD, DD and TD as well as equibiaxial tension. The comparison demonstrates that the four hardening curves are accurately predicted by the SY2009, CQN and CQN_Chen functions. The Hill48 function can only accurately characterize the uniaxial tensile hardening curve along the RD under isotropic hardening. Equation (8) is used to quantitatively evaluate the prediction accuracy of different models for the four hardening curves, as depicted in Figure 10. Consistent with the above analysis, the predicted errors of four hardening curves are zero for the SY2009, CQN and CQN_Chen models. The maximum predicted error of the Hill48 function is about 0.035 for uniaxial tension along TD. Meanwhile, the predicted errors of the Hill48 function are less than zero for the hardening curves of uniaxial tension along the DD and TD as well as equibiaxial tension, which indicates that the predicted yield stresses by the Hill48 function is underestimated under these stress states. The CQN_Chen function accurately describes the hardening behavior for equibiaxial tension and uniaxial tension along the RD, DD and TD.
E r r o r = ( σ p r e d σ exp ) σ exp
where σ p r e d and σ exp present the yield stresses of the experiment and prediction, respectively.

4. Conclusions

The uniaxial tensile tests along distinct loading directions are conducted to investigate the anisotropic hardening behaviors for SS304. The CQN_Chen function is used to characterize the plastic evolution under uniaxial tension along the RD, DD and TD as well as equibiaxial tension for SS304. The convex domain is verified by GINCA method for the CQN_Chen yield surface. The Hill48, SY2009 and CQN models are selected as the comparison to evaluate the accuracy of CQN_Chen function in describing anisotropic hardening. The predicted values of the four models are compared with the experimental data, including the 3D plane stress yield surface, 2D yield locus under biaxial loading space, the predicted uniaxial tensile yield stresses along different loading directions and the predicted hardening curves. The results shows that the anisotropy of SS304 is not obvious. The CQN_Chen function can accurately model the anisotropic evolution under uniaxial tension along the RD, DD and TD as well as equibiaxial tension with increasing plastic strain. Meanwhile, the yield surface difference are distinguished for BCC and FCC metals. The CQN_Chen model is recommended to characterize plastic evolving behavior for BCC and FCC metals under proportional loading condition.

Author Contributions

Conceptualization, X.G. and Y.L.; Methodology, X.G., S.W. and Y.L.; Software, X.G. and S.W.; Validation, S.W.; Formal analysis, S.W.; Investigation, S.W., Z.X., J.Z., X.W. and H.M.A.R.; Resources, X.G. and Y.L.; Writing—original draft, S.W.; Writing—review & editing, Y.L.; Supervision, Y.L.; Project administration, Y.L.; Funding acquisition, X.G. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 52075423 and U2141214) and the Fundamental Research Funds for the Central Universities (Grant No. xtr012019004 and zrzd2017027).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hou, Y.; Myung, D.; Park, J.K.; Min, J.Y.; Lee, H.R.; El–Aty, A.A.; Lee, M.G. A Review of Characterization and Modelling Approaches for Sheet Metal Forming of Lightweight Metallic Materials. Materials 2023, 16, 836. [Google Scholar] [CrossRef] [PubMed]
  2. Hill, R. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. A–Math Phys. Eng. Sci. 1948, 193, 281–297. [Google Scholar]
  3. Logan, R.W.; Hosford, W.F. Upper–bound anisotropic yield locus calculations assuming <111>–pencil glide. Int. J. Mech. Sci. 1980, 22, 419–430. [Google Scholar]
  4. Hosford, W.F. Comments on anisotropic yield criteria. Int. J. Mech. Sci. 1985, 27, 423–427. [Google Scholar] [CrossRef]
  5. Barlat, F.; Lege, D.J.; Brem, J.C. A six–component yield function for anisotropic materials. Int. J. Plast. 1991, 7, 693–712. [Google Scholar] [CrossRef]
  6. Barlat, F.; Brem, J.C.; Yoon, J.W.; Chung, K.; Dick, R.E.; Lege, D.J.; Pourboghrat, F.; Choi, S.H.; Chu, E. Plane stress yield function for aluminum alloy sheets–part 1: Theory. Int. J. Plast. 2003, 19, 1297–1319. [Google Scholar] [CrossRef]
  7. Barlat, F.; Aretz, H.; Yoon, J.W.; Karabin, M.E.; Brem, J.C.; Dick, R.E. Linear transfomation–based anisotropic yield functions. Int. J. Plast. 2005, 21, 1009–1039. [Google Scholar] [CrossRef]
  8. Cazacu, O.; Plunkett, B.; Barlat, F. Orthotropic yield criterion for hexagonal closed packed metals. Int. J. Plast. 2006, 22, 1171–1194. [Google Scholar] [CrossRef]
  9. He, Z.B.; Zhang, K.; Zhu, H.H.; Lin, Y.L.; Fu, M.W.; Yuan, S.J. An anisotropic constitutive model for forming of aluminum tubes under both biaxial tension and pure shear stress states. Int. J. Plast. 2022, 152, 103259. [Google Scholar] [CrossRef]
  10. Cazacu, O.; Barlat, F. Generalization of Drucker’s Yield Criterion to Orthotropy. Math. Mech. Solids 2001, 6, 613–630. [Google Scholar] [CrossRef]
  11. Drucker, D.C. Relation of Experiments to Mathematical Theories of Plasticity. J. Appl. Mech. 1949, 16, 349–357. [Google Scholar] [CrossRef]
  12. Cazacu, O.; Barlat, F. A criterion for description of anisotropy and yield differential effects in pressure–insensitive metals. Int. J. Plast. 2004, 20, 2027–2045. [Google Scholar] [CrossRef]
  13. Gao, X.S.; Zhang, T.T.; Zhou, J.; Graham, S.M.; Hayden, M.; Roe, C. On stress–state dependent plasticity modeling: Significance of the hydrostatic stress, the third invariant of stress deviator and the non-associated flow rule. Int. J. Plast. 2011, 27, 217–231. [Google Scholar] [CrossRef]
  14. Aretz, H.; Barlat, F. New convex yield functions for orthotropic metal plasticity. Int. J. Non–Linear Mech. 2013, 51, 97–111. [Google Scholar] [CrossRef]
  15. Yoshida, F.; Hamasaki, H.; Uemori, T. A user–friendly 3D yield function to describe anisotropy of steel sheet. Int. J. Plast. 2013, 45, 119–139. [Google Scholar] [CrossRef]
  16. Hu, Q.; Li, X.F.; Han, X.H.; Li, H.; Chen, J. A normalized stress invariant–based yield criterion: Modeling and validation. Int. J. Plast. 2017, 99, 248–273. [Google Scholar] [CrossRef]
  17. Lou, Y.S.; Yoon, J.W. Anisotropic yield function based on stress invariants for BCC and FCC metals and its extension to ductile fracture criterion. Int. J. Plast. 2018, 101, 125–155. [Google Scholar] [CrossRef]
  18. Xu, T.Y.; Li, F.G.; Wang, X.L. A user–friendly anisotropic yield function for modeling anisotropy of BCC and FCC sheet metals. J. Mater. Eng. Perform. 2023, 32, 2370–2391. [Google Scholar] [CrossRef]
  19. Choi, Y.; Walter, M.E.; Lee, J.K.; Han, C.-S. Observations of anisotropy evolution and identification of plastic spin parameters by uniaxial tensile tests. J. Mech. Mater. Struct. 2006, 1, 301–323. [Google Scholar] [CrossRef]
  20. Khan, A.S.; Suh, Y.S.; Kazmi, R. Quasi–static and dynamic loading responses and constitutive modeling of titanium alloys. Int. J. Plast. 2004, 20, 2233–2248. [Google Scholar] [CrossRef]
  21. Khan, A.S.; Kazmi, R.; Farrokh, B. Multiaxial and non–proportional loading responses, anisotropy and modeling of Ti-6Al-4V titanium alloy over wide ranges of strain rates and temperatures. Int. J. Plast. 2007, 23, 931–950. [Google Scholar] [CrossRef]
  22. Khan, A.S.; Pandey, A.; Gnaupel-Herold, T.; Mishra, R.K. Mechanical response and Texture evolution of AZ31 alloy at large strains for different strain rates and temperatures. Int. J. Plast. 2011, 27, 688–706. [Google Scholar] [CrossRef]
  23. Lou, X.Y.; Li, M.; Boger, R.K.; Agnew, S.R.; Wagoner, R.H. Hardening evolution of AZ31B Mg sheet. Int. J. Plast. 2007, 23, 44–86. [Google Scholar] [CrossRef]
  24. Stoughton, T.B.; Yoon, J.W. Anisotropic hardening and non–associated flow in proportional loading of sheet metals. Int. J. Plast. 2009, 25, 1777–1817. [Google Scholar] [CrossRef]
  25. Lee, E.H.; Stoughton, T.B.; Yoon, J.W. A yield criterion through coupling of quadratic and non–quadratic functions for anisotropic hardening with non–associated flow rule. Int. J. Plast. 2017, 99, 120–143. [Google Scholar] [CrossRef]
  26. Hosford, W.F. A generalized isotropic yield criterion. J. Appl. Mech. 1972, 39, 607–609. [Google Scholar] [CrossRef]
  27. Chen, Z.; Wang, Y.; Lou, Y. User-friendly anisotropic hardening function with non–associated flow rule under the proportional loadings for BCC and FCC metals. Mech. Mater. 2022, 165, 104190. [Google Scholar] [CrossRef]
  28. Hou, Y.; Min, J.Y.; Stoughton, T.B.; Lin, J.P.; Carsley, J.E.; Carlson, B.E. A non-quadratic pressure-sensitive constitutive model under non-associated flow rule with anisotropic hardening: Modeling and validation. Int. J. Plast. 2020, 135, 102808. [Google Scholar] [CrossRef]
  29. Hou, Y.; Min, J.Y.; Lin, J.P.; Lee, M.G. Modeling stress anisotropy, strength differential, and anisotropic hardening by coupling quadratic and stress-invariant-based yield functions under non–associated flow rule. Mech. Mater. 2022, 174, 104458. [Google Scholar] [CrossRef]
  30. Hu, Q.; Yoon, J.W.; Manopulo, N.; Hora, P. A coupled yield criterion for anisotropic hardening with analytical description under associated flow rule: Modeling and validation. Int. J. Plast. 2021, 36, 102882. [Google Scholar] [CrossRef]
  31. Zhang, C.; Lou, Y. Characterization and modelling of evolving plasticity behaviour up to fracture for FCC and BCC metals. J. Mater. Process. Technol. 2023, 317, 117997. [Google Scholar] [CrossRef]
  32. Zhou, Y.; Hu, Q.; Chen, J. A concise analytical framework for describing asymmetric yield behavior based on the concept of shape functions. Int. J. Plast. 2023, 164, 103593. [Google Scholar] [CrossRef]
  33. Hu, Q.; Yoon, J.W. Analytical description of an asymmetric yield function (Yoon2014) by considering anisotropic hardening under non–associated flow rule. Int. J. Plast. 2021, 140, 102978. [Google Scholar] [CrossRef]
  34. Yoon, J.W.; Lou, Y.S.; Yoon, J.; Glazoff, M.V. Asymmetric yield function based on the stress invariants for pressure sensitive metals. Int. J. Plast. 2014, 56, 184–202. [Google Scholar] [CrossRef]
  35. Hu, Q.; Yoon, J.W.; Chen, J. Analytically described polynomial yield criterion by considering both plane strain and pure shear states. Int. J. Plast. 2023, 162, 103514. [Google Scholar] [CrossRef]
  36. Lou, Y.S.; Zhang, C.; Zhang, S.J.; Yoon, J.W. A general yield function with differential and anisotropic hardening for strength modelling under various stress states with non–associated flow rule. Int. J. Plast. 2022, 158, 103414. [Google Scholar] [CrossRef]
  37. Hou, Y.; Du, K.; El-Aty, A.A.; Lee, M.G.; Min, J.Y. Plastic anisotropy of sheet metals under plane strain loading: A novel non-associated constitutive model based on fourth–order polynomial functions. Mater. Des. 2022, 223, 111187. [Google Scholar] [CrossRef]
  38. Lou, Y.S.; Yoon, J.W. Lode-dependent anisotropic-asymmetric yield function for isotropic and anisotropic hardening of pressure–insensitive materials. Part I: Quadratic function under non–associated flow rule. Int. J. Plast. 2023, 166, 103647. [Google Scholar] [CrossRef]
  39. Clausmeyer, T.; Svendsen, B. Comparison of two models for anisotropic hardening and yield surface evolution in bcc sheet steels. Eur. J. Mech. A/Solids 2015, 54, 120–131. [Google Scholar] [CrossRef]
  40. Min, J.Y.; Carsley, J.E.; Lin, J.P.; Wen, Y.Y.; Kuhlenkötter, B. A non–quadratic constitutive model under non–associated flow rule of sheet metals with anisotropic hardening: Modeling and experimental validation. Int. J. Mech. Sci. 2016, 119, 343–359. [Google Scholar] [CrossRef]
  41. Barlat, F.; Lian, J. Plastic behaviour and stretchability of sheet metals (Part I): A yield function for orthotropic sheet under plane stress conditions. Int. J. Plast. 1989, 5, 51–56. [Google Scholar] [CrossRef]
  42. Li, H.; Hu, X.; Yang, H.; Li, L. Anisotropic and asymmetrical yielding and its distorted evolution: Modeling and applications. Int. J. Plast. 2017, 82, 127–158. [Google Scholar] [CrossRef]
  43. Hao, S.; Dong, X.H. Interpolation–based anisotropic yield and hardening models. Eur. J. Mech. A/Solids 2020, 83, 104047. [Google Scholar] [CrossRef]
  44. Fu, J.W.; Xie, W.W.; Zhou, J.M.; Qi, L.H. A method for the simultaneous identification of anisotropic yield and hardening constitutive parameters for sheet metal forming. Int. J. Mech. Sci. 2020, 181, 105756. [Google Scholar] [CrossRef]
  45. Li, Z.G.; Yang, H.F.; Liu, J.G.; Liu, F. An improved yield criterion characterizing the anisotropic and tension–compression asymmetric behavior of magnesium alloy. J. Magnes. Alloys 2022, 10, 569–584. [Google Scholar] [CrossRef]
  46. Lou, Y.S.; Zhang, S.J.; Yoon, J.W. A reduced Yld2004 function for modeling of anisotropic plastic deformation of metals under triaxial loading. Int. J. Mech. Sci. 2019, 161, 105027. [Google Scholar] [CrossRef]
  47. Wang, S.C.; Shang, H.C.; Zhang, Z.; Lou, Y.S. Multi–scale numerical investigation of deep drawing of 6K21 aluminum alloy by crystal plasticity and a stress-invariant based anisotropic yield function under non-associated flow rule. J. Manuf. Process. 2023, 102, 736–755. [Google Scholar] [CrossRef]
  48. Hou, Y.; Min, J.Y.; Guo, N.; Lin, J.P.; Carsley, J.E.; Stoughton, T.B.; Traphöner, H.; Clausmeyer, T.; Tekkaya, A.E. Investigation of evolving yield surfaces of dual-phase steels. J. Mater. Process. Technol. 2021, 287, 116314. [Google Scholar] [CrossRef]
  49. Mu, Z.K.; Zhao, J.; Meng, Q.D.; Sun, H.L.; Yu, G.C. Anisotropic hardening and evolution of r–values for sheet metal based on evolving non–associated Hill48 model. Thin-Walled Struct. 2022, 171, 108791. [Google Scholar] [CrossRef]
  50. Du, K.; Huang, S.H.; Hou, Y.; Wang, H.B.; Wang, Y.X.; Zheng, W.T.; Yuan, X.G. Characterization of the asymmetric evolving yield and flow of 6016–T4 aluminum alloy and DP490 steel. J. Mater. Sci. Technol. 2023, 133, 209–229. [Google Scholar] [CrossRef]
  51. Mamros, E.M.; Mayer, S.M.; Banerjee, D.K.; Iadicola, M.A.; Kinsey, B.L.; Ha, J.J. Plastic anisotropy evolution of SS316L and modeling for novel cruciform specimen. Int. J. Mech. Sci. 2022, 234, 107663. [Google Scholar] [CrossRef]
  52. Yang, H.C.; Zhang, W.; Zhuang, X.C.; Zhao, Z. Anisotropic plastic flow of low/medium carbon steel plates in different loading conditions: Characterization of the r–value. J. Mater. Process. Technol. 2023, 321, 118159. [Google Scholar] [CrossRef]
  53. Lee, H.R.; Lee, M.G.; Park, N. Effect of evolutionary anisotropic hardening on the prediction of deformation and forming load in incremental sheet forming simulation. Thin-Walled Struct. 2023, 193, 111231. [Google Scholar] [CrossRef]
  54. Gawad, J.; Banabic, D.; Bael, A.V.; Comsa, D.S.; Gologanu, M.; Eyckens, P.; Houtte, P.V.; Roose, D. An evolving plane stress yield criterion based on crystal plasticity virtual experiments. Int. J. Plast. 2015, 75, 141–169. [Google Scholar] [CrossRef]
  55. Cazacu, O. New yield criteria for isotropic and textured metallic materials. Int. J. Solids Struct. 2018, 139–140, 200–210. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of convexity analysis of yield surface under plane stress space by GINCA.
Figure 1. Schematic diagram of convexity analysis of yield surface under plane stress space by GINCA.
Materials 16 06828 g001
Figure 2. Mechanical curves of uniaxial tension along the RD, DD and TD for SS304: (a) load–stroke curve; (b) engineering stress–engineering strain curve.
Figure 2. Mechanical curves of uniaxial tension along the RD, DD and TD for SS304: (a) load–stroke curve; (b) engineering stress–engineering strain curve.
Materials 16 06828 g002
Figure 3. The hardening curves of unaixial tension along theRD, DD and TD for SS304.
Figure 3. The hardening curves of unaixial tension along theRD, DD and TD for SS304.
Materials 16 06828 g003
Figure 4. The limit values of the parameter c of CQN_Chen function computed by GINCA method under biaxial loading space for SS304.
Figure 4. The limit values of the parameter c of CQN_Chen function computed by GINCA method under biaxial loading space for SS304.
Materials 16 06828 g004
Figure 5. The evolution of 3D plane stress yield surface of CQN_Chen function at different λ ¯ for SS304.
Figure 5. The evolution of 3D plane stress yield surface of CQN_Chen function at different λ ¯ for SS304.
Materials 16 06828 g005aMaterials 16 06828 g005b
Figure 6. Comparison of 2D yield surfaces predicted by different models with experimental data for SS304.
Figure 6. Comparison of 2D yield surfaces predicted by different models with experimental data for SS304.
Materials 16 06828 g006
Figure 7. The 2D yield surface evolution predicted by the CQN_Chen function at 0 λ ¯ 0.6 with an interval of 0.02 for SS304.
Figure 7. The 2D yield surface evolution predicted by the CQN_Chen function at 0 λ ¯ 0.6 with an interval of 0.02 for SS304.
Materials 16 06828 g007
Figure 8. Comparison of uniaxial tensile yield stresses along distinct loading directions predicted by four models with experimental values.
Figure 8. Comparison of uniaxial tensile yield stresses along distinct loading directions predicted by four models with experimental values.
Materials 16 06828 g008
Figure 9. Comparison of hardening curves predicted by different models with experimental data under equibiaxial tension and uniaxial tension along the RD, DD and TD for SS304: (a) Hill48 function; (b) SY2009; (c) CQN; (d) CQN_Chen.
Figure 9. Comparison of hardening curves predicted by different models with experimental data under equibiaxial tension and uniaxial tension along the RD, DD and TD for SS304: (a) Hill48 function; (b) SY2009; (c) CQN; (d) CQN_Chen.
Materials 16 06828 g009
Figure 10. Errors of hardening curves predicted by different functions under equibiaxial tension and uniaxial tension along the RD, DD and TD for SS304.
Figure 10. Errors of hardening curves predicted by different functions under equibiaxial tension and uniaxial tension along the RD, DD and TD for SS304.
Materials 16 06828 g010
Table 1. The parameters of the Hockett–Sherby/Hollomon hardening law for SS304 uniaxial tensile hardening curves.
Table 1. The parameters of the Hockett–Sherby/Hollomon hardening law for SS304 uniaxial tensile hardening curves.
Uniaxial Tension σ ( λ ¯ ) = A B exp ( C λ ¯ b ) + K λ ¯ n Fitted Error
A (MPa)B (MPa)CbK (MPa)n
RD465.291148.28067.9084.1721499.1450.8788.1 × 10−4
DD459.826151.99061.7474.0681481.6430.8669.1 × 10−4
TD417.364103.482107.6774.4911577.8810.8705.1 × 10−4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, X.; Wang, S.; Xu, Z.; Zhou, J.; Wan, X.; Rayhan, H.M.A.; Lou, Y. Plastic Evolution Characterization for 304 Stainless Steel by CQN_Chen Model under the Proportional Loading. Materials 2023, 16, 6828. https://doi.org/10.3390/ma16216828

AMA Style

Gao X, Wang S, Xu Z, Zhou J, Wan X, Rayhan HMA, Lou Y. Plastic Evolution Characterization for 304 Stainless Steel by CQN_Chen Model under the Proportional Loading. Materials. 2023; 16(21):6828. https://doi.org/10.3390/ma16216828

Chicago/Turabian Style

Gao, Xiang, Songchen Wang, Zhongming Xu, Jia Zhou, Xinming Wan, Hasib Md Abu Rayhan, and Yanshan Lou. 2023. "Plastic Evolution Characterization for 304 Stainless Steel by CQN_Chen Model under the Proportional Loading" Materials 16, no. 21: 6828. https://doi.org/10.3390/ma16216828

APA Style

Gao, X., Wang, S., Xu, Z., Zhou, J., Wan, X., Rayhan, H. M. A., & Lou, Y. (2023). Plastic Evolution Characterization for 304 Stainless Steel by CQN_Chen Model under the Proportional Loading. Materials, 16(21), 6828. https://doi.org/10.3390/ma16216828

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop