Next Article in Journal
3D/4D Printing Application for Shape Memory Materials
Previous Article in Journal
New PEDOT Derivatives Electrocoated on Silicon Nanowires Protected with ALD Nanometric Alumina for Ultrastable Microsupercapacitors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Constitutive Modeling of AA7050-T7451 with Anisotropic Stress Transformation

Department of Mechanical Engineering, University of Jinan, Jinan 250022, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(17), 5998; https://doi.org/10.3390/ma15175998
Submission received: 10 August 2022 / Revised: 23 August 2022 / Accepted: 26 August 2022 / Published: 30 August 2022

Abstract

:
The mechanical properties of anisotropic materials are generally characterized based on the orthotropy or transverse isotropy. However, the two-dimensional plane stress problems cannot comprehensively characterize the anisotropy of materials. In this study, based on the theory of elasticity and the transformation of the three-dimensional space coordinate system, combined with the projection relationship of the Cauchy stress tensor of an arbitrary section, the transformation relationship of the elastic modulus, shear modulus, and stress–strain between the orthogonal and load coordinate systems are obtained. The orthotropic Johnson-Cook (JC) constitutive model of AA7050-T7451 aluminum alloy is modified by fitting, and the constitutive relationship at any spatial angle is theoretically calculated by combining the obtained spatial coordinate transformation matrix. The generated spatial constitutive model is verified and modified through experiments. The results demonstrate that the theoretical mechanical properties and the modified spatial constitutive model can accurately reflect the effect of the spatial angle on the material stress distribution.

1. Introduction

Aluminum alloys are the most extensively used materials in industrial and civil applications and weapon manufacture, mainly due to their excellent properties, including high strength-to-weight ratio, low density, high-temperature resistance, and excellent corrosion resistance [1,2]. AA7050-T7451 aluminum alloy, which includes certain anisotropy because it is a typical poly crystalline alloy material, is typically used in aerospace applications for critical aircraft components [3,4,5]. Under different material pretreatment conditions, due to the different existence modes of micro-atoms (such as hydrogen) in the material, the grains in the material will produce different micro-effects, such as dislocation, twinning and other phenomena [6], so that the micro-grain of the material has complicated changes. Finally, the materials show different mechanical properties. Constitutive models have been effectively employed to characterize the mechanical properties of materials during plastic deformation [7,8]. Constitutive material modeling directly affects the accuracy of the theoretical study and finite element analyses [9]. Therefore, it is crucial to propose valid constitutive formulations that can describe the anisotropic characteristics for predicting the mechanical properties of materials.
The differences in the mechanical properties of anisotropic materials are generally manifested by the elastic modulus, Poisson’s ratio, and strength, in different directions [10,11,12]. Therefore, classical elastoplastic mechanics based on homogeneous, continuous, and isotropic assumptions cannot correctly reflect anisotropic mechanical behavior. The elastoplastic constitutive strength theory and the numerical analysis of anisotropic materials have been extensively researched. Plunkett et al. [13] proposed yield functions describing the anisotropic behavior of orthotropic sheet metal under tension and compression. Oana et al. [14] established a new isotropic criterion that well describes the yielding of randomly oriented face-centered polycrystalline metallic materials and extended this criterion to describe the orthotropy using generalized stress invariants. To describe the plastic anisotropy of metallic alloys, Cardoso et al. [15] developed a yield function that delivered an accurate prediction of directional stresses. Kazimierz et al. [16] studied a mathematical method for describing orthotropic material cutting based on the fracture theory and included the work of separation (fracture toughness) in addition to the material plasticity and friction. Meng et al. [17] established a modified JC constitutive model considering the forming direction, based on the dynamic impact shear test of AA7050-T7451 aluminum alloy in the typical forming direction (rolling direction (RD), transverse direction (TD), and normal direction (ND)).
The above studies were based on the orthotropic characteristics (RD, TD, ND) of materials because of the complexity of anisotropic mechanical behavior. In order to further study the anisotropy of materials, some scholars characterized the anisotropy of different orientations in a specific plane. Zhang et al. [18] accurately described the tensile stress–strain curves of anisotropic (along a certain angle to the rolling direction) aluminum alloy sheets using the three-parameter Barlat yield criterion. Aboubakr et al. [19] investigated the effect of the rolling orientation on the mechanical, in-plane anisotropy, and gamma ray shielding properties of Al-Cu-Li-Mg-X alloy. Based on the three-parameter Barlat yield criterion, Liu et al. [20] proposed a new five-parameter yield criteria in which two transverse shear stresses were added.
Although extensive research has been conducted on the anisotropy of materials, and yield criteria for describing their mechanical behavior have been proposed, the spatial conditions of orthotropy and transverse isotropy are yet to be resolved. The two-dimensional plane stress problem cannot comprehensively characterize the anisotropy of materials because three-dimensional characterization of the anisotropy cannot be realized.
Based on the theory of elasticity and the conversion of the space coordinate system, this study aims to obtain a theoretical method for comprehensively characterizing the anisotropy by introducing three-dimensional space angle variables. Quasi-static loading tests and dynamic impact compression/shear tests are performed on specimens in the typical forming direction (ND, TD, and RD) and space angle (RD-30°-ND, RD-45°-ND, RD-60°-ND). The quasi-static and high-strain-rate mechanical responses at different rates, directions, and temperatures are analyzed, and the proposed theoretical method is verified. Finally, the least squares method is used for fitting and modifying the constitutive material parameters, and a JC constitutive model that includes the spatial angle is established by combining the coordinate transformation matrix. The modified material constitutive model will be applied to the theoretical research and finite element simulation of material mechanical properties.

2. Theoretical Research on Anisotropy

2.1. Coordinate-System Conversion Matrix

The main form of metal plastic deformation is the crystal slip, which proceeds along a certain crystal plane and direction [21]. The strongest interatomic binding force has the least probability of disorder in the arrangement of atoms; hence, crystal slip preferentially occurs in the crystal plane with the highest atomic density and the crystal orientation with the densest atomic arrangement [22]. According to plastic mechanics, until the stress acting on the slip surface and direction reaches the critical value, the slip will proceed along a specific direction. AA7050-T7451 aluminum alloy is a typical face-centered cubic lattice, which consists of four groups of preferred slip surfaces [23,24,25,26]. As there are three preferred sliding directions for each slip surface, AA7050-T7451 has 12 slip systems and 384 slip states [27,28]. Based on the above analysis, it can be concluded that the simple two-dimensional plane stress problem cannot accurately characterize the anisotropy. In order to comprehensively characterize material anisotropy, a constitutive model that can characterize the mechanical parameters in each direction of the three-dimensional space must be established.
In the mechanical properties of anisotropic materials, the differences are mainly manifested by the elastic modulus, shear modulus, and stress. When anisotropic materials are subjected to non-orthotropic loads, the assumptions of orthotropy and transverse isotropy can no longer satisfy the descriptions of their stress parameters. Therefore, a load coordinate system must be established to characterize the mechanical parameters, under the load direction. An orthogonal (material) coordinate system was established with the x, y, and z axes of the TD, RD, and ND, respectively. A load coordinate system was created by considering the load direction as the z’ axis; thus, the three axes of it are not along the elastic direction (TD, RD and ND). The positional relationship between the axes of the orthogonal and load coordinate systems is shown in Figure 1; α, β, and γ represent the angle between the x, y, and z axes of the load coordinate system (oxyz′) and orthogonal coordinate system (oxyz), respectively, while li, mi, and ni represent the directional cosine between the coordinate axes. To characterize the mechanical properties in any direction of the three-dimensional space, the transformation relationship between the load and material coordinate systems must be obtained, and the orthotropic mechanical parameters in the material coordinate system and coordinate transformation matrix must be combined [29].
According to the theory of elasticity, the stress vector Fn of any oblique section can be uniquely determined by nine stress components in the elastic body, also known as the Cauchy stress tensor [30]. Figure 2 illustrates the stress components (positive stress components σxx, σyy and σzz, and the shear stress components) of the Cauchy tensor that completely describe the stress state at an arbitrary point in an elastic body. The normal stress component represents the tensile or compressive stress on the corresponding section. As per the equivalent law of the shearing stress, the shear stress can be expressed by τxy, τyz and τzx [31,32]. The expression for the stress tensor σ is as follows:
σ = σ x x σ x y σ x z σ y x σ y y σ y z σ z x σ z y σ z z = σ x x τ x y τ z x τ x y σ y y τ y z τ z x τ y z σ z z
Therefore, the stress state of any cross-section can be characterized by the stress parameters in the orthogonal state. The intermediate variable is given by the positional relationship between the cross-section and orthogonal coordinate system, i.e., the transformation relationship between the load and orthogonal coordinate systems. Using the above transformation matrix, the stress parameters obtained in the orthogonal state can be transformed into the stress state of an arbitrary cross-section, which can be used to characterize the mechanical parameters in three-dimensional space [33,34,35].
Considering the normal direction of an oblique section ABC as the z′ axis of the new coordinate system, the stress components Xv, Yv, and Zv on the section along the original coordinate axis can be expressed as follows:
X v = σ x l 3 + τ x y m 3 + τ x z n 3 Y v = τ y x + σ y m 3 + τ y z n 3 Z v = τ z x l 3 + τ z y m 3 + σ z n 3
Projecting Xv, Yv, and Zv onto the z′ axis, the normal stress σ z of ABC can be obtained as follows:
σ z = X v l 3 + Y v m 3 + Z v n 3 = σ x l 3 2 + σ y m 3 2 + σ z n 3 2 + 2 τ x y l 3 m 3 + 2 τ y z m 3 n 3 + 2 τ z x n 3 l 3
When the tetrahedral differential body OABC tends to be infinitely small, the x′ and y′ axes are in the oblique section, ABC. Therefore, the projection of Xv, Yv, and Zv on the x′ and y′ axes indicates the shear stress of the oblique differential body section.
τ x z = τ z x = X v l 1 + Y v m 1 + Z v n 1 = σ x l 1 l 3 + σ y m 1 m 3 + σ z n 1 n 3 + τ x y l 3 m 1 + l 1 m 3 + τ y z m 1 n 3 + m 3 n 1 + τ z x l 1 n 3 + l 3 n 1
τ y z = τ z y = X v l 12 + Y v m 2 + Z v n 2 = σ x l 2 l 3 + σ y m 2 m 3 + σ z n 2 n 3 + τ x y l 2 m 3 + l 3 m 2 + τ y z m 2 n 3 + m 3 n 2 + τ z x l 2 n 3 + l 3 n 2
The stress tensor of the oblique section, ABC, in the new coordinate system, oxyz′, can be expressed as
σ = σ x τ x y τ x z τ y x σ y τ y z τ z x τ z y σ z
The transformation matrix of the stress tensor from the original coordinate system oxyz to the new coordinate system oxyz′ can be obtained by combining Equations (2)–(6).
σ = A σ A T
where [A] = l 1 m 1 n 1 l 2 m 2 n 2 l 3 m 3 n 3 is the coordinate transformation matrix.
Therefore, the stress in the load coordinate system can be obtained by combining the mechanical parameters of the orthogonal coordinate system and the coordinate transformation matrix. Thus, the breakthrough from orthotropy to anisotropy characterization in three-dimensional space can be realized.

2.2. Mechanical Properties of the Anisotropy

The transformation of the stress state between the coordinate systems can be realized using the transformation matrix between the load and the orthogonal coordinate systems. As previously mentioned, the difference in the mechanical characteristics of anisotropic materials is mainly manifested by the elastic modulus, shear modulus, and stress–strain parameters; therefore, the transformation relationships of the strain, elastic modulus, and shear modulus between the orthogonal coordinate system and load coordinate system can be further analyzed using the transformation matrix.
In the material coordinate system oxyz, the stress–strain relationship of anisotropic materials can be expressed as follows:
ε = S σ
According to the theory of elasticity and the transformation relationship of the space coordinate system, the stress–strain relationship in the load coordinate system can be obtained as follows:
ε = T S σ = T S T T σ = S σ
where   ε   is the strain vector in the load coordinate system, T is the flexibility conversion matrix, S is the anisotropic flexibility matrix, σ   is the stress vector in the material coordinate system, σ   is the stress vector in the load coordinate system, S = T S T T is the equivalent anisotropic flexibility matrix in the load coordinate system, T and   S are expressed by Equations (10) and (11), respectively.
T   =   l 1 2 m 1 2 n 1 2 l 1 m 1 m 1 n 1 n 1 l 1 l 2 2 m 2 2 n 2 2 l 2 m 2 m 2 n 2 n 2 l 2 l 3 2 m 3 2 n 3 2 l 3 m 3 m 3 n 3 n 3 l 3 2 l 1 l 2 2 m 1 m 2 2 n 1 n 2 l 1 m 2   +   l 2 m 1 m 1 n 2   +   m 2 n 1 n 1 l 2   +   n 2 l 1 2 l 2 l 3 2 m 2 m 3 2 n 2 n 3 l 2 m 3   +   l 3 m 2 m 2 n 3   +   m 3 n 2 n 2 l 3   +   n 3 l 2 2 l 2 l 3 2 m 3 m 1 2 n 3 n 1 l 3 m 1   +   l 1 m 3 m 3 n 1   +   m 1 n 3 n 3 l 1   +   n 1 l 3
S = s 11 s 12 s 13 s 14 s 15 s 16 s 22 s 23 s 24 s 25 s 26 s 33 s 34 s 35 s 36 s 44 s 45 s 46 s 55 s 56 s 66
The elastic modulus E and shear elastic modulus G of the tetragonal system in an arbitrary direction can be obtained as follows [36]:
E 1 = R 1 4 + R 2 4 S 11 + R 3 4 S 33 + R 1 2 + R 2 4 R 3 2 2 S 13 + S 44 + R 1 2 R 2 2 2 S 12 + S 66
G 1 = 2 S 11 R 1 1 R 1 2 + R 2 2 1 R 2 2 + 2 S 33 R 3 2 1 R 3 2 4 S 12 R 1 2 R 2 2 4 S 13 R 1 2 R 3 2 + R 2 2 R 3 2 + 1 2 S 44 2 R 1 2 4 R 2 2 R 3 2 R 2 2 4 R 1 2 R 3 2 + 1 2 S 66 1 R 3 2 4 R 1 2 R 2 2
where Ri is the direction cosine of the load coordinate system relative to the orthogonal coordinate system, and Sij is the flexibility coefficient of the material.
The experimental data (elastic modulus, Poisson’s ratio, shear modulus, and stress–strain) in different directions of anisotropic material specimens were obtained through correlative experiments. The theoretical values corresponding to the load coordinate system, which can be experimentally verified, were obtained by combining the transformation matrix. Subsequently, the constitutive equation that can comprehensively characterize the anisotropic mechanical characteristics was stablished and corrected.

3. Experimental

3.1. Specimen Preparation

Pre-stretched AA7050-T7451 sheets were used in this study. In order to obtain the mechanical parameters of the material coordinate system, the specimens were sampled along the typical forming directions (RD, TD, and ND) of the sheet. A specific angle plane was selected for sampling to verify the theoretically derived results; the sampling direction is depicted in Figure 3. In Figure 3a, the specimen directions for impact compression and shear are respectively defined along their loading directions, namely, the axis direction of the cylinder and the thickness direction of the shear specimen. As shown in Figure 3b, the angle plane was selected based on the TD-RD plane at a certain angle γ (γ = 0°, 30°, 45°, 60°, 90°), which was defined as the RD-γ-TD plane. Cylindrical specimens were processed by wire-cutting in the RD-γ-TD plane; specimens at 0° and 90° constituted the ND and RD samples, respectively. Quasi-static compression tests were carried out at 25 °C (room temperature), and specimens sized Ø10 × 20 mm were sampled in the RD, TD, and ND, respectively. Impact compression and impact shear tests were performed using samples sized Ø5 × 4 mm and 40 × 23.5 ×2 mm (L × D × H), respectively.

3.2. Experimental Setup and Principle

Dynamic impact compression and shear tests were performed using a split Hopkinson pressure bar (SHPB) device, which is a typical method for obtaining the dynamic mechanical properties [37]. The elastic incident and reflected waves were measured by an input bar strain gauge. Further, the elastic incident wave was measured by an output bar strain gauge and was amplified by a dynamic strain gauge. The waveform data were collected and stored using a digital data recorder. Finally, the dynamic–mechanical property data were obtained according to the one-dimensional stress wave theory [38,39]. As the ultrahigh strain rate in the actual cutting process cannot be achieved by the SHPB device during the impact shearing process, a dynamic impact shearing “bar-tube” device was used to enhance the SHPB. The SHPB test equipment and “bar-tube” device are depicted in Figure 4. An intelligent temperature-controlled electric heating furnace was used for heating the specimens at a maximum temperature of 800 °C and temperature control precision of ±1 °C.
At room temperature (25 °C), quasistatic compression tests were carried out on specimens (Ø10 × 20 mm) in different directions (RD, ND, TD) at a strain rate of 1.6 × 10−3 s−1. Dynamic impact shear tests at strain rates of 1.0 × 104, 2.0 × 104, and 3.0 × 104 s−1 were conducted on shear specimens with different orientations (RD, ND, TD). Further, dynamic impact compression tests were performed on cylindrical specimens (Ø5 × 4 mm) in different directions (RD, ND, TD, and RD-γ-TD) at strain rates of 0.1 × 104, 0.2 × 104, 0.4 × 104, and 0.6 × 104 s−1, respectively. In order to obtain the mechanical properties at high temperature, impact compression tests were conducted at a strain rate of 0.4 × 104 s−1, at 100, 200, 300, 350, 400, and 450 °C, respectively. Each group of experiments was repeated thrice, and the experimental results were averaged to ensure accuracy.

4. Results and Discussion

4.1. Establishment of a Anisotropic Constitutive Model

As displayed in Figure 5, based on the selection of the angle plane in this study, the load coordinate system under angle γ was established as per Figure 1 and Figure 3b. According to the positional relationship between the orthogonal and load coordinate systems, the cosine of each direction was obtained and substituted into the coordinate system transformation matrix [A]. As depicted in Equation (14), the stress relationship expression between the load and orthogonal coordinate systems can be obtained by substituting the transformation matrix [A] of the obtained coordinate system in Equation (7).
σ = σ x x τ x y c o s γ τ x y s i n γ + τ x y c o s γ τ y x c o s γ σ y y c o s 2 γ σ y s i n γ c o s γ + τ y z c o s γ τ y x s i n γ + τ x z c o s γ σ y y + τ y z c o s γ c o s γ σ y y s i n 2 γ + σ z z c o s 2 γ + 2 τ z y c o s γ s i n γ
The normal stress at 0° (ND) and 90° (RD) is σzz and σyy, respectively, in the above equation. The stress in the orthogonal coordinate system is directly related to the accuracy of the theoretical stress value in the load coordinate system. Therefore, the constitutive equation in the orthogonal coordinate system, which can accurately express the stress–strain relationship of each direction in this system, must be established and modified.
The constitutive equation of a material describes the stress–strain relationship of the material under different loading conditions [40,41]. The empirical Johnson–Cook constitutive model was applied to describe the stress–strain relationship of ferrous and nonferrous metals, under conditions of high-strain rate and considerable deformation. In addition, the structure of the JC constitutive model is simple and convenient; it introduces the strain strengthening, strain-rate strengthening, and thermal softening parameters of the material [42], and can be expressed as follows:
σ = A + B ε n 1 + C ln ε ˙ 1 T T r T m T r m
where A, B, C, n, and m represent the initial yield stress, strain hardening modulus, hardening index, strain-rate sensitivity coefficient, and thermal softening coefficient, respectively; T, Tr, and Tm represent the deformation temperature, room temperature, and melting temperature of the material, respectively; σ, ε, and ε ˙ are the stress, strain, and reference strain rate, respectively [43].
For anisotropic materials, the constitutive model should be constructed in the typical forming direction and fitted in the orthogonal coordinate system. If the modified JC constitutive model can accurately characterize the stress–strain relationship of the material in the orthogonal coordinate system, the stress–strain constitutive relationship in the load coordinate system at an arbitrary angle γ can be obtained by combining Equations (14) and (15):
σ γ = σ R D s i n 2 γ + σ N D c o s 2 γ + 2 τ c o s γ s i n γ

4.2. Anisotropic Mechanical Property Analysis of AA7075-T7451

In order to modify the constitutive model and verify the anisotropy in three-dimensional space, dynamic impact compression tests were performed at different strain rates and temperatures for the specimens, in different directions. The stress–strain curves of the dynamic impact compression tests in the 0.1 × 104–0.6 × 104 s−1 strain-rate range were obtained, as depicted in Figure 6. When the strain was less than 0.02, the material was in the elastic deformation stage, and the stress increased rapidly because of the minor deformation caused by work hardening. With the increase in strain, thermal softening effect occurred in the material, and along with work hardening caused a significant decrease in the stress rate. When the work hardening and thermal softening effect gradually attained equilibrium, the stress exhibited a stable change stage. With the increase in strain rate, the maximum corresponding strain increased from 0.08 to 0.35. Meanwhile, the maximum stress in the stable stage also exhibited an upward trend, from 653 to 847 MPa.
At the same strain rate, the stresses at different angles γ, differed in the elastic, plastic deformation, and stable stages. As the material stress fluctuation was less in the stable deformation stage, the stress corresponding to the initial strain in the stable stage was selected [44]. The influence of angle γ on the stress was analyzed by comparing the stability stress under different loading conditions. Comparing the stability stresses at different strain rates, the variation trend of the stress with angle γ was as follows: The stress decreased with the increase in angle γ and reached a minimum when γ = 45°. As γ continued to increase, the stress increased after reaching a minimum value. Therefore, in the 0–90° range, the stress initially decreased and then increased with γ.
Figure 7 illustrates the stress–strain curves of AA7050-T7451 at different temperatures, for a strain rate of 0.4 × 104 s−1. When the temperature increased from room temperature to 200 °C, the height of the specimen decreased by only 0.72 mm, and the stress decreased from 755 to 733 MPa; therefore, the effect of the temperature on the stress of the material was less when the temperature was below 200 °C. When the temperature reached 300 °C, the stress decreased sharply, and the flow stress was only 60–70% of that at 200 °C. The specimen height decreased by 0.85 mm, compared to that at 200 °C, and the specimens exhibited obvious damage. The thermal activation increases the atomic kinetic energy in the alloy, when the deformation temperature is increased to a certain value, resulting the rapid softening. Further, the critical shear stress of the alloy slip systems decreases, and the dislocation slip can easily overcome the hindrance effect at lower stress [45]. Therefore, the flow stress of AA7050-T7451 decreases with the increase in temperature.
Different from the mechanical curves under different strain rates, the initial strain of the stable deformation stage at different temperatures was approximately 0.2; therefore, the stress value at a strain of 0.2 was selected as the stability stress (σ0.2) to analyze the variation trend, under different deformation temperatures. From Figure 7b, it can be observed that the variation in stress with temperature is consistent in the RD direction. When the temperature was below 300 °C, the stress of 0–90° specimens initially decreased and then increased, at the same temperature. The stress of the 45° specimen was the least, consistent with that at room temperature. Thus, the steady-state stress in each direction exhibited obvious anisotropy during the dominant stage of work hardening. Loading the materials at angles of 30°, 45°, and 60° reduced the stress, particularly, the minimum stress required for the 45° specimen deformation. The main reason is because the grain slip system in the material is anisotropic, and the preferred slip system direction is that with the highest atomic density or densest atomic arrangement. Therefore, it can be concluded that AA7050-T7451 material has the densest atomic arrangement in the direction of 45°, resulting in preferential slip orientation in this direction. The thermal softening effect was dominant at temperatures beyond 300 °C. With the increase in deformation temperature, the kinetic energy of the atoms in the material increased, and the crystal slip directions were no longer obvious.
In conclusion, the stress of AA7050-T7451 material is not only related to the strain rate, strain, and temperature, but also to the angle of its forming direction. Therefore, it is necessary to construct and modify the constitutive equation to include the spatial angle function, for describing the mechanical properties of the material.

4.3. Modification of the Constitutive Model

4.3.1. Constitutive Model in the Orthogonal Coordinate System

In order to introduce the spatial angle function and obtain the constitutive model of the anisotropy in three-dimensional space, the stress–strain curves in the RD, TD, and ND were initially fitted and modified. The constitutive equation that can accurately characterize the stress relationship in the orthogonal coordinate system was obtained. The parameters of the JC model A, B, C, and n were fitted and modified using the experimental stress–strain curve along with least squares method, as shown in Table 1.
As the single thermal softening coefficient m cannot satisfy the mechanical characteristic curve of AA7050-T7451 at high temperature, the thermal softening effect term Kt was polynomial fitted for accurately describing the rheological properties of the material at high temperature [46].
K t = σ / A + B ε n 1 + C ln ε ˙
Using the stress–strain curve of the high-temperature impact compression test, Kt at different temperatures can be obtained, i.e., Kt can be calculated by substituting the stress with a strain rate of 0.4 × 104 s−1 and strain of 0.1 at different temperatures in Equation (16). If the number of polynomials is high, the expression will be complicated, and if the number of times is less, the fitting effect will be poor. Therefore, fifth-degree polynomial fitting with a simple expression and good fitting effect was used to fit the result, as shown in Figure 8. When the temperature was below 300 °C, the thermal softening coefficient decreased gradually with the increase in temperature, decreasing by approximately 0.13 only; in addition, the thermal softening coefficients in the TD, RD, and ND directions were considerably different. Sharp softening occurred at 300 °C, and the thermal softening coefficient dropped sharply to approximately 0.45; although the thermal softening coefficients in all the directions were different, the difference was not significant. After 300 °C, the thermal softening coefficient exhibited a small downward trend with the change in temperature, decreasing by 0.05 only, for a temperature increase of approximately 150 °C. Thus, the fitting curves can better reflect the mechanical properties of AA7050-T7451 material at high temperature.
Therefore, by fitting and modifying the parameters of the JC constitutive model, a constitutive equation that can accurately describe the mechanical properties of AA7050-T7451 was obtained, under different strain rates, temperatures, and forming directions (ND, RD, and TD):
σ N D = 365 + 516 ε 0.27 1 + 0.042 ln ε ˙ K t N D σ R D = 342 + 452 ε 0.26 1 + 0.029 ln ε ˙ K t R D σ T D = 357 + 490 ε 0.3 1 + 0.04 ln ε ˙ K t T D
K t = N D : 0.9 0.002 t + 3.23 × 10 5 t 2 2 × 10 7 t 3 + 5.9 × 10 10 t 4 5.1 × 10 13 t 5 R D : 1.2 0.016 t + 2.22 × 10 4 t 2 1.22 × 10 6 t 3 + 2.8 × 10 9 t 4 2.3 × 10 12 t 5 T D : 0.98 0.006 t + 9.13 × 10 5 t 2 5.43 × 10 7 t 3 + 1.29 × 10 10 t 4 1.07 × 10 13 t 5

4.3.2. Verification of the Elastic Modulus and Shear Modulus

In order to verify the feasibility of the theoretical method, the elastic modulus, shear modulus, and Poisson’s ratio of the typical forming direction (orthogonal coordinate system) were obtained, through the dynamic impact compression and shear tests. The angles between the orthogonal and load coordinate systems were α = 0, β = 0, and γ. The anisotropic flexibility matrix of AA7050-T7451 aluminum alloy material was calculated as follows:
S = 0.0135 0.0128 0.0106 0.0152 0.0143 0.01187 0.0171 0.0162 0.0134 0.04115 0.0455 0.0446
According to the flexibility matrix [S] of AA7050-T7451 and combining Equations (12) and (13), the change trends of the theoretical values of the elastic modulus and shear modulus with the change in γ were obtained, as shown in Figure 9. As γ increased, the theoretical elastic modulus and shear modulus values initially decreased and then increased, and the minimum value appeared at 45°. This trend was approximately the same as that of the stress change in the impact compression test, shown in Figure 6 and Figure 7. The change trends of the theoretical values of the elastic modulus and shear modulus were the same as that of the experimental values, verifying the accuracy of the theory and the method for characterizing the three-dimensional space anisotropy.

4.3.3. Modification of the Spatial Constitutive Model

By comparing the experimental and theoretical values of the elastic modulus and shear modulus, respectively, the accuracy of the theoretical method for characterizing the three-dimensional mechanical parameters was preliminarily verified. To completely characterize the anisotropy of materials in three dimensions, the stress state should be characterized and verified to establish a constitutive relationship containing the spatial angle. In order to further verify the accuracy of the theoretical method, the theoretical value was verified using the results of the dynamic impact compression test. The stress obtained by dynamic impact compression is the normal stress σ z , in the load coordinate system. The stresses in the stable stage depicted in the stress–strain curves at different strain rates in different directions were selected and verified. The verification results are presented in Figure 10.
In Figure 10, the verification results demonstrate that under different strain rates and temperature conditions, the change trend of the obtained theoretical values under different loading conditions is consistent with the test results. A certain gap existed between the theoretical value of the stress and the experimental value at various angles. Although the error was large, the ratio between the theoretical and experimental values was within a range of 0.85–0.90 because the theoretically calculated stress values did not consider the effects of the material preforming processes on the anisotropy. Therefore, a correction factor k needs to be included, which accurately characterizes the mechanical properties of the anisotropy of AA7050-T7451 aluminum alloy in three-dimensional space. The final expression for the stress–strain constitutive equation in the load coordinate system as follows:
σ γ = k σ R D s i n 2 γ + σ N D c o s 2 γ + 2 τ c o s γ s i n γ
The stress–strain curves at arbitrary angles γ, α = 0, and β = 0 in three-dimensional space were predicted by combining Equations (18)–(20), and they are displayed in Figure 11. The accuracy of the constitutive equation containing the spatial angle was verified by predicting the strain curves at different temperatures and strain rates; the error was within 5%.

4.3.4. Simulation Modeling of Dynamic Shock Compression Process

In order to further verify the accuracy of the constructed model, the deform module finite element geometric modeling is carried out according to the size and parameters of the SHPB test device. As shown in Figure 12, in order to reduce the amount of calculation of simulated deformation, the finite element model only retains the incident rod, the specimen and the transmission rod.
The finite element model does not consider the deformation of the device, and the incident rod and transmission rod are defined as rigid bodies. According to the parameters of the JC constitutive model and the material characteristic parameters of 7050-T7451, the mechanical response behavior of the specimen under load is set. Due to the complex changes of field variables in the material during dynamic impact compression, the tetrahedral constant strain element of DEFORM software is used to realize the adaptive meshing of regional field variables. Set the transmission rod to be fixed, and set the end surface of the test piece in contact with the transmission rod to be fixed. The initial incident velocity of the incident rod is set according to the impact rod velocity corresponding to the actual strain rate to strike the specimen. At the same time, the ambient temperature and heat conduction parameters of the specimen are set according to the material and temperature test parameters. Considering the calculation accuracy and efficiency, the friction factor is set to 0.4, the incremental steps are set to 10, and the analysis steps are set to 100.
Through the software post-processing module, the simulated data are processed and analyzed to obtain the changes of stress field, strain field and temperature field under different loading conditions. Figure 13 shows the equivalent stress nephogram of samples with different angles under the temperature of 200 °C and strain rate of 4000 s−1 when the strain is 0.2 during compression deformation. The results show that the maximum stress appears at the edge of the end face of the specimen, which is consistent with the edge spalling of the main failure form of the test specimen.
Figure 14 shows a comparison curve of test values, simulation values and calculated values of materials under different forming angles. By comparing the curves, it can be concluded that the maximum errors of the calculation and simulation results in the plastic deformation stage are 13.3%, 9.2%, and the maximum errors in the elastic stage are 16.6%. The constitutive relation mainly describes the plastic deformation behavior of materials under different loading conditions; thus, the prediction error in the elastic stage is large. However, the error is still within 20%. Therefore, considering the existence of uncontrollable errors in experimental measurement and simulation, the spatial modified constitutive model has high prediction accuracy. It can accurately describe the dynamic compression mechanical behavior of AA7050-T7451 under different forming angles in space.

5. Conclusions

In this work, theoretical expressions for E, G, and the stress σ, which can characterize the mechanical properties of the spatial and orthogonal coordinates at arbitrary angles α, β, and γ, were obtained and verified. The theoretical variations of E, G and σ with spatial angles γ were consistent with the experimental results. A coordinate transformation matrix [A] was obtained using projection relation and the theoretical expressions of σ between orthogonal and load coordinate systems. It was verified in the Results section that the stresses of TD, ND, RD under high strain rate and temperature can be accurately forecasted by the modified orthogonal constitutive model.
Combining the [A] with a modified orthogonal constitutive model, a spatial JC constitutive equation that includes the spatial angle γ was proposed for describing the dynamic mechanical properties of materials at different spatial angles (γ). The ratio between the theoretical and experimental stress under different loading conditions was within the 0.85–0.90 range because the influence of the preforming processes on the anisotropy was ignored. Therefore, the spatial JC constitutive equation was modified by introducing a correction factor k (0.86), and the results demonstrated and verified the accuracy of the theory and the method for characterizing the three-dimensional space anisotropy.

Author Contributions

Conceptualization, Z.W. and Y.P.; Data curation, N.X. and Y.P.; Formal analysis, N.X., Y.P. and Y.Z.; Methodology, N.X.; Project administration, X.F.; Software, Y.Z.; Supervision, X.F.; Writing—original draft, Z.W.; Writing—review & editing, X.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (52175408), Key projects of Natural Science Foundation of Shandong (ZR2020KE022), and the Scientific Research Leader Studio Project of Jinan (CN) (2019GXRC054).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Komarasamy, M.; Alagarsamy, K.; Ely, L.; Mishra, R.S. Characterization of 3″ through-thickness friction stir welded 7050-T7451 Al alloy. Mater. Sci. Eng. A 2018, 716, 55–62. [Google Scholar] [CrossRef]
  2. Krishnan, P.K.; Christy, J.V.; Arunachalam, R.; Mourad, A.-H.I.; Muraliraja, R.; Al-Maharbi, M.; Murali, V.; Chandra, M.M. Production of aluminum alloy-based metal matrix composites using scrap aluminum alloy and waste materials: Influence on microstructure and mechanical properties. J. Alloys Compd. 2019, 784, 1047–1061. [Google Scholar] [CrossRef]
  3. Mishra, S.; Kulkarni, K.; Gurao, N.P. Effect of crystallographic texture on precipitation induced anisotropy in an aluminium magnesium silicon alloy. Mater. Des. 2015, 87, 507–519. [Google Scholar] [CrossRef]
  4. Bhattacharyya, J.J.; Bittmann, B.; Agnew, S.R. The effect of precipitate-induced backstresses on plastic anisotropy: Demonstrated by modeling the behavior of aluminum alloy, 7085. Int. J. Plast. 2019, 117, 3–20. [Google Scholar] [CrossRef]
  5. Linul, E.; Movahedi, N.; Marsavina, L. The temperature and anisotropy effect on compressive behavior of cylindrical closed-cell aluminum-alloy foams. J. Alloys Compd. 2018, 740, 1172–1179. [Google Scholar] [CrossRef]
  6. Safyari, M.; Moshtaghi, M.; Kuramoto, S.; Hojo, T. Influence of microstructure-driven hydrogen distribution on environmental hydrogen embrittlement of an Al–Cu–Mg alloy. Int. J. Hydrog. Energy 2021, 46, 37502–37508. [Google Scholar] [CrossRef]
  7. Sheikh-Ansari, M.H.; Aghaie-Khafri, M. Constitutive modeling of semisolid deformation for the assessment of dilatant shear bands. Appl. Math. Model. 2019, 70, 128–138. [Google Scholar] [CrossRef]
  8. Ryś, M.; Egner, H. Energy equivalence based constitutive model of austenitic stainless steel at cryogenic temperatures. Int. J. Solids Struct. 2018, 164, 52–65. [Google Scholar] [CrossRef]
  9. Abedini, A.; Butcher, C.; Nemcko, M.J.; Kurukuri, S.; Worswick, M.J. Constitutive characterization of a rare-earth magnesium alloy sheet (ZEK100-O) in shear loading: Studies of anisotropy and rate sensitivity. Int. J. Mech. Sci. 2017, 128–129, 54–69. [Google Scholar] [CrossRef]
  10. Lhadi, S.; Berbenni, S.; Gey, N.; Richeton, T.; Germain, L. Micromechanical modeling of the effect of elastic and plastic anisotropies on the mechanical behavior of β-Ti alloys. Int. J. Plast. 2018, 109, 88–107. [Google Scholar] [CrossRef]
  11. Borik, M.A.; Borichevskij, V.P.; Bublik, V.T.; Kulebyakin, A.V.; Lomonova, E.E.; Milovich, F.O.; Myzina, V.A.; Ryabochkina, P.A.; Sidorova, N.V.; Tabachkova, N.Y. Anisotropy of the mechanical properties and features of the tetragonal to monoclinic transition in partially stabilized zirconia crystals. J. Alloys Compd. 2019, 792, 1255–1260. [Google Scholar] [CrossRef]
  12. Marvi-Mashhadi, M.; Lopes, C.S.; Llorca, J. Effect of anisotropy on the mechanical properties of polyurethane foams: An experimental and numerical study. Mech. Mater. 2018, 124, 143–154. [Google Scholar] [CrossRef]
  13. Plunkett, B.; Cazacu, O.; Barlat, F. Orthotropic yield criteria for description of the anisotropy in tension and compression of sheet metals. Int. J. Plast. 2008, 24, 847–866. [Google Scholar] [CrossRef]
  14. Cazacu, O. New yield criteria for isotropic and textured metallic materials. Int. J. Solids Struct. 2018, 139–140, 200–210. [Google Scholar] [CrossRef]
  15. Cardoso, R.P.R.; Adetoro, O.B. A generalisation of the Hill’s quadratic yield function for planar plastic anisotropy to consider loading direction. Int. J. Mech. Sci. 2017, 128–129, 253–268. [Google Scholar] [CrossRef]
  16. Orlowski, K.A.; Ochrymiuk, T.; Sandak, J.; Sandak, A. Estimation of fracture toughness and shear yield stress of orthotropic materials in cutting with rotating tools. Eng. Fract. Mech. 2017, 178, 433–444. [Google Scholar] [CrossRef]
  17. Meng, Y.; Wang, X.; Zongcheng, H.; Fu, X. The Influence of Forming Directions and Strain Rate on Dynamic Shear Properties of Aerial Aluminum Alloy. Appl. Sci. 2018, 8, 520. [Google Scholar] [CrossRef]
  18. Zhang, X.; Liu, Y.; Wei, Y.; Yun, X. Finite element simulation for the punch and stress-strain relationship of anisotropic Aluminum alloy. J. Dalian Jiaotong Univ. 2008, 29, 44–47. [Google Scholar]
  19. Medjahed, A.; Moula, H.; Zegaoui, A.; Derradji, M.; Henniche, A.; Wu, R.; Hou, L.; Zhang, J.; Zhang, M. Influence of the rolling direction on the microstructure, mechanical, anisotropy and gamma rays shielding properties of an Al-Cu-Li-Mg-X alloy. Mater. Sci. Eng. A 2018, 732, 129–137. [Google Scholar] [CrossRef]
  20. Liu, Y.; Xu, D.; Xu, J. Updated 5-parameters barlat-lian yield criteria. J. Huazhong Univ. Sci. Technol. 2008, 36, 129–132. [Google Scholar]
  21. Pilthammar, J.; Sigvant, M.; Kao-Walter, S. Introduction of elastic die deformations in sheet metal forming simulations. Int. J. Solids Struct. 2017, 151, 76–90. [Google Scholar] [CrossRef]
  22. Feng, B.; Bronkhorst, C.A.; Addessio, F.L.; Morrow, B.M.; Li, W.H.; Lookman, T.; Cerreta, E.K. Coupled nonlinear elasticity, plastic slip, twinning, and phase transformation in single crystal titanium for plate impact loading. J. Mech. Phys. Solids 2019, 127, 358–385. [Google Scholar] [CrossRef]
  23. Cyr, E.D.; Brahme, A.; Mohammadi, M.; Mishra, R.K.; Inal, K. A new crystal plasticity framework to simulate the large strain behaviour of aluminum alloys at warm temperatures. Mater. Sci. Eng. A 2018, 727, 11–28. [Google Scholar] [CrossRef]
  24. Tayon, W.A.; Nygren, K.E.; Crooks, R.E.; Pagan, D.C. In-situ study of planar slip in a commercial aluminum-lithium alloy using high energy X-ray diffraction microscopy. Acta Mater. 2019, 173, 231–241. [Google Scholar] [CrossRef]
  25. Muhammad, W.; Sabat, R.K.; Brahme, A.P.; Kang, J.; Mishra, R.K.; Shalchi-Amirkhiz, B.; Hirsch, J.; Inal, K. Deformation banding in a precipitation hardened aluminum alloy during simple shear deformation. Scr. Mater. 2019, 162, 300–305. [Google Scholar] [CrossRef]
  26. Kireeva, I.V.; Chumlyakov, Y.I.; Pobedennaya, Z.V.; Vyrodova, A.V.; Kuksgauzen, I.V.; Kuksgauzen, D.A. Orientation and temperature dependence of a planar slip and twinning in single crystals of Al0.3CoCrFeNi high-entropy alloy. Mater. Sci. Eng. A 2018, 737, 47–60. [Google Scholar] [CrossRef]
  27. Turner, T.J.; Miller, M.P.; Barton, N.R. The influence of crystallographic texture and slip system strength on deformation induced shape changes in AA 7050 thick plate. Mech. Mater. 2002, 34, 605–625. [Google Scholar] [CrossRef]
  28. Wei, Y.L.; Godfrey, A.; Liu, W.; Liu, Q.; Huang, X.; Hansen, N.; Winther, G. Dislocations, boundaries and slip systems in cube grains of rolled aluminium. Scr. Mater. 2011, 65, 355–358. [Google Scholar] [CrossRef]
  29. Cantara, A.M.; Zecevic, M.; Eghtesad, A.; Poulin, C.M.; Knezevic, M. Predicting elastic anisotropy of dual-phase steels based on crystal mechanics and microstructure. Int. J. Mech. Sci. 2019, 151, 639–649. [Google Scholar] [CrossRef]
  30. Gao, K.; Harrison, J.P. Generation of random stress tensors. Int. J. Rock Mech. Min. Sci. 2016, 94, 18–26. [Google Scholar] [CrossRef]
  31. Biel, A.; Stigh, U. Comparison of J-integral methods to experimentally determine cohesive laws in shear for adhesives. Int. J. Adhes. Adhes. 2019, 94, 64–75. [Google Scholar] [CrossRef]
  32. Ribeiro, A.B.; Mascia, N.T. Numerical and experimental study of shear stress behavior of NBR and ASTM standard test specimens for FRP-wood bonds. Compos. Struct. 2019, 224, 111066. [Google Scholar] [CrossRef]
  33. Bruning, T.; Karakus, M.; Nguyen, G.D.; Goodchild, D. An experimental and theoretical stress-strain-damage correlation procedure for constitutive modelling of granite. Int. J. Rock Mech. Min. Sci. 2019, 116, 1–12. [Google Scholar] [CrossRef]
  34. Yan, B.; Regueiro, R.A. Definition and symmetry of averaged stress tensor in granular media and its 3D DEM inspection under static and dynamic conditions. Int. J. Solids Struct. 2019, 161, 243–266. [Google Scholar] [CrossRef]
  35. Ndubuaku, O.; Liu, X.; Martens, M.; Cheng, J.J.R.; Adeeb, S. The effect of material stress-strain characteristics on the ultimate stress and critical buckling strain of flat plates subjected to uniform axial compression. Constr. Build. Mater. 2018, 182, 346–359. [Google Scholar] [CrossRef]
  36. Chandola, N.; Cazacu, O.; Revil-Baudard, B. Prediction of plastic anisotropy of textured polycrystalline sheets using a new single-crystal model. Comptes. Rendus. Méc. 2018, 346, 756–769. [Google Scholar] [CrossRef]
  37. Saleh, M.; Kariem, M.M.; Luzin, V.; Toppler, K.; Li, H.; Ruan, D. High strain rate deformation of ARMOX 500T and effects on texture development using neutron diffraction techniques and SHPB testing. Mater. Sci. Eng. A 2018, 709, 30–39. [Google Scholar] [CrossRef]
  38. Pervier, M.L.A.; Lerma, B.G.; Moncholi, E.P.; Hammond, D.W. A new test apparatus to measure the adhesive shear strength of impact ice on titanium 6Al-4V alloy. Eng. Fract. Mech. 2019, 214, 212–222. [Google Scholar] [CrossRef]
  39. Daryadel, S.S.; Mantena, P.R.; Kim, K.; Stoddard, D.; Rajendran, A.M. Dynamic response of glass under low-velocity impact and high strain-rate SHPB compression loading. J. Non-Cryst. Solids 2015, 432, 432–439. [Google Scholar] [CrossRef]
  40. Eivani, A.R.; Vafaeenezhad, H.; Nikan, O.; Zhou, J. Modeling high temperature deformation characteristics of AA7020 aluminum alloy using substructure-based constitutive equations and mesh-free approximation method. Mech. Mater. 2019, 129, 104–112. [Google Scholar] [CrossRef]
  41. Ye, Z.B.; Huang, R.Y.; Li, Y.C.; Lv, L.; Zhao, K.; Zhang, Y.L.; Ma, J.; Lin, J.J. Steel fiber-reinforced concrete under impact loading dynamic constitutive equation. Constr. Build. Mater. 2018, 190, 1049–1055. [Google Scholar] [CrossRef]
  42. Wang, B.; Liu, Z. Shear localization sensitivity analysis for Johnson–Cook constitutive parameters on serrated chips in high speed machining of Ti6Al4V. Simul. Model. Pract. Theory 2015, 55, 63–76. [Google Scholar] [CrossRef]
  43. Dhananchezian, M.; Kumar, M.P.; Sornakumar, T. Cryogenic Turning of AISI 304 Stainless Steel with Modified Tungsten Carbide Tool Inserts. Mater. Manuf. Process. 2011, 26, 781–785. [Google Scholar] [CrossRef]
  44. Luong, H.; Hill, M.R. The effects of laser peening and shot peening on high cycle fatigue in 7050-T7451 aluminum alloy. Mater. Sci. Eng. A 2010, 527, 699–707. [Google Scholar] [CrossRef]
  45. Mello, A.W.; Nicolas, A.; Lebensohn, R.A.; Sangid, M.D. Effect of microstructure on strain localization in a 7050 aluminum alloy: Comparison of experiments and modeling for various textures. Mater. Sci. Eng. A 2016, 661, 187–197. [Google Scholar] [CrossRef]
  46. Harzallah, M.; Pottier, T.; Senatore, J.; Mousseigne, M.; Germain, G.; Landon, Y. Numerical and experimental investigations of Ti-6Al-4V chip generation and thermo-mechanical couplings in orthogonal cutting. Int. J. Mech. Sci. 2017, 134, 189–202. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Positional relationship and directional cosine between the axes of the orthogonal and load coordinate systems.
Figure 1. Positional relationship and directional cosine between the axes of the orthogonal and load coordinate systems.
Materials 15 05998 g001
Figure 2. Stress component of an arbitrary elastomer section.
Figure 2. Stress component of an arbitrary elastomer section.
Materials 15 05998 g002
Figure 3. (a) Sampling in the typical forming direction. (b) Sampling at different forming angles (γ = 0°, 30°, 45°, 60°, 90°).
Figure 3. (a) Sampling in the typical forming direction. (b) Sampling at different forming angles (γ = 0°, 30°, 45°, 60°, 90°).
Materials 15 05998 g003
Figure 4. (a) SHPB device. (b) Dynamic “Bar-Tube” shearing device.
Figure 4. (a) SHPB device. (b) Dynamic “Bar-Tube” shearing device.
Materials 15 05998 g004
Figure 5. Relationship between the orthogonal and load coordinate systems of AA7050-T7451.
Figure 5. Relationship between the orthogonal and load coordinate systems of AA7050-T7451.
Materials 15 05998 g005
Figure 6. Stress–strain curves of AA7050-T7451 in different directions. (a) 0.1 × 104 s1; (b) 0.2 × 104 s1; (c) 0.4 × 104 s1; (d) 0.6 × 104 s1.
Figure 6. Stress–strain curves of AA7050-T7451 in different directions. (a) 0.1 × 104 s1; (b) 0.2 × 104 s1; (c) 0.4 × 104 s1; (d) 0.6 × 104 s1.
Materials 15 05998 g006
Figure 7. High-temperature mechanical characteristic curve of AA7050-T7451 (strain rate: 0.4 × 104 s1). (a) Stress–strain curves at different temperatures (RD). (b) Variation in the stability stress with γ at different temperatures.
Figure 7. High-temperature mechanical characteristic curve of AA7050-T7451 (strain rate: 0.4 × 104 s1). (a) Stress–strain curves at different temperatures (RD). (b) Variation in the stability stress with γ at different temperatures.
Materials 15 05998 g007
Figure 8. Fitting curves of the thermal softening coefficient for different directions with a strain rate of 0.4 × 104 s1.
Figure 8. Fitting curves of the thermal softening coefficient for different directions with a strain rate of 0.4 × 104 s1.
Materials 15 05998 g008
Figure 9. Elastic modulus and shear modulus with the change in γ: (a) elastic modulus, (b) shear modulus.
Figure 9. Elastic modulus and shear modulus with the change in γ: (a) elastic modulus, (b) shear modulus.
Materials 15 05998 g009
Figure 10. Comparison and verification of the theoretical and experimental values of the stress in the load coordinate system. (a) Test verification at different strain rates. (b) Test verification at different temperatures.
Figure 10. Comparison and verification of the theoretical and experimental values of the stress in the load coordinate system. (a) Test verification at different strain rates. (b) Test verification at different temperatures.
Materials 15 05998 g010
Figure 11. Stress–strain prediction curves under different loading conditions, (a) Predicted stress–strain curve at a strain rate of 0.2 × 104 s1. (b) Predicted stress–strain curves of the 45° specimens at different temperatures.
Figure 11. Stress–strain prediction curves under different loading conditions, (a) Predicted stress–strain curve at a strain rate of 0.2 × 104 s1. (b) Predicted stress–strain curves of the 45° specimens at different temperatures.
Materials 15 05998 g011
Figure 12. Schematic diagram of dynamic impact compression finite element model.
Figure 12. Schematic diagram of dynamic impact compression finite element model.
Materials 15 05998 g012
Figure 13. Comparison of simulated equivalent stress nephogram and damage forms of workpieces: (a) ε = 0.2, RD-30°-TD; (b) ε = 0.2, RD-45°-TD; (c) ε = 0.2, RD-60°-TD; (d) comparison diagram of sample damage.
Figure 13. Comparison of simulated equivalent stress nephogram and damage forms of workpieces: (a) ε = 0.2, RD-30°-TD; (b) ε = 0.2, RD-45°-TD; (c) ε = 0.2, RD-60°-TD; (d) comparison diagram of sample damage.
Materials 15 05998 g013
Figure 14. Comparison of test value, simulation value and model value of samples at different angles: (a) RD-30°-TD; (b) RD-45°-TD; (c) RD-60°-TD.
Figure 14. Comparison of test value, simulation value and model value of samples at different angles: (a) RD-30°-TD; (b) RD-45°-TD; (c) RD-60°-TD.
Materials 15 05998 g014
Table 1. Correction values of the JC constitutive parameters A, B, C, and n.
Table 1. Correction values of the JC constitutive parameters A, B, C, and n.
Forming DirectionsA (MPa)B (MPa)Cn
ND 3655160.0420.26
TD 3574900.040.3
RD 3424520.0290.26
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Fu, X.; Xu, N.; Pan, Y.; Zhang, Y. Spatial Constitutive Modeling of AA7050-T7451 with Anisotropic Stress Transformation. Materials 2022, 15, 5998. https://doi.org/10.3390/ma15175998

AMA Style

Wang Z, Fu X, Xu N, Pan Y, Zhang Y. Spatial Constitutive Modeling of AA7050-T7451 with Anisotropic Stress Transformation. Materials. 2022; 15(17):5998. https://doi.org/10.3390/ma15175998

Chicago/Turabian Style

Wang, Zhenda, Xiuli Fu, Nianwei Xu, Yongzhi Pan, and Yijia Zhang. 2022. "Spatial Constitutive Modeling of AA7050-T7451 with Anisotropic Stress Transformation" Materials 15, no. 17: 5998. https://doi.org/10.3390/ma15175998

APA Style

Wang, Z., Fu, X., Xu, N., Pan, Y., & Zhang, Y. (2022). Spatial Constitutive Modeling of AA7050-T7451 with Anisotropic Stress Transformation. Materials, 15(17), 5998. https://doi.org/10.3390/ma15175998

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