# Damage Location and Quantification Indices of Shear Structures Based on Changes in the First Two or Three Natural Frequencies

Authors: Hien HoThu; Akira Mita
Journal: International Scholarly Research Notices
 (2014)

---

## Abstract

This study proposes damage detection algorithms for multistory shear structures that only need the first two or three natural frequencies. The methods are able to determine the location and severity of damage on the basis of damage location indices (DLI) and damage quantification indices (DQI) consisting of the changes in the first few squared natural frequencies of the undamaged and damaged states. The damage is assumed to be represented by a reduction in stiffness. This stiffness reduction causes a shift in the natural frequencies of the structure. The uncertainty associated with system identification methods for obtaining natural frequencies is also carefully considered. The methods are accurate and cost-effective means only requiring the changes in the natural frequencies.


---

## Body

## 1. Introduction

Structural health monitoring (SHM) systems are garnering attention as a way of maintaining building structures subject to natural hazards such as large earthquakes and strong winds [1]. SHM systems play an important role in assessing the health of a structure because they can determine the location and the severity of damage. The obstacle we often face when installing an SHM system in a building is the trade-off between the numbers of sensors and the accuracy of damage detection. A large number of sensors are costly and entail a large effort for wiring and installation. Such complicated and expensive SHM systems are not feasible for most buildings [2]. Therefore, a good SHM system should have few sensors yet obtaining enough information about the health of a structure.Doebling et al. [3] graded SHM systems into four levels of ability, as follows:(i)
Level 1: determining that damage is present in the structure;
(ii)
Level 2: determination of the location of the damage;
(iii)
Level 3: quantification of the severity of the damage;
(iv)
Level 4: prediction of the remaining service life of the structure.The first three levels are most often related to structural dynamic testing and modeling issues. Level 4 is not addressed in the structural vibration or modal analysis literature. Hence, most of damage detection methods aim to classify damage into the first three levels.Damage detection algorithms based on the modal properties of a structure, such as modal frequencies, mode shapes, curvature mode shapes, and modal flexibilities, have been studied in the SHM field for decades. However, most algorithms have difficulties in identifying the precise location and magnitude of the damage. Their accuracy and reliability are not considered sufficient, if not completely inadequate [2]. The key to making a successful damage detection method is thus using a few modal properties of a structure to identify Levels 1, 2, and 3.Zhao and DeWolf [4] presented a sensitivity study comparing the use of natural frequencies, mode shapes, and modal flexibilities for monitoring. Based on the fact that natural frequencies are sensitive indicators of structural integrity, the relationship between frequency changes and structural damage was discussed in a review by Salawu [5]. These studies showed that sensitivity analysis of the natural frequencies can be a valuable tool in SHM.A damage detection method based on natural frequencies only needs two vibration sensors (or even one acceleration sensor) to obtain the modal frequencies. It is known that, of the various characteristics, the natural frequencies are the least contaminated by measurement noise and can generally be measured with good accuracy [6]. Messina et al. [7] suggested a standard error of 0.15% as a benchmark figure for natural frequencies measured in the laboratory with the impulse hammer technique. Some researches [8–14] have achieved a standard deviation of less than 1% for the first few modes of natural frequencies. This low level of error suggests that a damage detection method using only information on the natural frequencies would have acceptable accuracy.Although many previous studies concluded that the frequency cannot provide spatial information about structural changes, multiple frequency shifts may provide spatial information about structural damage in situations where many natural frequencies can be measured [15]. However, there are only a few natural frequencies that can be measured in most buildings. Therefore, a natural frequency based method would need to work with only a few frequencies if it were to be practical.The purpose of this study is to devise a damage detection method to identify the existence, location, and amount of damage to multistory shear structures by using new damage indices consisting of the changes in the first two or three natural frequencies.


## 2. Sensitivity of Squared Natural Frequency Changes to Structural Damage

A multistory shear structure (N-story) can be modeled as a one-dimensional lumped mass shear model, as shown in Figure 1. Most of the damage to a structure, such as cracks, fatigue, corrosion, and loosening of bolted joints, manifests itself as a stiffness reduction.Figure 1
Simplified structural model withN degrees of freedom.The characteristic equation for such a structure is written as(1)
[
K
-
ω
r
2
M
]
∅
r
=
0
,
where ω
r and ∅
r are the rth frequency and mode shape vector, respectively,K is the stiffness matrix (N
×
N),(2)
K
=
[
k
1
+
k
2
-
k
2
0
-
k
2
k
2
+
k
3
-
k
3
0
⋱
⋱
…
0
⋱
⋮
⋱
-
k
i
k
i
+
k
i
+
1
-
k
i
+
1
⋮
⋱
0
⋯
⋱
⋱
0
-
k
N
-
1
k
N
-
1
+
k
N
-
k
N
0
-
k
N
k
N
]
,andM is the mass matrix (N
×
N),
(3)
M
=
[
m
1
0
0
m
2
⋱
…
0
⋱
⋮
m
i
⋮
⋱
0
⋯
⋱
m
N
-
1
0
0
m
N
]
.Solving (1) in the manner of (
8
) in [15], the general expression for the ith story stiffness can be obtained as
(4)
k
i
=
ω
r
2
(
ϕ
i
r
-
ϕ
(
i
-
1
)
r
)
∑
j
=
i
N
m
j
ϕ
j
r
,
where ϕ
i
r and ϕ
(
i
-
1
)
r are the rth mode shapes at the ith and (
i
-
1
)th stories, respectively.In addition, the squared natural circular frequency of therth mode is expressed in [2] as follows:
(5)
ω
r
2
=
K
r
M
r
=
ϕ
r
T
K
ϕ
r
ϕ
r
T
M
ϕ
r
,
where ω
r is the rth natural circular frequency, K
r is the rth modal stiffness, M
r is the rth modal mass, K and M are the stiffness and mass matrices, and ϕ
r is the rth mode shape vector.Although the yielded structure is undamaged and may also have reduced stiffness sometime later, in this study, we will assume that the structural damage can be directly expressed as a stiffness reduction and that the structural mass remains unchanged. The difference between the squared natural frequencies of therth mode of the undamaged structure and the same structure with damage to the ith story is expressed as
(6)
ω
r
2
-
ω
r
(
i
)
2
=
ϕ
r
T
K
ϕ
r
ϕ
r
T
M
ϕ
r
-
ϕ
r
(
i
)
T
K
i
ϕ
r
(
i
)
ϕ
r
(
i
)
T
M
ϕ
r
(
i
)
,
where the subscript (i) denotes the ith damaged story state.Here,ω
r and ω
r
(
i
) are the rth frequencies of the undamaged and ith story damaged state, K
i is the global stiffness matrix with the stiffness reduction considered to be a reflection of damage to the ith story, and is the rth mode shape vector of the ith damaged story state.In [16], Morita et al. gave an equation to obtain the difference in the squared natural frequency of the rth mode in (6) by neglecting the difference in mode shape arising from damage. That difference becomes
(7)
ω
r
2
-
ω
r
(
i
)
2
=
ϕ
r
T
Δ
K
i
ϕ
r
ϕ
r
T
M
ϕ
r
,
where Δ
K
i is the change in the stiffness matrix with the ith damaged story. When the ith story stiffness is reduced, only k
(
i
-
1
)
(
i
-
1
), k
(
i
-
1
)
i, k
i
(
i
-
1
) and k
i
i are altered in the stiffness matrix. Δ
K
i is written as
(8)
Δ
K
i
=
[
0
…
0
⋮
k
(
i
-
1
)
(
i
-
1
)
k
(
i
-
1
)
i
k
i
(
i
-
1
)
k
i
i
⋮
0
…
0
]
or
(9)
Δ
K
i
=
[
0
…
0
⋮
Δ
k
i
-
Δ
k
i
-
Δ
k
i
Δ
k
i
⋮
0
…
0
]
,
where Δ
k
i is the change in stiffness of the ith story.The squared natural frequency change ratio is determined by dividing both sides of (7) by ω
r
2:
(10)
Δ
i
(
r
)
=
ω
r
2
-
ω
r
(
i
)
2
ω
r
2
=
ϕ
r
T
Δ
K
i
ϕ
r
ω
r
2
ϕ
r
T
M
ϕ
r
=
ϕ
r
T
Δ
K
i
ϕ
r
ϕ
r
T
K
ϕ
r
(11)
⟺
Δ
i
(
r
)
=
1
ϕ
r
T
K
ϕ
r
Δ
k
i
(
ϕ
i
r
-
ϕ
(
i
-
1
)
r
)
2
,
where Δ
k
i is the change in stiffness of the ith story and ϕ
i
r and ϕ
(
i
-
1
)
r are the rth mode shapes at the ith and (
i
-
1
)th story.From (11), Δ
k
i is the extent of damage to the ith story and (
ϕ
i
r
-
ϕ
(
i
-
1
)
r
)
2 depends on the location of the ith damaged story. As damage normally causes a decrease in the natural frequencies, Δ
i
(
r
) is positive when damage occurs. This frequency change ratio Δ
i
(
r
) may vary depending on the damage location and quantification. Note that, for each natural frequency, there are some locations where the frequency is most sensitive to the damage, while there are other locations where the damage has little influence on the frequency.Similarly, the frequency change ratio of thesth mode, Δ
i
(
s
), is obtained as
(12)
Δ
i
(
s
)
=
1
ϕ
s
T
K
ϕ
s
Δ
k
i
(
ϕ
i
s
-
ϕ
(
i
-
1
)
s
)
2
.The sum of these changes is(13)
Δ
i
(
r
)
+
Δ
i
(
s
)
=
Δ
k
i
[
1
ϕ
r
T
K
ϕ
r
(
ϕ
i
r
-
ϕ
(
i
-
1
)
r
)
2
0000000
+
1
ϕ
s
T
K
ϕ
s
(
ϕ
i
s
-
ϕ
(
i
-
1
)
s
)
2
]
.Thus, the ratio of the changes in the two natural frequencies of therth and sth modes can be written as
(14)
Δ
i
(
r
)
Δ
i
(
r
)
+
Δ
i
(
s
)
=
(
1
/
ϕ
r
T
K
ϕ
r
)
(
ϕ
i
r
-
ϕ
(
i
-
1
)
r
)
2
(
1
/
ϕ
r
T
K
ϕ
r
)
(
ϕ
i
r
-
ϕ
(
i
-
1
)
r
)
2

+

(
1
/
ϕ
s
T
K
ϕ
s
)
(
ϕ
i
s
-
ϕ
(
i
-
1
)
s
)
2
.
(15)
⟹
Δ
i
(
r
)
Δ
i
(
r
)
+
Δ
i
(
s
)
=
1
1
+
(
ϕ
r
T
K
ϕ
r
/
ϕ
s
T
K
ϕ
s
)
(
(
ϕ
i
s
-
ϕ
(
i
-
1
)
s
)
/
(
ϕ
i
r
-
ϕ
(
i
-
1
)
r
)
)
2From (4), we can use the rth mode or sth mode to obtain the stiffness of the ith story:
(16)
ω
r
2
(
∅
i
r
-
∅
(
i
-
1
)
r
)
∑
j
=
i
N
m
j
∅
j
r
=
ω
s
2
(
∅
i
s
-
∅
(
i
-
1
)
s
)
∑
j
=
i
N
m
j
∅
j
s
(17)
⟹
ϕ
i
s
-
ϕ
(
i
-
1
)
s
ϕ
i
r
-
ϕ
(
i
-
1
)
r
=
ω
s
2
ω
r
2
∑
j
=
i
N
m
j
∅
j
s
∑
j
=
i
N
m
j
∅
j
r
.Substituting (17) into (15), we get
(18)
Δ
i
(
r
)
Δ
i
(
r
)
+
Δ
i
(
s
)
=
1
1

+

(
ϕ
r
T
K
ϕ
r
/
ϕ
s
T
K
ϕ
s
)
(
ω
s
4
/
ω
r
4
)
(
∑
j
=
i
N
m
j
∅
j
s
/
∑
j
=
i
N
m
j
∅
j
r
)
2
.This ratio depends on the location of damagei and the number of degrees of freedom, in the following coefficient: (
∑
j
=
i
N
m
j
∅
j
s
/
∑
j
=
i
N
m
j
∅
j
r
)
=
(
m
i
∅
i
s
+
m
i
+
1
∅
(
i
+
1
)
s
+
⋯
+
m
N
∅
N
s
)
/
(
m
i
∅
i
r
+
m
i
+
1
∅
(
i
+
1
)
r
+
⋯
+
m
N
∅
N
r
) for i
=
1
,
2
,
…
,
N.


## 3. Proposed Approach to Damage Detection

Let us formulate a simple algorithm for detecting structural damage, its location, and its amount. Note that the existence of damage is determined by a change in the first two or three natural frequencies. In previous, some researches [8–14] have achieved a standard deviation of less than 1% for the first few modes of natural frequencies. In particular, we assume that the squared natural frequencies have a standard deviation of 1% and damage has occurred when the squared natural frequency changes by 2% or more.


### 3.1. Damage Location Indices

The previous section derived the ratio of the changes in two squared natural frequencies of therth and sth modes. Like (18), this value depends on the location of the damage. The damage location index (DLI) is defined as
(19)
DL
I
r
s
=
Δ
i
(
r
)
Δ
i
(
r
)
+
Δ
i
(
s
)
.The number of obtained frequencies is usually smaller than the number of degrees of freedomN. When two of the first two or three frequencies are used, the DLI values are written as DLI
12, DL
I
23, and so on:
(20)
DLI
12
=
Δ
i
(
1
)
Δ
i
(
1
)
+
Δ
i
(
2
)
DL
I
23
=
Δ
i
(
2
)
Δ
i
(
2
)
+
Δ
i
(
3
)
.TheDLI
12 value depends on the number of damaged stories i, and it can be used as a damage location indicator in the N
-story shear structure. When taller building has many stories, DLI
23 and/or DL
I
34 and so forth, can be used in addition to DLI
12 to detect the damaged story.


### 3.2. Uncertainty of Damage Location Indices

The modal parameters are often sensitive to various environmental conditions such as temperature, humidity, and excitation amplitude. The effect of environmental conditions or excitation amplitude is treated as “noise” in a simulation, so we should obtain a confidence interval on the modal parameters.Denoting byσ
ω
2 the standard deviation of the squared natural frequency, we will discuss the reliability of our method on the basis of the theory presented in [17]. We will assume that the squared natural frequencies have a standard deviation of 1%.The standard deviation of the difference between the squared natural frequencies is(21)
σ
(
ω
r
2
-
ω
r
(
i
)
2
)
=
(
σ
ω
r
2
)
2
+
(
σ
ω
r
(
i
)
2
)
2
.Moreover, the standard deviation of the changes in therth frequency Δ
(
r
) is given by
(22)
σ
Δ
(
r
)
=
Δ
(
r
)
(
σ
(
ω
r
2
-
ω
r
(
i
)
2
)
ω
r
2
-
ω
r
(
i
)
2
)
2
+
(
σ
ω
r
2
ω
r
2
)
2
.The standard deviation of the changes in thesth frequency σ
Δ
(
s
) is calculated similarly to σ
Δ
(
r
) in (22).The standard deviation of the sum of the changes in the first two frequencies is(23)
σ
(
Δ
(
r
)
+
Δ
(
s
)
)
=
(
σ
Δ
(
r
)
)
2
+
(
σ
Δ
(
s
)
)
2
.The standard deviation of theDLI can be calculated as
(24)
σ
DLI
=
DLI
(
σ
Δ
(
r
)
Δ
(
r
)
)
2
+
(
σ
(
Δ
(
r
)
+
Δ
(
s
)
)
Δ
(
r
)
+
Δ
(
s
)
)
2
.Statistically, we expect theDLI to be reliable within 〈
DLI
±
2
σ
DLI
〉. Hence, the confidence interval is 95%, and the DLI can identify the damage location in a shear structure with satisfactory reliability.


### 3.3. Damage Quantification Indices

Now let us formulate a simple damage quantification technique. The change in the natural frequency can be used to interpolate the amount of damage at the detected location [18]. Equation (11) also shows that the change in the natural frequency is proportional to the stiffness reduction. In addition, Zhao and DeWolf [4] said that the change in natural frequency has a different sensitivity level for each location. That is, some stories show less of change in the natural frequency when the damage occurs. As [18], we can use the first frequency change Δ
(
1
) to interpolate the extent of damage. However, these values are too small for some free-end stories that are less sensitive to damage. Similarly, the second frequency change Δ
(
2
) also has less sensitive locations. In these cases, it is difficult to interpolate the extent of damage from the change in only one frequency.TheDLI defined in Section 3.1 was used to detect the location of a damaged story in the shear structure. Normally, DLI
12 is composed of the changes in the first two or three natural frequencies. Equation (7) in [16] can be used to approximate the extent of damage from the change in any one of the obtained frequencies. A more reliable value can be found by averaging the estimates of all of the modes used in the analysis.Because the sensitivity to a stiffness reduction depends on the modes, it would be a good idea to introduce the following averaging method. The damage quantification index (DQI) is defined by the changes in the first two or three frequencies that were used to determine the damaged story:
(25)
DQI
=
∑
j
=
1
2
(
3
)
w
j
Δ
k
i
(
j
)
∑
j
=
1
2
(
3
)
w
j
=
w
1
Δ
k
i
(
1
)
+
w
2
Δ
k
i
(
2
)
(
+
w
3
Δ
k
i
(
3
)
)
w
1
+
w
2
(
+
w
3
)
,
where Δ
k
i
(
j
) is the stiffness reduction of the damaged story ith, interpolated by the change in square frequency of the jth mode (
Δ
i
(
j
)
), and w
j is a weight with the variance, given by w
j
=
1
/
(
σ
Δ
(
j
)
)
2.The reliability band ofΔ
k
i
(
j
) is obtained from σ
Δ
i
(
j
) in (22) and the standard deviation of DQI, (
σ
DQI
), is also calculated.


## 3.1. Damage Location Indices

The previous section derived the ratio of the changes in two squared natural frequencies of therth and sth modes. Like (18), this value depends on the location of the damage. The damage location index (DLI) is defined as
(19)
DL
I
r
s
=
Δ
i
(
r
)
Δ
i
(
r
)
+
Δ
i
(
s
)
.The number of obtained frequencies is usually smaller than the number of degrees of freedomN. When two of the first two or three frequencies are used, the DLI values are written as DLI
12, DL
I
23, and so on:
(20)
DLI
12
=
Δ
i
(
1
)
Δ
i
(
1
)
+
Δ
i
(
2
)
DL
I
23
=
Δ
i
(
2
)
Δ
i
(
2
)
+
Δ
i
(
3
)
.TheDLI
12 value depends on the number of damaged stories i, and it can be used as a damage location indicator in the N
-story shear structure. When taller building has many stories, DLI
23 and/or DL
I
34 and so forth, can be used in addition to DLI
12 to detect the damaged story.


## 3.2. Uncertainty of Damage Location Indices

The modal parameters are often sensitive to various environmental conditions such as temperature, humidity, and excitation amplitude. The effect of environmental conditions or excitation amplitude is treated as “noise” in a simulation, so we should obtain a confidence interval on the modal parameters.Denoting byσ
ω
2 the standard deviation of the squared natural frequency, we will discuss the reliability of our method on the basis of the theory presented in [17]. We will assume that the squared natural frequencies have a standard deviation of 1%.The standard deviation of the difference between the squared natural frequencies is(21)
σ
(
ω
r
2
-
ω
r
(
i
)
2
)
=
(
σ
ω
r
2
)
2
+
(
σ
ω
r
(
i
)
2
)
2
.Moreover, the standard deviation of the changes in therth frequency Δ
(
r
) is given by
(22)
σ
Δ
(
r
)
=
Δ
(
r
)
(
σ
(
ω
r
2
-
ω
r
(
i
)
2
)
ω
r
2
-
ω
r
(
i
)
2
)
2
+
(
σ
ω
r
2
ω
r
2
)
2
.The standard deviation of the changes in thesth frequency σ
Δ
(
s
) is calculated similarly to σ
Δ
(
r
) in (22).The standard deviation of the sum of the changes in the first two frequencies is(23)
σ
(
Δ
(
r
)
+
Δ
(
s
)
)
=
(
σ
Δ
(
r
)
)
2
+
(
σ
Δ
(
s
)
)
2
.The standard deviation of theDLI can be calculated as
(24)
σ
DLI
=
DLI
(
σ
Δ
(
r
)
Δ
(
r
)
)
2
+
(
σ
(
Δ
(
r
)
+
Δ
(
s
)
)
Δ
(
r
)
+
Δ
(
s
)
)
2
.Statistically, we expect theDLI to be reliable within 〈
DLI
±
2
σ
DLI
〉. Hence, the confidence interval is 95%, and the DLI can identify the damage location in a shear structure with satisfactory reliability.


## 3.3. Damage Quantification Indices

Now let us formulate a simple damage quantification technique. The change in the natural frequency can be used to interpolate the amount of damage at the detected location [18]. Equation (11) also shows that the change in the natural frequency is proportional to the stiffness reduction. In addition, Zhao and DeWolf [4] said that the change in natural frequency has a different sensitivity level for each location. That is, some stories show less of change in the natural frequency when the damage occurs. As [18], we can use the first frequency change Δ
(
1
) to interpolate the extent of damage. However, these values are too small for some free-end stories that are less sensitive to damage. Similarly, the second frequency change Δ
(
2
) also has less sensitive locations. In these cases, it is difficult to interpolate the extent of damage from the change in only one frequency.TheDLI defined in Section 3.1 was used to detect the location of a damaged story in the shear structure. Normally, DLI
12 is composed of the changes in the first two or three natural frequencies. Equation (7) in [16] can be used to approximate the extent of damage from the change in any one of the obtained frequencies. A more reliable value can be found by averaging the estimates of all of the modes used in the analysis.Because the sensitivity to a stiffness reduction depends on the modes, it would be a good idea to introduce the following averaging method. The damage quantification index (DQI) is defined by the changes in the first two or three frequencies that were used to determine the damaged story:
(25)
DQI
=
∑
j
=
1
2
(
3
)
w
j
Δ
k
i
(
j
)
∑
j
=
1
2
(
3
)
w
j
=
w
1
Δ
k
i
(
1
)
+
w
2
Δ
k
i
(
2
)
(
+
w
3
Δ
k
i
(
3
)
)
w
1
+
w
2
(
+
w
3
)
,
where Δ
k
i
(
j
) is the stiffness reduction of the damaged story ith, interpolated by the change in square frequency of the jth mode (
Δ
i
(
j
)
), and w
j is a weight with the variance, given by w
j
=
1
/
(
σ
Δ
(
j
)
)
2.The reliability band ofΔ
k
i
(
j
) is obtained from σ
Δ
i
(
j
) in (22) and the standard deviation of DQI, (
σ
DQI
), is also calculated.


## 4. Performance of Proposed Method

A four-story shear structure was measured to show the feasibility of the proposed method. It was modeled as a one-dimensional lumped mass shear model, as shown in Figure1. The damping ratio for each mode was chosen to be 3%. The sampling frequency 200 Hz is commonly used in major monitoring systems for real buildings. However, as we can only need lower natural frequencies, the sampling frequency can be as low as 50 Hz provided that an appropriate antialiasing filter is used. In this study, the data sampling frequency was chosen 200 Hz. The excitation loading is white-noise. Although we used white-noise in this paper, the excitation can be any, such as winds, microtremors, and earthquakes.The stories of the structure had the same mass,m
i
=
1000 tons, but their stiffness was different: k
1
=
k
=
1.3
×
10
3 MN/m, k
2
=
0.9
k, k
3
=
0.8
k, and k
4
=
0.7
k. The first two natural frequencies of the undamaged structure were 1.89 Hz and 5.19 Hz. Four cases of damage (damage to the 1st, 2nd, 3rd, and 4th stories) were studied. The damage was simulated by reducing the stiffness of each story by 5%, 10%, 15%, 20%, 25%, and 30%.Figure2 shows the sensitivity of the changes in the first two natural frequencies to the location of the damaged story. From Figure 2, it can be seen that the influence of the damage location identified by the first natural frequency decreases on higher stories; the results mean that the first natural frequency is the most sensitive to damage on the 1st story. By contrast, the second natural frequency is the least sensitive to damage on the 2nd story.Sensitivity to changes in square natural frequencies (Δ
(
1
) and Δ
(
2
)) of 4-story structure.
(a)
(b)Figure3 uses different colors to plot DLI
12 values and their confidence interval of 5% to 30% stiffness reduction for each story. The DLI
12 values of these cases are indicated by the bold lines in the middle of each color; their reliability band considering 2
σ
DLI is indicated by the two thin lines above and below each thick line (the standard deviation of the squared natural frequencies is assumed to be 1%).Figure 3
DL
I
12 values and their intervals within the range 〈
DLI
±
2
σ
DLI
〉 of 4-story structure.From Figure3, we can see the following.(i)
TheDLI
12 values are stable at each damaged story for any level of damage; the location of the damaged story can be detected by examining the correlation between the unknown data and the data of the story on which damage was detected.
(ii)
The variances of2
σ
DLI depend on the stiffness reduction; when the reduction in stiffness is smaller, these values become larger. So the extent of damage, as indicated by the stiffness reduction, must be big enough to detect the location of damage without any mistakes.In this simulation, two modes were enough to detect the damaged story. After determining the location of the damage by using theDLI values presented in Figure 3, the quantification of damage on that story can be calculated from the DQI definition with its standard deviation. Figure 4 compares the DQI with its uncertainty against real damage to a 4-story structure that reduced the stiffness of each story by 5, 15, and 25%.Figure 4
Comparison ofDQI with uncertainty and actual damage of the 4-story structure.


## 5. Application to Tall Buildings

Next, we verified our method on an eight-story structure. The damping ratio was chosen to be 3% and the data sampling frequency was 200 Hz; the excitation loading is white-noise. The mass of each story was 1000 tons. The stiffness of the first story was assumed to bek
1
=
k
=
1.3
×
10
3 MN/m, and the stiffness of the other stories was k
2
=
0.95
k, k
3
=
0.9
k, k
4
=
0.85
k, k
5
=
0.8
k, k
6
=
0.75
k, k
7
=
0.7
k, and k
8
=
0.65
k.The first three natural frequencies of the undamaged structure were calculated to be 0.99, 2.82, and 4.58 Hz. Table1 summarizes these three cases in this simulation. The first three natural frequencies were obtained from acceleration data of the sensor on top of the simulated structure by applying the subspace identification method described in [19]. The resulting DLI values are displayed in Table 2.Table 1
Damage scenarios of eight-story structure.
Damage case number
Location of damage (damaged story)
Quantification of damages (percentage of the stiffness reduction)
1
1st
23%
2
4th
14%
3
7th
6%Table 2
Natural frequencies of 3 damage cases of eight-story structure.
Damage case number
The 1st natural frequency,ω
1
[Hz]
The 2nd natural frequency,ω
2
[Hz]
The 3rd natural frequency,ω
3
[Hz]
DLI12
DLI23
1
0.97
2.75
4.47
0.51
0.52
2
0.98
2.82
4.50
0.95
0.04
3
0.99
2.76
4.49
0.15
0.55Eight cases of damage (from the 1st to 8th story) were simulated. For each case, the damage was simulated by reducing the stiffness of each story from 5 to 25%, and these 21 levels of damage were reflected in the changes in the natural frequencies. TheDLI
12 values and their reliability band were calculated by assuming that the squared frequency had a standard deviation of 1% (Figure 5). There is some overlap of the 95% confidence intervals. The DLI
12 values could detect the location of the damage in three sections of the structure. The DLI
12 values of the three damage cases are plotted as the dash dotted lines in Figure 5.Figure 5
DL
I
12 values and their intervals for 8-story structure.The next step is using the 2nd and 3rd frequency changes to getDLI
23 and their reliability band. Figures 6, 7, and 8 plot the DLI
23 values.Figure 6
DL
I
23 and DQI values and their intervals for the blue sections.Figure 7
DL
I
23 and DQI values and their intervals for the green section.Figure 8
DL
I
23 and DQI values and their intervals for the purple section.After determining the location of the damage by using theDLI values, the amount of damage on that story was calculated by interpolating the DQI values. The DLI and DQI values together with their uncertainty in the three damage cases are displayed as the red crosses in Figures 6 to 8.Damage indices are dependent on the mass and stiffness distribution. Thus, we need a prior knowledge of them to apply our method. If design drawing is not available, some system identification tools may be needed. This is the limitation of this method to be applied to a real building.


## 6. Conclusion

The presented method determines the location and amount of damage to shear structures by using the changes in the first few natural frequencies only. Natural frequencies decrease as a result of damage and the damage indices (DLI and DQI) use these changes to reveal the location and amount of damage. The DLI uses two natural frequency changes between undamaged and damaged states and is stable for each damaged story. After identifying the location of the damage, the DQI of the damaged story is used to quantify the damage.As we need only a few natural frequencies, two vibration sensors are enough to obtain the modal frequencies, one on the ground detecting an input and the other on the roof detecting an output. If the input lasts long and the spectrum is flat, we may identify those parameters using the output data without input information. Thus, in such a case, only one sensor is needed.The uncertainty associated with system identification methods for obtaining natural frequencies was also carefully considered, and the confidence intervals of theDLI values were acceptable with high accuracy. The method was also shown to be able to detect the location of damage to tall buildings.



---

*Source: 459097-2014-08-27.xml*