1. Introduction
It has been long observed that plastic deformations induce anisotropy in initially isotropic materials [
1,
2], which has been a challenge for the traditional plastic constitutive theory to describe. The study of the evolution law of yield surface relies on the establishment of a plastic constitutive model, which serves as a basis for the analysis of macroscopic test data and verification of the applicability and validity of the model by multiaxial tests. As loading applied to the material increases, further deformation takes place, the stress state reaches the subsequent yield point, and there exists a subsequent yield surface for three-dimensional loading in the stress space. The evolution forms of a subsequent yield surface are decomposed into expansion and movement of the yield surface in classical plasticity theory, which correspond to isotropic hardening and kinematic hardening, respectively. Isotropic hardening and kinematic hardening can be characterized by a yield surface radius and a back stress in the von Mises theory. Until now, researchers have conducted many studies on plastic behaviors of metals [
3,
4,
5,
6,
7].
Chaboche [
8,
9,
10] proposed a more general cyclic plasticity constitutive model on the basis of Armstrong and Frederick (A-F model) [
11], which can better reflect the cyclic hardening and cyclic softening behaviors of materials under proportional loading. The experiments on various metals from Phillips [
12], Wu [
13], and Khan [
14] showed distortion of the yield surface characterized by a region of high curvature along the direction of loading and a region of flattening on the opposite side. Moreover, the distortion of the subsequent yield surface varies for different materials. Hu Guijuan and Zhang Keshi [
15,
16,
17] systematically studied the evolution behaviors of the yield surfaces for copper, steel, and aluminum under tension-shear loading. It was shown that the shape of the measured yield surface is related to the offset strain of the yield point. The subsequent yield surface defined by a small offset strain develops an obvious “sharp point”. The shape of the yield surface is close to the cylindrical surface with increases of offset strain, which is consistent with the circle yield surface described by Mises yield criterion. Typical distortional subsequent yield surfaces under pretension and pretorsion loadings are shown in
Figure 1 and
Figure 2 [
17] for 45 steel.
Since the initial and subsequent yield surfaces described by the von Mises theory are cylindrical surfaces in principal stress space, the traditional models cannot describe the cross effect (the sharp point in pre-deformation direction and the relatively flat bottom at its opposite direction) of a yield surface accurately. Many researchers have studied the anisotropic evolution of a yield surface by the macroscopic method [
18,
19,
20,
21] and the mesoscopic method [
22,
23], respectively. In classical plasticity theory, the anisotropy is reflected by the translation and distortion of the yield surface. Francois [
24] modified the classical von Mises yield criterion by replacing the usual stress deviator with distortional stress tensor
and described the yield surface as an “egg shape” to reduce the difference between simulation and experiment. Shi Baodong et al. [
25] improved the theoretical model by supplementing the distortion factor with the classical plastic constitutive theory. Fu Qiang et al. [
26] established a slip element model to simulate the subsequent yield surface evolution based on the mechanism of slip deformation. Chen Yuan et al. [
27] proposed the mathematical models for various deformation mechanisms of slip, twinning, and detwinning to capture the Bauschinger effect under cyclic loading. Wen et al. [
28] used the self-consistent model to determine the yield characteristics of polycrystalline materials under different loading conditions. The crystal plasticity method is powerful enough to simulate the inhomogeneous deformation and anisotropy induced by pre-deformation [
16,
29,
30,
31]. The crystal plasticity constitutive model follows continuum mechanics and considers the slip mechanisms at grain-level by introducing resolved shear stresses along different slip orientations and the nonlinear hardening of back stress. Zhang Keshi et al. [
16,
17] adopted the Voronoi polycrystalline representative volume element to compute the subsequent yield surfaces of pure copper and aluminum under different proportional pre-loadings and tension-compression cyclic loading. The results showed that the anisotropic distortion of the yield surface can be characterized with the polycrystalline crystal plastic model. Due to the higher computational cost of polycrystalline analysis of subsequent yield behavior of metals, it is difficult to apply this approach directly in engineering structural analysis. A constitutive model presenting the anisotropic evolution of yield surface under the hypotheses of continuous and uniform deformation will be more accurate and efficient in engineering application.
The anisotropic distortion of the subsequent yield surface has an obvious difference in metals, including different work hardening aluminum alloy [
18,
32]. Taking the subsequent yield of T2 pure copper behavior into account, a macroscopic constitutive model reflecting the distortional yield surface induced by plastic deformation is put forward under the hypothesis of continuous uniform in material. The main research work includes the following: (1) a distortional hardening multiplier in conjunction with an anisotropic degradation factor are introduced into the isotropic hardening term of the yield function, and an anisotropic distortional yield surface model considering cross effect is proposed; (2) the yield surfaces simulated by the anisotropic plastic constitutive model are compared with the experimental results for T2 pure copper under pre-deformation, the simulations of the Chaboche cyclic plastic model and the crystal plasticity model. The model parameters are calibrated and the reasonability of the anisotropic model is verified.
2. Theoretical Model
The strain tensor can be decomposed into elastic
and plastic
parts under a small strain assumption as:
According to Hooke’s law, the relationship of stress and strain is expressed as:
where the fourth-order tensor
represents the material stiffness. The plastic constitutive model concentrates on the calculation of plastic deformation, and the incremental theory considers the influence of loading history on present deformation. The plastic strain rate
of classical plastic theory in regards to rate-independent model is given by:
where
is the accumulated equivalent plastic strain rate,
is the loading potential function, and a second-order tensor
is the plastic flow direction, which is calculated by the partial differential of the loading potential function according to Equation (22).
For a rate-dependent model, the equivalent plastic strain rate
and plastic strain rate
can be expressed as:
where the symbol < > is the MacCauley bracket,
are the parameters reflecting the material viscosity.
is the loading potential function, which is same as the yield function on the premise of the associative flow. The second-order tensor above
is the plastic flow direction, as noted above.
The yield function is the foundation in plastic hardening and flow analysis, thus a proper yield function should be developed and chosen to characterize most features of the subsequent yield surface due to the diversity yield behavior of the numerous materials. Based on the classical von Mises criterion, the yield function of mixed hardening mode is defined by:
where
,
are the deviatoric stress and back stress, respectively,
represents the yield stress corresponding to zero plastic strain, and
is the isotropic hardening function. For mixed hardening, the subsequent yield surfaces evolve as the expansion (i.e., isotropic hardening) and translation (i.e., kinematic hardening due to the evolution of back-stress) constitutive model.
The most common kinematic rule is the linear one due to Prager [
3], given as:
A practical way to describe a nonlinear kinematic hardening—giving rise to a correct modelization of cyclic loops—consists of introducing an evanescent memory of the plastic strain path, as initially proposed by Armstrong and Frederick [
4]. This nonlinear kinematic effect appears in the internal stress equation:
where the first term corresponds to the Prager’s linear rule and the second term to the evanescent strain memory, which will gradually weaken and disappear as the cumulative plastic strain increases; the cumulated plastic length is defined, for example, by:
Chaboche et al. [
8] proposed a very good continuous description of cyclic loading of materials, obeying a combination of isotropic and nonlinear kinematic rules by introducing many strain memory items. The back stress is decomposed into several individual variables, and the general form is:
where
denotes the number of back stress components. Considering the description ability and conciseness of the model,
is adopted in this paper.
and
are the material constants responsible for the evolution of the back stress.
The isotropic hardening function
of Equation (6) is given as:
where
is used to describe the saturation value of an isotropic hardening function,
is a parameter describing the relationship between yield surface radius and plastic strain, and
is the accumulated equivalent plastic strain.
Subjected to the initial conditions:
Vincent [
19] established a general cyclic plastic model by introducing 25 distortional variables to replace the classical back stress components. Considering the rotation of principal stress axes in the loading direction, this model is suitable for exhibiting the ratcheting behavior under multi-axial loading. Based on the contribution of Feigenbaum [
20], the distortional hardening model of metal was proposed by implementing the fourth-order anisotropic evolution tensor into the classical Hill orthogonal anisotropic yield criterion in thermodynamic framework. The yield function is defined as:
where
,
are deviatoric stress and back stress of yield function, respectively,
is the fourth-order unit stress deviatoric tensor, and
is the fourth-order anisotropic distortional hardening tensor.
Francois [
24] replaced the deviatoric stress tensor with the distortional stress tensor in the yield function of the classical mixed hardening mode, and the distortional yield function can be defined as:
where
,
are deviatoric stress and back stress of yield function, respectively,
represents the yield stress corresponding to zero plastic strain, and
is the isotropic hardening function.
is the distortion stress tensor calculated by deviatoric stress
, back stress
, and isotropic hardening function
.
Introducing the anisotropic distortional function
into the yield function of classical mixed hardening model by combining the isotropic hardening and the kinematic hardening, which characterize, respectively, the expansion and translation of subsequent yield surface, the anisotropic yield function can be expressed as:
where
,
are deviatoric stress and back stress of yield function.
represents the yield stress corresponding to zero plastic strain.
is the anisotropic distortional function determined by anisotropic distortional factor
and the anisotropic degradation factor
. In this presentation, the nonlinear kinematic hardening rate of yield surface, which determines the evolution of back stress (center of yield surface), is described as Equation (10), and the isotropic hardening of yield surface
is given by Equation (11). The values for the corresponding coefficients
and
are indicated in
Table 1.
The anisotropic distortional function
in Equation (15) can be summarized as:
where
is the back stress during reloading, and it is constant when the stress state is within the yield surface. The material yields with the increase of reloading load and back stress
evolves with the following Equation (10). The fourth-order tensor
is the back stress corresponding to the unloading stress followed by preloading in subsequent yield surface test, which is related to preloading path and preloading load.
remains the constant during elastic deformation.
where
is the anisotropic degradation factor, which describes the evolution from anisotropy to quasi-isotropy with regard to yield function. It can be calibrated by the evolution test data of yield surface under different offset strains. If the anisotropic degradation factor
is very big, the subsequent yield surface presents only the proportional expansion of the initial anisotropic yield surface, i.e., the anisotropy of yield keeps stable. On the other hand, the initial anisotropic yield surface evolves toward the isotropic yield surface when the anisotropic degradation factor
is a smaller value, i.e., the anisotropy of yield degrades. An appropriate anisotropic degradation factor
is attempted based on the comparison of the evolution progress between the subsequent yield surface simulated and the experimental results for different given offset strains (confer the data in
Figure 3). The anisotropic distortional function
varies from
to 1 with the increase of reloading load. When the back stress reaches a certain critical value, which means
and
, the anisotropic distortional yield surface evolves into a yield surface without distortion.
The shape of subsequent yield surface under reloading tends to be circular in
stress space with the increase of the offset strain [
15,
17]. The subsequent yield surfaces of T2 pure copper under different offset strains [
16] are shown in
Figure 3. The results show that the anisotropy of the subsequent yield surfaces degrades gradually, and the subsequent yield surfaces approach a circle.
Noting the cross effect and anisotropic characteristics of metal subsequent yield surface test in
Figure 3, the anisotropic distortional factor
in anisotropic distortional function
can be given as:
where
are the parameters of the cross effect, which can be calibrated by the characteristic of distortional subsequent yield surface corresponding to the small offset strain under tension and torsion loading (cf.
Section 3.4). Anisotropic distortional yield surface degenerates into the circular mixed hardening yield surface without cross effect in the condition of
.
is the angle between the preloading direction and the reloading direction, and it is defined by the deviatoric stress and back stress:
Defining the two flow tensors as:
Substituting Equation (20) into (19), Equation (19) can be simplified as:
Following the orthogonal flow rule, the plastic flow direction of anisotropic yield surface can be described by the partial derivative of Equation (15):
Because H is 1 and
under preloading and unloading [given in Equation (16)], the plastic flow direction is orthogonal to the cylindrical yield surface during the two stages, i.e.,
. Similarly,
and
are established when anisotropy disappears in the process of reloading, i.e.,
. The yield surface evolves from anisotropy to quasi-isotropy gradually in the condition of
under reloading, as displayed in
Figure 3. Combining Equations (16), (20), and (22), the plastic flow direction
can be expressed as:
Based on the definition of anisotropic distortional factor
in Equation (18),
can be obtained by:
Due to the definition of
in Equation (20), which is irrelevant to
, the expression can be deduced as:
Integrating Equations (23), (24), and (26), the plastic flow direction
can be summarized as:
where the stress tensors
and
are normalized according to back stress tensor
and stress tensor
in Equation (20).
is a fourth-order tensor which can be defined as:
4. Simulation of Anisotropic Subsequent Yield Surface
The subsequent yield surface of T2 pure copper under pretension loading exhibited significant anisotropy [
34], as shown in
Figure 10. According to the test method described in
Section 3, the initial yield surface with the offset strain of 5 × 10
−5 was nearly circular when the material yielded firstly under pretension loading [cf. the point A in
Figure 10a,b]. T2 pure copper began to harden when the equivalent pre-tensile strains increased to 1% and 5% [cf. the point B and point C in
Figure 10a,b], and the subsequent yield surface [cf. the blue triangle and red rectangular in
Figure 10a] was larger than the initial yield surface [cf. the black circle in
Figure 10a]. The subsequent yield surface was characterized by a region of high curvature along the direction of preloading and a region of flattening on the opposite side. The yield surface calculated by the anisotropic model in this paper [cf. the blue, red, and black curves in
Figure 10a] were in good agreement with the experimental results, which indicates that this model can describe the anisotropic distortional characteristics of the subsequent yield surface.
The simulation results of the anisotropic distortional yield surface model were compared with the experimental data and the simulation results of polycrystalline plastic model [
16] to discuss the subsequent yield and hardening process of metal materials. In particular, the yield stresses determined by the very small offset strain showed significant dispersion, which gradually decreased as the offset strain increased. Thus the typical offset strains of 1 × 10
−4, 5 × 10
−4, and 1 × 10
−3 were adopted to investigate the subsequent yield behavior.
The subsequent yield surfaces determined by unloading and proportional reloading after three directions of pre-deformations (tensile, torsion, and tension-torsion) are shown in
Figure 11,
Figure 12 and
Figure 13 respectively, which includes yield stresses corresponding to the typical offset strains of 1 × 10
−4 [
Figure 11a,
Figure 12a and
Figure 13a], 5 × 10
−4 [
Figure 11b,
Figure 12b and
Figure 13b], 1 × 10
−3 [
Figure 11c,
Figure 12c and
Figure 13c]. Since the classical Chaboche plastic constitutive model only includes kinematic hardening and isotropic hardening, the subsequent yield surface calculated by Chaboche model without distortional hardening is a circle with the center of back stress, as shown in
Figure 11,
Figure 12 and
Figure 13.
The anisotropic distortional yield surface expanded continuously with the increase of offset strain, which was consistent with the experimental data and the simulation results of the crystal plastic model. The cross effect of the anisotropic yield surface was relatively prominent with the offset strain of 1 × 10
−4. The “sharp point” appeared in the preloading direction, and the yield surface was flat in the opposite direction. The cross effect of the distortional yield surface was gradually weakened when the offset strain increased from 1 × 10
−4 to 5 × 10
−4 and eventually reached 1 × 10
−3, and the anisotropic distortional yield surface was gradually approaching the isotropic circular yield surface computed by the Chaboche model (cf. the purple solid lines in
Figure 11,
Figure 12 and
Figure 13), which complied with the experimental data.
Comparing the results of the crystal plastic constitutive model (cf. the green lines in
Figure 11,
Figure 12 and
Figure 13) with the results of anisotropic distortional yield surface model shows that both of them could describe the expansion, shift, and distortion of subsequent yield surface, but the results of the anisotropic model in this paper were generally better than those of the crystal plastic constitutive model. The simulation results of the anisotropic model were very close to the experimental data under pre-tension, which was superior to the simulation results of crystal plasticity model, as presented in
Figure 11. The results of the anisotropic model were slightly higher than the experimental values in the case of pre-torsion and pre-tension-torsion, and the crystal plastic model was more accurate along the torsional direction. However, the yield stress perpendicular to the preloading direction for the crystal plastic model in
Figure 12 and
Figure 13 was not as precise as that of the anisotropic model.
The parameters of kinematic hardening and isotropic hardening parts of anisotropic yield model in
Table 1 are the same as those in Chaboche model. The diameter of the anisotropic distortional yield surface (cf. the red solid lines in
Figure 11,
Figure 12 and
Figure 13) was smaller than that of the circular yield surface simulated by Chaboche model under small offset strains, and the ratio of the two diameters along the loading direction was
, as shown in
Figure 11,
Figure 12 and
Figure 13. Therefore, the subsequent yield surface derived from the Chaboche model was the same as the subsequent yield surface simulated by the anisotropic yield model when anisotropy gradually degenerated to isotropy. The diameter of the circular yield surface simulated by Chaboche model could be reduced to fit the experimental curve by adjusting the parameters of the model, and the simulation results were obviously better than the curves with the parameters adopted in this paper. However, the non-circular characteristic of the experimental yield surface could not be well captured by Chaboche model.
Due to the different numbers of the experimental data of the three pre-loading conditions, the error of each yield surface determined by different offset strain under various pre-deformation conditions is defined based on the benchmark of the experimental yield stresses as:
where the
and
represent, respectively, the mean errors of the anisotropic model and the crystal plastic model with the i-th offset strain and the j-th pre-deformation type, the N is the total number of reloading directions, the k indicates the k-th reloading direction, and
expresses the 2-norm (Euclid norm) of the stress tensor.
The mean value of every error for the anisotropic model and the crystal plastic model are given respectively by:
where the
and
are the mean error of the anisotropic model and the crystal plastic model, and the total numbers of i and j are both three.
The standard deviation is calculated with fundamental statistics knowledge by the error and the mean error. The errors and the standard deviation of the yield surface with the anisotropic model and the crystal plastic model under every offset strain and pre-formation are listed in
Table 2. There was an obvious difference in the predicted accuracy between various offset strains and pre-formations conditions. In terms of the pre-torsion loading, the mean error and the standard deviation of the anisotropic model were higher than those of the crystal plastic model, which was attributed to the shear hardening model in slip process. On the whole, however, the mean error and the standard deviation of the anisotropic model, 11.54% and 8.23%, were both smaller than those of the crystal plastic model, as shown in the last column of
Table 2.