Next Article in Journal
A Study on a Cast Steel Reinforced with WC–Metal Matrix Composite
Previous Article in Journal
Durability of Nano-Reinforced Recycled Aggregate Concrete under Load and Chloride Ingress
Previous Article in Special Issue
Analysis of Cutting Stability of a Composite Variable-Section Boring Bar with a Large Length-to-Diameter Ratio Considering Internal Damping
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Training Deep Neural Networks with Novel Metaheuristic Algorithms for Fatigue Crack Growth Prediction in Aluminum Aircraft Alloys

by
Muhammad Hamza Zafar
1,†,
Hassaan Bin Younis
2,†,
Majad Mansoor
3,†,
Syed Kumayl Raza Moosavi
2,
Noman Mujeeb Khan
1 and
Naureen Akhtar
4,*
1
Department of Electrical, Capital University of Science and Technology, Islamabad 44000, Pakistan
2
School of Electrical and Electronics Engineering, National University of Sciences and Technology, Islamabad 44000, Pakistan
3
Department of Automation, University of Science and Technology of China, Hefei 230027, China
4
Department of Engineering Sciences, University of Agder, 4879 Grimstad, Norway
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 18 June 2022 / Revised: 7 August 2022 / Accepted: 29 August 2022 / Published: 6 September 2022
(This article belongs to the Special Issue Manufacturing and Fatigue Properties of Materials)

Abstract

:
Fatigue cracks are a major defect in metal alloys, and specifically, their study poses defect evaluation challenges in aluminum aircraft alloys. Existing inline inspection tools exhibit measurement uncertainties. The physical-based methods for crack growth prediction utilize stress analysis models and the crack growth model governed by Paris’ law. These models, when utilized for long-term crack growth prediction, yield sub-optimum solutions and pose several technical limitations to the prediction problems. The metaheuristic optimization algorithms in this study have been conducted in accordance with neural networks to accurately forecast the crack growth rates in aluminum alloys. Through experimental data, the performance of the hybrid metaheuristic optimization–neural networks has been tested. A dynamic Levy flight function has been incorporated with a chimp optimization algorithm to accurately train the deep neural network. The performance of the proposed predictive model has been tested using 7055 T7511 and 6013 T651 alloys against four competing techniques. Results show the proposed predictive model achieves lower correlation error, least relative error, mean absolute error, and root mean square error values while shortening the run time by 11.28%. It is evident through experimental study and statistical analysis that the crack length and growth rates are predicted with high fidelity and very high resolution.

1. Introduction

There is a significant importance of fracture mechanics in mechanical engineering. Most structural failures in automobiles, aircraft, bridges, and tanks are a result of fatigue crack propagation in materials [1]. Structural failure of a 58-year-old Grumman G73T Turbo Mallard Seaplane is one of the examples that occurred in 2005; the right hand wing of which shed on a domestic flight in the USA. In 2002, a Chinese Boeing 777 crashed mid-air as a result of a structural failure as well. The crash of a Japanese Boeing 747 SR100 was also due to structural failure. Such accidents pose a serious threat to automobile manufacturing companies to survive, as these failures lead towards human loss, as well as enormous economic deficit. This remains to be a matter of serious concern. Crack propagation under fatigue has been comprehensively explored by researchers and professionals in this field and many of them tried to predict it using various methodologies with the help of theoretical and experimental knowledge [2,3].
In engineering, fatigue is a process of material failure at a stage that is far below the actual material strength when subjected to cyclic loading [4]. The traditional models attempt to predict the fatigue crack growth (FCG) rate d a d N with respect to the stress intensity factor (SIF), Δ K [5]. A typical graph is shown in Figure 1. In this figure, a represents the length of the crack, while N represents the number of cycles that leads to failure of structure. The FCG curve is divided into three regions. The crack growth rate is much slower in Region I. In Region II, the graph shows a linear and stable growth followed by a rapid growth in Region III [6]. Paris and Erdogan [7] presented a pioneer model in their study using an extensive pool of data, which is also known as Paris’ law. It presented the accurate modeling of an FCG curve, especially in Region II, which is also called the Paris region. Forman et al. [8] claimed that the stress ratio effect was ignored in the Paris model. The model was improved in their study by considering the stress ratio and crack instability effect. It is stated that this model is more flexible compared to the classical Paris model as it is based upon a large amount of data. Priddle et al. [9] presented a method that improves Paris’ law. However, this model failed to address SIF at threshold, K t h , average SIF, and stress ratio effects [10]. Kujawski et al. [11] designed a parameter to relate stress ratio and FCG rate, both for long and short cracks in the Paris region. All the above mentioned models cannot handle non-linearity in Region II and III of the graph [12]. A short literature summary has been provided in Table 1.
In comparison to the extensive literature on empirical techniques, machine learning-based models have lesser literature available as this is a relatively modern technique [19]. Jang et al. [20] forecasted the fatigue life of chains used in cranes to lift heavy duty cargo loads using machine learning. Specifically, a logistic regression-based model was used. The experimental data were extracted using strain gauges. Fatigue life prediction in overhead wires was predicted using artificial neural networks (ANNs) by Nowel and Nowel [21]. Yan et al. [22] analyzed the failure of a bridge as a result of an overloading of vehicles using machine learning methodology. This approach gave more accurate results as compared to conventional finite element analysis-based techniques. In [23], the researchers used a radial basis function neural network (RBF-NN) to predict the fatigue crack growth rate in the airplane industry. The problem of vibrations in electrical transmission lines caused by wind, earthquakes, or snow loading was addressed by Pestana et al. using ANNs [24]. Younis [25] carried his work forward by predicting fatigue crack length using radial-based function-NN (RBF-NN). Moreover, optimized neural networks were used to predict the FCG rate in aircraft aluminum structures by Younis et al. [12]. The hill climbing optimized neural network produced better results among the other two by having better relevance with experimental data.

2. Fracture Mechanics

Fracture mechanics is the study of crack propagation in materials. All engineering structures are prone to geometrical discontinuities. The size and shapes of such features are crucial as they pose serious concerns about the strength of the component [26]. Structures exposed to cyclic loading fail far below their strength due to fatigue. It refers to the progressive degradation of material due to cyclic loading when subjected to applied force such as stress, strain, and torque [27]. Fatigue crack growth is discussed in context to different crack regions as shown in Figure 2 for 2324 T39 and 6013 T651 materials. Crack initiation is the first phase of crack propagation, followed by stable crack propagation and unstable crack propagation, respectively. In the first region, small cracks initiate that are of the order 10–25 µm or even less. This region is also called the threshold region. In Region II or the Paris region, the crack growth rate is of the order 10 6 10 3   mm / cycle . After crossing Region II, the crack propagates exponentially; that is, of the order 10 3   mm / cycle . The Paris region is a main concern of our research, which is explained by Paris and Erdogan [7] in the form of the crack growth equation. It is also called Paris’ law. The stress intensity factor K illustrates load all over the crack tip. The Δ K d can be utilized to correlate the R-ratio effects [27]. Crack propagation is the function of SIF range and Δ K during cycling loading. It is mathematically presented as (1):
d a d N = C Δ K m
where a is the crack length, d a d N is the crack growth rate for load cycle N . C and m are constants that are material dependent. SIF range can be represented as (2):
Δ K = K m a x K m i n

3. Proposed Technique

3.1. Dynamic Levy Flight Chimp Optimization Algorithm

In the following section, a new chimp optimization algorithm with Levy Flight random movement (DLFCO) is put forward in order to train DNN. The presented procedure updates the weights and biases of the DNN in order to decrease the cost mapping. DLFCO is an optimization algorithm based on the populace that has the application in the hunting and preying behavior of the chimps to update the weights and biases with the Levy flight random walk in order to obtain global maxima in the exploitation phase. Fusion of LF with the ChoA improves the exploration and exploitation capability of the algorithm [28]. This extremely exploratory attitude enables it to locate the overall best answer, that in turn efficiently trains the DNN model. This trained framework forecasts the fatigue crack growth for the aluminum alloys in aircraft applications. The preceding segment elucidates the chimp optimization and sine–cosine merged chimp optimization algorithm.

3.1.1. Chimp Optimization Algorithm (ChoA)

This model obtains its motivation by the intellect of the chimpanzee and breeding actions in cluster hunting. Four techniques are implemented to replicate the attitude, which are attack, chase, barrier, and driver as show in Figure 3. The scientific equations of the driver and attacker of the presented approach are mentioned in (3) and (4), respectively.
D = C · α p r e y m · α c h i m p  
α c h i m p n + 1 = α p r e y a · d
where m, c, and d are the co-efficient vectors. These factors may be rationalized by (5)–(7).
a = 2 · l · r 1 l
c = 2 · r 2
m = c h a o t i c   v a l u e
where r 1 and r 2 represent random numbers in the range of 0–1 and l is a constant which decreases in the direction of a line from 2.5 till 0 all along the iterations and m is a disordered vector. To mathematically simulate this system, the four best results with optimal capability are chosen which are attacker, barrier, chaser, and driver. The remaining inhabitants will upgrade their location with help from the data provided by the best four results. Their statistical demonstration is given by Equations (8), (9), (10), and (11), respectively.
d a t t a c k = C 1 · α a t t a c k e r m 1 · x n
d b a r r i e r = C 2 · α b a r r i e r m 2 · x n
d c h a s e r = C 3 · α c h a s e r m 3 · x n
d d r i v e r = C 4 · α d r i v e r m 4 · x n
The succeeding point of the chimps are restructured consequently using (12)–(15):
x 1 = α a t t a c k e r a 1 · d a t t a c k e r
x 2 = α b a r r i e r a 2 · d b a r r i e r
x 3 = α c h a s e r a 3 · d c h a s e r
x 4 = α d r i v e r a 4 · d d r i v e r
Using the above equations, the positions are upgraded through:
x 1 = x 1 + x 2 + x 3 + x 4 4
Subsequently, after upgrading the positions, (17) is applied.
α c h i m p n + 1 = α p r e y x · d , < 0.5 c h a o t i c , < 0.5

3.1.2. Levy Flight (LF)

Levy flight is the random walk used for position updating in a search of global maximum position in the area of search. LF is a non-Gaussian distribution based on arbitrary numbers. The Levy flight random walk is presented in Figure 4. The Gaussian method is a relatively stable procedure, containing summation of the Gaussian variables and Gaussian spread. The probability density function is shown in (18).
L x = 1 π 0 e γ p α cos p x d p
The Levy spread has two constraints, i.e., α , γ , and is proportional as well, regarding x = 0 . Where γ is the scaling aspect, γ > 0 and α is between [0, 2] [29]. The factor α defines the spreading figure to attain a number of probability distribution characters; specifically, the end of dispersion is governed by the factor α . Figure 4 provides the chaotic map of the Levy flight showing long and short jumps in the search space.

3.1.3. Improved ChoA with LF

The established ChoA keeps its agents posted in the direction of its prey, i.e., the best position on the basis of the driver, chaser, barrier, and attacker positions. On the other hand, the search neurons of the algorithm can become stuck in local minima in some circumstances. So, the premature conjunction issue may pose a threat. In a few circumstances, the acknowledged ChoA is unable to perform the transition from the exploration to exploitation phase and back smoothly. In this section, LF is applied to speak about the above stated limitations. The DLFCO offers more deep searching arrays that lead towards an efficient global search. This arrangement can solve immobility issues as well. In addition, the agent features must be upgraded in DLFCO through all iterations. The chimp’s location can be updated using (19) and (20).
x c h i m p = x 1 + x 2 + x 3 + x 4 4 + r a n d 0   L e v y W 1   if   | a |   <   0.5
  x c h i m p = x 1 + x 2 + x 3 + x 4 4 + r a n d 0   L e v y W 2   if   | a |   >   0.5
where ⊕ designates multiplication in the hierarchical entry, W1 and W2 are decreasing over the iterations, in that order; these equations are defined in a manner that W1 shows a descending angle and W2 shows an ascending angle, and both of the parameters show an arbitrary performance. The pseudo code of the proposed technique is presented in Figure 5.

3.2. Deep Neural Network (DNN)

An artificial neural network (ANN) is a brain-like structure in which learning is undertaken by experience. It consists of three layers, i.e., input layer, hidden layer, and output layer. The input layer acts as a channel to pass information to the hidden layer by multiplying with weights W i j , which are weights from the input layer to the hidden layer. The relationship of input and output data of the hidden layer is presented in Equation (22). The hidden layer applies to the predefined mathematical model which is called the activation function. The output of the hidden layer is passed onto the output layer after multiplying it with weights W j k . The relationship between the output and hidden layer is presented by Equation (23).
Figure 5. Pseudocode of DLFCO.
Figure 5. Pseudocode of DLFCO.
Materials 15 06198 g005
The number of weights and biases are directly proportional to the number of neurons in the hidden layer. The selection of the number of neurons in the hidden layer depend upon the non-linearity of the dataset. A high number of neurons could result in overfitting, while a lower number of neurons could lead to inefficient training of the ANN. Therefore, in order to deal with a highly non-linear dataset, multiple hidden layers need to be added between the input layer and output layer, which is called the deep neural network.
The deep neural network (DNN) is a multi-layered structure like the ANN, but it consists of more than one hidden layer. In this work, a layered DNN is implemented which contains one input layer, two hidden layer neurons, and one output layer, as shown in Figure 6. Six neurons are selected in the hidden layer.

3.3. Activation Function

For active mapping of inputs and outputs, a function is required in the hidden layer neurons called the activation function. Each hidden layer neuron has the activation function which is chosen on the basis of application. There are different types of activation functions in classification and regression applications. For classification, the sigmoid function is commonly used as in Equation (21).
a ˘ i = 1 1 + e x i     S i g m o i d   F u n c t i o n
As the fatigue crack growth rate is a regression problem and has a highly non-linear nature, the radial basis function neural network is chosen as the activation function, which is given by (22).
v j y j = exp i = 1 r y i j y ^   i j 2 μ R B F
Figure 6. Structure of DNN.
Figure 6. Structure of DNN.
Materials 15 06198 g006

3.4. Cost Function

The cost function or fitness function is the function that tells the difference between actual data and predicted data. This is the function which needs to be minimized during training of the ANN by updating the weights and biases. Various types of cost functions are presented in Equations (23)–(26).
MAE = i = 0 n y y n     M e a n   A b s o l u t e   E r r o r
MSE = i = 0 n y y 2 n     M e a n   S q u a r e d   E r r o r
RMSE = i = 0 n y y 2 n     R o o t   M e a n   S q u a r e d   E r r o r
NRMSE = RMSE X o b s , m a x X o b s , m i n     N o r m a l i z e d   R o o t   M e a n   S q u a r e d   E r r o r
The normalized root mean square error (NRMSE) is the cost function which is selected for this problem and is presented in Equation (26). The specifications of the DNN used in this network are presented in Table 2.

3.5. Neural Network Training Using DLFCO

Updating weights and biases in the neural network is called the training of the neural network. The back propagation (BP) algorithm is one of the simplest algorithms to train neural networks, but it takes lots of iterations for effective training, which makes BPNN expensive computationally. For effective training of weights and biases, the metaheuristic optimization algorithm is used in this work, which calculates the cost function and updates weights and biases according to the mathematical model of the optimization algorithm.
As shown in Figure 7, firstly, 50 sets of weights and biases are initialized over the search space randomly. These sets of weights and biases are the particle’s position. Then, the DNN structure with two hidden layers having six neurons in each layer are initialized. After that, fitness is calculated for every particle using NMSE, which defines the best particle’s position for the current iteration in the population. After calculation of the best particle, the particle’s position is updated using the optimization algorithm. This phase of optimization is called the exploration phase.

4. Results

This section first discusses the data collection for the fatigue crack growth rate, evaluation parameters for the prediction results, and prediction of the growth rate.

4.1. Data Collection

The data extraction has been undertaken for dissimilar aircraft structures of aluminum with the range of R-ratio from 0.1 to 0.7 after comprehensive experimentation of alloy specimens in the Fracture Technology Associates Laboratory, USA. The experimentation setup includes servo-controlled, hydraulically actuated, closed-loop assembly that is furthermore linked to a standalone PC for information gathering and governing of the factors. Test conditions which are undertaken for experimentation are summarized in Table 3.

4.2. Evaluation Parameter

In order to assess diverse techniques, further error directories are applied as well. The grade of spreading in outcomes may be corroborated through normalized root mean square error (NRMSE), given by Equation (27). To specify deviancy of extrapolation, mean absolute error (MAE) and mean absolute percentage error (MAPE) is given by Equations (28) and (29), correspondingly. Furthermore, the relationship between real and forecasted value can be calculated with the help of R-square ( R 2 ), which is mentioned in Equation (30).
NRMSE = 1 T 1 N i = 1 N T i P i 2 × 100 %
MAPE = 1 N i = 1 N T i P i P i × 100 %
NMAE = 1 N i = 1 N T i P i P i
R 2 = i = 1 N T i T i ¯ ) · P i P i ¯ 2 i = 1 N T i T i ¯ · i = 1 N P i P i ¯
where T i represents the actual value, P i is the forecasted value, N is the overall sum of illustrations. T i ¯ is the mean of the actual output value, P i ¯ is the mean of forecasted output value.

4.3. Prediction Results

The results are presented in the graphical, statistical, and qualitative measures. For aircraft alloys, namely, 7055 T7511, 6013 T651, and 2324 T39, prediction results and accuracy of prediction with respect to the sample are compared in Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 The results show the rate increase as the independent parameters progress non-linearly. The corresponding algorithm performances and statistical analysis is presented in Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9.

5. Discussion

To verify the prediction performance of shape and growth rate evolution, the testing set has been employed. The aluminum alloys under experimental data are 6013 T651, 7055 T7511, and 2324 T39. The 6013 T651 is medium strength aerospace alloy that offers enhanced formidability and corrosion resistance. The 2324 T39 alloy exhibits higher strength with better fracture toughness. The 7055 T7511 alloy is predominantly employed in compression-dominated structures due to improvised compressive and tensile strengths. The experiments conducted at the Fracture Technology Associates Laboratory utilize aluminum with an R-ratio range of 0.1–0.7 in several experiments. The stress was applied using a servo-controlled, hydraulically actuated, closed-loop assembly equipped with HIL data acquisition and parameter control. Strain gauge transducers measure the displacement that is utilized to measure the length and growth rates secondarily [30].
The performance of the proposed framework has been compared with three highly effective and recently studied optimization algorithms. It includes particle swarm optimization algorithms that provide standard behavior of metaheuristic optimization algorithms, which have the capability to distinguish among sub-optimum solutions. The second algorithm is the grey wolf optimization algorithm (GWO). The improved information sharing mechanism among GWO solutions enhances the accuracy towards optimum solutions. The comparison is made in terms of testing error, training errors, and algorithm run time. A more comprehensive statistical analysis is made for MFA-DNN, GWO-DNN, PSO-DNN, and DLFCO-DNN in Table 8 for statistical indices such as MAE, RMSE, correlation ( R 2 ), and mean RE [31].
Figure 9, Figure 11 and Figure 13 exhibit the trained DLFCO-DNN estimate RE of predicted values at various samples with high redundancy. The lack of accuracy in competing techniques is primarily associated with over- and/or under-fitting of machine learning algorithms. The stagnation and model fitting issues arise when the predictive values show high variance and low bias due to the noise present in the data. It is avoidable through accurate selection of parameters and preprocessing of data outliers and feature mapping during the pre-processing of the data. Table 4, Table 5, Table 6, Table 7, Table 8 and Table 9 show a summary of the results for d a / d N rate prediction using DLFCO-DNN on three aluminum alloys. The testing and training values are compared alongside run time. The best run time has been achieved by a simple criterion that defines the boundary of the small-scale yielding regime, to avoid invalid use of the LEFM parameter, ΔK, as the characterizing parameter for the fatigue crack growth rate [32].
While in uniform amplitude loading, the FCG graph—crack growth rate d a / d N against SIF (ΔK) in log–log measure—normally comprises three portions. Region I denotes the earlier crack propagation. In this region, the crack normally grows in the order of 10 6   mm / cycle span. This section of the curve is subjective to the micro-structural appearances, the R-ratio (lowest useful stress divided by maximum useful stress), plus the ecological circumstances. In general, it is believed that there subsists a SIF at threshold Δ K t h . Under this point, fatigue crack propagation does not take place [33]. Speedy and unsteady fatigue crack growth takes place in Region III of the curve before fracture. A crack progression rate of the order of   10 3   mm / cycle is ordinary in the case of metals and is asymptotic when crack growth reaches the threshold of fracture toughness K c of the material [34]. This is even though it is mostly dominated by the stress ratio effect, substance micro-cracks, and the specimen width. Moreover, crack growth in the third region of the curve is usually overlooked as a result of trivial leftover fatigue life while arriving in the region. Stable crack growth is observed in Region II of the FCG curve, typically in the range of 10 6     mm / cycle   to   10 3     mm / cycle in the case of metals [35].
It is pertinent to mention that in the case of a hole in the crack route, the crack growth rate in Region I reduces considerably, while in Region II, it elevates [36]. It is noticeable in the trajectory of crack propagation; in a straight line crack progression, KI escalates, but upon the change in direction, KII escalates [37]. The preceding graphs of SIF reveal that all simulated specimen testing data exhibit almost similar numbers for KI and KII throughout crack progression intervals. These insignificant deviations in the final values are dominated by the quantity of phases employed in every package and the mathematical methodology applied. The quantity of phases applied in ANSYS software for model simulation range between 11 and 20, with the total amount of percentage increase organized with the help of software. In FRANC2D/L, the steps are between 32 and 92 with a 0.1 to 0.13 increment per step.
Figure 8, Figure 10 and Figure 12 show comparisons of MFA-DNN, GWO-DNN, PSO-DNN, and DLFCO-DNN fatigue life simulations with experimental outcomes for 6013 T651, 7055 T7511, and 2324 T39. The simulated fatigue crack growth life using the proposed method comes to an impeccable agreement, while the experimental results are provided by Gomes and Miranda [38]. Moreover, the present study outcomes came out to be much more accurate for prediction of fatigue life as linked to the arithmetical results achieved by means of the software. The outcomes of such a study are carefully compared with Quera2D instead of BemCracker2D. Such results, moreover, establish the impact of the opening location regarding the fatigue life of the samples in various areas of impact with Δ K stress intensity range on a log scale [39]. Care must be taken while taking into account the use of Δ K , which is usable for high strength alloys as well [32]. The deterioration caused by the operational, environmental, and climate conditions quantitation should be improvised [40].

6. Conclusions

This study proposes a new framework for the fatigue crack growth rate in aluminum alloys usually utilized in aircraft manufacturing, specifically 7055 T7511, 2324 T39, and 6013 T651. A comprehensive modeling and analysis was undertaken to accurately predict the fatigue crack growth in aluminum aircraft alloys. A hybrid DLFCO-DNN framework has been proposed to utilize experimental data for training the predictive models. Existing models based on physical models lack robustness and accuracy demonstrated by the proposed model. The quantitative and statistical analysis made using statistical indices re-affirms the superior performance of DLFCO-DNN in terms of training, testing, and run time of the algorithms. Moreover, the accuracy in prediction as compared to recent works is up to 4–8 times higher; meanwhile, the DLFCO-DNN algorithm complexity allows for hardware implementation on less costly controllers such as c2000 for real-time functionality. In light of improvised results, it is safe to conclude that the proposed novel framework hybrid fatigue crack growth rate prediction is highly effective.

Author Contributions

Conceptualization, M.H.Z., H.B.Y. and M.M.; Data curation, M.H.Z., H.B.Y. and M.M.; Formal analysis, S.K.R.M.; Funding acquisition, N.A.; Investigation, S.K.R.M.; Methodology, M.H.Z., H.B.Y. and M.M.; Project administration, N.A.; Software, M.H.Z., H.B.Y., and M.M.; Supervision, N.A.; Validation, N.M.K.; Visualization, N.M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Drury, C.G. Human Factors in Aircraft Maintenance; State University of New York at Buffalo Dept of Industrial Engineering: New York, NY, USA, 2001. [Google Scholar]
  2. Zhao, K.; Xue, H.; Yang, F.; Zhao, L. Probability Prediction of Crack Growth Rate of Environmentally Assisted Cracks of Nickel-Based Alloys Based on Latin Hypercube Sampling. Int. J. Press. Vessel. Pip. 2019, 172, 391–396. [Google Scholar] [CrossRef]
  3. Mohanty, J.R.; Verma, B.B.; Ray, P.K.; Parhi, D.R. Prediction of mode-I overload-induced fatigue crack growth rates using neuro-fuzzy approach. Expert Syst. Appl. 2010, 37, 3075–3087. [Google Scholar] [CrossRef]
  4. Cheng, Z.Q.; Xiong, J.J.; Tan, W. Fatigue crack growth and life prediction of 7075-T62 aluminium-alloy thin-sheets with low-velocity impact damage under block spectrum loading. Int. J. Fatigue 2022, 155, 106618. [Google Scholar] [CrossRef]
  5. Yarullin, R.; Shlyannikov, V.; Amato, D.; Citarella, R. Mixed mode surface crack growth in aluminium alloys under complex stress state. Procedia Struct. Integr. 2022, 39, 364–378. [Google Scholar] [CrossRef]
  6. Jones, R.; Kovarik, O.; Cizek, J.; Ang, A.; Lang, J. Crack growth in conventionally manufactured pure nickel, titanium and aluminum and the cold spray additively manufactured equivalents. Addit. Manuf. Lett. 2022, 3, 100043. [Google Scholar] [CrossRef]
  7. Paris, P.; Erdogan, F. A critical analysis of crack propagation laws. J. Basic Eng. 1963, 85, 528–533. [Google Scholar] [CrossRef]
  8. Forman, R.G.; Kearney, V.; Engle, R. Numerical analysis of crack propagation in cyclic-loaded structures. J. Basic Eng. 1967, 89, 459–463. [Google Scholar] [CrossRef]
  9. Priddle, E. High cycle fatigue crack propagation under random and constant amplitude loadings. Int. J. Pressure Vessels Piping 1976, 4, 89–118. [Google Scholar] [CrossRef]
  10. Liu, C.; Gong, Y.; Gong, Y.; Li, W.; Liu, Z.; Hu, N. Mode II fatigue delamination behaviour of composite multidirectional laminates and the stress ratio effect. Eng. Fract. Mech. 2022, 264, 108321. [Google Scholar] [CrossRef]
  11. Kujawski, D. A new (ΔK+ Kmax) 0.5 driving force parameter for crack growth in aluminum alloys. Int. J. Fatigue 2001, 23, 733–740. [Google Scholar] [CrossRef]
  12. Younis, H.B.; Kamal, K.; Sheikh, M.F.; Hamza, A. Prediction of fatigue crack growth rate in aircraft aluminum alloys using optimized neural networks. Theor. Appl. Fract. Mech. 2022, 117, 103196. [Google Scholar] [CrossRef]
  13. Hu, D.; Su, X.; Liu, X.; Mao, J.; Shan, X.; Wang, R. Bayesian-based probabilistic fatigue crack growth evaluation combined with machine-learning-assisted GPR. Eng. Fract. Mech. 2020, 229, 106933. [Google Scholar] [CrossRef]
  14. Ma, X.; He, X.; Tu, Z. Prediction of fatigue–crack growth with neural network-based increment learning scheme. Eng. Fract. Mech. 2021, 241, 107402. [Google Scholar] [CrossRef]
  15. Mortazavi, S.; Ince, A. An artificial neural network modeling approach for short and long fatigue crack propagation. Comput. Mater. Sci. 2020, 185, 109962. [Google Scholar] [CrossRef]
  16. Wang, H.; Zhang, W.; Sun, F.; Zhang, W. A comparison study of machine learning based algorithms for fatigue crack growth calculation. Materials 2017, 10, 543. [Google Scholar] [CrossRef]
  17. Rovinelli, A.; Sangid, M.D.; Proudhon, H.; Guilhem, Y.; Lebensohn, R.A.; Ludwig, W. Predicting the 3D fatigue crack growth rate of small cracks using multimodal data via Bayesian networks: In-situ experiments and crystal plasticity simulations. J. Mech. Phys. Solids 2018, 115, 208–229. [Google Scholar] [CrossRef]
  18. Nejad, R.M.; Shariati, M.; Farhangdoost, K. Prediction of fatigue crack propagation and fractography of rail steel. Theor. Appl. Fract. Mech. 2019, 101, 320–331. [Google Scholar] [CrossRef]
  19. Shi, J.; Wang, J.; Macdonald, D.D. Prediction of crack growth rate in Type 304 stainless steel using artificial neural networks and the coupled environment fracture model. Corros. Sci. 2014, 89, 69–80. [Google Scholar] [CrossRef]
  20. Jang, D.W.; Lee, S.; Park, J.W.; Baek, D.C. Failure detection technique under random fatigue loading by machine learning and dual sensing on symmetric structure. Int. J. Fatigue 2018, 114, 57–64. [Google Scholar] [CrossRef]
  21. Nowell, D.; Nowell, P. A machine learning approach to the prediction of fretting fatigue life. Tribol. Int. 2020, 141, 105913. [Google Scholar] [CrossRef]
  22. Yan, W.; Deng, L.; Zhang, F.; Li, T.; Li, S. Probabilistic machine learning approach to bridge fatigue failure analysis due to vehicular overloading. Eng. Struct. 2019, 193, 91–99. [Google Scholar] [CrossRef]
  23. Younis, H.B.; Kamal, K.; Sheikh, M.F.; Hamza, A.; Zafar, T. Prediction of fatigue crack growth rate in aircraft aluminum alloys using radial basis function neural network. In Proceedings of the 2018 Tenth International Conference on Advanced Computational Intelligence (ICACI), Xiamen, China, 29–31 March 2018; IEEE: New York, NY, USA, 2018. [Google Scholar]
  24. Reinke, G.; Badibanga, R.K.; Pestana, M.S.; de Almeida Ferreira, J.L.; Araujo, J.A.; da Silva, C.R. Failure analysis of aluminum wires in all aluminum alloy conductors-AAAC. Eng. Fail. Anal. 2020, 107, 104197. [Google Scholar] [CrossRef]
  25. Younis, H.B.; Kamal, K.; Sheikh, M.F.; Iqbal, S. Prediction of fatigue crack length in aircraft aluminum alloys using radial basis function neural network. In Proceedings of the 2021 International Bhurban Conference on Applied Sciences and Technologies (IBCAST), Islamabad, Pakistan, 12–16 January 2021; IEEE: New York, NY, USA, 2021. [Google Scholar]
  26. Karpenko, O.; Oterkus, S.; Oterkus, E. Titanium alloy corrosion fatigue crack growth rates prediction: Peridynamics based numerical approach. Int. J. Fatigue 2022, 162, 107023. [Google Scholar] [CrossRef]
  27. Kujawski, D. Correlating R-ratio Effects on FCG Behavior Using ΔKd Function. Theor. Appl. Fract. Mech. 2022, 118, 103244. [Google Scholar] [CrossRef]
  28. Zafar, M.H.; Khan, N.M.; Mansoor, M.; Khan, U.A. Towards green energy for sustainable development: Machine learning based MPPT approach for thermoelectric generator. J. Clean. Prod. 2022, 351, 131591. [Google Scholar] [CrossRef]
  29. Mansoor, M.; Mirza, A.F.; Duan, S.; Zhu, J.; Yin, B.; Ling, Q. Maximum energy harvesting of centralized thermoelectric power generation systems with non-uniform temperature distribution based on novel equilibrium optimizer. Energy Convers. Manag. 2021, 246, 114694. [Google Scholar] [CrossRef]
  30. Yang, Q.; Yang, X.; Huang, W.; Shi, Y.; Shi, D. Small fatigue crack propagation rate and behaviours in a powder metallurgy superalloy: Role of stress ratio and local microstructure. Int. J. Fatigue 2022, 160, 106861. [Google Scholar] [CrossRef]
  31. Mansoor, M.; Mirza, A.F.; Long, F.; Ling, Q. An Intelligent Tunicate Swarm Algorithm Based MPPT Control Strategy for Multiple Configurations of PV Systems Under Partial Shading Conditions. Adv. Theory Simul. 2021, 4, 2100246. [Google Scholar] [CrossRef]
  32. Marques, B.; Borges, M.F.; Antunes, F.V.; Vasco-Olmo, J.M.; Díaz, F.A.; James, M.N. Limitations of small-scale yielding for fatigue crack growth. Eng. Fract. Mech. 2021, 252, 107806. [Google Scholar] [CrossRef]
  33. Ando, K.; Kim, M.H.; Nam, K.W. Analysis on peculiar fatigue fracture behaviour of shot peened metal using new threshold stress intensity factor range equation. Fatigue Fract. Eng. Mater. Struct. 2021, 44, 306–316. [Google Scholar] [CrossRef]
  34. Nakai, Y.; Kikuchi, S.; Osaki, K.; Kawabata, M.O.; Ameyama, K. Effects of rolling reduction and direction on fatigue crack propagation in commercially pure titanium with harmonic structure. Int. J. Fatigue 2021, 143, 106018. [Google Scholar] [CrossRef]
  35. Zhao, T.; Zhang, J.; Jiang, Y. A study of fatigue crack growth of 7075-T651 aluminum alloy. Int. J. Fatigue 2008, 30, 1169–1180. [Google Scholar] [CrossRef]
  36. Shlyannikov, V.N. Creep–fatigue crack growth rate prediction based on fracture damage zone models. Eng. Fract. Mech. 2019, 214, 449–463. [Google Scholar] [CrossRef]
  37. Henaff, G.; Odemer, G.; Benoit, G.; Koffi, E.; Journet, B. Prediction of creep–fatigue crack growth rates in inert and active environments in an aluminium alloy. Int. J. Fatigue 2009, 31, 1943–1951. [Google Scholar] [CrossRef]
  38. Gomes, G.; Miranda, A.C. Analysis of crack growth problems using the object-oriented program bemcracker2D. Frat. Ed. Integrità Strutt. 2018, 12, 67–85. [Google Scholar] [CrossRef]
  39. Zvyagin, A.V.; Luzhin, A.A.; Smirnov, N.N.; Shamina, A.A.; Shamin, A.Y. Stress intensity factors for branching cracks in space structures. Acta Astronaut. 2021, 180, 66–72. [Google Scholar] [CrossRef]
  40. Tamimi, M.F.; Khandel, O.; Soliman, M. A framework for quantifying fatigue deterioration of ship structures under changing climate conditions. Ships Offshore Struct. 2022, 1–16. [Google Scholar] [CrossRef]
Figure 1. Crack growth rate against stress intensity range.
Figure 1. Crack growth rate against stress intensity range.
Materials 15 06198 g001
Figure 2. (a) Alloy 2324 T39; (b) alloy 6013 T651.
Figure 2. (a) Alloy 2324 T39; (b) alloy 6013 T651.
Materials 15 06198 g002
Figure 3. ChoA position updating.
Figure 3. ChoA position updating.
Materials 15 06198 g003
Figure 4. Levy flight movement.
Figure 4. Levy flight movement.
Materials 15 06198 g004
Figure 7. Flow chart for training of DNN using DLFCO.
Figure 7. Flow chart for training of DNN using DLFCO.
Materials 15 06198 g007
Figure 8. Alloy 7055 T7511 crack growth rate prediction.
Figure 8. Alloy 7055 T7511 crack growth rate prediction.
Materials 15 06198 g008
Figure 9. Alloy 7055 T7511 crack growth rate relative error.
Figure 9. Alloy 7055 T7511 crack growth rate relative error.
Materials 15 06198 g009
Figure 10. Alloy 6013 T651 crack growth rate prediction.
Figure 10. Alloy 6013 T651 crack growth rate prediction.
Materials 15 06198 g010
Figure 11. Alloy 6013 T651 crack growth rate relative error.
Figure 11. Alloy 6013 T651 crack growth rate relative error.
Materials 15 06198 g011
Figure 12. Alloy 2324 T39 crack growth rate prediction.
Figure 12. Alloy 2324 T39 crack growth rate prediction.
Materials 15 06198 g012
Figure 13. Alloy 2324 T39 crack growth rate relative error.
Figure 13. Alloy 2324 T39 crack growth rate relative error.
Materials 15 06198 g013
Table 1. Brief summary of the literature work.
Table 1. Brief summary of the literature work.
SN. ReferenceTechnique(s)Summary
1.Hu. Dianyin et al. [13]Gaussian Process Regression (GPR)Bayesian-based calibration method is presented to improve computational efficacy and accuracy for FCG life prediction. In comparison with experimental data, predicted results are bounded within a factor of ±2.0
2.Ma. Xinran et al. [14]Increment Learning Scheme A fully connected neural network is presented using the Increment Learning Scheme for crack growth under tension on the specimens of 7B04 aluminum and TA15 titanium alloys. The proposed method is superior to conventional fitting techniques. The minimum achieved MSE is 0.10 for aluminum alloy for R = 0.06
3.Mortazavi et al. [15]Radial Basis Function Neural Network (RBF-NN)Radial basis function-based neural network approach is presented to predict FCG behavior. The proposed technique has been verified through experimentation on aluminum and titanium alloys. This model handles non-linearity in a better way for both long and short cracks
4.Wang et al. [16]Extreme Learning Machine, RBF-NN, and GA-BPThree different machine learning algorithms (MLAs) are compared for FCG prediction. The algorithms are tested on datasets of different alloys of aluminum and titanium. MLAs are compared with each other and with the classical K* approach. ELM outperforms its counterparts with minimum MSE, i.e., 3.12 × 10−8
5.Rovinelli et al. [17]Bayesian Networks This study utilizes in situ experimental and crystal plasticity simulations study multimodal data with Bayesian networks to predict the fatigue crack growth rate in 3D for small cracks. A non-local driving force is postulated to construct a data-driven probabilistic crack propagation framework.
6.Nejad et al. [18]Rockwell micro-hardness experimentPrediction of fatigue crack propagation and fractography of rail steel is made using Machine Learning algorithms. The modified Paris model is used to estimate fatigue crack growth rates. Fatigue crack growth and hardness tests are carried out by fractography studies on fractured specimens. A three-dimensional boundary element method is used for fatigue crack growth study under a stress field.
Table 2. Parameter descriptions of model used for prediction.
Table 2. Parameter descriptions of model used for prediction.
ParameterValue
Number of Layers4
Number of Hidden Layers2
Number of Neurons in Hidden Layer6
Number of Weights and Biases59
Number of Iterations50
Activation FunctionRBF
Cost FunctionNRMSE
Training AlgorithmDLFCOA
Table 3. Test conditions of alloy hardware experiments in the lab.
Table 3. Test conditions of alloy hardware experiments in the lab.
Specimen TypeM (T)
Notch Lenght (2 a 0 )10.2 mm
Width (W)101.6 mm
Thickness (B)2.2–3.1 mm (Variable)
Test Frequency15 Hz
d a / d N 1 × 10 12 1 × 10 4 m/cycle
Orientation L-T
Lab Air Temp.24 °C
Relative Humidity (R.H)50–55%
Table 4. Alloy 7055 T7511 algorithm performance results.
Table 4. Alloy 7055 T7511 algorithm performance results.
TechTesting MSETraining MSERun Time (s)
DLFCO-DNN1.33 × 10−105.394 × 10−1198.43 ms
MFA-DNN4.55 × 10−109.103 × 10−11101.45 ms
GWO-DNN8.21 × 10−108.547 × 10−10108.21 ms
PSO-DNN1.67 × 10−92.409 × 10−10112.04 ms
Table 5. Alloy 7055 T7511 statistical analysis results.
Table 5. Alloy 7055 T7511 statistical analysis results.
TechMAERMSER2Mean RE
DLFCO-DNN7.11 × 10−59.34 × 10−80.12130.0512
MFA-DNN2.98 × 10−47.97 × 10−80.19760.9563
GWO-DNN3.09 × 10−52.65 × 10−60.20151.1453
PSO-DNN3.85 × 10−41.01 × 10−70.27681.4300
Table 6. Alloy 6013 T651 algorithm performance results.
Table 6. Alloy 6013 T651 algorithm performance results.
TechTesting MSETraining MSERun Time (s)
DLFCO-DNN1.01 × 10−102.44 × 10−1199.01 ms
MFA-DNN2.44 × 10−108.01 × 10−11102.30 ms
GWO-DNN7.01 × 10−104.40 × 10−10109.22 ms
PSO-DNN4.21 × 10−91.99 × 10−10114.23 ms
Table 7. Alloy 6013 T651 statistical analysis results.
Table 7. Alloy 6013 T651 statistical analysis results.
TechMAERMSER2Mean RE
DLFCO-DNN7.14 × 10−57.33 × 10−63.48 × 10−40.0141
MFA-DNN4.26 × 10−41.21 × 10−44.22 × 10−40.2114
GWO-DNN5.65 × 10−45.85 × 10−55.04 × 10−40.5468
PSO-DNN9.23 × 10−41.99 × 10−45.47 × 10−40.5090
Table 8. Alloy 2324 T39 results.
Table 8. Alloy 2324 T39 results.
TechTesting MSETraining MSERun Time (s)
DLFCO-DNN2.71 × 10−105.04 × 10−1196.01 ms
MFA-DNN3.97 × 10−101.21 × 10−10106.30 ms
GWO-DNN6.05 × 10−95.43 × 10−10112.24 ms
PSO-DNN9.01 × 10−81.99 × 10−9115.03 ms
Table 9. Alloy 2324 39 results.
Table 9. Alloy 2324 39 results.
TechMAERMSER2Mean RE
DLFCO-DNN1.0913 × 10−41.0971 × 10−50.95210.0911
MFA-DNN8.4101 × 10−48.4721 × 10−50.97810.4023
GWO-DNN1.7400 × 10−41.6021 × 10−50.99010.2581
PSO-DNN8.6701 × 10−48.4431 × 10−50.99210.4221
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zafar, M.H.; Younis, H.B.; Mansoor, M.; Moosavi, S.K.R.; Khan, N.M.; Akhtar, N. Training Deep Neural Networks with Novel Metaheuristic Algorithms for Fatigue Crack Growth Prediction in Aluminum Aircraft Alloys. Materials 2022, 15, 6198. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15186198

AMA Style

Zafar MH, Younis HB, Mansoor M, Moosavi SKR, Khan NM, Akhtar N. Training Deep Neural Networks with Novel Metaheuristic Algorithms for Fatigue Crack Growth Prediction in Aluminum Aircraft Alloys. Materials. 2022; 15(18):6198. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15186198

Chicago/Turabian Style

Zafar, Muhammad Hamza, Hassaan Bin Younis, Majad Mansoor, Syed Kumayl Raza Moosavi, Noman Mujeeb Khan, and Naureen Akhtar. 2022. "Training Deep Neural Networks with Novel Metaheuristic Algorithms for Fatigue Crack Growth Prediction in Aluminum Aircraft Alloys" Materials 15, no. 18: 6198. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15186198

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

Article Metrics

Back to TopTop