Relevance feature selection of modal frequency-ambient condition pattern recognition in structural health assessment for reinforced concrete buildings

Abstract

Modal frequency is an important indicator for structural health assessment. Previous studies have shown that this indicator is substantially affected by the fluctuation of ambient conditions, such as temperature and humidity. Therefore, recognizing the pattern between modal frequency and ambient conditions is necessary for reliable long-term structural health assessment. In this article, a novel machine-learning algorithm is proposed to automatically select relevance features in modal frequency-ambient condition pattern recognition based on structural dynamic response and ambient condition measurement. In contrast to the traditional feature selection approaches by examining a large number of combinations of extracted features, the proposed algorithm conducts continuous relevance feature selection by introducing a sophisticated hyperparameterization on the weight parameter vector controlling the relevancy of different features in the prediction model. The proposed algorithm is then utilized for structural health assessment for a reinforced concrete building based on 1-year daily measurements. It turns out that the optimal model class including the relevance features for each vibrational mode is capable to capture the pattern between the corresponding modal frequency and the ambient conditions.

Keywords

Bayesian inference feature selection maximum likelihood model class selection structural health monitoring

Introduction

The goal of structural health monitoring (SHM) is to assess the health status of a structure based on structural responses and ambient conditions measurement.¹ Modal frequency, which is related to structural stiffness, is an important indicator in SHM. Previous studies have shown that this indicator is substantially affected by the fluctuation of ambient conditions, such as temperature and relative humidity.^2–6 In order to depict the relationship between modal parameters and ambient conditions, a number of long-term monitoring systems were operated for different types of structures, including bridge,^7–9 reinforced concrete buildings,^10–12 and other types.^13–15 In these studies, it has been shown that the fluctuation of the structural dynamical properties was significant due to the variation of ambient conditions, and the modal frequencies exhibited strong correlation with temperature and relative humidity. If the modal frequency-ambient condition pattern is not appropriately captured, there will be substantial bias in the structural health assessment results. Toward the goal for reliability enhancement of structural health assessment under changing ambient conditions, pattern recognition and machine-learning approaches have received tremendous attention.^16,17 A number of machine-learning algorithms have been proposed, including principal component analysis,¹⁸ support vector machine,¹⁹ nonlinear principal component analysis,²⁰ support vector machine multi-class clustering algorithm,²¹ and kernel-based algorithm.²²

Due to the fact that structural health assessment exhibits significant level of uncertainty,^23–26 a number of system identification approaches have been proposed.^27–32 In particular, the Bayesian inference–based approach has attracted special attention as it provides a rigorous solution for uncertainty quantification and it is applicable for different problems in pattern recognition and system identification,^33–35 such as modal/model updating,^36–38 robust signal processing, and sensor configuration.^39–43 In this article, a novel machine-learning approach is proposed to automatically select the relevance features in modal frequency-ambient condition pattern recognition based on structural dynamic response and ambient condition measurement.

In the traditional feature selection approaches,^44,45 a set of model class candidates is first constructed, and each candidate contains a subset of the components of the full extracted feature vector. Model plausibility evaluation is implemented for each candidate, and it is obvious that the total number of candidates for evaluation is huge although the number of the components of the extracted feature vector is moderate. For instance, when the number of components of the full extracted feature vector is 10, 2¹⁰-1 candidates in total are required for consideration given that at least one feature of the full extracted feature vector should be included in a candidate. Therefore, optimal model class determination in the traditional approaches requires large computational effort. In contrast, the proposed approach is capable of conducting continuous relevance feature selection, which is done by introducing a sophisticated hyperparameterization on the weight parameter vector controlling the relevancy of different features in the prediction model. The proposed algorithm is then utilized for structural health assessment for a reinforced concrete building based on 1-year daily measurements.

The structure of this article is outlined as follows. Dataset and feature extraction are first presented. Then, the approach for relevance feature selection is developed. Finally, the proposed approach is utilized for modal frequency-ambient condition pattern recognition for the monitored structure based on measurement.

Dataset and feature extraction

Description of dataset

The structure considered in this study is a 22-story reinforced concrete residential building, namely, the East Asia Hall. Figure 1 shows the side view and its layout plan. The floor plan of the building is in asymmetric L-shape, with the building height and the length of the two spans being 64.70, 51.90, and 61.75 m, respectively.

Figure 1.

(a) Side view and (b) layout of the East Asia Hall.

The monitoring period is 1 year between 2 May 2008 and 1 May 2009. During the monitoring period, the structure was subjected to four typhoons of 15 days. In order to focus on the pattern between the structural properties and ambient conditions under normal weather conditions, the 15-day typhoon attacking period is excluded and the rest of 350 datasets are utilized in the analysis. Two accelerometers, operated based on standard exploration geophone mass-spring systems with 50 V/g sensitivity, were installed on the 18th floor with two orthogonal directions (directions 1 and 2) shown in Figure 1(b).

In order to minimize the spatial temperature difference of the monitored structure, acceleration was measured at 11:00 p.m. every day. Based on the 10-min acceleration time history of each day, the Bayesian spectral density approach^46–48 is applied to identify the modal frequencies of the building. The Bayesian spectral density approach is a well-developed probabilistic approach for modal identification.^46–48 The modal identification is performed by utilizing the spectral density obtained from the structural response as the data to estimate the uncertain modal parameters. The efficiency and flexibility of the Bayesian spectral density approach on modal identification had been demonstrated through successful applications with in-field measurements.^11,12,38 Herein, each set of the 10-min acceleration measurement is partitioned into four segments with an equal time duration and then the averaged spectrum is obtained for modal identification. Since the squared modal frequencies of first three modes are considered, the range $[0, 4] s^{- 2}$ that covers the concerned squared modal frequencies is used in the modal identification.

Figure 2 shows the identified squared modal frequencies of the first three modes ( $ω_{1}^{2}, ω_{2}^{2}, ω_{3}^{2}$ ) and measurement of the temperature (T) and the relative humidity (H). The ranges of temperature and relative humidity are $13 \circ C \leq T \leq 35.0 \circ C$ and $32.0 % \leq H \leq 95.0 %$ , respectively. It can be seen that high temperature associates with high relative humidity, while low temperature associates with scattering relative humidity, reflecting the weather pattern influenced by monsoon and subtropical climate.

Figure 2.

Identified squared modal frequencies of the first three modes ( $ω_{1}^{2}, ω_{2}^{2}, ω_{3}^{2}$ ) and measurement of the temperature (T) and the relative humidity (H).

Feature extraction

Extracted features are constructed based on the previous research of the relationship between modal frequency and ambient conditions. Watson and Rajapakse⁴⁹ correlated structural properties with temperature based on a polynomial function; Xia et al.⁵ proposed a linear relationship between modal frequencies and temperature along with humidity; Rincón et al.⁵⁰ attempted to investigate the influence of structure due to the changing of relative humidity; Yuen and Kuok^11,12 proposed a second-order polynomial function for modal frequencies with respect to both temperature along with relative humidity; Moser and Moaveni⁵¹ and Moaveni and Behmanesh⁵² utilized high-order (up to fourth) polynomials for correlating modal frequencies and temperature only. In summary, most of the existing works considered low-order polynomials for modal frequencies with respect to temperature only or temperature and relatively humidity; on the other hand, very few works considered high-order polynomials (up to fourth) for frequencies with respect to temperature only, neglecting the cross combinations of temperature and relatively humidity. In this study, in order to comprehensively take different possible patterns into consideration, the extracted features include not only higher order functions of the temperature and the relative humidity but also the cross combinations of them. Finally, the feature vector x containing totally 15 components is given as follows

x = {[1, T_{n}, H_{n}, T_{n}^{2}, T_{n} H_{n}, H_{n}^{2}, T_{n}^{3}, T_{n}^{2} H_{n}, T_{n} H_{n}^{2}, H_{n}^{3}, T_{n}^{4}, T_{n}^{3} H_{n}, T_{n}^{2} H_{n}^{2}, T_{n} H_{n}^{3}, H_{n}^{4}]}^{T}

(1)

where the normalized temperature $T_{n} = T / T_{\max}$ and the normalized relative humidity $H_{n} = H / H_{\max}$ are obtained by considering the rescale of T and H by the maximum value $T_{\max} = \max (T)$ and $H_{\max} = \max (H)$ , respectively. Figure 3 shows the correlation coefficients between the squared modal frequencies and different feature components for the three modes. It is observed that all the modal frequencies possess positive correlation with each component of the feature vector. Toward the goal for reliability enhancement of the prediction model charactering modal frequency-ambient condition pattern, it is aimed to select the suitable relevance features based on the available dataset. In next section, a novel Bayesian inference–based learning approach is proposed. In contrast to the traditional approach requiring selecting relevance features by examining a large number of combinations of extracted features,^44,45 the proposed approach conducts continuous relevance feature selection by introducing a sophisticated hyperparameterization on the weight parameter vector controlling the relevancy of different features in the prediction model.

Figure 3.

Correlation coefficients between the squared modal frequency and each feature component for the three modes: (a) upper: mode 1; (b) middle: mode 2; and (c) lower: mode 3.

Relevance feature selection

Consider the mth vibrational mode of the structure. The modal frequency-ambient condition pattern given the feature vector with $N_{b}$ ( $N_{b} = 15$ ) components for the dataset with N ( $N = 350$ ) points is given as follows

Y_{m} = X b_{m} + ε_{m}

(2)

where $Y_{m} = [ω_{m, 1}^{2}, \dots, ω_{m, N}^{2}]^{T} \in R^{N}$ is the mth mode output vector; $X = [X_{1}, \dots, X_{N_{b}}] \in R^{N \times N_{b}}$ is the design/input matrix characterized by the feature vector of equation (1); $b_{m} \in R^{N_{b}}$ is the mth mode weight parameter vector to be trained; and $ε_{m} \in R^{N}$ is the mth mode residual vector following Gaussian distribution $N (ε_{m} | 0, σ_{m}^{2} I_{N})$ with the corresponding mth mode prediction-error variance $σ_{m}^{2}$ . Then, the likelihood function for $D_{m} = {X, Y_{m}}$ of the mth mode is given as follows⁴⁸

p (D_{m} | b_{m}, σ_{m}^{2}) = N (Y_{m} | X b_{m}, σ_{m}^{2} I_{N})

(3)

In order to conduct continuous relevance feature selection for the mth mode, a sophisticated hyperparameterization on the weight parameter vector $b_{m}$ is introduced. The automatic relevance determination (ARD) prior is adopted for hyperparameterization,^53–55 which is defined as a zero-mean Gaussian distribution

p (b_{m} | α_{m}) = N (b_{m} | 0, A_{m}^{- 1} (α_{m}))

(4)

where $A_{m} (α_{m}) = diag {α_{m, 1}, α_{m, 2}, \dots, α_{m, N_{b}}}$ is the precision matrix (inverse of the covariance matrix) and $α_{m} = [α_{m, 1}, α_{m, 2}, \dots, α_{m, N_{b}}]^{T}$ is the hyperparameter vector parameterizing the prior probability density function (PDF) of the weight parameter vector $b_{m}$ .

Using Bayes’ theorem, the posterior PDF for $b_{m}$ is given by⁵⁶

\begin{matrix} p (b_{m} | D_{m}, α_{m}, σ_{m}^{2}) & \propto p (D_{m} | b_{m}, σ_{m}^{2}) p (b_{m} | α_{m}) \\ = N (b_{m} | {\hat{b}}_{m}, {\hat{Σ}}_{m}) \end{matrix}

(5)

where

{\hat{b}}_{m} = σ_{m}^{- 2} {\hat{Σ}}_{m} X^{T} Y_{m}

(6)

{\hat{Σ}}_{m} = {(σ_{m}^{- 2} X^{T} X + A_{m})}^{- 1}

(7)

The posterior mean vector ${\hat{b}}_{m}$ and covariance matrix ${\hat{Σ}}_{m}$ , as well as the relevancy of each feature in the design matrix, continuously depends on the hyperparameter vector $α_{m}$ . This conclusion can be validated by considering two extreme cases for the value of the ith component of the hyperparameter vector $α_{m, i}$ . If $α_{m, i} \to 0$ , then the prior PDF $p (b_{m, i} | α_{m, i})$ is Gaussian distribution with very small precision (very large variance), corresponding to the case that the posterior PDF is dominated by the likelihood function so $b_{m, i} \to b_{m, i}^{mle}$ where $b_{m, i}^{mle}$ is the maximum likelihood estimate of $b_{m, i}$ , and its associated feature $X_{i}$ is retained in the model for modal frequency-ambient condition pattern of the mth mode in equation (2). On the other hand, if $α_{m, i} \to \infty$ , then the prior PDF $p (b_{m, i} | α_{m, i})$ is Gaussian distribution with very high precision (i.e. very small variance), corresponding to the case that the posterior PDF is dominated by the prior PDF so $b_{m, i} \to 0$ and its associated feature $X_{i}$ is irrelevant to the model in equation (2). In the remaining of this section, aiming for searching the optimal model class with the relevance features for each vibrational mode of the structure, an efficient strategy is proposed for optimization of the hyperparameter vector $α_{m}$ along with the prediction-error variance $σ_{m}^{2}$ .

Optimization of hyperparameter vector and prediction-error variance

According to Bayes’ theorem, the posterior PDF of the hyperparameter vector $α_{m}$ and the prediction-error variance $σ_{m}^{2}$ is given by $p (α_{m}, σ_{m}^{2} | D_{m}) \propto p (D_{m} | α_{m}, σ_{m}^{2}) p (α_{m}, σ_{m}^{2})$ , where $p (D_{m} | α_{m}, σ_{m}^{2})$ is the evidence conditional on $α_{m}$ and $σ_{m}^{2}$ , and $p (α_{m}, σ_{m}^{2})$ is the prior PDF of $α_{m}$ and $σ_{m}^{2}$ . As the maximum a posteriori estimation can be well approximated by the maximum evidence estimation, the optimal values for $α_{m}$ and $σ_{m}^{2}$ can be determined by considering the following optimization⁵⁷

[{\hat{α}}_{m}, {\hat{σ}}_{m}^{2}] = \underset{α_{m} \in R^{N_{b}}, σ_{m}^{2} \in R^{+}}{\arg max} L_{m} (α_{m}, σ_{m}^{2})

(8)

where $L_{m} (α_{m}, σ_{m}^{2}) \equiv \ln p (D_{m} | α_{m}, σ_{m}^{2})$ is the logarithm of the evidence of the mth vibrational mode, which can be evaluated as follows⁵⁷

\begin{matrix} L_{m} (α_{m}, σ_{m}^{2}) = \ln \int p (D_{m} | α_{m}, σ_{m}^{2}) p (b_{m} | α_{m}) d b_{m} \\ = - \frac{1}{2} [N \ln 2 π + \ln | C_{m} | + Y_{m}^{T} C_{m}^{- 1} Y_{m}] \end{matrix}

(9)

with the matrix $C_{m} \in R^{N \times N}$ given by

C_{m} = σ_{m}^{2} I_{N} + X^{T} A_{m}^{- 1} X

(10)

where $I_{N}$ is the $N \times N$ identity matrix. The optimization in equation (8) is $(N_{b} + 1)$ -dimensional, which is infeasible when the number of components of the feature vector $N_{b}$ is large. Here, separation optimization strategy is proposed for efficiency enhancement. Consider the following decomposition for the matrix $C_{m}$

\begin{matrix} C_{m} = σ_{m}^{2} I_{N} + \sum_{k = 1}^{N_{b}} α_{m, k}^{- 1} X_{k} X_{k}^{T} \\ = C_{m, - i} + α_{m, i}^{- 1} X_{i} X_{i}^{T} \end{matrix}

(11)

where $C_{m, - i} = σ_{m}^{2} I_{N} + \sum_{k \neq i} α_{m, k}^{- 1} X_{k} X_{k}^{T}$ represents $C_{m}$ excluding the ith term in the sum $α_{m, i}^{- 1} X_{i} X_{i}^{T}$ . Then, $| C_{m} |$ and $C_{m}^{- 1}$ can be rewritten as follows⁵⁸

| C_{m} | = | C_{m, - i} | | 1 + α_{m, i}^{- 1} X_{i} C_{m, - i}^{- 1} X_{i}^{T} |

(12)

C_{m}^{- 1} = C_{m, - i}^{- 1} - \frac{C_{m, - i}^{- 1} X_{i} X_{i}^{T} C_{m, - i}^{- 1}}{α_{m, i} + X_{i}^{T} C_{m, - i}^{- 1} X_{i}}

(13)

Let $α_{m, - i}$ denote the hyperparameter vector $α_{m}$ excluding $α_{m, i}$ . By plugging equations (12) and (13) into equations (9) and (10), $L_{m} (α_{m}, σ_{m}^{2})$ can be decomposed into two parts

L_{m} (α_{m}, σ_{m}^{2}) = L_{m} (α_{m, - i}, σ_{m}^{2}) + l_{m} (α_{m, i}, σ_{m}^{2})

(14)

The first term $L_{m} (α_{m, - i}, σ_{m}^{2})$ is the log-evidence removing $α_{m, i}$ (and thus $b_{m, i}$ and $X_{i}$ ) from the model of the mth mode

L_{m} (α_{m, - i}, σ_{m}^{2}) = - \frac{1}{2} [N \ln 2 π + \ln | C_{m, - i}^{- 1} | + Y_{m}^{T} C_{m, - i}^{- 1} Y_{m}]

(15)

The second term $l (α_{m, i}, σ_{m}^{2})$ , related to $α_{m, i}$ and $σ_{m}^{2}$ , is given by

l_{m} (α_{m, i}, σ_{m}^{2}) = \frac{1}{2} [\ln α_{m, i} - \ln (α_{m, i} + S_{m, i}) + \frac{Q_{m, i}^{2}}{α_{m, i} + S_{m, i}}]

(16)

with

S_{m, i} = X_{i}^{T} C_{m, - i}^{- 1} X_{i}

(17)

Q_{m, i} = X_{i}^{T} C_{m, - i}^{- 1} Y_{m}

(18)

As the ith hyperparameter $α_{m, i}$ has been isolated from the log-evidence $L_{m} (α_{m}, σ_{m}^{2})$ , the optimal value of $α_{m, i}$ can be obtained by taking the partial derivative of $l_{m} (α_{m, i}, σ_{m}^{2})$ in equation (16) with respect to $α_{m, i}$

\frac{\partial l_{m} (α_{m, i}, σ_{m}^{2})}{\partial α_{m, i}} = \frac{α_{m, i}^{- 1} S_{m, i}^{2} - (Q_{m, i}^{2} - S_{m, i})}{2 {(α_{m, i} + S_{m, i})}^{2}}

(19)

Thus, the optimal value of the ith hyperparameter ${\hat{α}}_{m, i}$ is achieved when $\partial l_{m} (α_{m, i}, σ_{m}^{2}) / \partial α_{m, i} = 0$ under the constrain that $α_{m, i} \geq 0$ (because $α_{m, i}$ is the diagonal element of the precision matrix $A_{m}$ ). Thus, there are two possible values for ${\hat{α}}_{m, i}$

{\hat{α}}_{m, i} = {\begin{matrix} \infty, & if Q_{m, i}^{2} \leq S_{m, i} \\ S_{m, i}^{2} / (Q_{m, i}^{2} - S_{m, i}), & if Q_{m, i}^{2} > S_{m, i} \end{matrix}

(20)

Based on the updated value of $α_{m}$ , the optimal prediction-error variance is re-estimated by considering the derivative of $L_{m} (α_{m}, σ_{m}^{2})$ with respect to $σ_{m}^{2}$ ^53,54

{\hat{σ}}_{m}^{2} = \frac{{(Y_{m} - X {\hat{b}}_{m})}^{T} (Y_{m} - X {\hat{b}}_{m})}{N - N_{b} + \sum_{i = 1}^{N_{b}} {\hat{α}}_{m, i} {\hat{Σ}}_{m, i}}

(21)

where ${\hat{b}}_{m}$ is obtained by equation (6); ${\hat{Σ}}_{m}$ is obtained by equation (7) and ${\hat{Σ}}_{m, i}$ is the ith diagonal element of the posterior covariance matrix ${\hat{Σ}}_{m}$ .

Summary of the proposed approach

The procedures of the proposed algorithm are shown in Figure 4 and are as follows. For pattern recognition of the mth mode, the algorithm starts from a single feature and proceeds by adding, pertaining and/or deleting features.

Initialize $σ_{m}^{2}$ . The default value is 20% of the sample variance of the mth mode output $Y_{m} = [ω_{m, 1}^{2}, \dots, ω_{m, N}^{2}]^{T}$ .

Select the initial feature $X_{i}$ . The correlation coefficients between the output $Y_{m}$ and each feature are calculated. The initial feature is selected as the one which possesses the maximum absolute correlation coefficient.

Initialize $α_{m}$ . Since only one of the initial feature $X_{i}$ is included, each component of $α_{m, - i}$ is infinity and, hence, $C_{m, - i} = σ_{m}^{2} I_{N}$ . The hyperparameter associated with the initial feature $X_{i}$ can be estimated using equation (20)

α_{m, i} = \frac{S_{m, i}^{2}}{Q_{m, i}^{2} - S_{m, i}} = \frac{X_{i}^{T} X_{i}}{{(X_{i}^{T} Y_{m})}^{2} / (X_{i}^{T} X_{i}) - σ_{m}^{2}}

(22)

Compute ${\hat{b}}_{m} = σ_{m}^{- 2} {\hat{Σ}}_{m} X^{T} Y_{m}$ and ${\hat{Σ}}_{m} = (σ_{m}^{- 2} X^{T} X + A_{m})^{- 1}$ by equations (6) and (7) based on $A_{m}$ and $σ_{m}^{2}$ .

Update $S_{m, j} = X_{j}^{T} C_{m, - j}^{- 1} X_{j}$ and $Q_{m, j} = X_{j}^{T} C_{m, - j}^{- 1} Y_{m}$ by equations (17) and (18) for $j = 1, \dots, N_{b}$ . Select a candidate feature $X_{j}$ , update its associated hyperparameter $α_{m, j}$ as follows

If $Q_{m, j}^{2} \leq S_{m, j}$ , then remove $X_{j}$ and set $α_{m, j} = \infty$ ;

If $Q_{m, j}^{2} > S_{m, j}$ and $α_{m, j} < \infty$ , then update $α_{m, j} = S_{m, j}^{2} / (Q_{m, j}^{2} - S_{m, j})$ (implying that $X_{j}$ has been already included in the model in current iteration);

If $Q_{m, j}^{2} > S_{m, j}$ and $α_{m, j} = \infty$ , then include $X_{j}$ and update $α_{m, j} = S_{m, j}^{2} / (Q_{m, j}^{2} - S_{m, j})$ .

Update ${\hat{σ}}_{m}^{2}$ by equation (21) with the updated ${\hat{Σ}}_{m}$ and ${\hat{b}}_{m}$ .

Repeat Steps 4–6 until convergence achieves so the final learning results are obtained.

Figure 4.

The procedures of the proposed algorithm.

Modal frequency-ambient condition pattern recognition

Table 1 shows the weight parameter identification results of the optimal models of the three modes, the optimal values (outside the parentheses), and the standard deviations (inside the parentheses), and the sign “–” indicates that the weight parameter is not included in the optimal model, representing that the associated feature is irrelevant. It is observed that the optimal models of the three modes possess different relevant features shown as follows

\begin{matrix} ω_{1}^{2} = 1.92 + 0.04 H_{0} - 0.931 T_{0}^{2} H_{0} \\ + 1.014 T_{0}^{3} H_{0} + 0.255 T_{0}^{2} H_{0}^{2} \end{matrix}

(23)

ω_{2}^{2} = 2.582 + 0.266 T_{0} - 0.187 T_{0}^{2} H_{0} + 0.358 T_{0}^{2} H_{0}^{2}

(24)

\begin{matrix} ω_{3}^{2} = 3.332 - 0.033 T_{0} H_{0} - 0.931 T_{0}^{2} H_{0} \\ - 0.015 T_{0}^{4} + 0.332 T_{0}^{2} H_{0}^{2} \end{matrix}

(25)

Table 1.

Weight parameter identification results of the optimal models of the three modes.

Mode	$b_{m, 1}$	$b_{m, 2}$	$b_{m, 3}$	$b_{m, 4}$	$b_{m, 5}$	$b_{m, 6}$	$b_{m, 7}$	$b_{m, 8}$
m = 1	1.920 (0.015)	−	0.040 (0.037)	−	−	−	−	−0.931 (0.127)
m = 2	2.582 (0.025)	0.226 (0.069)	−	−	−	−	−	−0.187 (0.126)
m = 3	3.332 (0.013)	−	−	−	−0.033 (0.057)	−	−	−
Mode	$b_{m, 9}$	$b_{m, 10}$	$b_{m, 11}$	$b_{m, 12}$	$b_{m, 13}$	$b_{m, 14}$	$b_{m, 15}$	$σ_{m}^{2} (\times 10^{- 4})$
m = 1	−	−	−	1.014 (0.144)	0.255 (0.073)	−	−	8.68
m = 2	−	−	−	−	0.358 (0.087)	−	−	6.02
m = 3	−	−	−0.015 (0.012)	−	0.332 (0.061)	−	−	7.34

– indicates that the weight parameter is not included in the optimal model.

Both the temperature and relative humidity exist in all the three optimal models, showing that these two ambient factors influence the modal frequencies of different modes and they should be considered simultaneously in modal frequencies prediction. This conclusion is consistent with that by Xia et al.⁵ and Yuen and Kuok.^11,12 On the other hand, the optimal models for different modes possess different combinations of the two ambient factors. The reasons are given as follows. In practical monitoring of a target structure under changing ambient and operational conditions, fluctuations of modal frequencies are induced by multiple environmental factors, not only temperature and relative humidity but also other factors, such as wind load,⁵⁹ rainfall, support reaction,⁶⁰ and soil–structure interaction,^61,62 and the degrees of participation of different environmental factors to the modal parameters of different modes are possibly not identical. On the other hand, it is worth noting that although only the temperature and relative humidity are utilized in the study, the temperature and relative humidity are the primary effects while other factors (wind load, rainfall, support reaction, soil–structure interaction) are the secondary effects in modal frequencies prediction. This conclusion is supported by checking the normalized residual plot of the three modes shown in Figure 5 that most of the normalized residuals satisfy $| {\bar{ε}}_{m} | \leq 1.96, m = 1, 2, 3$ , where the threshold “1.96” corresponds to that there is 95% in probability for a standard normal distributed sample falling into the interval [−1.96, 1.96]. Therefore, it can be concluded that the optimal models given in equations (23)–(25) are capable to properly capture the wandering pattern of the modal frequencies under normal weather conditions.

Figure 5.

Normalized residual plot of the three modes ( ${\bar{ε}}_{1}, {\bar{ε}}_{2}, and {\bar{ε}}_{3}$ are the normalized residuals of modes 1, 2, and 3, respectively).

Note that in Bayesian spectral density approach,^46–48 not only identified values but also associated uncertainties of modal frequencies are available. It is worth noting that unequal uncertainty values for modal frequencies of different days may lead to heterogeneous error pattern, requiring the development of the heterogeneous Bayesian inference and learning framework.⁶³

Figures 6 and 7 show the squared modal frequencies versus the temperature and relative humidity of the three modes, respectively. The modal frequency–temperature curves of Figure 6 are plotted with the relative humidity evaluated at its sample mean (67.9%), and the modal frequency–relative humidity curves of Figure 7 are plotted with the temperature with its sample mean (24.7°C). The increasing trends of these curves satisfy the correlation analysis results of Figure 3 for different modes. Finally, the design modal frequency contours of the optimal values and the associated uncertainties for the three modes are given in Figure 8. The left three subplots are for the optimal values of modal frequencies given the temperature and the relative humidity for modes 1, 2, and 3, respectively, and the right three subplots are for the corresponding associated uncertainties for modes 1, 2, and 3, respectively. Note that different modes associate with different uncertainty contours because they possess different optimal models. In addition, it can be observed that the uncertainty level of the corner regions in the uncertainty contours is higher than that of the central region. This is due to the non-uniform distribution of the input measurement that most of the measured data points locate in the center region while very few locate in the corner region, especially for the top right and bottom right regions. On the other hand, the non-uniform distribution of the input measurement also reflects the ambient condition pattern of the city for the monitored structure that, high temperature is with intermediate level of relative humidity rather than high or low level of relative humidity.

Figure 6.

The squared modal frequencies versus the temperature of the three modes.

Figure 7.

The squared modal frequencies versus the relative humidity of the three modes.

Figure 8.

Design modal frequency contours of the optimal values (left three subplots) and the associated uncertainties (right three subplots) for the three modes.

Conclusion

In this article, a novel machine-learning algorithm is proposed to automatically select relevance features in modal frequency-ambient condition pattern recognition based on structural dynamic response and ambient condition measurement. In the traditional feature selection approaches, 3 × (2¹⁵-1) candidates in total are required to obtain the optimal model classes for the three modes based on the 15-component full extracted feature vector, which requires large computational effort. On the other hand, in the proposed approach, the optimal model classes for modal frequency predictions of different vibrational modes of the structure are determined by first introduction of the hyperparameter in the ARD prior for the weight parameter vector, which controls the relevancy of different features in the prediction model, then conducting hyperparameter optimization by maximum evidence estimation. During the optimization process, the irrelevance components of the weight parameter vector as well as the corresponding features are automatically pruned out while the relevance terms are retained. The proposed algorithm is utilized for structural health assessment for a reinforced concrete building based on 1-year daily measurements. It turns out that the optimal model classes of different vibrational modes of the structure are capable to capture the modal frequency-ambient condition patterns of different modes. In addition, the design modal frequency contours of the optimal values and the associated uncertainties are given, which provides an efficient tool for modal frequency prediction in future analysis.

Footnotes

Academic Editor: Jun Li

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (51508201), the China Postdoctoral Science Foundation, and the State Key Laboratory of Subtropical Building Science, South China University of Technology (2016ZB26).

References

Farrar

Worden

An introduction to structural health monitoring. Philos T R Soc A 2007; 365: 303–315.

Sohn

Dzwonczyk

Straser

. An experimental study of temperature effect on modal parameters of the Alamosa Canyon Bridge. Earthq Eng Struct D 1999; 28: 879–897.

Bamu

Zingoni

Damage, deterioration and the long-term structural performance of cooling-tower shells: a survey of developments over the past 50 years. Eng Struct 2005; 27: 1794–1800.

Lee

Yun

CB.

Damage diagnosis of steel girder bridges using ambient vibration data. Eng Struct 2006; 28: 912–925.

Xia

Hao

Zanardo

. Long term vibration monitoring of an RC slab: temperature and humidity effect. Eng Struct 2006; 28: 441–452.

Kim

Park

Lee

BJ.

Vibration-based damage monitoring in model plate-girder bridges under uncertain temperature conditions. Eng Struct 2007; 29: 1354–1365.

Liu

DeWolf

JT.

Effect of temperature on modal variability of a curved concrete bridge under ambient loads. J Struct Eng 2007; 133: 1742–1751.

Catbas

Susoy

Frangopol

DM.

Structural health monitoring and reliability estimation: long span truss bridge application with environmental monitoring data. Eng Struct 2008; 30: 2347–2359.

Deraemaeker

Reynders

De Roeck

. Vibration-based structural health monitoring using output-only measurements under changing environment. Mech Syst Signal Pr 2008; 22: 34–56.

10.

Clinton

Bradford

Heaton

. The observed wander of the natural frequencies in a structure. B Seismol Soc Am 2006; 96: 237–257.

11.

Yuen

Kuok

SC.

Ambient interference in long-term monitoring of buildings. Eng Struct 2010; 32: 2379–2386.

12.

Yuen

Kuok

SC.

Modeling of environmental influence in structural health assessment for reinforced concrete buildings. Earthq Eng Eng Vib 2010; 9: 295–306.

13.

Chen

Lin

. Theoretical and experimental modal analysis of the Guangzhou New TV Tower. Eng Struct 2011; 33: 3628–3646.

14.

Bao

Xia

. Data fusion-based structural damage detection under varying temperature conditions. Int J Struct Stab Dy 2012; 12: 1250052.

15.

Xia

Chen

Weng

. Temperature effect on vibration properties of civil structures: a literature review and case studies. J Civil Struct Health Monit 2012; 2: 29–46.

16.

Worden

Farrar

Manson

. The fundamental axioms of structural health monitoring. Proc Roy Soc Lond A Mat 2007; 463: 1639–1664.

17.

Farrar

Worden

Structural health monitoring: a machine learning perspective. Hoboken, NJ: John Wiley & Sons, 2012.

18.

Yan

Kerschen

De Boe

. Structural damage diagnosis under varying environmental conditions—part I: a linear analysis. Mech Syst Signal Pr 2005; 19: 847–864.

19.

Hua

Fan

. Correlating modal properties with temperature using long-term monitoring data and support vector machine technique. Eng Struct 2005; 27: 1762–1773.

20.

. Modal identification of bridges under varying environmental conditions: temperature and wind effects. Struct Control Hlth 2010; 17: 495–512.

21.

Mechbal

Uribe

Rebillat

A probabilistic multi-class classifier for structural health monitoring. Mech Syst Signal Pr 2015; 60: 106–123.

22.

Santos

Figueiredo

Silva

MFM

. Machine learning algorithms for damage detection: kernel-based approaches. J Sound Vib 2016; 363: 584–599.

23.

Beck

JL.

Bayesian system identification based on probability logic. Struct Control Hlth 2010; 17: 825–847.

24.

Lam

Yin

Statistical detection of multiple cracks on thin plates utilizing dynamic response. Eng Struct 2010; 32: 3145–3152.

25.

Yin

Lam

Chow

HM.

A Bayesian probabilistic approach for crack characterization in plate structures. Comput-Aided Civ Inf 2010; 25: 375–386.

26.

Yuen

HQ.

Real-time system identification: an algorithm for simultaneous model class selection and parametric identification. Comput-Aided Civ Inf 2015; 30: 785–801.

27.

Zhao

XY.

Noise smoothing for structural vibration test signals using an improved wavelet thresholding technique. Sensors 2012; 12: 11205–11220.

28.

Thatoi

Das

Parhi

DR.

Review of techniques for fault diagnosis in damaged structure and engineering system. Adv Mech Eng 2012; 4: 327569.

29.

Błachowski

Zhong

. Rank-revealing QR decomposition applied to damage localization in truss structures. Struct Control Hlth. Epub ahead of print 18 February 2016. DOI: 10.1002/stc.1849.

30.

Wavelet based multi-step filtering method for bridge health monitoring using GPS and accelerometer. Smart Struct Syst 2013; 11: 331–348.

31.

Błachowski

A degree of dispersion-based damage localization method. Struct Control Hlth 2016; 23: 176–192.

32.

Spencer

A test method for damage diagnosis of suspension bridge suspender cables. Comput-Aided Civ Inf 2015; 30: 771–784.

33.

Yuen

Kuok

SC.

Bayesian methods for updating dynamic models. Appl Mech Rev 2011; 64: 010802.

34.

Yin

Lam

HF.

Dynamic analysis of finite-length circular cylindrical shells with a circumferential surface crack. J Eng Mech 2012; 139: 1419–1434.

35.

Lam

Peng

SK.

Development of a practical algorithm for Bayesian model updating of a coupled slab system utilizing field test data. Eng Struct 2014; 79: 182–194.

36.

Lam

Wong

MT.

The Bayesian methodology for the detection of railway ballast damage under a concrete sleeper. Eng Struct 2014; 81: 289–301.

37.

Yuen

Katafygiotis

LS.

Modal decomposition using multi-channel response measurements. Probabilist Eng Mech 2014; 37: 60–73.

38.

Kuok

Yuen

KV.

Investigation of modal identification and modal identifiability of a cable-stayed bridge with Bayesian framework. Smart Struct Syst 2016; 17: 445–470.

39.

Yuen

Kuok

SC.

Online updating and uncertainty quantification using nonstationary output-only measurement. Mech Syst Signal Pr 2016; 66: 62–77.

40.

Lei

Lin

A decentralized control algorithm for large-scale building structures. Comput-Aided Civ Inf 2012; 27: 2–13.

41.

Lei

Chen

Zhou

An algorithm based on two-step Kalman filter for intelligent structural damage detection. Struct Control Hlth 2015; 22: 694–706.

42.

Yuen

KV.

Novel outlier-resistant extended Kalman filter for robust online structural identification. J Eng Mech 2014; 141: 04014100.

43.

Yuen

Kuok

SC.

Efficient Bayesian sensor placement algorithm for structural identification: a general approach for multi-type sensory systems. Earthq Eng Struct D 2015; 44: 757–774.

44.

Yuen

HQ.

Peak ground acceleration estimation by linear and nonlinear models with reduced order Monte Carlo simulation. Comput-Aided Civ Inf 2011; 26: 30–47.

45.

Yuen

KV.

Seismic attenuation relationship with homogeneous and heterogeneous prediction-error variance models. Earthq Eng Eng Vib 2014; 13: 1–11.

46.

Yuen

Katafygiotis

Beck

JL.

Spectral density estimation of stochastic vector processes. Probabilist Eng Mech 2002; 17: 265–272.

47.

Yuen

Beck

JL.

Updating properties of nonlinear dynamical systems with uncertain input. J Eng Mech 2003; 129: 9–20.

48.

Yuen

KV.

Bayesian methods for structural dynamics and civil engineering. Hoboken, NJ: John Wiley & Sons, 2010.

49.

Watson

Rajapakse

RKND

. Seasonal variation in material properties of a flexible pavement. Can J Civil Eng 2000; 27: 44–54.

50.

Rincón

Hernández-López

de Valle-Moreno

. Environmental influence on point anodes performance in reinforced concrete. Constr Build Mater 2008; 22: 494–503.

51.

Moser

Moaveni

Environmental effects on the identified natural frequencies of the Dowling Hall Footbridge. Mech Syst Signal Pr 2011; 25: 2336–2357.

52.

Moaveni

Behmanesh

Effects of changing ambient temperature on finite element model updating of the Dowling Hall Footbridge. Eng Struct 2012; 43: 58–68.

53.

MacKay

DJC

. Bayesian interpolation. Neural Comput 1992; 4: 415–447.

54.

MacKay

DJC

. Bayesian nonlinear modeling for the prediction competition. ASHRAE Tran 1994; 100: 1053–1062.

55.

Tipping

ME.

Bayesian inference: an introduction to principles and practice in machine learning. In: Bousquet

von Luxburg

Rätsch

(eds) Advanced lectures on machine learning. Berlin: Springer, 2004, pp.41–62.

56.

Yuen

HQ.

A novel probabilistic method for robust parametric identification and outlier detection. Probabilist Eng Mech 2012; 30: 48–59.

57.

Bishop

CM.

Pattern recognition and machine learning. Berlin: Springer, 2006.

58.

Harville

DA.

Matrix algebra from a statistician’s perspective. Berlin: Springer, 1997.

59.

Siringoringo

Fujino

System identification of suspension bridge from ambient vibration response. Eng Struct 2008; 30: 462–477.

60.

Kulprapha

Warnitchai

Structural health monitoring of continuous prestressed concrete bridges using ambient thermal responses. Eng Struct 2012; 40: 20–38.

61.

Stewart

Fenves

GL.

System identification for evaluating soil–structure interaction effects in buildings from strong motion recordings. Earthq Eng Struct D 1998; 27: 869–885.

62.

Todorovska

MI.

Soil-structure system identification of Millikan Library North–South response during four earthquakes (1970–2002): what caused the observed wandering of the system frequencies?

B Seismol Soc Am 2009; 99: 626–635.

63.

Yuen

KV.

Ground motion prediction equation development by heterogeneous Bayesian learning. Comput-Aided Civ Inf. Epub ahead of print 13 June 2016. DOI: 10.1111/mice.12215.