Modeling Material Flow Behavior during Hot Deformation Based on Metamodeling Methods

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.Themerits 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.


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][5][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][9][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.

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: where  is the number of regression functions,   (⋅) is a regression function,   is the coefficient for   (⋅), x is the design point, and (⋅) is the random process function.The correlation function is defined by the following equation: where Cov(⋅) is the correlation function, x 1 and x 2 are two design points, and  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].

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 (x 1 , x 2 , . . ., x  ) are chosen first from the design points (x 1 , x 2 , . . ., x  ) based on some criteria.Then the basis functions are constructed by using these center points as where x is the design point, x  is the center point for   (⋅), and ‖x −x  ‖ is the Euclidian distance between the two points.The function (‖ ⋅ ‖) 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: where ĝ is the estimated output value,   (⋅) is a basis function,  0 is a constant, and   is the coefficient for   (⋅).

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 (x = [ 1 ,  2 , . . .,   ]) and their interactions such as A multivariate polynomial is constructed as the weighted sum of these basis functions: where ĝ is the estimated output value,   (⋅) is a basis function,  0 is a constant, and   is the coefficient for   (⋅).

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.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.

Model Construction and Discussion
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 (), strain rate ( ε ), and strain () whereas flow stress () is the output.In addition, the parameters of temperature and strain rate are processed by natural logarithm.

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 ( ave ) and the maximum ( max ) values of absolute relative errors () from each testing are employed as quantitative criteria to evaluate the performance of models and the applicability of modeling methods: where  is the experimental value,  is the predicted value, and  is the number of testing points.As shown in Table 2, the detailed information on  ave and  max from the developed models is provided.Except few specific conditions, the values of  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 True stress  (MPa)  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.

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  ave , the standard deviation () of  and the correlation coefficient () from all the testing points are employed as criteria for comparison of these four metamodeling methods: where , , and  are the average values of   ,   , and   , respectively.
From Table 3, it can be found that Kriging model has the smallest  ave , which indicates the highest prediction accuracy  among the developed models.Meanwhile, the listed values of  help to illustrate the prediction stability of models and Kriging model performs the best.Furthermore, the values of  provide information on the strength of linear relationship between the predicted and experimental flow stress.As shown in Figure 5, the fitting results of  clearly demonstrate that the coincidence degree between the experimental and predicted value from Kriging model is higher than that from other models.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 -axis is relative error and -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  ave , , , and the error distribution, Kriging method is proved to be the best one for modeling the material flow behavior.

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.

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  ave (6.56%) and  (0.0480).It is indicated that the prediction results from Kriging model match well with the experimental data.

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  ave and  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  ave 5.53% and the  0.0360.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.

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.

3. 1 .
Experiments and Modeling Preparation.The material used in the present experiment is 6013 aluminum alloy, with

Figure 1 :Figure 2 :
Figure 1: The architecture of the ANN model.

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

Figure 6 :
Figure 6: Distributions of the prediction errors.

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

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

Table 3 :
The predictability comparison of the models.