Next Article in Journal
DBOS: A Dialog-Based Object Query System for Hospital Nurses
Next Article in Special Issue
Performance of Damage Identification Based on Directional Wavelet Transforms and Entopic Weights Using Experimental Shearographic Testing Results
Previous Article in Journal
Effects of Wearable Devices with Biofeedback on Biomechanical Performance of Running—A Systematic Review
Previous Article in Special Issue
Identification of Building Damage from UAV-Based Photogrammetric Point Clouds Using Supervoxel Segmentation and Latent Dirichlet Allocation Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Damage Identification and Quantification in Beams Using Wigner-Ville Distribution

Department of Fundamentals of Machinery Design, Faculty of Mechanical Engineering, Silesian University of Technology, 44-100 Gliwice, Poland
Sensors 2020, 20(22), 6638; https://doi.org/10.3390/s20226638
Submission received: 10 November 2020 / Revised: 16 November 2020 / Accepted: 18 November 2020 / Published: 19 November 2020
(This article belongs to the Special Issue Sensors for Structural Damage Identification)

Abstract

:
The paper presents the novel method of damage identification and quantification in beams using the Wigner-Ville distribution (WVD). The presented non-parametric method is characterized by high sensitivity to a local stiffness decrease due to the presence of damage, comparable with the sensitivity of the wavelet-based approaches, however the lack of selection of the parameters of the algorithm, like wavelet type and its order, and the possibility of reduction of the boundary effect make this method advantageous with respect to the mentioned wavelet-based approaches. Moreover, the direct relation between the energy density resulting from the application of WVD to modal rotations make it possible to quantify damage in terms of its width and depth. The results obtained for the numerical modal rotations of a beam presented in this paper, simulating the results of non-destructive testing achievable with the shearography non-destructive testing method, confirm high accuracy in localization of a damage as well as quantification of its dimensions. It was shown that the WVD-based method is suitable for detection of damage represented by the stiffness decrease of 1% and can be identified and quantified with a high precision. The presented results of quantification allowed extracting information on damage width and depth.

1. Introduction

Modern maintenance and monitoring demands require detection, identification, and quantification of structural damage in the possible earliest stage of its initiation and propagation. For this purpose, numerous methods that use primarily non-destructive testing (NDT) techniques were developed in the last decades, which found application in damage identification and structural health monitoring (SHM) in mechanical, civil, military, and aerospace engineering, where enabling of a possibility of early damage detection and identification is of a crucial importance [1,2,3,4,5,6,7,8]. Many of the NDT and SHM approaches used for the structural assessment, especially for the large structures, are based on the analysis of structural modal parameters due to the possibility of easy and inexpensive measurements and effectiveness in damage identification.
Numerous early approaches in modal-based damage identification were based on the analysis of changes in natural frequencies [9,10,11,12] or whole frequency response functions [13,14,15], however, when the modal curvaturefors started to be taken into consideration, the sensitivity to damage presence was satisfactory enough to detect and identify small damage [16]. Nevertheless, the attempts to improve the sensitivity of the damage identification methods in the last decades resulted in the development of numerous algorithms that basically relied on the post-processing of raw vibration data. Many of these approaches were based on beams with a damage, which, in spite of their simplicity, still have numerous applications in mechanical and civil engineering.
The developed algorithms used to post-process raw modal curvatures used numerous mathematical formulations. In [17], Ratcliffe used finite difference approximation of the Laplacian operator to detect and localize damage in beams, which allowed identification of very small damage represented by a thickness reduction of 0.5%. The authors of [1] introduced curvature damage factor (CDF) in order to consider damage in multiple modal curvatures and presented the results of identification of a damage both for theoretical formulation of a beam and a real bridge. Messina et al. [18] introduced the Multiple Damage Location Assurance Criterion (MDLAC) to detect and localize damage in beam and frame structures based on the changes in natural frequencies, which was further improved by the authors of [19] by using a genetic algorithm. Lu et al. [20] identified damage in beams by determining differences in modal curvatures of intact and damaged beams using CDF proposed in [1]. The authors of [21] proposed a two-step damage identification algorithm which uses a damage index of strain statistical moment to localize damage in beam-type structures. Another damage index for damage localization was proposed in [22], which was based on the normalized uniform load surface curvature of tested structures. Silva and Maia [23] proposed an estimation of damage location using the errors determined from constitutive models of beams and demonstrated successful localization of the simulated damage starting from 10% of thickness reduction. Brasiliano et al. [24] proposed a baseline approach of damage detection and localization based on the residual error method, which implies subtraction of displacements of mode shapes of a damaged structure from the displacements of an intact one. A similar approach was also used by the authors of [25,26,27,28]. The improved version of this approach was reported by the authors of [29,30], where the adaptive baseline concept was proposed. These approaches, however, did not guaranteed proper localization and identification of small damage, which is essential in demanding NDT and SHM applications.
Numerous further attempts were focused on the application of wavelet transforms (WTs) due to their high sensitivity to small disturbances in the modal curvatures. The continuous wavelet transform (CWT) was successfully used by Gentile and Messina [31] for detection and localization of open cracks in beams. This approach was confirmed by Douka et al. [32] in the similar study of identification of cracks in beams, which was additionally improved by proposing an intensity factor law for the quantification of the crack size. CWT was successfully used by Rucka [33,34] to localize cracks with up to 5% of thickness reduction of a tested beam. Han et al. [35] adapted the wavelet packet transform for identification of damage in beam-type structures based on the proposed energy rate index. The approach based on CWT was applied to identification of cracks in a bridge beam-type structure [36], which justified its sensitivity to damage in practical applications. A similar study was performed by Bayissa et al. [37], where the authors used CWT for damage identification based on numerical data and then applied the proposed approach for identification of cracks in bridge structures. Montanari et al. [38] investigated the influence of minimizing the number of the measurement points of the mode shapes on the damage detectability using CWT and proposed an optimal sampling procedure for effective damage identification. The authors of [39] proposed the damage identification algorithm in beams using the discrete wavelet transform (DWT) which allows for reducing computational complexity of the post-processing procedures. This concept was enriched by Ren and Sun [40] by taking into consideration the metrics related to the wavelet entropy in order to enhance the damage detectability. Khorram et al. [41] applied CWT for identification and quantification of structural damage, however, the estimation of the damage size using WT-based approaches always depends on an applied wavelet, thus, cannot be generalized. The problem of dependency of the resulting damage signatures on the applied wavelet was investigated by Janeliukstis et al. [42,43], where the authors tested numerous wavelets and introduced a catalog of wavelets according to the introduced metric of the damage estimate reliability indicating sensitivity to damage, which partially resolved the problem. In order to overcome the mentioned problem with selection of parameters of applied wavelets in WT-based post-processing algorithms, several other methods were developed. For instance, the authors of [44,45] applied the fractal dimension (FD) concept, but the sensitivity to small damage was still unsatisfactory, i.e., the best results reported in these studies were obtained for 10% of a thickness reduction. Moreover, the FD-based approach is more demanding in terms of required number of samples in the modal curvatures [46].
The recent studies of structural damage identification in beams considered more sophisticated algorithms focused on the improvement of the damage detection and identification abilities. Although the WT-based post-processing algorithms reveal higher effectiveness in damage identification, besides the necessity of selection of an appropriate wavelet, they suffer from significant influence to numerical or experimental noise in modal curvatures, which require additional processing of modal curvatures or improving the algorithms by adding additional steps based on artificial intelligence or soft computing methods [47]. Several algorithms that allowed for overcoming the problem of environmental noise were developed the authors of [48,49], where they proposed a hybrid approach based on Teager energy operator and CWT able to identify damage under very noisy conditions and validated this approach in experimental studies. In another study reported by the same team of the authors [50], yet another method was proposed, which was based on the complex-wavelet approach allowing identification of cracks under the noisy conditions, however, the smallest damage was of 20% thickness reduction and its performance for small damage is unknown. The other ways of the detectability improvement reported to-date include a combination of WT with the principal component analysis [51] or optimization algorithms [52,53], but, in general, they need additional tuning of parameters of the algorithms, which make them difficult to generalize and implement in practical problems.
One of the solutions being considered to resolve the deficiencies of WT-based algorithms with tuning is to find an appropriate algorithm for a post-processing procedure which can allow retaining similar sensitivity to structural damage as the mentioned ones but eliminate the data dependency during selection of parameters of such algorithms. Since WT is an integral transform, it is evident to look for similar non-parametric transforms. One of such transforms is the short-time Fourier transform (STFT), which is one of the simplest integral time-frequency transforms (note that in the context of the considered class of problems we will further consider such a type of transform as the space-frequency one since in the modal-based damage identification problems the input signals are not the time series, but the distance variable). In several studies [54,55], the application of STFT for damage identification was reported, however, the direct application of this transform to the detection and localization of damage is not effective due to the lack of possibility to obtain good resolution both in space and frequency domains, which was confirmed in some studies [37,40].
An alternative approach is to select an integral transform which allows retaining the lack of necessity of selection and tuning of parameters and simultaneously retaining good resolution in the space-frequency domain. One of such a type of transform is the Wigner-Ville distribution (WVD), being the fundamental integral non-parametric transform from the Cohen’s class, similarly as STFT [56]. In spite of STFT, WVD is the quadratic transform with an optimal resolution and instantaneous power density spectrum in both considered domains [57]. Due to these properties, WVD found wide applicability primarily in the problems of technical diagnostics of rotary machinery and their condition monitoring [56,57,58,59,60], and in detecting defects in various applications [61,62,63]. To the best of the author’s knowledge, there were only several applications of WVD to the structural damage identification problems reported to date. Gillich and Praisach in their study [64] used WVD only for extraction the information on natural frequencies from the vibration signals, while Dai and He [65] applied WVD for the determination of damage position in the NDT technique based on ultrasonic guided waves.
The aim of this paper is to investigate the capability of the newly presented non-parametric damage identification algorithm based on WVD as well as to evaluate its performance in terms of the accuracy of damage identification and quantification of the damage based on the modal rotations obtained from numerical simulations. The simulations were prepared in such a way to simulate high-resolution data achievable from the shearography NDT method. The analysis of modal rotations is superior with respect to widely used modal curvatures as the input data, since it reveals higher sensitivity to damage, which was shown in several studies (see e.g., [24,66]), including the previous studies by the author of this paper [67,68]. In this study, this superiority will be shown for the newly presented post-processing algorithm. Moreover, the theoretical investigation on superiority of the WVD-based approach over STFT and WT is provided in the paper and the validation is performed by presenting the results of the case studies with the simulated damage in a beam.

2. Theoretical Background

The simplest transform used for the space-frequency analysis is the mentioned STFT proposed by Gabor [69], represented by the following integral equation:
S ( τ , ω ) = s ( x ) w ( x τ ) e i ω x d x ,
where s ( x ) is the transformed signal in the space domain, w ( τ ) is the window function (usually the Gaussian window is assumed), τ is the space shifting factor of the window w , ω is the frequency, and S ( τ , ω ) is the complex function representing the result of transform in the form of phase and magnitude in the space-frequency domain.
As it was already mentioned, STFT operates with a fixed window in the space-frequency domain and changing the widow size in one domain results in a corresponding change in another domain. This property is connected with the uncertainty principle
Δ x Δ ω 1 2 ,
where Δ x and Δ ω are the bounds limiting the window in space and frequency domains, respectively.
Due to this limitation, STFT cannot provide appropriate resolution of S ( τ , ω ) in the space-frequency domain, thus cannot be effectively used for damage identification problems. To overcome this limitation, numerous researchers used CWT proposed in its present form by Goupillaud, Grossman and Morlet [70], and given by the following integral equation:
W ( a , b ) = w ( a ) s ( x ) ψ ( x b a ) d x ,
where a and b are the scaling and translation parameters, w ( a ) is the weighting function usually set to 1 / a for energy conservation purposes, and ψ ( x b a ) is the wavelet kernel, which defines the basis function of the transform.
Due to the possibilities of scaling and translation of the basis function, CWT has the property of the multiresolution analysis, which allows representing of W ( a , b ) in the space-scale domain in such a way that all the components of s ( x ) are well-localized in this domain, which makes it possible to effectively perform damage identification (proven in the numerous studies cited in Section 1). The main deficiency of application of CWT in such a class of problems is the necessity of selection of the basis function ψ , which has a crucial influence on the resulting representation of W ( a , b ) , and thus, on the effectiveness of damage identification. Besides the question of the variable damage identification effectiveness, the difficulties occur when CWT results are used for damage quantification, since the magnitudes of the resulting wavelet coefficients depend on a selected basis function.
To solve this problem, it is proposed to apply WVD introduced by Ville [71] for applications in the problems related analysis of electrical signals. WVD is given by the following integral equation:
W V ( τ , ω ) = f ( τ + x 2 ) f * ( τ x 2 ) e 2 i π ω x d x ,
where f ( τ + x 2 ) f * ( τ x 2 ) is the instantaneous auto-correlation function, and W V ( τ , ω ) is the Wigner-Ville energy density spectrum of a transformed signal.
As one can observe, WVD is non-parametric, but retains most of properties useful from the point of view of damage identification applications. Moreover, as it can be observed from Equations (1) and (3), STFT and CWT has a similar origin and the difference between these transforms lies in the basis function, while the transformed signal is given in both cases in the explicit form. In contrast to this, in WVD the signal is not given in an explicit form, but is represented by the auto-correlation function, and enters twice into the integration operation, therefore WVD is a bi-linear transform [72] which belongs to the Cohen’s class of integral transforms. The WVD holds the properties of scaling and translation, similarly to CWT, which makes it useful for the numerous applications, including the damage identification problems investigated in this paper. These operations, although different than in CWT, where they are controlled by the parameters a and b, can be related to each other in both transforms (see [73] for more details). In spite of STFT and CWT, WVD is a non-linear (quadratic) transform which has the consequences in the resulting energy density spectra, namely, the appearance of cross-terms, which, in the case of damage identification problems, is the undesired property, since the cross-terms may mask the identified damage signatures. The WVD is directly connected with another distribution from the Cohen’s class—the Choi-Williams distribution (CWD) —given by the following equation:
C W ( τ , ω ) = e u 2 4 x 2 / σ 4 π x 2 / σ f ( τ u + x 2 ) f * ( τ u x 2 ) e 2 i π ω x d u   d x ,
where C W ( τ , ω ) is the Choi-Williams energy density spectrum of a transformed signal, σ is the temporal scaling parameter, and u is the ambiguity domain parameter. From (5), one can extract the kernel function of the CWD in the following form:
Φ C W ( x , u ) = e ( 2 π x u ) 2 σ .
For the case when σ , Φ C W ( x , u ) 1 , and (5) reduces to (4) [72].
Numerous comparative studies on the time-frequency representations obtained using the mentioned transforms can be found in [72,73]. The problem with cross-terms is also presented in the earlier works of one of the authors [74,75]. Nevertheless, the properties of WVD, including primarily the lack of parameters to tune and direct connection with the energy of the transformed signal, make it a good candidate for damage identification and quantification, which is the subject of the author’s studies presented in the next sections.

3. Analyzed Structure and Damage Scenarios

In this study, the vibration of the Euler-Bernoulli beam was simulated with the finite element analysis to obtain the first three bending natural modes of vibration according to the studies reported in [76]. The aluminum beam model with the length of 400 mm, the width of 40 mm, and the thickness of 3 mm, and the following material properties: Young’s modulus of 67.8 GPa and the density of 2700 kg/m3, was discretized in 2263 equally spaced finite elements. The definition of such a dense finite element mesh is related to a possibility of obtaining high-resolution data using shearography, the optical NDT technique used by the author in damage identification studies in the previous works (see e.g., [77,78]). Each element was of 2 nodes, having two displacement and two rotation values (see Figure 1). The modal displacements and rotations were determined for the intact beam as well as for the damaged beam, considering the first damage in the form of a slot of the width of 5 mm and located at the distance of 284 mm and the second damage in the form of a slot of the width of 3 mm and located at the distance of 200 mm from the beam end (see Figure 2). For the damage identification study, eight damage scenarios were defined by combining the thickness reduction of the two slots, the respective values of which are shown in Table 1. Both displacements and rotations were collected as the results of simulation (see Figure 3).
In order to perform the damage quantification study, the same numerical model was used. A single damage in the form of slot with the width of 5 mm, located at the distance of 284 mm from the end of the beam, was simulated within seven scenarios with variable depth starting from 1% and ending with 20% of the thickness reduction, which are presented in Table 2.

4. Results of Damage Identification

In this study, the investigation of an application of WVD for damage identification and quantification was performed. The selection of WVD among the other available TFDs is due to the fact that WVD is the most energy-concentrated distribution within the distributions of the Cohen’s class [79]. According to the results of the initial studies, the application of CWD presented in Section 2 does not result in satisfactory results for small values of the temporal scaling parameter σ , typically used for calculations, due to much less sensitivity to the simulated damage as well as due to the appearance of cross-terms.
Both baseline-free and baseline approaches were presented in order to determine the sensitivity of the proposed damage identification method to simulated damage. Then, the analysis of sensitivity to damage when the displacements and rotations of a tested beam are taken into consideration as an input to the proposed processing procedures was investigated to show the superiority of the chosen modal rotations.
Although WVD is the non-parametric transform, as was already mentioned before, the influence of a frequency range onto the results of identified damage was observed. The increasing frequency range influencing on the increase of the width of the damage signature band in the resulting energy density distributions as well as on the appearance and widening of the side bands is known in the literature as the boundary effect. Several examples of the damage identification results for the damage scenario 3I and mode 1 are presented in Figure 4. To increase the visual distinguishability of the damage signatures the obtained energy values were raised to the power of 0.01. This value was selected based on empirical analysis, and it was found that the distinguishability of the damage signatures improved asymptotically with the decrease of a power. The assumed value seems to be appropriate, since its further decreasing does not affect any visual improvement of the damage signatures.
From the results presented in Figure 4, it is clearly visible that for the frequency ranges of 0–20 Hz and higher the damage localization and identification are not possible due to blurring the damage signature, while for the range of 0–10 Hz one can precisely localize the simulated damage and even identify its boundaries. Moreover, the boundary effect for this range of frequency does not appear, which is a great advantage of this approach with respect to WT-based damage identification procedures, where this effect can mask damage signatures near the boundaries as well as significantly decrease distinguishability of the detected damage, as it was reported in numerous studies [31,33,43,80]. One can observe that in the case of TFDs the frequency range is equivalent to the length of the support of the applied wavelet in the WT-based approaches in the light of the increasing the width of the distortion on the boundaries due to the boundary effect. Considering the above observations, the range of 0–10 Hz was assumed in the next studies.

4.1. Baseline-Free Approach

The damage scenarios presented in Table 1 were used to investigate the detectability of the simulated damage in the form of slots for their various depths and locations. The obtained results for the particular damage scenarios and modes are presented in Figure 5.
From the obtained results, it can be noticed that only the second damage in the damage scenario 5I was not detected, while in the rest of the considered damage scenarios all the damage sites were well detectable and properly localized. It can be observed that the magnitude of energy density on the results presented in Figure 5 varies depending on the analyzed mode, which is especially well visible for the damage scenarios 6I–8I, where the damage signature for the second damage (located in the middle of the beam) is barely detectable for the second mode. This variability proves the proportional dependency of the resulting energy with the magnitude of vibration in the location of damage.

4.2. Baseline Approaches

In order to examine the possibilities of improvement of the ability of damage identification using the proposed WVD-based method, two baseline approaches were considered in this study. In general, the baseline approach was considered using of the modes of the intact beam as the reference data, however, the subtraction of the modal data for the damaged beam from the reference data can be performed in two ways: subtraction of mode shapes and then transforming the resulting differences and transforming the mode shapes for the intact and damaged beams separately and subtracting the resulting energy density spectra. According to the theoretical fundamentals of WVD [72], these operations are not interchangeable due to the appearance of the additional terms during subtraction of modes for the intact and damaged beams, thus the results of these baseline approaches will be different.

4.2.1. Baseline Approach Based on the Transform of Differences

The study was performed based on the results of subtraction of the mode shapes from the damage scenarios presented in Table 1 and the mode shapes on the intact beam. The results of the performed calculations are presented in Figure 6.
As it can be observed, the applied baseline-free approach makes it possible to identify all the simulated damage, including the damage scenario 5I, where the second damage became detectable. Moreover, the observed damage signatures are characterized by significantly higher distinguishability compared to the results of the baseline-free approach. However, in contrast to the results obtained using the baseline-free approach, the damage signatures observable in the obtained results using the presented baseline approach are blurred, which makes it impossible to identify the boundaries of the simulated slots.

4.2.2. Baseline Approach Based on the Differences of the Transformed Modes

In this approach, the modes were subjected to transform using the proposed WVD-based method separately for the intact and damaged beams using the data for the defined damage scenarios presented in Table 1. The results of the calculations are presented in Figure 7.
Analyzing the results presented in Figure 7, one can notice that all the damage signatures in all the considered scenarios are well detectable. Comparing these results with the results obtained using the baseline approach based on transformation of mode differences, it can be stated that the performance of the following approach is better, since the distinguishability of the second damage, especially for the second mode in particular damage scenarios, is much better when the latter approach was used. Moreover, it can be observed that the boundaries of the simulated slots are well identifiable in the case of using the second baseline approach, which were undetectable in the most of damage scenarios when the first baseline approach was used. The superiority of the second baseline approach over the first one is justified theoretically, due to the lack of additional disturbances caused by transforming the differences of modes, which was discussed in the preamble to Section 4.2.

4.3. Displacements vs. Rotations

In order to justify the selection of data for this study, namely, that modal rotations were selected instead of modal displacement, an additional study was performed. This justification was performed for the damage scenario 8I, where both simulated slots are well detectable with using modal rotations. The tests were performed for the considered baseline and baseline-free approaches, where the baseline approach based on the subtraction of the mode shapes before the transform was marked as “Baseline 1”, while the approach based on subtraction of the transformed mode shapes was marked as “Baseline 2”. The results for the modal displacements are presented in Figure 8.
It can be observed that, in general, the damage distinguishability is much lower compared to the results obtained for the modal rotations (cf. Figure 5, Figure 6 and Figure 7). The baseline-free approach, similar to the case of the modal rotations used as the input data, is characterized by the lowest sensitivity to damage. However, in contrast to the results obtained using modal rotations (see Figure 5), the simulated damage remains undetectable when the modal displacements were used. In the case of the baseline approaches, the simulated damage become to be barely detectable, excluding the results obtained for the second mode (both slots are undetectable for the “Baseline 1” approach, and the second slot is undetectable for the “Baseline 2” approach). Moreover, for both baseline approaches used for the third mode, additional disturbances appeared at the distance of 262 mm in the energy density spectra presented in Figure 8. These disturbances are probably the result of the application of WVD algorithm and can be considered as the cross-terms (see Section 2 for more details), which is undesirable due to the false damage signature appearance. Comparing the results obtained for modal displacements and rotations, one can conclude that the modal rotations provide better distinguishability of damage, since it is possible to obtain the derivatives of the mode shapes in a direct way, which is confirmed in the numerous independent studies (see e.g., [24,81]).

5. Results of Damage Quantification

The quantification of damage was performed considering two dimensions of the simulated slots: their width and depth. From the previously presented results, it is possible to observe that for all the approaches considered in this study the difference in the width of the simulated slots is noticeable. In order to validate the width of the considered slots quantitatively, the analysis of the obtained results was performed. The exemplary magnifications of the damage signatures of the resulting energy density spectra obtained for the damage scenario 8I (mode 1) for all the considered approaches is presented in Figure 9. The widths of the identified slots were determined separately for the particular modes and damage scenarios and it was found that the determined width does not depend on the considered mode, neither on the damage scenario (i.e., depth of the slot). The results show that in all cases the identified damage signature is overestimated with respect to the true damage width, and this overestimation is of 1.53 mm and 0.71 mm for the slots 1 and 2, respectively. This implies the conclusion that the quantification error in this case depends on the damage width. However, the independency of the estimated width on the considered mode and the depth of the damage makes the proposed WVD-based method advantageous for further applications.
In order to quantify the damage depth, another dataset with elevated damage depths according to the scenarios presented in Table 2 was analyzed. The quantification was performed using all the approaches presented in this paper, namely, the baseline-free and two baseline ones. Since small variabilities in the values within the damage signatures are observable, both on the width of the damage signatures (due to the differences in energy density values for the detected boundary of a slot and the interior of the slot) and their length (due to the influence of the frequency), it is essential to consider that the energy density values corresponded with the boundaries of the detected damage. For this purpose, the values of the energy density from the locations of detected boundaries of the damage were averaged for a given frequency, and then the maximum value was taken from the collected values. The results of the depth quantification are presented in Table 3.
The obtained values of the energy density for the investigated damage scenarios were subjected to the regression procedure in such a way that for particular approaches and particular modes the energy density values were plotted in a function on the increasing damage depth. This operation allowed to confirm the proportionality of the increase of the determined energy density values with respect to the considered damage depths. The cubic polynomial function was used for the regression due to the best fitting performance and the connection with the character of influence of a thickness reduction on the resulting structural stiffness. The estimation of the crack depth relationship with the local stiffness reduction based on polynomials is widely used in the literature (see e.g., [82,83,84,85]). The exemplary plot and the results of regression are presented in Figure 10. Note that the energy density for the intact beam was determined in the same way as for the damaged beams.
It was observed that the regression performance of the applied function allowed to obtain the exact relationship between the damage depth and the amount of the energy density. For all the considered approaches and all the considered modes, the coefficient of determination was of R2 = 1, which can be explained by the numerical character of the considered data. The obtained relationships allow for concluding that there exists a direct relation between the damage depth and the resulting energy density in the energy density spectra obtained using the WVD-based method. Following this, it is possible to quantify the damage depth based on comparison of the energy density values from the damage signature and damage-free regions. However, as it can be observed from the results presented in Table 3, the obtained energy density values for damage signatures are mode-dependent, which is a result of the various magnitudes of modal rotations for the particular modes. This means, that the operation of the quantification of the damage depth need to be performed for each considered mode separately.

6. Conclusions

The new method of damage identification and quantification based on the Wigner-Ville distribution has been presented to provide reliable information about location and size of single and multiple damage sites in beams. The presented algorithm is non-parametric, which means that the step of selection the parameters of processing typical for the WT-based approaches (such as selection of the wavelet function and its order) is omitted with retaining the sensitivity to damage comparable with the WT-based approach. Moreover, the proposed WVD-based method does not produce the boundary effect typical for WT-based approaches, and thus, allows identifying and quantifying damage also in the boundary regions. Three damage identification approaches were presented in this study: the baseline-free approach and two baseline approaches, which use the information of a condition of an intact structure. All the considered approaches reveal excellent sensitivity to damage, however the baseline approaches are characterized by considerably better distinguishability of the damage with respect to the baseline-free one, which is obvious due to the damage emphasizing by subtraction of the modes of damaged and intact structures. It was shown that the proposed method reveals high sensitivity to damage, i.e., for the baseline-free approach the detectability limit was of 0.1 mm of a thickness reduction, and for the baseline methods this limit was of 0.028 mm of a thickness reduction, which corresponds to 3.33% and 0.93%, respectively. It was also proven that using the modal rotations provides much better distinguishability of damage with respect to modal displacements, which was confirmed for all three considered approaches, both baseline and baseline-free. The performed damage identification studies proved that it is possible to identify the boundaries of a damage, while the results of the quantification studies show that the dimensions of a damage, namely its width and depth, can be accurately determined using the WVD-based method. The estimated damage width was insignificantly higher than the true width (1.53 mm for the 5 mm damage and 0.71 mm for the 3 mm damage), which seems to be acceptable from the point of view of practical applications, while the estimated damage depth was exactly the same as the true damage depth, which was proven by the exact matching regressions of the damage depth vs. energy density for all considered cases.
The above-presented advantages of the proposed WVD-based method creates a possibility of its application in practical problems and possible hardware implementation for the structural damage identification and SHM purposes. The performed studies were based on the data obtained from the numerical modal analysis of a beam in the form of modal rotations, and the influence of the measurement noise need to be investigated in order to evaluate the performance of this method and its validity for the application in practical problems. Nevertheless, the obtained results create a promising research direction for development of reliable non-parametric tools for NDT and SHM problems in mechanical and civil engineering. The proposed method can be applied to the experimental results obtained using the shearography NDT method, which among the high resolution and sensitivity to local structural changes can provide measurements of modal rotations in an explicit way.

Funding

This study was funded by the Polish National Agency for Academic Exchange within the Bekker programme under the agreement no. PPN/BEK/2019/1/00048.

Acknowledgments

The author would like to express his gratitude to Hernâni Lopes from the Polytechnic Institute of Porto for sharing numerical data for performing the studies and Cristóvão Manuel Mota Soares from the University of Lisbon for supervision of the research internship financed by the Polish National Agency for Academic Exchange.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Abdel Wahab, M.M.; De Roeck, G. Damage detection in bridges using modal curvatures: Application to a real damage scenario. J. Sound Vib. 1999, 226, 217–235. [Google Scholar] [CrossRef]
  2. Mehnen, J.; Tinsley, L.; Roy, R. Automated in-service damage identification. CIRP Ann. 2014, 63, 33–36. [Google Scholar] [CrossRef]
  3. Katunin, A.; Dragan, K.; Dziendzikowski, M. Damage identification in aircraft composite structures: A case study using various non-destructive testing techniques. Compos. Struct. 2015, 127, 1–9. [Google Scholar] [CrossRef]
  4. Rezaei, M.M.; Behzad, M.; Moradi, H.; Haddadpour, H. Modal-based damage identification for the nonlinear model of modern wind turbine blade. Renew. Energy 2016, 94, 391–409. [Google Scholar] [CrossRef]
  5. Xu, Y.-L.; Zhang, C.-D.; Zhan, S.; Spencer, B.F. Multi-level damage identification of a bridge structure: A combined numerical and experimental investigation. Eng. Struct. 2018, 156, 53–67. [Google Scholar] [CrossRef]
  6. Wu, B.; Wu, G.; Yang, C.; He, Y. Damage identification method for continuous girder bridges based on spatially-distributed long-gauge strain sensing under moving load. Mech. Syst. Signal Process. 2018, 104, 415–435. [Google Scholar] [CrossRef]
  7. Gong, Y.; Shao, H.; Luo, J.; Li, Z. A deep transfer learning model for inclusion defect detection of aeronautics composite materials. Compos. Struct. 2020, 252, 112681. [Google Scholar] [CrossRef]
  8. Hou, R.; Xia, Y. Review on the new development of vibration-based damage identification for civil engineering structures: 2010–2019. J. Sound Vib. 2020, 491, 115741, in press. [Google Scholar]
  9. Lee, Y.-S.; Chung, M.-J. A study on crack detection using eigenfrequency test data. Comput. Struct. 2000, 77, 327–342. [Google Scholar] [CrossRef]
  10. Kim, J.-T.; Stubbs, N. Crack detection in beam-type structures using frequency data. J. Sound Vib. 2003, 259, 145–160. [Google Scholar] [CrossRef]
  11. Xiang, J.; Liang, M.; He, Y. Experimental investigation of frequency-based multi-damage detection for beams using support vector regression. Eng. Fract. Mech. 2014, 131, 257–268. [Google Scholar] [CrossRef]
  12. Altunışık, A.C.; Okur, F.Y.; Kahya, V. Structural identification of a cantilever beam with multiple cracks: Modeling and validation. Int. J. Mech. Sci. 2017, 130, 74–89. [Google Scholar] [CrossRef]
  13. Sampaio, R.P.C.; Maia, N.M.M.; Silva, J.M.M. Damage detection using the frequency-response-function curvature methods. J. Sound Vib. 1999, 226, 1029–1042. [Google Scholar] [CrossRef]
  14. Araújo dos Santos, J.V.; Mota Soares, C.M.; Mota Soares, C.A.; Maia, N.M.M. Structural damage identification in laminated structures using FRF data. Compos. Struct. 2005, 67, 239–249. [Google Scholar] [CrossRef]
  15. Dilena, M.; Limongelli, M.P.; Morassi, A. Damage localization in bridges via the FRF interpolation method. Mech. Syst. Signal Process. 2015, 52–53, 162–180. [Google Scholar] [CrossRef]
  16. Pandey, A.K.; Biswas, M.; Samman, M.M. Damage detection from changes in curvature mode shapes. J. Sound Vib. 1991, 145, 321–332. [Google Scholar] [CrossRef]
  17. Ratcliffe, C.P. Damage detection using a modified Laplacian operator on mode shape data. J. Sound Vib. 1997, 204, 505–517. [Google Scholar] [CrossRef]
  18. Messina, A.; Williams, E.J.; Contursi, T. Structural damage detection by a sensitivity and statistical-based method. J. Sound Vib. 1998, 216, 791–808. [Google Scholar] [CrossRef]
  19. Gomes, H.M.; Silva, N.R.S. Some comparisons for damage detection on structures using genetic algorithms and modal sensitivity method. Appl. Math. Model. 2008, 32, 2216–2232. [Google Scholar] [CrossRef]
  20. Lu, X.B.; Liu, J.K.; Lu, Z.R. A two-step approach for crack identification in beam. J. Sound Vib. 2013, 332, 282–293. [Google Scholar] [CrossRef]
  21. Wang, D.; Xiang, W.; Zhu, H. Damage identification in beam type structures based on statistical moment using a two step method. J. Sound Vib. 2014, 333, 745–760. [Google Scholar] [CrossRef]
  22. Sung, S.H.; Jung, H.J.; Jung, H.Y. Damage detection for beam-like structures using the normalized curvature of a uniform load surface. J. Sound Vib. 2013, 332, 1501–1519. [Google Scholar] [CrossRef]
  23. Silva, T.A.N.; Maia, N.M.M. Detection and localisation of structural damage based on the error in the constitutive relations in dynamics. Appl. Math. Model. 2017, 46, 736–749. [Google Scholar] [CrossRef]
  24. Brasiliano, A.; Doz, G.N.; de Brito, J.L.V. Damage identification in continuous beams and frame structures using the Residual Error Method in the Movement Equation. Nucl. Eng. Des. 2004, 227, 1–17. [Google Scholar] [CrossRef]
  25. Lestari, W.; Qiao, P. Damage detection of fiber-reinforced polymer honeycomb sandwich beams. Compos. Struct. 2005, 67, 365–373. [Google Scholar] [CrossRef]
  26. Ciambella, J.; Vestroni, F. The use of modal curvatures for damage localization in beam-type structures. J. Sound Vib. 2015, 340, 126–137. [Google Scholar] [CrossRef]
  27. He, M.; Yang, T.; Du, Y. Nondestructive identification of composite beams damage based on the curvature mode difference. Compos. Struct. 2017, 176, 178–186. [Google Scholar] [CrossRef]
  28. Gorgin, R. Damage identification technique based on mode shape analysis of beam structures. Structures 2020, 27, 2300–2308. [Google Scholar] [CrossRef]
  29. Colombo, L.; Sbarufatti, C.; Giglio, M. Definition of a load adaptive baseline by inverse finite element method for structural damage identification. Mech. Syst. Signal Process. 2019, 120, 584–607. [Google Scholar] [CrossRef]
  30. Colombo, L.; Oboe, D.; Sbarufatti, C.; Cadini, F.; Russo, S.; Giglio, M. Shape sensing and damage identification with iFEM on a composite structure subjected to impact damage and non-trivial boundary conditions. Mech. Syst. Signal Process. 2021, 148, 107163. [Google Scholar] [CrossRef]
  31. Gentile, A.; Messina, A. On the continuous wavelet transforms applied to discrete vibrational data for detecting open cracks in damaged beams. Int. J. Solids Struct. 2003, 40, 295–315. [Google Scholar] [CrossRef]
  32. Douka, E.; Loutridis, S.; Trochidis, A. Crack identification in beams using wavelet analysis. Int. J. Solids Struct. 2003, 40, 3557–3569. [Google Scholar] [CrossRef]
  33. Rucka, M.; Wilde, K. Application of continuous wavelet transform in vibration based damage detection method for beams and plates. J. Sound Vib. 2006, 297, 536–550. [Google Scholar] [CrossRef]
  34. Rucka, M. Damage detection in beams using wavelet transform on higher vibration modes. J. Theor. Appl. Mech. 2011, 49, 399–417. [Google Scholar]
  35. Han, J.-G.; Ren, W.-X.; Sun, Z.-S. Wavelet packet based damage identification of beam structures. Int. J. Solids Struct. 2005, 42, 6610–6627. [Google Scholar] [CrossRef] [Green Version]
  36. Zhu, X.Q.; Law, S.S. Wavelet-based crack identification of bridge beam from operational deflection time history. Int. J. Solids Struct. 2006, 43, 2299–2317. [Google Scholar] [CrossRef] [Green Version]
  37. Bayissa, W.L.; Haritos, N.; Thelandersson, S. Vibration-based structural damage identification using wavelet transform. Mech. Syst. Signal Process. 2008, 22, 1194–1215. [Google Scholar] [CrossRef]
  38. Montanari, L.; Spagnoli, A.; Basu, B.; Broderick, B. On the effect of spatial sampling in damage detection of cracked beams by continuous wavelet transform. J. Sound Vib. 2015, 345, 233–249. [Google Scholar] [CrossRef]
  39. Grabowska, J.; Palacz, M.; Krawczuk, M. Damage identification by wavelet analysis. Mech. Syst. Signal Process. 2008, 22, 1623–1635. [Google Scholar] [CrossRef]
  40. Ren, W.-X.; Sun, Z.-S. Structural damage identification by using wavelet entropy. Eng. Struct. 2008, 30, 2840–2849. [Google Scholar] [CrossRef]
  41. Khorram, A.; Bakhtiari-Nejad, F.; Rezaeian, M. Comparison studies between two wavelet based crack detection methods of a beam subjected to a moving load. Int. J. Eng. Sci. 2012, 51, 204–215. [Google Scholar] [CrossRef]
  42. Janeliukstis, R.; Rucevskis, S.; Wesolowski, M.; Chate, A. Multiple damage identification in beam structure based on wavelet transform. Procedia Eng. 2017, 172, 426–432. [Google Scholar] [CrossRef]
  43. Janeliukstis, R.; Rucevskis, S.; Wesolowski, M.; Chate, A. Experimental structural damage localization in beam structure using spatial continuous wavelet transform and mode shape curvature methods. Measurement 2017, 102, 253–270. [Google Scholar] [CrossRef]
  44. Qiao, P.; Cao, M. Waveform fractal dimension for mode shape-based damage identification of beam-type structures. Int. J. Solids Struct. 2008, 45, 5946–5961. [Google Scholar] [CrossRef] [Green Version]
  45. Radzieński, M.; Krawczuk, M.; Palacz, M. Improvement of damage detection methods based on experimental modal parameters. Mech. Syst. Signal Process. 2011, 25, 2169–2190. [Google Scholar] [CrossRef]
  46. Dessi, D.; Camerlengo, G. Damage identification techniques via modal curvature analysis: Overview and comparison. Mech. Syst. Signal Process. 2015, 52–53, 181–205. [Google Scholar] [CrossRef]
  47. Fan, W.; Qiao, P. Vibration-based damage identification methods: A review and comparative study. Struct. Health Monit. 2011, 10, 83–111. [Google Scholar] [CrossRef]
  48. Cao, M.; Xu, W.; Ostachowicz, W.; Su, Z. Damage identification for beams in noisy conditions based on Teager energy operator-wavelet transform modal curvature. J. Sound Vib. 2014, 333, 1543–1553. [Google Scholar] [CrossRef]
  49. Cao, M.; Radzieński, M.; Xu, W.; Ostachowicz, W. Identification of multiple damage in beams based on robust curvature mode shapes. Mech. Syst. Signal Process. 2014, 46, 468–480. [Google Scholar] [CrossRef]
  50. Cao, M.-S.; Xu, W.; Ren, W.-X.; Ostachowicz, W.; Sha, G.-G.; Pan, L.-X. A concept of complex-wavelet modal curvature for detecting multiple cracks in beams under noisy conditions. Mech. Syst. Signal Process. 2016, 76–77, 555–575. [Google Scholar] [CrossRef]
  51. Shahsavari, V.; Chouinard, L.; Bastien, J. Wavelet-based analysis of mode shapes for statistical detection and localization of damage in beams using likelihood ratio test. Eng. Struct. 2017, 132, 494–507. [Google Scholar] [CrossRef]
  52. Katunin, A.; Przystałka, P. Detection and localization of delaminations in composite beams using fractional B-spline wavelets with optimized parameters. Eksploat. Niezawodn. 2014, 15, 391–399. [Google Scholar]
  53. Jiang, Y.; Wang, N.; Zhong, Y. A two-step damage quantitative identification method for beam structures. Measurement 2021, 168, 108434. [Google Scholar] [CrossRef]
  54. Xu, Y.F.; Chen, D.-M.; Zhu, W.D. Damage identification of beam structures using free response shapes obtained by use of a continuously scanning laser Doppler vibrometer system. Mech. Syst. Signal Process. 2017, 92, 226–247. [Google Scholar] [CrossRef] [Green Version]
  55. Xu, Y.F.; Chen, D.-M.; Zhu, W.D. Modal parameter estimation using free response measured by a continuously scanning laser Doppler vibrometer system with application to structural damage identification. J. Sound Vib. 2020, 485, 115536. [Google Scholar] [CrossRef]
  56. Padovese, L.R. Hybrid time–frequency methods for non-stationary mechanical signal analysis. Mech. Syst. Signal Process. 2004, 18, 1047–1064. [Google Scholar] [CrossRef]
  57. Wu, J.-D.; Huang, C.-K. An engine fault diagnosis system using intake manifold pressure signal and Wigner-Ville distribution technique. Expert Syst. Appl. 2011, 38, 536–544. [Google Scholar] [CrossRef]
  58. Zou, J.; Chen, J. A comparative study on time–frequency feature of cracked rotor by Wigner-Ville distribution and wavelet transform. J. Sound Vib. 2004, 276, 1–11. [Google Scholar] [CrossRef]
  59. Kim, B.S.; Lee, S.H.; Lee, M.G.; Ni, J.; Song, J.Y.; Lee, C.W. A comparative study on damage detection in speed-up and coast-down process of grinding spindle-typed rotor-bearing system. J. Mater. Proces. Technol. 2007, 187–188, 30–36. [Google Scholar] [CrossRef]
  60. Climente-Alarcon, V.; Antonino-Daviu, J.A.; Riera-Guasp, M.; Puche-Panadero, R.; Escobar, L. Application of the Wigner-Ville distribution for the detection of rotor asymmetries and eccentricity through high-order harmonics. Electr. Power Syst. Res. 2012, 91, 28–36. [Google Scholar] [CrossRef]
  61. Rajshekhar, G.; Gorthi, S.S.; Rastogi, P. Detection of defects from fringe patterns using a pseudo-Wigner-Ville distribution based method. Opt. Laser. Eng. 2012, 50, 1059–1062. [Google Scholar] [CrossRef]
  62. Zeng, W.; Wang, H.; Tian, G.; Hu, G. Application of laser ultrasound imaging technology in the frequency domain based on Wigner-Ville algorithm for detecting defect. Opt. Laser Technol. 2015, 74, 72–78. [Google Scholar] [CrossRef]
  63. Ramapriya, D.M.; Gradoni, G.; Creagh, S.C.; Tanner, G.; Moers, E.; Arteaga, I.A. Nearfield acoustical holography—A Wigner function approach. J. Sound Vib. 2020, 486, 115593. [Google Scholar] [CrossRef]
  64. Gillich, G.-R.; Praisach, Z.-I. Modal identification and damage detection in beam-like structures using the power spectrum and time-frequency analysis. Signal Process. 2014, 96, 29–44. [Google Scholar] [CrossRef]
  65. Dai, D.; He, Q. Structure damage localization with ultrasonic guided waves based on a time-frequency method. Signal Process. 2014, 96, 21–28. [Google Scholar] [CrossRef]
  66. Spagnoli, A.; Montanari, L.; Basu, B.; Broderick, B. Nonlinear damage identification in fiber-reinforced cracked composite beams through time-space wavelet analysis. Procedia Mater. Sci. 2014, 3, 1579–1584. [Google Scholar] [CrossRef] [Green Version]
  67. Reis Lopes, H.M.; Araújo dos Santos, J.V.; Mota Soares, C.M.; Miranda Guedes, R.J.; Pires Vaz, M.A. A numerical-experimental method for damage location based on rotation fields spatial differentiation. Comput. Struct. 2011, 89, 1754–1770. [Google Scholar] [CrossRef]
  68. Lopes, H.; Ribeiro, J. Structural health monitoring in composite automotive elements. In New Advances in Vehicular Technology and Automotive Engineering; Carmo, J.P., Ribeiro, J.E., Eds.; InTechOpen: London, UK, 2012; pp. 285–302. [Google Scholar]
  69. Gabor, D. Theory of communication. J. Inst. Electr. Eng. 1946, 93, 429–457. [Google Scholar] [CrossRef]
  70. Goupillaud, P.; Grossman, A.; Morlet, J. Cycle-octave and related transforms in seismic signal analysis. Geoexploration 1984, 23, 85–102. [Google Scholar] [CrossRef]
  71. Ville, J. Théorie et applications de la notion de signal analytique. Cables Transm. 1948, 2, 61–74. [Google Scholar]
  72. Cohen, L. Time-Frequency Analysis; Prentice Hall PTR: Upper Saddle River, NJ, USA, 1995. [Google Scholar]
  73. Boashash, B. (Ed.) Time-Frequency Signal Analysis and Processing, 2nd ed.; Elsevier: London, UK, 2016. [Google Scholar]
  74. Katunin, A. Application of time-frequency distributions in diagnostic signal processing problems: A case study. Diagnostyka 2016, 17, 95–103. [Google Scholar]
  75. Katunin, A. Localization of damage in beam-like structures applying time-frequency distributions to modal shapes of vibration. Diagnostyka 2016, 17, 53–58. [Google Scholar]
  76. Minnini, M.; Gabriele, S.; Lopes, H.; Araújo dos Santos, J.V. Damage identification in beams using speckle shearography and an optimal spatial sampling. Mech. Syst. Signal Process. 2016, 79, 47–64. [Google Scholar] [CrossRef]
  77. Katunin, A.; Lopes, H.; Araújo dos Santos, J.V. Identification of multiple damage using modal rotation obtained with shearography and undecimated wavelet transform. Mech. Syst. Signal Process. 2019, 116, 725–740. [Google Scholar] [CrossRef]
  78. Araújo dos Santos, J.V.; Katunin, A.; Lopes, H. Vibration-based damage identification using wavelet transform and a numerical model of shearography. Int. J. Struct. Stab. Dyn. 2019, 19, 1950038. [Google Scholar] [CrossRef]
  79. Qian, S. Introduction to Time-Frequency and Wavelet Transforms; China Machine Press: Beijing, China, 2005. [Google Scholar]
  80. Montanari, L.; Basu, B.; Spagnoli, A.; Broderick, M.B. A padding method to reduce edge effects for enhanced damage identification using wavelet analysis. Mech. Syst. Signal Process. 2015, 52–53, 264–277. [Google Scholar] [CrossRef]
  81. Bońkowski, P.A.; Bobra, P.; Zembaty, Z.; Jędraszak, B. Application of rotation rate sensors in modal and vibration analyses of reinforced concrete beams. Sensors 2020, 20, 4711. [Google Scholar] [CrossRef]
  82. Mazanoglu, K. A novel methodology using simplified approaches for identification of cracks in beams. Lat. Am. J. Solids Struct. 2015, 12, 2460–2479. [Google Scholar] [CrossRef] [Green Version]
  83. Liu, W.; Barkey, M.E. Nonlinear vibrational response of a single edge cracked beam. J. Mech. Sci. Technol. 2017, 31, 5231–5243. [Google Scholar] [CrossRef]
  84. Cung, L.T.; Dao, T.D.; Nguyen, P.C.; Bui, T.D. A model-based approach for estimation of the crack depth on a massive metal structure. Meas. Control 2018, 51, 182–191. [Google Scholar] [CrossRef]
  85. Chettah, M.; Lassoued, R. Comparative study of direct and inverse problems of cracked beams. MATEC Web Conf. 2018, 149, 02015. [Google Scholar]
Figure 1. The two-node Euler-Bernoulli finite element marked displacements, w 1 and w 2 , and rotations, θ 1 and θ 2 .
Figure 1. The two-node Euler-Bernoulli finite element marked displacements, w 1 and w 2 , and rotations, θ 1 and θ 2 .
Sensors 20 06638 g001
Figure 2. The geometry of the simulated beam with indication of the damage location.
Figure 2. The geometry of the simulated beam with indication of the damage location.
Sensors 20 06638 g002
Figure 3. Displacements and rotations for the investigated modes.
Figure 3. Displacements and rotations for the investigated modes.
Sensors 20 06638 g003
Figure 4. The influence of the frequency range on the damage detectability: the ranges of 0–10 Hz (a), 0–20 Hz (b), 0–40 Hz (c), 0–70 Hz (d), 0–100 Hz (e), and 0–200 Hz (f).
Figure 4. The influence of the frequency range on the damage detectability: the ranges of 0–10 Hz (a), 0–20 Hz (b), 0–40 Hz (c), 0–70 Hz (d), 0–100 Hz (e), and 0–200 Hz (f).
Sensors 20 06638 g004
Figure 5. The results of the damage identification for the considered damage scenarios using a baseline-free approach.
Figure 5. The results of the damage identification for the considered damage scenarios using a baseline-free approach.
Sensors 20 06638 g005aSensors 20 06638 g005b
Figure 6. The results of the damage identification for the considered damage scenarios using baseline approach based on subtraction of the mode shapes before the transform.
Figure 6. The results of the damage identification for the considered damage scenarios using baseline approach based on subtraction of the mode shapes before the transform.
Sensors 20 06638 g006
Figure 7. The results of the damage identification for the considered damage scenarios using baseline approach based on subtraction of the transformed mode shapes.
Figure 7. The results of the damage identification for the considered damage scenarios using baseline approach based on subtraction of the transformed mode shapes.
Sensors 20 06638 g007aSensors 20 06638 g007b
Figure 8. The results of the damage identification using modal displacements as the input data for the considered approaches.
Figure 8. The results of the damage identification using modal displacements as the input data for the considered approaches.
Sensors 20 06638 g008aSensors 20 06638 g008b
Figure 9. The exemplary results of the quantification of the damage width.
Figure 9. The exemplary results of the quantification of the damage width.
Sensors 20 06638 g009
Figure 10. The exemplary result of the cubic polynomial regression of the quantification results for the mode 1 of the baseline-free approach.
Figure 10. The exemplary result of the cubic polynomial regression of the quantification results for the mode 1 of the baseline-free approach.
Sensors 20 06638 g010
Table 1. The damage scenarios for the damage identification.
Table 1. The damage scenarios for the damage identification.
Scenario1I2I3I4I5I6I7I8I
Slot 1, mm0.10.2230.30.4090.4090.4090.4090.409
Slot 2, mm 0.0280.10.20.3
Table 2. The damage scenarios for the damage quantification.
Table 2. The damage scenarios for the damage quantification.
Scenario1Q2Q3Q4Q5Q6Q7Q
Depth, mm0.030.060.090.150.30.450.6
Table 3. The results of quantification of the damage depth.
Table 3. The results of quantification of the damage depth.
Scenario1Q2Q3Q4Q5Q6Q7Q
Baseline-free approachMode 11.423181.444101.466141.513911.659191.853572.11880
Mode 20.316380.359090.405640.511670.872231.427812.28861
Mode 30.018320.019450.020720.023730.034920.053890.08567
Baseline 1 approachMode 18.328·10−40.003460.008090.024330.119530.334660.75081
Mode 29.535·10−40.003960.009270.027880.137110.384360.86364
Mode 34.108·10−51.709·10−44.01·10−40.001210.006000.017020.03882
Baseline 2 approachMode 10.019860.040770.062810.110590.255860.450240.71548
Mode 20.039170.081870.128420.234450.595021.150602.01139
Mode 39.995·10−40.002130.003390.006410.017590.036560.06835
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Katunin, A. Damage Identification and Quantification in Beams Using Wigner-Ville Distribution. Sensors 2020, 20, 6638. https://doi.org/10.3390/s20226638

AMA Style

Katunin A. Damage Identification and Quantification in Beams Using Wigner-Ville Distribution. Sensors. 2020; 20(22):6638. https://doi.org/10.3390/s20226638

Chicago/Turabian Style

Katunin, Andrzej. 2020. "Damage Identification and Quantification in Beams Using Wigner-Ville Distribution" Sensors 20, no. 22: 6638. https://doi.org/10.3390/s20226638

APA Style

Katunin, A. (2020). Damage Identification and Quantification in Beams Using Wigner-Ville Distribution. Sensors, 20(22), 6638. https://doi.org/10.3390/s20226638

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