Next Article in Journal
Rapid and Resilient LoRa Leap: A Novel Multi-Hop Architecture for Decentralised Earthquake Early Warning Systems
Previous Article in Journal
OS-PSO: A Modified Ratio of Exponentially Weighted Averages-Based Optical and SAR Image Registration
Previous Article in Special Issue
A Water Environment-Based Simulated Method for Ultrasonic Testing of Slag Inclusion Weld Defects Based on Improved VMD
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Plastic Damage Warning and Real-Time Sensing System Based on Cointegration Theory

School of Mechanical Engineering, Dalian University of Technology, Dalian 116024, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(18), 5961; https://doi.org/10.3390/s24185961
Submission received: 24 August 2024 / Revised: 10 September 2024 / Accepted: 12 September 2024 / Published: 13 September 2024
(This article belongs to the Special Issue Sensing and Imaging for Defect Detection)

Abstract

:
Structural damage can affect the long-term operation of equipment. Real-time damage warning for structures can effectively avoid accidents caused by structural damage. In this paper, a real-time warning method of structural plastic damage based on the cointegration theory is proposed. This method calculates the cointegration relationship between the strain signals at different measuring points, and the stability of the strain signal relationships is also evaluated. The problem of inaccurate detection caused by the error of strain measurement and environmental influence can be eliminated by the comprehensive judgment of strain between asymmetrical measuring points. A real-time damage sensing system is developed in this paper. In order to improve the real-time and practicability of the system, this paper proposes and determines the residual warning coefficient by analyzing the proportion of the strain residuals exceeding the residual threshold. The research on this sensing system has certain value for the engineering application of damage monitoring methods.

1. Introduction

The health of structural parts directly affects the service life of equipment [1,2,3]. Under the long-term action of loads, structural parts will suffer a certain degree of damage, which will affect their service life and even cause engineering accidents. Therefore, the accurate monitoring and warning of structural damage can provide a timely reference for equipment maintenance [4,5], which can effectively avoid losses caused by structural damage.
Structural damage monitoring has always been a hot topic in the field of engineering research [6,7]. Common structural damage mainly includes cracks, holes, plastic deformation, structural fractures and other forms. Researchers have conducted relevant research on different types of structural damage using different types of sensors and signal processing methods [8]. Detection methods based on eddy current [9], acoustic emission [10], piezoelectric ultrasonic guided wave [11], electrical impedance [12], visual recognition [13], fiber bragg grating [14] and other technologies each have their respective advantages in the field of structural health monitoring. In recent years, ultrasonic detection technology and machine learning methods have been widely introduced into structural health monitoring [15]. Among them, ultrasonic detection is widely used to monitor cracks in structures [16,17]. Liu et al. studied a wireless sensor based on a nonlinear ultrasonic modulation method to achieve the real-time online monitoring of fatigue cracks [18]. Kim et al. proposed a reference-free ultrasonic guided-wave crack detection method using the polarization characteristics of piezoelectric wafers [19]. With the continuous optimization and development of machine learning methods, they have been applied to monitor crack damages in structures [20,21]. Wang et al. predicted the growth law of fatigue cracks based on the particle filter model [22]. Chen et al. identified cracks in bridge structures based on a deep learning algorithm and realized the integrated application of the algorithm [23]. Zhou et al. used laser ultrasound and the YOLOv5 framework to monitor and identify structural cracks in interference fits [24]. The continuous development of sensor technology has also provided new methods for monitoring crack damage. Xu et al. studied a new layer sensor that monitors the crack extension of steel bridges by different signals at different fatigue stages [25]. Pang et al. studied a sensor based on microstrip antennas that can detect cracks under layers [26]. Strain can accurately reflect changes in structural stress, and the method of identifying damage through strain in a structure has certain advantages. Yu et al. studied full-optical strain FBG with tunable DFB laser demodulation for monitoring tensile and fatigue damage [27]. Pauw et al. used embedded optical fiber sensors to monitor the fatigue damage of structures [28]. Different detection technologies have their own advantages, and the monitoring of structural damage through the fusion of multiple sensors is also being studied by scholars [29]. Qi et al. monitored crack growth during fatigue by integrating strain sensors and piezoelectric guided waves [30]. Wang et al. used a weighted adaptive Kalman filtering-based method to integrate optical fiber and guided-wave information to identify structural fatigue cracks [31]. In summary, most studies on structural damage monitoring focus on monitoring macroscopic damage such as cracks. In actual engineering applications, the plastic deformation of materials often occurs before cracks initiate, and even before obvious cracks appear, the structure has already been destroyed. Therefore, the monitoring and warning of structural damage before crack initiation are of great significance for early maintenance of structures.
Based on the cointegration theory, an efficient and real-time plastic damage warning method is proposed in this paper, which can accurately identify the plastic damage of structures. The method identifies structural plastic damage through the change in the cointegration residual of the strain series between the measuring points. By analyzing the proportion of residuals exceeding the residual threshold, this paper proposes a warning coefficient. And the stability of the real-time alarming of the damage sensing system is improved by comparing the residual and warning coefficient. This method further realizes the application of damage monitoring theory in practical engineering.

2. Cointegration Theory

Time series can be divided into two types: unstable and stable time series. Unstable time series refers to a series in which the mean, variance and other statistical data change with time, while stable time series refers to a series in which the statistical data do not change with time. Most of the strain data for structural damage monitoring are unstable series that change with time, and their linear combination may become a series that does not change with time.
For an unstable time series y t , after d times of difference, it becomes a stationary time series, which is integrated of order d, denoted by y t ~ I d . For an unstable variable Y = y 1 , y 2 , , y n T , its components are all integrated of order 1. If there is a set of linear combination coefficients β = β 1 , β 2 , β n , which makes the previous Y become a stationary series, then y 1 , y 2 , , y n satisfies the cointegration relationship, and the equation is shown below.
β T Y = β 1 y 1 + β 2 y 2 + + β n y n   ~   I
where β = β 1 , β 2 , β n , βTY is residual covariance, and it is a Gaussian noise series with a mean equal to zero.
From the above formula, we can see that the main purpose of cointegration is to seek a long-term stable equilibrium relationship between two sets of unstable time series. If two data series become stable after the first-order difference calculation and are also stable after linear combination, it means that there is a cointegration relationship between them.
An unstable series is more likely to cause pseudo-regression. The significance of the cointegration test is to determine whether the relationship described by its regression equation is pseudo-regression. Therefore, before solving the cointegration residual, it is necessary to test the stationarity of the time series. Commonly used methods include the ADF (Augmented Dickey–Fuller) test and E-G (Engel–Granger) test.

2.1. Augmented Dickey–Fuller (ADF) Method

The ADF method is used to test the stability of time series. If there is a time series y t , its ADF model is as follows:
y t = ρ y t 1 + a + i = 1 p 1 b i y i 1 + ε t
where y i 1 is the difference component of y t and p and bi are coefficients. The a represents the constant term. The ρ is the lag order, and ε t is the residual series. The coefficient ρ is used to determine the stationarity of the time series y t . When each series or its difference series is integrated of order 1, the requirements of the cointegration equation are met, and the next step is to continue the E-G test to complete the subsequent solution.

2.2. Engel–Granger (E-G) Method

The theoretical basis of the E-G test is to perform a unit root test on the residuals of the regression equation. If the dependent variable can be described by various linear combinations of the independent variables, it means that there is a relatively stable equilibrium relationship between the two. The part that cannot be described forms a residual series. There is no need for serial correlation between the residual series, so the residuals must be stationary. Therefore, the purpose of the E-G method is to verify whether its residual series is stable. The steps of this method are as follows.
Assuming that x t and y t are both integrated of order 1, the residual difference between them is calculated as follows.
ε t = y t β x t α
After the ADF test, if the residual series ε t is stable, it means that there is a cointegration relationship between the original series x t and y t , where α and β are the corresponding cointegration coefficients.
The existence of a cointegration relationship means that the two sets of time series have a stable relationship in the long term, that is, their long-term equilibrium states are interdependent. The cointegration residual is a stable series obtained after eliminating unstable factors. If the cointegration residual changes suddenly, it means that the cointegration relationship has changed, which can be used as a key indicator for damage warning.

3. Process of Structural Damage Warning Method

By loading the sample, the strain signals at different measuring points are tested. The residual threshold of damage warning is calculated by the strain signals of the measuring points in the undamaged state. The residual threshold is calculated according to the following method:
R T U = μ + 3 σ R T L = μ 3 σ
where the RTU is the upper limit of the residual threshold, the RTL is the lower limit of the residual threshold, μ is the residual mean and σ is the residual variance. The flowchart of the damage warning method is shown in Figure 1. First, the stability of the strain data is tested by the ADF method. If the two sets of test strains are unstable, this method is used to evaluate the stability of the corresponding difference series. If the first-order difference of the strain series is stable, the series is integrated of order 1. It meets the requirements of the cointegration equation, and then the E-G test is performed on it. If the E-G result shows that the residual is stable, the residual equation can be constructed. After calculating the residuals of the two series, the ADF stability test is performed on the residual series. If the residuals are unstable, it means that the stability of the two series is destroyed, which also means that structural damage has occurred. If the residuals are stable, the strain series will continue to be calculated to realize the warning of structural damage.

4. Test System and Experiment Analysis

4.1. Test System and Samples

As shown in Figure 2, the test system consists of a sample, testing machine, load control terminal, strain collecting instrument and sensing system. The load control terminal controls the testing machine to apply the load to the sample. The strain collecting instrument transmits the collected strain data to the sensing system. The sensor system analyzes the residual of the strain series of multiple channels through the co-arrangement theory and finally realizes the alarm of sample damage.
Figure 3 shows the sample, which is made of Q345D structural steel. Sample 1 is the stress concentration sample, and sample 2 is the standard sample. Strain gauges are pasted in the center and on both sides of samples 1 and 2. The detailed dimensions of the structure and strain gauge pasting are shown in Figure 3a. The thickness of the samples is 5 mm. The actual pictures of the two samples are shown in Figure 3b.

4.2. Analysis of Experimental Results

By applying tensile force to the sample, the strain signal of the sample measuring point is tested. Each strain channel is synchronously sampled with a sampling rate of 200 Hz. The strain test results of different measuring points of the two samples are shown in Figure 4 and Figure 5, respectively. The signal of the entire test process is divided into three stages: the no-load stage, the elastic deformation stage and the plastic deformation stage. As shown in Figure 4a, 0–45 s is the no-load stage, and the strain signal tested at this time is the noise of the system. Next, 45–82 s is the elastic deformation stage. At 45 s, the test machine starts to load the sample with a loading rate of 0.5 kN/s. The enlarged curve of the no-load stage and the elastic deformation stage is shown in Figure 4b. It can be seen from the figure that the strain signals of each measuring point from 45 to 82 s show a linear growth trend. It can be judged by Hooke’s law that this stage belongs to the elastic deformation stage. After 82 s of continuous loading, the sample slowly enters the plastic deformation stage. It can be seen from Figure 4a that the strain of each measuring point increases rapidly after entering the plastic zone. The enlarged view of Figure 4b shows that the strain growth rate of measuring point 2 located at the stress concentration position is greater than the strain growth rate of the symmetrically arranged measuring points 1 and 3.
The strain changes of sample 2 during the entire loading process are also analyzed based on the strain test data in Figure 5, and the loading parameters are consistent with those of sample 1. In Figure 5a, the test machine starts loading at 27 s, the strain and load of the sample are linearly related before 60 s and the sample is in the elastic deformation area. It can be clearly seen from Figure 5a,b that after 60 s, with the continuous loading of the load, the strains of the three measuring points of sample 2 gradually change and increase. After 60 s, the sample begins to slowly undergo plastic deformation, and the strain also increases rapidly until the strain gauge fails and the structure is damaged. Similarly, the strain growth rate of measuring point 2 at the center position is greater than the strain growth rate of symmetrical points 1 and 3 when entering the plastic deformation stage. For engineering structural parts, they can work normally in the elastic deformation stage and can still return to the initial state after unloading. Once the structure undergoes plastic deformation, it will cause irreversible damage to the structure and cannot return to the initial state after unloading. Irreversible structural damage will occur after the structural part enters the plastic stage, so accurate monitoring of it can provide timely warning of structural damage.
According to the cointegration theory, the strain series needs to be integrated of order 1 in order to conduct a stability assessment. The strain and force changes of the structural parts are linear in elastic deformation. Figure 6 is a first-order difference series of the strain signal at each measuring point. Sample 1 is in the no-load and elastic deformation stages from 0–82 s, and the first-order difference series of the two stages shows the same law. The strain first-order difference series of sample 2 also maintains the same law of the first-order difference in the no-load and elastic deformation stages from 0 to 60 s. After calculating the difference series, the stability of the difference series needs to be checked.
The strain signals of the first 25 s of sample 1 and sample 2 are selected as the stability test data. The ADF test is performed on each strain series and the first-order difference series of the strain, respectively, and the verification results are shown in Table 1 and Table 2. Among them, y1, y2 and y3 are the strain series corresponding to different measuring points of the sample, and Δy1, Δy2 and Δy3 are the first-order difference series of the corresponding strains. By comparing the ADF check value with the ADF critical value at the 1% and 5% significance levels, if the ADF check value of the series is less than the critical value, it means that the series is stable. If the ADF check value of the series is greater than the critical value, the series is unstable. It can be seen from the ADF test results in Table 1 and Table 2 that the ADF test values of series y1, y2 and y3 are all greater than the critical values of 1% and 5%, so the strain series are all unstable. The first-order difference series Δy1, Δy2 and Δy3 are all stable after verification. The strain series of sample 1 and sample 2 are both unstable after verification, and their first-order series are stable. Therefore, the strain series of the two samples are all integrated of order 1, which meets the cointegration requirements.
The threshold of the residual is calculated using the selected data. The cointegration equation of the strain series between different measuring points is constructed by the E-G method, and then the residual is calculated by the cointegration equation. The parameters of the cointegration equation of the strain series of each measuring point are shown in Table 3 and Table 4. Strain series #1–#2 represents the combination of data from measuring points 1 and 2. The thresholds of the strain residual of sample 1 and sample 2 are calculated by Formula (4), and its parameters are shown in Table 5 and Table 6. In addition, the threshold mean value calculated by the three groups of data is calculated in the table, and it is used to judge whether the limit is exceeded.
The ADF method is used to test the stability of the first-order different residual series. From the test results of Table 7 and Table 8, by comparing with the 1% and 5% critical values, the ADF test results of the residuals between each series are all less than the critical value, satisfying the stable condition. And it can be seen that the cointegration residuals between different strain series are all stable. Therefore, there is a cointegration relationship between the strains of the three measurement points of sample 1 and sample 2.
The corresponding residual series is calculated for the strains of the three stages. By comparing the residual calculation results with the residual threshold, the warning of plastic deformation damage is realized. Figure 7 and Figure 8 are the residual calculation results of the strain series between the measuring points of sample 1 and sample 2, respectively. For sample 1, the stage of 0–82 s corresponds to the no-load and elastic deformation stages, and its curve is enlarged, as shown in Figure 7b. It can be seen from the enlarged residual diagram that the residual calculation results of the two groups #1–#2 and #2–#3 do not exceed the RTU and RTL before the plastic deformation stage, showing a stable state. With the continuous loading of force, the structural parts began to slowly enter the plastic deformation stage after 82 s. At this time, the relationship between strain and force at each measuring point becomes a nonlinear relationship. Since measuring points 1 and 3 are arranged on the symmetrical side of the stress concentration part of sample 1, the cointegration residual of #1–#3 changes little after plastic deformation. The residual verification results show that the strain cointegration relationship between the stress concentration measuring point and the symmetrical measuring point is very sensitive to plastic deformation. After plastic deformation occurs, the residual value will deviate from the original stable state. As the degree of plastic deformation increases, the residual value gradually increases. The change in the structural damage state can be clearly seen through the change in the residual value.
The residual calculation results of sample 2 are analyzed. As shown in Figure 8, before 60 s of loading, the residual results calculated for the sample in the no-load and elastic deformation stages do not exceed the RTU and RTL. As shown in Figure 8b, after entering the plastic deformation stage, the residual calculation results increased significantly, and the residual stability is destroyed at this time. After the material enters the plastic deformation, the strain relationship changes at the center measuring point and the symmetrical measuring point are very obvious. However, for the symmetrically arranged #1–#3 combination, the residual changes are small just after entering the plastic deformation. In the later stage of plastic deformation, the residual changes significantly. The structural damage warning through residual values needs to be calculated under the action of dynamic loading. After 160 s and 130 s in Figure 7 and Figure 8, the loading stops, the strains are all stable values of the mean and the residual cannot be judged at this time.
Accurate warning of structural damage through the cointegration method requires calculation and analysis based on the relatively stable relationship between the strains of the two measuring points. For the measuring points of the structural parts, the strains of the asymmetric points of the structure should be selected to evaluate the stability of the cointegration relationship. Within the elastic deformation range of the structure, the strains of the two points are linearly related, and their residuals show a certain stability. After the structure enters the plastic deformation, the strains of the asymmetric measuring points will obviously destroy their original stable relationship, and then the warning of structural damage is performed by comparison with the residual threshold.

5. Real-Time Damage Sensing System

In order to ensure the stability of the system alarm, the warning coefficient q is proposed in this paper. The warning coefficient is obtained by analyzing the data. In addition, the warning coefficient should be determined according to the sampling data and sampling frequency. The warning coefficient is the upper warning limit of the proportion of data exceeding the residual threshold in the sampled data. A too-large warning coefficient will lead to no warning after damage occurs in the structure. A too-small warning coefficient may be disturbed by unstable factors such as noise in the actual load, resulting in a false alarm. In Figure 9 and Figure 10, the regions in which the residuals exceed the residual thresholds are marked. In 3000 data samples, the proportion of residuals that exceed the threshold is calculated, as shown in Table 9. According to the data analysis in the table, the smallest proportion is 0.19. In practical applications, the residual threshold setting of the damage sensing system needs to be less than 0.19. Since the sampling data of the real-time sensing system is phased sampling, in order to improve the stability of the system alarm, the warning coefficient of the two samples is set to 0.05.
The warning of structural damage is inseparable from the real-time calculation and interactive display of a sensing system. In this paper, an intelligent damage sensing system is developed by combining the hybrid programming of LabVIEW 2015 and MATLAB R2023b. The workflow of the sensing system is shown in Figure 11a. The strain-collecting instrument collects the strain data of the sample in real time and transmits the data to the sensing system. And the sensing system packages the strain data of three measuring points into n sampling data and transmits them to the MATLAB module. According to the analysis conclusion of this paper, the strain series of measuring points 1 and 2 (#1–#2) and measuring points 2 and 3 (#2–#3) are combined. In order to improve the stability of the alarming function of the sensing system, the system counts the number of m1 (#1–#2) and m2 (#2–#3) of the calculated residuals exceeding the threshold in n sampling data in real time. The accuracy of the residual calculated for small sampling data will decrease, and too much sampling data will affect the real-time alarming of the sensing system. Through testing and analysis, the n of the sensing system is 3000. The MATLAB module calculates the corresponding residuals based on the two sets of data sampled from #1–#2 and #2–#3. Then, the proportion of the two sets of residuals in n data is counted in LabVIEW. If the proportion of any set exceeding the threshold is greater than q, the result will be reported to the sensor system for damage alarming. Through the continuous sampling and analysis of the sensing system, the warning of structural damage can be achieved in real time. Figure 11b,c are the interface displays of the system damage warning for sample 1 and sample 2, respectively.
Figure 12 shows the monitoring value of the warning coefficient of the structural damage sensing system. In the no-load and elastic deformation stages, the damage warning coefficient is 0, and the system does not issue an alarm. In Figure 12, all 3000 data samples correspond to a sampling stage. From the beginning of the damage to structural destruction, the damage warning coefficient basically increases gradually. Once the structure is destroyed, the stress change in the structure is a constant value, and the damage warning coefficient becomes 0. It is worth noting that in the third warning stage of sample 1, both damage coefficients become 0, indicating that after the damage occurs in the first stage, the strain in the structure produces a stable fluctuation, causing the residual to return to a stable state. After continuous loading, the warning coefficient of the structure is greater than the warning coefficient of the second stage, and the warning of damage can still be achieved. It can be seen from the curve of the warning coefficient that under the condition of a sampling number of 3000, a sampling frequency of 200 Hz and a warning coefficient of 0.05, the damage warning accuracy of sample 1 is 75%, and the damage warning accuracy of sample 2 is 100%. Therefore, different structural forms will have a certain impact on the accuracy of the warning system. A comparison of existing damage monitoring studies is presented in Table 10.

6. Conclusions

Based on the cointegration theory, this paper proposes a structural damage warning method. This method realizes the warning of structural damage by calculating the cointegration residual between two measuring points in the structure, avoiding the uncertainty of the warning of the single strain threshold method. The strain cointegration residuals of stress-concentrated samples and standard samples are analyzed. The residual analysis results show that the residual of the strain series between two asymmetric measuring points is more sensitive to the plastic deformation damage of the structure. The calculated residual is compared with the threshold value in the residual steady state to realize the identification of structural damage. In order to improve the engineering practicality of the method, this paper combines LabVIEW and MATLAB to develop a real-time structural damage sensing system. The system calculates the residual warning coefficients of the two groups of strain series #1–#2 and #2–#3 in real time and comprehensively considers the two prediction results to warn of structural damage. The test results show that under the conditions of a sampling number of 3000, a sampling frequency of 200 Hz and a warning coefficient of 0.05, the system alarm accuracy for the stress concentration sample is 75%, and for the standard sample it is 100%. Due to different structural characteristics, in actual applications, the system prediction accuracy for different types of structures will vary. The proposed method further realizes the engineering application of structural damage monitoring.

Author Contributions

Conceptualization, Q.G. and J.H.; methodology, Q.G. and Y.W.; software, X.W. and C.W.; writing—original draft, Q.G.; writing—review and editing, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China under grant number 52275236, Liaoning Province major science and technology grant number 2022HJ1/10400031, and Liaoning Province science and technology plan joint plan grant number 2023JH2/101700286.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

Thanks to the authors of this paper for their work on the research content. The authors are thankful for the experimental platform and test conditions provided by the project.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Achenbach, J.D. Structural Health Monitoring—What is the prescription? Mech. Res. Commun. 2009, 36, 137–142. [Google Scholar] [CrossRef]
  2. Li, S.; Liu, F.; Peng, G.; Cheng, F.; Zhao, B.; Ji, M. A lightweight SHM framework based on adaptive multisensor fusion network and multigeneration knowledge distillation. IEEE Trans. Instrum. Meas. 2022, 71, 3526919. [Google Scholar] [CrossRef]
  3. Kralovec, C.; Schagerl, M. Review of structural health monitoring methods regarding a multi-sensor approach for damage assessment of metal and composite structures. Sensors 2020, 20, 826. [Google Scholar] [CrossRef] [PubMed]
  4. Moradi, M. MEMS multisensor intelligent damage detection for wind turbines. IEEE Sens. J. 2015, 15, 1437–1444. [Google Scholar] [CrossRef]
  5. Etxaniz, J. Ultrasound-based structural health monitoring methodology employing active and passive techniques. Eng. Fail. Anal. 2023, 146, 107077. [Google Scholar] [CrossRef]
  6. Arcadius Tokognon, C.; Gao, B.; Tian, G.Y.; Yan, Y. Structural health monitoring framework based on internet of things: A survey. IEEE Internet Things J. 2017, 4, 619–635. [Google Scholar] [CrossRef]
  7. Janeliukstis, R.; Mironovs, D. Smart composite structures with embedded sensors for load and damage monitoring—A review. Mech. Compos. Mater. 2021, 57, 131–152. [Google Scholar] [CrossRef]
  8. Broer, A.A.R.; Benedictus, R.; Zarouchas, D. The need for multi-sensor data fusion in structural health monitoring of composite aircraft structures. Aerospace 2022, 9, 183. [Google Scholar] [CrossRef]
  9. Ha, N.; Lee, H.-S.; Lee, S. Development of a wireless corrosion detection system for steel-framed structures using pulsed eddy currents. Sensors 2021, 21, 8199. [Google Scholar] [CrossRef]
  10. Joseph, R.; Giurgiutiu, V. Analytical and experimental study of fatigue-crack-growth ae signals in thin sheet metals. Sensors 2020, 20, 5835. [Google Scholar] [CrossRef]
  11. Na, W.S.; Baek, J. A review of the piezoelectric electromechanical impedance based structural health monitoring technique for engineering structures. Sensors 2018, 18, 1307. [Google Scholar] [CrossRef] [PubMed]
  12. Migot, A.; Saaudi, A.; Giurgiutiu, V. Delamination depth detection in composite plates using the lamb wave technique based on convolutional neural networks. Sensors 2024, 24, 3118. [Google Scholar] [CrossRef] [PubMed]
  13. Castro-Toscano, M.J.; Rodriguez-Quinonez, J.C.; Sergiyenko, O.; Flores-Fuentes, W.; Ramirez-Hernandez, L.R.; Hernandez-Balbuena, D.; Lindner, L.; Rascon, R. Novel sensing approaches for structural deformation monitoring and 3D measurements. IEEE Sens. J. 2021, 21, 11318–11328. [Google Scholar] [CrossRef]
  14. Soman, R. Sequential niching particle swarm optimization algorithm for localization of multiple damage locations using fiber bragg grating sensors. NDT E Int. 2024, 143, 103069. [Google Scholar] [CrossRef]
  15. Khani, M.M.; Vahidnia, S.; Ghasemzadeh, L.; Ozturk, Y.E.; Yuvalaklioglu, M.; Akin, S.; Ure, N.K. Deep-learning-based crack detection with applications for the structural health monitoring of gas turbines. Struct. Health Monit. 2019, 19, 1440–1452. [Google Scholar] [CrossRef]
  16. Jiao, P.; Egbe, K.-J.I.; Xie, Y.; Matin Nazar, A.; Alavi, A.H. Piezoelectric sensing techniques in structural health monitoring: A state-of-the-art review. Sensors 2020, 20, 3730. [Google Scholar] [CrossRef]
  17. Qing, X.; Li, W.; Wang, Y.; Sun, H. Piezoelectric transducer-based structural health monitoring for aircraft applications. Sensors 2019, 19, 545. [Google Scholar] [CrossRef]
  18. Liu, P.; Lim, H.J.; Yang, S.; Sohn, H.; Lee, C.H.; Yi, Y.; Kim, D.; Jung, J.; Bae, I. Development of a “stick-and-detect” wireless sensor node for fatigue crack detection. Struct. Health Monit. 2017, 16, 153–163. [Google Scholar] [CrossRef]
  19. Kim, S.B.; Sohn, H. Instantaneous reference-free crack detection based on polarization characteristics of piezoelectric materials. Smart Mater. Struct. 2007, 16, 2375–2387. [Google Scholar] [CrossRef]
  20. Yoon, J. Deep neural network-based structural health monitoring technique for real-time crack detection and localization using strain gauge sensors. Sci. Rep. 2022, 12, 20204. [Google Scholar] [CrossRef]
  21. Azimi, M.; Eslamlou, A.; Pekcan, G. Data-driven structural health monitoring and damage detection through deep learning: State-of-the-art review. Sensors 2020, 20, 2778. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, L.; Zhang, C.; Tao, C.; Ji, H.; Yang, Y.; Qiu, J. Prediction of multiple fatigue crack growth based on modified paris model with particle filtering framework. Mech. Syst. Signal Process. 2023, 190, 110124. [Google Scholar] [CrossRef]
  23. Chen, W.; He, Z.; Zhang, J. Online monitoring of crack dynamic development using attention-based deep networks. Autom. Constr. 2023, 154, 105022. [Google Scholar] [CrossRef]
  24. Zhou, Y.; Chen, X.; Li, H.; Hao, S.; Yan, X.; Zhao, Y. Automatic identification and location method for shaft-hole interference fit crack based on deep learning. Mech. Syst. Signal Process. 2023, 192, 110194. [Google Scholar] [CrossRef]
  25. Xu, W.; Cui, C.; Luo, C.; Zhang, Q. Fatigue crack monitoring of steel bridge with coating sensor based on potential difference method. Constr. Build. Mater. 2022, 350, 128868. [Google Scholar] [CrossRef]
  26. Pang, Q.; Dong, G.; Yang, X. Metal crack detection sensor based on microstrip antenna. IEEE Sens. J. 2023, 23, 8375–8384. [Google Scholar] [CrossRef]
  27. Yu, J.; Li, C.; Qiu, X.; Chen, H. A full-optical strain fbg sensor for in-situ monitoring of fatigue stages via tunable DFB laser demodulation. Opt. Quantum Electron. 2021, 53, 156. [Google Scholar] [CrossRef]
  28. De Pauw, B.; Hinderdael, M.; Moonens, M.; De Baere, D.; Geernaert, T.; Berghmans, F.; Guillaume, P. Fatigue failure monitoring of 316L stainless steel coupons using optical fibre based distributed strain sensing. Smart Mater. Struct. 2019, 28, 105054. [Google Scholar] [CrossRef]
  29. Yu, Y.; Liu, X.; Yan, J.; Wang, Y.; Qing, X. Real-time life-cycle monitoring of composite structures using piezoelectric-fiber hybrid sensor network. Sensors 2021, 21, 8213. [Google Scholar] [CrossRef]
  30. Qi, C.; Weixi, Y.; Jun, L.; Heming, G.; Yao, M. A research on fatigue crack growth monitoring based on multi-sensor and data fusion. Struct. Health Monit. 2021, 20, 848–860. [Google Scholar] [CrossRef]
  31. Wang, Y.; He, M.; Sun, L.; Wu, D.; Wang, Y.; Qing, X. Weighted adaptive kalman filtering-based diverse information fusion for hole edge crack monitoring. Mech. Syst. Signal Process. 2022, 167, 108534. [Google Scholar] [CrossRef]
  32. Feng, L.; Qian, X. Enhanced real-time crack monitoring and updating in welded structural components with limited measurement data. Mech. Syst. Signal Process. 2023, 198, 110403. [Google Scholar] [CrossRef]
  33. Malik, S.; Kontsos, A. A novel information entropy approach for crack monitoring leveraging nondestructive evaluation sensing. Mech. Syst. Signal Process. 2024, 214, 111207. [Google Scholar] [CrossRef]
  34. Zhu, W.; Xu, Z.; Xiang, Y.; Liu, C.; Deng, M.; Qiu, X.; Sun, D.; Xuan, F. Nonlinear ultrasonic detection of partially closed cracks in metal plates using static component of lamb waves. NDT E Int. 2021, 124, 102538. [Google Scholar] [CrossRef]
  35. Fan, X.; He, Y.; Chen, T. Research on fatigue crack quantitative monitoring based on eddy current sensor with an interactive induction coil layout. Smart Mater. Struct. 2023, 32, 035013. [Google Scholar] [CrossRef]
  36. Qiu, L.; Lin, X.; Wang, Y.; Yuan, S.; Shi, W. A mechatronic smart skin of flight vehicle structures for impact monitoring of light weight and low-power consumption. Mech. Syst. Signal Process. 2020, 144, 106829. [Google Scholar] [CrossRef]
Figure 1. Flowchart of plastic damage warning method based on cointegration theory.
Figure 1. Flowchart of plastic damage warning method based on cointegration theory.
Sensors 24 05961 g001
Figure 2. Test system diagram.
Figure 2. Test system diagram.
Sensors 24 05961 g002
Figure 3. Strain test structure and strain gauge installation position: (a) sample structure dimensions and strain gauge installation position; (b) sample actual picture (unit: mm).
Figure 3. Strain test structure and strain gauge installation position: (a) sample structure dimensions and strain gauge installation position; (b) sample actual picture (unit: mm).
Sensors 24 05961 g003
Figure 4. Strain test results of measuring points 1, 2 and 3 of sample 1: (a) strain at different measuring points during the entire loading process; (b) strain enlarged diagram of the no-loading stage and the plastic deformation stage.
Figure 4. Strain test results of measuring points 1, 2 and 3 of sample 1: (a) strain at different measuring points during the entire loading process; (b) strain enlarged diagram of the no-loading stage and the plastic deformation stage.
Sensors 24 05961 g004
Figure 5. Strain test results of measuring points 1, 2 and 3 of sample 2: (a) strain at different measuring points during the entire loading process; (b) strain enlarged diagram of the no-loading stage and the plastic deformation stage.
Figure 5. Strain test results of measuring points 1, 2 and 3 of sample 2: (a) strain at different measuring points during the entire loading process; (b) strain enlarged diagram of the no-loading stage and the plastic deformation stage.
Sensors 24 05961 g005
Figure 6. First-order difference series of strain at measuring points 1, 2 and 3: (a) strain difference series of sample 1; (b) strain difference series of sample 2.
Figure 6. First-order difference series of strain at measuring points 1, 2 and 3: (a) strain difference series of sample 1; (b) strain difference series of sample 2.
Sensors 24 05961 g006
Figure 7. Calculation results of residual at the measuring point of sample 1: (a) calculation results of residual during the whole process; (b) enlarged residual during the no-load stage and the elastic deformation stage.
Figure 7. Calculation results of residual at the measuring point of sample 1: (a) calculation results of residual during the whole process; (b) enlarged residual during the no-load stage and the elastic deformation stage.
Sensors 24 05961 g007
Figure 8. Calculation results of residual at the measuring point of sample 2: (a) calculation results of residual during the whole process; (b) amplified residual during the unloaded stage and the elastic deformation stage.
Figure 8. Calculation results of residual at the measuring point of sample 2: (a) calculation results of residual during the whole process; (b) amplified residual during the unloaded stage and the elastic deformation stage.
Sensors 24 05961 g008
Figure 9. The regions of sample 1 in which the residual exceeds the residual thresholds: (a) the labeled figure of region 1; (b) the labeled figure of region 2.
Figure 9. The regions of sample 1 in which the residual exceeds the residual thresholds: (a) the labeled figure of region 1; (b) the labeled figure of region 2.
Sensors 24 05961 g009
Figure 10. The regions of sample 2 in which the residual exceeds the residual thresholds: (a) the labeled figure of region 1; (b) the labeled figure of region 2.
Figure 10. The regions of sample 2 in which the residual exceeds the residual thresholds: (a) the labeled figure of region 1; (b) the labeled figure of region 2.
Sensors 24 05961 g010
Figure 11. Structural damage sensing system: (a) system workflow diagram; (b) damage warning result of sample 1; (c) damage warning result of sample 2.
Figure 11. Structural damage sensing system: (a) system workflow diagram; (b) damage warning result of sample 1; (c) damage warning result of sample 2.
Sensors 24 05961 g011
Figure 12. Warning coefficient of structural damage sensing system: (a) warning coefficient of sample 1; (b) warning coefficient of sample 2.
Figure 12. Warning coefficient of structural damage sensing system: (a) warning coefficient of sample 1; (b) warning coefficient of sample 2.
Sensors 24 05961 g012
Table 1. ADF check results of strain series and first-order difference series of sample 1.
Table 1. ADF check results of strain series and first-order difference series of sample 1.
SeriesADF Test Value1% Critical Value5% Critical ValueStability
y1−1.0398−2.5690−1.9416Unstable
y1−99.1729−2.5690−1.9416Stable
y2−1.0791−2.5690−1.9416Unstable
y2−101.1411−2.5690−1.9416Stable
y30.8419−2.5690−1.9416Unstable
y3−95.9264−2.5690−1.9416Stable
Table 2. ADF check results of strain series and first-order difference series of sample 2.
Table 2. ADF check results of strain series and first-order difference series of sample 2.
SeriesADF Test Value1% Critical Value5% Critical ValueStability
y1−0.1069−2.5690−1.9416Unstable
y1−121.2684−2.5690−1.9416Stable
y2−0.2294−2.5690−1.9416Unstable
y2−123.6009−2.5690−1.9416Stable
y3−1.1644−2.5690−1.9416Unstable
y3−111.5629−2.5690−1.9416Stable
Table 3. Cointegration equation parameters for different strain series of sample 1.
Table 3. Cointegration equation parameters for different strain series of sample 1.
Strain Series β α
#1–#21.02171.7590
#1–#30.618345.1140
#2–#30.605542.4706
Table 4. Cointegration equation parameters for different strain series of sample 2.
Table 4. Cointegration equation parameters for different strain series of sample 2.
Strain Series β α
#1–#21.526422.3674
#1–#3−0.1026−126.8029
#2–#3−0.0400−95.5583
Table 5. Residual mean, variance and upper and lower limits of different strain series of sample 1.
Table 5. Residual mean, variance and upper and lower limits of different strain series of sample 1.
Strain SeriesμσRTURTL
#1–#2−1.3747 × 10−40.43731.3118−1.3121
#1–#3−0.00480.50021.4958−1.5054
#2–#3−0.00460.49431.4785−1.4876
Mean valueN/AN/A1.4287−1.4350
Table 6. Residual mean, variance and upper and lower limits of different strain series of sample 2.
Table 6. Residual mean, variance and upper and lower limits of different strain series of sample 2.
Strain SeriesμσRTURTL
#1–#2−7.7404 × 10−40.74472.2332−2.2348
#1–#35.5210 × 10−40.42541.2769−1.2658
#2–#36.5017 × 10−40.41791.2543−1.2530
Mean valueN/AN/A1.5881−1.5845
Table 7. ADF test results of strain residual series of sample 1.
Table 7. ADF test results of strain residual series of sample 1.
Strain SeriesADF Test Value1% Critical Value5% Critical ValueStability
#1–#2−71.4768−2.5821−1.9416Stable
#1–#3−119.5111−2.5821−1.9416Stable
#2–#3−120.6987−2.5821−1.9416Stable
Table 8. ADF test results of strain residual series of sample 2.
Table 8. ADF test results of strain residual series of sample 2.
Strain SeriesADF Test Value1% Critical Value5% Critical ValueStability
#1–#2−106.9975−2.5821−1.9416Stable
#1–#3−119.1516−2.5821−1.9416Stable
#2–#3−116.0702−2.5821−1.9416Stable
Table 9. The number and proportion of residuals exceeding the threshold in different regions.
Table 9. The number and proportion of residuals exceeding the threshold in different regions.
SamplesRegions#1–#2Proportion#1–#3Proportion#2–#3Proportion
Sample 1Region 15700.19 006520.22
9770.33 0011410.38
Region 210860.36 0022530.75
21720.72 0011130.37
Sample 2Region 100008140.27
7090.24 0000
Region 228570.95 13040.43 32761.09
15180.51 16490.55 23570.79
Table 10. Comparison of studies on damage monitoring.
Table 10. Comparison of studies on damage monitoring.
StudiesSensorMonitoring Damage TypeMonitoring SystemReal Time
Qian et al. [32]Strain sensorCrackNoNo
Kontsos et al. [33]Piezoelectric sensorCrack initiationNoNo
Xiang et al. [34]Piezoelectric sensorFatigue crackNoNo
Chen et al. [35]Eddy current sensorFatigue crackNoYes
Yuan et al. [36]Piezoelectric sensorImpactYesYes
This paperStrain sensorEarly plastic damageYesYes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, Q.; Huo, J.; Wang, Y.; Wang, X.; Wang, C. Structural Plastic Damage Warning and Real-Time Sensing System Based on Cointegration Theory. Sensors 2024, 24, 5961. https://doi.org/10.3390/s24185961

AMA Style

Gao Q, Huo J, Wang Y, Wang X, Wang C. Structural Plastic Damage Warning and Real-Time Sensing System Based on Cointegration Theory. Sensors. 2024; 24(18):5961. https://doi.org/10.3390/s24185961

Chicago/Turabian Style

Gao, Qiang, Junzhou Huo, Youfu Wang, Xiaotian Wang, and Chongru Wang. 2024. "Structural Plastic Damage Warning and Real-Time Sensing System Based on Cointegration Theory" Sensors 24, no. 18: 5961. https://doi.org/10.3390/s24185961

APA Style

Gao, Q., Huo, J., Wang, Y., Wang, X., & Wang, C. (2024). Structural Plastic Damage Warning and Real-Time Sensing System Based on Cointegration Theory. Sensors, 24(18), 5961. https://doi.org/10.3390/s24185961

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