Research on a Lamb Wave and Particle Filter-Based On-Line Crack Propagation Prognosis Method

Prognostics and health management techniques have drawn widespread attention due to their ability to facilitate maintenance activities based on need. On-line prognosis of fatigue crack propagation can offer information for optimizing operation and maintenance strategies in real-time. This paper proposes a Lamb wave-particle filter (LW-PF)-based method for on-line prognosis of fatigue crack propagation which takes advantages of the possibility of on-line monitoring to evaluate the actual crack length and uses a particle filter to deal with the crack evolution and monitoring uncertainties. The piezoelectric transducers (PZTs)-based active Lamb wave method is adopted for on-line crack monitoring. The state space model relating to crack propagation is established by the data-driven and finite element methods. Fatigue experiments performed on hole-edge crack specimens have validated the advantages of the proposed method.


Introduction
Reliability and availability are key problems for safety-critical systems such as aircraft, wind turbines, bridges and nuclear plants [1,2]. However, conventional maintenance frameworks may involve longer downtime and greater cost without considering the actual status of the systems [3]. In recent years, prognostics and health management (PHM) techniques, which consider the actual system condition via diagnostic techniques and the future condition through prognosis methods, have drawn widespread attention due to their ability to enable maintenance activities based on need [4][5][6]. Accurate prognosis for the degradation state and failure time of the critical structure plays an important role in the PHM technique, which leads to an increase of reliability and availability. Moreover, the life cycle cost will be reduced by undertaking maintenance activities only as necessary and minimizing downtime and spare part storage.
Fatigue cracks are commonly regarded as a principal failure mode for various structural and mechanical systems [7,8]. In recent years, a lot of attention has been paid to the methods which combine fatigue crack propagation models with Bayes' theorem for fatigue crack propagation prognosis, including stochastic filters [9][10][11][12][13][14][15] and the Bayesian inference [16][17][18]. Within these methods, the uncertainties during fatigue crack propagation are taken into account. The measurement information of the actual crack propagation state is used to update the result obtained by the crack propagation model to achieve a more accurate one. As to stochastic filters, the Kalman filter (KF) [19] offers the optimal solution to linear problems under Gaussian uncertainty assumption. Nevertheless, most realistic cases such as the process of fatigue crack propagation are nonlinear with non-Gaussian uncertainty. To tackle these cases, various types of approximation are developed for the KF [20]. The extended Kalman filter (EKF) utilizes the first-order Taylor expansion of the nonlinear model to solve the fatigue crack propagation problem [9]. The unscented Kalman filter (UKF) is proposed to deal with the fatigue crack propagation problem using unscented transformation [10]. However, the Gaussian uncertainty assumption is still needed for the EKF and UKF. As an improvement of the KF, the particle filter (PF) is capable of handling the prognosis problem of nonlinear and non-Gaussian processes without restrictive assumptions based on Monte Carlo methods [21]. In recent years, several PF based methods have been reported for fatigue crack propagation prognosis. Shin et al. [11] adopted the PF to deal with the prognosis problem of fatigue crack propagation. The visual inspections performed by an optical magnifier were used as the measurements of the crack lengths. Corbetta et al. [12] proposed a kind of stochastic dynamic state space model for the PF based method, which utilized the crack lengths obtained with a caliper as the measurements. Compare and Zio [13] proposed a PF-based method for predictive maintenance, in which simulated crack lengths were employed for validation. Sun et al. [14] analyzed the sources of uncertainties in fatigue crack propagation prognosis using Virkler's data which was measured by a zoom stereomicroscope and explored a PF-based algorithm for uncertainty management. All these studies indicate that the PF has the potential for fatigue crack propagation prognosis under uncertainties.
However, most literatures published so far have only used simulation or off-line NDT results, which have limitations for on-line application. On-line crack monitoring is capable of offering convenient and quick inspections of crack damages with sensors. Timely detection and prognosis of crack damages can maximize the operational availability and safety by optimizing operation and maintenance strategies in real-time. Recently, on-line crack monitoring methods have been gradually combined with PF-based methods to realize on-line crack propagation prognosis. For example, Chen et al. [15] proposed a PF-based method for the machine condition prediction, in which the vibration feature extraction method is applied for crack monitoring. However, vibration feature extraction-based methods are insensitive to small damage or damage growth [22]. In general, research on prognosis methods integrating on-line crack monitoring with the PF is still lacking, as well as experimental verifications. With the development of the structural health monitoring (SHM) technology, different kinds of methods have been developed for on-line crack monitoring [23]. Among them, The PZTs-based active Lamb wave (LW) technique is one of the most appealing and effective methods [24][25][26], which has the merits including the ability of traveling a long distance, the capacity to access hidden components, as well as sensitivity to small crack damages [27,28].
Aiming at realizing on-line fatigue crack propagation prognosis, a LW-PF-based method is proposed to combine the PZTs-based active Lamb wave method with the PF. The PZT sensor array is used to actuate and acquire Lamb wave signals in the structure. The cross-correlation damage index extracted from the monitored Lamb wave signal is employed to capture signal characteristics and quantify the actual crack length. Each time a new damage index is available, the PF utilizes this damage index to estimate the posterior probability density function (pdf) of the crack length with the crack propagation state space model. This state space model is derived from a stochastic Paris's law and an active Lamb wave-based measurement equation, whose parameters and uncertainty are determined by data driven methods. On the basis of the obtained posterior pdf, the prognosis of the crack propagation is performed and the failure cycle is calculated afterwards. The proposed method is evaluated with the fatigue experiments of 6 hole-edge crack specimens, and the posterior estimation and prognostic value of the crack length are discussed, as well as the failure cycle.
The rest of the paper is organized as follows: in Section 2, the proposed on-line LW-PF-based crack propagation prognosis method is explained. First, the state equation of the fatigue crack propagation is presented. Then, the PZTs-based active Lamb wave method is introduced. The modeling process of the measurement equation is proposed with the cross-correlation damage indices extracted from experimental Lamb wave signals. At the end, the LW-PF based prognosis method is presented in detail. In Section 3, the proposed method is implemented and validated on 6 hole-edge crack specimens. Finally, the conclusions are given in Section 4.

On-Line LW-PF-Based Crack Propagation Prognosis Method
To apply the PF-based method, essentially a state space model describing fatigue crack propagation is needed, which consists of a state equation and a measurement equation. The state equation describes the evolution of the crack length, and the measurement equation governs the relationship between the damage index and the crack length.

Physical Model Based State Equation with Data Driven Parameters
In this section, the method of establishing the state equation is proposed, in which the crack propagation mechanism is based on Paris' law with the parameters obtained by data driven methods.

Stochastic Paris' Law Based State Equation
The general state equation of fatigue crack propagation can be expressed by Equation (1): where k is the discrete time index, x k is the crack length at time k, f p¨q is a nonlinear function representing the evolution of the crack length from time k´1 to time k, ω k´1 is a random variable denoting the uncertainty during fatigue crack propagation. Paris' law [29] has been widely used to describe fatigue crack propagation, giving the crack propagation rate as Equation (2): where N is the number of loading cycles, C and m are material constants, ∆K is the stress intensity factor (SIF) range determined by the crack length and the fatigue load.
To take the uncertainty during fatigue crack propagation into consideration, the crack propagation rate can be modified as Equation (3) according to Yang and Manning's model [30]: where X(N) is a stationary random process denoting the uncertainty during fatigue crack propagation. For simplification, X(N) can be reduced to a lognormal random variable X [31] and expressed as exppωq, where ω follows the normal distribution Np0, σ 2 ω q. As a result, Equation (3) can be discretized as the state Equation (4): where ∆N is the discrete step of loading cycles.

Parameters of the State Equation
The parameters ∆K, C, m, and σ 2 ω in Equation (4) should be determined in advance for prognosis. Since different kinds of structures have different parameters, the parameters of a certain kind of structure need to be decided specifically. A data driven method is proposed to obtain the parameters as the process shown in Figure 1, where the variable S denotes the number of the specimens.

Computation of the ∆K with the Finite Element Method
First, the finite element method (FEM) is adopted to compute the SIF range ∆K. Given the structure, the crack with specific crack length and position can be simulated in the finite element model by introducing discontinuities between elements. The ∆K is then calculated through J-integral after applying the boundary conditions and the load range ∆P " P max´Pmin . To improve the numerical results, singular elements and fine mesh are applied around the crack tip [32]. By repeating analyses for cracks with different lengths assuming the crack position is constant, a set of SIF ranges can be acquired. A third-order polynomial approximation is then applied for the relationship between the crack length and ∆K as expressed in Equation (5): where ta n , n " 0, 1, 2, 3u are the polynomial coefficients, x is the crack length.

Calculating Material Constants
Fatigue experiments are performed for a set of this kind of structures to determine the material constants C and m. The logarithmic transformation of Equation (2) is expressed as Equation (6): where ln(C) and m are the intercept and slope considering the linear relationship between lnpdx{dNq and lnp∆Kq.
During the experiments, S specimens labeled as T j , j " 1, . . . , S ( are tested and measurements of the crack lengths tx j i , i " 1, ..., Mu are performed M times for each specimen. Their corresponding loading cycles are recorded as tN j i , i " 1, ..., Mu. It should be noted that x j 0 is the crack length first measured, and N j 0 is the corresponding loading cycle. The crack propagation rate at crack length x j i of the specimen T j can be calculated approximately using the experimental data as shown in Equation (7): Taking these crack propagation rates tpdx{dNq j i , i " 1, ..., Mu of each specimen with corresponding SIF ranges t∆Kpx j i q, i " 1, ...., Mu, the two parameters lnpC j q and m j of the specimen T j are calculated using linear regression [31]. In this paper, the mean value of lnpC j q and the mean value of m j are adopted as the material constants in the state equation.

Estimation of the Crack Propagation Uncertainty
The logarithm transformation of Equation (4) is shown as Equation (8). Empirically, the difference between the logarithm of the crack propagation rate obtained from the experimental observations and the one from Paris' law is used to approximate the crack propagation uncertainty as illustrated in Equation (9): where the material constants C and m are determined previously, the difference value e j i is supposed to be the crack propagation uncertainty that disturbs the crack propagation rate. Therefore, the variance σ 2 ω of the random variable ω is calculated as the variance of all the difference values from all the specimens as shown in Equation (10): where Var denotes the variance of the difference values. After these processes, the state equation is obtained as Equation (11):

On-Line Crack Monitoring by the Active Lamb Wave Method
The PZTs-based active Lamb wave method is adopted for on-line crack monitoring. Lamb waves are elastic waves that propagate in plate-like structures. A typical configuration used for crack monitoring is the pitch-catch way as shown in Figure 2. The PZT sensor array is arranged on the structure. After the Lamb wave is excited in the structure by a PZT, interaction of the Lamb wave with the crack can influence the Lamb wave propagation. By comparing the received signals under healthy and cracked conditions from other PZTs, the crack can be estimated.
To evaluate the variations between the baseline signal collected when the structure is healthy and the monitored signal when the crack propagates in the structure, many damage indices have been proposed [22]. In this paper, the cross-correlation damage index is adopted to quantify the crack length as expressed in Equation (12): where H(t) is the baseline signal, D(t) is the monitored signal during the crack propagation process, t 0 and t 1 are the sampling start time and stop time. This damage index reflects the phase changes between the signals caused by the increasing crack length.

Active Lamb Wave-Based Measurement Equation
The measurement equation is defined as Equation (13): where the function gp¨q represents the relationship between the damage index y k and the crack length x k , ν is a random variable denoting the measurement uncertainty. The proposed active Lamb wave-based measurement equation is established by the data driven method as shown in Figure 3. During the fatigue experiment of each specimen T j , j " 1, ..., S, the baseline signal is acquired from the initial state. With the crack propagation in the specimen T j , damage indices are obtained at M specific crack lengths tx i , i " 1, ..., Mu, which are denoted as ty j i , i " 1, ..., Mu. The mean value of the damage indices at the crack length x i of all the specimens can be calculated as Equation (14): After that, a third order polynomial as shown in Equation (15) is employed to govern the relationship between the crack length and the mean value of the damage indices: where tb n , n " 0, 1, 2, 3u are the polynomial coefficients. The measurement uncertainty is assumed to be normally distributed, which denotes ν " Np0, σ 2 ν q. With the damage indices obtained from all the specimens, the variance pσ 2 v q i of the damage indices at crack length x i can be approximated as Equation (16): where Var is the variance of the damage indices. As a result, the maximum one in the set pσ 2 v q i , i " 1, ..., M ( is adopted as σ 2 v of the measurement uncertainty to illustrate the capability the proposed method. Hence, the measurement equation is obtained as Equation (17):

On-Line LW-PF-Based Crack Propagation Prognosis Method
Assuming the current time is k, the objective of the prognosis method is to predict the crack length x k`d , where k + d is the future time of interest. In this paper, the PF is employed to incorporate the damage index of the actual crack state for a more precise prognostic result, which gives the conditional pdf ppx k`d |y 1:k q with the on-line monitored damage indices y 1:k " y j , j " 1, ..., k ( . This conditional pdf represents the probability distribution of x k`d and contains all the available measurement information to the current time k. Then the estimation of the crack lengthx k`d can be expressed as Equation (18) It should be noted that the latest damage index in ppx k`d |y 1:k q can be obtained at time k is y k . There are two procedures in the proposed LW-PF method: the first one is to estimate the posterior pdf ppx k |y 1:k q of the current crack length x k according to the on-line monitored damage indices y 1:k . The second procedure is to predict the future crack length x k`d on the basis of the posterior pdf ppx k |y 1:k q.

Estimation of the Posterior Pdf with On-Line Monitored Damage Indices
From a Bayesian perspective, the posterior pdf ppx k |y 1:k q of the crack length x k is recursively calculated with on-line monitored damage indices in two stages: prediction and update.
Suppose the posterior pdf ppx k´1 |y 1:k´1 q at time k´1 is obtained. The prior pdf ppx k |y 1:k´1 q can be calculated by Equation (19) in the prediction step, where ppx k |x k´1 q is the transition pdf defined by Equation (11). Once the damage index y k is obtained, it is used to calculate the posterior pdf ppx k |y 1:k q through Bayes' theorem as expressed in Equation (20): Update : ppx k |y 1:k q " ppy k |x k qppx k |y 1:k´1 q ppy k |y 1:k´1 q where ppy k |x k q is the likelihood function defined by the measurement equation, ppy k |y 1:k´1 q is a normalized constant [21]. If the posterior pdf ppx k |y 1:k q is obtained, the posterior estimation of the crack length x k can be calculated as Equation (21) x k ppx k |y 1:k qdx k However, Equations (19) and (20) do not have analytical solutions in the most cases. The PF makes the evaluation feasible by resorting to Monte Carlo methods. The basic idea of the PF is to approximate the posterior pdf ppx k |y 1:k q by means of a set of particles tx piq k , i " 1, ..., N s u with their normalized weights t r w piq k , i " 1, ..., N s u: where δ is the Dirac delta function, N s is the number of particles, tx piq k , i " 1, ..., N s u are particles sampled from the pdf ppx k |y 1:k q. Since the posterior pdf is difficult to sampled from, instead, these particles are sampled from a known and easily sampled pdf qpx k |y 1:k q called the importance density function [21]. Hence, the posterior estimation of x k can be expressed as Equation (23), and the corresponding weight w piq k is defined as Equation (24): x k ppx k |y 1:k q qpx k |y 1:k q qpx k |y 1:k qdx k (23) If the importance density function is chosen to factorize such that: qpx k |y 1:k q " qpx k´1 |y 1:k´1 qqpx k |x k´1 , y 1:k q the non-normalized weight w piq k can be derived in a iterative form as shown in Equation (26) [20], and the weight normalization is given as Equation (27): where, ppy k |x piq k q is the likelihood value corresponding to the particle x piq k given the damage index y k . Due to the normal distribution assumption of the measurement uncertainty, this likelihood value is calculated as Equation (28): A common choice of the importance density function is the transition pdf, i.e., qpx piq k |y 1:k q " ppx piq k |x piq k´1 q. In this sense, the calculation of the weights can be reduced to Equation (29): One of the main problems of this PF is the degeneracy phenomenon. After a few iterations, all but one particle will have negligible weight. The effective sample size N e f f expressed in Equation (30) is introduced to evaluate the degradation degree [21]. To deal with this problem, a resampling procedure is performed when N e f f is less than a fixed threshold N th . The resampling eliminates particles with small weights, copying which have large weights, and setting all the weights to 1{N s : As a result, the flow chart to estimate the posterior pdf of the crack length is illustrated in Figure 4.

On-Line Prognosis of Fatigue Crack Propagation
After the posterior pdf ppx k |y 1:k q is obtained, the prognosis procedure is performed. During this procedure, weight updates are no longer taken since no new damage indices can be collected. The condition pdf ppx k`d |y 1:k q is calculated with the posterior pdf ppx k |y 1:k q under the hypothesis of the Markov processes of order one as shown in Equation (31): Replacing ppx k |y 1:k q as its approximation with the particles and corresponding weights gives Equation (32): These integrals can be evaluated by extending the particle x piq k using the state equation as shown in Equation (33). The conditional probability ppx k`d |y 1:k q is approximated by the particles tx piq k`d , i " 1, ..., N s u and the unchanged weights t r w piq k , i " 1, ..., N s u as expressed in Equation (34): After that, the prognostic result of the crack lengthx k`d can be evaluated as Equation (35):

Estimation of the Failure Cycle for the Structure
As to the fatigue crack propagation problem, the failure time of the structure may be represented by the failure cycle N f , which indicates the loading cycle when the cracked structure becomes unusable. Usually, a threshold crack length x th is defined to evaluate whether this structure is still usable. Once the crack length exceeds this threshold, the cracked structure is regarded as failed. As mentioned above, each particle with its weight represents a possible crack propagation path. The failure cycle is estimated by collecting the loading cycle when each path reaches the threshold x th . The pdf of the failure cycle can be expressed as Equation (36): where N piq f is the failure cycle of the i th crack propagation path. Then the failure cycle is evaluated as Equation (37): As illustrated in Figure 5, the flow chart shows the process to obtain the crack propagation prognosis and the estimation of the failure cycle based on the obtained posterior pdf.

Experimental Setup
The experiments of 6 hole-edge crack specimens were performed to verify the effectiveness of the proposed method. The specimens were made of 2 mm thick YL12 aluminum plate and labeled from T 1 to T 6 as shown in Figure 6a.  Table 1 lists their mechanical properties. For each specimen, a 3 mm notch was machined at the edge of the through hole to initiate the crack and control the crack propagation direction during the fatigue test. Two PZTs were attached to the specimen and used as the actuator and the sensor respectively, whose positions are illustrated in Figure 6b. The experimental setup is illustrated in Figure 7. The material test system MTS810 was used for applying the fatigue load. The multi-channel PZT array scanning system developed by the authors' group [33] was employed to perform the active Lamb wave based monitoring. A 5-cycle tone-burst signal with the center frequency of 290 kHz and˘10 V amplitude was used as the excitation signal, which is expressed in Equation (38). Lamb wave signals were sampled at 10 MHz: where N c = 5, f c is the central frequency, A is the amplitude. In this validation research, a tensile test was conducted at first for the specimen T 6 . The ultimate tensile load of the specimen T 6 was obtained as 45 kN. Referring to this result, a sinusoidal load with peak value P max = 15 kN was chosen for the fatigue experiments of the specimens T 1 to T 5 , with the frequency of 10 Hz and the load ratio R = P min /P max = 0.1.

Fatigue Crack Propagation and Damage Indices
The signal under the initial state of each specimen was collected as the baseline signal. Figure 8 illustrates a set of typical Lamb wave signals obtained as the fatigue crack propagation of the specimen T 1 . It is observed that the amplitudes and phases of the signals are influenced by the increasing crack lengths. As illustrated in Figure 9, the fatigue crack propagation results of the specimens T 1 -T 5 are recorded at 2 mm crack length intervals. Obvious uncertainties can be seen in the results. The corresponding damage indices are presented in Figure 10. It can be observed that the damage indices of each specimen are capable of representing the crack propagation process, and likewise measurement uncertainty exists.

On-Line Crack Propagation Estimation and Prognosis of the Specimen T 5
In this validation research, the specimen T 5 was used as the target structure and deemed as unknown. For the on-line crack propagation prognosis of the specimen T 5 , the information of the crack propagation of this kind of structure should be extracted to establish the state space model. Thus, the fatigue experiments of the specimens T 1 -T 4 were conducted ahead to establish the state space model of the specimen T 5 .
3.3.1. Establishing the State Space Model for the Specimen T 5 As the process described in Section 2, the expression of ∆K was obtained as Equation (39): ∆K " 0.016x 3´0 .728x 2`1 9.322x`466.826 The material parameters rlnpC j q, m j s, j " 1, ..., 4 ( of the specimens T 1 -T 4 were calculated as shown in Table 2, where the subscript j represents the specimen T j . The mean values ln(C) =´32, and m = 3.897 were chosen as the material constants C and m in the state equation. The variance of the crack propagation uncertainty was approximated as σ 2 ω " 0.166 2 with the method mentioned above. The state equation is expressed as Equation (40): To obtain the measurement equation, the mean values of the damage indices were calculated by Equation (14) excluding the singular point. As shown in Figure 11, the 3rd order polynomial fit is applied for the mean values, which gives Equation (41): y "´2.962ˆ10´5¨x 3`0 .00156x 2´0 .00656x`0.00798 (41) Figure 11. The relationship of damage index versus the crack length.
The variance of the measurement uncertainty was approximated from the experimental data as expressed in Equation (16). The maximum value pσ 2 ν q max " 0.075 2 was adopted in this paper. As a result, the measurement equation was obtained as Equation (42): where the variable ν subjects to normal distribution Np0, 0.075 2 q.

On-Line Crack Propagation Estimation and Prognosis
The particle number and the resampling threshold of the PF algorithm was set to N s " 2000, N th " 0.8N s respectively. The loading cycle step was chosen as ∆N = 50 cycles. The initial crack length is 3 mm. The N s random samples from the distribution Np3, 0.6 2 q were adopted as the initial particles to obtain a diverse particle set. When the crack length exceeds the predefined threshold x th " 31 mm, the structure is considered to be in a very critical state where the crack tip is close to the structure boundary, and is defined as structural failure.
The specimen T 5 was regarded as unknown and monitored by the active Lamb wave method. The damage indices extracted from the specimen T 5 were used as the measurements of the actual crack lengths to estimate the posterior pdf of the crack length. The damage index was sequentially collected on-line. Once a new damage index was obtained, it was used for updating the weight. Figure 12 illustrates the posterior estimation and prognosis procedures with the proposed LW-PF based method. After the posterior pdf of the crack length is obtained using the available damage indices, the prognosis procedure is conducted based on the obtained posterior pdf. Each particle with the corresponding weight represents a possible crack propagation path. Taking all the possible paths, the prognostic values of the crack length in the future can be calculated, as well as the failure cycle. As shown in Figure 13, the posterior estimations of the crack length considering all the damage indices are closer to the experiment results comparing with the results calculated directly by Paris' law. This is because the LW-PF-based method takes advantage of the on-line monitored information with the crack propagation model. Table 3 gives the crack lengths estimated at N = 7848, N = 17,281, N = 23,286 cycles. As to the experimental results, it can be found that the errors of the LW-PF-based method are smaller than those obtained by Paris' law.  Moreover, the crack propagation paths predicted at loading cycles N = 7848, N = 17,281, N = 23,286 are shown in Figure 14. From these figures, it can be observed that the prognostic values of the crack lengths gradually approach the experimental results. Considering the prognostic value of the crack length at the future loading cycles N = 26,620, the results obtained at different cycles are shown in Table 4. The prognostic value of the proposed method becomes closer to the experimental value as more damage indices are available. Furthermore, the standard deviation of the prognostic pdf represented by the particles and corresponding weights decreases, which indicates the capability of uncertainty reduction of the proposed method.  The crack length threshold is defined as x th " 31 mm. Once the crack length exceeds this threshold, the cracked structure is regarded as failed and the estimation of the failure cycle is obtained as expressed in Equation (37). The relative error of the failure cycle is specified as the percentage of the experimental failure cycle as shown in Equation (43): Relative error " Estimated failure cycle-Experimental failure cycle Experimental failure cycleˆ1 00% (43) As illustrated in Table 5, the estimations of the failure cycle at N = 7848, 17,281, 23,286 cycles are obtained, respectively. From the relative error of the estimated failure cycle, it can be found that the results have been improved by using the proposed method. As more damage indices become available, the prognostic estimation becomes more accurate, which shows the advantages of the proposed method.

Cross Validation Performed on the Specimen T 3
For better validating the proposed method, a leave-one-out cross validation approach was employed by picking a different specimen as the target structure. In this validation, the specimen T 3 was picked as the validation specimen instead, while the experimental data of the specimens T 1 , T 2 , T 4 , T 5 was used for calculating the state space model parameters. By the same procedure, the state space model parameters of the specimen T 3 were obtained as shown in Table 6. The parameters of the PF algorithm, including the particle number, the resampling threshold, and the loading cycle step were set as the same values of the specimen T 5 . Figure 15 shows the posterior estimation of the crack propagation path considering all the damage indices. Similarly, the posterior estimations are closer to their experimental results comparing with the results using the Paris' law. The estimations of the LW-PF based method at load cycles N = 7266, N = 16,514, N = 22,588 are shown in Table 7.  The crack propagation paths of the specimen T 3 predicted at loading cycles N = 7266, N = 14,514, N = 22,588 are obtained as shown in Figure 16. In addition, Table 8 presents the prognostic crack length of the specimen T 3 at future loading cycles N = 26,091. On the other hand, the estimations of the failure cycle obtained at loading cycles N = 7266, N = 16,514, N = 22,588 are shown in Table 9.    The validation experiment of the specimen T 3 have also shown the advantage of the proposed method. By integrating the on-line monitored damage indices, the posterior estimations are more accurate than those obtained by Paris' Law. Moreover, the prognostic crack propagation paths are improved as more damage indices are available, as well as the estimations of the failure cycle.

Conclusions
This paper proposes a LW-PF-based method for on-line crack propagation prognosis. In this method, the active Lamb wave-based method is employed for on-line crack monitoring. The cross-correlation damage indexes extracted from Lamb wave signals are used for measuring the actual crack length. To implement the PF algorithm, the parameters of the physical model-based state equation and the measurement equation are proposed to be predetermined through finite element analysis and data driven methods. The PF gives the probabilistic result for the crack propagation prognosis and the failure cycle estimation by integrating the on-line crack length measurements. The validation experiments performed on the hole-edge crack specimens have shown the capability of the proposed method.