Research Article  Open Access
Ehsan Darvishan, Gholamreza Ghodrati Amiri, Pedram Ghaderi, "A Study on Instantaneous TimeFrequency Methods for Damage Detection of Nonlinear MomentResisting Frames", Shock and Vibration, vol. 2014, Article ID 523675, 14 pages, 2014. https://doi.org/10.1155/2014/523675
A Study on Instantaneous TimeFrequency Methods for Damage Detection of Nonlinear MomentResisting Frames
Abstract
Most of the civil structures exhibit nonlinear hysteresis behavior during earthquakes. However, detection of damage in these structures is a challenging issue due to successive change in structural characteristics during a seismic excitation. The current paper presents a promising approach for damage detection of nonlinear moment frames. First, several instantaneous timefrequency methods including HilbertHuang transform, direct quadrature, Teager energy operator, and higherorder energy operator are investigated as signal processing tools and the most appropriate method is selected using an outlier analysis. Next, a procedure is proposed based on timefrequency analysis in conjunction with clustering to find damage extension in moment frames under a seismic excitation using frequency, amplitude, and energy damage measures. A probabilistic approach is implemented to investigate capability of the procedure for different ground motion records using incremental dynamic analysis. Results show that frequency is not an appropriate feature to detect damage in nonlinear structures.
1. Introduction
During an earthquake many civil structures undergo damages. Therefore, it is important to ensure serviceability and safety of the structures after a seismic excitation. Today, nondestructive evaluation techniques (NDE) are widely utilized for health monitoring and damage detection of structures. At present, current practical methods are based on visual inspection, CT scanning, ultrasonic, stress waves, acoustic emission, and so forth [1]. However, these local damage detection methods are effective only for small structures or structural members. Recently, global vibrationbased techniques have attracted researchers since they can solve the problem of large and complicated structures. These methods include power spectrum, Fourier transform, spectrum analysis, and cepstrum analysis [2]. The main drawback of these procedures is that they cannot trace timevarying nature of signals due to their fixed functions. Instead, timefrequency methods are effective tools for analyzing nonstationary signals which makes them appropriate for health monitoring and fault detection of civil, mechanical, and aerospace systems [3].
Doebling et al. [4] and also Fan and Qiao [5] gave a detailed review of vibrationbased damage detection methods. An experimental investigation was conducted on a threestory shear building by Xu and Chen [6]. They studied applicability of empirical mode decomposition (EMD) to detect sudden change in structural stiffness. Yang et al. [7] tried to extract damage spikes from vibration response of the ASCE benchmark building using EMD. Rezaei and Taheri [8] proposed a new damage index for damage detection of beams using EMD. Chen [9] experimentally investigated the capability of EMD for structural health monitoring of large civil structures. HilbertHuang transform (HHT) was utilized as a signal processing tool to detect existing damages in an actual benchmark problem by Liu et al. [10]. Lin and Chu [11] characterized the acoustic emission signals due to cracks in welded parts of an offshore structure by application of HHT. Hsu et al. [12] proposed a new damage index to detect damage in a steel structure at early state. Xu et al. [13] employed HHT to identify modal properties of Di Wang tower in China. Bouchikhi et al. [14] compared EMD and Bspline EMD associated with HHT and THT for multicomponent nonstationary signals.
However, despite extensive investigations on structural damage detection, most researchers have focused on damage detection of elastic structures by monitoring stiffness change due to damage and researches on nonlinear hysteretic structures are scarce in the literature. Saadat et al. [15] proposed an intelligent parameter varying procedure to identify damage in a nonlinear massspring system. Xia [16] proposed a new approach to detect timevarying characteristics of a nonlinear structure with hysteresis behavior.
This study focuses on detection of damage in nonlinear moment frames due to seismic excitations. This paper is organized into two parts. In the first part, capability of several wellknown timefrequency procedures to extract signals is investigated using acceleration response of the 3story SAC steel structure [17] under a seismic excitation. An outlier analysis is performed to investigate the adverse effect of outliers on statistical parameters of extracted frequencies. The second part presents a promising approach to detect damage due to hysteretic behavior of moment connections of a structure under seismic excitations. Three damage measures including frequency, amplitude, and energy are compared. A probabilistic approach is implemented to quantify damage under various ground motion records.
2. Signal Processing Procedures
2.1. Empirical Mode Decomposition
In 1996, Huang et al. [18] proposed empirical mode decomposition (EMD) for nonlinear and nonstationary data. EMD is an adaptive procedure which relies on two basic assumptions: the number of extrema and zerocrossings must be equal or differ at most by one for all the data; mean value of envelope defined by local maxima and minima must be zero at each point. Extraction of intrinsic mode functions (IMFs) is performed via sifting process described below.
First, determine local extrema of the data. Plot an upper envelope by using a cubic spline passing all the extrema, . Repeat the same procedure for minima to find . Calculate their mean, , as
The difference between signal and is designated as the first component, :
Now, the first component is considered as the signal and the above procedure is repeated times to satisfy the basic conditions:
Thus, will be the first IMF:
A stopping criterion is computed at each repeat as follows:
When the stopping criterion is satisfied the first IMF, , is obtained. To obtain other IMFs the residue, , is considered as the signal and the procedure is repeated:
When the residue is constant or a monotonic signal or no more sifting is possible the procedure is stopped. Finally, we can reconstruct the original data by summing up
2.2. HilbertHuang Transform
HilbertHuang transform (HHT) implements the benefits of EMD and Hilbert spectral analysis (HSA) to extract a signal to its intrinsic mode functions [19]. Hilbert transform of a signal, , is an integral transform that can be thought of as the convolution with the Cauchy kernel where PV is the Cauchy principal value. Frequency content of a signal, , can be obtained by providing the conjugate part to give an analytic complex signal as
Thus where a(t) is the instantaneous amplitude and is the instantaneous phase. The instantaneous frequency is the time derivative of phase function
2.3. Teager Energy Operator
Teager energy operator (TEO) is a powerful technique for timefrequency analysis of nonstationary signals. TEO tracks the modulation energy and determines the instantaneous frequency (FM) and instantaneous amplitude (AM) of a signal [20].
For a continuous TEO, is defined as where and are the first and second time derivatives of , respectively. In the discrete case, the time derivatives are substituted by time differences. Thus, the discrete time TEO is given by
According to (13), TEO only requires three samples to calculate energy. Therefore, it provides high resolution and is able to capture energy variations.
Maragos et al. [20] proposed discrete energy separation algorithm (DESA1) for the instantaneous frequency and amplitude of an AMFM signal as follows:
Or
The main drawback of TEO is its sensitivity. If the value of data points changes abruptly the discrete time differentiation yields unrealistic results. Thus, TEO is appropriate for bandwidth signals.
2.4. HigherOrder Energy Operator
Higherorder energy operator (HEO) is an extension of TEO to a large class of operators [21]. These generalized operators can improve the sampling rate precision and computation effort and show finer frequency resolution in the presence of noise. HEO, , implements four partial th, th, th, and th derivatives of a signal as follows:
in which conditions of and must be satisfied. The combination of , , , yields TEO.
In discrete case, continuous time derivatives are obtained by substituting two sample symmetric differences
Thus, the symmetric higherorder discrete operator is obtained byeven: odd:
By including higherorder derivatives, more sampling points are included which can lead to less sensitivity of energy operator against frequency fluctuations but less instantaneous results:
2.5. Direct Quadrature Method
Direct quadrature (DQ) method is defined based on a normalization scheme. An iterative normalization is applied on the original signal to separate AM and FM parts [22]. First, local maxima of absolute values of data, , are determined. Next, a cubic spline is passed from the local maxima points to create envelope of the signal, . The first normalized data is obtained by
The normalized data is considered as the original data and the same procedure is used repeatedly to obtain the th normalized data as
The normalization process stops when all the values of are equal to or less than unity. is now considered as the empirical FM part of the data
The empirical AM part is defined by
Also, the original signal can be reconstructed by
3. Modeling and Analysis
The SAC 3story building in Los Angeles was selected for analysis. This building is designed as a standard office building as a part of SAC task 5.4.3 [17]. Figure 1(a) depicts the building in plan and elevation. Perimeter moment frames are indicated in Figure 1(b) with thick lines. Design of the structure is performed based on UBC94 [23] and provisions of FEMA 267 [24] are also considered. Grade A36 steel is utilized in conjunction with Gr. 50 steel for beams and girders. Seismic mass of 65.5 kipssec^{2}/ft is assumed for stories and 70.9 kipssec^{2}/ft for the roof. It is assumed that the building is located on stiff soil [23].
(a)
(b)
The open source finite element program, Opensees [25], is used for modeling. The building is modeled as a centerline 3D moment frame. Rigid diaphragm is defined for each story. No cyclic degradation is defined. “Steel01” material is utilized with “fiber” sections to model momentrotation relationship of moment connections. 5% strain hardening is defined. “Displacementbased beamcolumn” element is used to model nonlinear behavior of elements.  effect is considered to account for secondorder displacements during the analysis.
To improve the robustness of the procedure several sensors are installed on the structure. Four sensors are mounted on the corners of each story for each direction. Figure 2 illustrates the orientation of principal directions (i.e., and ) of the building and position of the sensors.
Two horizontal components of “EL Centro Array number 1” record of “Imperial Valley” earthquake are utilized for nonlinear dynamic analyses. Both components are simultaneously imposed on the structure in two orthogonal directions for analysis.
4. Selection of Signal Processing Method
Timefrequency methods are robust procedures to detect instantaneous changes in frequency of a system. However, time dependence of extracted frequencies can increase complexity in identification of the structure characteristic. To clarify this issue, the frequency of vibration of the structure is investigated in healthy state. Four instantaneous timefrequency methods including two wellknown HHT and TEO procedures and two newer DQ and HEO methods are implemented and the results are compared.
For this reason, a nonlinear dynamic analysis is performed. A low scale factor is selected for the acceleration record to ensure elastic behavior of the structure (maximum drift ratio below 0.01). Acceleration output of sensor 17 is recorded during the nonlinear dynamic analysis (Figure 3). Sampling frequency of 800 Hz is assigned for the sensor. EMD is used to extract the acceleration records of the sensor to its IMFs, as depicted in Figure 4. Only the first IMF is used for further process since it contains most of the frequency content. Finally, the aforementioned timefrequency methods are applied to extract instantaneous frequencies and amplitudes (Figure 5). Statistical variables of the frequencies are listed in Table 1. All the calculations are performed using MATLAB software package [26].

(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(a)
(b)
(c)
(d)
As can be easily seen from Figure 5, large fluctuation in frequencies is observed which is associated with numerous spikes. An outlier analysis is carried out to further investigate the outlier frequencies.
5. Outlier Analysis
An outlier is an unusual point that is markedly inconsistent from rest of the data. Outliers can appear accidentally in a set of data, but in most cases this deviated data is due to measurement error or natural deviations in the dataset [27]. A common practice for the former case is to discard outliers. However, statistical tools should be cautiously implemented for the latter case. Outliers can influence the estimation of statistical parameters of the observation. Therefore, it is important to distinguish the origin of outliers.
5.1. Univariate Analysis
This section describes the results of a univariate analysis of frequencies. The following statistical features were considered to compare the procedures including root mean square (RMS), range. and : where is data sample. Table 1 compares the aforementioned extracted frequencies.
Comparison of mean and median values shows the deviation of mean values due to presence of outliers that are significantly deviant from the rest of measurements. Therefore, median values are more appropriate for comparison.
HHT yields the largest median frequency equal to 2.84. Standard deviation of the HHT is large due to high values of outliers. Although the range of frequencies for HHT is less than TEO and DQ, the number of outliers with large values is high which leads to large RMS and kfactor.
In contrast to HHT, HEO has low median frequency and the lowest range, RMS, and Kfactor. The extremely low RMS and kfactors show that although the number of outliers is not less than the others, the outliers with very large values are less. It is notable that although HEO shows poor accuracy to predict the median frequency in comparison to TEO, the low standard deviation of the results shows that using more data samples for numerical differentiation (17) can reduce the fluctuations of the results. However, this can adversely affect the instantaneous nature of the HEO and leads to poor estimations. DQ has the lowest median frequency among all. Also, DQ has the highest range, RMS, and factor. The combination of numerous outliers with considerable deviation from the rest of data set has resulted in large difference between median and mean values and also large standard deviation. This can be due to normalization process (22) which may disassemble the frequency content.
5.2. Multivariate Analysis
To study the effect of outliers for all the sensors simultaneously, a multivariate outlier analysis is performed using a discordancy test. Mahalanobis squared distance (MSD) is designated as a discordancy test variant. MSD for a pdimensional problem with n samples is defined as where is the mean vector of data samples, is the observed outlier, and is the covariance matrix. and are pdimensional vectors, and is a square matrix. Figure 6 plots MSD values for the assumed signal processing procedures.
(a)
(b)
(c)
(d)
It is clear that HHT has the lowest values of MSD. Also the number of MSDs with large values is the least among all. The most outspread results belong to TEO. HEO shows a slight improvement in comparison to TEO. However, the number of MSDs with large values is still remarkable. Although DQ has lower values of MSD compared to TEO and HEO, large discordance is observed specially for the first 5000 samples (first 6 second of seismic analysis). According to the above results, HHT is selected for further postprocessing since it provides more appropriate results.
Since the outliers are observed in the undamaged case of the structure, they cannot be considered as damage spikes as a result of sudden change in structural stiffness. Therefore, we assumed that these outliers cannot be useful to detect damage and eliminated them since they adversely affect the postprocessing results.
5.3. Sampling Rate Effect
A sensitivity analysis is carried out to study the sampling rate effect for HHT procedure. The previous statistical variables are used to measure dispersion of data. Table 2 refers to statistical characteristics of extracted frequency for sampling rate of 25 to 800 Hz.

As shown in Table 2, the median frequencies are almost constant. But standard deviation increases by increase in sampling rate which affects the mean frequency. Increase in dispersion due to sampling rate increase is also confirmed by RMS, range, and kfactor which is due to presence of spikes (outliers) in frequencies. These spikes can appear as a result of numerical differentiation of phase function which is sensitive to sampling time step (11). However, difference in results of 400 Hz and 800 Hz is not sensible. Generally speaking, frequency results are less sensitive to sampling rates over 400 Hz.
On the other hand, the lower values of spikes for lower sampling frequencies (i.e., 25 and 50 Hz) are due to lower Nyquist frequency that limits the maximum frequency representable in the data sample, not as a result of more precise extraction. Thus, 800 Hz sampling frequency was designated for postprocessing. However, the outliers were eliminated as mentioned before.
6. Damage Detection Methodology
This section investigates the capability of HHT to detect damage in a nonlinear moment frame. First, the moment connections of the structures are modeled to present nonlinear behavior. Second, they are modeled to behave elastically. An IDAbased [28] approach is utilized for analysis. Seven intensity levels are defined for analysis covering the range of elastic (maximum drift ratio of 0.01) to near collapse of the structure (maximum drift ratio of 0.1) in accordance with FEMA350 [29]. Table 3 lists the details of the intensity levels and the corresponding response of the nonlinear structure.

It is worth noting that no damage occurs in intensity level 1. Response acceleration of the structure under the aforementioned intensity levels is recorded by all the sensors for both elastic and nonlinear cases.
Three features are implemented to detect damage including frequency, amplitude, and energy. Frequency and amplitude are directly extracted from the HHT. Energy index is defined as sum of square of the first IMF [30]: where t_{0} is the time length of the signal.
We did not filter the frequencies since filtering may disassemble the frequency content of the response. But the outlier frequencies are discarded, as aforementioned. According to Figure 7, Fourier spectrum of response of the healthy state shows that vibration frequency of the structure is mainly between 1 and 7 Hz. Thus, the frequencies outside this range are cut off.
To quantify changes in the frequencies and amplitudes, fuzzy cmeans clustering [31] is employed. Fuzzy cmeans allows each of the data to belong to several clusters by using rationale of fuzzy sets. The aim is to minimize the leastsquared error of data, x_{i}, and center of cluster, c_{j}: where is weighting exponent and is degree of membership of for cluster . Partitioning is an iterative process. In each iteration, the memberships and center of clusters are updated as follows:
When the iteration threshold is reached, the iterations are terminated.
Before clustering data, a silhouette statistic [32] is employed to find the appropriate number of clusters. For a given sampling point i, is the average dissimilarity of to all other data points for any cluster . If is defined by the minimum of all , the silhouette width for the ith point is
By finding the average of silhouette widths for all sampling points the average silhouette width is obtained by
Larger silhouette width indicates better clustered data. For this case using two clusters yields the best results. For instance, the average silhouette width of frequencies for different numbers of clusters is plotted in Figure 8.
In this study, the value of center of the first cluster for frequencies and amplitudes is designated as frequency and amplitude feature indexes, respectively. Also, the energy value is considered as energy feature index. Figure 9 shows the major components of the procedure.
To detect damage, the sum of residues between nonlinear and elastic feature indexes for all the sensors is considered as damage index. Table 4 lists the damage indexes for different intensity levels. As shown in the table, damage indexes of all features are zero for intensity level 1 implying the undamaged state for the nonlinear structure. Frequency index rapidly increases for intensity levels 2 and 3 since the structure does not tolerate severe nonlinearity. The frequency index remains approximately constant for higher intensity levels. It can be due to hysteresis behavior of moment connections which results in successive changes in stiffness of members and consequently successive changes in vibration frequency of the structure. When large numbers of beamcolumn connections undergo successive frequency changes under a moderate or severe seismic action (e.g., intensity levels 4 to 7) the overall frequency change of the structure will be complicated. In this study no drastic change in frequency index was observed for highintensity levels. Thus, frequency may not be an appropriate damage index to estimate extent of damage. Amplitude index shows more uniform behavior. No residue is observed for intensity level 1. However, intensity level 2 shows a nonnegligible residue since the structure behaves slightly nonlinearly. The rate of increase in amplitude index increases for higher levels of damage. Energy index shows the most rational results. The energy index increases more rapidly after 0.01 maximum drift.

Figure 10(a) shows the behavior of normalized values of the indexes versus maximum drift ratios for and directions. It can be easily seen that amplitude and energy indexes show better correlation with maximum drift ratios than frequencies.
(a)
(b)
Also, feature indexes are plotted versus spectral acceleration (S_{a}) as an intensity measure in Figure 10(b) to study the ability of feature indexes to be implemented as a damage measure (DM) [28] in IDA analysis. Acceptable slope change of the amplitude and energy indexes is observed for both directions according to the figure. In contrast, frequency curve shows a negative slope at higher intensity levels. This ability is investigated in more detail in the next section.
7. Ground Motion Variability
The procedure in the previous section was tested for only one ground motion record. However, large variability is often observed in characteristics of different ground motions. In attempt to investigate the damage of the structure under different seismic actions the above procedure was assessed for a set of ground motion records. FEMAP695 [33] farfield record set was selected for analysis. This record set includes twentytwo pairs of horizontal components of ground motions that is claimed to be adequate to consider recordtorecord variability. An IDA analysis is carried out. The analysis is stopped when 10% maximum story drift ratio is reached for each ground motion record. Figure 11 refers to IDA curves of the features in direction. The measures are normalized to provide an ability to compare the features.
(a)
(b)
(c)
Slope of frequency IDA curves significantly increases for high level of spectral acceleration implying the lower rate of frequency change for these intensity levels. Nevertheless, IDA curves show relatively erratic behavior in some cases. Abrupt slope changes occur even in low damage measures. Regarding the figure, amplitudes present more smooth behavior. Slope of amplitude IDA curves gradually decreases to reach the maximum damage measure. However, no sensible slope change is observed after amplitude index of about 0.2, on average, making estimation of extent of damage difficult. Slope of IDA curves for energy index shows a smooth transition from elastic to collapse state. In contrast to amplitude feature, one can observe a moderate slope reduction for energy feature over 0.2. These results are consistent with the previous results.
As expected, large recordtorecord variability exists in IDA curves of all the feathers. For interpretation of randomness included in the records, a probabilistic approach is employed. IDA curves are summarized to form 16%, 50%, and 84% crosssectional fractiles [28], as depicted in Figure 12.
(a)
(b)
(c)
Although the IDA curves of frequency behave irregularly for some records, the median curve shows more uniform behavior. However, it is notable that the median curve cannot show the frequency tendency properly for normalized values of more than about 0.9 since the frequency indexes stay constant or decrease. Still, energy seems to be a more reliable feature to evaluate extent of damage.
To quantify the damage extent, 5 limit states are defined. Table 5 compares median values of DM in and directions. All three median values follow almost similar paths in both directions. This is more obvious for amplitude and energy index. Frequency median values in direction show more rapid increase in S_{a} from 0.5 to 1.0 which may be a result of dispersion of IDA curves in this range since the other features do not show this deviation in direction.

To evaluate sufficiency of the damage measures, dispersion of estimations for different limit states is investigated. It is clear that more dispersion in the results reduces the confidence in estimation of the features. Dispersion of the damage measures is estimated using the conditional factional standard deviation defined by half the distance between 16% and 84% values for a certain damage measure in logarithmic scale [34]. Table 6 shows that frequency presents the largest dispersion. Amplitude and energy indexes show less dispersion which is almost constant for different limit states.

8. Conclusion
This study aims to find a promising method to detect damage in nonlinear moment frames. The capability of several timefrequency procedures is evaluated to find a suitable signal processing method. Then a probabilistic approach is employed to find the extension of damage under various seismic actions. Based on the results, the following major conclusions can be drawn.(i)Despite the fact that several timefrequency procedures are proposed recently, HHT is still a robust procedure to extract frequencies of a signal.(ii)Large dispersion was observed in estimation of damage features due to large recordtorecord variability. Thus, a probabilistic approach can be useful to quantify damage.(iii)Energy variable yields more appropriate results than frequency and amplitude for damage detection of nonlinear moment frames.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The authors gratefully appreciate Samaneh Ghorbani from University of Mohaghegh Ardabili for her constructive comments.
References
 P. Cawley, “Nondestructive testing—current capabilities and future directions,” Journal of Materials Design and Applications, Proceedings of the Institution of Mechanical Engineers L, vol. 215, no. 4, pp. 213–223, 2001. View at: Google Scholar
 C. R. Farrar, S. W. Doebling, and D. A. Nix, “Vibrationbased structural damage identification,” Philosophical Transactions of the Royal Society A, vol. 359, no. 1778, pp. 131–149, 2001. View at: Publisher Site  Google Scholar
 P. C. Chang, A. Flatau, and S. C. Liu, “Review paper: health monitoring of civil infrastructure,” Structural Health Monitoring, vol. 2, no. 3, pp. 257–267, 2003. View at: Publisher Site  Google Scholar
 S. W. Doebling, C. R. Farrar, and M. B. Prime, “A summary review of vibrationbased damage identification methods,” Shock and Vibration Digest, vol. 30, no. 2, pp. 91–105, 1998. View at: Google Scholar
 W. Fan and P. Qiao, “Vibrationbased damage identification methods: a review and comparative study,” Structural Health Monitoring, vol. 10, no. 1, pp. 83–111, 2011. View at: Publisher Site  Google Scholar
 Y. L. Xu and J. Chen, “Structural damage detection using empirical mode decomposition: experimental investigation,” Journal of Engineering Mechanics, vol. 130, no. 11, pp. 1279–1288, 2004. View at: Publisher Site  Google Scholar
 J. N. Yang, Y. Lei, S. Lin, and N. Huang, “HilbertHuang based approach for structural damage detection,” Journal of Engineering Mechanics, vol. 130, no. 1, pp. 85–95, 2004. View at: Publisher Site  Google Scholar
 D. Rezaei and F. Taheri, “Damage identification in beams using empirical mode decomposition,” Structural Health Monitoring, vol. 10, no. 3, pp. 261–274, 2011. View at: Publisher Site  Google Scholar
 J. Chen, “Application of empirical mode decomposition in structural health monitoring: some experience,” Advances in Adaptive Data Analysis, vol. 1, no. 4, pp. 601–621, 2009. View at: Publisher Site  Google Scholar
 J. Liu, X. Wang, S. Yuan, and G. Li, “On HilbertHuang transform approach for structural health monitoring,” Journal of Intelligent Material Systems and Structures, vol. 17, no. 89, pp. 721–728, 2006. View at: Publisher Site  Google Scholar
 L. Lin and F. Chu, “Feature extraction of AE characteristics in offshore structure model using HilbertHuang transform,” Measurement, vol. 44, no. 1, pp. 46–54, 2011. View at: Publisher Site  Google Scholar
 W.K. Hsu, D.J. Chiou, C.W. Chen, M.Y. Liu, W.L. Chiang, and P.C. Huang, “Sensitivity of initial damage detection for steel structures using the HilbertHuang transform method,” Journal of Vibration and Control, 2012. View at: Publisher Site  Google Scholar
 Y. L. Xu, S. W. Chen, and R. C. Zhang, “Modal identification of Di Wang building under Typhoon York using the HilbertHuang transform method,” Structural Design of Tall Buildings, vol. 12, no. 1, pp. 21–47, 2003. View at: Google Scholar
 A. Bouchikhi, A.O. Boudraa, S. Benramdane, and E.H. Diop, “Empirical mode decomposition and some operators to estimate instantaneous frequency: a comparative study,” in Proceedings of the 3rd International Symposium on Communications, Control, and Signal Processing (ISCCSP '08), pp. 608–613, March 2008. View at: Publisher Site  Google Scholar
 S. Saadat, M. N. Noori, G. D. Buckner, T. Furukawa, and Y. Suzuki, “Structural health monitoring and damage detection using an Intelligent Parameter varying (IPV) technique,” International Journal of NonLinear Mechanics, vol. 39, no. 10, pp. 1687–1697, 2004. View at: Publisher Site  Google Scholar
 Y. Xia, System identification and damage detection of nonlinear structures [Ph.D. dissertation], University of California, Irvine, Calif, USA, 2011.
 A. Gupta and H. Krawinkler, “Seismic demands for performance evaluation of steel moment resisting frame structures (SAC Task 5.4.3),” John A. Blume Earthquake Engineering Research Center Report 132, Stanford University, 1999. View at: Google Scholar
 N. E. Huang, S. R. Long, and Z. Shen, “The mechanism for frequency downshift in nonlinear wave evolution,” Advances in Applied Mechanics, vol. 32, pp. 59–117, 1996. View at: Publisher Site  Google Scholar
 N. E. Huang, Z. Shen, S. R. Long et al., “The empirical mode decomposition and the Hubert spectrum for nonlinear and nonstationary time series analysis,” Proceedings of the Royal Society A, vol. 454, no. 1971, pp. 903–995, 1998. View at: Google Scholar
 P. Maragos, J. F. Kaiser, and T. F. Quatieri, “On amplitude and frequency demodulation using energy operators,” IEEE Transactions on Signal Processing, vol. 41, no. 4, pp. 1532–1550, 1993. View at: Publisher Site  Google Scholar
 F. Salzenstein, A.O. Boudraa, and J.C. Cexus, “Generalized higherorder nonlinear energy operators,” Journal of the Optical Society of America A, vol. 24, no. 12, pp. 3717–3727, 2007. View at: Publisher Site  Google Scholar
 N. E. Huang, Z. Wu, S. R. Long, K. C. Arnold, X. Chen, and K. Blank, “On instantaneous frequency,” Advances in Adaptive Data Analysis, vol. 1, no. 2, pp. 177–229, 2009. View at: Publisher Site  Google Scholar
 Uniform Building Code (UBC), “Uniform building code,” in Proceedings of the International Conference of Building Officials (ICBO '94), Whittier, Calif, USA, 1994. View at: Google Scholar
 FEMA267, Interim Guidelines, Inspection, Evaluation, Repair, Upgrade and Design of Welded Moment Resisting Steel Structures, The SAC Joint Venture for the Federal Emergency Management Agency, Washington, DC, USA, 1995.
 F. McKenna, G. L. Fenves, and M. H. Scott, OpenSees: Open System for Earthquake Engineering Simulation, Pacific Earthquake Engineering Center, University of California, Berkeley, Calif, USA, 2006.
 MATLAB User’s Guide, vol. 4, The MathWorks Inc., Natick, Mass, USA, 1998.
 V. J. Hodge and J. Austin, “A survey of outlier detection methodologies,” Artificial Intelligence Review, vol. 22, no. 2, pp. 85–126, 2004. View at: Publisher Site  Google Scholar
 D. Vamvatsikos and C. Allin Cornell, “Incremental dynamic analysis,” Earthquake Engineering and Structural Dynamics, vol. 31, no. 3, pp. 491–514, 2002. View at: Publisher Site  Google Scholar
 FEMA, Recommended Seismic Design Criteria for New Steel MomentFrame Buildings, FEMA350, SAC Joint Venture, Federal Emergency Management Agency, Washington, DC, USA, 2000.
 N. Cheraghi and F. Taheri, “A damage index for structural health monitoring based on the empirical mode decomposition,” Journal of Mechanics of Materials and Structures, vol. 2, no. 1, pp. 43–62, 2007. View at: Google Scholar
 J. C. Dunn, “A fuzzy relative of the ISODATA process and its use in detecting compact wellseparated clusters,” Journal of Cybernetics, vol. 3, no. 3, pp. 32–57, 1973. View at: Publisher Site  Google Scholar
 P. J. Rousseeuw, “Silhouettes: a graphical aid to the interpretation and validation of cluster analysis,” Journal of Computational and Applied Mathematics, vol. 20, pp. 53–65, 1987. View at: Google Scholar
 FEMA, Quantification of Building Seismic Performance Factors, FEMA P695, Applied Technology Council for the Federal Emergency Management Agency, Washington, DC, USA, 2009.
 C. A. Cornell, F. Jalayer, R. O. Hamburger, and D. A. Foutch, “Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines,” Journal of Structural Engineering, vol. 128, no. 4, pp. 526–533, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Ehsan Darvishan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.