Next Article in Journal
Numerical Analysis and Experimental Verification of Damage Identification Metrics for Smart Beam with MFC Elements to Support Structural Health Monitoring
Next Article in Special Issue
Integration with 3D Visualization and IoT-Based Sensors for Real-Time Structural Health Monitoring
Previous Article in Journal
Beamforming Using Exact Evaluation of Leakage and Ergodic Capacity of MU-MIMO System
Previous Article in Special Issue
Remaining Useful Life Estimation for Engineered Systems Operating under Uncertainty with Causal GraphNets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structural Response Prediction for Damage Identification Using Wavelet Spectra in Convolutional Neural Network

by
Edisson Alberto Moscoso Alcantara
,
Michelle Diana Bong
and
Taiki Saito
*
Department of Architecture and Civil Engineering, Toyohashi University of Technology, Aichi 441-8580, Japan
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(20), 6795; https://doi.org/10.3390/s21206795
Submission received: 31 August 2021 / Revised: 9 October 2021 / Accepted: 10 October 2021 / Published: 13 October 2021
(This article belongs to the Special Issue Sensors in Structural Health Monitoring and Smart Structural Control)

Abstract

:
If damage to a building caused by an earthquake is not detected immediately, the opportunity to decide on quick action, such as evacuating the building, is lost. For this reason, it is necessary to develop modern technologies that can quickly obtain the structural safety condition of buildings after an earthquake in order to resume economic and social activities and mitigate future damage by aftershocks. A methodology for the prediction of damage identification is proposed in this study. Using the wavelet spectrum of the absolute acceleration record measured by a single accelerometer located on the upper floor of a building as input data, a CNN model is trained to predict the damage information of the building. The maximum ductility factor, inter-story drift ratio, and maximum response acceleration of each floor are predicted as the damage information, and their accuracy is verified by comparing with the results of seismic response analysis using actual earthquakes. Finally, when an earthquake occurs, the proposed methodology enables immediate action by revealing the damage status of the building from the accelerometer observation records.

1. Introduction

Earthquakes in the proximity of structurally vulnerable buildings could cause damage of varying intensities. Damage of different risk levels is often difficult to classify rapidly, making it difficult to accurately determine the structural safety of a building. For example, according to the National Institute of Civil Defense of Peru, during the Pisco earthquake on 15 August 2007, in the five main regions of Peru (Lima is included, which is the capital of Peru), 136,149 dwellings, 1278 educational buildings, and 126 health buildings collapsed or were damaged, and their use was classified as restricted or unsafe [1]. However, this report was released almost two months after the earthquake, during which time, all activities in the affected areas of the main regions had to be suspended, including the construction of temporary dwellings.
Resilient cities are goals that countries are building towards to increase the capacity for learning from past disasters for better future protection and to improve risk reduction measures [2]. In particular, as part of this concept, there is a need to develop modern structures for which we can quickly obtain the structural safety information after an earthquake for resuming economic and social activities in order to minimize social disruption and mitigate the effects of future earthquakes [3]. In order to promote and disseminate knowledge to increase social resilience and reduce earthquake risk, experts from academia and industry gathered in 2019 for a workshop focused on state-of-the-art risk-reduction strategies. It identified a need for research in the area of structural health monitoring (SHM) to assess the integrity and performance of engineering structures in order to quickly detect damage after an earthquake and enable decision making [4]. SHM is a field where it is possible to obtain the real-time structural responses and successful fast post-earthquake damage detection of monitored buildings, bridges, cultural heritage structures, dams, base-isolated buildings, etc. [5,6]. For instance, Goulet et al. proposed a methodology that updates the prediction of the damage state of uninspected monitored buildings as the model learns from collected data of the damage state of inspected buildings [7]. This proposal was validated in a city with 1000 buildings. Furthermore, Sivasuriyan et al. reviewed a large number of studies on the practical implementation and operations of SHM in multi-story buildings, as well as damage evaluation of monitored buildings, and discussed the structural response by considering static and dynamic analysis using numerical simulations such as finite element analysis (FEA) [8].
In the field of SHM, there are several types of sensors to measure and diagnose the static and dynamic properties of the monitored buildings. Antunez et al. demonstrated that optical fiber sensors can be useful in the static and dynamic monitoring of large raw earth masonry structures common in cultural, historical, and architecturally recognized buildings around the world [9]. Piezoelectric sensors are another type of monitoring device, and Roghaei et al. proposed a method to identify stress and deformation using an array of sensors mounted in certain locations [10]. They verified the proposed method using a three-story steel building and confirmed that continuous monitoring and analysis of sensor signals can help the building manager to apply warning alarms and call for evacuation. However, the most common monitoring control sensor is the accelerometer. For instance, Wang et al. developed a method to evaluate the story damage index (SDI) based on the modal frequency and mode shape obtained from the records of earthquake response of a building [11]. Furthermore, an approximate story damage index (ASDI) was developed without considering the information of the floor mass to identify the extent of damage to the story. Although it was possible to verify the damage index by some numerical simulations and the experimental data analysis established previously, it was necessary to calculate the modal frequency and mode shapes from the post-earthquake structural responses of each story and to compare with the values of the building before the earthquake. It is worth pointing out that a large number of sensors will require a high investment. For this reason, Xu et al. estimated the maximum drift and time histories of relative displacement in all stories of multi-degree-of-freedom (MDOF) structures considering only one accelerometer, verifying the effectiveness of the method by taking into account the robustness, installation location, and truncation error [12].
The machine learning method, which predicts the structural responses using a learning model specific to the structure, may provide higher accuracy by updating the model after each earthquake. According to study [13], there are two approaches for damage identification: model-driven methods and data-driven methods. In a model-driven approach, usually, a high-fidelity physical model of the structure is used to establish a comparison metric between the model and the measured data from the real structure to distinguish the damage condition from the normal condition. In a data-driven approach, a structural model is used as a statistical representation of the system, and the main algorithms developed for this purpose are those in the field of pattern recognition or, more broadly, machine learning. A convolutional neural network (CNN) is a tool for solving the problem of pattern recognition related to image and video recognition, classification, natural language processing, and others. Oh et al. studied a method of predicting the time histories of displacement of building structures from the measured acceleration responses on each floor based on a CNN, considering that the time series of acceleration structural response is similar to pixel-based image data (every acceleration value corresponds to one pixel), which is the basic input data in CNN [14]. The validation of their proposed method was from a numerical process using the ASCE benchmark model and an experimental test on a reinforced concrete (RC) frame structure. However, the structural model and dynamic responses used in the studies exhibited linear behavior. Tsuchimoto et al. proposed a rapid safety evaluation of multi-story buildings using sparse acceleration measurements [15]. Their proposed method predicts the maximum story drift ratio, and ultimately classifies the damage into three classes, namely “Safe”, “Restricted Use”, and “Unsafe” from a damage-sensitive feature (comparison between linear and nonlinear acceleration measurement responses) and ground acceleration as input data. Subsequently, Tsuchimoto et al. modified the previous method for high-rise buildings and validated considering an experimental test of a large-scale structure (1/3-scale 18-story steel building tested on the shaking table at E-Defense in Japan) [16].
There are two main characteristics observed on the ground motion records due to earthquakes. The first is the non-stationary characteristics in which the intensity of the ground motion varies with time; they are represented by the acceleration, velocity, and displacement. The second is the non-stationary characteristics in which the frequency content of the ground motion varies with time; they depend on several parameters such as magnitude, source and path effects, local site conditions, etc. [17]. Time–frequency distribution analysis is a method of obtaining a two-dimensional spectral function (there are several types of functions according to resources and needs) from a one-dimensional signal (ground motion or time–history structural response) that reflects the time and frequency of the original signal and is suitable to analyze the changes in the linear and nonlinear structural responses with only one function. For instance, Tao et al. used the matching pursuit decomposition algorithm to analyze the time–frequency distribution of the ground motion and verify the effect on the dynamic response of a nonlinear structure, and finally, this method reveals the effect of the ground motion on the nonlinear structural response [18]. Moreover, Cao et al. demonstrated the effect of energy concentration on the structural nonlinear response by using the wavelet transform to obtain a local spectrum and change the energy distribution over time for several earthquake records [19]. Spanos et al. analyzed the undamaged and damaged condition of a 20-story steel frame building using the harmonic wavelet transform applied to structural responses to obtain the variation of the effective natural frequencies due to the influence of the nonlinearity developed during the seismic event [20]. Balafas and Kiremidjian used the continuous wavelet transform of the input and output acceleration measurements to extract damage sensitive features for seismic damage estimation in civil structures [21]. Noh et al. proposed an extraction method of three damage-sensitive features using wavelet transform spectrum for structural damage diagnosis and applied them to experimental data of a reinforced concrete bridge column and a four-story steel moment-resisting frame structure [22]. In general, time–frequency distributions are two-dimensional spectral functions that can be used as input data for a CNN to predict dynamic issues related to structural engineering. For example, Xu et al. proposed a methodology to recognize and classify different types of vibrational events (digging, walking, vehicles passing, and damaging) [23]. First, they denoise the unknown signal and use the short-time Fourier transform (STFT) to obtain the time–frequency spectra and input them to the CNN for automatic feature extraction and classification. The proposed method used the support vector machine method to compare the obtained recognition rates of vibration events over 90% with the previous soft-max classifier. Dokht et al. used a CNN and STFT to consider a dataset of over 4900 earthquakes recorded over 3 years in Canada to classify between earthquake and noise signals. They also used another CNN and wavelet spectrum to classify and separate P from S waves and estimate their approximate arrival times [24]. Their results achieved an average accuracy of nearly 99% for both networks. Mousavi et al. proposed a detector based on a deep neural network (CNN belong to this field) called CNN-RNN Earthquake Detector (CRED), which is a network that combines a CNN and a recurrent neural network (RNN), specifically the bidirectional long-short-term-memory (LSTM) method, to learn the time-frequency characteristics of the dominant phases in an earthquake signal from three-component data recorded at a single station, having an accuracy of 99.95% [25]. In addition, Liao et al. proposed an identification method for a structural seismic response using a wavelet spectrum as input data in a CNN to distinguish the responses during an earthquake event under serviceability conditions [26]. Linear and nonlinear behaviors are considered in the research. According to previous studies, the CNN method in the SHM field has advantages over other methods in terms of higher accuracy by updating the model after each earthquake, flexibility to combine different methodologies, wide application areas, etc., however, it requires a large database of known data to train the model.
Previous studies have not fully investigated how to define the damage level of each floor of a structure from the time–frequency distribution of the observation data of a single sensor. The Japan Structural Consultants Association (JSCA), an organization of building structural engineers in Japan, uses three parameters of safety criteria used on the assessment of a building: absolute acceleration, ductility ratio, and story drift ratio [27]. Acceleration is related to damage in nonstructural components, and ductility and story drift ratio are related to damage in structural components. It is worth pointing out that the use of only one sensor implies a low-cost investment. This study proposes a methodology to predict the absolute acceleration, ductility ratio, and story drift ratio on each floor under earthquake conditions using machine learning. In the beginning, the earthquake responses of a model building are calculated under the scaled earthquake records with several intensities (scale factors). The level of intensity is established to obtain a range of linear and nonlinear behavior of the building. Then, wavelet spectra are developed from the structural response accelerations on the upper floor of the building. The wavelet spectra are the input data of a CNN model to predict the absolute acceleration, ductility ratio, and story drift ratio on each floor, which correspond to the damage of the nonstructural and structural components of the building.
This paper contains sections as follows: In Section 2, the basis and methods of the structural response prediction for damage identification are described, including the structural model of the case study, wavelet spectrum, convolutional neural network, input ground motion, and scale factor of records. Next, the application of the methodology is carried out by two processes: training and validation. The results and the comparison of the prediction and reference values of the case study are shown in Section 3. In Section 4, a summary and discussion of the research results are presented.

2. Structural Response Prediction Method

2.1. Structural Model and Structural Responses for Damage Identification

In this study, a lumped mass model (LMM) is considered as the structural model of the building, which takes into account the concentrated mass and the hysteresis model in each story of a low- to mid-rise building as shown in Figure 1.
The structural responses (displacement, acceleration, etc.) of each story of the LMM under the ground motion acceleration are obtained by a time history response analysis using the STERA_3D software [28]. The process is shown in Figure 2.
The maximum ductility ratio (ductility ratio from now on) indicates the amount of inelastic deformation over the yielding threshold as defined in Figure 3. This parameter is related to damage in the structural components of a building. Damage identification based on the ductility ratio is based on the performance-based guideline developed by JSCA [27] as follows: a ductility ratio <1.0 means no damage, a ductility ratio ≥1.0 but <2.0 means minor damage, a ductility ratio ≥2.0 but <3.0 means significant damage, a ductility ratio ≥3.0 but <4.0 means severe damage, and a ductility ratio ≥4.0 means collapse. These values are shown in Table 1. Notice that the ductility ratio is always greater than 1, however, in this study, ratios less than 1 are obtained as well to differentiate between the elastic and inelastic behavior.
The maximum story drift ratio (story drift ratio from now on) represents the maximum relative displacement that a certain story reaches that is associated with the damage of structural components as defined in Figure 4. A larger story drift (relative displacement) after the yielding stage corresponds to a larger extent of damage. Damage identification based on the story drift ratio is based on the performance-based guideline developed by JSCA [27], as follows: a story drift ratio <1/300 means no damage, a story drift ratio ≥1/300 but <1/150 means minor damage, a story drift ratio ≥1/150 but <1/100 means significant damage, a story drift ratio ≥1/100 but <1/75 is severe damage, and a story drift ratio ≥1/75 means collapse. These values are shown in Table 1.
The maximum absolute acceleration (acceleration from now on) indicates the intensity that a certain story is subjected to and is associated with the damage of nonstructural components. Damage identification based on the acceleration is based on the performance-based guideline developed by JSCA [27], as follows: an acceleration <250 gal means no damage, acceleration ≥250 gal but <500 gal means minor damage, acceleration ≥500 gal but <1000 gal means significant damage, an acceleration ≥1000 gal but <1500 gal means severe damage, and an acceleration ≥1500 gal means collapse. These values are shown in Table 1.
Note that the damage condition after severe damage in all cases is considered collapse condition. Besides, no damage and minor damage represent a building that is safe for use, significant damage represents a building that can have restricted use, and severe damage represents a building that is unsafe for use, that is, a value greater than minor damage is a restricted or unsafe condition, which is a parameter used for evacuating the building.

2.2. Wavelet Spectrum

Various transformation functions are used to extract the characteristics of a signal. For example, the Fourier transform can be used to obtain the frequency components of a signal, but it cannot capture the changes over time. On the other hand, if the frequency component varies with time, there are methods such as using the instantaneous frequency or the short-time Fourier transform, both of which have the property that the resolution of time and frequency is constant. However, in actual analysis, it is often the case that low frequency components change slowly over time, while high frequency components change rapidly over time. In the wavelet transform, the optimal time and frequency resolution for each component can be obtained by changing the time resolution according to the frequency of the signal component (see Figure 5 and Figure 6, respectively).
The continuous wavelet transform (CWT) of a signal s(t) is given by Equation (1):
w a ,   b = 1 a s t ψ * t b a d t
where the function ψ t is the mother wavelet (Morlet wavelet [29] used in this study), and “a” and “b” are dilation (scale) and translation (position) parameters, respectively [20]. The symbol (*) denotes complex conjugation. Therefore, the wavelet transform permits transformation from a signal to a spectrum (wavelet spectrum) in two dimensions (time and frequency) with coefficients (scales) that represent the intensity of the signal, in the time-domain and frequency-domain. The wavelet spectrum shows the highest intensity of the wave on the time-domain and frequency-domain only in one graph (Figure 7b). As a reference, Figure 7a shows the acceleration wave, and Figure 7c shows a 3D graph of the wavelet spectrum.
This is a powerful tool for extracting the characteristics of the waveform signals such as response acceleration, velocity, and displacement. Thus, the wavelet spectrum of the acceleration response waveform obtained from the accelerometer installed in the building is computed in this study and used as an input to the CNN model.

2.3. Convolutional Neural Network (CNN)

An image is processed by a computer as a grayscale image (image from now on) represented by an arrangement of numbers. For example, in Figure 8, the right matrix contains numbers between 0 and 255, each of which corresponds to the pixel brightness in the left image [30].
The convolution of the input image is performed by applying a set of weights, also known as a kernel or filter, as shown in Figure 9 [31].
In the CNN method, images are used as input data, and for every input data set, the features of the input data are extracted by the convolution of the kernels. However, this convolution step loses information that might exist on the border of the image because they are only captured when the kernel slides (the kernel has to start and finish its process on the image borders) [32]. For this reason, the size of the input image is reduced as shown in Figure 9 (from input size: 3 × 3 to output size: 2 × 2). In order to obtain the same size as the original input, it is possible to apply the “same padding”, also called “zero-padding”, method (used in this study), which means the input is filled with zeros along its border as shown in Figure 10.
Then, every resultant matrix is evaluated by a nonlinear activation function to allow for the learning of more complex models. The nonlinear activation function (activation function from now on) used in this study is the rectified linear unit (ReLU), defined as the function Y = max(X, 0) [33], as shown in Figure 11.
Finally, the new input data, the feature maps, are obtained. The process from the input data to the feature maps using the previous definitions is called the typical convolutional layer (see Figure 12).
Usually, the typical convolutional layer is followed by a pooling layer to reduce the number of operations since the number of parameters increases as the network processes more kernels. A type of pooling layer is the “maximum pooling” or “max pooling” process, which takes the maximum value sliding along the feature map [30], as shown in Figure 13.
The pooling layer is required for image classification. It adjusts the features’ robustness to noise and disorder by reducing the resolution of the previous feature maps [33]. However, in this study, the CNN models with and without the maximum pooling layer were trained, and the CNN model without the maximum pooling layer converged on the output prediction more effectively. Therefore, the pooling layer is not used in the proposed CNN model.
Usually, a hierarchical architecture is used in advance to propose the number of convolutional layers for the CNN architecture model [34]. In this study, after training the CNN models with different numbers of convolutional layers, 17 convolutional layers are finally used in the proposed CNN model, as shown in Figure 14.
Subsequently, the last convolutional layer is fully connected to the 1D layer or the flattening layer [33] (matrix of one column) with the number of stories as shown in Figure 14. In order to optimize the convergence and measure the error between the predicted and reference output, “Adam” [35] and mean squared error (MSE) are used as the optimizer function and the loss function. Equation (2) defines MSE, where ypred is the prediction output, yref is the reference output, and N is the number of samples.
MSE = 1 N · i = 1 N y p r e d y r e f 2 ,
The flattening matrix contains structural responses for the damage identification, which can be the ductility ratio, story drift ratio, and acceleration (see Section 2.1 for their definitions). Figure 14 shows the CNN scheme used in this research.
Table 2 shows the architecture of the CNN for the structural response prediction method. This was finalized by extensive analysis of trained CNNs in advance. In Table 2, “No. kernels” is the number of filters or kernels assigned in each layer. Ten different kernels are used for the first layer and eight kernels are used for the other layers. Two types of kernel initializer are used in this study. “He_Normal” is used for the first four convolutional layers and “glorot_uniform” is used for the rest of the others. The kernel size is 10 × 10 for the first convolutional layer and 3 × 3 for the rest. The “same padding” and ReLU activation function are used in all convolutional layers.
Figure 15 shows the convergence curve of the CNN model using the CNN architecture shown in Table 2, where “Loss” is the value of the loss function, and “Number of epochs” is the number of training iterations over the input data [31].
Firstly, the CNN model is trained with known input and output data. This is called the “training process”. Subsequently, new unknown input data are used to validate the trained CNN model by comparing the output data (structural responses for damage identification) with the reference structural responses. This process is called the “validation process” and the MSE function is used to evaluate the error.

2.4. Case Study and Input Ground Motion

2.4.1. Case Study

The case study is a building of five stories with the following considerations (see Table 3 for more details):
  • The fundamental period is considered the following: T1 = 0.025 H (H: total height of the building). The height of each story (h) is considered to be 4.0 m, then, H is 20 m and T1 is 0.5 s.
  • LMM is used for the model of the building, and the bilinear hysteresis model (see Figure 5) is used to represent the nonlinear relationship between shear force and story drift for each story.
  • The structural responses for damage identification (ductility ratio, story drift ratio, and acceleration) under earthquake ground motions are calculated by STERA 3D software [28].
In order to build the bilinear hysteresis model, the yielding shear force (Qi) is calculated to be equal to the design shear force under the horizontal seismic load according to Japanese code. Moreover, the story stiffness (ki) is calculated so that the first mode shape becomes a triangular shape. Table 4 shows the parameters used in this study to define the bilinear hysteresis model in each story. The post-yield stiffness ratio (k2/k1, see Figure 16) is 0.1 for each story.

2.4.2. Input Ground Motion

Table 5 shows the 25 earthquake ground motions considered in this study. Every earthquake contains two directions (E–W and N–S). As a consequence, the total number of records used is 50. As mentioned, there are two processes in the CNN method—the training and the validation processes. For this reason, the records are subdivided into two groups. The number of records for the training is 40 (20 earthquakes) and the number of records for the validation is 10 (5 earthquakes). This obeys the split ratio recommended for typical CNN procedures (80% training records and 20% validation records). The earthquakes are selected randomly to avoid extracting the same characteristics between different records. Figure 17 shows the acceleration response spectrum of the 50 records scaled to have the same values at the fundamental period of the structure (T1 = 0.5 s) as Sa (T1) = 100 gal.

2.4.3. Scale Factor of Records

The linear and nonlinear behavior of the structure is obtained by using different intensities of earthquake ground motions. Thus, the records are scaled to include a wide range of earthquake intensity. In order to evaluate the range of the scale factors, an incremental dynamic analysis with the structural responses for damage identification is conducted by taking into account the variation of the Peak Ground Acceleration (PGA), and the ordinate of the response acceleration spectrum evaluated on the fundamental period of the structure (Sa(T1)).
Figure 18 presents the incremental structural responses for damage identification (ductility ratio, story drift ratio, and acceleration) in each story for the input ground motion “El Centro 1940” (Figure 18a) and “Northridge” (Figure 18b) using the same scale factor applied to Sa(T1) such that the minimum scale factor produces Sa(T1) = 100 gals and the maximum scale factor produces Sa(T1) = 1500 gals. Figure 18a shows the structural response under the maximum PGA of El Centro up to 500 gals. As shown in Figure 18b, the maximum PGA of Northridge must be around 1000 gals to achieve the same degree of response. Furthermore, the PGA of the threshold of the nonlinear behavior is around 150 gals in Figure 18a (El Centro) and 250 gals in Figure 18b (Northridge). On the other hand, the relationship between the responses and Sa(T1) is roughly the same in Figure 18a (El Centro 1940) and 18b (Northridge). Therefore, the Sa(T1) is more stable for characterizing the structural response of the structure. For this reason, the scale factor is based on Sa(T1) such that the minimum scale factor produces Sa(T1) = 100 gals, and the maximum scale factor produces Sa(T1) = 1500 gals and Sa(T1) = 1000 gals to train and validate the CNN model, respectively. Figure 19 shows the Acceleration Response Spectra of the “Loma Prieta” input ground motion considering the minimum and maximum scale factor and the original value.

2.5. Machine Learning Methodology

The methodology for predicting the structural responses for damage identification was as follows:
  • The wavelet spectrum was obtained from the time–history acceleration response on the upper floor of the building. The frequency range was from 0.1/T1 to 5/T1, where T1 is the fundamental period of the case study structure (T1 = 0.5s), which is from 0.2 Hz to 10 Hz. This covered the high and low frequencies produced during high mode vibrations and nonlinear frequencies.
  • There were two sets of scale factors for the training and validation of CNN processes.
  • The training scale factor set was the minimum scale factor, which produces Sa(T1) = 100 gal, to the maximum, which produces Sa(T1) = 1500 gal, at increments of 50 gal.
  • The validation scale factor set was the minimum scale factor, which produces Sa(T1) = 100 gal, to the maximum, which produces Sa(T1) = 1000 gal, at increments of 25 gal.
  • There were 1160 structural analyses conducted for the training process by considering 40 records with 29 scale factors, while there were 370 structural analyses conducted for the validation process by considering 10 new records with 37 scale factors. Therefore, 1530 structural analyses carried out were used in this study.
The application of the methodology to predict the structural responses for damage identification was conducted as follows:
TRAINING PROCESS
  • STEP 01: 40 training records are scaled with 29 scale factors per record. As a result, 1160 scaled records are generated.
  • STEP 02: 1160 structural analyses are carried out for the structural model of the case study. As a result, 1160 absolute acceleration data on the upper floor are obtained. Additionally, the responses for damage identification (ductility ratio, story drift, and acceleration) are computed from the structural analyses for validating and calibrating the CNN model.
  • STEP 03: 1160 wavelet spectra are obtained from the absolute acceleration of the previous step. The wavelet spectra are the input data for training the CNN model.
  • STEP 04: The CNN model is trained for each structural response for damage identification (ductility ratio, story drift ratio, and acceleration).
VALIDATION PROCESS
  • STEP 01: 10 validation records are scaled with 37 scale factors per record. As a result, 570 scaled records are generated.
  • STEP 02: 370 structural analyses are carried out for the structural model of the case study. As a result, 370 absolute acceleration data on the upper floor are obtained. Additionally, the responses for damage identification (ductility ratio, story drift, and acceleration) are computed as reference outputs to validate the prediction.
  • STEP 03: 370 wavelet spectra are obtained from the absolute acceleration of the previous step. The wavelet spectra are the input data for predicting the structural response for the damage identification using the trained CNN model.
  • STEP 04: Prediction outputs are calculated using the CNN model for each structural response for damage identification (ductility ratio, story drift ratio, and acceleration).
  • STEP 05: The reference and prediction outputs are compared.

3. Prediction and Validation of the Case Study

An example of the analysis results is shown in Figure 20. Figure 20a shows the ductility ratio results under the scaled Petrolia California E–W records, comparing the prediction (horizontal axis) and the reference (vertical axis). In the figure, the straight line represents the perfect prediction. The points represent the results of each story and scale factor defined in Section 2.1 and Section 2.4. Additionally, Figure 20a shows the regions that define the damage condition. The green, yellow, orange, and red regions represent the no damage, minor damage, significant damage, and severe damage conditions, respectively. The collapse condition is considered for any value greater than the severe damage condition. The dashed red rectangle encloses the region for any value that is greater than the minor damage condition and means that the use of the building is restricted or unsafe (condition for evacuating the building). Figure 20b shows an example of the prediction and reference values of each story for a scale factor that produces Sa(T1) = 900 gal.
Figure 21 shows the results of the ductility ratio, story drift ratio, and acceleration for the validation process under the scaled Petrolia California N–S records. The regions that define the damage condition are also shown in the figure. As seen in Figure 21b, the story drift ratios do not reach the significant damage, severe damage, and collapse condition. Likewise, the restricted or unsafe use condition is not reached. Figure 22 shows the prediction and reference values of the ductility ratio, story drift ratio, and acceleration on each floor considered under the same record for a scale factor that produces Sa(T1) = 875 gal.
The coefficient of correlation (r) is used to measure the accuracy of the CNN model in this study, and it is defined as shown in Equation (3):
r = 1 N · i N y p r e d ,   i y p r e d ¯ y r e f ,   i y r e f ¯ 1 N · i N y p r e d ,   i y p r e d ¯ 2 · 1 N · i N y r e f ,   i y r e f ¯ 2
where ypred is the prediction output by the CNN model, yref is the reference output by the structural analysis, y p r e d ¯ and y r e f ¯ are the mean of ypred and yref, respectively, and N is the number of samples. Table 6 shows the r-values for the validation process. The average values of the r-values of the ductility ratio, story drift ratio, and acceleration are 0.905, 0.846, and 0.829, respectively. In particular, the accuracy of the estimation of the ductility ratio is the highest.
Two new ratios are introduced, the damage condition ratio (DCR) and the restricted or unsafe use ratio (RUUR), to examine the accuracy of the prediction of structural damage. The damage condition ratio (DCR) is defined as the ratio of the number of the predicted values and the number of reference values inside the damage condition region as shown in Equation (4). Likewise, the restricted or unsafe use ratio (RUUR) is defined as the number of the predicted values and the number of reference values inside the restricted or unsafe region as shown in Equation (5).
DCR   = No .   of   Predicted   values   inside   the   damage   condition   region No .   of   Reference   values   inside   the   damage   condition   region · 100 %
RUUR = No .   of   Predicted   values   inside   the   restricted   or   unsafe   use   region No .   of   Reference   values   inside   the   restricted   or   unsafe   use   region · 100 %
Figure 23 shows the comparison of the DCR and RUUR for the ductility ratio. In general, the DCR of no damage and collapse condition are larger and more accurate than others. In most cases, RUUR has high precision—greater than 80%. Notice that DCR and (or) RUUR for some records is not reached because the structural response is not over the limit for being measured.
Figure 24 shows the comparison of the DCR and RUUR for the story drift ratio. In general, the DCR of no damage and minor damage condition are larger and more accurate than other conditions. Few data reach DCR of severe damage and collapse conditions.
Figure 25 shows the comparison of the DCR and RUUR for the acceleration. In general, the DCR of significant damage condition is larger and more accurate than others. Few data reach DCR of severe damage and collapse conditions. In most cases, RUUR has high precision—greater than 90%.

4. Conclusions and Discussion

In this study, a method is proposed with which to estimate the damage of a building by applying a machine learning method from the acceleration response at the upper floor of the building. The results of this research are summarized as follows:
  • The maximum ductility factor, inter-story drift ratio, and maximum response acceleration of each floor were predicted via a CNN model using the acceleration record at the upper floor of the building.
  • The wavelet spectrum of the acceleration record of the upper floor of the building was used as the input of the CNN model to account for the non-stationarity of both the amplitude and frequency of the building response.
  • A CNN model was trained for the linear to nonlinear response of a building by inputting two horizontal components of 20 different earthquake ground motions with varying scales. The trained CNN model was then validated by inputting the two-directional horizontal components of five different earthquake motions to the building with different scales.
  • The correlation coefficients between the predicted values and the reference values by the CNN model exceeded 0.8 for all response values, confirming the high accuracy of the model.
  • The damage information evaluated by the CNN model was classified according to the target performance of the building as “no damage”, “minor damage”, “significant damage”, and “severe damage”. Furthermore, new ratios, DCR and RUUR, are proposed to examine the accuracy of the prediction of structural damage.
Using this method, it is possible to estimate the degree of damage to a building immediately after an earthquake using only the record of accelerometers installed on the upper floor of the building. The results will be useful for countermeasures after an earthquake, such as evacuation and decisions on the continued use of the building.

Author Contributions

Conceptualization, T.S., E.A.M.A. and M.D.B.; methodology, E.A.M.A. and M.D.B.; software, T.S., E.A.M.A. and M.D.B.; formal analysis, E.A.M.A.; investigation, E.A.M.A. and M.D.B.; writing—original draft preparation, E.A.M.A.; writing—review and editing, T.S.; supervision, T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. National Institute of Civil Defense (INDECI) of Peru. Emergency Report No. 53—COEN-SINADECI; INDECI: Pisco, Peru, 2007. [Google Scholar]
  2. United Nations World Conference on Disaster Reduction. In Proceedings of the Report of the World Conference on Disaster Reduction, Kobe, Japan, 18–22 January 2005.
  3. Tilio, L.; Murgante, B.; Di Trani, F.; Vona, M.; Masi, A. Resilient City and Seismic Risk: A Spatial Multicriteria Approach. Comput. Sci. Appl.—ICCSA 2011, 6782, 410–422. [Google Scholar]
  4. Freddi, F.; Galasso, C.; Cremen, G.; Dall’Asta, A.; Di Sarno, L.; Giaralis, A.; Gutiérrez-Urzúa, F.; Málaga-Chuquitaype, C.; Mitoulis, S.A.; Petrone, C.; et al. Innovations in earthquake risk reduction for resilience: Recent advances and challenges. Int. J. Disaster Risk Reduct. 2021, 60, 102267. [Google Scholar] [CrossRef]
  5. Limongelli, M.P. Seismic Structural Health Monitoring: From Theory to Successful Applications; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  6. Gkoktsi, K.; Giaralis, A. A compressive MUSIC spectral approach for identification of closely-spaced structural natural frequencies and post-earthquake damage detection. Probabilistic Eng. Mech. 2020, 60, 103030. [Google Scholar] [CrossRef]
  7. Goulet, J.A.; Michel, C.; Der Kiureghian, A. Data-driven post-earthquake rapid structural safety assessment. Earthq. Eng. Struct. Dyn. 2015, 44, 549–562. [Google Scholar] [CrossRef] [Green Version]
  8. Sivasuriyan, A.; Vijayan, D.S.; Górski, W.; Wodzynski, Ł.; Vaverková, M.D.; Koda, E. Practical Implementation of Structural Health Monitoring in Multi-Story Buildings. Buildings 2021, 11, 263. [Google Scholar] [CrossRef]
  9. Antunes, P.; Lima, H.; Varum, H.; André, P. Optical fiber sensors for static and dynamic health monitoring of civil engineering infrastructures: Abode wall case study. Measurement 2012, 45, 1695–1705. [Google Scholar] [CrossRef]
  10. Roghaei, M.; Zabihollah, A. An Efficient and Reliable Structural Health Monitoring System for Buildings after Earthquake. APCBEE Procedia 2014, 9, 309–316. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, J.F.; Lin, C.-C.; Yen, S.-M. A story damage index of seismically excited buildings based on modal frequency and mode shape. Eng. Struct. 2007, 29, 2143–2157. [Google Scholar] [CrossRef]
  12. Xu, K.; Mita, A. Estimation of maximum drift of multi-degree-of-freedom shear structures with unknown parameters using only one accelerometer. Struct. Control. Health Monit. 2021, 28, e2799. [Google Scholar] [CrossRef]
  13. Worden, K.; Manson, G. The application of machine learning to structural health monitoring. Phil. Trans. R. Soc. A 2007, 365, 515–537. [Google Scholar] [CrossRef]
  14. Oh, B.P.; Park, H. Seismic response prediction method for building structures using convolutional neural network. Struct. Control. Health Monit. 2020, 27, e2519. [Google Scholar] [CrossRef]
  15. Tsuchimoto, K.; Narazaki, Y.; Hoskere, V.; Spencer, B. Rapid postearthquake safety evaluation of buildings using sparse acceleration measurements. Struct. Health Monit. 2021, 20, 1822–1840. [Google Scholar] [CrossRef]
  16. Tsuchimoto, K.; Narazaki, Y.; Spencer, B., Jr. Development and Validation of a Post-Earthquake Safety Assessment System for High-Rise Buildings Using Acceleration Measurements. Mathematics 2021, 9, 1758. [Google Scholar] [CrossRef]
  17. Huerta-Lopez, C.; Shin, Y.; Powers, E.; Roesset, J. Time-Frequency Analysis of Earthquake Records. In Proceedings of the 12th World Conference on Earthquake Engineering, Auckland, New Zealand, 30 January–4 February 2000. [Google Scholar]
  18. Tao, D.; Lin, J.; Lu, Z. Time-Frequency Energy Distribution of Ground Motion and Its Effect on the Dynamic Response of Nonlinear Structures. Sustainability 2019, 11, 702. [Google Scholar] [CrossRef] [Green Version]
  19. Cao, H.; Friswell, M. The effect of energy concentration of earthquake ground motions on the nonlinear response of RC structures. Soil Dyn. Earthq. Eng. 2009, 29, 292–299. [Google Scholar] [CrossRef]
  20. Spanos, P.D.; Giaralis, A.; Politis, N.P.; Roesset, J.M. Numerical Treatment of Seismic Accelerograms and of Inelastic Seismic Structural Responses Using Harmonic Wavelets. Comput.-Aided Civ. Infrastruct. Eng. 2007, 22, 254–264. [Google Scholar] [CrossRef] [Green Version]
  21. Balafas, K.; Kiremidjian, A. Development and validation of a novel earthquake damage estimation scheme based on the continuous wavelet transform of input and output acceleration measurements. Earthq. Eng. Struct. Dyn. 2015, 44, 501–522. [Google Scholar] [CrossRef]
  22. Noh, H.Y.; Nair, K.K.; Lignos, D.G.; Kiremidjian, A.S. Use of Wavelet-Based Damage-Sensitive Features for Structural Damage Diagnosis Using Strong Motion Data. J. Struct. Eng. 2011, 137, 1215–1228. [Google Scholar]
  23. Xu, C.; Guan, J.; Bao, M.; Lu, J.; Ye, W. Pattern recognition based on time-frequency analysis and convolutional neural networks for vibrational events in φ-OTDR. Opt. Eng. 2018, 57, 016103. [Google Scholar] [CrossRef]
  24. Dokht, R.; Kao, H.; Visser, R.; Smith, B. Seismic Event and Phase Detection Using Time–Frequency Representation and Convolutional Neural Networks. Seismol. Res. Lett. 2019, 90, 481–490. [Google Scholar] [CrossRef]
  25. Mousavi, S.; Zhu, W.; Sheng, Y.; Beroza, G. CRED: A Deep Residual Network of Convolutional and Recurrent Units for Earthquake Signal Detection. Sci. Rep. 2019, 9, 1–14. [Google Scholar] [CrossRef]
  26. Liao, W.; Chen, X.; Lu, X.; Huang, Y.; Tian, Y. Deep Transfer Learning and Time-Frequency Characteristics-Based Identification Method for Structural Seismic Response. Front. Built Environ. 2021, 7, 627058. [Google Scholar] [CrossRef]
  27. The guide to safe buildings. In JSCA Performance-Based Seismic Design; JSCA: Tokyo, Japan, 2000. (In Japanese)
  28. Saito, T. Structural Earthquake Response Analysis, STERA_3D Version 10.8. Available online: http://www.rc.ace.tut.ac.jp/saito/software-e.html (accessed on 1 October 2020).
  29. Torrence, C.; Compor, G.P. A Practical Guide to Wavelet Analysis. Bull. Am. Meteorol. Soc. 1998, 79, 61–78. [Google Scholar] [CrossRef] [Green Version]
  30. Mueller, J.P.; Massaron, L. Deep Learning. In For Dummies, 1st ed.; IDG: Needham, MA, USA, 2019. [Google Scholar]
  31. Goodfellow, I.; Bengio, Y.; Courville, A. Deep Learning; MIT Press: Cambridge, MA, USA, 2016. [Google Scholar]
  32. Albawi, S.; Mohammed, A.; Al-Zawi, S. Understanding of a Convolutional Neural Network. In Proceedings of the 2017 International Conference on Engineering and Technology (ICET), Antalya, Turkey, 21–23 August 2017. [Google Scholar]
  33. Namatēvs, I. Deep Convolutional Neural Networks: Structure, Feature Extraction and Training. Inf. Technol. Manag. Sci. 2017, 20, 40–47. [Google Scholar] [CrossRef] [Green Version]
  34. Gerón, A. Hands-On Machine Learning with Scikit-Learn, Keras, and TensorFlow, 2nd ed.; O’Reilly: Newton, MA, USA, 2019. [Google Scholar]
  35. Kingma, D.P.; Ba, J. Adam: A Method for Stochastic Optimization. In Proceedings of the International Conference Learn. Represent. (ICLR), San Diego, CA, USA, 5–8 May 2015. [Google Scholar]
Figure 1. Structural model.
Figure 1. Structural model.
Sensors 21 06795 g001
Figure 2. The general process of structural analysis.
Figure 2. The general process of structural analysis.
Sensors 21 06795 g002
Figure 3. Definition of ductility ratio.
Figure 3. Definition of ductility ratio.
Sensors 21 06795 g003
Figure 4. Definition of story drift ratio.
Figure 4. Definition of story drift ratio.
Sensors 21 06795 g004
Figure 5. Types of dilation of the mother wavelet function.
Figure 5. Types of dilation of the mother wavelet function.
Sensors 21 06795 g005
Figure 6. Translation of wavelets over time.
Figure 6. Translation of wavelets over time.
Sensors 21 06795 g006
Figure 7. (a) Acceleration wave; (b) 2D wavelet spectrum; (c) 3D wavelet spectrum.
Figure 7. (a) Acceleration wave; (b) 2D wavelet spectrum; (c) 3D wavelet spectrum.
Sensors 21 06795 g007
Figure 8. Digital image by an arrangement of pixels represented as numbers.
Figure 8. Digital image by an arrangement of pixels represented as numbers.
Sensors 21 06795 g008
Figure 9. Convolution process of a part of an image by matrix multiplication. The symbol (*) means the convolution operator.
Figure 9. Convolution process of a part of an image by matrix multiplication. The symbol (*) means the convolution operator.
Sensors 21 06795 g009
Figure 10. Padding same. The input image is filled with zeros along its border so that the output size is the same as the original input size. The symbol (*) means the convolution operator.
Figure 10. Padding same. The input image is filled with zeros along its border so that the output size is the same as the original input size. The symbol (*) means the convolution operator.
Sensors 21 06795 g010
Figure 11. ReLU activation function.
Figure 11. ReLU activation function.
Sensors 21 06795 g011
Figure 12. Typical convolutional layer, note that the input data is a representation of the wavelet coefficient’s matrix.
Figure 12. Typical convolutional layer, note that the input data is a representation of the wavelet coefficient’s matrix.
Sensors 21 06795 g012
Figure 13. Maximum pooling process.
Figure 13. Maximum pooling process.
Sensors 21 06795 g013
Figure 14. Convolutional neural network scheme. Structural response: ductility ratio, story drift ratio, or acceleration.
Figure 14. Convolutional neural network scheme. Structural response: ductility ratio, story drift ratio, or acceleration.
Sensors 21 06795 g014
Figure 15. Converge curve of the trained CNN.
Figure 15. Converge curve of the trained CNN.
Sensors 21 06795 g015
Figure 16. Bilinear hysteresis model for each story.
Figure 16. Bilinear hysteresis model for each story.
Sensors 21 06795 g016
Figure 17. Acceleration response spectrum of 50 records scaled to have the same values at the fundamental period T1 = 0.5 s.
Figure 17. Acceleration response spectrum of 50 records scaled to have the same values at the fundamental period T1 = 0.5 s.
Sensors 21 06795 g017
Figure 18. (a) Incremental structural responses for damage identification in each story with “El Centro 1940” input ground motion; (b) Incremental structural responses for damage identification in each story with “Northridge” input ground motion.
Figure 18. (a) Incremental structural responses for damage identification in each story with “El Centro 1940” input ground motion; (b) Incremental structural responses for damage identification in each story with “Northridge” input ground motion.
Sensors 21 06795 g018aSensors 21 06795 g018b
Figure 19. Acceleration response spectra of the “Loma Prieta” input ground motion. The red line is with the maximum scale factor such that it produces Sa(T1) = 1500 gal, the green line is with the minimum scale factor such that it produces Sa(T1) = 100 gal, and the black dashed line considers the original input ground motion.
Figure 19. Acceleration response spectra of the “Loma Prieta” input ground motion. The red line is with the maximum scale factor such that it produces Sa(T1) = 1500 gal, the green line is with the minimum scale factor such that it produces Sa(T1) = 100 gal, and the black dashed line considers the original input ground motion.
Sensors 21 06795 g019
Figure 20. Example of the ductility ratio results (Petrolia California E–W record); (a) Comparison between prediction and reference values (points) and damage condition regions; (b) Prediction and reference ductility ratio of each story for a scale factor that produces Sa(T1) = 900 gal.
Figure 20. Example of the ductility ratio results (Petrolia California E–W record); (a) Comparison between prediction and reference values (points) and damage condition regions; (b) Prediction and reference ductility ratio of each story for a scale factor that produces Sa(T1) = 900 gal.
Sensors 21 06795 g020
Figure 21. Comparison between reference and prediction of the Petrolia California N–S record for the validation process of (a) Ductility ratio; (b) Story drift ratio, and (c) Acceleration.
Figure 21. Comparison between reference and prediction of the Petrolia California N–S record for the validation process of (a) Ductility ratio; (b) Story drift ratio, and (c) Acceleration.
Sensors 21 06795 g021
Figure 22. Prediction and reference values on each floor of the Petrolia California N–S record and scale factor that produces Sa(T1) = 875 gal for the validation process of (a) Ductility ratio; (b) Story drift ratio, and (c) Acceleration.
Figure 22. Prediction and reference values on each floor of the Petrolia California N–S record and scale factor that produces Sa(T1) = 875 gal for the validation process of (a) Ductility ratio; (b) Story drift ratio, and (c) Acceleration.
Sensors 21 06795 g022
Figure 23. DCR and RUUR for ductility ratio of the validation process.
Figure 23. DCR and RUUR for ductility ratio of the validation process.
Sensors 21 06795 g023
Figure 24. DCR and RUUR for story drift ratio of the validation process.
Figure 24. DCR and RUUR for story drift ratio of the validation process.
Sensors 21 06795 g024
Figure 25. DCR and RUUR for the acceleration of the validation process.
Figure 25. DCR and RUUR for the acceleration of the validation process.
Sensors 21 06795 g025
Table 1. Proposal of damage condition according to the structural response for damage identification.
Table 1. Proposal of damage condition according to the structural response for damage identification.
Damage Condition No DamageMinor DamageSignificant DamageSevere DamageCollapse
Ductility ratio<1≥1.0 but <2.0≥2.0 but <3.0≥3.0 but <4.0≥4.0
Story drift ratio<1/300≥1/300 but <1/150≥1/150 but <1/100≥1/100 but < 1/75≥1/75
Acceleration (gal)<250≥250 but <500≥500 but <1000≥1000 but < 1500≥1500
Table 2. CNN architecture for the structural response prediction method.
Table 2. CNN architecture for the structural response prediction method.
No. of LayerNo. KernelsKernel SizePaddingKernel InitializerActivation Function
Convolutional Layer 011010 × 10SameHe NormalReLU
Convolutional Layer 0283 × 3SameHe NormalReLU
Convolutional Layer 0383 × 3SameHe NormalReLU
Convolutional Layer 0483 × 3SameHe NormalReLU
Convolutional Layer 0583 × 3Sameglorot_uniformReLU
Convolutional Layer 0683 × 3Sameglorot_uniformReLU
Convolutional Layer 0783 × 3Sameglorot_uniformReLU
Convolutional Layer 0883 × 3Sameglorot_uniformReLU
Convolutional Layer 0983 × 3Sameglorot_uniformReLU
Convolutional Layer 1083 × 3Sameglorot_uniformReLU
Convolutional Layer 1183 × 3Sameglorot_uniformReLU
Convolutional Layer 1283 × 3Sameglorot_uniformReLU
Convolutional Layer 1383 × 3Sameglorot_uniformReLU
Convolutional Layer 1483 × 3Sameglorot_uniformReLU
Convolutional Layer 1583 × 3Sameglorot_uniformReLU
Convolutional Layer 1683 × 3Sameglorot_uniformReLU
Convolutional Layer 1783 × 3Sameglorot_uniformReLU
Table 3. Structural configuration of the case study.
Table 3. Structural configuration of the case study.
DescriptionNomenclature (Units)Value
Number of storiesN5
Story heighth (m)4
Building heightH (m) = h × N20
WidthB (m)30
Area of floorA = B2 (m2)900
Weight per floor areaw (kN/m2)12
Weight of floorW (kN)1080
Fundamental PeriodT1 (s)0.5
Table 4. Parameters of the bilinear hysteresis model used in the case study.
Table 4. Parameters of the bilinear hysteresis model used in the case study.
Storyki (kN/mm)Qi (kN)
587587.87
4157954.15
32091240.66
22441460.87
12611620.00
Table 5. Earthquake ground motions.
Table 5. Earthquake ground motions.
No. Name LocationStation ReferenceMagnitudeDate
Training Process *
01Anza_01USA33.706N, 116.235W/Ground Floor: South WingMw = 5.212/06/2005
02Anza_02USA33.706N, 116.235W/Roof: Center Hallway of S. WingMw = 5.212/06/2005
03El Centro 1940USAImperial Valley EarthquakeMw = 6.918/05/1940
04Kobe 1995JapanGreat Hanshin Earthquake/Kobe Marine ObservatoryMw = 6.917/01/1995
05Loma Prieta_01USA36.974N, 121.952W/Capitola—Fire StationMs = 7.118/10/1989
06Loma Prieta_02USA36.973N, 121.572W/Gilroy #1—Gavilan CollegeMs = 7.118/10/1989
07Loma Prieta_03USA36.987N, 121.536W/Gilroy #3—Gilroy Sewage PlantMs = 7.118/10/1989
08Loma Prieta_04USA37.046N, 121.803W/Corralitos—Eureka Canyon Rd.Ms = 7.118/10/1989
09Loma Prieta_05USA37.118N, 121.550W/Coyote Lake DamMs = 7.118/10/1989
10Loma Prieta_06USA37.255N, 122.031W/Saratoga—Aloha Ave.Ms = 7.118/10/1989
11Northridge_01USA34.068N, 118.439W/Los Angeles—UCLA GroundsMw = 6.717/01/1994
12Northridge_02USA34.236N, 118.439W/Arleta—NordHoff Ave. Fire Station Mw = 6.717/01/1994
13Northridge_03USA34.387N, 118.530W/Newhall—LA County Fire StationMw = 6.717/01/1994
14Petrolia_01USA40.325N, 124.287W/PetroliaMw = 7.025/04/1992
15Petrolia_02USA40.503N, 124.100W/Rio Dell—101/Painter St. OverpassMw = 7.025/04/1992
16Petrolia Aftershock_01 USA40.325N, 124.287W/Petrolia/04/26/92, 07:41:40 UTCMs = 6.626/04/1992
17Petrolia Aftershock_02 USA40.325N, 124.287W/Petrolia/04/26/92, 11:18:25 UTCMs = 6.626/04/1992
18Petrolia Aftershock_03 USA40.026N, 124.069W/Shelter Cove—AirportMs = 6.626/04/1992
19Whittier_01USA34.037N, 118.178W/Los Angeles—Obregon ParkMl = 6.101/10/1987
20Whittier_02USA34.160N, 118.534W/Tarzana—Cedar Hill NurseryMl = 6.101/10/1987
Validation Process *
21Palm SpringsUSA33.962N, 116.509W/Desert Hot SpringsMl = 6.108/07/1986
22Petrolia CaliforniaUSA40.325N, 124.287W/PetroliaMl = 5.917/08/1991
23Taft 1952USAKern County, California EarthquakeMw = 7.321/07/1952
24Tohoku 1978JapanMiyagi Earthquake/Recorded at Tohoku UniversityMs = 7.712/06/1978
25WestmorlandUSA33.037N, 115.623W/WestmorlandMl = 6.026/04/1981
* See Section 2.5 for details.
Table 6. Coefficient of correlation for the validation process.
Table 6. Coefficient of correlation for the validation process.
No.RecordDuctility RatioStory Drift RatioAcceleration
01Palm Springs E–W0.9530.9470.928
02Palm Springs N–S0.8950.9170.951
03Petrolia California E–W0.9330.8450.873
04Petrolia California N–S0.9720.9260.956
05Taft 1952 E–W0.8720.7710.417
06Taft 1952 N–S0.8060.8480.870
07Tohoku 1978 E–W0.8330.4660.562
08Tohoku 1978 N–S0.9430.8900.797
09Westmorland E–W0.9250.9690.985
10Westmorland N–S0.9160.8830.950
Average0.9050.8460.829
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Moscoso Alcantara, E.A.; Bong, M.D.; Saito, T. Structural Response Prediction for Damage Identification Using Wavelet Spectra in Convolutional Neural Network. Sensors 2021, 21, 6795. https://doi.org/10.3390/s21206795

AMA Style

Moscoso Alcantara EA, Bong MD, Saito T. Structural Response Prediction for Damage Identification Using Wavelet Spectra in Convolutional Neural Network. Sensors. 2021; 21(20):6795. https://doi.org/10.3390/s21206795

Chicago/Turabian Style

Moscoso Alcantara, Edisson Alberto, Michelle Diana Bong, and Taiki Saito. 2021. "Structural Response Prediction for Damage Identification Using Wavelet Spectra in Convolutional Neural Network" Sensors 21, no. 20: 6795. https://doi.org/10.3390/s21206795

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