Abstract
Zn–Ni electrophosphate coating is one of the most commonly used materials in industrial applications. The corrosion resistance of this coating is very important in order to achieve the minimum corrosion current of the Zn–Ni electrophosphate coating. This study described a new reliability simulation framework to determine the corrosion behavior of coating using a gene artificial neural network (ANN) to estimate the corrosion current of the coating. The input parameters of the model are temperature, pH of electroplating bath, current density, and Ni2+ concentration, and corrosion current defined as output. The effectiveness and accuracy of the model were checked by utilizing the absolute fraction of variance (R2 = 0.9999), mean absolute percentage error (MAPE = 0.0171), and root mean square error (RMSE= 0.0002). This is determined using the genetic algorithm (GA) and the optimum practice condition.
Introduction
Hot-dip galvanization is an efficient safety measure of steel against the corrosion of atmospheric corrosion, but it cannot offer adequate protection of the steel substrates against the corrosion substance produced in the atmosphere (; ; ). In general, the surface modification treatment is applied to galvanized steel for greater protection of zinc-coated steel from corrosion. Phosphating, which has been used extensively in a wide variety of industries, is one of the major processes of chemical conversion for applications such as corrosion prevention, paint primers, wear prevention, metal forming lubricants, electrical insulation, and even decoration (; ). The highest operating temperatures of the most traditional phosphating baths mentioned in the literature are 90–98°C. A high level of energy consumption is the biggest problem for the use of a high temperature bath, which is a huge crisis nowadays. In addition, due to the development of a scale that contributes to excessive heating of the bath solution, the application and retaining of the heating coils is difficult, and thus a repeated replacement is needed (). The overheating of bath solution, allowing an early conversion of primary phosphate into tertiary phosphate before the metal is treated, is another issue in this case. These reasons increase the free acidity of the bath and thus delay the precipitation of the coating of phosphates. The low-temperature phosphating processes have become more promising and successful due to the high cost of energy in the high-temperature bath. However, the phosphating bath at low temperatures is slow, and certain techniques must be used to accelerate. In order to facilitate the phosphating process, chemical, electrochemical, and mechanical methods are used. The chemical methods are most widely used in many industries. However, the chemical accelerators such as nitrite are poisonous and have a detrimental effect on the environment and human health, and it is desirable to be eliminated from the bath composition and replaced with some replacements such as electrochemical accelerator (; ).
In many industries, the chemical processes are most commonly used. However, chemical accelerators such as nitrite are toxic and have a harmful effect on the environment. So, it is required to restrict the use of these hazardous chemicals by replacing the electrochemical accelerator (; ; ; ; Zhu et al., 2021). The electrochemical methods using cathodic and anodic treatments for accelerating phosphating processes become popular these days for their low operating temperature, high speed of coating, and improved coating properties. So, experimental parameters will affect the structural and chemical properties as well as electroplating efficiency during electrodeposition of an electrophosphating coating. Typical methods and measurements make it very difficult to analyze the results of these deposition parameters and to consider the relations between the performance characteristics.
The use of the artificial neural network (ANN) and the genetic algorithm is an engineering technique to improve process variables. The ANN is one of the important non-linear strategies to simulate the complicated behavior of materials (; Zhu et al., 2021). This method can learn and simulate the experiential information through weight and bias in the neuron, but cannot determine the optimal condition and may be hindered at the local minimum.
Genetic algorithms (GA) are used to find the optimum solution(s) for a computational problem that maximizes or minimizes a certain function. Genetic algorithms are also called as evolutionary computation, which mimic biological reproductive and selection processes in order to find the “fittest” lst solutions. Artificial intelligence (AI) techniques such as ANN and gene expression programming (GEP) are important in many areas of physical, chemical, and mathematical science (). The ANN method offers a modern technique for predicting material deterioration in different conditions. The ANN is a technology for artificial intelligence to simulate human brain biological processes (; ). The system consists of one-way signal flow channels of interconnected operators. Samples are obtained with a distributed code that forms a non-linear, feasible framework. It is also environmentally self-adapting to react rationally to the different inputs (; ). ANN models have few drawbacks of local minima, and to prevent local minima, which leads to false convergence of the ANN model, a derivative-free optimization algorithm should be added in the ANN algorithm training process (). There are a lot of studies about the applications of GEP in the literature for different engineering problems (). The combination of regression strategies and systematic design of experiment is an efficient alternative approach for providing the experimental data in a new popular model approach such as GEP. ANN and GEP can capture complex interactions among input/output variables in a system without using prior knowledge about the nature of these interactions. Isleem et al. reported a database based on the ANN for predicting the hardening and softening behavior of glass fiber-reinforced polymers (). Finke et al. have developed anticorrosion behavior of coating based on the ANN (). Hu et al. reported on the corrosion behavior of Mo–V–Cr–Ni steel using an ANN model to predict polarization curves in deep-sea environments without the need for experiments ().
To the best of our knowledge, very few studies are reported about ANN and GEP together to compare the expected results and explain the research protocol (; ; ; ; ). The purpose of this study was to assess the oxidation behavior of Ni base alloys using particle swarm optimization (PSO)-ANN and GEP models.
Due to their flexibility in general regression processes, the present study was conducted using combined GA–ANN to determine an appropriate process for zinc–nickel electrophosophate coating on galvanized steel conditions.
Materials Production and Processing
The thickening of zinc covering on a steel layer of galvanized stain was 10 μm, and the dimension of hot-dip galvanized steel with a dimension (50 × 20 × 2 mm) was used in this experiment. It was ultrasonically degassed for 12 min in acetone, and then washed with deionized water for 20 s in 4% sulfuric acid. A dip into a titanium phosphate solution (1 g/L) for 1 min was used to enable the samples. Finally, prepared specimens were immersed in a cell-contained phosphating solution under different coating conditions.
Constant direct current was used for electrophosphating cells using a galvanostat. Table 1 indicates the chemical composition and operating conditions of the bath (; ).
TABLE 1
| Parameters | Amount | |
|---|---|---|
| Bath composition | Zinc oxide (g/L) | 2.04 |
| Orthophosphoric acid (ml/L) | 11 | |
| Sodium fluoride (g/L) | 0.3 | |
| Nickel nitrate | 1 | |
| Operating conditions | pH | 1.6–2.2 |
| Temperature (°C) | 25–35 | |
| Ni2+ (g/L) | 0–0.5 | |
| Current density (mA/Cm2) | 5–30 |
Experimental conditions.
In this analysis, anodes of stainless steel were used. Two sets of galvanized steel substrates have been mounted on both ends. The area under the review was one square centimeter, and the remainder was sealed. The cell for corrosion measurement comprises three electrodes, namely, the operating electrode, stainless steel plate, and reference saturated calomel electrode (SCE). The potentiodynamic polarization experiments were performed with a scanning speed of 1 mV/S ranging from - 300 mV (relating to the open circuit potential) to +300 mV. The corrosion current density (icorr) was evaluated using the linear polarization method.
Artificial Neural Network
The ANN, in particular the structure of the brain, is influenced by biological neural networks. An ANN model includes elements such as inputs and works in a manner that creates a relation between the elements in an approximation to what a brain is doing in other words, it works by learning the connection between input and output data (). A connected neuron formula in terms of weight is reported as Eq. 1:
Herein, b represents the bias of neurons, p represents the number of elements, and wi represents the weight of the input vector (ai). Each neuron receives sum of the weight inputs with bias and is used from the activation function to validate its output signal, as Eq. 2:
where represents the transfer function of the neuron. Some of the most common transfer functions are log-sigmoid (Logsig), tan-sigmoid (Tansig), and linear transfer function (Purelin), followed by Eqs 2–4.
Figure 1 schematically shows these transfer function. Feed-forward back-propagation is one of the most common strategies for training of the ANN. This approach enables us to present effective solutions for engineering applications, especially in material science (; ; ).
FIGURE 1
Detailed explanations of the functions of hidden and output layers are descripted in the work by
GA
As Darwin’s theory explains, the GA is a random global method of search and optimization which mimics natural biological evolution. Using a population of individuals, the GA examines all regions of the solution space (
PSO
PSO is one of the most popular population-based stochastic optimization algorithms (
FIGURE 2

Pseudo-code of standard PSO.
The aim of the current study was to determine the most appropriate values for the weights and bias of NN (i.e., optimized NN) based on the PSO algorithm (
GEP
According to the literature, the use of AI, for example, the adaptive neuron fuzzy inference system (ANFIS), and the ANN in combination with optimization techniques, for example, the GA, simulated annealing (SA), are powerful techniques for modification and simulation of engineering phenomena (
ETs are combined from constants, functions, operators, and input variables with their nodes and terminals. A mathematical function is responsible to connect each other ETs as shown in Figure 3 (
FIGURE 3

Typical representation of linear string and their GEP model.
Figure 4 illustrates the GEP operators, schematically. As shown, the generation of chromosomes in a random manner is the first stage of GEP. In the second stage, the validation of each individual population is investigated by the use of fitness functions as the criteria. The generation of new output is performed by the usage of reproduction, crossover, and mutation operators in GEP.
FIGURE 4

Schematical representation of GEP flowchart (
During the regeneration stages, functions with lower performance are omitted after the adoption of the algorithm. Selection and combination of most appropriate data from parents were carried out using a swamp randomly selected section of two trees. Mutation is responsible to investigate the non-local properties during simulation (
Results and Discussion
In this study, a feed-forward neural network with a back-propagation algorithm is used to simulate the process. The output of each neuron in the feed-forward technique is simply linked with the next layer. This strategy consists of one or more hidden layers and has one output layer. To increase the performance of evolution, all input data were normalized in the range 0–1. Hidden layers control the network by using a non-linear transfer function to learn linear or non-linear behavior between inputs and outputs. Normalizing is performed using Eq. 6 (
where Xmax and Xmin are the highest and lowest limits, respectively. In this work, the back-propagation used by a network that has an input layer with its operating condition factor neurons (Table 1) and an output layer that has one neuron (icorr); 47 experimental result patterns used to train the ANN model are seen in Table 2, and the remainder is used for research. As shown in Table 3, by consideration of regression as the criteria, “Logsig” and “Purelin” are the best functions for hidden and output layers, respectively.
TABLE 2
| No | Input | Output | |||
|---|---|---|---|---|---|
| Temperature (°C) | pH | Current density (mA/Cm2) | Ni2+ (g/L) | icorr (mA/Cm2) | |
| 1 | 25 | 1.60264 | 30.2635 | 0 | 11.0039 |
| 2 | 25.9151 | 1.60562 | 30.0475 | 0.075731 | 10.3807 |
| 3 | 25.7118 | 1.61354 | 30.1443 | 0.075731 | 10.3923 |
| 4 | 26.2202 | 1.62157 | 29.6782 | 0.075731 | 10.195 |
| 5 | 26.6267 | 1.62708 | 29.6468 | 0.074311 | 9.78467 |
| 6 | 27.1352 | 1.63259 | 29.3698 | 0.106615 | 9.37438 |
| 7 | 28.0502 | 1.64337 | 29.1583 | 0.138208 | 8.97185 |
| 8 | 28.2535 | 1.65404 | 29.1583 | 0.159744 | 8.77836 |
| 9 | 28.5586 | 1.66471 | 28.9736 | 0.18128 | 8.58486 |
| 10 | 28.9653 | 1.67825 | 28.5774 | 0.223997 | 7.97719 |
| 11 | 29.5753 | 1.69179 | 28.2959 | 0.244113 | 7.36953 |
| 12 | 29.7785 | 1.69718 | 28.2959 | 0.243758 | 7.16826 |
| 13 | 30.1852 | 1.70785 | 28.269 | 0.286475 | 6.97477 |
| 14 | 30.7953 | 1.71875 | 27.7805 | 0.296888 | 6.36322 |
| 15 | 32.219 | 1.72976 | 27.5779 | 0.327417 | 5.54264 |
| 16 | 32.5238 | 1.74319 | 27.4587 | 0.370844 | 5.144 |
| 17 | 32.5238 | 1.7511 | 27.4587 | 0.370844 | 5.15565 |
| 18 | 33.0323 | 1.76177 | 27.0893 | 0.370134 | 4.96215 |
| 19 | 33.4389 | 1.77256 | 26.7944 | 0.379837 | 4.55963 |
| 20 | 33.9474 | 1.78598 | 5.17569 | 0.401373 | 4.16099 |
| 21 | 34.2524 | 1.79677 | 5.95284 | 0.433676 | 3.75846 |
| 22 | 34.2524 | 1.80732 | 5.95284 | 0.433676 | 3.77399 |
| 23 | 34.8624 | 1.82327 | 6.29077 | 0.477103 | 3.58826 |
| 24 | 34.8624 | 1.83118 | 6.29077 | 0.477103 | 3.5999 |
| 25 | 35.2692 | 1.92124 | 9.86549 | 0.497929 | 3.10482 |
| 26 | 34.8624 | 1.92904 | 26.0379 | 0.497929 | 3.32549 |
| 27 | 34.8624 | 1.93696 | 26.0379 | 0.497929 | 3.33713 |
| 28 | 34.2524 | 1.95256 | 5.95284 | 0.433676 | 3.77848 |
| 29 | 34.2524 | 1.96048 | 5.95284 | 0.433676 | 3.79012 |
| 30 | 34.2524 | 1.97631 | 5.43026 | 0.433676 | 3.81342 |
| 31 | 34.2524 | 1.9895 | 5.43026 | 0.433676 | 3.83283 |
| 32 | 33.9474 | 1.99994 | 5.17569 | 0.401373 | 4.05738 |
| 33 | 33.9474 | 2.03688 | 5.17569 | 0.401373 | 4.11173 |
| 34 | 33.4389 | 2.04996 | 26.7944 | 0.379837 | 4.34017 |
| 35 | 33.0323 | 2.06028 | 27.0893 | 0.370134 | 4.77375 |
| 36 | 33.0323 | 2.07875 | 27.0893 | 0.370134 | 4.80092 |
| 37 | 32.3206 | 2.10755 | 27.4587 | 0.33783 | 5.26168 |
| 38 | 32.219 | 2.11798 | 27.4587 | 0.327417 | 5.48623 |
| 39 | 31.8121 | 2.12315 | 27.5779 | 0.327417 | 5.70302 |
| 40 | 31.2018 | 2.13095 | 27.5779 | 0.327417 | 5.92369 |
| 41 | 30.4901 | 2.15435 | 28.1498 | 0.286475 | 6.58571 |
| 42 | 29.5753 | 2.175 | 28.2959 | 0.244113 | 7.45287 |
| 43 | 28.2535 | 2.18773 | 28.9736 | 0.191338 | 8.30838 |
| 44 | 27.6438 | 2.18991 | 29.1807 | 0.106615 | 9.14837 |
| 45 | 28.9653 | 2.1906 | 28.5774 | 0.223997 | 7.89421 |
| 46 | 25.9151 | 2.20252 | 30.0475 | 0.075731 | 10.2129 |
Experimental data descriptions for GEP models (
TABLE 3
| No | Neurons | Function | R2 | RMSE | MAPE |
|---|---|---|---|---|---|
| 1 | 3-3-1 | logsig-purelin | 0.9999 | 0.0002 | 0.0170 |
| 2 | 3-3–1 | tansig-purelin | 0.998 | 0.0011 | 0.0315 |
| 3 | 3-5–1 | logsig-purelin | 0.9957 | 0.0026 | 0.0614 |
| 4 | 3-5–1 | tansig-purelin | 0.986 | 0.0025 | 0.0564 |
| 5 | 3-7–1 | logsig-purelin | 0.9808 | 0.0123 | 0.0654 |
| 6 | 3-7–1 | tansig-purelin | 0.9848 | 0.0019 | 0.0543 |
| 7 | 3-9–1 | logsig-purelin | 0.9905 | 0.021 | 0.0678 |
| 8 | 3-9–1 | tansig-purelin | 0.9955 | 0.0044 | 0.0682 |
Amount of regressions for different ANN structures.
RMSE: Root mean square error; MAPE: Mean absolute percentage error.
In Figure 7, it can be observed that different ANN models with different neurons in the hidden layer provide higher R2 of 0.9999, lower RMSE of 0.0002, and MAPE of 0.0171. Hence, this 3-3-1 topology is used in further analysis in this study, as shown in Figure 3.
The comparison between experimental and expected data is shown in Figure 5. The correlation coefficient for the training and evaluation processes is also 0.9999 and 0.9985, respectively. Neural network results confirmed high precision in the experimental data.
FIGURE 5

Result of ANN models: (A)R2, (B) RMSE, and (C) MAPE.
The close linear pattern between the ANN predicts the values for the performance parameters, and the experimental dataset shows the adjacency of the model (
where the level of sensitivity related to input is represented by Si (%), and N = 47 is the number of datasets.
Changes in the hardness of each input variable are shown in Figure 6. As shown, the magnitude of operating conditions such as pH, temperature, Ni2+, and current density has more positive impact on the corrosion current density (icorr) of products.
FIGURE 6

Schematical representation of the ANN model.
The corrosion current of Ni–Zn electrophosphate coating is a function of temperature. The corrosion current decreased as the temperature increased.
Only tiny needlelike crystals shaped on a substrate surface at room temperature do not cover the entire region. In general, the development of phosphate coating includes nucleation and growth phases. As the nucleation stage occurs at low temperatures, the tiny needles are formed, but the growth kinetics at this temperature is too low, and thus the growth process cannot be completed. As can be observed, the kinetic energy of the growth stage is high enough to complete the growth process by raising the temperature to 43°C, thereby increasing the corrosion current (
The corrosion current increases slightly when the pH increases from 1.8 to 2.2 and increases greatly up to 2.2. Increasing of corrosion current with increase of pH is related to the porosity of coating and thus decreasing of resistance of corrosion.
As the current density is increased to 22 mA/cm2, further increase in the current density is directly dependent on corrosion current. This increase in the current density results from a high over potential, and this leads to an increased nucleation rate. The increasing value of the current density also leads to Ni–Zn electrophosphate coating (
The presence of nickel ions in the bath decreases the size of phosphate coatings and improves the corrosion behavior of phosphate coatings. The improvement of corrosion resistance of phosphate coating in the presence of Ni2+ is due to the fact that nickel ions in the phosphatization bath is a place for nucleation of phosphate crystals which acts and reduces the size of phosphate coating seeds.
Genetic algorithms have been used for the choosing of a single point crossover and the roulette wheel. Each individual was developed using an ANN model to determine the fitness function. The crossover probability, mutation probability, starting population, and generation size have been chosen to be 0.8, 0.2, 50, and 100, respectively. According sensitivity analysis and GA simulation predicted the combination sintering temperature 8.35°C, pH 9.93, current density 8.64 mA/cm2, and Ni2+ 9.09 g/L to give the optimum corrosion current density (icorr) of 1.61 mA/cm2. Experiments in various setups were performed to test the results, and the results are described in Table 3. The experimentally obtained optimum fitness function was determined by number 1 closer to GA’s predicted conditions. The validity of the simulated results is conclusively confirmed.
The practical parameters of each GEP model are abbreviated in Table 4. The maximum number of genes chosen is equal to 4 to prevent an increase in complexity. Since 1 gene has a very weak performance, the number of genes chosen is equal to 3. The head size is changed continuously from 10 at first to 11. MAPE is more appropriate than RMSE for fitness function.
TABLE 4
| Operating conditions | Sample.1 | Sample.2 |
|---|---|---|
| pH | 1.99994 | 2.06028 |
| Temperature (°C) | 33.9474 | 33.0323 |
| Ni2+ (g/L) | 0.401373 | 0.370134 |
| Current density (mA/Cm2) | 5.17569 | 27.0893 |
Proposed optimum condition by GA.
Table 5 abbreviates the statistical characteristics of GEP models. Accordingly, R2 values for the training phase change between 0.9913–0.9981. As shown in Figure 7 the minimum value of MAPE and RMSE are 0.0435 and 0.1040, respectively, for training data.
TABLE 5
| Parameters | Values |
|---|---|
| No. of chromosomes | 30 |
| Head size | 10–11 |
| No. of genes | 3–4 |
| Linking function | Addition, multiplication |
| Fitness function error type | MAPE |
| Constant per gene | 1 |
| Rate of mutation | 0.0044 |
| Rate of inversion | 0.0015 |
| Recombination rate (RR) for one or two points | 0.0277 |
| Rate of gene recombination | 0.0155 |
| Rate of gene transposition | 0.0155 |
Different parameters of GEP.
RMSE: Root mean square error; MAPE: Mean absolute percentage error.
FIGURE 7

(A) Training data and (B) testing data using regression analysis of predicted and experimental hardness.
In Table 6, the best GEP model is GEP-2 with linking function addition, head size 10, and number of genes about 4. Changing the linking function to multiplication with the same head size and number of genes in the GEP-6 model shows the best model is GEP-2 with 30 chromosomes. Table 7 shows the differed chromosomes, head size, and genes number for GEP-2 models. The R2, RMSE, and MAPE of different models of GEP-2 are shown in Figure 8.
TABLE 6
| Model | Linking function | Head size | Number of genes | Type of function | R2 | RMSE | MAPE |
|---|---|---|---|---|---|---|---|
| GEP-1 | Addition | 10 | 3 | +, −, ,/ | 0.9970 | 0.1304 | 0.0545 |
| GEP-2 | Addition | 10 | 4 | +, −,,/,ln,, cos | 0.9981 | 0.1040 | 0.0435 |
| GEP-3 | Addition | 11 | 3 | , , tan, arcsin | 0.9966 | 0.1448 | 0.0606 |
| GEP-4 | Addition | 11 | 3 | , −, , cos, sin | 0.9931 | 0.1977 | 0.0827 |
| GEP-5 | Multiplication | 10 | 3 | +, −, ,/ | 0.9851 | 0.2919 | 0.1222 |
| GEP-6 | Multiplication | 10 | 4 | +, −, ,/, ln, cos | 0.9930 | 0.1993 | 0.0834 |
| GEP-7 | Multiplication | 11 | 3 | , , tan, arcsin | 0.9937 | 0.1893 | 0.0792 |
| GEP-8 | Multiplication | 11 | 3 | , −, , cos, sin | 0.9913 | 0.2245 | 0.0939 |
Statistics of GEP models.
TABLE 7
| Parameter | Value |
|---|---|
| Training sample | 35 |
| Testing sample | 11 |
| R2 training | 0.9981 |
| RMSE training | 0.1040 |
| MAPE training | 0.0435 |
| R2 testing | 0.9876 |
| RMSE testing | 0.3721 |
| MAPE testing | 0.1658 |
Statistical data and parameters for selected GEP model.
FIGURE 8

Significances of input variables in mechanical properties.
The best GEP model is Model 1 as shown in Table 7, according to R2 values. In Figure 9, this model expression tree is provided for the four subgenes of the proposed model.
FIGURE 9

Result of GEP models: (A)R2, (B) RMSE, and (C) MAPE.
As shown in Figure 10, the expression tree of the GEP-2 and formulation of this shows
FIGURE 10

Expression tree of the GEP-2.
Conclusion
The ANN model with one hidden layer with three neurons in the hidden layer is a useful method for the prediction of operating conditions to the corrosion current density (icorr). The combined GA–ANN algorithm was an effective model for optimizing icorr parameters to the minimum of operating conditions. Sensitivity analysis showed that Ni2+ concentration and temperature are the most significant parameter, and Ni2+ concentration is the most important parameter.
This study proposed the most appropriate models by the employment of GEP for the estimation of operating conditions to the minimum corrosion current density (icorr). One type of function, that is, monotonic and non-monotonic, has been used within GEP. Then, by consideration of the lowest value of MAPE as the criteria, the best models were selected. The combination of sintering temperature 8.35°C, pH 9.93, current density 8.64 mA/cm2, and Ni2+ 9.09 g/L predicted by GA simulation and sensitivity analysis would result in an optimal corrosion current density (icorr) of 1.61 mA/cm2. In summary, the GEP models had a relatively high accuracy for the prediction of friction factor of multi-stage counter-current fluidized bed reactor, by consideration of the process parameters, operating conditions to the minimum corrosion current density (icorr).
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
We consciously assure that this manuscript is our own original work, which has not been previously published elsewhere. It is not currently being considered for publication elsewhere and also reflects its analysis in a truthful and in complete manner. We have properly given credits to the meaningful contributions of co-authors. The results are appropriately placed in the context of prior and existing research. All sources used are properly disclosed (correct citation). All authors have been personally and actively involved in substantial work leading to the manuscript, and will take public responsibility for its content.
Author contributions
MZ: conceptualization, data curation, formal analysis, funding acquisition, investigation, methodology, resources, software, validation, visualization, writing—original draft. HA: writing—formal analysis and writing—revising draft. PG: formal analysis and writing—revising draft. NS: conceptualization, data curation, formal analysis, investigation, methodology, resources, supervision, validation, visualization, writing—original draft, and writing—review and editing. GS: conceptualization, data curation, formal analysis, funding acquisition, investigation, methodology, project administration, resources, software, supervision, validation, visualization, writing—original draft, and writing—review and editing.
Acknowledgments
The authors are thankful to Bam University of Medical Sciences.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AsadiP.GiviM. K. B.RastgooA.AkbariM.ZakeriV.RasouliS. (2012). Predicting the Grain Size and Hardness of AZ91/SiC Nanocomposite by Artificial Neural Networks. Int. J. Adv. Manuf Technol.63, 1095–1107. 10.1007/s00170-012-3972-z
2
AshtianiH. R. R.ShahsavariP. (2016). A Comparative Study on the Phenomenological and Artificial Neural Network Models to Predict Hot Deformation Behavior of AlCuMgPb alloy. J. Alloys Compd.687, 263–273. 10.1016/j.jallcom.2016.04.300
3
AtesH. (2007). Prediction of Gas Metal Arc Welding Parameters Based on Artificial Neural Networks. Mater. Des.28, 2015–2023. 10.1016/j.matdes.2006.06.013
4
BeneaL.WengerF.PonthiauxP.CelisJ. P. (2009). Tribocorrosion Behaviour of Ni-SiC Nano-Structured Composite Coatings Obtained by Electrodeposition. Wear266, 398–405. 10.1016/j.wear.2008.04.018
5
CanakciA.VarolT.OzsahinS. (2015). Artificial Neural Network to Predict the Effect of Heat Treatment, Reinforcement Size, and Volume Fraction on AlCuMg alloy Matrix Composite Properties Fabricated by Stir Casting Method. Int. J. Adv. Manuf Technol.78, 305–317. 10.1007/s00170-014-6646-1
6
CoşkunM. İ.Karahanİ. H. (2018). Modeling Corrosion Performance of the Hydroxyapatite Coated CoCrMo Biomaterial Alloys. J. Alloys Compd.745, 840–848. 10.1016/j.jallcom.2018.02.253
7
DarbandG. B.AfsharA.AliabadiA. (2016). Zn-Ni Electrophosphating on Galvanized Steel Using Cathodic and Anodic Electrochemical Methods. Surf. Coat. Technol.306, 497–505. 10.1016/j.surfcoat.2015.12.089
8
DarbandG. B.AfsharA.RabaniM. (2016). Effect of Treatment Time and Temperature on Microstructure and Corrosion Behavior of Zn-Ni Electrophosphate Coating. J. Alloys Compd.688, 596–604. 10.1016/j.jallcom.2016.07.032
9
DashtianK.Zare-DorabeiR. (2017). Synthesis and Characterization of Functionalized Mesoprous SBA-15 Decorated with Fe3O4 Nanoparticles for Removal of Ce(III) Ions from Aqueous Solution: ICP-OES Detection and central Composite Design Optimization. J. Colloid Interf. Sci.494, 114–123. 10.1016/j.jcis.2017.01.072
10
DattaS.ChattopadhyayP. P. (2013). Soft Computing Techniques in Advancement of Structural Metals. Int. Mater. Rev.58, 475–504. 10.1179/1743280413y.0000000021
11
EvisZ.ArcakliogluE. (2011). Artificial Neural Network Investigation of Hardness and Fracture Toughness of Hydroxylapatite. Ceramics Int.37, 1147–1152. 10.1016/j.ceramint.2010.10.037
12
FerreiraC. (2006). Gene Expression Programming: Mathematical Modeling by an Artificial Intelligence. Berlin/Heidelberg, Germany: Springer.
13
FinkeA.EscobarJ.MunozJ.PetitM.TechnologyC. (2021). Prediction of Salt spray Test Results of Micro Arc Oxidation Coatings on AA2024 Alloys by Combination of Accelerated Electrochemical Test and Artificial Neural Network. Surf. Coat. Technol.421, 127370. 10.1016/j.surfcoat.2021.127370
14
GandomiA. H.RokeD. A. (2015). Assessment of Artificial Neural Network and Genetic Programming as Predictive Tools. Adv. Eng. Softw.88, 63–72. 10.1016/j.advengsoft.2015.05.007
15
GolafshaniE. M.RahaiA.SebtM. H. (2015). Artificial Neural Network and Genetic Programming for Predicting the Bond Strength of GFRP Bars in concrete. Mater. Struct.48, 1581–1602. 10.1617/s11527-014-0256-0
16
HoseinianF. S.FaradonbehR. S.AbdollahzadehA.RezaiB.Soltani-MohammadiS. (2017). Semi-autogenous Mill Power Model Development Using Gene Expression Programming. Powder Technol.308, 61–69. 10.1016/j.powtec.2016.11.045
17
HosseinzadehF.SarpoolakiH.HashemiH. (2014). Precursor Selection for Sol-Gel Synthesis of Titanium Carbide Nanopowders by a New Intuitionistic Fuzzy Multi-Attribute Group Decision-Making Model. Int. J. Appl. Ceram. Technol.11, 681–698. 10.1111/ijac.12108
18
HouT.-H.SuC.-H.LiuW.-L. (2007). Parameters Optimization of a Nano-Particle Wet Milling Process Using the Taguchi Method, Response Surface Method and Genetic Algorithm. Powder Technol.173, 153–162. 10.1016/j.powtec.2006.11.019
19
HuQ.LiuY.ZhangT.GengS.WangF. (2019). Modeling the Corrosion Behavior of Ni-Cr-Mo-V High Strength Steel in the Simulated Deep Sea Environments Using Design of experiment and Artificial Neural Network. J. Mater. Sci. Technol.35, 168–175. 10.1016/j.jmst.2018.06.017
20
HuangY.AnC.ZhangQ.ZangL.ShaoH.LiuY.et al (2021). Cost-effective Mechanochemical Synthesis of Highly Dispersed Supported Transition Metal Catalysts for Hydrogen Storage. Nano Energy80, 105535. 10.1016/j.nanoen.2020.105535
21
IsleemH. F.TayehB. A.AlaloulW. S.MusaratM. A.RazaA. (2021). Artificial Neural Network (ANN) and Finite Element (FEM) Models for GFRP-Reinforced Concrete Columns under Axial Compression. Materials14, 7172. 10.3390/ma14237172
22
JafariM. M.SoroushianS.KhayatiG. R. (2017). Hardness Optimization for Al6061-MWCNT Nanocomposite Prepared by Mechanical Alloying Using Artificial Neural Networks and Genetic Algorithm. J. Ultrafine Grained Nanostructured Mater.50, 23–32. 10.7508/JUFGNSM.2017.01.04
23
JafariS.MahiniS. S. (2017). Lightweight concrete Design Using Gene Expression Programing. Construction Building Mater.139, 93–100. 10.1016/j.conbuildmat.2017.01.120
24
JamshidiM.GhaediM.DashtianK.HajatiS.BazrafshanA. A. (2016). Sonochemical Assisted Hydrothermal Synthesis of ZnO: Cr Nanoparticles Loaded Activated Carbon for Simultaneous Ultrasound-Assisted Adsorption of Ternary Toxic Organic Dye: Derivative Spectrophotometric, Optimization, Kinetic and Isotherm Study. Ultrason. Sonochem.32, 119–131. 10.1016/j.ultsonch.2016.03.004
25
Jiménez-ComeM. J.MartínM. D. L. L.MatresV.BaladésJ. D. M. (2020). The Use of Artificial Neural Networks for Modelling Pitting Corrosion Behaviour of EN 1.4404 Stainless Steel in marine Environment: Data Analysis and New Developments. Corrosion Rev.38. 10.1515/corrrev-2019-0095
26
Karahanİ. H.ÖzdemirR. (2015). A Comparison for Grain Size Calculation of Cu-Zn Alloys with Genetic Programming and Neural Networks. Acta Physica Pol. A128, B427–B431. 10.12693/aphyspola.128.b-427
27
Karahanİ. H.OzdemirR.ErkaymanB. (2013). A Comparison of Genetic Programming and Neural Networks; New Formulations for Electrical Resistivity of Zn-Fe Alloys. Appl. Phys. A.113, 459–476. 10.1007/s00339-013-7544-3
28
KennedyJ.EberhartR. (1995). “Particle Swarm Optimization,” in Proceedings of the IEEE International Conference on Neural Networks, Perth, WA, Australia, 27 Nov.-1 Dec. 1995.
29
KhalajG.NazariA.KhoieS. M. M.KhalajM. J.PouraliakbarH. (2013). Chromium Carbonitride Coating Produced on DIN 1.2210 Steel by Thermo-Reactive Deposition Technique: Thermodynamics, Kinetics and Modeling. Surf. Coat. Technol.225, 1–10. 10.1016/j.surfcoat.2013.02.030
30
KimJ. H.ReddyN. S.YeomJ. T.HongJ. K.LeeC. S.ParkN.-K. (2009). Microstructure Prediction of Two-phase Titanium alloy during Hot Forging Using Artificial Neural Networks and FE Simulation. Met. Mater. Int.15, 427–437. 10.1007/s12540-009-0427-7
31
LaugierS.RichonD. (2003). Use of Artificial Neural Networks for Calculating Derived Thermodynamic Quantities from Volumetric Property Data. Fluid phase equilibria210, 247–255. 10.1016/s0378-3812(03)00172-9
32
LiJ. G.YaoY. X.GaoD.LiuC. Q.YuanZ. J. (2008). Cutting Parameters Optimization by Using Particle Swarm Optimization (PSO). Appl. Mech. Mater.10, 879–883. 10.4028/www.scientific.net/AMM.10-12.879
33
LiY.MacdonaldD. D.YangJ.QiuJ.WangS. (2020). Point Defect Model for the Corrosion of Steels in Supercritical Water: Part I, Film Growth Kinetics. Corrosion Sci.163, 108280. 10.1016/j.corsci.2019.108280
34
LinB.-L.LuJ.-T.KongG. (2008). Effect of Molybdate post-sealing on the Corrosion Resistance of Zinc Phosphate Coatings on Hot-Dip Galvanized Steel. Corrosion Sci.50, 962–967. 10.1016/j.corsci.2007.12.002
35
LiuG.JiaL.KongB.GuanK.ZhangH. (2017). Artificial Neural Network Application to Study Quantitative Relationship between Silicide and Fracture Toughness of Nb-Si Alloys. Mater. Des.129, 210–218. 10.1016/j.matdes.2017.05.027
36
LuZ.ZhangF.-Y.FuH.DingH.ChenL.-Q. (2021). Rotational Nonlinear Double-Beam Energy Harvesting. Smart Mater. Structures31, 210–218. 10.1088/1361-665X/ac4579
37
MansouriI.GholampourA.KisiO.OzbakkalogluT. (2018). Evaluation of Peak and Residual Conditions of Actively Confined concrete Using Neuro-Fuzzy and Neural Computing Techniques. Neural Comput. Applic29, 873–888. 10.1007/s00521-016-2492-4
38
MansouriI.KisiO. (2015). Prediction of Debonding Strength for Masonry Elements Retrofitted with FRP Composites Using Neuro Fuzzy and Neural Network Approaches. Composites B: Eng.70, 247–255. 10.1016/j.compositesb.2014.11.023
39
MartíP.ShiriJ.Duran-RosM.ArbatG.De CartagenaF. R.Puig-BarguésJ. (2013). Artificial Neural Networks vs. Gene Expression Programming for Estimating Outlet Dissolved Oxygen in Micro-irrigation Sand Filters Fed with Effluents. Comput. Electron. Agric.99, 176–185. 10.1016/j.compag.2013.08.016
40
MozaffariS.LiW.ThompsonC.IvanovS.SeifertS.LeeB.et al (2017). Colloidal Nanoparticle Size Control: Experimental and Kinetic Modeling Investigation of the Ligand-Metal Binding Role in Controlling the Nucleation and Growth Kinetics. Nanoscale9, 13772–13785. 10.1039/c7nr04101b
41
NouiriM.BekrarA.JemaiA.NiarS.AmmariA. C. (2018). An Effective and Distributed Particle Swarm Optimization Algorithm for Flexible Job-Shop Scheduling Problem. J. Intell. Manuf29, 603–615. 10.1007/s10845-015-1039-3
42
OkhovatA.MousaviS. M. (2012). Modeling of Arsenic, Chromium and Cadmium Removal by Nanofiltration Process Using Genetic Programming. Appl. Soft Comput.12, 793–799. 10.1016/j.asoc.2011.10.012
43
RashidiA. M.HayatiM.RezaeiA. (2012). Application of Artificial Neural Network for Prediction of the Oxidation Behavior of Aluminized Nano-Crystalline Nickel. Mater. Des.42, 308–316. 10.1016/j.matdes.2012.06.011
44
SuH.-Y.LinC.-S. (2014). Effect of Additives on the Properties of Phosphate Conversion Coating on Electrogalvanized Steel Sheet. Corrosion Sci.83, 137–146. 10.1016/j.corsci.2014.02.002
45
TsaiC.-Y.LiuJ.-S.ChenP.-L.LinC.-S. (2010). A Two-step Roll Coating Phosphate/molybdate Passivation Treatment for Hot-Dip Galvanized Steel Sheet. Corrosion Sci.52, 3385–3393. 10.1016/j.corsci.2010.06.020
46
VarolT.CanakciA.OzsahinS. (2013). Artificial Neural Network Modeling to Effect of Reinforcement Properties on the Physical and Mechanical Properties of Al2024-B4c Composites Produced by Powder Metallurgy. Composites Part B: Eng.54, 224–233. 10.1016/j.compositesb.2013.05.015
47
VettivelS. C.SelvakumarN.LeemaN. (2013). Experimental and Prediction of Sintered Cu-W Composite by Using Artificial Neural Networks. Mater. Des.45, 323–335. 10.1016/j.matdes.2012.08.056
48
WuX.LiuZ.YinL.ZhengW.SongL.TianJ.et al (2021). A Haze Prediction Model in Chengdu Based on LSTM. Atmosphere12, 1479. 10.3390/atmos12111479
49
XuC.RangaiahG. P.ZhaoX. S. (2015). Application of Artificial Neural Network and Genetic Programming in Modeling and Optimization of Ultraviolet Water Disinfection Reactors. Chem. Eng. Commun.202, 1415–1424. 10.1080/00986445.2014.952813
50
XuY.-Y.LinB.-L. (2007). Effect of Silicate Pretreatment, post-sealing and Additives on Corrosion Resistance of Phosphated Galvanized Steel. Trans. Nonferrous Met. Soc. China17, 1248–1253. 10.1016/s1003-6326(07)60257-x
51
YusupN.ZainA. M.HashimS. Z. M. (2012). Overview of PSO for Optimizing Process Parameters of Machining. Proced. Eng.29, 914–923. 10.1016/j.proeng.2012.01.064
52
ZadshakoyanM.PourmostaghimiV. (2013). Genetic Equation for the Prediction of Tool-Chip Contact Length in Orthogonal Cutting. Eng. Appl. Artif. Intelligence26, 1725–1730. 10.1016/j.engappai.2012.10.016
53
ZakiM. R.VarshosazJ.FathiM. (2015). Preparation of agar Nanospheres: Comparison of Response Surface and Artificial Neural Network Modeling by a Genetic Algorithm Approach. Carbohydr. Polym.122, 314–320. 10.1016/j.carbpol.2014.12.031
54
ZeraatiM.KhayatiM. (2018). Optimization of Micro Hardness of Nanostructure Cu-Cr-Zr Alloys Prepared by the Mechanical Alloying Using Artificial Neural Networks and Genetic Algorithm. J. Ultrafine Grained Nanostructured Mater.51, 183–192. 10.22059/JUFGNSM.2018.02.11
55
ZhangL.WangL.ZhangY.WangD.GuoJ.ZhangM.et al (2021). The Performance of Electrode Ultrafiltration Membrane Bioreactor in Treating Cosmetics Wastewater and its Anti-fouling Properties. Environ. Res.206, 112629. 10.1016/j.envres.2021.112629
56
ZhangQ.ZhouC.XiaoW.NelsonP. C. (2007). “Improving Gene Expression Programming Performance by Using Differential Evolution,” in Proceedings of the Sixth International Conference on Machine Learning and Applications (ICMLA 2007), Cincinnati, OH, USA, December 2007 (Piscataway, New Jersey, United States: IEEE), 31–37. 10.1109/icmla.2007.62
57
ZhangX.TangY.ZhangF.LeeC.-S. (2016). A Novel Aluminum-Graphite Dual-Ion Battery. Adv. Energ. Mater.6, 1502588. 10.1002/aenm.201502588
58
ZhuH.ZhuJ.ZhangZ.ZhaoR. (2021). Crossover from Linear Chains to a Honeycomb Network for the Nucleation of Hexagonal Boron Nitride Grown on the Ni (111) Surface. The J. Phys. Chem. C48, 26542–26551. 10.1021/acs.jpcc.1c09334
Summary
Keywords
artificial neural network, Zn–Ni electrophosphate coating, corrosion current, modeling, genetic algorithm, gene expression programming
Citation
Zeraati M, Abbasi H, Ghaffarzadeh P, Chauhan NPS and Sargazi G (2022) Application of Artificial Neural Networks for Corrosion Behavior of Ni–Zn Electrophosphate Coating on Galvanized Steel and Gene Expression Programming Models. Front. Mater. 9:823155. doi: 10.3389/fmats.2022.823155
Received
26 November 2021
Accepted
04 February 2022
Published
21 March 2022
Volume
9 - 2022
Edited by
Norhazilan Md Noor, University of Technology Malaysia, Malaysia
Reviewed by
Sudagar J, VIT-AP University, India
Indra Bhushan Singh, Council of Scientific and Industrial Research (CSIR), India
Updates

Check for updates
Copyright
© 2022 Zeraati, Abbasi, Ghaffarzadeh, Chauhan and Sargazi.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ghasem Sargazi, g.sargazi@gmail.com; Narendra Pal Singh Chauhan, narendrapalsingh14@gmail.com
This article was submitted to Computational Materials Science, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.