Next Article in Journal
Information-Theoretic Analysis of a Family of Improper Discrete Constellations
Next Article in Special Issue
Low Computational Cost for Sample Entropy
Previous Article in Journal
Image Segmentation Based on Statistical Confidence Intervals
Previous Article in Special Issue
Automated Multiclass Classification of Spontaneous EEG Activity in Alzheimer’s Disease and Mild Cognitive Impairment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy-Based Structural Health Monitoring System for Damage Detection in Multi-Bay Three-Dimensional Structures

Department of Civil Engineering, National Chiao Tung University, Hsinchu 30010, Taiwan
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(1), 49; https://doi.org/10.3390/e20010049
Submission received: 28 November 2017 / Revised: 3 January 2018 / Accepted: 5 January 2018 / Published: 11 January 2018

Abstract

:
In this paper, a structural health monitoring (SHM) system based on multi-scale cross-sample entropy (MSCE) is proposed for detecting damage locations in multi-bay three-dimensional structures. The location of damage is evaluated for each bay through MSCE analysis by examining the degree of dissimilarity between the response signals of vertically-adjacent floors. Subsequently, the results are quantified using the damage index (DI). The performance of the proposed SHM system was determined in this study by performing a finite element analysis of a multi-bay seven-story structure. The derived results revealed that the SHM system successfully detected the damaged floors and their respective directions for several cases. The proposed system provides a preliminary assessment of which bay has been more severely affected. Thus, the effectiveness and high potential of the SHM system for locating damage in large and complex structures rapidly and at low cost are demonstrated.

1. Introduction

Structural health monitoring (SHM) has attracted considerable attention among engineers, because structures are inevitably subject to internal or external factors that affect their service lives. SHM enables early detection of any damage caused by factors such as the environment or faulty construction, thus facilitating timely maintenance or repair jobs. Over the past two decades, signal-processing techniques have mainly been applied in SHM methods to analyze the measured displacement, velocity or acceleration signals of structures to obtain dynamic characteristics such as basic vibration frequency (natural frequency) and damping. These characteristics contribute to the diagnosis of damage, as well as its condition and possible location [1,2].
In 1999, Wahab and De Roeck [3] used the change in dynamic parameters between the undamaged and damaged conditions of simply-supported and continuous beams to localize damage. In 2000, Maeck et al. [4] identified the degree and location of damage in reinforced concrete beams according to their modal characteristics for dynamic stiffness analysis. Subsequently, vibration-based SHM algorithms and their application limitations were examined and summarized by Chang et al. [5]. The feasibility of using measured modal parameters for damage detection of steel towers through Bayesian probability theory was examined by Lam and Yang [6]. Furthermore, Chen et al. [7] identified bridge bearing damage by using beam modal information as the input to a neural network.
German physicist Clausius [8] first introduced the concept of entropy in 1865 to evaluate the uncertainty of events in a thermodynamic system. Decades later, Shannon [9] proposed the influential Shannon entropy in the field of information theory. Moreover, Kolmogorov [10] defined the notion of entropy for a new class of dynamical systems; Sinai [11] then introduced a definition of entropy that can be applied to all dynamical systems. Jointly, the entropy defined in these two studies was called KS entropy, which was used to measure the complexity of measured time series in D-dimensional dynamic systems. Thereafter, research demonstrated that although KS entropy can be effectively applied to low-dimensional chaotic systems, it cannot be applied to experimental data. The calculation results are affected by various levels of noise that may be involved [12].
In 1991, Pincus [13] developed an analytical method called “approximate entropy” (ApEn) by modifying the KS entropy formula. This method permitted the quantification of regularity in a time series as a single number. The performance of ApEn was tested using clinical data, and the test results demonstrated the predictive and diagnostic capabilities of ApEn. An and Ou [14] proposed the mean curvature difference method based on the ApEn theory and successfully located the damage in the shear frame structures. A modification of ApEn, called “sample entropy” (SampEn), was proposed by Richman and Moorman [15] in 2000. The advantage of SampEn is that the entropy value obtained is not affected by the length of the time series. Moreover, higher relative consistency can be achieved under different parameters such as the threshold r, sample length m and signal length N. Richman and Moorman also argued that comparing a dataset with itself is meaningless; therefore, SampEn does not count self-matches and is free of any bias caused in the entropy estimation.
In 2002, multi-scale entropy (MSE) was proposed by Costa et al. [16,17] to measure the entropy in physiological time series. Traditional entropy-based algorithms that use only a single time scale for analysis may yield misleading results; therefore, a coarse-graining procedure was proposed to obtain more accurate results during entropy calculation. The MSE method was verified by applying it to heartbeat time series of healthy subjects, subjects with congestive heart failure and subjects with atrial fibrillation. The results showed that on multiple scales, healthy heartbeats had the highest entropy values; thus, diagnosing patients through the MSE method is feasible. MSE analysis has been used not only in the medical field, but also in the area of mechanical engineering. The damage condition of roller bearings has been successfully diagnosed through the MSE analysis of vibration signals [18,19].
Similar to ApEn, Cross-ApEn was introduced by Pincus and Singer [20] to analyze the degree of asynchrony of two related time series. Subsequently, Richman and Moorman [15] proposed Cross-SampEn as a superior alternative for the same purpose, because unmatched templates result in undefined probabilities in Cross-ApEn. The performance of both methods was tested using cardiovascular time series, demonstrating that the results obtained using Cross-SampEn had a higher level of relative consistency than those obtained using Cross-ApEn.
In 2013, Fabris et al. [21] applied the SampEn and Cross-SampEn algorithms to electroglottogram and microphone signals. The healthy patients and those with throat or vocal disorders can be identified by quantifying the degree of asynchrony between time series. Subsequently, SHM systems based on the MSE and MSCE algorithms introduced in [15,16,17] in the field of physiology have been proposed to identify damage locations and directions on a single-bay structure [22,23]. The vertical MSCE analysis was performed to identify the damaged floor, and planar MSCE analysis was performed to identify the damage directions. The resulting MSCE curves indicated that a higher degree of synchronicity between two signals yields lower entropy values. Furthermore, time series with high complexity have high entropy values, indicating damage. Considering the higher entropy values obtained for damaged floors, vertical and planar damage index (DI) values were proposed for efficiently quantifying damage. A comparison of healthy and damaged signals revealed positive DI values for damaged floors, whereas negative DI values were observed for healthy floors. Moreover, the parameters of sample entropy and wavelet transformation were optimized to detect the possible crack on a cantilever beam [24], and a cross-entropy optimization technique was also applied to identify the damage of the shear structure component [25].
The aim of the present study was to implement the vertical MSCE analysis coupled with the vertical DI on a large and complex three-bay bi-axial numerical model to identify the damaged floors and damaged bays of the structure. The remainder of this paper is organized as follows: The proposed SHM system is described in Section 2. In Section 3, a numerical evaluation of a three-bay, seven-story steel structure is presented. On the basis of the numerical results, the performance of the MSCE and DI analyses is discussed in Section 4. Finally, Section 5 provides the conclusions.

2. SHM Algorithm

2.1. SampEn

In this section, SampEn is first introduced to understand more clearly the methods used for the SHM system [13]. As a statistical method for analyzing time series, SampEn estimates the entropy value of a measured time series to quantify the complexity of a system. The results of SampEn, an unbiased refinement of ApEn, are not affected by the time series length or calculation parameters.
For a time series { X i } = { x 1 , , x i , x N } with length N, a vector of m data points u m ( i ) = { x i , x i + 1 , , x i + m 1 } , 1 i N m + 1 can be defined as the template. The combination of all templates with length m is represented by the template space T of the signal; for example, [ x i , x i + 1 , , x i + m 1 ] represents the i-th template of the time series. Various N m + 1 templates may constitute the time series. The template space T, which is the combination of all N m + 1 templates with length m, is expressed as follows:
T = x 1 x 2 x m x 2 x 3 x m 1 x N m + 1 x N m + 2 x N
Let d i j be the maximum distance between two templates i and j and r be a predetermined threshold.
d i j = max { | x ( i + k ) x ( j + k ) | : 0 k m 1 }
Next, the number of similarities n i m ( r ) between templates u m ( i ) and u m ( j ) can be calculated as follows:
n i m ( r ) = j = 1 N m d [ u m ( i ) , u m ( j ) ]
where similarity is defined as follows:
d [ u m ( i ) , u m ( j ) ] = 1 d i j r 0 d i j > r
The distance between templates is calculated through Equation (2) and then substituted into Equation (4) to define the similarity between the two. The two templates are determined to be similar when the distance d i j does not exceed the threshold r. By contrast, the two templates are dissimilar when d i j exceeds r. Different pattern templates can be substituted for comparisons with template i. The degree of sample similarity U i m ( r ) can then be calculated as follows:
U i m ( r ) = n i m ( r ) ( N m + 1 )
The average degree of sample similarity can be calculated as follows after obtaining the degree of sample similarity:
U m ( r ) = 1 N m i = 1 N m U i m ( r )
Here, U m ( r ) represents the average degree of similarity between all templates of length m in the template space T. Subsequently, a new template space is created by assembling templates with length m + 1. The average degree of similarity U m + 1 ( r ) of the new template space is calculated by repeating the aforementioned steps. Consequently, the SampEn values of the time series with parameters m, r and N can be obtained as follows:
S E ( m , r , N ) = ln U m + 1 ( r ) U m ( r )

2.2. MSE

MSE can extract much more information of a time series, compared with single scale-based entropy methods. Briefly, a time series is subjected to a coarse-graining procedure to construct multiple time series at different time scales [15,16]. The procedure is described as follows: A discrete time series x 1 , x 2 , , x N with length N is segmented into multiple time series with length τ , where τ ( τ = 1 , 2 , , N ) is the scale factor. Subsequently, a new time series { y j ( τ ) } is constructed by deriving the arithmetic mean of each set of data values according to the following equation:
y j ( τ ) = 1 τ i = ( j 1 ) τ + 1 j τ x i , 1 j N τ
The length of each coarse-grained time series is N / τ , meaning that at scale 1, the coarse-grained time series is the original time series. Moreover, the length of the coarse-grained time-series decreases as τ increases. After the process is completed, SampEn is used to calculate the entropy values for each coarse-grained time series { y j ( τ ) } . The obtained sample entropy values are the MSE of the time series and are plotted as a function of the scale factor ( f ( τ ) = S E ) .

2.3. Cross-SampEn

With a similar procedure to SampEn, Cross-SampEn is used to evaluate the degree of dissimilarity between two time series derived from the same system [20]. The estimation of Cross-SampEn can be summarized as follows. Consider two individual time series { X i } = { x 1 , , x i , , x N } and { Y j } = { y 1 , , y j , , y N } with length N. The signals are segmented into the following templates of length m: u m ( i ) = { x i , x i + 1 , , x i + m 1 } , 1 i N m + 1 and v m ( j ) = { y j , y j + 1 , , y j + m 1 } , 1 j N m + 1 . The template space T x is presented as follows:
T x = x 1 x 2 x m x 2 x 3 x m 1 x N m + 1 x N m + 2 x N
Similarly, template space T y is expressed as follows:
T y = y 1 y 2 y m y 2 y 3 y m 1 y N m + 1 y N m + 2 y N
The number of similarities between u m ( i ) and v m ( j ) , defined as n i m ( r ) , is calculated under the following criterion:
d [ u m ( i ) , v m ( j ) ] r , 1 j N m
The similarity probability of the templates is evaluated using the following equations:
U i m ( r ) ( v u ) = n m ( r ) ( N m )
Then, the average probability of similarity of length m is calculated as follows:
U m ( r ) ( v u ) = 1 ( N m ) i = 1 N m U i m ( r ) ( v u )
where U m ( r ) ( v u ) is the degree of dissimilarity between the two time series when m points are segmented. Subsequently, new template spaces T x and T y are created by assembling templates with length m + 1 . The steps above are repeated to obtain the average similarity probability U m + 1 ( r ) ( v u ) , and the Cross-SampEn values can be then derived as follows:
C S E ( m , r , N ) = ln U m + 1 ( r ) ( v u ) U m ( r ) ( v u )

2.4. MSCE and DI

The SHM system proposed in this study relies on MSCE analysis for identifying structural damage. Because an extended three-bay model was analyzed in this study, it was also of interest to identify the damaged bay. However, diagnosing the location of damage through a simple observation of the obtained MSCE curves is typically difficult. Therefore, a DI based on [22,23] was proposed for rapidly and efficiently diagnosing the damaged floor, axis and bay in the structure.
For the three-bay structure described in the following section, the SHM process is applied to each bay separately to detect the possible damage, while the interaction between each bay has been included in the measured signals. Thus, the following procedure is repeated three times. Two groups of curves representing the healthy and damaged conditions of the structure are analyzed. For a biaxial structure with N floors, the MSCE curves for the x- and y-axes under the healthy condition can be expressed as matrices:
M S C E u n d a m a g e d = H 1 x H 2 x H N x M S C E u n d a m a g e d H 1 y H 2 y H N y
Similarly, the MSCE curves for the x- and y-axes under the damaged condition are expressed as follows:
M S C E d a m a g e d = D 1 x D 2 x D N x M S C E d a m a g e d D 1 y D 2 y D N y
H and D represent the MSCE curves for the healthy and damaged conditions, respectively. The superscripts x and y represent the analyzed axes, and the subscript depicts the analyzed floor. For example, H N x is the x-axis MSCE of the signal in the analyzed N-th floor and the signal of the floor beneath it under healthy conditions. This can be further expressed as follows:
H N a x i s = { C S E H N a x i s 1 , C S E H N a x i s 2 , C S E H N a x i s 3 , , C S E H N a x i s τ }
Similarly, D N a x i s can be arranged as follows:
D N a x i s = { C S E D N a x i s 1 , C S E D N a x i s 2 , C S E D N a x i s 3 , , C S E D N a x i s τ }
where C S E denotes the Cross-SampEn value in each element, the superscript denotes the scale factor τ , the first superscript denotes the health condition, the second subscript denotes the analyzed floor and the third subscript denotes the analyzed axis. Subsequently, the following formulae can be used to calculate the dual-axis DI per floor on a specific bay:
D I N x = k = 1 τ ( C S E D N x q C S E H N x q ) D I N y = k = 1 τ ( C S E D N y q C S E H N y q )
The DI is evaluated by calculating the difference between the MSCE values of the damaged and healthy structures. Each bay of the structure has two DI values per floor: one on the x-axis and the other on the y-axis. For a specific floor, a positive DI indicates that the floor has sustained damage, whereas a negative DI indicates no damage on the floor, where the structure is under a more stable condition compared to the original healthy state. As verified in previous research [22], the proposed method can sustain the velocity changes for an approximately 10–20% noise level.

3. Numerical Simulation

To verify the feasibility of the SHM system proposed in this study, SAP2000 software (SAP2000 v9) was used to construct and analyze the three-bay, seven-story model, which is an extension of the areal benchmark structure commonly used at the National Center for Research on Earthquake Engineering (NCREE), for numerical simulation. The dimensions and characteristics of the numerical model are outlined as follows: The model was a steel structure comprising seven stories and three bays on the x-axis and a single bay on the y-axis. The height of each story was 1.06 m, and the widths of the bays on the x- and y-axes were 1.32 and 0.92 m, respectively. The columns were steel plates measuring 75 × 50 mm 2 . The beams were steel plates of 70 × 100 mm 2 . All sides of the structure were fitted with steel bracing chosen as L-shaped steel angles measuring 65 × 65 × 6 mm. Apart from the self-weight of the structure, an additional 500-kg mass was added per bay on each story to simulate the actual characteristics of a structure.
Time history and modal analyses were performed on the model. To perform the time history analysis, white noise signals were first generated with a power intensity of 1 MW and selected as the input accelerations to excite the numerical model as a simulation of ambient vibrations. The direct integration time-history analysis was applied to solve the structural response for 150 s. The response signals of the undamaged scenario under ambient vibrations shown in Figure 1 were used as the reference of the SHM database for damage detection. Detailed setting of the finite element model is listed in Table 1.

Damage Database

The numerical model of the steel structure was fitted with diagonal braces on all directions of every floor, implying that for the three-bay model, each floor had six braces on the strong (X) axis and four braces on the weak (Y) axis. Structural damage was simulated by removing the installed braces symmetrically; that is, two braces were removed on the x-axis to simulate damage on the strong axis, and two bracings were removed on the y-axis to simulate damage on the weak axis.
After the SHM database was constructed by performing the time history and modal analyses on the numerical model for each damage condition, biaxial velocity response data were extracted from the center of each floor per bay. Figure 2 shows the data extraction points of the structure. The various damage conditions along with their respective cases are presented in Table 2. The cases comprised various combinations of single-story, two-story or multistory damage, paired with single- or multi-bay and single- or multi-direction damage. The scenarios for the damage database were selected with a broad spectrum of diverse damage locations. In total, the damage conditions were classified into 12 categories and 26 cases. Numbers in the case names indicate damaged floors; X and Y represent the damaged axes; and L, C and R (left, center and right) denote the damaged bay. For example, Case 9 (6X-L & 6Y-R) represents a case involving damage on the sixth floor, x-axis, left bay and sixth floor, y-axis, right bay. Figure 3 illustrates the numerical model, with the dotted lines representing the braces removed to exemplify the damaged bracing for Case 9. A list of the modal analysis for all damage cases is shown in Table 3 to reflect the change of the global behavior of the numerical model. As depicted, the fundamental frequencies of the first and second modes drop accompanying the increase of damage level in both the X and Y directions. Significant changes are caused by damages located in the lower floors.

4. Results and Discussion

For the undamaged case and every damage condition, the velocity signals for the long (X) and short (Y) axes were extracted per bay from the center of each floor. Subsequently, the signals of two vertically-adjacent floors under the same damage condition were processed through Cross-SampEn at multiple scales (MSCE) to evaluate the dissimilarity between floors. After a series of optimization searches [23], where different combinations of parameters were considered for the best performance on damage detection accuracy, the parameters required for the calculation of Cross-SampEn, such as the template length m, threshold r and signal length N, were selected as 4, 0.10 × the standard deviation (SD) of the time series and 30,000, respectively. In addition, Cross-SampEn was calculated across 25 scales ( τ = 25).
After analyzing the undamaged case and the 26 damage cases through the MSCE method, the results were compared and quantified per bay using the DI. When a specific floor sustains damage, the structure experiences a loss of stiffness, which causes a rise in the complexity of the extracted response signals for that specific floor. Time series with high complexity have high entropy values, indicating damage. Thus, analysis of the changes in complexity or entropy of healthy and damaged signals results in positive DI values for damaged floors and negative DI values for healthy floors. Small positive values can be excluded by a predetermined threshold value of one, which was selected by experimentation. Thus, the damaged floor and direction can be easily detected.
Figure 4a–f presents the MSCE diagrams obtained for the healthy condition in the X- and Y-directions for the left, center and right bays; the diagrams of the X-direction are shown on the left side, and those of the Y-direction are on the right side. In these figures, G-1F denotes the curve for the first floor, 1F-2F denotes the curve for the second floor, 2F-3F denotes the curve for the third floor, and so forth. The Cross-SampEn for Channel G-1F was calculated from the velocity signals of the ground and the first floor.

4.1. Single-Story, Single-Bay, Multidirectional Damage: Case 4 (4XY-L)

Figure 5 illustrates the MSCE diagrams for Case 4. In the X-direction, the rise of the curve for the fourth floor is evident in the MSCE diagrams of all bays, whereas the remaining curves are in similar positions as in the healthy condition. In the Y-direction, the rise of the fourth-floor curve and the dropping of the remaining curves are also noticeable in the diagrams of all bays. Furthermore, Figure 5b shows a large gap between the fourth-floor curve and the remaining curves at scale factors ranging from 8–25. These results indicate damage on the fourth floor in both directions, with a higher degree of damage possibly occurring on the left bay.
The DI analysis results shown in Figure 6 indicate that the fourth floor sustained damage in the X- and Y-directions. Figure 6c,d enables a closer evaluation of the DI values of the fourth floor. In both directions, the bar representing the left bay is higher than the remaining bars, indicating that this bay could have sustained more damage. In the Y-direction, the DI values are shown to be higher, and the differences among bays are shown to be more significant than those in the X-direction. Because the analyzed model involved three bays in the X-direction and a single bay in the Y-direction, the removal of bracings in the short (Y) axis could foster a larger impact on its stiffness; therefore, more significant differences could be observed.

4.2. Single-Story, Multi-Bay, Multidirectional Damage: Case 8 (3X-R & 3Y-C)

The results for the MSCE analysis of Case 8 are shown in Figure 7. In the X-direction, the MSCE diagrams for the three bays are very similar to one another. The figure reveals that the curves for the third floor are in higher positions than they are in the diagram for the healthy condition. Furthermore, at scales ranging from 4–24, the third-floor curves are shown to maintain significantly higher positions than the remaining curves. In the Y-direction, the curves for the third floor are also shown to rise relative to their original positions in the diagram for the undamaged condition, with the difference being slightly higher for the MSCE diagram for the center bay (Figure 7d). Case 8 was therefore classified as involving damage in both directions of the third floor.
The DI diagrams for Case 8 are illustrated in Figure 8. In the X-direction, a positive DI value is shown for the third floor, with the right bay bar revealing the highest DI. Although the second floor also had a positive DI, it was excluded because the value fell under the threshold. Because damage happened on the third floor, an effect on the adjacent second floor is natural due to the production of a rigid body response. The DI diagram for the Y-direction also reveals a positive DI for the third floor. In this case, the center bay bar is shown to be the highest. These results, which can be examined more closely in Figure 8c,d, confirm the MSCE analysis findings, with a further indication of the possible damaged bay.

4.3. Two-Story, Single-Bay, Single-Direction Damage: Case 11 (1Y-R & 5Y-R)

The results obtained from the MSCE analysis of Case 11 are presented in Figure 9. In the X-direction, the MSCE curves for all bays remain almost identical to those of the undamaged condition; therefore, no damage occurred in the X-direction. Figure 9b illustrates the MSCE curves for the left bay in the Y-direction, revealing that that the curve for the first floor is at a higher position than the remaining curves at most scales. Furthermore, at scales of 5–15, a rise in complexity can be observed for the fifth floor. The changes in the curves for the first and fifth floors are even more evident in the MSCE diagrams for the center and right bays (Figure 9d,f). Therefore, the first and fifth floors were damaged in the Y-direction. Additionally, the curve for the fifth floor in the right bay displays higher entropy values than in the center bay for scales ranging from 10–25, suggesting a higher degree of damage to the right bay.
The DI analysis results are shown in Figure 10, indicating no damage in the X-direction. Moreover, in the Y-direction, the first and fifth floors are demonstrated to have positive DI values. Figure 10c shows a close-up of the DI values obtained for the first and fifth floors in the Y-direction; the bars representing the right bay exhibit the highest values. Therefore, the first and fifth floors were damaged in the Y-direction, with the right bay sustaining a higher degree of damage. Hence, the MSCE analysis results are confirmed for this case.

4.4. Two-Story, Multi-Bay, Single-Direction Damage: Case 16 (2X-L & C, 6X-C & R)

Figure 11 illustrates the MSCE curves obtained from the analysis of Case 16. The curves for the second and sixth floors are shown to clearly separate and maintain a wide gap from the remaining curves in the MSCE diagrams obtained from analysis of the three bays in the X-direction. The curves for the second and sixth-floor reach their peaks at scales 5 and 7, respectively. In the Y-direction, no evident changes are shown in the diagrams for any bay; therefore, damage was determined to have occurred on only the second and sixth floor in the X-direction.
The DI diagrams for Case 16 are shown in Figure 12. In the X-direction, positive DI values are shown for the second and sixth floors, thus resulting in these floors being classified as damaged. The DI value of the second floor is higher than that in the sixth floor because the damage on the second floor can relatively affect the global behavior of the structure compared with similar damage conditions. In this particular case, the second floor was damaged on the left and center bays, whereas the sixth floor was damaged on the center and right bays. The damaged bay could not be identified by the DI during the analysis of damage to a single floor in adjacent bays; however, the results could still indicate that the second and sixth floors received more damage in the left and right bays, respectively. Figure 12c enables a closer observation of these results. Finally, the DI diagram for the Y-direction indicates no damage on this axis.

4.5. Multistory, Multi-Bay, Single-Direction Damage: Case 23 (3X-L & 4X-C & 5X-R)

The MSCE diagrams for Case 23 are illustrated in Figure 13. For the three bays in the X-direction, an increase in complexity is shown for the third and fourth floors. The fifth-floor curve is shown to remain in an almost identical position as in the diagram for the healthy condition, whereas the remaining curves are illustrated to decrease relative to their positions under the healthy condition. Although the fifth floor was incorrectly classified as undamaged, the third and fourth floors were determined as damaged in the X-direction. To determine the damaged bays, further examination was required through DI analysis. No evident changes were observed in the Y-direction.
The DI diagrams for this case are presented in Figure 14, revealing positive DI values for the third and fourth floors in the X-direction. The resulting diagram, a close-up of which is provided in Figure 14c, suggests that the third floor sustained more damage in the left bay, whereas the bar is slightly higher for the center bay on the fourth floor. The fifth floor was not identified as damaged. In the Y-direction, small positive DI values are shown in the figure; however, these values were determined to be negligible because they fell under the predetermined threshold of one, thereby confirming the MSCE analysis results.

4.6. Multistory, Multi-Bay, Multidirectional Damage: Case 26 (7XY-R & 4Y-L & 6Y-C)

Figure 15 illustrates the MSCE curves for Case 26. In the X-direction, the curves for the seventh floor for all bays are shown to increase and reach their peaks at scale 5. However, the curves for the remaining floors are shown to drop or maintain their positions as in the diagrams for the healthy condition. Therefore, the seventh floor was determined to have sustained damage in the X-direction. In the Y-direction, the MSCE diagram for the left bay (Figure 15b), reveals that the curve for the fourth-floor increases dramatically. The trend for the first floor also increases, albeit slightly erratically. A slight increase in the curves for the sixth and seventh floors can be observed in the MSCE diagram for the center bay (Figure 15d). Moreover, a change is observed in the fourth-floor curve, but to a lower degree than that in the left bay MSCE diagram for the left bay. Regarding the MSCE diagram for the right bay (Figure 15f), the curves for the first, sixth and seventh floors increase significantly. On the basis of these results, the first, fourth, sixth and seventh floors sustained damage; however, this was false for the first floor.
Figure 16 illustrates the DI diagrams obtained for this case, with a close-up presented in Figure 16c,d. Positive DI values can be observed for the seventh floor in the X-direction, with the bar representing the right bay having the highest DI value. The first, fourth, sixth and seventh floors were determined to have sustained damage in the Y-direction. However, slightly erratic curves are shown in the figure for the first floor, which may have caused its misclassification. The DI diagrams suggest that the fourth and seventh floors received more damage in the left and right bays, respectively. For the sixth floor, the bar representing the right bay indicates the highest value, possibly because the damage to the center bay also affected the right bay because the bracing was shared by both. The results obtained in the DI analysis for this case demonstrate the accuracy of the proposed system.

4.7. General Discussion

In this study, 26 damage cases classified into 12 categories were examined to verify the feasibility of MSCE and DI analyses for damage detection in a complex three-bay structure. The complete results are summarized in Table 4. The DI values are lower when there are more damaged floors as the structural complexity is redistributed after different damage conditions. These were further analyzed in two stages. First, the DI results obtained for the identification of the damaged floors and damage directions were quantified through a precision and recall analysis. For the precision and recall analyses, the DI results for the X- and Y-directions were first classified into four categories: true positives (TP), representing damaged floors that have been correctly identified; false positives (FP), representing floors that have been misclassified as damaged; true negatives (TN), representing undamaged floors that have been correctly classified; and false negatives (FN), representing damaged floors that have been misclassified as undamaged. Precision and recall can then be calculated as follows:
Precision = T P T P + F P Recall = T P T P + F N

4.7.1. Damage Location

High precision denotes few false positives, which represents the percentage of detecting the damage location correctly, and a high recall indicates few false negatives, which indicates the reliability of not missing the possible damage. The combined results for both directions are summarized in Table 5. The derived precision and recall were 83% and 87%, respectively. These results indicate that 83% of the floors and their respective directions classified by the DI as damaged were true positive. As most of the damages were simulated along the X-direction, better performance was achieved in the X-direction for 100% than the Y-direction, which is 73% in the study. Moreover, 87% of all actual damaged floors and their respective directions were correctly classified as damaged. The high-percentage of accuracy proves that the system can reliably detect any level of damage with only a small probability of misdiagnoses, which are false negatives. The results have demonstrated the capacity of the proposed SHM system to locate the damaged floor and direction of a large, three-bay numerical model.

4.7.2. Damaged Bay Identification

A separate analysis of the identification accuracy of the damaged bays was then performed; therefore, the results pertaining to the damaged bays were not considered for the precision and recall analyses. Because the removal of bracings on a specific bay would inevitably affect those adjacent to it, the calculated DI values would be fairly close. The combined results for both directions are shown in Table 6, indicating an average identification accuracy of 75%. Better performance was achieved in the Y-direction, as it only has one bay in the Y-direction. For the complicated damage combinations in the X-direction, the accuracy dropped to 71%. Furthermore, the proposed system is limited when more than one bay is damaged. Therefore, these results are merely a suggestion as to which bay might have been affected more severely. With the support of the damaged bay identification, the proposed method can be easily extended to any complex structure.

5. Conclusions

In this study, the feasibility of detecting damage in a complex three-bay, seven-story numerical model by SHM methods was examined. Through the proposed SHM system, damage locations can be rapidly and effectively detected by only measuring the velocity response data of the model. The ambient vibrations from the center of each floor before and after damage occurs were first recorded, and subsequently, the complexity of the signals can be analyzed by the MSCE method. The reliability and viability of the proposed SHM system were examined through the numerical analysis of 26 damage cases in 12 categories representing several degrees of damage severity. The results of the analyses were examined in two stages. First, the results pertaining to the damaged floor and direction indicate that 83% of the floors and their respective directions were truly damaged. Furthermore, 87% of all the actually damaged floors and their respective directions were correctly classified as damaged by the DI. Subsequently, the identification of the damaged bays was analyzed, and an identification accuracy of 75% was obtained. Identification of the damaged bays through the proposed SHM system is limited, especially when multi-bay damage exists on a single floor.
The obtained results verify the feasibility and further potential of the proposed SHM system for the detection and localization of damage in large and complex structures.

Acknowledgments

This work is supported by the Ministry of Science and Technology (MOST) (105-2221-E-009-056-) at Taipei, Taiwan.

Author Contributions

Tzu-Kang Lin conceived of and put forward the research ideas. Ana Gabriela Laínez carried out the numerical analysis and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Friswell, M.; Penny, J.; Garvey, S. Parameter subset selection in damage location. Inverse Probl. Eng. 1997, 5, 189–215. [Google Scholar] [CrossRef]
  2. Doebling, S.W.; Farrar, C.R.; Prime, M.B. A summary review of vibration-based damage identification methods. Shock Vib. Dig. 1998, 30, 91–105. [Google Scholar] [CrossRef]
  3. Wahab, M.A.; 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]
  4. Maeck, J.; Wahab, M.A.; Peeters, B.; De Roeck, G.; De Visscher, J.; De Wilde, W.; Ndambi, J.M.; Vantomme, J. Damage identification in reinforced concrete structures by dynamic stiffness determination. Eng. Struct. 2000, 22, 1339–1349. [Google Scholar] [CrossRef]
  5. Chang, P.C.; Flatau, A.; Liu, S. Health monitoring of civil infrastructure. Struct. Health Monit. 2003, 2, 257–267. [Google Scholar] [CrossRef]
  6. Lam, H.F.; Yang, J. Bayesian structural damage detection of steel towers using measured modal parameters. Earthq. Struct. 2015, 8, 935–956. [Google Scholar] [CrossRef]
  7. Chen, Z.; Fang, H.; Ke, X.; Zeng, Y. A new method to identify bridge bearing damage based on Radial Basis Function Neural Network. Earthq. Struct. 2016, 11, 841–859. [Google Scholar] [CrossRef]
  8. Clausius, R. On several convenient forms of the fundamental equations of the mechanical theory of heat. Ann. Phys. 1865, 201, 352–400. [Google Scholar]
  9. Shannon, C.E. A mathematical theory of communication, Part I, Part II. Bell Syst. Tech. J. 1948, 27, 623–656. [Google Scholar] [CrossRef]
  10. Kolmogorov, A.N. New metric invariant of transitive dynamical systems and endomorphisms of Lebesgue spaces. Dokl. Russ. Acad. Sci. 1958, 119, 861–864. [Google Scholar]
  11. Sinai, Y.G. On the notion of entropy of a dynamical system. Dokl. Akad. Nauk. SSSR 1959, 124, 768–771. [Google Scholar]
  12. Pincus, S.M.; Gladstone, I.M.; Ehrenkranz, R.A. A regularity statistic for medical data analysis. J. Clin. Monit. Comput. 1991, 7, 335–345. [Google Scholar] [CrossRef]
  13. Pincus, S.M. Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef] [PubMed]
  14. An, Y.-H.; Ou, J.-P. Structural damage localisation for a frame structure from changes in curvature of approximate entropy feature vectors. Nondestruct. Test. Eval. 2014, 29, 80–97. [Google Scholar] [CrossRef]
  15. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol.-Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed]
  16. Costa, M.; Goldberger, A.L.; Peng, C.K. Multiscale entropy analysis of complex physiologic time series. Phys. Rev. Lett. 2002, 89, 068102. [Google Scholar] [CrossRef] [PubMed]
  17. Costa, M.; Goldberger, A.L.; Peng, C.K. Multiscale entropy analysis of biological signals. Phys. Rev. E 2005, 71, 021906. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, L.; Xiong, G.; Liu, H.; Zou, H.; Guo, W. Bearing fault diagnosis using multi-scale entropy and adaptive neuro-fuzzy inference. Expert Syst. Appl. 2010, 37, 6077–6085. [Google Scholar] [CrossRef]
  19. Liu, H.; Han, M. A fault diagnosis method based on local mean decomposition and multi-scale entropy for roller bearings. Mech. Mach. Theory 2014, 75, 67–78. [Google Scholar] [CrossRef]
  20. Pincus, S.; Singer, B.H. Randomness and degrees of irregularity. Proc. Natl. Acad. Sci. USA 1996, 93, 2083–2088. [Google Scholar] [CrossRef] [PubMed]
  21. Fabris, C.; De Colle, W.; Sparacino, G. Voice disorders assessed by (cross-) sample entropy of electroglottogram and microphone signals. Biomed. Signal Process. Control 2013, 8, 920–926. [Google Scholar] [CrossRef]
  22. Lin, T.K.; Liang, J.C. Application of multi-scale (cross-) sample entropy for structural health monitoring. Smart Mater. Struct. 2015, 24, 085003. [Google Scholar] [CrossRef]
  23. Lin, T.K.; Tseng, T.C.; Lainez, A.G. Three-dimensional structural health monitoring based on multiscale cross-sample entropy. Earthq. Struct. 2017, 12, 673–687. [Google Scholar]
  24. Wimarshana, B.; Wu, N.; Wu, C. Crack identification with parametric optimization of entropy & wavelet transformation. Struct. Monit. Maint. 2017, 4, 33–52. [Google Scholar]
  25. Guan, X.; Wang, Y.; He, J. A probabilistic damage identification method for shear structure components based on Cross-Entropy optimizations. Entropy 2017, 19, 27. [Google Scholar] [CrossRef]
Figure 1. The time history of the velocity response of the undamaged scenario.
Figure 1. The time history of the velocity response of the undamaged scenario.
Entropy 20 00049 g001
Figure 2. Signal extraction points.
Figure 2. Signal extraction points.
Entropy 20 00049 g002
Figure 3. Three-dimensional view of the numerical model. The dotted braces represent the damaged bracing for Case 9 (6X-L & 6Y-R).
Figure 3. Three-dimensional view of the numerical model. The dotted braces represent the damaged bracing for Case 9 (6X-L & 6Y-R).
Entropy 20 00049 g003
Figure 4. MSCE diagrams for undamaged condition: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Figure 4. MSCE diagrams for undamaged condition: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Entropy 20 00049 g004aEntropy 20 00049 g004b
Figure 5. MSCE diagrams for Case 4: damage on fourth floor, left bay, X- and Y-directions: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Figure 5. MSCE diagrams for Case 4: damage on fourth floor, left bay, X- and Y-directions: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Entropy 20 00049 g005
Figure 6. Damage index (DI) diagrams for Case 4: (a) X-direction; (b) Y-direction; (c) close up of DI diagrams for the fourth floor in Case 4: X direction; (d) Y-direction.
Figure 6. Damage index (DI) diagrams for Case 4: (a) X-direction; (b) Y-direction; (c) close up of DI diagrams for the fourth floor in Case 4: X direction; (d) Y-direction.
Entropy 20 00049 g006
Figure 7. MSCE diagrams for Case 8: damage on the third floor, right bay, X-direction and third floor, center bay, Y-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Figure 7. MSCE diagrams for Case 8: damage on the third floor, right bay, X-direction and third floor, center bay, Y-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Entropy 20 00049 g007
Figure 8. DI diagrams for Case 8: (a) X-direction; (b) Y-direction; close-up of DI diagrams for the third floor in Case 8: (c) X-direction; (d) Y-direction.
Figure 8. DI diagrams for Case 8: (a) X-direction; (b) Y-direction; close-up of DI diagrams for the third floor in Case 8: (c) X-direction; (d) Y-direction.
Entropy 20 00049 g008
Figure 9. MSCE diagrams for Case 11: damage on first floor, right bay, Y-direction and fifth floor, right bay, Y-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Figure 9. MSCE diagrams for Case 11: damage on first floor, right bay, Y-direction and fifth floor, right bay, Y-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Entropy 20 00049 g009
Figure 10. DI diagrams for Case 11: (a) X-direction; (b) Y-direction; (c) close-up of DI diagrams for the first and fifth floors in Case 11 (Y-direction).
Figure 10. DI diagrams for Case 11: (a) X-direction; (b) Y-direction; (c) close-up of DI diagrams for the first and fifth floors in Case 11 (Y-direction).
Entropy 20 00049 g010
Figure 11. MSCE diagrams for Case 16: damage on second floor, left and center bays, X-direction and sixth floor, center and right bays, X-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Figure 11. MSCE diagrams for Case 16: damage on second floor, left and center bays, X-direction and sixth floor, center and right bays, X-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Entropy 20 00049 g011
Figure 12. DI diagrams for Case 16: (a) X-direction; (b) Y-direction; (c) close-up of DI diagrams for the second and sixth floors in Case 16 (X-direction).
Figure 12. DI diagrams for Case 16: (a) X-direction; (b) Y-direction; (c) close-up of DI diagrams for the second and sixth floors in Case 16 (X-direction).
Entropy 20 00049 g012
Figure 13. MSCE diagrams for Case 23: damage on the third floor, left bay, X-direction; fourth floor, center bay, X-direction; fifth floor, right bay, X-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Figure 13. MSCE diagrams for Case 23: damage on the third floor, left bay, X-direction; fourth floor, center bay, X-direction; fifth floor, right bay, X-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Entropy 20 00049 g013
Figure 14. DI diagrams for Case 23: (a) X-direction; (b) Y-direction; (c) close-up of DI diagrams for the third and fourth floors in Case 23 (X-direction).
Figure 14. DI diagrams for Case 23: (a) X-direction; (b) Y-direction; (c) close-up of DI diagrams for the third and fourth floors in Case 23 (X-direction).
Entropy 20 00049 g014
Figure 15. MSCE diagrams for Case 26: damage on seventh floor, right bay, X- and Y-directions; fourth floor, left bay, Y-direction; sixth floor, center bay, Y-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Figure 15. MSCE diagrams for Case 26: damage on seventh floor, right bay, X- and Y-directions; fourth floor, left bay, Y-direction; sixth floor, center bay, Y-direction: (a) left bay, X-direction; (b) left bay, Y-direction; (c) center bay, X-direction; (d) center bay, Y-direction; (e) right bay, X-direction; (f) right bay, Y-direction.
Entropy 20 00049 g015
Figure 16. DI diagrams for Case 26: (a) X-direction; (b) Y-direction; (c) seventh floor, X-direction; (d) fourth, sixth and seventh floor, Y-direction.
Figure 16. DI diagrams for Case 26: (a) X-direction; (b) Y-direction; (c) seventh floor, X-direction; (d) fourth, sixth and seventh floor, Y-direction.
Entropy 20 00049 g016
Table 1. Numerical setting of the finite element model.
Table 1. Numerical setting of the finite element model.
ParameterSetting
Constitutive equationLinear elastic
MaterialA36 steel
DOF of the foundationFixed in 6 DOF
Finite element for column, beam and bracingBeam
Integration methodDirect-integration time-history analysis
Duration150 s
Ambient vibrationWhite noise of 1 mW power
Table 2. Damage cases for analysis.
Table 2. Damage cases for analysis.
Case NumberDamage GroupDamaged Floor, Direction and Bay
1Single-story, single-bay, single-direction5X-L
23Y-C
37Y-R
4Single-story, single-bay, multidirectional4XY-L
56XY-C
6Single-story, multi-bay, single-direction2X-L & 2X-C
75Y-L & 5Y-C & 5Y-R
8Single-story, multi-bay, multidirectional3X-R & 3Y-C
96X-L & 6Y-R
10Two-story, single-bay, single-direction3X-L & 6X-L
111Y-R & 5Y-R
12Two-story, single-bay, multidirectional4X-C & 7Y-C
132XY-R & 3XY-R
14Two-story, multi-bay, single-direction5X-R & 7X-L
152Y-C & 4Y-R
162X-L & C, 6X-C & R
17Two-story, multi-bay, multidirectional4X-R & 2Y-L
186XY-R & 7XY-L
19Multistory, single-bay, single-direction3X-L & 4X-L & 6X-L
201Y-R & 4Y-R & 7Y-R
21Multistory, single-bay, multidirectional4X-L & 5Y-L & 6Y-L
221XY-C & 3XY-C & 5XY-C
23Multistory, multi-bay, single-direction3X-L & 4X-C & 5X-R
246Y-L & 2Y-C & 7Y-R
25Multistory, multi-bay, multidirectional1X-R & 2X-R & 1Y-L
267XY-R & 4Y-L & 6Y-C
Table 3. Modal analysis of the numerical model.
Table 3. Modal analysis of the numerical model.
Case No.Damage CaseFrequency (Hz)
Mode 1(X)Mode 2(X)Mode 1(Y)Mode 2(Y)
Undamaged (UN)6.519.522.6911.56
15X-L6.4118.842.6911.56
23Y-C6.5019.532.6511.55
37Y-R6.5017.972.6811.02
44XY-L6.1218.692.6011.32
56XY-C6.4418.582.6811.15
62X-L & 2X-C5.7818.932.6911.57
75Y-L & 5Y-C & 5Y-R6.2717.412.448.88
83X-R & 3Y-C6.3219.522.6511.55
96X-L & 6Y-R6.4317.022.6610.48
103X-L & 6X-L6.2918.812.6911.57
111Y-R & 5Y-R5.8817.552.499.21
124X-C & 7Y-C6.3219.092.6911.36
132XY-R & 3XY-R5.4118.232.3910.86
145X-R & 7X-L6.4218.552.6911.56
152Y-C & 4Y-R6.319.162.5511.02
162X-L & C, 6X-C & R5.6716.212.6911.57
174X-R & 2Y-L5.9118.662.5711.00
186XY-R & 7XY-L6.4315.502.6610.19
193X-L & 4X-L & 6X-L6.1318.522.6911.57
201Y-R & 4Y-R & 7Y-R5.7816.132.459.33
214X-L & 5Y-L & 6Y-L6.2917.562.599.96
221XY-C & 3XY-C & 5XY-C5.9218.072.5710.54
233X-L & 4X-C & 5X-R6.1218.462.7011.57
246Y-L & 2Y-C & 7Y-R6.4916.442.609.94
251X-R & 2X-R & 1Y-L5.6718.092.5610.40
267XY-R & 4Y-L & 6Y-C6.2715.942.5910.52
Table 4. Classification results of DI analysis.
Table 4. Classification results of DI analysis.
Case NumberDamage GroupDamage CaseDamage Index (X-Direction)Damage Index (Y-Direction)
1Single-story,5X-LOKOK
2single-bay,3Y-COKOK
3single-direction7Y-ROK1F 1
4Single-story, single-4XY-LOKOK
5bay, multi-directional6XY-C6F 2 OK
6Single-story, multi-2X-L & 2X-C2F 2 OK
7bay, single-direction5Y-L & 5Y-C & 5Y-ROK5F 2
8Single-story, multi-3X-R & 3Y-COKOK
9bay, multidirectional6X-L & 6Y-ROK1F 1
10Two-story, single-bay,3X-L & 6X-LOKOK
11single-direction1Y-R & 5Y-ROKOK
12Two-story, single-bay,4X-C & 7Y-COKOK
13multidirectional2XY-R & 3XY-ROK 1 F 1
14Two-story,5X-R & 7X-LOKOK
15multi-bay,2Y-C & 4Y-ROK1F 1 & 2F 2
16single-direction2X-L & C, 6X-C & R2F 2 & 6F 2 OK
17Two-story, multi-4X-R & 2Y-LOK1F 1
18bay, multidirectional6XY-R & 7XY-LOK1F 1 & 2F 1
19Multistory, single-3X-L & 4X-L & 6X-L6F 3 OK
20bay, single-direction1Y-R & 4Y-R & 7Y-ROKOK
21Multistory, single-4X-L & 5Y-L & 6Y-LOK1F 1 & 2F 1
22bay, multidirectional1XY-C & 3XY-C & 5XY-C1F 3 & 5F 3 1F 3 & 3F 2
23Multistory, multi-bay,3X-L & 4X-C & 5X-R5F 3 OK
24single-direction6Y-L & 2Y-C & 7Y-ROK1F 1 & 2F 3
25Multistory, multi-bay,1X-R & 2X-R & 1Y-L1F 3 & 2F 3 OK
26multidirectional7XY-R & 4Y-L & 6Y-COK1F 1 & 6F 2
1 Indicates that the floor has been misclassified as damaged; 2 indicates that the damaged bay has not been successfully identified; 3 indicates that the damaged floor has not been detected.
Table 5. Precision and recall analysis results (X- and Y-directions).
Table 5. Precision and recall analysis results (X- and Y-directions).
DirectionTrue PositivesFalse PositivesTrue NegativesFalse NegativesPrecisionRecall
X2501516100%81%
Y3011139273%94%
Total5511290883%87%
Table 6. Identification accuracy of damaged bays (X- and Y-directions).
Table 6. Identification accuracy of damaged bays (X- and Y-directions).
DirectionDamage InstancesCorrectly Identified BayAccuracy
X342471%
Y342779%
Total685175%

Share and Cite

MDPI and ACS Style

Lin, T.-K.; Laínez, A.G. Entropy-Based Structural Health Monitoring System for Damage Detection in Multi-Bay Three-Dimensional Structures. Entropy 2018, 20, 49. https://doi.org/10.3390/e20010049

AMA Style

Lin T-K, Laínez AG. Entropy-Based Structural Health Monitoring System for Damage Detection in Multi-Bay Three-Dimensional Structures. Entropy. 2018; 20(1):49. https://doi.org/10.3390/e20010049

Chicago/Turabian Style

Lin, Tzu-Kang, and Ana Gabriela Laínez. 2018. "Entropy-Based Structural Health Monitoring System for Damage Detection in Multi-Bay Three-Dimensional Structures" Entropy 20, no. 1: 49. https://doi.org/10.3390/e20010049

APA Style

Lin, T. -K., & Laínez, A. G. (2018). Entropy-Based Structural Health Monitoring System for Damage Detection in Multi-Bay Three-Dimensional Structures. Entropy, 20(1), 49. https://doi.org/10.3390/e20010049

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