A New Regression Model for the Prediction of the Stress–Strain Relations of Different Materials

Experimental flow stress–strain data under different stress states are often used to calibrate the plastic constitutive model of anisotropic metal materials or identify the appropriate model that is able to reproduce their plastic deformation behavior. Since the experimental stress–strain data are discrete, they need to be mathematically returned to a continuous function to be used to describe an equivalent hardening increment. However, the regression results obtained using existing regression models are not always accurate, especially for stress–strain curves under biaxial stress loading conditions. Therefore, a new regression model is proposed in this paper. The highest-order term in the recommended form of the new model is quadratic, so the functional relationships between stress–strain components can be organized into explicit expressions. All the experimental data of the uniform deformation stage can be substituted into the new model to reasonably reproduce the biaxial experimental stress–strain data. The regression results of experimental data show that the regression accuracy of the new model is greatly improved, and the residual square sum SSE of the regression curves of the new model reduced to less than 50% of the existing three models. The regression results of stress–strain curves show significant differences in describing the yield and plastic flow characteristics of anisotropic metal materials, indicating that accurate regression results are crucial for accurately describing the anisotropic yielding and plastic flow behaviors of anisotropic metal materials.


Introduction
The application of high-strength steel, aluminum alloy, magnesium alloy and titanium alloy materials has been increasing in various manufacturing industries, such as automobiles, aerospace, bicycles, etc. [1][2][3].These metal materials generally exhibit significant anisotropy, including anisotropic yielding and anisotropic plastic flow behaviors [4][5][6].Therefore, various experimental data from different directions and different stress states has been utilized to calibrate the plastic constitutive models of these anisotropic metal materials or identify the appropriate model which are able to reproduce the plastic deformation behavior of them [7][8][9].
Uniaxial tensile test results in three directions, namely the rolling direction (RD), transverse directions (TD) and the direction 45 • from the RD of the material, respectively, and the flow stress-strain curves of biaxial tensile tests, plane-strain tension and disk compression experiments, etc., were commonly used to calibrate the plastic constitutive models of sheet metal or identify the proper model [7,[10][11][12][13][14].For instance, the plastic constitutive model of a 1.2 mm thick dual phase steel sheet with a tensile strength of 780 MPa was calibrated using the results of uniaxial tensions in RD, TD and 45 • directions and equibiaxial tension [11], and its appropriateness was identified using the results of biaxial tensile tests that were performed using the cruciform specimen.The plastic constitutive model of an AA6022-T4 sheet with 1 mm thickness was calibrated via uniaxial tension, plane-strain tension and disk compression experiments [12].Similarly, the plastic constitutive models of aluminum sheets of AA7075 with 1.5 mm thickness in two different tempers [13], a zinc-coated low carbon steel sheet (SPCD) with a thickness of 0.66 mm [14], a 590-MPa-class high-strength steel of 1.2 mm in thickness (590HSS) [15] were all determined and calibrated.
For anisotropic thin-walled metal tubes, the information of the flow stress-strain curves under different directions and different stress states are often provided via axial uniaxial tensile tests, hoop tensile tests [1,2], free bulging tests [16,17], biaxial controllable loading bulging tests [9, 18,19], pure shearing tests [20], shear-tension tests [21,22], and so on.These flow stress-strain curves are often used to calibrate the plastic constitutive model of thin-walled metal, or to identify their aappropriateness.
When using experimental stress-strain data from different directions and different stress states to calibrate the plastic constitutive model of a material or determine its applicability, each set of experimental data must correspond to the same equivalent hardening state, that is, each set of data must correspond to the same level of plastic work [23,24].Since the experimental stress-strain data are discrete, all of these experimental stress-strain data should be reproduced as continuous functions [25].The regression precision will affect the specific expression of the plastic constitutive model, and furthermore will affect the theoretical prediction accuracy of the yielding and plastic flow behaviors of materials under complex stress states.Therefore, accurate regression of experimental data is crucial.
Nonlinear regression models have been widely used to describe the nonlinear behavior of materials.Currently, experimental stress-strain data are commonly reproduced using nonlinear regression models, such as power law, the second-order function model, the fourth-order function model, and so on [25][26][27][28][29][30][31].But with the development of technology, many new materials have emerged, and many flow stress-strain data under biaxial stress state can be obtained using the new proposed test method.For the experimental stressstrain data of some new materials and some flow stress-strain data under a biaxial stress state, the regression accuracy of the existing model is not high enough [14].Therefore, proposing new regression models for these experimental data will be a fundamental issue for improving the accuracy of plastic constitutive models.
In this paper, a new regression model is proposed, and analytical and numerical methods for determining the pending parameters of the new model are introduced.Then, the experimental flow stress-strain data of three kinds of thin-walled metal materials with obvious anisotropy are used to analyze the regression feasibility and universality of the new model.The regression results of other models were also analyzed for comparison.Finally, the effects of the regression results of different models on the yielding and plastic flow properties were analyzed through the experimental results of a 0.7 mm thick cold-rolled IF steel sheet.Consequently, it was found that the new regression model in the recommended form can well reproduce most of the experimental flow stress-strain data of the three kinds of thin-walled metal materials, whether it is in the uniaxial tensile stress state or in the biaxial stress state.The regression results of different models have a significant impact on the yielding and plastic flow properties of the material in the case of the 0.7 mm thick cold-rolled IF steel sheet.

The Expression and Analysis of the New Regression Model
In order to improve the regression accuracy of experimental stress-strain data, a new regression model is proposed.
The general expression of the new model is (2) the numerator is a constant, i.e., a 1 , a 2 , • • • , a n = 0, a n+1 is a constant; (3) the order of the numerator is lower than the denominator's, i.e., n ≤ m; and (4) the order of the numerator is higher than the denominator's, i.e., m ≤ n, respectively.
The regression accuracy of the new general model is analyzed through regressing the typical biaxial experimental data σ x /σ y = 1 : 2, and the regression results obtained through the four forms of the new general model are shown in Table A1 in the Appendix A. From Table A1, it can be seen that the accuracy of the corresponding regression curves is very good (R-square = 0.999, i.e., close to one, RMSE is less than 0.55, and SSE is no more than 5.51) when n and m take any combined data in Table A1 except for m = 2 or 3, n = 0 and n = 2, m = 0.This shows that the general model shown in Equation ( 1) is reasonable, and that proper combinations of m and n can be selected as needed.In addition, it can be seen that the difference between the regression results is very small when m, n ≥ 3. R-squares are all equal to 0.99, and SSEs are essentially the same.
The regression results when m and n are both less than 3 are shown in Figure 1.All the regression results are good, and the regression coefficients are shown in Table 1.The regression coefficients are too large for m = 2, n = 1 or 2; this is not conducive to the stability of the model, so it is not recommended to choose these two functions.Furthermore, from Table A1, it can be seen that the regression result of n = 2, m = 1 is more accurate than that of n =1, m = 1.Therefore, after the comprehensive balance of regression accuracy, considering the monotonicity of the first derivative, whether it can be expressed as an explicit function of stress and strain and calculation efficiency, the combination of n = 2, m = 1 is first recommended.At this point, the corresponding specific expression of the new model is The new model shown in Equation ( 2) only has four undetermined coefficients and the highest-order term in the new model is quadratic, so the functional relationships between stress-strain components can be organized into explicit expressions.The analytic expression of the strain using the stress is given as follows: Table 1.Expressions and corresponding regression coefficients in Figure 1.Therefore, after the comprehensive balance of regression accuracy, considering the monotonicity of the first derivative, whether it can be expressed as an explicit function of stress and strain and calculation efficiency, the combination of n = 2, m = 1 is first recommended.At this point, the corresponding specific expression of the new model is

Expression Regression Coefficients
The new model shown in Equation ( 2) only has four undetermined coefficients and the highest-order term in the new model is quadratic, so the functional relationships between stress-strain components can be organized into explicit expressions.The analytic expression of the strain using the stress is given as follows: Further, the new model shown in Equation ( 2) also has the advantage that the relationship between the stress increment dσ and the strain increment dε can be explicitly expressed alongside each other, For a biaxial stress state, the equations of regression in both directions would be expressed as where dσ x and dσ y are the stress increments in the direction x and direction y, respectively.ε x and ε y are the plastic strains in the direction x and direction y, respectively.dε x and dε y are the plastic strain increments in the direction x and direction y, respectively.a 1x , a 2x , a 3x and b 1x are the pending coefficients of the function in the direction x. a 1y , a 2y , a 3y and b 1y are the pending coefficients in the direction y.
Based on Equation (5), the expression of an equivalent increment on a biaxial stress state would be expressed as where d 2 w p is the second-order differentiation of plastic work, α is the stress ratio, and α = dσ x /dσ y .If dε x is chosen as the increment step, dε x can be calculated using the ε x , ε y , σ x and α of the current step through Equation (6).Further, dσ x , dσ y and dε y can be calculated and determined sequentially using Equation (5a), stress ratio α, and Equation (5b).So far, all the stress-strain components in the next step of the two directions under biaxial stress state have been obtained.

Analytical Method for Determining the Pending Parameters of the New Model
For the new model shown in Equation (2), the pending coefficients a 1 , a 2 , a 3 and b 1 can be determined using four different sets of experimental data.If the experimental data (ε i , σ i ) (where, i = 1, 2, 3, 4) are known, then Further, Equation ( 7) can be given as Substituting four known sets of experimental data, Equation ( 8) can be expressed in a matrix form: Defining the pending coefficients vector as c s = a 1 a 2 a 3 b 1 T , Equation (9) can be given as where If matrix A is reversible, i.e., det(A) = 0, the pending coefficients of the model can be obtained as follows: For the general expression of the new model shown in Equation (1), the pending coefficients a 1 , a 2 , • • • , a n+1 and b 1 , b 2 , • • • , b m can be determined using n + m + 1 different sets of experimental data.Suppose that the experimental data (ε i , σ i ) (where i = 1, 2, • • • , n + m + 1) are known, the pending coefficient vector where ×(n + m + 1), and is a (n + m + 1)-dimensional column vector.

Numerical Method for Determining the Pending Parameters of the New Model
As we all know, there are always some test errors in experimental data due to reasons such as test equipment, experimental conditions, and experimental environment.The experimental flow stress-strain curve is always fluctuating throughout the deformation process, especially for biaxial experimental stress-strain curves under controllable loading paths.As a feedback control is required, the fluctuation is bigger.Therefore, if analytical methods are used to determine the undetermined coefficients of the new model, it means that only four experimental points are used to determine the entire regression curve, and the regression accuracy will inevitably be affected by the selected points.Thus, the leastsquares algorithm is recommended for substituting all the experimental data points in the plastic deformation stage into the expression of the new model to determine the pending coefficients, in order to minimize the impact of test errors on the regression results and obtain more accurate regression results.
The general model shown in Equation ( 1) can be represented as where T is the pending coefficient vector.If there are n + m + 1 different sets of experimental data (ε i , σ i ) (where i = 0, 1, . . . . . ., n + m + 1) known, the optimum parameter vector P will be determined through minimizing the sum of squared residuals [32], i.e., P = argmin P {M(P)} M(P) is the objective function for solving the pending coefficient vector.Usually, the sum of squared residuals is chosen as the objective function: where g(P) = [g 1 (P), g 2 (P), • • • , g N (P)] T is the residual vector under the given experimental data, and N is a positive integer.The 34] algorithm improves the iteration accuracy and stability by introducing a damping term, so it is chosen to solve the parameter vector P, which can be expressed as where µ is the damping constant which is a positive constant, and I is an m 1 × m 1 identity matrix.P k+1 and P k are the iterative values of the parameter vector in steps k + 1 and k, respectively, g(P k ) is the residual vector of step k, and J(P k ) is a Jacobin matrix of step k, which can be expressed as When the value of µ is small enough, LM algorithm is similar to the Gauss-Newton algorithm.When µ is big enough, the LM algorithm is similar to the gradient descent method.The constant µ can be self-adaptively adjusted in real time; therefore, the LM algorithm has the advantages of both the Gauss-Newton algorithm and the gradient descent method.In addition, J T c k J c k + µI is an invertible matrix, so the stability of the algorithm can be guaranteed.
The specific steps to solve the pending coefficient vector P using the LM algorithm are as follows: Step 1 Determine the value of the constant matrix η, which is used to control the convergence error; the initial value µ 0 ; the scaling factor λ, which is used to adjust the damping constant; and the initial value of pending coefficient vector P 0 .
Step 2 Calculate g(P k ), J(P k ) and F(P k ).
Step 3 Calculate the iterative step.
Step 5 if Step 7 if h k ≤ η, the iteration is terminated, and P k+1 is the result.Otherwise, let k = k + 1, and go to step 2.

Research Program
In order to verify the following advantages of the new model in regressing the flow stress-strain data of metal materials with different stress states and properties (such as regression feasibility and universality, regression accuracy, and stability), as well as the effects of the regression results on the yielding property and plastic flow behavior, the following research program is proposed.
Experimental flow stress-strain data of three kinds of materials with obvious anisotropy were used to verify regression feasibility and universality in this paper.The first kind of experimental data used in this study are from the results of the biaxial tensile test using cruciform specimens of 0.7 mm thick cold-rolled IF steel sheets, involving the data of the stress states σ x /σ y = 4:1, 2:1, 4:3, 1:4, 1:2, 3:4, 1:1 and 1:0 [7].The r values at 0 • , 45 • and 90 • (transverse direction; TD) to the rolling direction (RD) are 2.27, 1.77 and 2.65, respectively.The RD and TD directions of the material are defined as the x-and y-axes, respectively.The second kind of experimental data are from the results of controlled hydro-bulging tests of 1.9 mm thick extruded 6061F aluminum alloy tubes with a 40 mm initial outer diameter, including the data of the stress states σ z /σ θ =6 : 8, 8:7 and −4:8 [28], where σ z and σ θ are the axial and circumferential stresses, respectively.The last kind of experimental data are from the results of uniaxial tensile test in different directions of hot-rolled magnesium alloy sheet [35].
In order to illustrate the feasibility and superiority of the new model, the regression results of it were compared with the results reproduced by the second-order function model (abbreviated as SOF model) and power functions.The expression of the SOF model is where σ y and ε y are the initial yielding stress and the corresponding strain, σ max and ε max are the maximum stress and the corresponding strain, ε A and σ A are the strain and stress of experimental point A, which is located between the initial yield point and the maximum stress point, such as the 5% strain point, or some other point with a different strain value, and X 1 , X 2 , X 3 are the three pending coefficients used to address the experimental data.The power function is where K is the strength coefficient, and n is the strain hardening exponent.
The regression results of the new model were compared with results of the SOF model and power function.The regression accuracy of these models was evaluated using the residual sum of squares SSE (SSE = ∑ σ theoretical − σ experimental 2 ), the mean absolute percentage error MAPE (MAPE = ∑ σ theoretical − σ experimental • 100/σ experimental /n * ) (where n * is the number of samples, and n * = 1 for this study), and residuals of each point δ (δ = σ theoretical − σ experimental ) from overall and local angles, respectively.The regression results of the SOF model obtained separately with three different intermediate locations of point A were analyzed due to the sensitivity of the point A used for the regression [30].The point A is located at ε A = 1/4ε max , ε A = 1/2ε max and ε A = 3/4ε max , respectively.Finally, the experimental results of the 0.7 mm thick cold-rolled IF steel sheet were used to evaluate the effects of regression accuracy on the yielding and plastic flow properties of the material.

Regression Results of the New Model and Other Models
The regression universality, regression accuracy, and stability of the new model were discussed using the experimental flow stress-strain data of different stress states and material properties.

Comparative Analysis of Regression Results of Different Models
The regression results of different models are comparatively analyzed through the experimental data of the uniaxial tensile stress state.Figure 2 shows the regression results of the cold-rolled IF steel sheet with different models and the corresponding discrete experimental data.The regression parameters of these models are shown in Tables 2 and 3.  Figure 3 shows the total fitting deviation dealing with uniaxial tensile stress-strain relationships.The maximum regular residual δ of the new model is 1.38 MPa, while most of regular residual δ obtained using the SOF model and power function is greater than the value.It is obvious that better results can be obtained with the new model.It can be seen that the regression result of the new model can well reproduce the experimental data, which are significantly better than those of other models.From Table 3, it can be seen that the three coefficients of the SOF model differ greatly, and X 1 is too small.Table 4 shows the overall error SSEs and MAPE of regression results using the new model; the SOF model with ε A = 1/4ε max , ε A = 1/2ε max and ε A = 3/4ε max ; and the power function.Regardless of SSE or MAPE, the result of the new model is the smallest, closer to the ideal case of SSE = 0, followed by the SOF model, with ε A = 1/4ε max ; the other two errors corresponding to the SOF model with ε A = 1/2ε max and ε A = 3/4ε max are much larger.The SSE of the regression result obtained using the power function is the maximum.The SSE of the new model is only 6.5% of the SOF model with ε A = 1/4ε max .The predicting accuracy of the new model is better than that of the SOF model, and much better than the power function.In addition, the highest-order term in the new model is quadratic, so the functional relationships between stress-strain components can be easily organized into explicit expressions.Figure 3 shows the total fitting deviation dealing with uniaxial tensile stress-strain relationships.The maximum regular residual δ of the new model is 1.38 MPa, while most of regular residual δ obtained using the SOF model and power function is greater than the value.It is obvious that better results can be obtained with the new model.Figure 3 shows the total fitting deviation dealing with uniaxial tensile stress-strain relationships.The maximum regular residual δ of the new model is 1.38 MPa, while most of regular residual δ obtained using the SOF model and power function is greater than the value.It is obvious that better results can be obtained with the new model.Furthermore, the results in Tables 3 and 4 show that the regression results of the SOF model were significantly affected by the value of point A. Among the three selected point A, the SSE is the smallest when εA = 1/4εmax, and the biggest when εA = 3/4εmax.Furthermore, the results in Tables 3 and 4 show that the regression results of the SOF model were significantly affected by the value of point A. Among the three selected point A, the SSE is the smallest when ε A = 1/4ε max , and the biggest when ε A = 3/4ε max .

Regression Results of Biaxial Stress State
Figure 4 shows the regression results of stress ratio σ x /σ y = 1 : 2 of cold-rolled IF steel sheets with different models, and the corresponding discrete experimental data.The corresponding regression parameters of these models are shown in Tables 5 and 6.From Table 6, it also can be seen that X 1 is too small, and X 1 and X 3 vary widely, which is not conducive to the stable calculation of the model.ε A has a significant effect on the coefficient X 2 , which will lead to different regression results.
coefficient X2, which will lead to different regression results.From Figure 4, it also can be seen that the regression accuracy of the SOF model is sensitive to the point A, which is an unstable regression model.Different points A are selected, the results vary widely.When the selected point A is unreasonable, the regression accuracy decreases significantly, and the SSE rises sharply.The regression results obtained using the power function significantly deviate from the experimental data.Moreover, there are no longer physical meanings to the regression results of many biaxial flow stress-strain data using power law, unlike the regression results of uniaxial tensile stress-strain data.
Table 7 shows the overall error SSEs and MAPE of the new model, the SOF model with   From Figure 4, it also can be seen that the regression accuracy of the SOF model is sensitive to the point A, which is an unstable regression model.Different points A are selected, the results vary widely.When the selected point A is unreasonable, the regression accuracy decreases significantly, and the SSE rises sharply.The regression results obtained using the power function significantly deviate from the experimental data.Moreover, there are no longer physical meanings to the regression results of many biaxial flow stress-strain data using power law, unlike the regression results of uniaxial tensile stress-strain data.
Table 7 shows the overall error SSEs and MAPE of the new model, the SOF model with ε A = 1/4ε max , ε A = 1/2ε max and ε A = 3/4ε max , and the power function along direction y.For the biaxial stress state, regardless of SSE or MAPE, the result of the new model is still the smallest, followed by the SOF model with ε A = 1/2ε max , which has a significant regression advantage.Furthermore, all the experimental points in the plastic deformation stage can be used to optimize the regression results of experimental curves through this new model, using numerical algorithms; thus, the experimental data are evenly distributed on both sides of the regression curve, as well as the calculation error of the plastic work up and down offset.The regression results are reasonable and stable.
Table 7. Overall error of biaxial stress-strain curves with σ x /σ y = 1 : 2 along direction y. Figure 5 presents the results of the deviation analysis.The residuals reproduced by the new model in the whole plastic deformation process are evenly distributed on both sides, with a maximum stress deviation of no more than 2.5 MPa.The predicting accuracy of this suggested model looks better than that of the SOF model, and much better than that of the power function, while the regular residuals reproduced by the SOF model with ε A = 1/4ε max , ε A = 1/2ε max and ε A = 3/4ε max are different.Compared with the other two values, the residual distribution is more uniform when ε A = 1/2ε max .The regression residuals are biased to the positive side when ε A = 3/4ε max ; the regression curve is always lower than the experimental results., which has a significant regression advantage.Furthermore, all the experimental points in the plastic deformation stage can be used to optimize the regression results of experimental curves through this new model, using numerical algorithms; thus, the experimental data are evenly distributed on both sides of the regression curve, as well as the calculation error of the plastic work up and down offset.The regression results are reasonable and stable.Figure 5 presents the results of the deviation analysis.The residuals reproduced by the new model in the whole plastic deformation process are evenly distributed on both sides, with a maximum stress deviation of no more than 2.5 MPa.The predicting accuracy of this suggested model looks better than that of the SOF model, and much better than that of the power function, while the regular residuals reproduced by the SOF model with  Figure 6 shows the regression results of various stress ratios of cold-rolled IF steel sheets with different models and the corresponding discrete experimental data.The matched regression parameters of these models are shown in Table 8.From Figure 6, it can be seen that there are significant differences in the experimental flow stress-strain curves under different stress states.However, the regression results of the new model all well reproduce the experimental data of these stress ratios.The SSEs of the new model and the power function are shown in Tables 2 and 3. Table 2 shows the SSEs in the direction x, and Table 3 shows the SSEs in the direction y.Table 9 shows the SSE of the Figure 6 shows the regression results of various stress ratios of cold-rolled IF steel sheets with different models and the corresponding discrete experimental data.The matched regression parameters of these models are shown in Table 8.From Figure 6, it can be seen that there are significant differences in the experimental flow stress-strain curves under different stress states.However, the regression results of the new model all well reproduce the experimental data of these stress ratios.The SSEs of the new model and the power function are shown in Tables 2 and 3. Table 2 shows the SSEs in the direction x, and Table 3 shows the SSEs in the direction y.Table 9 shows the SSE of the new model is only 10 −5 ∼ 10 −6 times that of the power function, and the regression results of the biaxial stress states have been significantly improved.The regression feasibility and universality of these models are briefly analyzed using the following experimental flow stress-strain data, including the experimental data of an extruded 6061F aluminum alloy tube obtained through controlled hydro-bulging tests, and the experimental data of hot-rolled magnesium alloy sheet measured via the uniaxial tension in different directions.

Models
Figure 7 displays the experimental flow stress-strain curves of an extruded 6061F aluminum alloy tube with stress ratios −4:8, 6:8 and 8:7 [28] measured through controlled hydro-bulging tests and the regression results of different models.It can be seen that there are significant differences in the experimental flow stress-strain curves under different stress states.However, they all can be well reproduced, using a numerical algorithm, through the new model.tension in different directions.
Figure 7 displays the experimental flow stress-strain curves of an extruded 6061F aluminum alloy tube with stress ratios −4:8, 6:8 and 8:7 [28] measured through controlled hydro-bulging tests and the regression results of different models.It can be seen that there are significant differences in the experimental flow stress-strain curves under different stress states.However, they all can be well reproduced, using a numerical algorithm, through the new model.Figure 8 shows the experimental flow stress-strain data of a hot-rolled magnesium alloy sheet measured via uniaxial tensile tests in different directions, and the corresponding results reproduced by different models.It can be seen that they can all be well reproduced, using a numerical algorithm, through the new model.The regression results of the new model are significantly better than those of the power function.Figure 8 shows the experimental flow stress-strain data of a hot-rolled magnesium alloy sheet measured via uniaxial tensile tests in different directions, and the corresponding results reproduced by different models.It can be seen that they can all be well reproduced, using a numerical algorithm, through the new model.The regression results of the new model are significantly better than those of the power function.

The Effect of the Regression Results of Different Models on the Yielding Property
For anisotropic materials, in order to obtain accurate FE simulation results of plastic deformation, the yielding property of the material should be accurately reproduced via the plastic constitutive model used in the FE analysis.Therefore, the effects of the regression results of different models on the yielding and plastic flow properties were analyzed through the experimental results of a 0.7 mm thick cold-rolled IF steel sheet.

Effect of Regression Accuracy on the Yield Loci
First, the initial yield points of different stress states were determined according to the equivalent strain calculated through the Mises function, which is 0.002.Assuming that the material is incompressible, the initial yield condition can be expressed as where ε x , ε y are the true strains along the directions x and y, respectively.

The Effect of the Regression Results of Different Models on the Yielding Property
For anisotropic materials, in order to obtain accurate FE simulation results of plastic deformation, the yielding property of the material should be accurately reproduced via the plastic constitutive model used in the FE analysis.Therefore, the effects of the regression results of different models on the yielding and plastic flow properties were analyzed through the experimental results of a 0.7 mm thick cold-rolled IF steel sheet.

Effect of Regression Accuracy on the Yield Loci
First, the initial yield points of different stress states were determined according to the equivalent strain calculated through the Mises function, which is 0.002.Assuming that the material is incompressible, the initial yield condition can be expressed as where x y ε ε are the true strains along the directions x and y, respectively.Further, experimental points on the subsequent yield loci can be obtained according to ( ) where ( ) Figure 9 shows the yield loci of the cold-rolled IF steel sheet, which were obtained with the regression results of the new model and the power function, respectively.Figure 9a-d are the yield loci corresponding to 0 0.01, 0.015, 0.02 p ε = and 0.03, respectively.The yield loci are obtained via the least-squares method based on the Hill48′s quadratic yield function [36]: Further, experimental points on the subsequent yield loci can be obtained according to x and σ x are the true stress for the steps i and i-1 along the direction x, σ are the true strains for the steps i and i-1 along the direction y.The expression of f −1 x or y (σ) see Equation (3), dσ i 0 , and dε i 0 are the stress increment and the strain increment from uniaxial tensile testing, respectively.
Figure 9 shows the yield loci of the cold-rolled IF steel sheet, which were obtained with the regression results of the new model and the power function, respectively.Figure 9a-d are the yield loci corresponding to ε p 0 = 0.01, 0.015, 0.02 and 0.03, respectively.The yield loci are obtained via the least-squares method based on the Hill48 s quadratic yield function [36]: where f σ ij represents the yield function of the material element, and A, B, C, D are the pending coefficients, which are determined using the experimental data.
It can be seen that there are significant differences in the yield loci of the cold-rolled IF steel sheet obtained based on the regression results of different models.The regression results have a significant impact on the yield loci of the cold-rolled IF steel sheet.The yield loci of the cold-rolled IF steel sheet obtained based on the new model are significantly smaller than those obtained based on the power function.The smaller the degree of deformation is, the greater the difference in the yield loci will be.The biggest difference between the two shapes is nearly 10%, as shown in Figure 10.

Effect of the Yield Loci on the Predicted Plastic Flow Direction
If the associated flow rule is adopted, then the flow direction of the material will be predicted by the yield function.As mentioned above, the regression results of different models have a significant impact on the yield loci, so will inevitably have a significant impact on the calculated flow direction, as shown in Figure 11.The regression results of different models directly affect the yield loci and plastic flow directions predicted using theory.Higher theoretical prediction accuracy for the yielding and plastic flow properties can be obtained based on a better regression result.

The Effect of the Regression Results of Different Models on the Plastic Flow Property
The influence of the regression results of different models on the plastic flow properties was also analyzed through the experimental results of 0.7 mm thick cold-rolled IF steel sheets.The plastic flow properties of thin-walled metal materials are usually quantitatively described using the ratio of the two in-plane strain increments [7].It can be seen that there are significant differences in the yield loci of the cold-rolled IF steel sheet obtained based on the regression results of different models.The regression results have a significant impact on the yield loci of the cold-rolled IF steel sheet.The yield loci of the cold-rolled IF steel sheet obtained based on the new model are significantly smaller than those obtained based on the power function.The smaller the degree of deformation is, the greater the difference in the yield loci will be.The biggest difference between the two shapes is nearly 10%, as shown in Figure 10.

Effect of the Yield Loci on the Predicted Plastic Flow Direction
If the associated flow rule is adopted, then the flow direction of the material will be predicted by the yield function.As mentioned above, the regression results of different models have a significant impact on the yield loci, so will inevitably have a significant impact on the calculated flow direction, as shown in Figure 11.The regression results of different models directly affect the yield loci and plastic flow directions predicted using theory.Higher theoretical prediction accuracy for the yielding and plastic flow properties can be obtained based on a better regression result.

Effect of the Yield Loci on the Predicted Plastic Flow Direction
If the associated flow rule is adopted, then the flow direction of the material will be predicted by the yield function.As mentioned above, the regression results of different models have a significant impact on the yield loci, so will inevitably have a significant impact on the calculated flow direction, as shown in Figure 11.The regression results of different models directly affect the yield loci and plastic flow directions predicted using theory.Higher theoretical prediction accuracy for the yielding and plastic flow properties can be obtained based on a better regression result.The real-time ratio of the two in-plane strain increments in the plastic deformation process can be obtained based on the regression results of the experimental flow stress-strain data of the two directions x and y.Assuming that thin-walled metal materials are uncompressible bodies, the ratio of the two in-plane strain increments can be written as

The Effect of the Regression Results of Different Models on the Plastic Flow Property
where ε x and ε y are the strain increments in the directions x and y, respectively.Based on Equation ( 4), the strain increments can be determined using the corresponding stress increments: where a ix , b jx and a iy , b jy (i = 1, 2, 3, j = 1) are the pending coefficients in the directions x and y, respectively.
By substituting Equation (24) into Equation ( 23), the real-time ratio of the two in-plane strain increments can be represented by the regression functions in two directions.For comparison purposes, the strain increments in the directions x and y are also determined based on the following power functions: where K x , K y and n x , n y , are the pending parameters of the power functions in the directions x and y.
Figure 12 shows that the ratio of the two in-plane strain increments varies with the strain value in the tensile direction in real time.The experimental data of β-value are also provided in Figure 12.In Figure 12, there is a big difference in the ratios β obtained based on the regression results of different models, and the ratios β calculated on the basis of the regression results of the new model are closer to the experimental results.But the β-values obtained based on the power function are far from the experimental results, which will lead to a large error in the description of the plastic flow property of the material.In addition, the ratio β varies with the hardening process, which is similar to the results in much of the literature [37].But the variation in the ratio β with the hardening process is also different for different loading paths.The ratio β increases with the hardening process for the stress ratios α = 1:4, 2:1 and 3:4, while in contrast, the ratio β decreases with the hardening process for α = 4:1, 1:2 and 4:3.This means that a poor predicting accuracy will be obtained if a constant β-value is used in the plastic constitutive model of FE simulations, like the value currently used in stamping simulations.An accurate regression model is crucial for describing the plastic flow properties of thin-walled metal materials.

Conclusions
In this paper, a new regression model was proposed, and the analytical and numerical methods for determining the pending parameters of the new model were introduced.Then, the experimental flow stress-strain data of three kinds of thin-walled metal materials with obvious anisotropy were used to analyze the regression feasibility and universality of the new model.The regression results of other models were also analyzed for comparison.Finally, the effects of the regression results of different models on the yielding and plastic flow properties were analyzed through the experimental results of 0.7 mm thick cold-rolled IF steel sheets; the following conclusions were obtained.

1.
The new regression model has a general expression, the regression results of which are all good, except for several combinations of m and n.The combination of n = 2, In Proceedings of the m = 1 is first recommended after comprehensive balancing of regression accuracy, due to the monotonicity of the first derivative, whether it can be expressed as an explicit function of stress and strain and calculation efficiency.

2.
The recommended form of the new model only has four pending coefficients, and the highest-order term in the recommended form of the new model is quadratic, so the functional relationships between stress-strain components can be organized into explicit expressions using this model, which is easy to use for the determination of the coefficients of plastic constitutive models.

3.
The new regression model in the recommended form can well reproduce most of the experimental flow stress-strain data of the three kinds of thin-walled metal materials, whether it is in the uniaxial tensile stress state or in the biaxial stress state.

4.
The regression results of different models have a significant impact on the yielding properties of 0.7 mm thick cold-rolled IF steel sheets.The yield loci of the cold-rolled IF steel sheets obtained based on the new model are smaller than those obtained based on the power function.The biggest difference between the two shapes is nearly 10%.

5.
The regression results of different models have a significant impact on the plastic flow properties of the material, in the case of 0.7 mm thick cold-rolled IF steel sheets.The ratios β, calculated on the basis of the regression results of the new model, are much closer to the experimental results., where n * is the number of measurements).The smaller the RMSE is, the better the regression result will be.R-square ranges from 0 to 1; the closer to one the R-square is, the better the regression result will be.R-square = 1 is an ideal case, i.e., the regression results are completely coincident with the experimental results.

where, a 1 ,
a 2 , • • • , a n , a n+1 and b 1 , • • • , b m−1 , b m are the pending coefficients, n and m are nonnegative integers, n ≥ 0, m ≥ 0, σ is the flow stress, and ε is the plastic strain corresponding to σ. n and m have different combinations.When n and m take different values, the forms of functions are different.They can be classified into four forms which are as follows.(1) Pure polynomial (i.e., the denominator is a constant), i.e., b 1 , • • • , b m−1 = 0, b m is a constant;

Figure 2 .
Figure 2. Regression results of cold-rolled IF steel sheets with different models, and the corresponding discrete experimental data.

Figure 2 .
Figure 2. Regression results of cold-rolled IF steel sheets with different models, and the corresponding discrete experimental data.

Figure 2 .
Figure 2. Regression results of cold-rolled IF steel sheets with different models, and the corresponding discrete experimental data.

Figure 4 .
Figure 4. Experimental stress-strain relationships of a cold-rolled IF steel sheet with a stress ratio of 1: 2 x y σ σ = , reproduced with different functions.

Figure 4 .
Figure 4. Experimental stress-strain relationships of a cold-rolled IF steel sheet with a stress ratio of σ x /σ y = 1 : 2, reproduced with different functions.
curve is always lower than the experimental results.

Figure 7 .
Figure 7. Experimental flow stress-strain curves of extruded 6061F aluminum alloy tubes measured via controlled hydro-bulging tests and the regression results of different models.(a) 4 : 8 α =

Figure 8 .
Figure 8. Uniaxial tensile stress-strain curves of hot-rolled magnesium alloy sheets in different directions, and the corresponding results reproduced by different models.(a) Strething in rolling direction.(b) Strething in transverse direction.
the true stress for the steps i and i-1 along the direction x, the true stress for the steps i and i-1 along the direction y, strains for the steps i and i-1 along the direction x, and strains for the steps i and i-1 along the direction y.The expression of ( ) increment and the strain increment from uniaxial tensile testing, respectively.

Figure 8 .
Figure 8. Uniaxial tensile stress-strain curves of hot-rolled magnesium alloy sheets in different directions, and the corresponding results reproduced by different models.(a) Strething in rolling direction.(b) Strething in transverse direction.
stress for the steps i and i-1 along the direction y, f −1 strains for the steps i and i-1 along the direction x, and f −1

Figure 9 .
Figure 9. Yield loci of the cold-rolled IF steel sheet which were obtained with the regression results of the new model and the power function ( 0 p ε is the plastic strain of the uniaxial tensile test). (a)

Figure 9 .
Figure 9. Yield loci of the cold-rolled IF steel sheet which were obtained with the regression results of the new model and the power function (ε p 0 is the plastic strain of the uniaxial tensile test).(a) ε p 0 = 0.01; (b) ε p 0 = 0.015; (c) ε p 0 = 0.02; (d) ε p 0 = 0.03; (e) Yield loci corresponding to different values of ε p 0 .

Figure 10 .
Figure 10.Shape difference in the yield loci obtained based on the regression results of the new model and the power function.

Figure 11 .
Figure 11.Predicted plastic flow direction in terms of the yield loci shown in Figure 9b.(a) Different models can have an impact on regression results; (b) Direction of flow.

Figure 10 .
Figure 10.Shape difference in the yield loci obtained based on the regression results of the new model and the power function.

Figure 10 .
Figure 10.Shape difference in the yield loci obtained based on the regression results of the new model and the power function.

Figure 11 .
Figure 11.Predicted plastic flow direction in terms of the yield loci shown in Figure 9b.(a) Different models can have an impact on regression results; (b) Direction of flow.

Figure 11 .
Figure 11.Predicted plastic flow direction in terms of the yield loci shown in Figure 9b.(a) Different models can have an impact on regression results; (b) Direction of flow.

Table 2 .
Regression parameters of the new model and power function for uniaxial tensile stressplastic strain curve.

Table 3 .
Regression coefficients of the SOF model for uniaxial tensile stress-plastic strain curve.

Table 4 .
Overall error of each model for uniaxial tensile stress-plastic strain curve.

Table 5 .
Regression parameters of the new model and power function for biaxial stress-strain curves with 1: 2

Table 6 .
Regression coefficients of the SOF model for biaxial stress-strain data of 1 : 2

Table 5 .
Regression parameters of the new model and power function for biaxial stress-strain curves with σ x /σ y = 1 : 2.

Table 6 .
Regression coefficients of the SOF model for biaxial stress-strain data of σ x /σ y = 1 : 2.
For the biaxial stress state, regardless of SSE or MAPE, the result of the new model is still the smallest, followed by the SOF model with

Table 7 .
Overall error of biaxial stress-strain curves with

Table 8 .
Regression parameters of the new model and power function for biaxial stress-strain curves with different stress ratios.

Table 9 .
SSEs of the new model and power function for biaxial stress-strain curves with different stress ratios.
: SSE is the residual sum of squares (SSE = ∑ σ theoretical − σ experimental 2 ).The smaller the SSE is, the better the regression result will be.SSE = 0 is an ideal case, i.e., the regression results are completely coincident with the experimental results.RMSE is the root-mean-square error (RMSE = ∑ σ theoretical − σ experimental 2 /n * Note