Next Article in Journal
The Effect of Cooling Rate from Solution Treatment on γ′ Reprecipitates and Creep Behaviors of a Ni-Based Superalloy Single-Crystal Casting
Next Article in Special Issue
β-Phase-Induced Quasi-Cleavage Fracture Mechanism by Dual-Phase High-Strength Titanium Alloy at Elevated Temperature
Previous Article in Journal
Evaluation of the Skin Permeation-Enhancing Abilities of Newly Developed Water-Soluble Self-Assembled Liquid Crystal Formulations Based on Hexosomes
Previous Article in Special Issue
Modeling the Mechanical Properties of Heat-Treated Mg-Zn-RE-Zr-Ca-Sr Alloys with the Artificial Neural Network and the Regression Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Microstructure and Mechanical Properties of Heat Treated ZE41-Ca-Sr Alloys for Integrated Computing Platform

1
Ningbo Branch of China Academy of Ordnance Science, Ningbo 315103, China
2
Key Laboratory of Solidification Control and Digital Preparation Technology, School of Materials Science and Engineering, Dalian University of Technology, Dalian 116024, China
*
Authors to whom correspondence should be addressed.
Crystals 2022, 12(9), 1237; https://doi.org/10.3390/cryst12091237
Submission received: 17 August 2022 / Accepted: 23 August 2022 / Published: 1 September 2022
(This article belongs to the Special Issue High-Performance Light Alloys 2022)

Abstract

:
The main objective of this study is to present a methodology to model the microstructure and mechanical properties of ZE41-xCa-ySr alloys for integrated optimization calculation of the heat treatment process of gearbox casting. Firstly, the models of microstructure and mechanical properties of ZE41-xCa-ySr alloys (0 ≤ x ≤ 2, 0 ≤ y ≤ 0.2) are developed using an artificial neural network (ANN) and multivariate regression. The dataset for ANN and regression models is generated by investigating the microstructures and mechanical properties of the ZE41-xCa-ySr alloys. The inputs for ANN and regression models are Ca and Sr contents, aging temperature and aging time. The outputs are grain size, ultimate tensile strength, elongation and microhardness. The optimal ANN model is obtained by testing the performance of different network architectures. In addition, multivariate regression models have been built based on the Least Squares method. Secondly, based on SiPESC software, an Integrated Computing Platform is constructed by combining the scripting language with the command line operation of simulation software, realizing the “process—microstructure—property” optimization calculation. Finally, based on the developed regression model, an Integrated Computing Platform batch called MATLAB achieves the heat treatment process optimization based on mechanical property prediction. The optimum aging temperature of the ZE41-0.17Ca-0.2Sr alloy is 322 °C, and the corresponding aging time is 11 h. Furthermore, the optimized results are validated by the ANN model, suggesting that ANN predicted results are in good agreement with optimized results. As a consequence, this work provides a new strategy for the research and development of Mg alloys, contributing to acceleration in the development of magnesium alloys.

1. Introduction

Magnesium (Mg) has the potential to provide 60–75% mass savings relative to steel or cast iron, and die casting can be used to economically produce large and thin-wall Mg alloy castings to replace steel subsystems [1]. Mg alloys containing rare-earth (RE) elements are considered to be promising alternatives to commercially available alloys, due to their good oxidation resistance and high strength [2]. In addition, mechanical properties of Mg alloys can be improved by employing relevant heat treatment and other engineering processes [3,4,5]. The processes of tailoring the alloy to deliver application-dependent property requirements play key roles in accelerating the conversion of improved Mg alloys into products.
In 2008, Integrated Computational Materials Engineering (ICME) was first proposed, which involved the tools for modeling and simulation, as well as design optimization. ICME aims to reduce the time to market of innovative products by exploiting the design and development of material and process path [6]. Since 2012, ICME has been gradually applied to the research and development of alloys and their components. This approach integrates the multiscale models and key experiments into the entire process of aluminum alloys development. Therefore, ICME promotes the research and development mode to change from “Traditional Design” to “Predictive Design”.
Designing materials for targeted property as required in ICME demands modeling and simulation, which includes both process–microstructure and microstructure–property relations. However, it is a difficult task to briefly summarize a range of multiscale modeling approaches for all process–microstructure and microstructure–property relations for all the alloys.
A few experimental research efforts were devoted to studying the strengthening mechanism, heat treatment technology and microstructure evolution of the Mg-Zn-RE-Zr alloys [7,8,9,10]. The relationship between process parameters and mechanical properties for Mg-Zn-RE-Zr alloys has only been studied empirically. The “composition/process-microstructure-property” relation that is foundational to the design of the alloy composition and process has not yet been studied. The few systematic efforts that have been made to analyze trends in experimental data as a basis for predictions have, in large part, been inconclusive, not the least of which is due to the lack of large amounts of organized data, and even more importantly the challenge of sifting through them in an efficient manner.
In recent decades, artificial neural networks (ANNs) and various regression methods have been powerful and flexible modeling tools that can lead to significant improvements in materials science, modeling complex problems and exploring the correlations between different materials’ properties [11,12,13,14]. Yang employed the ANN model with a backpropagation (BP) algorithm to explore correlations that exist between heat treatment processes and mechanical properties of A357 alloy [15]. Conduit developed an ANN to enable the prediction of individual materials properties both as a function of the composition and heat treatment routine [16]. In addition, linear and nonlinear regression models were applied to develop the input–output relationship in many casting processes [17,18]. Chen applied the nonlinear mathematical model to quantitatively analyze the effects of the heat treatment on Vickers hardness of the Al-Si-Mg alloy, and achieved the optimum heat treatment process using the sequential approximation optimization method [19]. Akhtar optimized process parameters in the turning of Aluminum Alloy 7075 using a Computer Numerical Control machine [20].
ZE41 magnesium alloy is one of the most popular magnesium alloys in the Mg-Zn-RE-Zr based alloy and has been widely used for the aircraft gearbox and generator housings on military helicopters [21,22,23]. In this study, we take ZE41-xCa-ySr alloy as a research object. Firstly, we focus on modeling and correlating the microstructure and mechanical properties of the heat-treated ZE41-Ca-Sr alloys, with alloying elements and aging process using the ANN approach and regression models. Secondly, based on the developed regression model facilitating Integrated Computing Platform, the heat treatment process and alloying composition of ZE41-xCa-ySr alloy are optimized by the Integrated Computing Platform. Finally, the optimized results are validated by the ANN model. This work aims to provide a new strategy for the research and development of Mg alloys, which is conducive to accelerating the development of magnesium alloys.

2. Experiment and Methods

2.1. Experimental Procedure

Figure 1 illustrates the flow chart of this work. In the present work, ZE41-xCa-ySr (Mg-4.2Zn-1.7RE-0.8Zr-xCa-ySr x = 0, 0.2 wt.%, y = 0, 0.2 wt.%) was used as the casting material due to its advantages such as excellent fluidity, good heat resistance and low wall-thickness effect. The casting ingots were produced by high-purity Mg, Zn, Ce-rich mischmetal (50 wt.% Ce, 28 wt.% La, 16 wt.% Nd, 4% wt.% Pr and 2 wt.% impurity), Mg-30Zr, Mg-20Ca and Mg-20Sr in an electric resistance furnace under a protected argon atmosphere at 730 °C. Then, the samples were subjected to different heat treatments. Experiments were designed and accomplished according to underage, peak age and overage conditions at different temperatures. The chemical compositions of the as-cast alloys were determined by the Inductive Coupled Plasma (ICP) method and the results were summarized in Table 1. The RE element consisted of Ce, La, trace Nd and Pr. The heat-treatment parameters (T1 conditions) used in this work were as follows: artificial aging temperatures (300 °C, 325 °C and 350 °C) and aging time (0–32) h. After heat treatments, the tensile specimens were machined with 30 mm in gauge length and 6 mm in diameter. Based on the standard GB/T 228.1-2010, tensile tests were conducted at room temperature. Ultimate tensile strength (UTS) and elongation to failure (El.) were measured. Furthermore, Vickers hardness tests were performed with a 1 kg load. Next, 10 indentations per sample were analyzed to improve precision. The average value was reported as the microhardness (HV). In addition, microstructural examinations were conducted using an optical microscope (OM). The specimens for OM observations were etched in a solution of 2.1 g picric acid, 5 mL acetic acid, 5 mL H2O and 35 mL ethanol. Average grain size was measured using the linear intercept method described in ASTM standard E 112-88.
Based on the previous study, the optical micrographs of the as-cast ZE41-xCa-ySr alloys are shown in Figure 2 [24]. It can be seen that the as-cast alloys consist of equiaxed α-Mg matrix and eutectic compounds distributed along grain boundaries. The average grain size of α-Mg phase decreases with increasing Ca and Sr. Figure 3 presents the optical micrographs of the peak-aged ZE41-xCa-ySr alloys, wherein the microstructures consist of primary α-Mg grains and interdendritic phases dispersing along the grain boundaries, as well as the intermetallic particle existing inside grains [25]. The average grain sizes for the peak-aged ZE41 and ZE41-0.2Ca-ySr (y = 0, 0.1 and 0.2 wt.%) alloys are (49.6 ± 2.7) μm, (37.5 ± 2.4) μm, (34.8 ± 2.1) μm, (32.5 ± 2.5) μm and (35.2 ± 2.1) μm, respectively, revealing a slight growth of α-Mg grains after peak age at 325 °C compared with the as-cast alloys. Therefore, the average grain size is selected as an output variable to express the microstructure.

2.2. ANN Modeling

An ANN is a mathematical model consisting of many highly interconnected processing elements organized into layers. ANN keeps knowledge with connection weights [26]. Input signals are processed in the course of the hidden layer, and the following outputs are computed in the output layer. Input–output pairs are presented to the ANN and weights are adjusted to minimize the error between the predicted output and actual value. The multilayered neural network (MLP) is utilized to establish an ANN model, which is used to predict the microstructure and mechanical properties of ZE41-xCa-ySr alloys. Since the Back Propagation (BP) algorithm is a representative method used to reduce the errors created by the gradient descent method, it is used to train the multilayer feed forward network. A three-layered ANN model is developed using the MATLAB R2018a. Figure 4 shows the network configuration of the ANN model. As it seen in Figure 2, the network structure includes three layers: an input layer (four neurons), one hidden layer and an output layer (four neurons). The input variables consist of Ca, Sr content, aging temperature (Ta) and aging time (ta). The outputs correspond to average grain size (D), ultimate tensile strength (UTS), elongation (El.) and microhardness (HV). The dataset has been generated from a series of experiments on samples, as discussed in Section 2.1, which are listed in Appendix A. In order to avoid overfitting in the ANN training, 120 groups of data are divided randomly into three subsets: 70% training set, 15% testing set and 15% validation set. Both input variables and output variables are fed into the neural network toolbox.
The ANN model training parameters are presented in Table 2. Based on the training parameters, the mathematical model of the ANN is shown in Figure 5. The data are normalized to reduce the order of magnitude difference between the various dimensions [27]. As activation transfer functions, the hyperbolic tangent “tan-sigmoid” function and the linear transfer function “Purelin” are used. Compared with the standard gradient descent algorithm, the Levenberg–Marquardt (LM) algorithm possesses fast convergence and a small mean square error [26]. As a result, a BP neural network with LM algorithm is adopted to train the ANN model. In addition, the number of hidden neurons is estimated according to the empirical formula:
M = n + m + a
where m and n are the number of neurons in the input and output layers, respectively, and a is a constant in the range of 1 to 10 [28,29]. According to the empirical formula, nine kinds of ANN architectures with different neurons in the hidden layer are tested to determine the optimal one. The correlation coefficient (R) is used to evaluate the linear correlation between the experimental and predicted results. However, it is unable to evaluate when the prediction is biased towards a local scope. Therefore, mean squared error (MSE) is also adopted to evaluate the ANN performance. The corresponding functions are as follows:
R = i = 1 n ( T i T ¯ ) ( Y i Y ¯ ) i = 1 n ( T i T ¯ ) 2 i = 1 n ( Y i Y ¯ ) 2
M S E = 1 n t m = 1 t i = 1 n [ T i ( m ) Y i ( m ) ] 2 × 100 %
where Ti is the experimental value, Yi is the predicted value and T ¯ and Y ¯ are the mean values of all the experimental and predicted results, respectively. n denotes the number of outputs and t denotes the amount of training sets. The convergence to MSE of 0.0167 is established in 1000 epochs.

2.3. Multivariate Regression Modeling

Multivariate regression analysis is a statistical tool that can be used to analyze the effects of parameters, develop the input–output relationship and obtain the optimal process parameter settings [17]. In this study, regression analysis is used to establish nonlinear regression models for the outputs, namely grain size (D), ultimate tensile strength (UTS), elongation (El.) and microhardness (HV).
The input variables are coded to reduce the difference in the scales of variables. Coding is performed for each input based on its minimum and maximum. The corresponding formula is as follows:
x i = 2 X i X i 0 Δ X i
where xi (i = 1, 2, 3, 4) is the coded input variable, Xi is the actual input variable, Xi0 is the value of Xi at the center level and ΔXi is the variation range in Xi. Regression models and statistical analysis are established based on the data in Appendix A. The fitted second-order polynomial regression model describing Y (output variables) expressed as a function of input variables is listed below:
Y ( x ) = b 0 + i = 1 4 b i x i + i = 1 4 j = i + 1 4 b i j x i x j + i = 1 4 b i i x i 2
where b0 is a constant, and bi, bii and bij (i, j = 1, 2, 3, 4) are the coefficients of linear, quadratic and cross product terms, respectively. Furthermore, the coefficients are calculated using the Least Squares method, which is a trial-and-error process. The statistical adequacy of regression models is determined by the coefficient of determination (R2) and Fisher’s criterion (F-test). The developed regression models can provide predictive models for integrated optimization calculation of heat treatment process of the ZE41-xCa-ySr alloy gearbox casting.

2.4. Integrated Computing Platform Building

The Integrated Computing Platform is constructed based on an integrated software platform SiPESC (Software Integration Platform for Engineering and Scientific Computation) [28]. SiPESC adopts the “microkernel and plugin” architecture, achieving the dynamic combination, extension and cooperative work of plugins. Therefore, the Integrated Computing Platform integrates multiple software by combining the scripting language with command line operation of simulation software, realizing the entire process simulation from CAD model input to “process-microstructure/defect-property” simulation. In addition, many algorithms are embedded in the Integrated Computing Platform, such as BFGS, SLP, SQP, MMA and GA algorithms, which are used to search for optimal parameter values.
In this work, the optimization calculation of composition and heat treatment processes for the ZE41-xCa-ySr alloy gearbox casting is implemented using the Integrated Computing Platform. Integrated Computing Platform batch calls MATLAB, achieving the optimization calculation of “heat treatment process/composition—microstructure/property” for the gearbox casting. However, before that, developed regression models need to be compiled into a “m.” script file for MATLAB. Figure 6 presents the architecture diagram showing “heat treatment process/composition-microstructure/property” integrated computing system. The optimization calculation of heat treatment and composition parameters for ZE41-xCa-ySr alloy gearbox casting involves the following steps:
  • Input the geometric model of ZE41 alloy gearbox casting.
  • Add design variables, such as Ca content (X1), Sr content (X2), aging temperature (X3) and aging time (X4). Moreover, set initial values and calculation ranges of separate variables, as shown in Table 3.
  • Add an external program for the Integrated Computing Platform to call. Here, it is MATLAB.
  • Set the target variable. Set the target variable to UTS in order to obtain the optimal design variables corresponding to the maximum UTS.
  • Select the optimization algorithm. The SQP (Sequential Quadratic Programming) method is selected as the optimization algorithm, which is an iterative method for nonlinear optimization.
  • Solve. The Integrated Computing Platform starts to make iterative optimization calculations and output each calculation result. Outputs include grain size (D), ultimate tensile strength (UTS), elongation (El.) and microhardness (HV).

3. Results and Discussion

3.1. Subsection

The aim of using the ANN model is to develop a prediction model which correlates between the input variables such as Ca, Sr content, aging temperature and aging time, and outputs such as D, UTS, El. and HV of the alloys. ANN Models are developed based on trial and error by adjusting the number of neurons in the hidden layer until maximum R and minimum MSE are obtained. The values of R and MSE for the different ANN models are shown in Figure 7, respectively. The results suggest that the ANN model with 12 neurons in the hidden layer obtains the maximum R value and the minimum MSE value, suggesting that 12 hidden-layer neurons of the ANN model exhibit favorable performance. As a result, the optimal ANN architecture is 4-12-4 (4 neurons in the input layer, 12 neurons in the hidden layer and 4 neurons in the output layer).
Figure 8 shows the performance of the optimal ANN model, indicating that the ANN model prediction exhibits good correlation with experimental results in training (R = 0.997), testing (R = 0.999), validation (R = 0.995) and whole (R = 0.996) datasets. The training convergence curve for the optimal ANN model is given in Figure 9. In the training stage of the ANN model, it is observed that the MSE decreases with the increasing number of iterations and reaches the error goal at about 120 epochs. Therefore, the well-trained ANN model has robust generalization ability.
Figure 10 presents the comparisons between the experimental and the predicted results for the entire datasets, suggesting that the ANN model has an adequate approximation to experimental values. Therefore, the microstructure and mechanical properties of the heat-treated ZE41-xCa-ySr alloy can be predicted with high accuracy using this ANN model.
Table 4 shows the predicted and experimental values of the BP neural network and percentage error. The results show that, for all the input parameters, the neural network projections agree very well with the experimental values. The percentage of error (P) between the experimental and ANN results is within 5%. These error levels are acceptable. As a consequence, the accuracy of the ANN model is validated. The percentage of error (P) is quantified as follows:
P e r c e n t a g e   o f   e r r o r   ( % ) = 100 ( T i Y i T i )
where Ti is the experimental value and Yi is the predicted value, respectively.

3.2. Model Development Using Regression Model

Equations (7)–(10) present the multivariate regression models and their adequacy validations. The coefficients of the terms suggest the influence of the variables on the response variables. The coefficient of determination (R2) is adopted to test the fit of regression models. For the model of D, the value of R2 is 0.747, which indicates that 74.7% of the variability in the response (D) can be explained by the regression model. In addition, the F-test is used to determine whether all the terms of regression models are statistically significant. If the calculated F ratio exceeds the critical F1−α, k−1, n−k value with degrees of freedom (k−1) and (n−k), the terms are significant at the α level of significance (α = 0.05, k is the number of terms, and n is the number of the experimental dataset). Consequently, from Equations (7)–(10), we can see that developed models quantify the effects of alloying elements and the aging process on the response variables, and the multivariate regression models are capable of making accurate predictions. The regression models are provided to the Integrated Computing Platform for “heat treatment process/composition-microstructure/property” integrated optimization calculation.
D = 38.466 + 0.599 x 1 3.251 x 2 + 1.418 x 3 0.448 x 4 + 0.621 x 1 · x 3 + 4.137 x 1 · x 4 0.866 x 2 · x 3 1.549 x 2 · x 4 + 6.239 x 3 2 0.617 x 4 2 R 2 = 0.747 , F   ratio = 2.398 > F 0.95 , 10 , 19 = 2.378
U S T = 193.052 + 3.231 x 1 + 3.22 x 2 2.209 x 3 + 11.777 x 4 + 2.929 x 1 · x 3 3.064 x 2 · x 3 1.577 x 2 · x 4 + 0.621 x 3 · x 4 11.143 x 3 2 22.021 x 4 2 R 2 = 0.8 , F   ratio = 3.383 > F 0.95 , 10 , 19 = 2.378
E l . = 3.584 0.069 x 1 0.269 x 2 + 0.205 x 3 0.579 x 4 + 0.233 x 2 · x 3 + 0.099 x 2 · x 4 0.204 x 3 · x 4 + 0.41 x 3 2 + 0.919 x 4 2 R 2 = 0.724 , F   ratio = 2.445 > F 0.95 , 9 , 20 = 2.393
H V = 40.869 + 32.815 x 1 31.78 x 2 1.012 x 3 + 2.401 x 4 + 32.553 x 1 · x 2 0.254 x 1 · x 3 0.748 x 1 · x 4 + 0.526 x 2 · x 3 0.189 x 2 · x 4 2.76 x 3 2 10.007 x 4 2 R 2 = 0.905 , F   ratio = 7.392 > F 0.95 , 11 , 18 = 2.377 .
where, x1 = 10X1 − 1, x2 = 10X2 − 1, x3 = 0.04X3 − 13, x4 = 0.0625X4 − 1.

3.3. Integrated Optimization Calculation

Based on the ICME paradigm, an Integrated Computing Platform for the “process-microstructure-property” simulation is established, which can realize the property prediction and process optimization. The Integrated Computing Platform first calls ProCAST software to implement “casting process-microstructure/defect” simulation, and then calls MATLAB software for “microstructure–property” calculation. Figure 11 shows the calculation interface of the Integrated Computing Platform. The Integrated Computing Platform automatically iterates 76 times. Table 5 lists partial iterations of integrated optimization calculation, suggesting that the optimized design variables correspond to the maximum UTS, HV and minimum D. Table 6 shows the optimal solution corresponding to the 61st iteration. Then, the optimized results are validated by the developed ANN model. The comparison of results is also shown in Table 5. Through analyzing the obtained results, it can be found that ANN predicted results coincide well with the integrated computing results. Therefore, the ANN model and multivariate regression models have high prediction accuracy. Furthermore, the optimization of the heat treatment process and composition based on property prediction is achieved through the Integrated Computing Platform.

4. Conclusions

  • The ANN model was developed using the BP algorithm. The optimal architecture (4-12-4) processed the maximum R value and the minimum MSE value. The ANN model was capable of predicting the microstructure and mechanical properties of heat-treated ZE41-xCa-ySr alloys with high reliability.
  • Multivariate regression analysis was employed to model the microstructure and mechanical properties of the heat-treated ZE41-xCa-ySr alloys. The adequacy of the models was tested by the coefficient of determination and Fisher’s criterion. All the nonlinear regression models were statistically adequate, which provided mathematical models for Integrated Computing Platform.
  • Based on SiPESC software, the Integrated Computing Platform was established by combining the scripting language with command line operation of the simulation software, realizing “process-microstructure/defect-property” simulation. An Integrated Computing Platform called MATLAB achieved the optimization calculation of “heat treatment process/composition—microstructure/property” for the ZE41-xCa-ySr alloy gearbox casting. The optimum aging temperature of the ZE41-0.17Ca-0.2Sr alloy is 322 °C, and the corresponding aging time is 11 h.

Author Contributions

Conceptualization, Y.F. and H.H.; methodology, C.L.; software, Y.S.; validation, Y.X., J.W. and X.Z.; formal analysis, Z.S.; writing—original draft preparation, Y.F.; writing—review and editing, H.H.; visualization, Y.X.; supervision, X.Z.; project administration, Y.F.; funding acquisition, Y.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Major Special Projects of the Plan “Science and Technology Innovation 2025” (No. 2019B10105, 2019B10086, 2019B10103, 2020Z096, 2020Z060, 2020Z048, 20221ZDYF020197), the Ningbo Natural Science Foundation (No. 202003N4340) and Zhejiang province military–civilian integration technology collaborative innovation platform.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The experimental results.
Table A1. The experimental results.
Coded Input VariablesActual Input VariablesActual Output Variables
x1x2x3x4X1: Ca (wt.%)X2: Sr (wt.%)X3: Ta (°C)X4: ta (h)Y1: D (μm)Y2: UTS (MPa)Y3: El. (%)Y4: HV
Center point (0)0.10.132516 ----
Range ΔXi0.20.25032 ----
High level (1)0.20.235032 ----
Low level (−1)003000 ----
No.
1−1−1−1−100300048.40130.206.2054.00 ± 1.22
2−1−1−1−0.7500300449.08169.305.3466.64 ± 2.13
3−1−1−1−0.375003001049.15163.904.4867.57 ± 1.54
4−1−1−11003003239.29170.705.1165.48 ± 2.01
5−1−10−0.687500325545.49170.505.7067.00 ± 1.87
6−1−10−0.375003251034.29189.003.8067.57 ± 1.94
7−1−101003253240.01168.703.2064.40 ± 2.34
8−1−11−0.62500350653.14153.304.8663.22 ± 2.01
9−1−11−0.500350847.83158.004.6364.25 ± 2.84
10−1−111003503246.93168.904.9661.72 ± 1.17
111−1−1−10.20300036.70131.605.6059.00 ± 2.18
121−1−1−0.750.20300449.59162.704.7767.10 ± 1.31
131−1−1−0.250.203001250.09158.704.6767.65 ± 1.94
141−1−110.203003247.64186.105.3064.39 ± 2.00
151−10−0.68750.20325541.86185.205.6466.90 ± 2.31
161−10−0.250.203251242.98194.403.4574.25 ± 1.94
171−1010.203253251.93170.403.1863.50 ± 1.76
181−11−0.6250.20350648.25171.204.7863.30 ± 2.74
191−11−0.3750.203501053.16175.205.1664.50 ± 2.03
201−1110.203503251.76175.564.5661.40 ± 1.79
2111−1−10.20.2300031.30144.104.9061.00 ± 2.11
2211−1−0.750.20.2300438.15187.704.1067.30 ± 1.84
2311−1−0.250.20.23001242.98173.303.4569.50 ± 2.31
2411−110.20.23003249.39179.303.0863.90 ± 1.71
25110−0.68750.20.2325541.19176.203.4966.90 ± 2.54
26110−0.250.20.23251246.01208.003.5077.10 ± 1.90
……
11711010.20.23253241.74173.805.1764.80 ± 1.57
118111−0.6250.20.2350640.15170.324.9165.10 ± 2.09
119111−0.3750.20.23501042.40167.404.7268.90 ± 2.27
12011110.20.23503245.14190.703.5463.20 ± 1.68
Table A2. Nomenclature.
Table A2. Nomenclature.
SymbolImplicationUnits
MSEMean squared error-
RThe correlation coefficient-
DGrain sizeμm
UTSUltimate tensile strengthMPa
El.Elongation%
HVMicrohardness-
taAging times
T aAging temperature°C

References

  1. Tresa, M.P. Weight loss with magnesium alloys. Science 2010, 328, 986–987. [Google Scholar]
  2. Kondori, B.; Benzerga, A.A. Modeling damage accumulation to fracture in a magnesium-rare earth alloy. Acta Mater. 2017, 124, 225–236. [Google Scholar] [CrossRef]
  3. Wang, Y.D.; Wu, G.H.; Liu, W.C.; Pang, S.; Zhang, Y.; Ding, W.J. Effects of chemical composition on the microstructure and mechanical properties of gravity cast Mg-xZn-yRE-Zr alloy. Mater. Sci. Eng. A 2014, 594, 52–61. [Google Scholar] [CrossRef]
  4. Yang, Y.; Zhang, K.; Li, X.G.; Li, Y.J.; Ma, M.L.; Shi, G.L.; Yuan, J.W. Microstructure and phase transformation of as-cast and annealed Mg-4Zn-1Y alloy containing quasi-crystal phase. Rare Met. 2015, 34, 239–244. [Google Scholar] [CrossRef]
  5. Hou, X.L.; Peng, Y.L.; Lv, P.; Cai, J.; Ji, L.; Guan, Q.F. Hot deformation behavior and microstructure evolution of a Mg-Gd-Nd-Y-Zn alloy. Rare Met. 2016, 35, 532–536. [Google Scholar] [CrossRef]
  6. Panchal, J.H.; Kalidindi, S.R.; McDowell, D.L. Key computational modeling issues in integrated computational materials engineering. Comput. Aided Des. 2013, 45, 4–25. [Google Scholar] [CrossRef]
  7. Qiu, X.; Yang, Q.; Cao, Z.Y.; Liu, Y.B.; Meng, J. Microstructure and mechanical properties of Mg-Zn-(Nd)-Zr alloys with different extrusion processes. Rare Met. 2016, 35, 841–849. [Google Scholar] [CrossRef]
  8. Luo, K.; Zhang, L.; Wu, G.H.; Liu, W.C.; Ding, W.J. Effect of Y and Gd content on the microstructure and mechanical properties of Mg-Y-RE alloys. J. Magnes. Alloy 2019, 7, 345–354. [Google Scholar] [CrossRef]
  9. Zhou, Y.Y.; Fu, P.H.; Peng, L.M.; Wang, D.; Wang, Y.X.; Hu, B.; Liu, M.; Anil, K.S.; Ding, W.J. Precipitation modification in cast Mg-1Nd-1Ce-Zr alloy by Zn addition. J. Magnes. Alloy 2019, 7, 113–123. [Google Scholar] [CrossRef]
  10. Wang, Y.D.; Wu, G.H.; Liu, W.C.; Pang, S.; Zhang, Y.; Ding, W.J. Influence of heat treatment on microstructure and mechanical properties gravity cast Mg-4.2Zn-1.5RE-0.7Zr magnesium alloy. Trans. Nonferrous Met. Soc. China 2013, 23, 3611–3620. [Google Scholar] [CrossRef]
  11. Ashtiani, H.R.R.; Shayanpoor, A.A. Processing Workability and Artificial Neural Network of AA1070 to the Prediction of Hot Flow Stress. Phys. Met. Metallogr. 2021, 122, 1426–1435. [Google Scholar] [CrossRef]
  12. Sharath, B.N.; Venkatesh, C.V.; Afzal, A.; Aslfattahi, N.; Aabid, A.; Baig, M.; Saleh, B. Multi ceramic particles inclusion in the aluminium matrix and wear characterization through experimental and response surface-artificial neural networks. Materials 2021, 14, 2895. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, D.W.; Ren, D.X.; Zhao, K.M.; Pan, S.; Guo, X.L. Effect of welding parameters on tensile strength of ultrasonic spot welded joints of aluminum to steel-By experimentation and artificial neural network. J. Manuf. Process. 2017, 30, 63–74. [Google Scholar] [CrossRef]
  14. Castin, N.; Pascuet, M.I.; Messina, L.; Domain, C.; Olsson, P.; Pasianot, R.C.; Malerba, L. Advanced atomistic models for radiation damage in Fe-based alloys: Contributions and future perspectives from artificial neural networks. Comput. Mater. Sci. 2018, 148, 116–130. [Google Scholar] [CrossRef]
  15. Yang, X.W.; Zhu, J.C.; Nong, Z.S.; He, D.; Lai, Z.H.; Liu, Y.; Liu, F.W. Prediction of mechanical properties of A357 alloy using artificial neural network. Trans. Nonferrous Met. Soc. China 2013, 23, 788–795. [Google Scholar] [CrossRef]
  16. Conduit, B.D.; Jones, N.G.; Stone, H.J.; Conduit, G.J. Design of a nickel-base superalloy using a neural network. Mater. Des. 2017, 131, 358–365. [Google Scholar] [CrossRef]
  17. Patel, M.P.G.; Mathew, R.; Krishna, P.; Parappagoudar, M.B. Investigation of Squeeze Cast Process Parameters Effects on Secondary Dendrite Arm Spacing Using Statistical Regression and Artificial Neural Network Models. Procedia Technol. 2014, 14, 149–156. [Google Scholar] [CrossRef]
  18. Bhatt, A.; Parappagoudar, M.B. Modeling and Analysis of Mechanical Properties in Structural Steel-DOE Approach. Arch. Foundry Eng. 2015, 15, 5–12. [Google Scholar] [CrossRef]
  19. Chen, L.W.; Zhao, Y.H.; Wen, Z.Q.; Tian, J.Z.; Hou, H. Modelling and Optimization for Heat Treatment of Al-Si-Mg Alloy Prepared by Indirect Squeeze Casting Based on Response Surface Methodology. Mater. Res. 2017, 20, 1274–1281. [Google Scholar] [CrossRef]
  20. Akhtar, M.N.; Sathish, T.; Mohanavel, V.; Afzal, A.; Arul, K.; Ravichandran, M.; Rahim, I.A.; Alhady, S.S.N.; Bakar, E.A.; Saleh, B. Optimization of Process Parameters in CNC Turning of Aluminum 7075 Alloy Using L27 Array-Based Taguchi Method. Materials 2021, 14, 4470. [Google Scholar] [CrossRef]
  21. Zhao, M.C.; Liu, M.; Song, G.L.; Atrens, A. Influence of Microstructure on Corrosion of As-cast ZE41. Adv. Eng. Mater. 2008, 10, 104–111. [Google Scholar] [CrossRef]
  22. Neil, W.C.; Forsyth, M.; Howlett, P.C.; Hutchinson, C.R.; Hinton, B.R.W. Corrosion of magnesium alloy ZE41—The role of microstructural features. Corros. Sci. 2009, 51, 387–394. [Google Scholar] [CrossRef]
  23. Fu, Y.; Wang, H.; Liu, X.T.; Hai, H. Effect of calcium addition on microstructure, casting fluidity and mechanical properties of Mg-Zn-Ce-Zr magnesium alloy. J. Rare Earths. 2017, 35, 503–509. [Google Scholar] [CrossRef]
  24. Fu, Y.; Wang, H.; Zhang, C.; Hai, H. Effects of minor Sr additions on the as-cast microstructure, fluidity and mechanical properties of Mg-4.2Zn-1.7RE-0.8Zr-0.2Ca (wt%) alloy. Mater. Sci. Eng. A 2018, 723, 118–125. [Google Scholar] [CrossRef]
  25. Fu, Y.; Liu, C.; Hai, H.; Xu, Y.D.; Zhu, X.R. Effect of ageing treatment on microstructures, mechanical properties and corrosion behavior of Mg-Zn-RE-Zr alloy micro-alloyed with Ca and Sr. China Foundry 2021, 18, 131–140. [Google Scholar] [CrossRef]
  26. Ozerdem, M.S.; Kolukisa, S. Artificial neural network approach to predict the mechanical properties of Cu-Sn-Pb-Zn-Ni cast alloys. Mater. Des. 2009, 30, 764–769. [Google Scholar] [CrossRef]
  27. Britto, A.S.F.; Raj, R.E.; Mabel, M.C. Prediction and optimization of mechanical strength of diffusion bonds using integrated ANN-GA approach with process variables and metallographic characteristics. J. Manuf. Process. 2018, 32, 828–838. [Google Scholar] [CrossRef]
  28. Soundararajan, R.; Ramesh, A.; Sivasankaran, S.; Sathishkumar, A. Modeling and Analysis of Mechanical Properties of Aluminium Alloy (A413) Processed through Squeeze Casting Route Using Artificial Neural Network Model and Statistical Technique. Adv. Mater. Sci. Eng. 2015, 2015, 714762. [Google Scholar] [CrossRef]
  29. Li, J.C.; Zhao, D.L.; Ge, B.F.; Yang, K.W.; Chen, Y.W. A link prediction method for heterogeneous networks based on BP neural network. Physica A 2018, 495, 1–17. [Google Scholar] [CrossRef]
Figure 1. The flow chart of this work.
Figure 1. The flow chart of this work.
Crystals 12 01237 g001
Figure 2. Optical micrographs of the as-cast ZE41-Ca-Sr alloys: (a) ZE41; (b) ZE41-0.2Ca; (c) ZE41-0.2Ca-0.1Sr; (d) ZE41-0.2Ca-0.2Sr; (e) ZE41-0.2Ca-0.4Sr [24].
Figure 2. Optical micrographs of the as-cast ZE41-Ca-Sr alloys: (a) ZE41; (b) ZE41-0.2Ca; (c) ZE41-0.2Ca-0.1Sr; (d) ZE41-0.2Ca-0.2Sr; (e) ZE41-0.2Ca-0.4Sr [24].
Crystals 12 01237 g002
Figure 3. Optical micrographs of the peak-aged alloys: (a) ZE41, (b) ZE41-0.2Ca, (c) ZE41-0.2Ca-0.1Sr, (d) ZE41-0.2Ca-0.2Sr, (e) ZE41-0.2Ca-0.4Sr and (f) partial enlarged drawing of the ZE41-0.2Ca-0.2Sr alloy [25].
Figure 3. Optical micrographs of the peak-aged alloys: (a) ZE41, (b) ZE41-0.2Ca, (c) ZE41-0.2Ca-0.1Sr, (d) ZE41-0.2Ca-0.2Sr, (e) ZE41-0.2Ca-0.4Sr and (f) partial enlarged drawing of the ZE41-0.2Ca-0.2Sr alloy [25].
Crystals 12 01237 g003
Figure 4. The network configuration of the ANN model.
Figure 4. The network configuration of the ANN model.
Crystals 12 01237 g004
Figure 5. The mathematical model of the ANN.
Figure 5. The mathematical model of the ANN.
Crystals 12 01237 g005
Figure 6. The architecture diagram showing “heat treatment process/composition-microstructure/property” integrated computing system.
Figure 6. The architecture diagram showing “heat treatment process/composition-microstructure/property” integrated computing system.
Crystals 12 01237 g006
Figure 7. The values of R and MSE for the different ANN models.
Figure 7. The values of R and MSE for the different ANN models.
Crystals 12 01237 g007
Figure 8. Analysis of the optimal ANN model predicted results in comparison with experimental results during training, validation and testing of the model. (a) training (R = 0.997), (b) testing (R = 0.999), (c) validation (R = 0.995) and (d) whole (R = 0.996) datasets.
Figure 8. Analysis of the optimal ANN model predicted results in comparison with experimental results during training, validation and testing of the model. (a) training (R = 0.997), (b) testing (R = 0.999), (c) validation (R = 0.995) and (d) whole (R = 0.996) datasets.
Crystals 12 01237 g008
Figure 9. Training convergence curve for the optimal ANN model.
Figure 9. Training convergence curve for the optimal ANN model.
Crystals 12 01237 g009
Figure 10. The comparisons between the experimental and the predicted results for the entire datasets. (a) D, (b) UTS, (c) El. and (d) HV.
Figure 10. The comparisons between the experimental and the predicted results for the entire datasets. (a) D, (b) UTS, (c) El. and (d) HV.
Crystals 12 01237 g010
Figure 11. The calculation interface of Integrated Computing Platform.
Figure 11. The calculation interface of Integrated Computing Platform.
Crystals 12 01237 g011
Table 1. Chemical compositions of the experimental alloys (wt.%).
Table 1. Chemical compositions of the experimental alloys (wt.%).
Nominal AlloysActual Composition
MgZnREZrCaSr
Mg-4.2Zn-1.7RE-0.8ZrBal.4.091.670.70--
Mg-4.2Zn-1.7RE-0.8Zr-0.2CaBal.4.141.610.760.18-
Mg-4.2Zn-1.7RE-0.8Zr-0.2Ca-0.2SrBal.4.131.720.720.220.21
Table 2. The ANN model training parameters.
Table 2. The ANN model training parameters.
ParametersANN Model
Number of layers3
The number of neurons on the layersInput: 4, Hidden: 4~12, Output: 4
Transfer functionsHidden layer: Tan-Sigmoid
Output layer: Purelin
Training methodLevenberg–Marquardt (LM)
Initial weights and biasesRandomly between −1 and 1
Target error value0.0167
Learning rateVariable learning rate
Table 3. Initial values and calculation ranges of separate variables.
Table 3. Initial values and calculation ranges of separate variables.
ValuesCa Content (X1)/wt.%Sr Content (X2)/wt.%Aging Temperature (X3)/°CAging Time (X4)/h
Initial values0.10.132516
Calculation ranges0–0.20–0.2300–3500–32
Table 4. Predicted and experimental values of the BP neural network and percentage error.
Table 4. Predicted and experimental values of the BP neural network and percentage error.
Group 1Group 2Group 3
InputsCa Content
(wt.%)
Sr
Content
(wt.%)
Ta (°C)ta (h)Ca Content
(wt.%)
Sr
Content
(wt.%)
Ta (°C)ta (h)Ca Content
(wt.%)
Sr
Content
(wt.%)
Ta (°C)ta (h)
Values0.1750.232580.20.2325120.40.432530
OutputsD/
μm
UTS
/MPa
El.
/%
HVD/
μm
UTS
/MPa
El.
/%
HVD/
μm
UTS
/MPa
El.
/%
HV
ANN model predicted values37.12183.33.9371.926.19202.53.4574.950.21163.23.1160.2
Experimental values39.03185.43.7770.025.45208.03.5077.148.72159.83.2061.2
Percentage of error between the experimental and ANN resultsPD
(%)
PUTS
(%)
PEl.
(%)
PHV
(%)
PD
(%)
PUTS
(%)
PEl.
(%)
PHV
(%)
PD
(%)
PUTS
(%)
PEl.
(%)
PHV
(%)
4.891.13−4.24−2.71−2.912.641.432.85−3.06−2.132.811.63
Table 5. Partial iterations of integrated optimization calculation.
Table 5. Partial iterations of integrated optimization calculation.
NoCa/wt%Sr/wt%Ta/°Cta/hD/
μm
UTS
/MPa
El.
/%
HV
10.10.11153305.1838.362193.9993.51636.728
20.10010.11153305.8238.364194.0023.51736.764
30.10.111633020.838.357194.0013.51636.693
40.10.111533020.82138.360193.9973.51736.723
50.10.1115330.0320.838.366193.9903.51736.725
60.100040.1114329.9720.79538.357194.0093.51636.738
70.100120.1116329.92120.78738.349194.0283.51536.754
570.178520.2321.80129.70536.022193.8173.16859.771
580.199920.2321.37329.02837.104193.2383.43459.014
590.174480.2321.88210.32235.879200.0193.16167.539
600.174390.2321.88410.23035.876200.0203.16067.488
610.174570.2321.88411.27135.879200.0253.16067.921
620.174390.198321.88410.29135.883200.0133.16167.503
630.174390.2321.88410.11335.877200.0173.16167.483
640.174390.2321.90512.20135.876200.0193.16167.488
650.174390.2318.18912.26136.092199.7983.12867.426
660.174390.2320.03612.29135.950199.9693.14267.472
670.174390.2321.10310.36135.899200.0133.15267.485
680.174390.2321.71610.73635.879200.0203.15867.487
690.174570.2321.71612.26735.883200.0253.15867.811
700.174390.199321.71612.27135.887200.0143.15967.503
710.174390.2321.71610.31135.881200.0183.15867.483
720.174390.2321.73810.29135.879200.0203.15967.488
730.173450.2321.71610.67135.865199.9933.15966.875
740.173920.2321.71610.31135.872200.0073.15967.181
750.174200.2321.71611.36135.877200.0143.15967.358
760.174390.2321.71610.36035.879200.0203.15867.487
Table 6. The optimal solution and results of model validation.
Table 6. The optimal solution and results of model validation.
Input Variables of the Optimal Solution (Maximum UTS, HV and Minimum D)Output Variables through Integrated Computing Platform
Ca/wt.%Sr/wt.%Ta/°Cta/hD/μmUTS/MPaEl./%HV
0.170.23221135.88200.033.1667.92
ANN model predicted results
36.41194.743.1368.79
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fu, Y.; Liu, C.; Song, Y.; Hao, H.; Xu, Y.; Shao, Z.; Wang, J.; Zhu, X. Modeling of Microstructure and Mechanical Properties of Heat Treated ZE41-Ca-Sr Alloys for Integrated Computing Platform. Crystals 2022, 12, 1237. https://doi.org/10.3390/cryst12091237

AMA Style

Fu Y, Liu C, Song Y, Hao H, Xu Y, Shao Z, Wang J, Zhu X. Modeling of Microstructure and Mechanical Properties of Heat Treated ZE41-Ca-Sr Alloys for Integrated Computing Platform. Crystals. 2022; 12(9):1237. https://doi.org/10.3390/cryst12091237

Chicago/Turabian Style

Fu, Yu, Chen Liu, Yunkun Song, Hai Hao, Yongdong Xu, Zhiwen Shao, Jun Wang, and Xiurong Zhu. 2022. "Modeling of Microstructure and Mechanical Properties of Heat Treated ZE41-Ca-Sr Alloys for Integrated Computing Platform" Crystals 12, no. 9: 1237. https://doi.org/10.3390/cryst12091237

APA Style

Fu, Y., Liu, C., Song, Y., Hao, H., Xu, Y., Shao, Z., Wang, J., & Zhu, X. (2022). Modeling of Microstructure and Mechanical Properties of Heat Treated ZE41-Ca-Sr Alloys for Integrated Computing Platform. Crystals, 12(9), 1237. https://doi.org/10.3390/cryst12091237

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop