The Effect of Hot Forming–Quenching and Heat Treatment Processes on the Mechanical Properties of AA6016 Aluminum Alloy Sheets

: This study explored the impact of Hot Forming–Quenching (HFQ) and heat treatment processes on the mechanical properties of AA6016 sheets. The experimental findings demonstrated that at high-temperature pre-straining (HT-PS) of 15%, the strength performance of the AA6016 sheet exhibited enhancement, with a progressive increase in both the heat treatment temperature and duration. Conversely, under HT-PS conditions of 3% and 7%, the heat treatment process exhibited a relatively modest impact on the mechanical properties of the AA6016 sheet. Differential scanning calorimetry (DSC) was employed to understand the influence of different process conditions on the precipitated phases. By comparing the precipitation peaks of the β ′′ phase at HT-PS of 3% and 15%, it was observed that the precipitation peak of the β ′′ phase decreased with an increase in HT-PS. This indicated that HT-PS promoted the precipitation of the β ′′ phase. In order to forecast the mechanical performance of the AA6016 sheets after applying various pre-straining and heat treatment parameters, two models were used: a backpropagation (BP) neural network and a genetic algorithm (GA)-BP neural network. These models were evaluated for their fitting and predictive capabilities. The research findings demonstrated that the GA-BP neural network model exhibited superior fitting and predictive accuracy compared to the BP neural network model.


Introduction
Under high-temperature conditions, aluminum alloy sheets exhibit a remarkable plastic deformation ability, which greatly facilitates the production of high-performance components.Among various aluminum alloy hot stamping techniques, Hot Forming-Quenching (HFQ) technology has gained significant attention.This technology involves rapidly transferring a fully solutionized aluminum alloy onto a cold die, promptly closing the die, and achieving in-die quenching [1].In a study conducted by Omer [2] et al., the influence of the HFQ process on the mechanical properties of AA5754 was investigated.The results demonstrated that the HFQ process was highly advantageous for forming parts and enhancing their performance.Similarly, Zheng [3] et al. explored the effects of the HFQ process on AA7075 and AA6082.The findings indicated that HFQ significantly reduced the post-treatment mechanical properties of highly quench-sensitive AA7075, but improved those of less quench-sensitive AA6082.This process enables dual enhancement of the formability and strength of aluminum alloys, aligning with the lightweighting trend in the automotive industry.As a result, it holds promising prospects for development in the automotive sector.aluminum alloys that take into account the HFQ and heat treatment processes.Therefore, the objective of this study is to examine the microstructure of AA6016 according to HT-PS parameters and heat treatment processes parameters.The experimental results serve as the basis for employing the backpropagation (BP) neural network model and a BP neural network model optimized with a genetic algorithm (GA-BP) to fit the data.Subsequently, the accuracy of fitting and the prediction of the two models are analyzed.

Materials
The experimental specimens were fabricated from 1.1 mm thick AA6016 rolled aluminum sheets, with their chemical composition provided by the material manufacturer, as shown in Table 1.The dimensions of the hot tensile specimens were depicted in Figure 1, featuring a total length of 110 mm and a gauge length of 15mm.Most of the studies mentioned above have primarily focused on investigating the mechanical properties of aluminum alloy sheets after undergoing cold deformation and heat treatment processes.However, limited research has been conducted on the effects of hot deformation, particularly the HFQ process and heat treatment process, on the properties of aluminum alloy sheets.Furthermore, there is a lack of research on constitutive models for aluminum alloys that take into account the HFQ and heat treatment processes.Therefore, the objective of this study is to examine the microstructure of AA6016 according to HT-PS parameters and heat treatment processes parameters.The experimental results serve as the basis for employing the backpropagation (BP) neural network model and a BP neural network model optimized with a genetic algorithm (GA-BP) to fit the data.Subsequently, the accuracy of fitting and the prediction of the two models are analyzed.

Materials
The experimental specimens were fabricated from 1.1 mm thick AA6016 rolled aluminum sheets, with their chemical composition provided by the material manufacturer, as shown in Table 1.The dimensions of the hot tensile specimens were depicted in Figure 1, featuring a total length of 110 mm and a gauge length of 15mm.

HT-PS and Heat Treatment Experiments
In order to accurately reflect the mechanical behavior of AA6016 during the hot stamping process and after heat treatment, an experimental plan for AA6016 was developed based on the deformation temperature during hot stamping and the heat treatment temperature and the holding time during the heat treatment.The specific process parameters for AA6016 are outlined in Table 2, and the designated process flow is illustrated in Figure 2. It should be noted that the heat treatment in this paper is intended to simulate the automotive paint baking process.

HT-PS and Heat Treatment Experiments
In order to accurately reflect the mechanical behavior of AA6016 during the hot stamping process and after heat treatment, an experimental plan for AA6016 was developed based on the deformation temperature during hot stamping and the heat treatment temperature and the holding time during the heat treatment.The specific process parameters for AA6016 are outlined in Table 2, and the designated process flow is illustrated in Figure 2. It should be noted that the heat treatment in this paper is intended to simulate the automotive paint baking process.
The specimens were subjected to HT-PS using the 200 kN Thermecmastor (Tokyo, Japan) thermal simulation testing machine.Thermocouples (Type K) were welded to the center of the gauge length section of the samples to measure the temperature.The welding of the thermocouples onto the gauge length section allowed for control over the heating time and cooling rate.The AA6016 sheet undergoes a solution treatment at 550 • C for 200 s, followed by air cooling to 450 • C (cooling rate of 10 • C/s).Subsequently, at 390 • C, the sheets were deformed at a strain rate of 0.01 s −1 to different pre-strain levels (3%, 7%, 15%, 18%).Simultaneously, the cooling process was conducted during the deformation of the sheets.After the pre-strain, all the samples were kept at room temperature for 2 days and then subjected to the heat treatment process following natural aging.It should be noted that each set of processes was tested at least twice.If there were significant differences between the results of the two tests, a third test was conducted to ensure the accuracy of the experimental data.A representative set of data was selected from each group of experiments to construct the stress-strain curves.The mechanical properties of the AA6016 sheet were obtained by calculating the deviation from the two most accurate sets of data chosen.The specimens were subjected to HT-PS using the 200 kN Thermecmastor (Tokyo, Japan) thermal simulation testing machine.Thermocouples (Type K) were welded to the center of the gauge length section of the samples to measure the temperature.The welding of the thermocouples onto the gauge length section allowed for control over the heating time and cooling rate.The AA6016 sheet undergoes a solution treatment at 550 °C for 200 s, followed by air cooling to 450 °C (cooling rate of 10 °C/s).Subsequently, at 390 °C, the sheets were deformed at a strain rate of 0.01 s −1 to different pre-strain levels (3%, 7%, 15%, 18%).Simultaneously, the cooling process was conducted during the deformation of the sheets.After the pre-strain, all the samples were kept at room temperature for 2 days and then subjected to the heat treatment process following natural aging.It should be noted that each set of processes was tested at least twice.If there were significant differences between the results of the two tests, a third test was conducted to ensure the accuracy of the experimental data.A representative set of data was selected from each group of experiments to construct the stress-strain curves.The mechanical properties of the AA6016 sheet were obtained by calculating the deviation from the two most accurate sets of data chosen.

Mechanical Testing before and after the Heat Treatment Process
Both the samples before and after the heat treatment process underwent quasi-static tensile testing at room temperature.This testing was conducted using a universal testing machine, specifically model 5982 (Boston, MA, USA).The crosshead speed during the testing was set at 0.8 mm/min.This testing procedure provided valuable data for evaluating the mechanical properties of the samples.

Microstructure Characterization
Samples were selected from the vicinity of the gauge length section of the tensile specimens.These samples underwent grinding and polishing processes until a mirror-like finish was achieved, without the need for etching.Subsequently, scanning electron microscopy (SEM) with a tungsten filament Vega scanning electron microscope (Zeiss, accelerating voltage of 10 kV) was utilized to perform microstructure analysis on the samples.SEM combined with backscattered electrons (BSE) was used to enhance the visibility and facilitate the statistical analysis of the second phase (iron-rich phase).All the specimens were scanned using SEM-BSE, capturing a total of 24 images each.Images numbered 8-16 represent the second phase (iron-rich phase) in the central region, while images numbered 1-7 and 17-24 represent the second phase (iron-rich phase) on both sides.Data statistics were calculated for all the captured images, employing threshold segmentation and batch processing with machine learning techniques.Additionally, the fractured surfaces of the tensile samples were cut and analyzed using the same SEM (accelerating voltage of 15 kV).

Mechanical Testing before and after the Heat Treatment Process
Both the samples before and after the heat treatment process underwent quasi-static tensile testing at room temperature.This testing was conducted using a universal testing machine, specifically model 5982 (Boston, MA, USA).The crosshead speed during the testing was set at 0.8 mm/min.This testing procedure provided valuable data for evaluating the mechanical properties of the samples.

Microstructure Characterization
Samples were selected from the vicinity of the gauge length section of the tensile specimens.These samples underwent grinding and polishing processes until a mirrorlike finish was achieved, without the need for etching.Subsequently, scanning electron microscopy (SEM) with a tungsten filament Vega scanning electron microscope (Zeiss, accelerating voltage of 10 kV) was utilized to perform microstructure analysis on the samples.SEM combined with backscattered electrons (BSE) was used to enhance the visibility and facilitate the statistical analysis of the second phase (iron-rich phase).All the specimens were scanned using SEM-BSE, capturing a total of 24 images each.Images numbered 8-16 represent the second phase (iron-rich phase) in the central region, while images numbered 1-7 and 17-24 represent the second phase (iron-rich phase) on both sides.Data statistics were calculated for all the captured images, employing threshold segmentation and batch processing with machine learning techniques.Additionally, the fractured surfaces of the tensile samples were cut and analyzed using the same SEM (accelerating voltage of 15 kV).
Samples were taken from the gauge length section of the specimens that underwent HT-PS.The samples were directly fabricated into circular discs with a diameter of 3 mm and a thickness of 1 mm.The precipitation behavior of AA6016 was analyzed using the Netzsch DSC 214 (Selb, Germany) differential scanning calorimetry (DSC) instrument.This method involved operating the 1100LF system within a temperature range of 50 • C to 550 • C, with a heating rate of 10 • C/min, under an argon atmosphere.

BP Neural Network Model
The essence of an artificial neural network is to recognize and predict data using nonlinear equations.After inputting data into the input layer, the initial data begin to pass through the neural network.The data inputted by multiple neurons undergo multiple iterations in the hidden layer and are eventually output in the output layer.The backpropagation (BP) neural network [20] has historically been one of the most widely used error backpropagation network structures.This algorithm had advantages such as a high classification accuracy, a fast self-learning speed, and strong parallel processing capabilities.BP primarily achieves the predetermined design or minimum value of the system error of the entire network by iteratively adjusting the network weights and thresholds [21].Its structure is shown in Figure 3.In this study, HT-PS, heat treatment temperature, heat treatment time, and strain were used as inputs, while stress, yield strength, tensile strength, and elongation were used as outputs.This allowed for fitting of the constitutive model of AA6016 by utilizing HT-PS and heat treatment parameters as the parameters.
Metals 2024, 14, x FOR PEER REVIEW 5 of 19 Samples were taken from the gauge length section of the specimens that underwent HT-PS.The samples were directly fabricated into circular discs with a diameter of 3 mm and a thickness of 1 mm.The precipitation behavior of AA6016 was analyzed using the Netzsch DSC 214 (Selb, Germany) differential scanning calorimetry (DSC) instrument.This method involved operating the 1100LF system within a temperature range of 50 °C to 550 °C, with a heating rate of 10 °C/min, under an argon atmosphere.

BP Neural Network Model
The essence of an artificial neural network is to recognize and predict data using nonlinear equations.After inputting data into the input layer, the initial data begin to pass through the neural network.The data inputted by multiple neurons undergo multiple iterations in the hidden layer and are eventually output in the output layer.The backpropagation (BP) neural network [20] has historically been one of the most widely used error backpropagation network structures.This algorithm had advantages such as a high classification accuracy, a fast self-learning speed, and strong parallel processing capabilities.BP primarily achieves the predetermined design or minimum value of the system error of the entire network by iteratively adjusting the network weights and thresholds [21].Its structure is shown in Figure 3.In this study, HT-PS, heat treatment temperature, heat treatment time, and strain were used as inputs, while stress, yield strength, tensile strength, and elongation were used as outputs.This allowed for fitting of the constitutive model of AA6016 by utilizing HT-PS and heat treatment parameters as the parameters.

Data Normalization
A large amount of experimental data were obtained from uniaxial tensile tests, and the range of the data was relatively large.Therefore, the data needed to be standardized and organized into a certain range for ease of calculation.In this paper, the method of maximum-minimum normalization was adopted for normalization, and the method for normalization of type A is as follows: In this method, Y represents the normalized data, mina is the minimum value in class a, and maxa is the maximum value.This method was used to normalize the temperature and stress data obtained from the uniaxial tensile tests.

Function Selection and Forward Propagation
The Sigmoid function was chosen as the transfer function in the model, and its mathematical formula is as follows:

Data Normalization
A large amount of experimental data were obtained from uniaxial tensile tests, and the range of the data was relatively large.Therefore, the data needed to be standardized and organized into a certain range for ease of calculation.In this paper, the method of maximum-minimum normalization was adopted for normalization, and the method for normalization of type A is as follows: In this method, Y represents the normalized data, mina is the minimum value in class a, and maxa is the maximum value.This method was used to normalize the temperature and stress data obtained from the uniaxial tensile tests.

Function Selection and Forward Propagation
The Sigmoid function was chosen as the transfer function in the model, and its mathematical formula is as follows: Take x 1 , x 2 , x 3 ∼ x n as the input variables, y as the output variable, and u m as the output of the hidden layer neurons, which, passing through the neurons, goes through the activation function; let v ij be the weight of the ith input variable and the jth hidden layer Metals 2024, 14, 599 6 of 18 neuron, and let θ u j be the jth neuron threshold of the hidden layer u: then, the expression of u j can be written as: Let w j be the weight of the jth neuron connected to y and θ y be the threshold of y, which is obtained by activating the Sigmoid function:

Genetic Algorithm Optimization of BP Neural Networks
Although BP neural networks have strong nonlinear analysis capabilities, they often encounter issues such as local error minimization and slow convergence, which can result in reduced predictive accuracy.Genetic algorithms (GAs) are adaptive parallel optimization algorithms that simulate biological genetic evolution mechanisms.They possess several advantages, including global search, high parallelism, and a strong generalization ability [22].By applying genetic algorithms to BP neural networks, we can enhance the accuracy of the network training convergence and leverage the global optimization properties of genetic algorithms.In the GA-BP neural network, the initial weights and thresholds of the BP neural network are improved using a genetic algorithm to obtain the global optimal region.Subsequently, the BP neural network is trained using these optimal weights and thresholds until the training process is completed while simultaneously searching for parameters corresponding to the optimal fitness [23].Figure 4 depicts a schematic diagram generated by the machine learning model.In this study, the genetic algorithm primarily serves to optimize the weights and thresholds of the neural network.
Take  ,  ,  ~ as the input variables, y as the output variable, and  as the output of the hidden layer neurons, which, passing through the neurons, goes through the activation function; let  be the weight of the ith input variable and the jth hidden layer neuron, and let  be the jth neuron threshold of the hidden layer u: then, the expression of  can be written as: Let  be the weight of the jth neuron connected to y and  be the threshold of y, which is obtained by activating the Sigmoid function:

Genetic Algorithm Optimization of BP Neural Networks
Although BP neural networks have strong nonlinear analysis capabilities, they often encounter issues such as local error minimization and slow convergence, which can result in reduced predictive accuracy.Genetic algorithms (GAs) are adaptive parallel optimization algorithms that simulate biological genetic evolution mechanisms.They possess several advantages, including global search, high parallelism, and a strong generalization ability [22].By applying genetic algorithms to BP neural networks, we can enhance the accuracy of the network training convergence and leverage the global optimization properties of genetic algorithms.In the GA-BP neural network, the initial weights and thresholds of the BP neural network are improved using a genetic algorithm to obtain the global optimal region.Subsequently, the BP neural network is trained using these optimal weights and thresholds until the training process is completed while simultaneously searching for parameters corresponding to the optimal fitness [23].Figure 4 depicts a schematic diagram generated by the machine learning model.In this study, the genetic algorithm primarily serves to optimize the weights and thresholds of the neural network.The multi-objective problem can be described by the following equation: In the equation, T min represents vector minimization, indicating that each sub-objective function f 1 (x), f 2 (x), f 3 (x), • • • f n (x) should be minimized to the greatest extent possible.In practical applications, there is often interdependence and conflict between the objectives and parameters.It becomes challenging to accurately describe the relationship between the target values and input values using precise functional relationships.Moreover, conflicts Metals 2024, 14, 599 7 of 18 frequently arise between multiple target values.Enhancing the performance of one subobjective may lead to a decrease in the performance of another sub-objective [24].In other words, it is not feasible to simultaneously achieve optimal values for multiple sub-objectives.Therefore, compromises and coordination must be made to optimize each sub-objective as much as possible, aiming to achieve the best overall optimization results for the system.(10,20,30) denotes the heat treatment durations.In Figure 5a, for HT-PS (HT-PS) of 3%, different heat treatments had a minimal effect on the strength of the AA6016 sheet but slightly impacted the elongation.Figure 5b showed that for HT-PS of 7%, particularly with a heat treatment temperature of 185 • C for 30 min, the strength of the AA6016 sheet increased.In Figure 5c, for an HT-PS of 15%, there was a noticeable increase in strength with an increase in the heat treatment temperature and time, but the elongation decreased significantly.

Mechanical Properties
In the equation,  represents vector minimization, indicating that each sub-objective function  (),  (),  (), ⋯  () should be minimized to the greatest extent possible.In practical applications, there is often interdependence and conflict between the objectives and parameters.It becomes challenging to accurately describe the relationship between the target values and input values using precise functional relationships.Moreover, conflicts frequently arise between multiple target values.Enhancing the performance of one sub-objective may lead to a decrease in the performance of another subobjective [24].In other words, it is not feasible to simultaneously achieve optimal values for multiple sub-objectives.Therefore, compromises and coordination must be made to optimize each sub-objective as much as possible, aiming to achieve the best overall optimization results for the system.(10,20,30) denotes the heat treatment durations.In Figure 5a, for HT-PS (HT-PS) of 3%, different heat treatments had a minimal effect on the strength of the AA6016 sheet but slightly impacted the elongation.Figure 5b showed that for HT-PS of 7%, particularly with a heat treatment temperature of 185 °C for 30 min, the strength of the AA6016 sheet increased.In Figure 5c, for an HT-PS of 15%, there was a noticeable increase in strength with an increase in the heat treatment temperature and time, but the elongation decreased significantly.The mechanical properties under different experimental conditions are shown in Table 3, and Figure 6 shows a relationship graph of the yield strength (YS), tensile strength (UTS), and elongation (EI) under different HT-PS conditions.It can be seen in Figure 6a that the YS and UTS for HT-PS3% and HT-PS7% remained relatively stable, while the average values for yield strength and tensile strength at HT-PS15% and 18% were noticeably higher.It can be seen in Figure 6b that the average values of EI decrease with an increase in HT-PS.Therefore, HT-PS has a positive correlation with strength and a negative correlation with elongation.Specifically, when HT-PS increased from 15% to 18%, the heat treatment temperature was set to 200 • C, and the heat treatment time was 30 min, the yield strength significantly increased from 229 MPa to 257 MPa, and elongation decreased from 19% to 9%.This indicated that higher HT-PS combined with the appropriate heat treatment temperature and time could significantly enhance the material's yield strength, but it could decrease the elongation dramatically.that the YS and UTS for HT-PS3% and HT-PS7% remained relatively stable, while the av-erage values for yield strength and tensile strength at HT-PS15% and 18% were noticeably higher.It can be seen in Figure 6b that the average values of EI decrease with an increase in HT-PS.Therefore, HT-PS has a positive correlation with strength and a negative correlation with elongation.Specifically, when HT-PS increased from 15% to 18%, the heat treatment temperature was set to 200 °C, and the heat treatment time was 30 min, the yield strength significantly increased from 229 MPa to 257 MPa, and elongation decreased from 19% to 9%.This indicated that higher HT-PS combined with the appropriate heat treatment temperature and time could significantly enhance the material's yield strength, but it could decrease the elongation dramatically.To further study the impact of the heat treatment process on the mechanical properties, mechanical property tests were conducted on specimens that had undergone HT-PS and natural aging treatment (without the heat treatment process).The experimental results are shown in Table 4.For low HT-PS (3%), an increase in the heat treatment temperature and time did not significantly enhance the yield strength and tensile strength.This demonstrated that the strengthening effect of the heat treatment process was limited at lower levels of HT-PS.For high HT-PS (15%), the strengthening effect of the heat treatment process depends on the heat treatment temperatures and times.The heat treatment effect became more pronounced when the temperature was higher than 185 °C and the heat treatment time was 30 min.Under HT-PS (3% and 15%) conditions, the heat treatment process had a minor impact on the elongation of the AA6016 sheet.To further study the impact of the heat treatment process on the mechanical properties, mechanical property tests were conducted on specimens that had undergone HT-PS and natural aging treatment (without the heat treatment process).The experimental results are shown in Table 4.For low HT-PS (3%), an increase in the heat treatment temperature and time did not significantly enhance the yield strength and tensile strength.This demonstrated that the strengthening effect of the heat treatment process was limited at lower levels of HT-PS.For high HT-PS (15%), the strengthening effect of the heat treatment process depends on the heat treatment temperatures and times.The heat treatment effect became more pronounced when the temperature was higher than 185 • C and the heat treatment time was 30 min.Under HT-PS (3% and 15%) conditions, the heat treatment process had a minor impact on the elongation of the AA6016 sheet.Based on the analysis of the tensile test results, it was concluded that the heat treatment process had a negligible impact on the elongation of the samples.To further validate this conclusion, all the fracture surfaces of the tensile specimens were sectioned.The cut surfaces were then scanned using SEM, and various fracture morphologies were analyzed.Four fracture surface images with significantly different elongation values were selected for comparison, as shown in Figure 7.As shown in Figure 7a,b, when the HT-PS was 3%, the fracture morphologies of the different heat-treated specimens were similar, exhibiting ductile fractures with deep dimples.Similarly, as depicted in Figure 7c,d, when the HT-PS was 15%, the fracture morphologies of the different heat-treated specimens were also similar but with shallower dimples.By comparing images with different pre-strains, it was evident that with increasing pre-strain, the number of dimples decreased, and the dimples became shallower.Therefore, the influence of HT-PS on the material fracture mode was greater than that of heat treatment.

Influence of the Manufacturing Process on Microstructure and Properties
Based on the analysis of the tensile test results, it was concluded that the heat treatment process had a negligible impact on the elongation of the samples.To further validate this conclusion, all the fracture surfaces of the tensile specimens were sectioned.The cut surfaces were then scanned using SEM, and various fracture morphologies were analyzed.Four fracture surface images with significantly different elongation values were selected for comparison, as shown in Figure 7.As shown in Figure 7a,b, when the HT-PS was 3%, the fracture morphologies of the different heat-treated specimens were similar, exhibiting ductile fractures with deep dimples.Similarly, as depicted in Figure 7c,d, when the HT-PS was 15%, the fracture morphologies of the different heat-treated specimens were also similar but with shallower dimples.By comparing images with different pre-strains, it was evident that with increasing pre-strain, the number of dimples decreased, and the dimples became shallower.Therefore, the influence of HT-PS on the material fracture mode was greater than that of heat treatment.The tensile test results indicated that the heat treatment process significantly improved the strength of the AA6016 sheet under high HT-PS conditions.Figure 8 shows SEM-BSE images of the second phase under partial processing conditions.According to the SEM-EDX spectra, it was observed that these micron-scale second phases were iron-rich phases, with the EDX results shown in Table 5.To further analyze the size and distribution of the micron-scale second phases (iron-rich phases) in Figure 8, quantitative statistics were calculated, as presented in Table 6.
The tensile test results indicated that the heat treatment process significantly improved the strength of the AA6016 sheet under high HT-PS conditions.Figure 8 shows SEM-BSE images of the second phase under partial processing conditions.According to the SEM-EDX spectra, it was observed that these micron-scale second phases were ironrich phases, with the EDX results shown in Table 5.To further analyze the size and distribution of the micron-scale second phases (iron-rich phases) in Figure 8, quantitative statistics were calculated, as presented in Table 6.In Table 6, it was observed that the quantified statistical results on the second phase at different positions on the same specimen were very close.Additionally, there was little In Table 6, it was observed that the quantified statistical results on the second phase at different positions on the same specimen were very close.Additionally, there was little difference in the proportion of the second phases between process 1 (3-200-20) and process 2 (15-200-30).The number density of the second phase in process 2 was slightly higher than that in process 1, and the average size of the second phase in process 2 was also slightly higher than that in process 1.However, there was a significant difference in the mechanical properties of the materials obtained from process 1 and process 2, despite the small difference in the statistical data on the micron-scale second phases.This indicated that there was no apparent relationship between the micron-scale second phases and the mechanical properties of the material, suggesting the need for further research on nanoprecipitates.
To study the effect of HT-PS on the precipitation phases in AA6016, DSC experiments were conducted, and the results are depicted in Figure 9.In the HT-PS samples, it was observed that the dissolution trough between 200 • C and 220 • C, arising from cluster reversion during natural aging, decreased or disappeared [25,26].This indicates that HT-PS treatment effectively inhibits natural aging [27].Moreover, as the HT-PS increased, the precipitation peak of the β ′′ phase in the samples shifted towards lower temperatures.The precipitation temperature of the β ′′ phase in the HT-PS 3% sample (248.24• C) was higher than that in the HT-PS 15% sample (242.85 • C), suggesting that a higher degree of HT-PS promotes the precipitation of the β ′′ phase.difference in the proportion of the second phases between process 1 (3-200-20) and process 2 (15-200-30).The number density of the second phase in process 2 was slightly higher than that in process 1, and the average size of the second phase in process 2 was also slightly higher than that in process 1.However, there was a significant difference in the mechanical properties of the materials obtained from process 1 and process 2, despite the small difference in the statistical data on the micron-scale second phases.This indicated that there was no apparent relationship between the micron-scale second phases and the mechanical properties of the material, suggesting the need for further research on nanoprecipitates.To study the effect of HT-PS on the precipitation phases in AA6016, DSC experiments were conducted, and the results are depicted in Figure 9.In the HT-PS samples, it was observed that the dissolution trough between 200 °C and 220 °C, arising from cluster reversion during natural aging, decreased or disappeared [25,26].This indicates that HT-PS treatment effectively inhibits natural aging [27].Moreover, as the HT-PS increased, the precipitation peak of the ′′ phase in the samples shifted towards lower temperatures.The precipitation temperature of the ′′ phase in the HT-PS 3% sample (248.24°C) was higher than that in the HT-PS 15% sample (242.85 °C), suggesting that a higher degree of HT-PS promotes the precipitation of the ′′ phase.The precipitation modeling is based on classic theory [28], and an incoherent interface of the precipitate/matrix is assumed.The stationary nucleation rate is given by: is the Zeldovich factor,  * is the condensation rate of solute atoms in nuclei,  is the number of nucleation sites per unit volume,  is the Boltzmann constant,  is temperature.∆ is the energy barrier for nucleation.According to the references, heterogeneous nucleation on dislocations releases the core energy of the dislocations, and hence the nucleation barrier is reduced [29,30].By analyzing the reduction in the nucleation barrier described by Equation ( 6), under the premise of the same material and the same precipitated phase, where ,  * , ,  , and  are all constant, only the nucleation barrier ∆ decreases.Therefore, the nucleation rate  of the precipitated phase will increase.The precipitation modeling is based on classic theory [28], and an incoherent interface of the precipitate/matrix is assumed.The stationary nucleation rate is given by: Z is the Zeldovich factor, β * is the condensation rate of solute atoms in nuclei, N is the number of nucleation sites per unit volume, k B is the Boltzmann constant, T is temperature.∆G nucl is the energy barrier for nucleation.According to the references, heterogeneous nucleation on dislocations releases the core energy of the dislocations, and hence the nucleation barrier is reduced [29,30].By analyzing the reduction in the nucleation barrier described by Equation ( 6), under the premise of the same material and the same precipitated phase, where Z, β * , N, k B , and T are all constant, only the nucleation barrier ∆G nucl decreases.Therefore, the nucleation rate J s of the precipitated phase will increase.Consequently, with an increase in HT-PS, the dislocation density also increases, indicating an increase in the nucleation rate of the precipitate phase.This increased nucleation rate results in the precipitation of a greater amount of the β ′′ phase.
The numerical method employed to track the evolution of the precipitates, including nucleation, growth, and coarsening, is the numerical Kampmann-Wagner (KWN) model [31].In this kinetic simulation model, a crucial parameter for calculating the nucleation rate of precipitates is the change in Gibbs free energy associated with the formation of critical nuclei.The total energy change resulting from nucleation can be expressed as follows: r represents the radius of the precipitate phase, γ represents the interfacial energy, and ∆G 0 bulk represents the volume free energy.Based on the aforementioned Equation ( 7), under the premise of the same material and the same precipitated phase, ∆G nucl and γ are constants.Analyzing the reduction in nucleation potential energy indicates that the corresponding precipitate radius decreases.Therefore, with an increase in HT-PS, the reduction in the nucleation barrier leads to a decrease in the precipitation radius.This suggests that higher HT-PS can promote the formation of finer precipitates.
During the HT-PS treatment process, an increase in the dislocation density within the crystals occurred [32].The dislocations introduced by straining acted as vacancy sinks.These vacancy sinks played a crucial role in the aggregation activity during natural aging, as the relaxation of over-saturated vacancies could significantly suppress the formation of Cluster (1) during natural aging.Serizawa et al. [33] reported two types of nanoclusters, denoted as Cluster (1) and Cluster (2), formed near room temperature and at 100 • C, respectively.The Cluster (1) structure had a wider range of Mg/Si ratios and was nearly insoluble at heat treatment temperature (180-200 • C).Cluster (2) could form above 100 • C, suppressing the formation of Cluster (1) during room temperature storage, and over-saturated vacancies inhibited natural aging, thus inhibiting the formation of Cluster (1).Cluster (2) had a Mg/Si ratio similar to the β ′′ phase, facilitating the formation of the β ′′ phase during the heat treatment process.In Al-Mg-Si alloys, the β ′′ phase was the primary strengthening phase [34,35].During heat treatment, dislocations could provide heterogeneous nucleation sites for Cluster (2) and accelerate β ′′ phase precipitation [34].Meanwhile, as the HT-PS increased, the dislocation density increased, and more dislocations led to a decrease in the nucleation barrier [29,30], combined with the previous Equations ( 8) and (9), resulting in an increase in the nucleation rate of the precipitates and a decrease in their radius with increasing HT-PS.Jia [12] studied and compared precipitates in samples with HT-PS of 7% and RT-PS of 7% (mainly the β ′′ phase), revealing that the precipitates in the hightemperature pre-strained samples had a higher volume fraction and a smaller average radius than those in the room-temperature pre-strained samples, resulting in a higher BHR.Therefore, with increasing HT-PS, the precipitates exhibited a higher volume fraction and a smaller average radius, leading to a significant increase in the mechanical properties of AA6016.
Research by Zandbergen [14] et al. indicated that increasing the heat treatment temperature and time promoted the precipitation of the β ′′ phase.Combined with the finding of the previous discussion, where the heat treatment process facilitated the formation of Cluster (2), and Cluster (2) promoted the precipitation of the β ′′ phase [34].Due to the hindrance of dislocation motion by the β ′′ phase.Therefore, greater force was required for dislocations to bypass these β ′′ phases, resulting in an improvement in the mechanical properties.Therefore, when the percentage of precipitates was higher and their average size was smaller, it was easier to impede dislocation motion, thus also facilitating the improvement of material strength.
In summary, the combined effect of HT-PS treatment and the heat treatment process on the mechanical properties of AA6016 was investigated.With increasing HT-PS treatment, as well as an increasing heat treatment temperature and time, a higher volume fraction and a smaller average radius of the β ′′ phase were precipitated.The increased presence of the β ′′ phase hindered dislocations, consequently enhancing the mechanical properties of AA6016.The overall mechanism is illustrated in Figure 10.
treatment, as well as an increasing heat treatment temperature and time, a higher volume fraction and a smaller average radius of the ′′ phase were precipitated.The increased presence of the ′′ phase hindered dislocations, consequently enhancing the mechanical properties of AA6016.The overall mechanism is illustrated in Figure 10.

Comparison of Neural Network Model Fitting Accuracy
Figure 11 depicts a comparison between the experimental and fitted stress-strain curves of AA6016.In order to objectively assess the prediction accuracy of the two models, random experimental data were selected for each model.As depicted in Figure 11, the fitting accuracy of the two models was quite similar.To further quantitatively evaluate the fitting performance, this study utilized statistical metrics such as the coefficient of determination (R2), mean absolute error (MAE), and root mean square error (RMSE).These metrics provided a comprehensive assessment of the model's performance: = ∑ ( −  )  (10) where  is the experimental value;  is the fitted value;  is is the mean of the experimental value; N is the amount of experimental data.The fitting values of the models were comprehensively evaluated using the coefficient of determination (R2), mean absolute error (MAE), and root mean square error (RMSE).Figure 11 depicts a comparison between the experimental and fitted stress-strain curves of AA6016.In order to objectively assess the prediction accuracy of the two models, random experimental data were selected for each model.As depicted in Figure 11, the fitting accuracy of the two models was quite similar.To further quantitatively evaluate the fitting performance, this study utilized statistical metrics such as the coefficient of determination (R2), mean absolute error (MAE), and root mean square error (RMSE).These metrics provided a comprehensive assessment of the model's performance: where σ i E is the experimental value; σ i P is the fitted value; σ E is is the mean of the experimental value; N is the amount of experimental data.racy.To further validate the accuracy of the GA-BP model in predicting the mechanical properties, yield strength, tensile strength, and elongation were used as the fitting results.Table 7 shows a comparison between the fitting results of AA6016 and the experimental data.From Table 7, it can be observed that the data fitted by the GA-BP model approximated the experimental data, significantly outperforming the BP model.The GA-BP artificial neural network model demonstrated good fitting accuracy.The fitting values of the models were comprehensively evaluated using the coefficient of determination (R2), mean absolute error (MAE), and root mean square error (RMSE).Figure 12 presents the statistical values for both models.The correlation coefficients of the GA-BP and BP models were 0.986 and 0.985, respectively, indicating that the two models had similar correlations.The GA-BP model had a slightly lower MAE of 2.147 compared to the BP model.Additionally, the RMSE of the GA-BP model was also lower than that of the BP model, indicating that the GA-BP model had an advantage in overall fitting accuracy.To further validate the accuracy of the GA-BP model in predicting the mechanical properties, yield strength, tensile strength, and elongation were used as the fitting results.Table 7 shows a comparison between the fitting results of AA6016 and the experimental data.From Table 7, it can be observed that the data fitted by the GA-BP model approximated the experimental data, significantly outperforming the BP model.The GA-BP artificial neural network model demonstrated good fitting accuracy.

Comparison of Model Prediction Accuracy
The trained neural network models have the capability to predict data, and in this study, two types of predictions were conducted for each model.Firstly, the models were utilized to predict the experimental data within the range of the training set (15-170-20, the data were not used as a training set to train the model).Secondly, the models were used to predict the data points outside the range of the training set (18-200-30), as illustrated in Figure 13.It was observed that the GA-BP model outperformed the BP model in terms of its prediction accuracy.To further evaluate the prediction performance of the two models, the aforementioned Equations ( 7)-( 9) were employed to provide a comprehensive assessment of their performance, as presented in Table 8.Table 8 showcases the statistical values of the predictions generated by the two models.The GA-BP model achieved a higher correlation coefficient compared to the BP model, indicating a stronger relationship between the predicted and actual values.Additionally, the GA-BP model exhibited smaller mean absolute error (MAE) and root mean square error (RMSE) values than the BP model, suggesting the better prediction accuracy of the GA-BP model.

Comparison of Model Prediction Accuracy
The trained neural network models have the capability to predict data, and in this study, two types of predictions were conducted for each model.Firstly, the models were utilized to predict the experimental data within the range of the training set (15-170-20, the data were not used as a training set to train the model).Secondly, the models were used to predict the data points outside the range of the training set (18-200-30), as illustrated in Figure 13.It was observed that the GA-BP model outperformed the BP model in terms of its prediction accuracy.To further evaluate the prediction performance of the two models, the aforementioned Equations ( 7)-( 9) were employed to provide a comprehensive assessment of their performance, as presented in Table 8.Table 8 showcases the statistical values of the predictions generated by the two models.The GA-BP model achieved a higher correlation coefficient compared to the BP model, indicating a stronger relationship between the predicted and actual values.Additionally, the GA-BP model exhibited smaller mean absolute error (MAE) and root mean square error (RMSE) values than the BP model, suggesting the better prediction accuracy of the GA-BP model.Figure 14 illustrates a comparison between the experimental and predicted values of the mechanical properties.In Figure 14a, data prediction is conducted for the experimental values within the training set range (15-170-20).The GA-BP model predicts a yield strength of 173.6 MPa, which is 5.4 MPa lower than the experimental value.On the other hand, the BP model predicts a yield strength that is 36.8MPa higher than the experimental value.The GA-BP model provides predicted values of tensile strength and elongation that are closer to the experimental values compared to the BP model.Overall, both models perform reasonably well in predicting the mechanical properties within the training set range, but the GA-BP model exhibits a better prediction accuracy.In Figure 14b, data prediction is carried out for the experimental values outside the training set range (18-200-20).The results indicate that the GA-BP model offers significantly better predictions for both yield strength and tensile strength compared to the BP model.However, both models show notable deviations in predicting elongation.In summary, the GA-BP model demonstrates a better prediction accuracy for both sets of results, particularly in terms of yield strength and tensile strength.Figure 14 illustrates a comparison between the experimental and predicted values of the mechanical properties.In Figure 14a, data prediction is conducted for the experimental values within the training set range (15-170-20).The GA-BP model predicts a yield strength of 173.6 MPa, which is 5.4 MPa lower than the experimental value.On the other hand, the BP model predicts a yield strength that is 36.8MPa higher than the experimental value.The GA-BP model provides predicted values of tensile strength and elongation that are closer to the experimental values compared to the BP model.Overall, both models perform reasonably well in predicting the mechanical properties within the training set range, but the GA-BP model exhibits a better prediction accuracy.In Figure 14b, data prediction is carried out for the experimental values outside the training set range (18-200-20).The results indicate that the GA-BP model offers significantly better predictions for both yield strength and tensile strength compared to the BP model.However, both models show notable deviations in predicting elongation.In summary, the GA-BP model demonstrates a better prediction accuracy for both sets of results, particularly in terms of yield strength and tensile strength.

Conclusions
This study investigated the mechanical properties of AA6016 sheets after undergoing the HFQ process and subsequent heat treatment.The following conclusions could be drawn:

Conclusions
This study investigated the mechanical properties of AA6016 sheets after undergoing the HFQ process and subsequent heat treatment.The following conclusions could be drawn: (1) At HT-PS of 3% and 7%, the heat treatment process has a minimal effect on the mechanical properties of the AA6016 aluminum alloy.At HT-PS of 15%, the mechanical strength of the AA6016 sheet improved with an increasing heat treatment temperature and time.A further increase in HT-PS from 15% to 18% could reduce the elongation significantly.
(2) Micron-scale second phases (iron-rich phases) underwent quantitative analysis using SEM-BSE.The results revealed that these micron-scale second phases exhibited no significant impact on the mechanical properties of the AA6016 sheets.The DSC results indicated that as the HT-PS increased from 3% to 15%, the precipitation peak of the β ′′ phase decreased.This suggests that HT-PS can promote the precipitation of the β ′′ phase and enhance the strength of the AA6016 sheet.(3) Through analysis of the nucleation rate formula and the Gibbs free energy formula in the KWN model, it can be concluded that with an increase in HT-PS, the precip-itated phase will obtain a higher volume fraction and a smaller mean radius, and AA6016 sheets will obtain a higher strength.(4) Comparing the fitting and prediction results of the BP neural network model and the GA-BP neural network model, it is evident that the GA-BP neural network model outperforms the other in terms of fitting and prediction accuracy.
A comprehensive assessment was conducted on the influence of the HFQ process and heat treatment techniques on the mechanical properties of AA6016 sheets.This provided a process design guideline for using HFQ and heat treatment processes on AA6016 in the automotive industry.

Figure 2 .
Figure 2. Schematic representation of heat treatment procedures for AA6016 alloy.

Figure 2 .
Figure 2. Schematic representation of heat treatment procedures for AA6016 alloy.

Figure 3 .
Figure 3. Structure of a BP neural network.

Figure 3 .
Figure 3. Structure of a BP neural network.

Figure 5
Figure5illustrates the stress-strain curves of AA6016 sheets under different heat treatment treatments after HT-PS.In the diagram, the first set of numbers(3,7,15,18) represents the HT-PS values, the second set of numbers (170, 185, 200) indicates the heat treatment temperatures, and the third set of numbers(10,20,30) denotes the heat treatment durations.In Figure5a, for HT-PS (HT-PS) of 3%, different heat treatments had a minimal effect on the strength of the AA6016 sheet but slightly impacted the elongation.Figure5bshowed that for HT-PS of 7%, particularly with a heat treatment temperature of 185 • C for 30 min, the strength of the AA6016 sheet increased.In Figure5c, for an HT-PS of 15%, there was a noticeable increase in strength with an increase in the heat treatment temperature and time, but the elongation decreased significantly.

Figure 5
Figure5illustrates the stress-strain curves of AA6016 sheets under different heat treatment treatments after HT-PS.In the diagram, the first set of numbers(3,7,15,18) represents the HT-PS values, the second set of numbers (170, 185, 200) indicates the heat treatment temperatures, and the third set of numbers(10,20,30) denotes the heat treatment durations.In Figure5a, for HT-PS (HT-PS) of 3%, different heat treatments had a minimal effect on the strength of the AA6016 sheet but slightly impacted the elongation.Figure5bshowed that for HT-PS of 7%, particularly with a heat treatment temperature of 185 °C for 30 min, the strength of the AA6016 sheet increased.In Figure5c, for an HT-PS of 15%, there was a noticeable increase in strength with an increase in the heat treatment temperature and time, but the elongation decreased significantly.

Figure 5 .
Figure 5. Stress-strain curve with (a) HT-PS 3%, (b) HT-PS 7%, (c) HT-PS 15% and 18%.The mechanical properties under different experimental conditions are shown in Table 3, and Figure6shows a relationship graph of the yield strength (YS), tensile strength (UTS), and elongation (EI) under different HT-PS conditions.It can be seen in Figure6a

Figure 9 .
Figure 9. DSC curves of samples under different pre-straining conditions.

Figure 9 .
Figure 9. DSC curves of samples under different pre-straining conditions.

Figure 10 .
Figure 10.Flowchart illustrating the influence of manufacturing processes on material microstructure.

Figure 10 .
Figure 10.Flowchart illustrating the influence of manufacturing processes on material microstructure.

4. 3 .
Comparative Analysis of Neural Network Model Fitting Accuracy and Prediction Accuracy 4.3.1.Comparison of Neural Network Model Fitting Accuracy

Figure 12 .
Figure 12.Statistical values for both models (BP and GA-BP models).

Table 2 .
HT-PS and heat treatment process.

Table 2 .
HT-PS and heat treatment process.

Table 3 .
Mechanical properties under different experimental conditions.

Table 3 .
Mechanical properties under different experimental conditions.

Table 6 .
Effect of heat treatment process on the second phase of HT-PS 3% and HT-PS 15%.

Table 6 .
Effect of heat treatment process on the second phase of HT-PS 3% and HT-PS 15%.

Table 7 .
Fitted values for both models (BP and GA-BP models).

Table 8 .
Predictive statistics for both models (BP and GA-BP models).

Table 8 .
Predictive statistics for both models (BP and GA-BP models).

Table 8 .
Predictive statistics for both models (BP and GA-BP models).