Print version ISSN 1516-1439
Mat. Res. vol.16 no.1 São Carlos Jan./Feb. 2013 Epub Nov 13, 2012
Guo-zheng Quan*; Gui-sheng Li; Yang Wang; Wen-quan Lv; Chun-tang Yu; Jie Zhou
School of Material Science and Engineering, Chongqing University, Chongqing 400044, China
In order to perform the numerical simulations of forging response and establish the processing parameters for as-extruded 7075 aluminum alloy, the compressive deformation behavior of as-extruded 7075 aluminum alloy were investigated at the temperatures of 573 K, 623 K, 673 K and 723 K and the strain rates of 0.01 s1, 0.1 s1, 1 s1 and 10 s1 on a Gleeble1500 thermo-mechanical simulator. Based on the analysis of the effect of strain, temperature and strain rate on flow stress, dynamic recrystallization (DRX) type softening characteristics of the stress-strain curve with single peak were identified. The traditional Arrhenius type model is in favor of the prediction for the flow stress at a fixed strain, and can not satisfy the need of the numerical simulations of various hot forming processes due to the lack of the effect of strain on flow stress. Lin et al. improved Arrhenius type model with a series of variable coefficients as functions of true strain (including activation energy of deformation Q, material constants n and a, and structure factor A) to predict the flow stress during the hot compression. The application has been demonstrated in this work for as-extruded 7075 aluminum alloy. The comparisons between the predicted and experimental results show that, for the worst case, the error in the flow stress estimate is 5.63%, and the max mean error is 3.6%. The developed model provides fast, accurate and consistent results, making it superior to the conventional Arrhenius type model. In further it can be used in computer code to model the forging response of 7075 aluminum alloy mechanical part members under the prevailing loading conditions.
Keywords: flow stress, constitutive equation, material constants, aluminum alloy
Aluminum alloys of the 7075 type (Al-Zn-Mg-Cu) have some advantageous properties that make them an excellent choice for a number of applications in aerospace component manufacture1. The understanding of metals and alloys behavior at hot deformation condition has a great importance for designers of hot metal forming processes (hot rolling, forging and extrusion) because of its effective role on metal flow pattern as well as the kinetics of metallurgical transformation, and the constitutive relation is often used to describe the plastic flow properties of the metals and alloys in a form that can be used in computer code to model the forging response of mechanical part members under the prevailing loading conditions2. In common, metals and alloys are very sensitive to deformation conditions, the shape of the workpiece and product, the shapes of the tools, friction, temperature, forming speed, etc4. In such parameters, the constitutive relationship is one of the most important factors for the analysis problem of metal forming process, and understanding their effects has been a difficult task due to their complex nature3-5.
A number of research groups have attempted to develop constitutive equations of materials and suggested their own formulations by putting the experimentally measured data into one single equation. Shida6 and Misaka7 suggested the models giving the stress of carbon steels as a function of the strain, strain rate, temperature and carbon content. Johnson and Cook8 developed a constitutive equation which assumes that the dependence of the stress on the strain, strain rate and temperature can be multiplicatively decomposed into three separate functions that include five constants to be determined by experimental data obtained for a specific material. Lee et al.9 established the flow stress constitution equation of 7075 aluminum alloy using a split Hopkinson bar (SHPB) experiment. Jamal et al.10 investigated the ANN constitutive model for high strain-rate deformation of Al 7075-T6. The traditional Arrhenius type model is in favor of the prediction for the flow stress at a fixed strain. Due to the lack of the effect of strain on flow stress, it can not satisfy the need of the numerical simulations of various hot forming processes. Lin et al.2 introduced a series of variable coefficients (including activation energy of deformation Q, material constants n and a, and structure factor A) as functions with strain into Arrhenius type model, and then the influence of strain on flow stress were created. Such improved constitutive model was remarkable for its higher accuracy and applied to 42CrMo steel successfully. However, there is still no available information indicating such model applied in as-extruded 7075 aluminum alloy1-8.
The object of this study is to construct the constitutive description equation for as-extruded 7075 aluminum alloy which represents the relationships of strain, temperature, strain rate and flow stress. The main work is to uncover and describe the general nature of the influence of strain on a series of variable coefficients (including activation energy of deformation Q, material constants n and a, and structure factor A) in Arrhenius type model for as-extruded 7075 aluminum alloy using a series of true stress-strain curves collected from compression tests at the temperature of 573~723 K and the strain rate of 0.01~10 s1. And then an improved Arrhenius type constitutive equation incorporating a series of polynomial functions for each coefficient was constructed for as-extruded 7075 aluminum alloy. The validity and effectiveness of the proposed equation were also demonstrated by the max mean error of 3.6% for a fixed strain rate.
2. Materials and Experimental Methods
The chemical compositions of as-extruded 7075 aluminum alloy (Al-Zn-Mg-Cu) used in this study was as follows (wt. %): Zn-5.5, Mg-2.2, Cr-2.2, Cu-1.7, Si-0.4, Fe-0.3, Mn-0.1, Al (balance). The homogenized ingot was scalped to diameter 10 mm and height 12 mm, and seventeen cylindrical specimens were prepared. A computer-controlled servo-hydraulic Gleeble1500 machine was used for the isothermal compression experiments. The specimens were resistance heated at a heating rate of 1 K/s and held at a fixed temperature for 180 s by thermo-coupled-feedback-controlled AC current, which decreased the anisotropy in flow deformation behavior effectively. Three thermocouples were welded at the mid-span of billet to provide accurate temperature control and measurement during testing. During the tests, the temperatures of specimens were controlled to within ± 0.1 K through three thermocouples adhered to the specimens at the gauge zero. The compression tests corresponding to a height reduction of 60% were carried out at four different temperatures of 573 K, 623 K, 673 K and 723 K and four different strain rates of 0.01 s1, 0.1 s1, 1 s1, 10 s1. During the compression process, the variations of stress and strain were monitored continuously by a personal computer equipped with an automatic data acquisition system. The true stress and true strain were derived from the measurement of the nominal stress-strain relationship according to the following formula: σT = σN(1 + εN), εT = ln (1 + εN), where σT is the true stress, σN is the nominal stress, εT is the true strain and εN is the nominal strain11.
3. Results and Discussion
3.1. Flow stress behavior of as-extruded 7075 Aluminum alloy
The true stressstrain curves of as-extruded 7075 Aluminum alloy compressed at different deformation conditions are shown in Figure 1a-d. The flow stress as well as the shape of the flow curve is sensitively dependent on temperature and strain rate. The general characteristics of the flow stress curves are similar at all deformation conditions. For every curve, after a rapid increase in the stress to a peak value, the flow stress decreases monotonically towards a steady state regime with a varying softening rate which typically indicates the onset of DRX. It reveals that 7075 aluminum alloy has an apparent DRX soften characteristic. Comparing these curves with one another, it is found that the stress level decreases with temperature increasing and strain rate decreasing because lower strain rate and higher temperature provide longer time for the energy accumulation and higher mobilities at boundaries which result in the nucleation and growth of dynamically recrystallized grains and dislocation annihilation2,11,12.
3.2. Constitutive equations for flow behavior of as-extruded 7075 aluminum alloy
So far, several empirical equations have been proposed to determine the deformation activation energy and hot deformation behavior of alloys. The most frequently used is Arrhenius equation which describes the behavior of the material deformed at different temperatures and strain rates. The effects of the temperatures and strain rate on the deformation behavior can be represented by Zener-Hollomon parameter, Z in an exponent-type equation. The hyperbolic law in the Arrhenius type equation gives better approximations between Z parameter and stress13-16.
in which, is the strain rate (s1), R is the universal gas constant (8.31 J.mol1.K1), T is the absolute temperature (K), Q is the activation energy of DRX (kJ.mol1), σ is the flow stress (MPa) for a given stain, A, α and n are the material constants, α = β/n. For the low stress level (α|σ|<0.8) and high stress (α|σ|>1.2), by substituting power law and exponential law of F(σ)into Equation 2, the relationships between flow stress and strain rate can be expressed as following equations respectively
Where B and B' are the material constants that are dependent of deformation temperatures.
It is commonly accepted that the effect of strain on stress hasn't been considered in Equation 1-4. Here the effects of deformation strain on stress are investigated by the consideration of the influence of strain on a series of variable coefficients (including activation energy of deformation Q, material constants n and a, and structure factor A) in Arrhenius type model. The following is taking the strain of 0.1 as an example.
3.2.1. Calculation of the material constants n, α and β
By taking the natural logarithm of both sides of Equation 3 and 4, respectively, gives
Then, , and . Substituting the values of the flow stress and corresponding strain rate at the strain of 0.1 into the logarithm Equations 5 and 6 gives the relationship between stress and strain rate as shown in Figure 2. Figure 2a, b shows the relationships of |σ| In|| and In|σ| In|| for ε = 0.1 at the temperatures of 573 K, 623 K, 673 K and 723 K. It can be seen that the stresses under different temperatures and strain rates can be approximated by a group of parallel and straight lines. In Figure 2a the average value of all the lines' slopes is equal to the inverse of n-value, thus n = 7.5565 MPa1 for ε = 0.1. Meanwhile in Figure 2b the average value of all the lines' slopes is equal to the inverse of b-value, thus b = 0.0751 MPa1 for ε = 0.1. Then α = β/n = 0.0087.
3.2.2. Calculation of the apparent activation energy of deformation Q
The apparent activation energy of deformation is a term defined as the energy that must be overcome in order for the nucleation and growth of new surface or grain boundary to occur. Activation energy may also be defined as the minimum energy required to start DRX, and it can be thought as the height of the potential energy barrier. The apparent activation energy of deformation is usually denoted by Q, which an important physical parameter is serving as an indicator of deformation difficulty degree in the plasticity deformation theory. For all the stress level (including low and high stress levels), Equation 2 can be represented as the following:
If is constant, there is a linear relationship between In sinh(α|σ|) and 1/T, and Equation 7 can be rewritten as
The stresses at the strain of 0.1 under different temperatures and strain rates can be identified for the present target stresses. The linear relationships between In sinh(α|σ|) and 1/T at different strain rates were fitted out as Figure 3a. The mean value of all the slope rates is accepted as Q/(Rn) value, furthermore Q value for ε = 0.1 are obtained as 96.1749 kJ.mol1.
3.2.3. Construction of constitutive equation at strain of 0.1
Taking the logarithm of both sides of Equation 7 gives:
By substituting the values of the flow stress and strain rate for all the tested temperatures into Equation 9, the linear relationships between In [sinh(α|σ|)] and In|| at different temperatures were fitted out as Figure 3b. The mean value of all the intercepts of In [sinh(α|σ|)] versus In|| plots is accepted as A value, furthermore A value for ε = 0.1 are obtained as 2.628 × 1012 s1. Substituting n, α, Q and A into Equation 7 gives the relationships between, T and σ as following.
Substituting Zener-Hollomon parameter Z = ||exp(Q/RT) into Equation 2, thus, the flow stress can be expressed as Equation 11.
Here, substituting Z = ||exp[(96.1749 ×103) 8.31T] , and A = 2.628 × 1012s1 into Equation 11, the constitutive equation for flow behavior of as-extruded 7075 aluminum alloy can be expressed as Equation 12.
3.2.4. Construction of constitutive equation in the strain range of 0~0.9
The above is taking the strain of 0.1 as an example, while the following is about the influence of strain (0~0.9) with interval of 0.1 on flow stress. A series of variable coefficients (including activation energy of deformation Q, material constants n, and α, and structure factor A) in Arrhenius type model at different deformation strains within the range of 0~0.9 and the interval of 0.1 were computed by the way as above. The relationships between Q, In A, n, α and true strain e for as-extruded 7075 aluminum alloy were polynomial fitted by the compensation of strain, as shown in Equation 13. The polynomial fitted results of Q, In A, n, α of as-extruded 7075 aluminum alloy are provided in Table 1. The variation of material constants are shown in Figure 4.
Substituting a, n, Q and A in Equation 13 into Equation 7 gives the relationships between , T and σ as Equation 14.
Thus, Z = ||exp[j(ε)/) 8.31T]. Furthermore, the constitutive equation for flow behavior of as-extruded 7075 aluminum alloy in a wide strain range of 0~0.1 can be expressed as Equation 15.
where f(ε), g(ε), h(ε) and j(ε) are the polynomial functions of A, α, n and Q at different true strains, and their expressions are as Equation 13 and Table 1.
3.3. Adaptability evaluation of the developed constitutive equation
In order to verify the developed constitutive equations of 7075 aluminum alloy, the comparison between the experimental and predicted results was carried out. Using the Equation 15, by applying the determined material constants of the constitutive equation and fixed the ε, , T, the flow stress can be calculated for strain rates between 0.01 s1~10 s1 and temperatures between 573K~723K. Figure 5 shows comparisons between the experimental and predicted results by the developed Equation 15 considering the compensation of strain at the temperatures of 573 K, 623 K, 673 K and 723 K, and the strain rates of 0.01 s1, 0.1 s1, 1 s1 and 10 s1. It can be seen that the calculation results are well in agreement with the experimental results not only at the strain-hardening stage but also at the strain-softening stage. In order to evaluate the detailed accuracy of the developed constitutive equation of as-extruded 7075 aluminum alloy, the error between the calculated flow stress (σC) and measured flow stress (σM) can be calculated as Equation 16. Table 2 lists the detailed evaluation results for flow stress at the true strains of 0.3, 0.5 and 0.7, and the strain rates of 0.01 s1, 0.1 s1, 1 s1 and 10 s1. It can be easily found that, for the worst case, the error in the flow stress estimate is 5.63%, and the max mean error for a fixed strain rate is 3.6%. The results show that applying improved Arrhenius model with variable parameters to predict the flow stress has higher accuracy degree. The cause lies in that the improved constitutive model with a series of variable coefficients as function of true strain, including activation energy of deformation Q, material constants n and a, and structure factor A. Research show that the proposed deformation constitutive equation gives an accurate and precise estimate of the flow stress for as-extruded 7075 aluminum alloy, and can be used to analyze the problems during metal forming process.
3.4. Potentiality in practice
With the development of industrial technology, it is more and more significant to evaluate and ensure the intrinsic workability of deforming workpiece corresponding to the microstructure evolution and deformation mechanism transformation17. One of the most potential applications is FEM, which is popularly used to simulate hot deformation processes and predict the problems in production. In the deforming workpiece the distribution of the strain, strain rate and temperature can be calculated by FEM, and constitutive relation is the basis of numerical simulation. The constitutive model could be applied effectively in computer codes for simulation by two main techniques. As for the first approach, the constitution equation could be applied in FEM by inserting an in-house code using FORTRAN language into the commercial software such as ABAQUS' Marc, etc. As for the second approach, the intensive data predicted by the developed Equation 15 are inputted into the metallic material warehouse of FEM software directly. By contrast, the first method is with higher efficiency and accuracy. However, whether the simulation result is feasible or not depends on the accuracy degree of constitutive models. This improved constitutive model is remarkable for its higher accuracy and can be applied to 7075 aluminum alloy successfully.
In the present study, the flow behavior of as-extruded 7075 aluminum alloy was investigated using hot compression tests in the temperature range of 573~723 K and the strain rate range of 0.01~10 s1. The main work is to construct the constitutive description equation for as-extruded 7075 aluminum alloy which represents the relationships of strain, temperature, strain rate and flow stress. The application of improved Arrhenius type model with a series of variable coefficients as functions of true strain (including activation energy of deformation Q, material constants a, and n, and structure factor A) to predict the flow stress during the hot compression of as-extruded 7075 aluminum alloy has been demonstrated in this work. The comparisons between the predicted and experimental results show that, for the worst case, the error in the flow stress estimate is 5.63%, and the max mean error is for a fixed strain rate 3.6%. The developed model provides fast, accurate and consistent results, making it advantageous regarding the conventional Arrhenius type model. In further it can be used in computer code to model the forging response of 7075 aluminum alloy mechanical part members under the prevailing loading conditions.
This work was supported by National Key Technologies R & D Program of China (2012ZX04010-081), Science and Technology Committee of Chongqing (cstc2009aa3012-1), Fundamental Research Funds for the Central Universities (CDJZR11130009), Sharing Fund of Chongqing University's Large-scale Equipment (2011063014).
1. Quan GZ, Liu KW and Zhou J. Dynamic softening behaviors of 7075 aluminum alloy. Transactions of nonferrous metals society of China. 2009; 19(s3):s537-s541. http://dx.doi.org/10.1016/S1003-6326(10)60104-5 [ Links ]
2. Lin YC, Chen MS and Zhang J. Modeling of flow stress of 40CrMo steel under hot compression. Materials science and engineering A. 2009; 499(1-2):88-92. http://dx.doi.org/10.1016/j.msea.2007.11.119 [ Links ]
3. Dan WJ, Zhang WG, Li SH and Lin ZQ. Finite element simulation on strain-induced martensitic transformation effects in TRIP steel sheet forming. Computational materials science. 2007; 39(3): 593-599. http://dx.doi.org/10.1016/j.commatsci.2006.08.013 [ Links ]
4. Kim SI, Lee Y and Byon SM. Study on constitutive relation of AISI 4140 steel subject to large strain at elevated temperatures. Materials Processing Technology. 2003; 140(1-3):84-89. http://dx.doi.org/10.1016/S0924-0136(03)00742-8 [ Links ]
5. Berbenni S, Favier V and Berveiller M. Impact of the grain size distribution on the yield stress of heterogeneous materials. Computational materials science. 2007; 24(1):114-142. [ Links ]
6. Shida S. Empirical formula of flow stress of carbon steels resistance to deformation of carbon steels at elevated temperature. 2nd repore. Journal of The JSTP. 1969; 10:610-617. [ Links ]
7. Misaka Y and Yoshimoto T. Formulation of mean resistance of deformation of plain carbon steel at elevated temperature. Journal of the Japan Society for Technology of Plasticity. 1967-1968; 8:414-422. [ Links ]
8. Johnson GR and Cook WH. A constitutive model and data for matals subjected to large strains, high strain rates and high temperature. Proceedings of the Seventh International Symposium on Ballistics. 1983; 541-547. [ Links ]
9. Lee WS, Sue WC, Lin CF and Wu CJ. The strain rate and temperature dependence of the dynamic impact properties of 7075 aluminum alloy. Materials Processing Technology. 2000; 100:116-122. http://dx.doi.org/10.1016/S0924-0136(99)00465-3 [ Links ]
10. Jamal AS and Janet T. ANN constitutive model for high strain-rate deformation of Al 7075-T6. Materials Processing Technology. 2007; 186:339-345. http://dx.doi.org/10.1016/j.jmatprotec.2006.11.228 [ Links ]
11. Quan GZ, Li GS, Chen T, Wang YX, Zhang YW and Zhou J. Dynamic recrystallization kinetics of 42CrMo steel during compression at different temperatures and strain rates. Materials Science and Engineering A. 2011; 528(13-14):4643-4651. http://dx.doi.org/10.1016/j.msea.2011.02.090 [ Links ]
12. Lei LM, Xu H, Wang MM, Wang L, Qin J, Li H et al. Effect of hot compressive deformation on the martensite transformation of Ti-10V-2Fe-3Al titanium alloy. Materials Science and Engineering A. 2011; 530:591-601. http://dx.doi.org/10.1016/j.msea.2011.10.028 [ Links ]
14. Shi H, McLaren AJ, Sellars CM, Shahani R and Bolingbroke R.Constitutive equation for high temperature flow stress of aluminum alloys. Materials Science and Technology. 1997; 13:210-216. http://dx.doi.org/10.1179/026708397790302421 [ Links ]
16. Han MY, Chen ZH, Zhang H and Xia WJ. Flow stress behavior of porous FVS0812 aluminum alloy during hot-compression. Mechanics Research Communications. 2006; 33:508-514. http://dx.doi.org/10.1016/j.mechrescom.2005.05.008 [ Links ]
17. Quan GZ, Kang BS, Ku TW and Song WJ. Identification for the optimal working parameters of Al-Zn-Mg-Cu alloy with the processing maps based on DMM. Advanced manufacturing technology. 2011; 56(9-12):1069-1078. [ Links ]
Received: May 6, 2012
Revised: June 11, 2012