Modeling Material Flow Behavior during Hot Deformation Based on Metamodeling Methods

Abstract

Modeling material flow behavior is an essential step to design and optimize the forming process. In this context, four popular metamodel types Kriging, radial basis function, multivariate polynomial, and artificial neural network are investigated as potential methods for modeling the flow behavior of 6013 aluminum alloy. Based on the experimental data from hot compression tests, the modeling performance of these four methods was tested and subsequently compared from different aspects. It is found that all the methods are capable of constructing models for describing the hot deformation behavior. The merits of Kriging method over other three methods are highlighted when the sample size for modeling is decreased. Furthermore, the applicability of Kriging method is validated while decreasing the sample uniformity with respect to temperature or strain rate. It is proved that Kriging method is competent in modeling the material flow behavior and is the most effective one among the four popular types of metamodeling method.

Gang Xiao 1, 2 and Qinwen Yang 2 and Luoxing Li 1, 2 and Zhengbing Xu 3

Academic Editor:Mohsen Asle Zaeem

1, State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body, Hunan University, Changsha 410082, China
2, College of Mechanical and Vehicle Engineering, Hunan University, Changsha 410082, China
3, Key Laboratory of Nonferrous Materials and New Processing Technology of Ministry of Education of China, Guangxi University, Nanning 530004, China

Received 27 July 2015; Accepted 30 August 2015; 20 September 2015

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Hot forming technology of materials is widely applied in practical manufacturing. In order to achieve the required microstructure and mechanical property of products, the design and optimization of thermomechanical process parameters have to be properly implemented [1, 2]. For this reason, the designers strive to study the effect of the parameters with respect to the kinetics of the metallurgical deformation. The flow behavior of materials is compactly affected by the factors of temperature, strain rate, and strain [3]. For exploring and describing the hot deformation behavior, many research groups have made use of regression method to establish the phenomenological and physically based constitutive models [4-6], which would be committed to provide a complete mathematical description of the flow stress. Lin et al. [5] developed a new phenomenological model for describing high-temperature flow behavior of inconel 718 superalloy. Saadatkia et al. [6] developed a physically based model of low and medium carbon steels to investigate the deformation behaviors. However, the deformation behavior of metal at elevated temperature is always associated with various metallurgical phenomena and thereby complicated in nature [7], such as work hardening, dynamic recovery, dynamic recrystallization, and flow instabilities. The effects of these factors on flow stress are so complex that the relationship among them is highly nonlinear. It can be found that the physical procedure is difficult to be systematically interpreted by the conventional constitutive model [8], which would reduce the prediction accuracy and limit the application range. Moreover, the development of modeling flow behavior in the conventional way is usually time-consuming and relies on amounts of experimental data.

Considering the disadvantages of the abovementioned method, artificial neural network (ANN) method is gradually utilized as an alternative approach for modeling material flow behavior to improve processing scheme [8-10]. In essence, ANN is one type of metamodeling method which can get rid of the constraints from describing physical mechanisms. Metamodeling method is widely accepted as a valuable and efficient technique to model various complex nonlinear relationships. In previous researches [11], some other types of metamodeling method were maturely developed, such as Kriging, radial basis function (RBF), and multivariate polynomial method. Each metamodel type has its associated fitting method and the corresponding characteristics [12]. Due to the complicated nature of the relationships between metamodeling methods and engineering problems, there is no conclusion on which method is definitely superior to the others [11]. In order to investigate the applicability of metamodeling methods and search for the most appropriate one in modeling flow behavior, the four popular methods Kriging, RBF, multivariate polynomial, and ANN were utilized and compared for modeling the flow stress of 6013 aluminum alloy during hot compression.

In this work, the hot plane strain compression tests were carried out on Gleeble-3500 thermomechanical simulator. Considering deformation heating and heat transfer, the measured flow stress data at relatively high strain rates were corrected. Based on the measured and corrected data, the feasibility of these methods is firstly tested in modeling the flow behavior. The most effective modeling method with satisfactory accuracy, stability, and efficiency was then selected on the basis of comparative study and detailed analysis. Finally, the performance of Kriging method that is considered as the best one was further assessed for this modeling process.

2. Metamodeling Methods

2.1. Kriging Method

After being proposed by Krige [13] and improved by others such as Matheron [14], Kriging method was gradually developed to be a popular metamodel type and systematically introduced into the area of computer experiment [15]. In Kriging method, the random output is assumed to be obtained from a linear combination of regression functions plus a random process factor as follows: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the number of regression functions, [figure omitted; refer to PDF] is a regression function, [figure omitted; refer to PDF] is the coefficient for [figure omitted; refer to PDF] , [figure omitted; refer to PDF] is the design point, and [figure omitted; refer to PDF] is the random process function. The correlation function is defined by the following equation: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the correlation function, [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are two design points, and [figure omitted; refer to PDF] is a structural parameter to be optimized. The correlation function could be defined in several different ways such as Gaussian, exponential, linear, spherical, and cubical [16].

2.2. Radial Basis Function Method

Radial basis function (RBF) method was initially developed as an exact interpolation technique for data in multidimensional space [17]. In RBF method, a series of center points ( [figure omitted; refer to PDF] ) are chosen first from the design points ( [figure omitted; refer to PDF] ) based on some criteria. Then the basis functions are constructed by using these center points as [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the design point, [figure omitted; refer to PDF] is the center point for [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] is the Euclidian distance between the two points. The function [figure omitted; refer to PDF] could have different forms such as linear, Gaussian, cubic, thin-plate, spline, and multiquadratic [17]. The relationship between inputs and outputs is constructed as a linear combination of radial basis functions: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the estimated output value, [figure omitted; refer to PDF] is a basis function, [figure omitted; refer to PDF] is a constant, and [figure omitted; refer to PDF] is the coefficient for [figure omitted; refer to PDF] .

2.3. Multivariate Polynomial Method

Multivariate polynomial method is one of the most fundamental metamodeling methods used in computer experiment. It is mostly known from response surface method [18], which employs quadratic polynomial in the field of engineering design optimization. In multivariate polynomial method, basis functions are built directly by using input variable components ( [figure omitted; refer to PDF] ) and their interactions such as [figure omitted; refer to PDF]

A multivariate polynomial is constructed as the weighted sum of these basis functions: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the estimated output value, [figure omitted; refer to PDF] is a basis function, [figure omitted; refer to PDF] is a constant, and [figure omitted; refer to PDF] is the coefficient for [figure omitted; refer to PDF] .

2.4. Artificial Neural Network Method

Artificial neural network (ANN) method is based on a particular set of nonlinear functions to build the relationship between inputs and outputs. Its architecture consists of input layer, output layer, and hidden layer, which are connected by the processing units called neurons [19], as shown in Figure 1. Each neuron in the input layer and output layer represents one independent variable, while the neurons in the hidden layers are only for computation purpose. During training procedure, the known input-output pairs are used to update the weights and biases [20]. The objective is to minimize the errors between the ANN outputs and the targets for the corresponding inputs.

Figure 1: The architecture of the ANN model.

[figure omitted; refer to PDF]

3. Model Construction and Discussion

3.1. Experiments and Modeling Preparation

The material used in the present experiment is 6013 aluminum alloy, with a nominal composition of 0.95% Mg, 0.75% Si, 0.9% Cu, and 0.35% Mn (mass fraction). The cubic specimens with a size of 20 mm × 15 mm × 10 mm were prepared for testing. The hot plane strain compression tests were carried out 25 times using Gleeble-3500 thermomechanical simulator at different temperatures and strain rates (listed in Table 1). All specimens were heated to the preset deformation temperature at a rate of 10 K/s and held for 120 s to eliminate thermal gradient. The specimens were compressed to a true strain of 0.8 and cooled down to room temperature in water. The experimental procedure is explained in Figure 2.

Table 1: Experimental conditions.

Level	1	2	3	4	5
Temperature (K)	613	653	693	733	773
Strain rate (s-1 )	0.001	0.01	0.1	1	10

Figure 2: Experimental procedure for the compression tests.

[figure omitted; refer to PDF]

It is widely accepted that the flow stress data measured directly from tests at relatively high strain rates are always affected by instantaneous temperature rise during deformation [21, 22]. According to the previous work of the authors [23], the temperature rise during hot plane strain compression depends on deformation heating integrated with heat transfer which does occur among the deformed region, the undeformed region, and the press indenter. This type of temperature rise considering deformation heating and heat transfer has been quantitatively described by a proposed model in this literature. Based on the published work, the measured flow stress data at strain rates of 10 and 1 s^-1 are corrected and then used for comparison and analysis of different modeling methods in this paper.

Based on the 25 experimental (including measured and corrected) flow stress curves, the true stress data corresponding to the true strain between 0.05 and 0.8 with the interval of 0.05 were adopted to construct models and test their performance. It can be found that a total of 400 sample points are located on the 25 curves which are considered as 25 cases in the modeling and testing processes. These data will be divided into training and testing datasets according to the requirements in different sections of this paper.

The schematic of modeling the flow behavior of 6013 aluminum alloy is shown in Figure 3. The inputs of the model are temperature ( [figure omitted; refer to PDF] ), strain rate ( [figure omitted; refer to PDF] ), and strain ( [figure omitted; refer to PDF] ) whereas flow stress ( [figure omitted; refer to PDF] ) is the output. In addition, the parameters of temperature and strain rate are processed by natural logarithm.

Figure 3: Schematic of modeling flow behavior.

[figure omitted; refer to PDF]

3.2. Model Construction and Applicability Assessment

The feasibility of metamodeling methods for modeling the flow behavior is investigated in this section. Three cases with totally different temperatures and strain rates are successively selected for model testing, as shown in Table 2. In each testing process, 24 cases (excluding the testing one) are utilized for model training based on the schematic given in Figure 3. The average ( [figure omitted; refer to PDF] ) and the maximum ( [figure omitted; refer to PDF] ) values of absolute relative errors ( [figure omitted; refer to PDF] ) from each testing are employed as quantitative criteria to evaluate the performance of models and the applicability of modeling methods: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the experimental value, [figure omitted; refer to PDF] is the predicted value, and [figure omitted; refer to PDF] is the number of testing points.

Table 2: Three testing cases to analyze the prediction accuracy of models.

Number	[figure omitted; refer to PDF] (K)	[figure omitted; refer to PDF] (s-1 )	Kriging	RBF	Polynomial	ANN
[figure omitted; refer to PDF] (%)	[figure omitted; refer to PDF] (%)	[figure omitted; refer to PDF] (%)	[figure omitted; refer to PDF] (%)	[figure omitted; refer to PDF] (%)	[figure omitted; refer to PDF] (%)	[figure omitted; refer to PDF] (%)	[figure omitted; refer to PDF] (%)
1	653	1	2.87	8.37	8.24	11.19	3.47	5.75	1.57	2.79
2	693	0.1	3.79	6.19	5.91	8.56	2.42	4.26	3.44	5.68
3	733	0.01	3.92	8.53	1.20	6.26	3.11	5.97	9.49	13.85

As shown in Table 2, the detailed information on [figure omitted; refer to PDF] and [figure omitted; refer to PDF] from the developed models is provided. Except few specific conditions, the values of [figure omitted; refer to PDF] are all less than 10%. Meanwhile, Figure 4 displays the comparisons between the predicted values and experimental curves. All the developed models can basically track the dynamic change of stress with the strain increasing at different temperatures and strain rates. Therefore, it can be considered that the four metamodeling methods Kriging, RBF, multivariate polynomial, and ANN are all capable of modeling the elevated temperature flow behavior of 6013 aluminum alloy. However, it is really difficult to find out which one is the most appropriate method based on the existing testing results. In order to make a reasonable selection from these modeling methods, further comparison and analysis are still needed.

Figure 4: Comparisons between the predicted values and experimental flow stress curves of three testing cases: (a) Kriging; (b) RBF; (c) Polynomial; (d) ANN.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

3.3. Comparison of Metamodeling Methods

Sample size refers to the number of data points in a dataset. When sample size is sufficiently large, an accurate model can usually be constructed. However, the experiments and simulations are usually expensive or time-consuming, which may lead to the insufficiency of samples. Most of the modeling methods do not show good performance when sample size is not large enough. This could be an indication that the samples are limited for modeling method to capture the general features of the problems. The method which can provide the most satisfactory modeling performance with limited samples will thus be considered as the most effective one.

In this section, five cases with totally different temperatures and strain rates are selected for each time of model testing to assess the developed models. That is, one-fifth of the data are evenly excluded from the training data for testing each time. There are 120 different combinations of five cases being utilized. The general performances of the metamodeling methods are compared based on these 120 times of testing. In addition to [figure omitted; refer to PDF] , the standard deviation ( [figure omitted; refer to PDF] ) of [figure omitted; refer to PDF] and the correlation coefficient ( [figure omitted; refer to PDF] ) from all the testing points are employed as criteria for comparison of these four metamodeling methods: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] are the average values of [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] , respectively.

From Table 3, it can be found that Kriging model has the smallest [figure omitted; refer to PDF] , which indicates the highest prediction accuracy among the developed models. Meanwhile, the listed values of [figure omitted; refer to PDF] help to illustrate the prediction stability of models and Kriging model performs the best. Furthermore, the values of [figure omitted; refer to PDF] provide information on the strength of linear relationship between the predicted and experimental flow stress. As shown in Figure 5, the fitting results of [figure omitted; refer to PDF] clearly demonstrate that the coincidence degree between the experimental and predicted value from Kriging model is higher than that from other models.

Table 3: The predictability comparison of the models.

Model	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]	[figure omitted; refer to PDF]
Kriging	5.68%	0.0485	0.9952
RBF	7.73%	0.0760	0.9814
Polynomial	9.95%	0.1098	0.9883
ANN	9.68%	0.1898	0.9721

Figure 5: Correlations between the predicted and experimental flow stress: (a) Kriging; (b) RBF; (c) Polynomial; (d) ANN.

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

(c) [figure omitted; refer to PDF]

(d) [figure omitted; refer to PDF]

In addition to the three criteria above related to the prediction error, the distributions of all the errors are compared among the models at a macro level. As 80 points are predicted in each testing, there are a total of 9600 values of errors collected from 120 times of testing. With the errors being classified into different groups according to their values, the distributions are reflected in Figure 6, while [figure omitted; refer to PDF] -axis is relative error and [figure omitted; refer to PDF] -axis is the percentages of error distribution. It can be found that 84.91% of predicted data from Kriging method locate in the relative error range of ±10%, larger than 72.96% from RBF, 66.29% from multivariate polynomial, and 71.53% from ANN. Based on the comparative analysis on [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and the error distribution, Kriging method is proved to be the best one for modeling the material flow behavior.

Figure 6: Distributions of the prediction errors.

[figure omitted; refer to PDF]

3.4. Assessing Performance of Kriging Method

In the experiments, the deformation conditions such as temperature and strain rate are evenly distributed. However, these conditions are not always evenly distributed in practical manufacturing. Uniformity is a measure to evaluate how uniform a point set is scattered in a space. In this task, the performance of Kriging method has been further assessed when the sample uniformity is decreased. The assessment is carried out by the way that one level of temperature and one level of strain rate are successively excluded from model training data for model testing.

3.4.1. Decreasing Sample Uniformity Related to Temperature

All the data related to temperature 733 K are excluded from model training data and then used for model testing. That is, the training data cover four levels of temperature and five levels of strain rate. The comparisons between the predicted values and experimental flow stress curves at 733 K are shown in Figure 7, which are quantified by the calculated values of [figure omitted; refer to PDF] (6.56%) and [figure omitted; refer to PDF] (0.0480). It is indicated that the prediction results from Kriging model match well with the experimental data.

Figure 7: Comparisons between the predicted values and experimental flow stress curves at 733 K.

[figure omitted; refer to PDF]

3.4.2. Decreasing Sample Uniformity Related to Strain Rate

With the information related to all the temperature levels, the data related to strain rate 0.01 s^-1 are excluded from model training data for model testing. The predictability of Kriging model is also analyzed based on the [figure omitted; refer to PDF] and [figure omitted; refer to PDF] from all the testing data. The predicted values from Kriging model are all close to or on the experimental flow stress curves, and the deviations are shown in Figure 8. The high prediction accuracy is also validated by the [figure omitted; refer to PDF] 5.53% and the [figure omitted; refer to PDF] 0.0360.

Figure 8: Comparisons between the predicted values and experimental flow stress curves at 0.01 s^-1 .

[figure omitted; refer to PDF]

From the analysis above, it is found that Kriging model can be constructed with high accuracy and stability even when the sample quality is decreased from the aspects of sample size and sample uniformity. The applicability of Kriging method has been fully validated in modeling the flow behavior of 6013 aluminum alloy during hot deformation.

4. Conclusions

In order to investigate the applicability of metamodeling methods and search for the appropriate method in modeling material flow behavior, four popular metamodel types Kriging, radial basis function, multivariate polynomial, and artificial neural network are investigated and compared for modeling the elevated temperature flow behavior of 6013 aluminum alloy. These four metamodeling methods are proved to be capable of modeling the flow behavior of 6013 aluminum alloy. When sample size for modeling is decreased by 20%, the superiority of Kriging model is revealed in prediction with satisfactory accuracy and stability. Furthermore, the excellent performance of Kriging model is again demonstrated while decreasing the sample uniformity related to temperature or strain rate. It is indicated that Kriging method can be taken as an appropriate option to model the material flow behavior, even if the sample size is not large enough or the samples are not evenly distributed.

Acknowledgments

This work was supported by National Natural Science Foundation of China (no. 51475156), National Key Project of Science and Technology of China (no. 2014ZX04002071), and the Opening Foundation of Key Laboratory for Nonferrous Metal and Featured Material Processing, Guangxi Zhuang Autonomous Region (no. GXKFJ14-08).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] C. X. Wang, F. X. Yu, D. Z. Zhao, X. Zhao, L. Zuo, "Hot deformation and processing maps of DC cast Al-15%Si alloy," Materials Science and Engineering A , vol. 577, pp. 73-80, 2013.

[2] R. S. Nalawade, A. J. Puranik, G. Balachandran, K. N. Mahadik, V. Balasubramanian, "Simulation of hot rolling deformation at intermediate passes and its industrial validity," International Journal of Mechanical Sciences , vol. 77, pp. 8-16, 2013.

[3] A. Rusinek, J. R. Klepaczko, "Shear testing of a sheet steel at wide range of strain rates and a constitutive relation with strain-rate and temperature dependence of the flow stress," International Journal of Plasticity , vol. 17, no. 1, pp. 87-115, 2001.

[4] Y. C. Lin, X.-M. Chen, "A critical review of experimental results and constitutive descriptions for metals and alloys in hot working," Materials & Design , vol. 32, no. 4, pp. 1733-1759, 2011.

[5] Y. C. Lin, K. K. Li, H. B. Li, J. Chen, X. M. Chen, D. X. Wen, "New constitutive model for high-temperature deformation behavior of inconel 718 superalloy," Materials & Design , vol. 74, pp. 108-118, 2015.

[6] S. Saadatkia, H. Mirzadeh, J. M. Cabrera, "Hot deformation behavior, dynamic recrystallization, and physically-based constitutive modeling of plain carbon steels," Materials Science and Engineering A , vol. 636, pp. 196-202, 2015.

[7] W. Liu, H. Zhao, D. Li, Z. Zhang, G. Huang, Q. Liu, "Hot deformation behavior of AA7085 aluminum alloy during isothermal compression at elevated temperature," Materials Science and Engineering A , vol. 596, pp. 176-182, 2014.

[8] G.-Z. Quan, W.-Q. Lv, Y.-P. Mao, Y.-W. Zhang, J. Zhou, "Prediction of flow stress in a wide temperature range involving phase transformation for as-cast Ti-6Al-2Zr-1Mo-1V alloy by artificial neural network," Materials & Design , vol. 50, pp. 51-61, 2013.

[9] V. Senthilkumar, A. Balaji, D. Arulkirubakaran, "Application of constitutive and neural network models for prediction of high temperature flow behavior of Al/Mg based nanocomposite," Transactions of Nonferrous Metals Society of China (English Edition) , vol. 23, no. 6, pp. 1737-1750, 2013.

[10] Y. Han, G. J. Qiao, J. P. Sun, D. N. Zou, "A comparative study on constitutive relationship of as-cast 904L austenitic stainless steel during hot deformation based on Arrhenius-type and artificial neural network models," Computational Materials Science , vol. 67, pp. 93-103, 2013.

[11] G. G. Wang, S. Shan, "Review of metamodeling techniques in support of engineering design optimization," Transactions of the ASME-Journal of Mechanical Design , vol. 129, no. 4, pp. 370-380, 2007.

[12] D. Zhao, D. Xue, "A comparative study of metamodeling methods considering sample quality merits," Structural and Multidisciplinary Optimization , vol. 42, no. 6, pp. 923-938, 2010.

[13] D. G. Krige A statistical approach to some mine valuations and allied problems at the witwatersrand [M.S. thesis] , University of the Witwatersrand, Johannesburg, South Africa, 1951.

[14] G. Matheron, "Principles of geostatistics," Economic Geology , vol. 58, no. 8, pp. 1246-1266, 1963.

[15] J. Sacks, W. J. Welch, T. J. Mitchell, H. P. Wynn, "Design and analysis of computer experiments," Statistical Science , vol. 4, no. 4, pp. 409-435, 1989.

[16] S. N. Lophaven, H. B. Nielsen, J. Søndergaard, "Dace: a Matlab kriging toolbox, version 2.0,", no. IMMREP-2002-12, Technical University of Denmark, Kongens Lyngby, Denmark, 2002.

[17] M. J. D. Powell, J. C. Mason, M. G. Cox, "Radial basis functions for multivariable interpolation: a review," Algorithms for Approximation , pp. 143-167, Clarendon Press, 1987.

[18] R. H. Myers, D. C. Montgomery Response Surface Methodology: Process and Product Optimization Using Designed Experiments , John Wiley & Sons, 1995.

[19] N. Haghdadi, A. Zarei-Hanzaki, A. R. Khalesian, H. R. Abedi, "Artificial neural network modeling to predict the hot deformation behavior of an A356 aluminum alloy," Materials & Design , vol. 49, pp. 386-391, 2013.

[20] H.-Y. Li, X.-F. Wang, D.-D. Wei, J.-D. Hu, Y.-H. Li, "A comparative study on modified Zerilli-Armstrong, Arrhenius-type and artificial neural network models to predict high-temperature deformation behavior in T24 steel," Materials Science and Engineering A , vol. 536, pp. 216-222, 2012.

[21] F. J. Humphreys, M. Hatherly Recrystallization and Related Annealing Phenomena , Elsevier, 2004., 2nd.

[22] J. Zhang, H. Di, X. Wang, Y. Cao, J. Zhang, T. Ma, "Constitutive analysis of the hot deformation behavior of Fe-23Mn-2Al-0.2C twinning induced plasticity steel in consideration of strain," Materials & Design , vol. 44, pp. 354-364, 2013.

[23] G. Xiao, L.-X. Li, T. Ye, "Modification of flow stress curves and constitutive equations during hot plane compression deformation of 6013 aluminum alloy," The Chinese Journal of Nonferrous Metals , vol. 24, no. 5, pp. 1268-1274, 2014.