Next Article in Journal
Raman and Conductivity Analysis of Graphene for Biomedical Applications
Next Article in Special Issue
Linear and Nonlinear Guided Wave Imaging of Impact Damage in CFRP Using a Probabilistic Approach
Previous Article in Journal
Design of Boiler Welding for Improvement of Lifetime and Cost Control
Previous Article in Special Issue
A Multi-Level Decision Fusion Strategy for Condition Based Maintenance of Composite Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrated Health Monitoring Method for Structural Fatigue Life Evaluation Using Limited Sensor Data

1
School of Reliability and System Engineering, Beihang University, Beijing 100191, China
2
School of Energy and Power Engineering, Beihang University, Beijing 100191, China
3
Siemens Corporation, Corporate Technology, 755 College Rd. E., Princeton, NJ 08540, USA
*
Author to whom correspondence should be addressed.
Materials 2016, 9(11), 894; https://doi.org/10.3390/ma9110894
Submission received: 29 August 2016 / Revised: 31 October 2016 / Accepted: 31 October 2016 / Published: 4 November 2016
(This article belongs to the Special Issue Advances in Structural Health Monitoring for Aerospace Structures)

Abstract

:
A general framework for structural fatigue life evaluation under fatigue cyclic loading using limited sensor data is proposed in this paper. First, limited sensor data are measured from various sensors which are preset on the complex structure. Then the strain data at remote spots are used to obtain the strain responses at critical spots by the strain/stress reconstruction method based on empirical mode decomposition (REMD method). All the computations in this paper are directly performed in the time domain. After the local stress responses at critical spots are determined, fatigue life evaluation can be performed for structural health management and risk assessment. Fatigue life evaluation using the reconstructed stresses from remote strain gauge measurement data is also demonstrated with detailed error analysis. Following this, the proposed methodology is demonstrated using a three-dimensional frame structure and a simplified airfoil structure. Finally, several conclusions and future work are drawn based on the proposed study.

1. Introduction

Structural health monitoring (SHM) is an effective tool to monitor the status of a target system for maintenance purposes [1,2]. Moreover, SHM has drawn a lot of attention in the SHM community due to flexibilities of smart sensors and advances in composite materials [3,4,5,6]. For most mechanical systems subjected to fatigue loads, fatigue crack propagation is one of the major failure mechanisms, and fatigue damage is mainly dependent on the applied stress history under realistic conditions. The practical challenge for fatigue life prediction in such cases is that the stress and fatigue load in the vicinity of the fatigue crack is very difficult to measure. Therefore, one of the critical aspects in fatigue life prediction analysis is the stress information. The knowledge of the stress information of a structure allows for reliable and accurate fatigue life prediction. It is frequently seen that the critical spots for fatigue damage usually have no direct sensor measurements. Using measurements from pre-installed sensors to infer the information about the critical spot—i.e., a fatigue crack-damaged location—without direct sensor measurements may be the only practical option in realistic situations. Therefore, extrapolation to a critical spot from measurements located in a remote spot is a mandatory step for fatigue prognostics and the remaining life prediction.
Recent advances of dynamical response reconstruction methods include frequency-domain methods and direct time-domain methods [7,8,9]. However, very few published works have been found for fatigue life evaluation using limited sensor data. Direct time-domain reconstruction of dynamical responses for acceleration and velocity using the empirical mode decomposition (EMD) method is demonstrated in [10]. For convenience, the reconstruction methods based on the EMD and the finite element modeling (FEM) of the structure is abbreviated to REMD in this paper. One of the objectives of the proposed reconstruction methodology is to provide accurate stress information for fatigue life prediction using sparse and remote strain gauges installed in the structure. The proposed reconstruction methodology can be readily integrated into existing structural health monitoring systems for prognostics.
Stress responses of a location of interest can be obtained whereby the strain responses for the location are reconstructed using the REMD method with remote strain gauge measurements. The reconstructed stress responses form a spectrum loading for the location of interest. In practice, the spectrum loading for a mechanical part is a complex and often stochastic sequence of load amplitudes. The calculations of fatigue life are different between crack-damaged and normal mechanical parts. For a component without fatigue cracks, the calculation of fatigue life based on the S-N curve [11] and Miner’s rule [12] can be adopted. In such cases, the following procedure is widely used:
  • Reduce the stochastic loading spectrum to a set of simple cyclic loadings using the rainflow counting method [13,14].
  • Perform statistical analysis for the results from the rainflow analysis to create a histogram of cyclic stress and form a fatigue damage spectrum.
  • For each stress level in the fatigue damage spectrum, calculate the degree of cumulative damage using the S-N curve (of the material).
  • Combine the damage contributions of each stress level using Miner’s rule.
Additionally, a number of interesting studies of fatigue life evaluation based on strain field measurements have been raised in recent years. A fast and reliable algorithm for automatic simulations of crack propagation in bi-dimensional and non-planar shell FE models is described in [15]. Giglio [16] performed some damage-tolerant behavior verification tests on full-scale panel Al-lithium alloy specimens instrumented with several strain gauges to obtain a strain map. He monitored the propagation of the crack caused by artificial damage until the progressive failure of the panel reached one or more stringers. Luo [17] investigated an energy-based multiscale damage criterion for a biaxial loading case and conducted numerical fatigue analysis which showed that the multiscale fatigue damage model can provide acceptable prediction of failure. A damage tolerance analysis is also implemented using an FEM-based approach combined with a 3D automatic crack propagation algorithm in [18]. Sbarufatti [19] applied sequential Monte-Carlo sampling to estimate the probabilistic residual life of a structural component subjected to fatigue crack propagation, while real-time estimation of crack length is provided through a committee of artificial neural networks, trained with finite element-simulated strain patterns. However, most of these approaches require cycle-counting techniques to transform the direct stress history to cycle history before the fatigue analysis can be performed. In engineering, loading history is usually complex or even far from being cyclic; generally, the conversion from random load history to the cycle load sequence significantly modifies the load history. Therefore, an integrated health monitoring approach for structural fatigue life evaluation with real-time spectrum loading is proposed in this paper. Fatigue crack growth models reported in [20,21,22] which changed the cycle-counting to small time counting can deal with stochastic spectrum loadings and the reconstructed stress responses can be directly used.
The objective of the study is to develop a novel method for fatigue life evaluation using limited measurements from usage monitoring system. The proposed method extends the previous EMD-based time-domain reconstruction method by using the finite element model to derive a strain transformation function in modal coordinates. EMD method with intermittency criteria is employed in the proposed method to decouple the remote strain measurements into modal coordinates, and the extrapolation of the strain and stress responses is made using the strain and stress transformation function. Mode superposition is used to obtain the strain and stress responses of the critical spot for fatigue life prediction. One potential benefit is that the fatigue damage can be obtained concurrently with structural dynamics analysis, which is critical for real-time decision making.
Figure 1 presents the overall fatigue life evaluation method. First, limited sensor data are measured from various sensors which are preset on the complex structure. Then the strain data at remote spots are used to obtain the strain responses at critical spots. All the computations in this paper are directly performed in the time domain. After the local stress responses at critical spots are determined, fatigue life evaluation can be performed for structural health management and risk assessment.
The paper is organized as follows. First, the EMD method for signal decomposition is briefly introduced, and the required signal filtering process using the intermittency criteria is discussed. Next, the transformation equation for strain and stress is derived using the FE model of the structure. Following this, the fatigue life prediction using the reconstructed stress is calculated and compared with theoretical solutions. A three-dimensional frame structure and a simplified airfoil structure are used to demonstrate and validate the proposed method. Finally, conclusions are drawn based on the current study.

2. Integrated Health Monitoring Method Development

The proposed integrated health monitoring methodology is demonstrated in this section. The overall process for the reconstruction involves several steps, as illustrated in [10]. Details of each of the calculation steps are discussed in this section.

2.1. Strain/Stress Reconstruction Method Based on Empirical Mode Decomposition

An essential step of the proposed dynamic response reconstruction method is to obtain the individual modal responses of the measured signals. The EMD method with intermittency criteria is adopted to decompose the measured signals into modal responses which are called intrinsic mode functions (IMFs) in EMD. In this section, the standard sifting process of the EMD method is briefly discussed.
Using those modal responses, the modal responses at inaccessible locations are computed with the mode shapes. Then modal superposition methodology is employed to calculate the dynamic responses in time domain after all the modal responses have been obtained. To ensure that each IMF contains only one frequency component, the sifting process with intermittency criteria is used. This method is very efficient in terms of computational cost, and very suitable for various dynamic response reconstructions based on the different types of sensor measurements. A limited number of measurement locations are required.
By using the EMD method with intermittency criteria discussed in [10] for different natural frequencies, all modal responses can be obtained. These IMFs have several characteristics: (1) Each IMF contains the intrinsic characteristics of the signal; (2) Once an IMF is obtained, the next IMF will not have the same frequency at the same time instant [23,24]; and (3) The first IMF for each IMF series is considered the approximation of modal response. Using the sifting process with intermittency criteria, the original signal expression can be written as Equation (1).
y ( t ) i = 1 m x i ( t ) + i = 1 n m f i ( t ) + r ( t )
where the original signal y ( t ) can be expressed as the summation of n IMFs and a residue term r ( t ) by applying the standard EMD method. x i ( t ) is the m modal response (that is also an IMF) for the ith mode. Terms f i ( t ) ( i = 1 , n m ) are other IMFs but not modal responses.
For a general structure, finite element model (FEM) can be used as the structural model to derive transformation equations. Considering a general FEM to describe the structure under analysis, the system dynamics equation can be expressed as
M X ¨ + C X ˙ + K X = F
where M, K and C are mass, stiffness, and damping matrices, respectively. X is the displacement vector and F is the load vector. For practical structures subject to stochastic excitations, F is unknown and directly solving Equation (2) to obtain the dynamical responses of a sensor from an inaccessible location is not possible. However, the mode superposition method allows for correlating displacement responses of two different degree of freedoms (DOFs) in the modal coordinates through the mode shape matrix. The mode shape matrix can readily be obtained by solving the eigenvalue problem of
[ Φ , λ ] = eig ( [ M - 1 K ] ) F
where Φ and λ are the eigenvectors and eigenvalues, respectively. Φ (Equation (4)) is also referred to as the mode shape matrix. λ corresponds to the natural frequencies of the structure, i.e., λ = ( 2 π f ) 2 where f is the vector of natural frequencies.
Φ = [ φ 11 φ n 1 φ 1 n φ n n ]
The physical meaning of Φ can be interpreted as follows: each column of Φ represents a mode and each component in the column represents the displacement contribution of a DOF in the structure.
The transformation equation for strain responses under ith mode between two elements indexed by e and u can be obtained by using Equation (6) as
B ( e ) Φ i ( e ) B ( u ) Φ i ( u ) = α B ( e ) δ i ( e ) α B ( u ) δ i ( u ) = η i ( e ) η i ( u )
where B ( k ) is the strain-displacement matrix for element k and Φ i is the ith column vector in the mode shape matrix Φ , η i ( k ) is the strain response vector associated with ith mode and kth element under modal coordinates, α is a scalar constant for a given time index t. The result of Equation (5) indicates that if the physical strain responses are measured at one location (represented by element e in the FE model), it can be decomposed to its modal responses, i.e., ε ( e ) i = 1 ... m η i ( e ) , the physical strain responses at a sensor inaccessible location (represented by element u in the FE model), can be reconstructed using the following transformation equation:
ε ( u ) ( t ) i = 1... m [ η i ( e ) ( t ) ( B ( e ) Φ i ( e ) B ( u ) Φ i ( u ) ) 1 ]
where i = 1 ... m denotes the participating modes and all other notations are defined as before. Once the physical strain responses are reconstructed using Equation (6), the stress responses can be readily calculated using the following constitutive equation:
σ ( t ) = c ε ( t )
where c is the material matrix.
Once the strain and stress responses of locations of interest are obtained using the proposed reconstruction method, the stress responses can be used in the fatigue crack growth models based on fracture mechanics or the concept of S-N curves. The fatigue life prediction integrating the strain and stress reconstruction method is presented next.

2.2. Fatigue Life Prediction Integrating Strain and Stress Reconstruction Method

The fatigue crack growth model used in this study is a small time scale model. Since the development of the fatigue crack growth model is beyond the scope of this paper, only a brief introduction is given here. Detailed derivation and model validation for the new material fatigue crack growth model can be found in the referred article [25].
The small time scale model is developed based on the geometric relationship between the crack tip opening displacement (CTOD) and the instantaneous crack growth kinetics. The geometric relationship between the CTOD and the instantaneous crack growth kinetics is shown in Figure 2. The schematic illustration in Figure 2 is for a through thickness crack in an infinite plate. Only the tip region is shown. As shown in Figure 2, the crack will extend a distance (da), after a small time increment (dt) and the crack tip will extend from O to O’. Considering the geometry of crack tips at two time points (t and t + dt), the crack growth rate da/dt for an infinitesimal crack growth is derived as it is shown in Equation (8), where θ is the crack tip opening angle (CTOA).
d a = cot ( θ ) 2 d δ = C d δ
The CTOD can be approximately expressed as Equations (9) and (10) using the plastic zone model proposed by Irwin [26]. Equation (9) is for elastic-perfect-plastic material behavior and ignores the hardening effect [27].
δ = 1 2 K 2 E σ y = λ σ 2 a
λ = π 2 E σ y
where E is the Young’s modulus. σy is the yield strength. It should be noted that σ in Equation (9) is the nominal stress. According to Figure 2, the crack is propagated from O to O’ during a small time increment. The crack length increment is da and the CTOD increment is . The CTOD increment , is expressed as Equation (11)
d δ = λ ( 2 σ a d σ + σ 2 d a )
Substituting Equation (11) into Equation (8) and dividing on both sides by a small time increment dt, the instantaneous crack growth rate is represented as
1 C λ a d a d t = 2 σ 1 C λ σ 2 d σ d t
The proposed methodology describes crack growth rate in terms of time scale instead of cycle, according to Equation (12). The crack length at any arbitrary time can be calculated by direct time integration. The above discussion is for the case when the crack starts to grow. The crack may not grow during the entire duration of the cyclic loading. A general expression considering the non-uniform crack growth is expressed as
a ˙ = H ( σ ˙ ) H ( σ σ r e f ) 2 C λ 1 C λ σ 2 σ ˙ σ a
where “ ˙ ” denotes time derivative throughout this paper. H is the Heaviside function. σref is the reference stress level where the crack starts to grow. Cracks do not grow during the entire duration of the cyclic loading. For example, crack does not grow during the unloading path due to the energy principle. Also, crack only starts to grow when the applied loading is beyond a certain stress level. The crack length needs to be calculated by performing the integration from the lower integration limit (e.g., the time when the reference stress level is reached) to the upper integration limit (e.g., the time when the maximum stress level is reached). The crack closure model is essentially contained in the calculation of reference stress level. In this case, the determination of the reference stress level is critical in order to perform the integration. The small time scale model has been validated using the existing experimental data for various materials under both constant loading and variable amplitude loadings [25]. There are several advantages for the developed small time scale model: (1) Comparing with Paris’s model, the small time scale model does not need cycle-counting (e.g., only using direct time domain integral) under random variable loadings. Paris’s model considers the average crack growth per cycle and cannot include the detailed mechanisms within one loading cycle (e.g., non-uniform crack growth kinetics within one cycle as shown in Equation (13); (2) Stress ratio effect has been included in the small time scale model since the direct stress state instead of stress range is used; and (3) One unique advantage of the small time scale model is that it can be seamlessly coupled with structural analysis and has great potential for concurrent structural dynamic analysis and multilevel (e.g., structural level and material level) fatigue life evaluation. This capability makes it an ideal model for real-time damage analysis and on-line decision-making [21].

3. Example Validation

Two numerical examples are presented here. The first example is a three-dimensional frame structure subject to random forces. The second example is a complex airfoil structure with 19,956 DOFs. These two examples demonstrate the reconstruction methodology integrating the fatigue damage life evaluation.

3.1. A Three-Dimensional Frame Structure Example

The proposed methodology can be directly applied to structural level reconstruction analysis. To demonstrate the basic idea, a 216-DOF spatial frame structure is used here. The finite element model diagram of the structure and its dimensions are shown in Figure 3. The structure is a four-story steel structure and it has extents of 2 m, 1 m, and 5 m in the x-axis direction, y-axis direction, and z-axis direction, respectively. Each story has a height of 1.25 m. The figure nodes are labeled using nearby numbers and elements are labeled with numbers in round brackets. Each of the elements has two nodes and 12 DOFs in total. Each node has six DOFs, namely the displacements of x, y, and z directions and the rotations along x, y, and z axes. All the nodes attached to the ground (i.e., z = 0) are prescribed and the structure has 216 DOFs. Required properties of the structural member element for FE modeling are listed in Table 1.
At each floor, stochastic forces in the x-axis direction are applied to simulate the ambient excitations. The surface geometry center of the element (78) is used as the location of the actual strain gauge. Sensor measurements are obtained by solving the dynamical equations of the finite element model and three different (2%, 5%, 10%) RMS noise terms are added to the deterministic results to represent the measurement uncertainty. The sampling frequency is 1000 Hz. To simulate the realistic measured signals, the noisy components are of Gaussian pulse processes with root mean square (RMS) setting to a percentage of the largest RMS of the measured signal. Then the noise components can be quantified by the noise level in this paper. For example, a 5% noise level is to generate noise components from Gaussian pulse processes with RMS setting to 5% of the largest RMS of the strain responses. Without loss of generality, the element (52) is arbitrarily chosen to represent the location of interest. For this element, strain and stress responses at the geometry center of the surface perpendicular to the y-axis direction can have maximum response amplitudes and this geometry center is the location of interest for reconstruction. Both the strain gauge measurement location and the reconstructed location are shown in Figure 3. For illustration purposes, only the normal strain and stress are considered to be significant. The reconstruction for shear strain and stress responses can be obtained using the same methodology.
The 60 s of synthesized strain gauge measurement data with 5% RMS noise terms are presented in Figure 4a. The Fourier spectra of the strain gauge measurement data indicate four significant frequencies as shown in Figure 4c. Those four frequencies are used to design band-pass filters for the modal responses extraction. Figure 5 presents the modal responses obtained by the EMD method with intermittency criteria. The reconstructed strain gauge responses for the location of interest in element (52) and the theoretical results under 2%, 5%, 10% noise level are compared in Figure 6. The reconstructed stress responses based on the reconstructed strain results under three noise levels are shown in Figure 7. The correlation between the reconstructed strain responses and the theoretical strain calculation results is 0.9789, 0.9788, 0.9785. The reconstructed responses under the three noise levels suggests that the reconstructed responses are very consistent, which indicates that the presented integrated health monitoring method is robust against noisy measurement signals.

3.2. A Simplified Airfoil Structure Model Example

The proposed methodology can be directly applied to structural level reconstruction analysis. To demonstrate the basic idea, a 19956-DOF simplified airfoil structure model is used here. The finite element model diagram of the structure and its dimensions are shown in Figure 8. The structure is a simplified airfoil without skin and it has extents of 3.49 m, 0.11 m, and 0.95 m in the x-axis direction, y-axis direction, and z-axis direction, respectively. The element type is solid185 so that it is automatically meshed with hexahedron or (reduced) tetrahedron elements in ANSYS. Each of the elements has eight nodes and 24 DOFs in total. Each node has three DOFs, namely the displacements of x, y, and z directions. All the nodes attached to the airframe (i.e., z = 0) are prescribed and the structure has 19956 DOFs. Required properties of the FE modeling in ANSYS are listed in Table 2.
Near the end of the airfoil (node 3224 in FE model), transient force with the amplitude of 1000 N in the y-axis direction is applied to simulate the ambient excitations for the first 0.01 s. The location K (element 1847 in FE model) in Figure 8 is used as the location of the actual strain gauge. Sensor measurements are presented by transient dynamic analysis in ANSYS and 5% RMS noise terms are added to the deterministic results to represent the measurement uncertainty. The sampling frequency is 1000 Hz. For the airfoil example, the fixed end of the airfoil (Loc. I1, I2, I3) are chosen as points of interest for reconstruction. The strain gauges are placed near the free end of the airfoil (Loc. K) where the sensors are easy to set. The proposed method provides a fundamental tool to reconstruct the stress information of critical spots using remote and sparse sensor measurements. The location of critical spots can be inferred based on engineering experience and numerical simulations. Without loss of generality, the locations I1, I2, I3 (element 1082, 2558, 3043 in FE model) in Figure 8 are arbitrarily chosen to represent the location of interest for reconstruction. Both the strain gauge measurement location and the reconstructed location are shown in Figure 8.

3.2.1. Strain and Stress Response Reconstruction

Three seconds of synthesized strain gauge measurement data with 5% RMS noise terms are presented in Figure 9a. The Fourier spectra of the strain gauge measurement data indicate a significant frequency component as shown in Figure 9b. The frequency is used to design the band-pass filter for the modal responses extraction. The reconstructed strain gauge responses for the locations of interest I1, I2, I3, and the theoretical results are compared in Figure 10. The reconstructed stress responses based on the reconstructed strain results are shown in Figure 11. The correlation between the reconstructed strain responses and the theoretical strain calculation results is 0.9849, 0.9867, and 0.9865 respectively.
The discrepancies on the right end side of Figure 10 and Figure 11 are caused by the Gibbs phenomenon which has also been referred to as boundary/end effect in the EMD method [28], which has been widely discussed as the intrinsic weakness of EMD. The solution and verification of the Gibbs phenomenon is also worthy of future research. This is also a limitation of the method. Some clarification is discussed below in the conclusions.
To investigate the effect of measurement noise level to the reconstructed results, several numerical cases with various noise levels are studied. The noise components are added to the strain responses and the results are used as the representative noisy strain measurement data. Correlation coefficient is used as a metric to evaluate the similarity between the theoretical responses and the reconstructed responses for bending stresses. RMS is set to taking value from 0% to 10% with 1% increment. At each of the RMS settings, bending stress responses at location I1 (as shown in Figure 8) are reconstructed based on the strain measurement at location K (as shown in Figure 8). The correlation coefficient between the reconstructed stress responses at location I1 and the theoretical stress responses (also with noise components) at location I1 is calculated. Results for all RMS settings are presented in Figure 12.
It can be seen that the performance is degraded with the increase of noise level. The overall performance of the reconstruction method is larger than 97.5% when noise levels are not larger than 10% RMS in this example. Next, fatigue life evaluation integrating the reconstructed stress responses is presented.

3.2.2. Fatigue Life Evaluation Using the Reconstructed Stress Responses

To demonstrate the fatigue life evaluation with reconstructed stress responses, a numerical crack growth analysis using the original and reconstructed signal with different noise levels in the airfoil structure model example are performed. The strain gauges are placed at the location K, while the location I1 is the point of interest. A mode I planar crack is considered in the demonstrated examples which is reduced to a single edge notch bending problem as shown in Figure 13.
The initial crack is assumed to be 1.5 mm at location I1 in Figure 13. The given periodic force which is composed of normal stress in Figure 11 is applied at location I. For clear demonstration, the crack size versus number of cycles (a-N) curves are used for fatigue crack growth trajectory demonstration in Figure 14. Figure 14a presents the overall fatigue crack growth evaluation results using both theoretical and proposed REMD methods. Four different noise levels (1%, 2%, 5%, 10% RMS noise) are used to verify the robustness of the proposed method. Results of 2532 to 2684 cycles are also shown in Figure 14b for clear presentation of fatigue crack growth behavior. The x-axis represents the cycles in the numerical simulation. The y-axis is the simulated crack length using the proposed small time scale fatigue crack growth model. As shown in Figure 14, all crack length curves agree well with each other. In the current study, the noisy components are Gaussian pulse process with root mean square (RMS) setting to a percentage of the largest RMS of the measured signal. Therefore, as the percentage of the noise increases, it may cause a larger peak of the strain response, which will lead to an increase of the crack growth rate. Quantitative error analysis is performed and shown in Figure 15. The x-axis represents the time points. The y-axis is the maximum absolute error of the crack length prediction normalized by the crack length. The error definition is
E r r o r = max | a ^ ( t ) a ( t ) | a e ( t )
where a ^ ( t ) is the crack length at time t using reconstructed stress responses. a ( t ) is the crack length at time t using theoretical stress responses. a e ( t ) is the final crack length of calculation. It is shown that the prediction error of crack length does not have a monotonic trend with respect to the introduced noise level. The maximum error is less than 2.5% with all the stress histories investigated in the current study. It should be noted that the proposed study focuses on the theoretical development of the stress reconstruction methodology. Further investigation is required for the experimental validation studies.

4. Conclusions

In this study, an integrated fatigue life evaluation method for structural health monitoring using limited sensor data is developed. Sparse and remote strain gauge measurements from the existing structural health monitoring (SHM) system can be directly used to reconstruct the strain and stress responses for critical spots without direct sensor measurements. The proposed method extends the previous empirical mode decomposition (EMD)-based time-domain reconstruction method by performing real-time fatigue life evaluation for structural health monitoring purposes. The strain/stress reconstruction method based on empirical mode decomposition (REMD) is employed in the proposed method to extrapolate the strain and stress responses using the derived strain and stress transformation function. Mode superposition is used to obtain the strain and stress responses of the critical spot for fatigue life prediction. The overall method is demonstrated using a three-dimensional frame structure and a more complex 19956-DOF simplified airfoil structure. Fatigue analysis using the reconstructed stress responses is performed to demonstrate the effectiveness of the overall method. Based on the current study, several conclusions are drawn:
  • In this study, a novel structural fatigue life evaluation method is proposed to provide prompt, informed fatigue damage predictions of the structures based on limited strain gauge measurement. Using this method, the evolution of fatigue damage within structures can be performed using the limited strain gauge measurement.
  • According to numerical analysis results, the proposed method can produce results which are very close to theoretical solutions considering a practical noisy measurement system. The reconstructed results have an overall correlation coefficient larger than 0.975 under 10% RMS noise settings.
  • Fatigue crack growth analysis using the reconstructed stress responses agrees well with the crack size predictions using the original stress responses. The crack size prediction using the reconstructed stress responses is not sensitive to the introduced noise.
However, limitations should be stated as follows: Firstly, in the current study, locations and damage severity of critical spots are assumed to be known from either previous field data or damage detection results. On the other hand, a crack will definitely alter the reconstructed dynamic responses when it is large enough to affect the structure parameters such as stiffness and natural frequency. From this point of view, the reconstructed responses containing damage information can be used as damage detection. However, damage detection (including locating and evaluating the severity of the damage) is beyond the scope of this paper and is listed as future work. Secondly, if the high damping C is introduced, the amplitude of the measured strain signal will reduce rapidly. Then, the mode response may be extracted incompletely from the measured strain signal which will lead to the uncertainty of the strain reconstruction method.
The current study focuses on theoretical development, and further study is thus required to experimentally validate the proposed methodology. Only mode I planar crack is considered in the demonstrated examples. Curvilinear or more general 3D cracks need further investigations, and the uncertainty of the loading on the fatigue life evaluation is also an interesting study for the future.

Acknowledgments

This study is supported by the National Natural Science Foundation of China (Grant No. 11404018) and the Collaborative Innovation Center for Advanced Aero-Engine (CICAAE).

Author Contributions

Jingjing He conceived and designed the experiments; Yibin Zhou performed the experiments and wrote the paper; Wei Zhang, Xuefei Guan and Yibin Zhou analyzed the data. Weifang Zhangand Yanrong Wang monitored and led the whole research project. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kulkarni, S.; Achenbach, J. Structural health monitoring and damage prognosis in fatigue. Struct. Health Monit. 2008, 7, 37–49. [Google Scholar] [CrossRef]
  2. Beral, B.; Speckmann, H. Structural health monitoring (SHM) for aircraft structures: A challenge for system developers and aircraft manufactures. In Structural Health Monitoring 2003: From Diagnostics & Prognostics to Structural Health Management: Proceedings of the 4th International Workshop on Structural Health Monitoring, Stanford University, Stanford, CA, September 15–17, 2003; Chang, F.K., Ed.; Stanford University: Stanford, CA, USA, 2003; pp. 12–29. [Google Scholar]
  3. Chang, P.C.; Flatau, A.; Liu, S. Review paper: Health monitoring of civil infrastructure. Struct. Health Monit. 2003, 2, 257–267. [Google Scholar] [CrossRef]
  4. Giurgiutiu, V.; Cuc, A. Embedded non-destructive evaluation for structural health monitoring, damage detection, and failure prevention. Shock Vib. Dig. 2005, 37, 83–105. [Google Scholar] [CrossRef]
  5. Katsikeros, C.E.; Labeas, G. Development and validation of a strain-based structural health monitoring system. Mech. Syst. Signal Process. 2009, 23, 372–383. [Google Scholar] [CrossRef]
  6. Farrar, C.R.; Worden, K. An introduction to structural health monitoring. Philos. Trans. R. Soc. A 2007, 365. [Google Scholar] [CrossRef] [PubMed]
  7. Law, S.; Li, J.; Ding, Y. Structural response reconstruction with transmissibility concept in frequency domain. Mech. Syst. Signal Process. 2011, 25, 952–968. [Google Scholar] [CrossRef]
  8. Kammer, D.C. Estimation of structural response using remote sensor locations. J. Guid. Control Dyn. 1997, 20, 501–508. [Google Scholar] [CrossRef]
  9. Ribeiro, A.; Silva, J.; Maia, N. On the generalisation of the transmissibility concept. Mech. Syst. Signal Process. 2000, 14, 29–36. [Google Scholar] [CrossRef]
  10. He, J.; Guan, X.; Liu, Y. Structural response reconstruction based on empirical mode decomposition in time domain. Mech. Syst. Signal Process. 2012, 28, 348–366. [Google Scholar] [CrossRef]
  11. Stephens, R.I.; Fuchs, H.O. Metal Fatigue in Engineering, 2nd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 2001. [Google Scholar]
  12. Miner, M.A. Cumulative damage in fatigue. J. Appl. Mech. 1945, 12, A159–A164. [Google Scholar]
  13. Matsuishi, M.; Endo, T. Fatigue of metals subjected to varying stress. Jpn. Soc. Mech. Eng. 1968, 68, 37–40. [Google Scholar]
  14. Downing, S.D.; Socie, D. Simple rainflow counting algorithms. Int. J. Fatigue 1982, 4, 31–40. [Google Scholar] [CrossRef]
  15. Colombo, D.; Giglio, M. A methodology for automatic crack propagation modelling in planar and shell Fe models. Eng. Fract. Mech. 2006, 73, 490–504. [Google Scholar] [CrossRef]
  16. Giglio, M.; Manes, A. Crack propagation on helicopter panel: Experimental test and analysis. Eng. Fract. Mech. 2008, 75, 866–879. [Google Scholar] [CrossRef]
  17. Luo, C.T.; Mohanty, S.; Chattopadhyay, A. Fatigue damage prediction of cruciform specimen under biaxial loading. J. Intell. Mater. Syst. Struct. 2013, 24, 2084–2096. [Google Scholar] [CrossRef]
  18. Colombo, D.; Giglio, M.; Manes, A. 3d fatigue crack propagation analysis of a helicopter component. Int. J. Mater. Prod. Technol. 2007, 30, 107–123. [Google Scholar] [CrossRef]
  19. Sbarufatti, C.; Corbetta, M.; Manes, A.; Giglio, M. Sequential monte-carlo sampling based on a committee of artificial neural networks for posterior state estimation and residual lifetime prediction. Int. J. Fatigue 2016, 83, 10–23. [Google Scholar] [CrossRef]
  20. Newman, J.C., Jr. A crack-closure model for predicting fatigue crack growth under aircraft spectrum loading. ASTM Int. 1981, 748, 53–84. [Google Scholar]
  21. He, J.; Lu, Z.; Liu, Y. New method for concurrent dynamic analysis and fatigue damage prognosis of bridges. J. Bridge Eng. 2012, 17, 396–408. [Google Scholar] [CrossRef]
  22. De Koning, A. A simple crack closure model for prediction of fatigue crack growth rates under variable-amplitude loading. ASTM Int. 1981, 743, 63–85. [Google Scholar]
  23. Bao, C.; Hao, H.; Li, Z.X.; Zhu, X. Time-varying system identification using a newly improved HHT algorithm. Comput. Struct. 2009, 87, 1611–1623. [Google Scholar] [CrossRef]
  24. Huang, N.E.; Shen, Z.; Long, S.R. A new view of nonlinear water waves: The hilbert spectrum 1. Annu. Rev. Fluid Mech. 1999, 31, 417–457. [Google Scholar] [CrossRef]
  25. Lu, Z.; Liu, Y. Small time scale fatigue crack growth analysis. Int. J. Fatigue 2009, 32, 1306–1321. [Google Scholar] [CrossRef]
  26. Janssen, M.; Zuidema, J.; Wanhill, R.J.H. Fracture Mechanics, 1st ed.; VSSD: Delft, The Netherlands, 2004; p. 378. [Google Scholar]
  27. Skorupa, M. Load interaction effects during fatigue crack growth under variable amplitude loading—A literature review. Part II: Qualitative interpretation. Fatigue Fract. Eng. Mater. Struct. 1999, 22, 905–926. [Google Scholar] [CrossRef]
  28. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. A 1998, 454, 903–995. [Google Scholar] [CrossRef]
Figure 1. Overall procedure of the fatigue life evaluation method.
Figure 1. Overall procedure of the fatigue life evaluation method.
Materials 09 00894 g001
Figure 2. Schematic representation of crack tip geometry profiles during fatigue crack growth.
Figure 2. Schematic representation of crack tip geometry profiles during fatigue crack growth.
Materials 09 00894 g002
Figure 3. FE model of the spatial frame structure.
Figure 3. FE model of the spatial frame structure.
Materials 09 00894 g003
Figure 4. Strain gauge measurement data and Fourier spectra of the data. (a) Strain gauge measurement data (0–60 s); and (b) concentrated on 4–6 s for clear presentation; and (c) Fourier spectra of the measurement data (0–60 s).
Figure 4. Strain gauge measurement data and Fourier spectra of the data. (a) Strain gauge measurement data (0–60 s); and (b) concentrated on 4–6 s for clear presentation; and (c) Fourier spectra of the measurement data (0–60 s).
Materials 09 00894 g004aMaterials 09 00894 g004b
Figure 5. Modal responses of the strain gauge measurement data obtained by EMD method with intermittency criteria. Model 1: modal response for the 1st mode; Model 2: modal response for the 2nd mode; Model 3: modal response for the 3rd mode; Model 4: modal response for the 4th mode.
Figure 5. Modal responses of the strain gauge measurement data obtained by EMD method with intermittency criteria. Model 1: modal response for the 1st mode; Model 2: modal response for the 2nd mode; Model 3: modal response for the 3rd mode; Model 4: modal response for the 4th mode.
Materials 09 00894 g005
Figure 6. Reconstructed and theoretical strain responses for the location of interest under (a) 2%; (b) 5%; (c) 10% noise level. Results are concentrated on 34–38 s for clear presentation.
Figure 6. Reconstructed and theoretical strain responses for the location of interest under (a) 2%; (b) 5%; (c) 10% noise level. Results are concentrated on 34–38 s for clear presentation.
Materials 09 00894 g006
Figure 7. Reconstructed and theoretical bending stress responses for the location of interest under (a) 2%; (b) 5%; (c) 10% noise level. Results are concentrated on 34–38 s for clear presentation.
Figure 7. Reconstructed and theoretical bending stress responses for the location of interest under (a) 2%; (b) 5%; (c) 10% noise level. Results are concentrated on 34–38 s for clear presentation.
Materials 09 00894 g007aMaterials 09 00894 g007b
Figure 8. FE model of the simplified airfoil structure.
Figure 8. FE model of the simplified airfoil structure.
Materials 09 00894 g008
Figure 9. Strain gauge measurement data and Fourier spectra of the data. (a) Strain gauge measurement data (0–3 s); and (b) Fourier spectra of the measurement data (0–3 s).
Figure 9. Strain gauge measurement data and Fourier spectra of the data. (a) Strain gauge measurement data (0–3 s); and (b) Fourier spectra of the measurement data (0–3 s).
Materials 09 00894 g009
Figure 10. Reconstructed and theoretical strain responses for (a) the location I1; (b) the location I2; (c) the location I3.
Figure 10. Reconstructed and theoretical strain responses for (a) the location I1; (b) the location I2; (c) the location I3.
Materials 09 00894 g010
Figure 11. Reconstructed and theoretical bending stress responses for (a) the location I1; (b) the location I2; (c) the location I3.
Figure 11. Reconstructed and theoretical bending stress responses for (a) the location I1; (b) the location I2; (c) the location I3.
Materials 09 00894 g011
Figure 12. Reconstruction performance measured in correlation coefficient under different noise levels.
Figure 12. Reconstruction performance measured in correlation coefficient under different noise levels.
Materials 09 00894 g012
Figure 13. Single edge notch bending solid segment.
Figure 13. Single edge notch bending solid segment.
Materials 09 00894 g013
Figure 14. Crack growth trajectories calculated using theoretical and reconstructed stress responses under different RMS noise levels. (a) Entire 3500 cycles crack growth trajectories; (b) crack growth trajectories (2532–2684 cycles).
Figure 14. Crack growth trajectories calculated using theoretical and reconstructed stress responses under different RMS noise levels. (a) Entire 3500 cycles crack growth trajectories; (b) crack growth trajectories (2532–2684 cycles).
Materials 09 00894 g014
Figure 15. Error analysis of crack length growth evaluation using theoretical and reconstructed stress responses under different RMS noise levels.
Figure 15. Error analysis of crack length growth evaluation using theoretical and reconstructed stress responses under different RMS noise levels.
Materials 09 00894 g015
Table 1. Properties of structural member elements.
Table 1. Properties of structural member elements.
PropertyValue
Cross-section area B (m2)0.01
Moment of inertia Iy (m4)8.3333 × 10−6
Moment of inertia Iz (m4)8.3333 × 10−6
Torsion constant J (m4)6.6667 × 10−5
Young’s modulus E (GPa)200
Poisson’s ratio ν0.3
Shear modulus G (GPa)E/(2 + 2ν)
Mass per unit volume ρ (kg/m3)7.8 × 10−3
Table 2. Properties of FE model in ANSYS.
Table 2. Properties of FE model in ANSYS.
PropertyValue
Materialaluminum 7075
Element typesolid185
Young’s modulus E (GPa)72
Poisson’s ratio ν0.33
Mass per unit volume ρ (kg/m3)2.81 × 103
Number of element14951

Share and Cite

MDPI and ACS Style

He, J.; Zhou, Y.; Guan, X.; Zhang, W.; Wang, Y.; Zhang, W. An Integrated Health Monitoring Method for Structural Fatigue Life Evaluation Using Limited Sensor Data. Materials 2016, 9, 894. https://doi.org/10.3390/ma9110894

AMA Style

He J, Zhou Y, Guan X, Zhang W, Wang Y, Zhang W. An Integrated Health Monitoring Method for Structural Fatigue Life Evaluation Using Limited Sensor Data. Materials. 2016; 9(11):894. https://doi.org/10.3390/ma9110894

Chicago/Turabian Style

He, Jingjing, Yibin Zhou, Xuefei Guan, Wei Zhang, Yanrong Wang, and Weifang Zhang. 2016. "An Integrated Health Monitoring Method for Structural Fatigue Life Evaluation Using Limited Sensor Data" Materials 9, no. 11: 894. https://doi.org/10.3390/ma9110894

APA Style

He, J., Zhou, Y., Guan, X., Zhang, W., Wang, Y., & Zhang, W. (2016). An Integrated Health Monitoring Method for Structural Fatigue Life Evaluation Using Limited Sensor Data. Materials, 9(11), 894. https://doi.org/10.3390/ma9110894

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