Next Article in Journal
Electronic Structure of NdFeCoB Oxide Magnetic Particles Studied by DFT Calculations and XPS
Next Article in Special Issue
Forming Limit Analysis of Thin-Walled Extruded Aluminum Alloy Tubes under Nonlinear Loading Paths Using an Improved M-K Model
Previous Article in Journal
Dispersion and Polishing Mechanism of a Novel CeO2-LaOF-Based Chemical Mechanical Polishing Slurry for Quartz Glass
Previous Article in Special Issue
Optimization of Processing Parameter and Mechanical Response Analysis of Advanced Heterogeneous Laminated Composites Using Ni/Al Foils by In Situ Reaction Synthesis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Homogeneous Anisotropic Hardening Model in Plane Stress State for Sheet Metal under Nonlinear Loading Paths

1
School of Mechanical Engineering, Dalian University of Technology, Dalian 116024, China
2
Department of Integrated Systems Engineering, The Ohio State University, Columbus, OH 43210, USA
3
National Key Laboratory for Precision Hot Processing of Metals, Harbin Institute of Technology, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(3), 1151; https://doi.org/10.3390/ma16031151
Submission received: 12 December 2022 / Revised: 14 January 2023 / Accepted: 25 January 2023 / Published: 29 January 2023
(This article belongs to the Special Issue Mechanical Behavior of Advanced Engineering Materials)

Abstract

:
In sheet metal forming, the material is usually subjected to a complex nonlinear loading process, and the anisotropic hardening behavior of the material must be considered in order to accurately predict the deformation of the sheet. In recent years, the homogeneous anisotropic hardening (HAH) model has been applied in the simulation of sheet metal forming. However, the existing HAH model is established in the second-order stress deviator space, which makes the calculation complicated and costly, especially for a plane stress problem such as sheet metal forming. In an attempt to reduce the computational cost, an HAH model in plane stress state is proposed, and called the HAH-2d model in this paper. In the HAH-2d model, both the stress vector and microstructure vector contain only three in-plane components, so the calculation is significantly simplified. The characteristics of the model under typical nonlinear loading paths are analyzed. Additionally, the feasibility of the model is verified by the stress–strain responses of DP780 and EDDQ steel sheets under different two-step uniaxial tension tests. The results show that the HAH-2d model can reasonably reflect the Bauschinger effect and the permanent softening effect in reverse loading, and the latent hardening effect in cross loading, while the predictive accuracy for cross-loading softening remains to be improved. In the future, the HAH-2d model can be further modified to describe more anisotropic hardening behaviors and applied to numerical simulations.

1. Introduction

With increasing requirements for lightweight materials in the aerospace, aviation, and automotive industries, integrated thin-walled components with complex shapes are finding broader applications [1,2]. During the forming of an integrated component, the material often undergoes complex nonlinear loading paths, even including multi-step pre-deformation. The deformation behavior of sheet metal often exhibits strong path dependence [3,4,5] and pronounced anisotropic hardening behaviors in nonlinear loading paths; these may be, for example, the Bauschinger effect, the permanent softening effect and the latent hardening effect [6,7].
The Bauschinger effect is the phenomenon wherein the yield stress under reverse loading decreases after a certain amount of pre-deformation [8,9]; the permanent softening effect is another phenomenon during reverse loading wherein the reloading flow stress remains lower than the flow stress of monotonic loading [10]. Meanwhile, there is a wide range of strain path changes between monotonic loading and reverse loading. For some materials, the reloading flow stress at an angle with respect to the first loading may overshoot the monotonic stress–strain curve; this phenomenon is referred to the latent hardening effect [11,12].
Practically, the anisotropic hardening behavior significantly affects the spring-back and forming limit [13,14,15], and the reduction of spring-back and the avoidance of necking or fracture are of paramount importance in the forming of an integrated component. In order to guide process formulation and reduce the cost of process development, it is necessary to predict the deformation process and optimize the forming process by simulation. At present, the isotropic hardening assumption is commonly used in general simulations of metal forming, i.e., the subsequent yield surface expands uniformly based on the initial yield surface. However, an isotropic hardening model cannot reasonably capture the anisotropic hardening behavior under nonlinear loading paths [16,17]. Therefore, establishing advanced anisotropic hardening models to accurately predict the plastic deformation behavior of thin-walled metals under nonlinear loading conditions is an important research direction [18].
A variety of anisotropic hardening models have been proposed [19,20], for example, the kinematic hardening model established by introducing a back stress tensor into the isotropic hardening model. The kinematic hardening model was first proposed by Prager to capture the Bauschinger effect, wherein the yield surface translates with the deformation but its shape remains the same [21]. Subsequently, some kinematic hardening models with a nonlinear relationship between the increments of back stress and plastic strain were proposed to more accurately capture the flow stress under reverse loading or cyclical loading conditions [22,23]. In fact, the typical hardening behavior of metals is often a mixture of isotropic hardening and kinematic hardening, i.e., the yield surface changes both in size and position. Hence, the mixed hardening mode was proposed [24,25].
As an alternative to the kinematic hardening model, Barlat et al. [26] proposed the homogeneous anisotropic hardening (HAH) model, which is a distortional hardening model without a back stress tensor. In the HAH model, a fluctuating component is added into the traditional yield function to change the shape of yield surface, and a microstructure deviator related to the loading history is introduced to reflect the state of the microstructure. In order to further describe the latent hardening, more variables were introduced into the model [27,28]. More importantly, the HAH model can be regarded as a theoretical framework for the anisotropic hardening model, in which any homogeneous yield function and hardening law can be used, and the yield surface can be distorted depending on the specific loading path [17]. Some researchers have successfully implemented the HAH model in numerical simulations, and the prediction accuracy has improved clearly [29,30,31].
More recently, in order to capture complex and varied anisotropic hardening behavior and improve prediction accuracy, more parameters were introduced, for example, the HAH-20 model (using, at most, 23 coefficients [17]) and the HEXAH model (using, at most, 15 coefficients [32]). On the other hand, in order to reduce the computational cost in numerical simulation, some efficient algorithms were implemented, including a multi-step return-mapping algorithm for the HAH model [33], a fully implicit numerical algorithm that can solve a complete set of residuals [34], and a fast and robust stress-update algorithm based on the general cutting-plane method [35].
It should be pointed out that the HAH model and these improvements were proposed based on the second-order stress deviator with six independent components, corresponding to the fully 3D (three dimensional) stress state. However, in the plastic forming of thin-walled components, the normal stress is usually much smaller than the in-plane principal stresses, which means that the stress condition can be simplified to the plane stress state with only three components. Therefore, shell elements and material constitutive models for plane stress states are widely used in numerical simulations of sheet metal forming. When the HAH model is selected to analyze the deformation of sheet metals under nonlinear loading paths, it is required to transform the in-plane stress components to a stress deviator, then calculate the strain increments and new stress deviator after a deformation increment using an iterative method, and lastly, transform the new stress deviation tensor into in-plane stress components for output. Obviously, the plane stress problem is transformed into a more complex 3D stress problem, which makes the calculation complicated and costly. Therefore, it is of great significance to establish an anisotropic hardening model in plane stress state to analyze the deformation of sheet metals subjected to complex loading.
In this paper, an anisotropic hardening model for plane stress problems such as thin-walled metal forming will be proposed based on the framework of the HAH model. For convenience, the proposed model is named the HAH-2d model in this work. In the meantime, the performance of the HAH-2d model in predicting the evolution of yield loci and the stress–strain curve will be analyzed under typical nonlinear loading paths. Moreover, the new model will be verified by experimental results from two-step uniaxial tension tests.

2. Fundamentals of the HAH Model

The original HAH model was proposed by Barlat et al. in the second-order stress deviator space [26]. The yield function, as well as the plastic potential, is as follows:
Φ ( s ) = [ ϕ q + ϕ h q ] 1 q = [ ϕ q ( s ) + f 1 q | h ^ : s | h ^ : s | | q + f 2 q | h ^ : s + | h ^ : s | | q ] 1 q = σ ¯ ( ε ¯ )
where, s and h ^ are the stress deviator and microstructure deviator, respectively, and “:” denotes the double dot product.
The yield function Φ ( s ) consists of a stable component ϕ and a fluctuating component ϕ h . Any yield function may be used as the stable component after being reduced to a homogeneous function of degree 1 with the form of ϕ ( s ) = σ ¯ . σ ¯ ( ε ¯ ) is a hardening law, where σ ¯ and ε ¯ are the equivalent stress and equivalent plastic strain, respectively. f 1 and f 2 are state variables related to deformation history leading to the distortion of the yield surface, and q is a constant exponent. If f 1 = f 2 = 0 , Equation (1) reduces to a conventional yield function, Φ ( s ) = ϕ ( s ) = σ ¯ ( ε ¯ ) .
The microstructure deviator h ^ , a normalized tensorial state variable, is proposed to capture the deformation history and reflect the microstructure evolution of the material. It defines an axis, not a direction; namely, h ^ and h ^ represent the same microstructural state. Its initial value h ^ 0 is defined as the normalized stress deviator s ^ 0 corresponding to the initial yield, as in the following equation, where the factor 8/3 is used for convenience [26,27].
h ^ 0 = s ^ 0 = s ^ 0 8 3 s i j 0 s i j 0
If the material is reloaded in a different stress state, h ^ will rotate gradually towards the new stress deviator s when cos χ 0 , or towards −s when cos χ < 0 . Additionally, cos χ determined by Equation (3) is the cosine of the angle between s and h ^ , representing the variation of the loading path [27,28]. Monotonic, reverse and cross-loading sequences are, respectively, represented by cos χ = 1, −1 and 0.
cos χ = s : h ^ s h ^
where, s = s : s = s i j s i j denotes the norm of the second-order tensor s.

3. HAH-2d Model in Plane Stress State

3.1. Stress Vector and Microstructure Vector

During the forming process of sheet metal, the material is mainly in plane stress state; the three in-plane stress components can be simply expressed as a stress vector σ = ( σ 11 ,   σ 22 ,   σ 12 ) T , while the other stress components such as σ 33 , σ 13 and σ 23 are equal to 0. In order to describe the microstructure evolution of the sheet metal under plane stress conditions, a normalized microstructure vector H ^ = ( H ^ 11 ,   H ^ 22 ,   H ^ 12 ) T is defined in terms of stress vector in this paper. Its initial value H ^ 0 also corresponds to the normalized stress vector σ ^ 0 leading to the first increment of plastic deformation.
H ^ 0 = σ ^ 0 = σ 0 σ 0 : σ 0
Considering that the plastic deformation of metals is driven by the stress deviator s, it is still necessary to discuss the variation of loading path based on stress deviator when analyzing plastic deformation processes under plane stress conditions. The second-order deviatoric tensors, s and h, corresponding to the stress vector σ and the microstructure vector H ^ , can be obtained through the following two equations, respectively.
s = ( σ 11 σ 12 0 σ 12 σ 22 0 0 0 0 ) σ 11 + σ 22 3 ( 1 0 0 0 1 0 0 0 1 )
h = ( H ^ 11 H ^ 12 0 H ^ 12 H ^ 22 0 0 0 0 ) H ^ 11 + H ^ 22 3 ( 1 0 0 0 1 0 0 0 1 )
Then, the double dot product s : h can be expressed as:
s : h = 1 3 ( 2 σ 11 H ^ 11 + 2 σ 22 H ^ 22 σ 11 H ^ 22 σ 22 H ^ 11 + 6 σ 12 H ^ 12 ) = σ T D H ^
where the matrix D = 1 3 ( 2 1 0 1 2 0 0 0 6 ) .
Therefore, the value of cos χ can be calculated through vectors σ and H ^ , as follows:
cos χ = σ T D H ^ σ T D σ H ^ T D H ^
Moreover, H ^ defines an axis in the vector space of plane stress. When the sheet metal undergoes loading-path changes, H ^ will remain the same if cos χ = ± 1 , or rotate gradually towards the new stress vector σ if 0 cos χ < 1 , or rotate towards σ if 1 cos χ < 0 . In this work, a possible evolution for H ^ is given as follows:
{   d H ^ = k ( σ ^ H ^ ) d ε ¯ ,   cos χ 0   d H ^ = k ( σ ^ H ^ ) d ε ¯ ,   cos χ < 0
H i + 1 = H ^ i + d H ^
H ^ i + 1 = H i + 1 ( H 11 i + 1 ) 2 + ( H 22 i + 1 ) 2 + ( H 12 i + 1 ) 2
in which σ ^ is the normalized stress vector of σ , k is a constant that controls the rate of rotation, and d ε ¯ is the increment of equivalent strain ε ¯ , which can be calculated according to the plastic work increment d W p :
σ ¯ d ε ¯ = σ d ε = d W p
d ε = ( d ε 11 ,   d ε 22 ,   d γ 12 ) T
where d ε is the plastic strain increment vector, d ε 11 and d ε 22 are strain increments along directions 1 and 2, and d γ 12 is the shear strain increment.
The evolution of H ^ is shown schematically in Figure 1. H ^ rotates in the plane determined by H ^ i and σ ^ , and finally, H ^ takes the direction of σ ^ or σ ^ .

3.2. Formulation of the HAH-2d Model

3.2.1. Distortional Yield Function

The distortional yield function of the HAH model in the plane stress space (HAH-2d model) is as follows:
Φ ( σ ) = 1 F L [ ϕ q + ϕ h q ] 1 q = 1 F L [ ϕ q ( σ ) + f 1 q | σ p | σ p | | q + f 2 q | σ p + | σ p | | q ] 1 q = σ ¯ ( ε ¯ )
σ p = s cos χ = σ T D H ^ H ^ T D H ^
The stable component ϕ can be any homogeneous yield function of degree 1 for plane stress. The fluctuating component ϕ h is adopted to cover the Bauschinger and permanent softening effects, and the multiplier component 1 / F L is used to capture the latent hardening effect. Correspondingly, Φ ( σ ) is also a homogeneous function of degree 1, i.e., Φ ( a σ ) = a Φ ( σ ) holds for any real number a. In addition, σ p is the projection of the stress deviator s in the direction of the microstructure deviator h.
It is important to note that the associated flow rule is applied in this paper, the yield function also serves as the plastic potential. In other words, the plastic strain increment vector d ε is defined as
d ε = Φ ( σ ) σ d λ
where d λ 0 is the plastic multiplier.

3.2.2. Bauschinger Effect and Permanent Softening Effect

A schematic diagram of yield loci in the σ 11 σ 22 plane with f 1 > 0 , f 2 > 0 and F L = 1 is shown in Figure 2. σ f and σ r are yield stress vectors on the current distortional yield locus along the positive and negative directions of H ^ , respectively, and σ iso is the yield stress vector on the isotropic yield locus. State variables g 1 and g 2 are defined as Equation (17) to represent the relationship between σ f , σ r and σ iso , and the parameters f 1 and f 2 in Equation (14) can be expressed as Equation (18).
{ g 1 = σ r / σ iso 1 g 2 = σ f / σ iso 1
f i = 6 4 ( g i q 1 ) 1 q ,   i = 1 ,   2
For a sheet metal exhibiting the Bauschinger effect, the normalized yield surface in the opposite direction of loading tends to contract during forward loading, and the value of f 1 should increase from 0 towards a finite value. As a result, the yield stress of reverse loading is reduced compared with the final forward-loading stress, which is the Bauschinger effect. During reverse loading, the flat yield surface at the side of reverse loading tends to recover the yield surface determined by the isotropic hardening function ϕ ( σ ) = σ ¯ . At the same time, the yield surface opposite to the reverse loading contracts, and f 2 starts to increase in the same way as f 1 . If the yield surface at the side of reverse loading cannot recover to the level of isotropic hardening, the permanent softening effect will occur during reverse loading correspondingly.
A possible evolution for g 1 and g 2 is given in Equation (19) which has been successfully applied in existing HAH models [26,27,28]. State variables g 3 and g 4 represent the maximum saturation values of g 2 and g 1 , respectively. Material coefficients k1 and k2 control the evolution rate of g 1 and g 2 , and k3 controls their lower bound; k 4 controls the minimum value of g 3 and g 4 , and k 5 controls their evaluation rate. Meanwhile, the value of k 4 is usually slightly less than 1, according to the permanent softening effect. When the initial values of g 1 and g 2 are both 1 and k 1 = k 2 = 0 , the ability to describe the Bauschinger effect will be suppressed, and when g 3 = g 4 = 1 and k 4 = 1 (or k 5 = 0 ), the ability to describe the permanent softening effect will be suppressed.
If   cos χ 0 d g 1 d ε ¯ = k 2 ( k 3 σ ¯ ( 0 ) σ ¯ ( ε ¯ ) g 1 ) d g 2 d ε ¯ = k 1 g 3 g 2 g 2 d g 4 d ε ¯ = k 5 ( k 4 g 4 ) If   cos χ < 0 d g 1 d ε ¯ = k 1 g 4 g 1 g 1 d g 2 d ε ¯ = k 2 ( k 3 σ ¯ ( 0 ) σ ¯ ( ε ¯ ) g 2 ) d g 3 d ε ¯ = k 5 ( k 4 g 3 )

3.2.3. Latent Hardening Effect

For some materials, the flow stress during reloading may overshoot the stress during monotonic loading under the condition of the same equivalent strain. The phenomenon of overshooting was explained as the latent hardening effect. Additionally, the degree of overshooting is usually maximum for cross loading, in which the angle between the two corresponding normalized stress deviators s ^ 1 and s ^ 2 is 90° ( s ^ 1 : s ^ 2 = 0 ). As the degree of plastic deformation increases in the second step, h ^ rotates gradually towards s ^ 2 or s ^ 2 .
Suppose s v and s v are the stress deviators vertical to h ^ on the distortional yield surface and the isotropic one, respectively. Additionally, the corresponding plane stress vectors are σ v and σ v . Then, a state variable g L can be defined as Equation (20), which represents that the distortional yield surface at cos χ = 0 expands to g L times the yield surface under monotonic loading, and its minimum value is 1.
g L = ϕ ( σ v ) ϕ ( σ v ) = s v s v
Because the degree of the latent hardening effect is different under different strain path changes, it is important to determine the yield surface expansion for arbitrary χ . The stress deviator s can be decomposed into s parallel to h ^ and s perpendicular to h ^ , as shown in Figure 3 [28]. Considering that the latent hardening effect is mainly caused by the vertical component s , a new deviator s is composed as
s = s / g L + s
Then, the expansion ratio of the yield surface for arbitrary χ can be represented by
F L = s s = g L ( g L 2 1 ) cos 2 χ + 1
Thus, the surface of ϕ ( σ ) / F L = σ ¯ ( ε ¯ ) expands faster than that of ϕ ( σ ) = σ ¯ ( ε ¯ ) under any strain path change except cos χ = ± 1 . Furthermore, the evolution equation of g L is as follows [27]:
d g L d ε ¯ = k L [ ( σ ¯ ( ε ¯ ) σ ¯ ( 0 ) σ ¯ ( ε ¯ ) ) ( L ( 1 cos 2 χ ) + cos 2 χ 1 ) + 1 g L ]
where L and k L are constants controlling the upper bound and evolution rate of g L , respectively. When the initial value of g L is 1 and L = 1, F L will be identical to 1, and the latent hardening effect is suppressed.

3.3. Coefficient Identification

It is recommended to identify the coefficients in the HAH-2d model in the following order: yield function ϕ ( σ ) , isotropic hardening law σ ¯ ( ε ¯ ) , and distortional hardening [27]. The coefficients of ϕ ( σ ) and σ ¯ ( ε ¯ ) can be determined independently in a conventional way; for example, the coefficients of ϕ ( σ ) can be calculated based on the yield stresses and anisotropy coefficients along different directions of the sheet, and the coefficients of σ ¯ ( ε ¯ ) can be obtained by fitting a monotonic stress–strain curve. Additionally, the nine coefficients related to the distortion, namely, q, k, k1, k2, k3, k4, k5, kL and L, can be identified using an optimization method in which the input experimental data should come from reverse-loading and cross-loading tests.
When the stress–strain data of forward-reverse loading are available, coefficients k1, k2, k3, k4 and k5 can be determined independently. Otherwise, the five coefficients have to be evaluated together with k, kL, L and q. In particular, the Bauschinger, permanent softening and latent hardening effects can be expressed simultaneously in the stress–strain curve from a two-step tensile test with 1 < cos χ < 0 , for example, re-tension in the direction orthogonal to the fist tensile direction ( cos χ = 0.5 ). Given that the orthogonal two-step tensile test is easy to carry out, all the coefficients associated with the distortion could be identified based on the stress–strain curve.

4. Hardening Behavior Predicted under Typical Nonlinear Loading Paths

In this section, the stress–strain curves and yield locus evolutions of a generic sheet material under several typical nonlinear loading paths were predicted by the HAH-2d model. The typical nonlinear loading paths in the range of 1 cos χ < 1 include reverse loading ( cos χ = 1 ), cross loading ( cos χ = 0 ) and two-step tensions conducted at 45° ( cos χ = 0.25 ) and 90° ( cos χ = 0.5 ) from the first loading direction.
The initial yield condition of the generic sheet material is assumed to satisfy the Mises isotropic yield criterion, see Equation (24), and the hardening satisfies the Swift hardening law [36], see Equation (25). The material has a Bauschinger effect, permanent softening effect and latent hardening effect simultaneously. The coefficients of the swift hardening law and coefficients associated with the distortion are listed in Table 1.
ϕ Mises ( σ ) = σ 11 2 + σ 22 2 σ 11 σ 22 + 3 σ 12 2 = σ ¯
σ ¯ = K ( ε ¯ + ε 0 ) n
where K, ε 0 and n are material coefficients.

4.1. Reverse Loading

The evolution of the normalized yield locus of the generic material, predicted by the HAH-2d model under a reverse loading, namely, uniaxial compression (UC) to true strain of 0.10, followed by uniaxial tension (UT) in the reverse direction, is shown in Figure 4. The microstructure vector H ^ maintains (−1, 0, 0) during the whole deformation. The normalized yield loci in σ 11 σ 22 plane are divided into left and right sides by the line σ 22 / σ 11 = 2.0 . This line is chosen because the two stress deviators corresponding to stress vectors ( 1 ,   2 ,   0 ) T and ( 1 ,   0 ,   0 ) T are orthogonal. During the compression step, the yield locus on the left side remains unchanged, while the yield locus on the right side contracts inward, which will result in a reduction in re-loading yield stress when tension occurs in the opposite direction. As the tensile strain increases in the second step, the yield locus on the left side contracts inward instead. The yield locus on the right side recovers towards the isotropic hardening yield, but not completely, because of the permanent softening effect.
The predicted stress–strain curve during reverse loading is given in Figure 5. The reloading yield stress after compression pre-deformation is significantly lower than the monotonic flow stress under the same deformation level, which is the Bauschinger effect. When k4 = 0.9, the stress will not reach the monotonic uniaxial tensile curve, which is the permanent softening effect. Additionally, the permanent softening phenomenon disappears when k4 = 1.0. Therefore, the HAH-2d model can describe the Bauschinger effect and permanent softening effect during reverse loading.
The sensitivity of the HAH-2d model to parameters k1, k2, and k3 should be considered when predicting the Bauschinger effect during reverse loading. The state of the yield surface associated with the Bauschinger effect is reflected by state variables g1 and g2, and their evolutions follow the same rule. Therefore, the effects of k1, k2, and k3 on the predicted responses, such as the evaluations of g1 and the stress–strain curves, were analyzed as shown in Figure 6. It can be found from Figure 6a,b that k1 affects the evaluation rate of g1 during reloading. The greater k1 is, the greater the increase rate of g1, and the faster the flow stress reaches the saturation state. Figure 6c shows that k2 affects the evaluation rate of g1 during preloading. However, the corresponding stress–strain curves of reloading under the condition of 30 ≤ k2 ≤ 50, k1 = 100 and k3 = 0.5 are almost coincident, as shown in Figure 6d. The reason is that there is little difference in g1 after a pre-strain of 0.1, and the difference is quickly covered by the evolution during reverse loading. For the same reason, the difference among the curves in Figure 6f under different k3 is not obvious too. Even so, it should be noted that the value of g1 during preloading and the reloading yield stress are significantly affected by k3; the smaller the value of k3, the smaller the value of g1 and the reloading yield stress.
Figure 7 shows the effects of k4 and k5 on the predicted responses of the permanent softening effect during reverse loading, where the state variable g4 represents the maximum saturation value of the reloading stress. It can be seen that the permanent softening effect is mainly affected by k4, and the smaller the value of k4, the lower the value of g4 and the saturation value of reloading stress. In the meantime, when using a value of k4 less than 1.0, the larger the value of k5, the faster the value of g4 approaches the saturation value, and the lower the reloading stress–strain curve will be after a certain pre-strain. Of course, if the pre-strain is large enough, the influence of k5 on the reloading stress–strain curve will be significantly reduced or may even disappear.

4.2. Cross Loading

It is known that the two stress deviators corresponding to σ ^ PT = ( 1 / 5 ,   2 / 5 ,   0 ) T and σ ^ UT = ( 1 ,   0 ,   0 ) T are orthogonal to each other, and the corresponding deformation types are plane strain tension and uniaxial tension, respectively. Therefore, the loading sequence of a plane strain tension (PT) followed by a uniaxial tension at an angle of 90° from the first tension direction is a typical cross loading path, and cos χ = 0 at the start of the second loading.
The evolution of the normalized yield locus of the generic material, predicted by the HAH-2d model during uniaxial tension after plane strain tension to equivalent strain of 0.10, is shown in Figure 8. In the plane strain tension process, there is H ^ = σ ^ PT , the yield locus above the line σ 22 = 0 remains unchanged, and the yield locus below the line σ 22 = 0 contracts inward because of the Bauschinger effect. In the second loading step, H ^ rotates gradually towards σ ^ UT with the increase in tensile strain, and the yield locus near σ ^ UT expands outward first and then returns to the isotropic hardening yield locus.
The uniaxial tensile stress–strain curves in the second loading step after different degrees of plane strain pre-tension are shown in Figure 9. Obvious stress overshooting is observed, and the higher the pre-deformation the greater the overshooting. When the tensile strain reaches about 0.06 in the second step, the tensile stress returns to the level of single uniaxial tension. It indicates that the HAH-2d model can capture the latent hardening in cross loading.
The latent hardening effect predicted by the HAH-2d model during cross loading is affected mainly by parameters k, L, and kL. Therefore, the sensitivity of the model to these three parameters was analyzed. Figure 10 shows the predicted responses of latent hardening effect under different values of k during a typical cross loading, PT0.1−UT. With the increase in the reloading strain, the value of χ gradually decreases from 90° to 0°, and the value of gL rapidly increases to a certain value and then gradually decreases to 1.0. As a result, the phenomenon of stress overshooting is observed. Meanwhile, the larger the value of k, the faster the value of χ and g4 decrease, and the faster the overshoot stress recovers.
Figure 11 shows the effects of L and kL on the predicted responses of latent hardening effect during cross loading. It can be seen from Figure 11a,b that the larger the value L, the larger gL can be obtained and the more significant the stress overshooting. On the other hand, the latent hardening effect is also influenced by the value of kL, as shown in Figure 11c,d. Additionally, the larger the value of kL, the faster the stress value reaches its peak, and the larger the maximum stress value is. However, the sensitivity of the model to kL is much lower than that to L. In addition, neither L nor kL affects the strain span of the phenomenon of stress overshooting.

4.3. Two-Step Uniaxial Tension

Stress–strain curves in two-step uniaxial tensions with tensile axes at 45° and 90° from each other were predicted by the HAH-2d model and are shown in Figure 12 and Figure 13, respectively. The curves with 45° express only latent hardening, while the curves with 90° express the Bauschinger effect, permanent softening and latent hardening at the same time. Therefore, it is better to optimize the coefficients of the HAH-2d model by using the stress–strain data of a two-step uniaxial tension with 90°, or the orthogonal two-step tension mentioned in Section 3.3.

5. Model Validation

5.1. Determination of Material Coefficients

In this paper, the HAH-2d model will be applied to describe the hardening behavior of DP780 and EDDQ steel sheets with 1.2 mm thickness that were tested by Ha et al. using two-step uniaxial tension tests [6]. As shown in Figure 14, the sub-specimens were cut along different directions at the center of the pre-deformed big-specimen [6]. Additionally, the experimental results from tests consisted first of tension in the rolling direction (RD) and second of tension at 90° (TD, the transverse direction), 60° and 45° from RD were used in this work. Coefficients of the Yld2000-2d anisotropic yield function (see Appendix A) and the Swift hardening law (see Equation (25)) for the two materials are listed in Table 2, where m = 6 was adopted as recommended for BCC metals [6].
To quantitatively evaluate the error between the predicted and experimental data, the root mean square error (RMSE) and the average absolute relative error (AARE) were used in this paper. The definition of RMSE and AARE can be expressed by Equations (26) and (27).
RMSE = 1 N 1 N ( σ P σ E ) 2
AARE   ( % ) = 1 N 1 N | σ P σ E σ E | × 100 %
where σ P is the predicted stress, σ E is the experimental stress, and N is the number of data.
The coefficients associated with the distortion of the HAH-2d model for each material were identified with an optimization method using the stress–strain data from an orthogonal two-step tensile test with 10% pre-strain in RD, as shown in Figure 15. In the optimization method, the objective is to minimize the RMSE. Considering that the DP780 steel expresses no permanent softening and latent hardening, k4 = 1.0, k5 = 0, kL = 0 and L = 1.0 were set directly. Similarly, as the EDDQ steel expresses no permanent softening, k4 = 1.0 and k5 = 0 were set. In the meantime, the value of q was set to 2.0 for convenience. The obtained coefficients of the HAH-2d models for DP780 and EDDQ are listed in Table 3. Additionally, the predicted stress–strain curves are plotted in Figure 15a,b, both of which agree well with the experimental data. Additionally, the values of RMSE for DP780 and EDDQ are 9.81 and 1.85, respectively. Incidentally, the reason for using the experimental data from reference [27] in Figure 15 is that the experimental data used in reference [27] are also from reference [6], but more points are given for this testing condition.

5.2. Comparison of Prediction and Experimental Results

In order to validate the HAH-2d model, the hardening curves in two-step uniaxial tensions of DP780 and EDDQ sheets, which consisted of 4% or 10% pre-strain in RD and re-tension at 90° (TD), 60° or 45° from RD, were predicted and compared with the experimental results. Incidentally, the two-step tension with tensile axes at 60° from each other can be called “pseudo cross loading”, because it is close to the ideal cross-loading condition, about 55° between the first and second tensile directions.
The predicted stress–strain curves of the HAH-2d model for DP780 under different two-step uniaxial tensions are compared with the experimental curves in Figure 16, and the corresponding prediction errors are listed in Table 4. It can be seen from Figure 16a that the HAH-2d model captures the transient hardening behavior of DP780 fairly well for re-tension in TD after a 4% pre-strain in RD, and the AARE value is only 1.28%. However, the prediction at 60° shown in Figure 16b slightly overestimates the re-loading yield stress and leads to a lower hardening rate. As a result, the values of RMSE and AARE increase to 37.85 and 3.44% under the condition of 10% pre-strain. On the other hand, as shown in Figure 16c, the phenomenon of cross-loading softening was observed in re-loading stress–strain curves at 45°, but the predicted re-loading yield stresses at 45° reach the monotonic curve immediately. Additionally, the high RMSE values of 24.49 and 47.13 are mainly due to the significant difference between the predicted stress and the experimental stress at the beginning of the second tension.
Figure 17a–c show the predicted reloading flow stress curves for EDDQ in different two-step uniaxial tensions with pre-tension in RD and re-tension at 90° (TD), 60° and 45° from RD. The predicted curves capture all the experimental results well, even though there is a little overestimation of the flow stress in the re-tension at 60° from RD. The corresponding prediction errors are listed in Table 5. The maximum values of RMSE and AARE, 9.86 and 2.00%, both occur in the condition of re-tension at 60° after a 10% pre-tension. This indicates that the HAH-2d model can capture anisotropic hardening responses such as flow stress overshooting and strain-hardening stagnation.

6. Conclusions

In this work, an HAH model in plane stress state was proposed, and the evolution of the yield loci and the stress–strain curve of a generic material under some typical nonlinear loading paths was analyzed using this model. Furthermore, the model was validated by experimental results from two-step uniaxial tension tests. The conclusions are as follows:
(1) The HAH-2d model was proposed in plane stress state. Both the loading stress vector and microstructure vector contain only three in-plane components. Compared with the HAH model in fully 3D space, the transformation between the plane stress vector and the second-order stress deviator and the complicated calculation in the stress deviator space can be avoided in the HAH-2d model. The computational cost in finite element applications is therefore reduced.
(2) The Bauschinger effect and permanent softening effect in reverse loading and the latent hardening effect in cross loading were investigated in this paper. The evolution of the yield surface and the stress–strain curve of the sheet metal with these effects can be predicted by the HAH-2d model.
(3) The HAH-2d model can reasonably predict experimental phenomena such as the reduction of re-loading yield stress in two-step loading with cos χ < 0 and stress overshooting in cross loading. However, the phenomenon of cross-loading softening observed in two-step loading with 0 cos χ < 1 cannot be captured by the model.

Author Contributions

Conceptualization, H.Z. and Z.H..; methodology, H.Z. and Y.L.; validation, Y.L. and K.C.; formal analysis, Y.L. and K.C.; investigation, H.Z., Y.L. and Z.H.; resources, Z.H. and S.Y.; data curation, H.Z. and Y.L.; writing—original draft preparation, H.Z., Y.L. and K.C.; writing—review and editing, H.Z., Y.L., K.C. and Z.H.; supervision, Z.H. and S.Y.; project administration, S.Y.; funding acquisition, Z.H. 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 numbers: 52075075, 51575131; and the National Key R&D Program of China, grant numbers: 2017YFB0304404.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The present work was financed by the National Natural Science Foundation of China (52075075, 51575131) and the National Key R&D Program of China (2017YFB0304404). The authors would like express their sincere appreciation to the funds.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Yld2000-2d anisotropic yield criterion [37] is defined as:
ϕ = | X 1 X 2 | m + | X 1 + 2 X 2 | m + | 2 X 1 + X 2 | m = 2 σ ¯ m
where σ ¯ is the equivalent stress, X i and X i (i = 1, 2) are the principal values of tensors X = L σ and X = L σ , respectively. σ is Cauchy stress tensor. L and L can be written in the following form
[ L ] = [ L 11 L 12 0 L 21 L 22 0 0 0 L 66 ] = 1 3 [ 2 α 1 α 1 0 α 2 2 α 2 0 0 0 α 7 ]
[ L ] = [ L 11 L 12 0 L 21 L 22 0 0 0 L 66 ] = 1 9 [ 8 α 5 2 α 3 2 α 6 + 2 α 4 4 α 6 4 α 4 4 α 5 + α 3 0 4 α 3 4 α 5 4 α 4 + α 6 8 α 4 2 α 6 2 α 3 + 2 α 5 0 0 0 9 α 8 ]
where α i (i = 1, 2…8) are pending coefficients.
The principal values of X are
{ X 1 = 1 2 ( X 11 + X 22 + ( X 11 X 22 ) 2 + 4 X 12 2 ) X 2 = 1 2 ( X 11 + X 22 ( X 11 X 22 ) 2 + 4 X 12 2 )
In the Yld2000-2d criterion, only the plane stress components, σ x x , σ y y and σ x y are relevant, where x and y indicate the RD and TD, respectively. For the exponent m in Equation (A1), it is recommended to use 6 for BCC material and 8 for FCC material.

References

  1. Yuan, S.; Fan, X. Developments and perspectives on the precision forming processes for ultra-large size integrated components. Int. J. Extrem. Manuf. 2019, 1, 022002. [Google Scholar] [CrossRef]
  2. Zheng, K.; Politis, D.J.; Wang, L.; Lin, J. A review on forming techniques for manufacturing lightweight complex—Shaped aluminium panel components. Int. J. Lightweight Mater. Manuf. 2018, 1, 55–80. [Google Scholar] [CrossRef]
  3. Bouvier, S.; Gardey, B.; Haddadi, H.; Teodosiu, C. Characterization of the strain-induced plastic anisotropy of rolled sheets by using sequences of simple shear and uniaxial tensile tests. J. Mater. Process. Tech. 2006, 174, 115–126. [Google Scholar] [CrossRef]
  4. Qin, J.; Holmedal, B.; Hopperstad, O.S. A combined isotropic, kinematic and distortional hardening model for aluminum and steels under complex strain-path changes. Int. J. Plast. 2018, 101, 156–169. [Google Scholar] [CrossRef]
  5. He, J.; Gu, B.; Li, Y.; Li, S. Forming limits under stretch-bending through distortionless and distortional anisotropic hardening. J. Manuf. Sci. Eng.: Trans. ASME 2018, 140, 121013. [Google Scholar] [CrossRef]
  6. Ha, J.; Lee, M.; Barlat, F. Strain hardening response and modeling of EDDQ and DP780 steel sheet under non-linear strain path. Mech. Mater. 2013, 64, 11–26. [Google Scholar] [CrossRef]
  7. Mánik, T.; Holmedal, B.; Hopperstad, O.S. Strain-path change induced transients in flow stress, work hardening and r-values in aluminum. Int. J. Plast. 2015, 69, 1–20. [Google Scholar] [CrossRef] [Green Version]
  8. Bauschinger, J. Changes of the elastic limit and the modulus of elasticity on various metals. Ziivilingenieur 1881, 27, 289–348. [Google Scholar]
  9. Hasegawa, T.; Yakou, T.; Karashima, S. Deformation behaviour and dislocation structures upon stress reversal in polycrystalline aluminium. Mater. Sci. Eng. 1975, 20, 267–276. [Google Scholar] [CrossRef]
  10. Hasegawa, T.; Yakou, T. The effect of strain reversal and thermal recovery on stress vs strain behavior in aluminum. Scr. Metall. 1980, 14, 1083–1087. [Google Scholar] [CrossRef]
  11. Schmitt, J.H.; Fernandes, J.V.; Gracio, J.J.; Viera, M.F. Plastic behaviour of copper sheets during sequential tension tests. Mater. Sci. Eng. A 1991, 147, 143–154. [Google Scholar] [CrossRef]
  12. Barlat, F.; Ferreira Duarte, J.M.; Gracio, J.J.; Lopes, A.B.; Rauch, E.F. Plastic flow for non-monotonic loading conditions of an aluminum alloy sheet sample. Int. J. Plast. 2003, 19, 1215–1244. [Google Scholar] [CrossRef]
  13. Ben Bettaieb, M.; Abed-Meraim, F. Localized necking in elastomer-supported metal layers: Impact of kinematic hardening. J. Manuf. Sci. Eng. Trans. ASME 2017, 139, 61008. [Google Scholar] [CrossRef]
  14. Min, J.; Guo, N.; Hou, Y.; Jiang, K.; Chen, X.; Carsley, J.E.; Lin, J. Effect of tension-compression testing strategy on kinematic model calibration and springback simulation of advanced high strength steels. Int. J. Mater. Form. 2020, 14, 435–448. [Google Scholar] [CrossRef]
  15. Yetna N’Jock, M.; Badreddine, H.; Labergere, C.; Yue, Z.; Saanouni, K.; Dang, V.T. An application of fully coupled ductile damage model considering induced anisotropies on springback prediction of advanced high strength steel materials. Int. J. Mater. Form. 2021, 14, 739–752. [Google Scholar] [CrossRef]
  16. Clausmeyer, T.; Güner, A.; Tekkaya, A.E.; Levkovitch, V.; Svendsen, B. Modeling and finite element simulation of loading-path-dependent hardening in sheet metals during forming. Int. J. Plast. 2014, 63, 64–93. [Google Scholar] [CrossRef]
  17. Barlat, F.; Yoon, S.; Lee, S.; Wi, M.; Kim, J. Distortional plasticity framework with application to advanced high strength steel. Int. J. Solids Struct. 2020, 202, 947–962. [Google Scholar] [CrossRef]
  18. Cao, J.; Banu, M. Opportunities and challenges in metal forming for lightweighting: Review and future work. J. Manuf. Sci. Eng. Trans. ASME 2020, 142, 110813. [Google Scholar] [CrossRef]
  19. Lee, M.G.; Kim, C.; Pavlina, E.J.; Barlat, F. Advances in sheet forming—Materials modeling, numerical simulation, and press technologies. J. Manuf. Sci. Eng. Trans. ASME 2011, 133, 61001. [Google Scholar] [CrossRef]
  20. Banabic, D.; Barlat, F.; Cazacu, O.; Kuwabara, T. Advances in anisotropy of plastic behaviour and formability of sheet metals. Int. J. Mater. Form. 2020, 13, 749–787. [Google Scholar] [CrossRef]
  21. Prager, W. A new method of analyzing stresses and strains in work-hardening plastic solids. J. Appl. Mech. 1956, 23, 493–496. [Google Scholar] [CrossRef]
  22. Frederick, C.; Armstrong, P. A mathematical representation of the multiaxial Bauschinger effect. Mater. High Temp. 2007, 24, 1–26. [Google Scholar] [CrossRef]
  23. Chaboche, J.L. Time-independent constitutive theories for cyclic plasticity. Int. J. Plast. 1986, 2, 149–188. [Google Scholar] [CrossRef]
  24. Bouhamed, A.; Jrad, H.; Said, L.B.; Wali, M.; Dammak, F. A non-associated anisotropic plasticity model with mixed isotropic-kinematic hardening for finite element simulation of incremental sheet metal forming process. Int. J. Adv. Manuf. Technol. 2018, 100, 929–940. [Google Scholar] [CrossRef]
  25. Grilo, T.J.; Valente, R.A.F.; Alves De Sousa, R.J. Modelling non-quadratic anisotropic yield criteria and mixed isotropic-nonlinear kinematic hardening: Analysis of forward- and backward-Euler formulations. Int. J. Mater. Form. 2014, 8, 533–547. [Google Scholar] [CrossRef]
  26. Barlat, F.; Gracio, J.J.; Lee, M.; Rauch, E.F.; Vincze, G. An alternative to kinematic hardening in classical plasticity. Int. J. Plast. 2011, 27, 1309–1327. [Google Scholar] [CrossRef]
  27. Barlat, F.; Ha, J.; Grácio, J.J.; Lee, M.; Rauch, E.F.; Vincze, G. Extension of homogeneous anisotropic hardening model to cross-loading with latent effects. Int. J. Plast. 2013, 46, 130–142. [Google Scholar] [CrossRef]
  28. Barlat, F.; Vincze, G.; Grácio, J.J.; Lee, M.G.; Rauch, E.F.; Tomé, C.N. Enhancements of homogenous anisotropic hardening model and application to mild and dual-phase steels. Int. J. Plast. 2014, 58, 201–218. [Google Scholar] [CrossRef]
  29. Choi, J.; Lee, J.; Bae, G.; Barlat, F.; Lee, M. Evaluation of springback for DP980 S rail using anisotropic hardening models. JOM 2016, 68, 1850–1857. [Google Scholar] [CrossRef]
  30. Chen, S.; Liao, J.; Xiang, H.; Xue, X.; Pereira, A.B. Pre-strain effect on twist springback of a 3D P-channel in deep drawing. J. Mater. Process. Technol. 2021, 287, 116224. [Google Scholar] [CrossRef]
  31. Choi, H.; Choi, S.; Kang, S.; Kim, C.; Lee, M. Prediction of anisotropic strengths of steel plate after prior bending-reverse bending deformation: Application of distortional hardening model. Int. J. Mech. Sci. 2021, 204, 106512. [Google Scholar] [CrossRef]
  32. Reyne, B.; Barlat, F. A new concept for continuum distortional plasticity. Int. J. Plast. 2022, 155, 103303. [Google Scholar] [CrossRef]
  33. Hc, A.; Jwya, B. Stress integration-based on finite difference method and its application for anisotropic plasticity and distortional hardening under associated and non-associated flow rules. Comput. Methods Appl. Mech. Eng. 2019, 345, 123–160. [Google Scholar] [CrossRef]
  34. Choi, H.; Choi, S.; Kang, S.C.; Lee, M.G. Fully Implicit stress update algorithm for distortion-based anisotropic hardening with cross-loading effect: Comparative algorithmic study and application to large-size forming problem. Appl. Sci. 2021, 11, 5509. [Google Scholar] [CrossRef]
  35. Lee, J.; Bong, H.J.; Ha, J. Efficient and robust stress integration algorithm for anisotropic distortional hardening law under cross-loading with latent hardening. Eur. J. Mech.-A/Solids 2022, 96, 104775. [Google Scholar] [CrossRef]
  36. Swift, H. Plastic instability under plane stress. J. Mech. Phys. Solids 1952, 1, 1–18. [Google Scholar] [CrossRef]
  37. 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]
Figure 1. Schematic diagram of the evolution of microstructure vector H ^ .
Figure 1. Schematic diagram of the evolution of microstructure vector H ^ .
Materials 16 01151 g001
Figure 2. Schematic diagram of yield loci in σ 11 σ 22 plane with f 1 > 0 , f 2 > 0 and F L = 1 .
Figure 2. Schematic diagram of yield loci in σ 11 σ 22 plane with f 1 > 0 , f 2 > 0 and F L = 1 .
Materials 16 01151 g002
Figure 3. Schematic of the decomposition and composition of stress deviator [28]. Redraw with reference to International Journal of Plasticity, Vol 58, Barlat, F.; Vincze, G.; Grácio, J.J.; Lee, M.G.; Rauch, E.F.; Tomé, C.N., Enhancements of homogenous anisotropic hardening model and application to mild and dual-phase steels, 201–218, Copyright Elsevier (2014), with permission from Elsevier.
Figure 3. Schematic of the decomposition and composition of stress deviator [28]. Redraw with reference to International Journal of Plasticity, Vol 58, Barlat, F.; Vincze, G.; Grácio, J.J.; Lee, M.G.; Rauch, E.F.; Tomé, C.N., Enhancements of homogenous anisotropic hardening model and application to mild and dual-phase steels, 201–218, Copyright Elsevier (2014), with permission from Elsevier.
Materials 16 01151 g003
Figure 4. Evolution of the normalized yield locus of the generic material predicted by the HAH-2d model in a uniaxial compression-tension loading path.
Figure 4. Evolution of the normalized yield locus of the generic material predicted by the HAH-2d model in a uniaxial compression-tension loading path.
Materials 16 01151 g004
Figure 5. True stress–strain curves of the generic material predicted by the HAH-2d model in a uniaxial compression-tension loading path.
Figure 5. True stress–strain curves of the generic material predicted by the HAH-2d model in a uniaxial compression-tension loading path.
Materials 16 01151 g005
Figure 6. Predicted responses of the Bauschinger effect during reverse loading; evaluation of g1 during (a) UC0.1−UT loading under different k1; (c) UC loading under different k2; and (e) UC loading under different k3; and stress–strain curves during UC0.1−UT loading under (b) different k1; (d) different k2 and (f) different k3.
Figure 6. Predicted responses of the Bauschinger effect during reverse loading; evaluation of g1 during (a) UC0.1−UT loading under different k1; (c) UC loading under different k2; and (e) UC loading under different k3; and stress–strain curves during UC0.1−UT loading under (b) different k1; (d) different k2 and (f) different k3.
Materials 16 01151 g006
Figure 7. Predicted responses of permanent softening effect during reverse loading; evaluation of g4 during UC loading under (a) different k4 and (c) different k5; and stress–strain curves during UC0.1−UT loading under (b) different k4 and (d) different k5.
Figure 7. Predicted responses of permanent softening effect during reverse loading; evaluation of g4 during UC loading under (a) different k4 and (c) different k5; and stress–strain curves during UC0.1−UT loading under (b) different k4 and (d) different k5.
Materials 16 01151 g007aMaterials 16 01151 g007b
Figure 8. Evolution of the normalized yield locus of the generic material predicted by HAH-2d model in plane strain tension–uniaxial tension.
Figure 8. Evolution of the normalized yield locus of the generic material predicted by HAH-2d model in plane strain tension–uniaxial tension.
Materials 16 01151 g008
Figure 9. Stress–strain curves predicted by the HAH-2d model in plane strain tension–uniaxial tension with different pre-strain levels.
Figure 9. Stress–strain curves predicted by the HAH-2d model in plane strain tension–uniaxial tension with different pre-strain levels.
Materials 16 01151 g009
Figure 10. Predicted responses of latent hardening effect during cross loading of PT0.1−UT under different k: (a) evaluation of χ, (b) evaluation of gL and (c) stress–strain curves.
Figure 10. Predicted responses of latent hardening effect during cross loading of PT0.1−UT under different k: (a) evaluation of χ, (b) evaluation of gL and (c) stress–strain curves.
Materials 16 01151 g010
Figure 11. Predicted responses of latent hardening effect during cross loading of PT0.1−UT; evaluation of gL under (a) different L and (c) different kL; and stress–strain curves under (b) different L and (d) different kL.
Figure 11. Predicted responses of latent hardening effect during cross loading of PT0.1−UT; evaluation of gL under (a) different L and (c) different kL; and stress–strain curves under (b) different L and (d) different kL.
Materials 16 01151 g011
Figure 12. Stress–strain curves predicted by the HAH-2d model in two-step uniaxial tensions with tensile axes at 45° from each other.
Figure 12. Stress–strain curves predicted by the HAH-2d model in two-step uniaxial tensions with tensile axes at 45° from each other.
Materials 16 01151 g012
Figure 13. Stress–strain curves predicted by the HAH-2d model in two-step uniaxial tensions with tensile axes at 90° from each other.
Figure 13. Stress–strain curves predicted by the HAH-2d model in two-step uniaxial tensions with tensile axes at 90° from each other.
Materials 16 01151 g013
Figure 14. Schematic illustration of the specimens used in the two-step uniaxial tension tests: (a) the sub-specimen at the center of the big-specimen and (b) directions of the second tension [6].
Figure 14. Schematic illustration of the specimens used in the two-step uniaxial tension tests: (a) the sub-specimen at the center of the big-specimen and (b) directions of the second tension [6].
Materials 16 01151 g014
Figure 15. Experimental and predicted stress–strain curves of re-tension in TD with 10% uniaxial pre-strain in RD for: (a) DP780 steel and (b) EDDQ steel. The experimental data are derived from reference [27]: International Journal of Plasticity, Vol 46, Barlat, F.; Ha, J.; Grácio, J.J.; Lee, M.; Rauch, E.F.; Vincze, G., Extension of homogeneous anisotropic hardening model to cross-loading with latent effects, 130–142, Copyright Elsevier (2013), with permission from Elsevier.
Figure 15. Experimental and predicted stress–strain curves of re-tension in TD with 10% uniaxial pre-strain in RD for: (a) DP780 steel and (b) EDDQ steel. The experimental data are derived from reference [27]: International Journal of Plasticity, Vol 46, Barlat, F.; Ha, J.; Grácio, J.J.; Lee, M.; Rauch, E.F.; Vincze, G., Extension of homogeneous anisotropic hardening model to cross-loading with latent effects, 130–142, Copyright Elsevier (2013), with permission from Elsevier.
Materials 16 01151 g015
Figure 16. Experimental and predicted stress–strain curves of re-tension at (a) 90°, (b) 60° and (c) 45° from RD with 4% or 10% uniaxial pre-strain in RD for DP780. The experimental data are derived from reference [6]: Mechanics of Materials, Vol 64, Ha, J.; Lee, M.; Barlat, F., Strain hardening response and modeling of EDDQ and DP780 steel sheet under non-linear strain path, 11–26, Copyright Elsevier (2013), with permission from Elsevier.
Figure 16. Experimental and predicted stress–strain curves of re-tension at (a) 90°, (b) 60° and (c) 45° from RD with 4% or 10% uniaxial pre-strain in RD for DP780. The experimental data are derived from reference [6]: Mechanics of Materials, Vol 64, Ha, J.; Lee, M.; Barlat, F., Strain hardening response and modeling of EDDQ and DP780 steel sheet under non-linear strain path, 11–26, Copyright Elsevier (2013), with permission from Elsevier.
Materials 16 01151 g016
Figure 17. Experimental and predicted stress–strain curves of re-tension at (a) 90°, (b) 60° and (c) 45° from RD with 4% or 10% uniaxial pre-strain in RD for EDDQ. The experimental data are derived from reference [6]: Mechanics of Materials, Vol 64, Ha, J.; Lee, M.; Barlat, F., Strain hardening response and modeling of EDDQ and DP780 steel sheet under non-linear strain path, 11–26, Copyright Elsevier (2013), with permission from Elsevier.
Figure 17. Experimental and predicted stress–strain curves of re-tension at (a) 90°, (b) 60° and (c) 45° from RD with 4% or 10% uniaxial pre-strain in RD for EDDQ. The experimental data are derived from reference [6]: Mechanics of Materials, Vol 64, Ha, J.; Lee, M.; Barlat, F., Strain hardening response and modeling of EDDQ and DP780 steel sheet under non-linear strain path, 11–26, Copyright Elsevier (2013), with permission from Elsevier.
Materials 16 01151 g017
Table 1. Coefficients of the generic material.
Table 1. Coefficients of the generic material.
CoefficientsSwift Hardening Law Coefficients Associated with the Distortion
K (MPa)ε0n qkk1k2k3k4k5LkL
Values5000.010.25 220100500.50.9202.0300
Table 2. Coefficients of yield function and hardening law for DP780 and EDDQ (Data from [6]).
Table 2. Coefficients of yield function and hardening law for DP780 and EDDQ (Data from [6]).
MaterialYld2000-2d Yield Function Swift Hardening Law
α1α2α3α4α5α6α7α8m K (MPa)ε0n
DP7800.9461.0221.0151.0001.0110.9681.0101.0066 12950.00080.142
EDDQ1.0141.1180.9310.8920.9040.8111.0290.9186 5380.00750.267
Table 3. Coefficients of HAH-2d models for DP780 and EDDQ.
Table 3. Coefficients of HAH-2d models for DP780 and EDDQ.
Materialkk1k2k3k4k5kLLq
DP78045.0135.239.60.4751.0001.02.0
EDDQ12.80773.9295.90.5511.00205.92.2262.0
Table 4. Prediction errors for DP780 under different two-step uniaxial tensions.
Table 4. Prediction errors for DP780 under different two-step uniaxial tensions.
ConditionsRD 4%−TDRD 4%–60°RD 10%–60°RD 4%–45°RD 10%–45°
RMSE20.1124.0737.8524.4947.13
AARE (%)1.282.143.441.942.58
Table 5. Prediction errors for EDDQ under different two-step uniaxial tensions.
Table 5. Prediction errors for EDDQ under different two-step uniaxial tensions.
ConditionsRD 4%−TDRD 4%–60°RD 10%–60°RD 4%–45°RD 10%–45°
RMSE5.598.659.864.745.40
AARE (%)1.531.782.001.331.13
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

Zhu, H.; Lin, Y.; Chen, K.; He, Z.; Yuan, S. A Homogeneous Anisotropic Hardening Model in Plane Stress State for Sheet Metal under Nonlinear Loading Paths. Materials 2023, 16, 1151. https://doi.org/10.3390/ma16031151

AMA Style

Zhu H, Lin Y, Chen K, He Z, Yuan S. A Homogeneous Anisotropic Hardening Model in Plane Stress State for Sheet Metal under Nonlinear Loading Paths. Materials. 2023; 16(3):1151. https://doi.org/10.3390/ma16031151

Chicago/Turabian Style

Zhu, Haihui, Yanli Lin, Kelin Chen, Zhubin He, and Shijian Yuan. 2023. "A Homogeneous Anisotropic Hardening Model in Plane Stress State for Sheet Metal under Nonlinear Loading Paths" Materials 16, no. 3: 1151. https://doi.org/10.3390/ma16031151

APA Style

Zhu, H., Lin, Y., Chen, K., He, Z., & Yuan, S. (2023). A Homogeneous Anisotropic Hardening Model in Plane Stress State for Sheet Metal under Nonlinear Loading Paths. Materials, 16(3), 1151. https://doi.org/10.3390/ma16031151

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