Seismic Damage Detection ofMoment Resisting Frame Structures Using Time-Frequency Features

To detect seismic damage of moment resisting frame (MRF) structures, a data-driven method using the fractal dimension (FD) of time-frequency feature (TFF) of structural seismic dynamic responses at measured stories is extended and refined. +e TFF is defined as the real part of Gabor wavelet transform of translational interstory displacement, and FD is used to give a quantitative value to describe the calculated TFF. Static condensation method is first used to reduce the degrees-of-freedom (DOFs) of MRF and to express the rotational displacements using translational displacements. For linearMRF, the FDs of TFFs at all stories are the same using the definition of TFF and modal superposition principle. For damaged MRF with plastic hinges at the ends of beams and columns, the force analogy method is implemented to establish transformation matrix from plastic hinge rotations to translational interstory inelastic displacements. Due to the sparseness of the transformation matrix, plastic hinges only generate interstory inelastic displacements, which are low-frequency contents, in the vicinity of plastic hinges. Correspondingly, the FDs of TFFs of interstory displacements with inelastic component are different from the FDs of TFFs of the interstory displacements that do not contain inelastic component. A numerical simulation on a 16-story MRF was conducted.+e simulation included 10 cases such as no damage or linear structure, plastic hinges in single-story beams, plastic hinges in single-story columns, plastic hinges in single-story beams and columns, and plastic hinges in multiple story beams and columns. +e robustness to measurement noise was also investigated. +e seismic damage detection results demonstrated that the proposed method was capable of locating the stories where the plastic hinges occurred.


Introduction
After big earthquakes occur, the national and local governments and owners, in order to take actions of postearthquake emergency rescue and reconstruction, need to know the safety and applicability of building structures damaged by strong ground motions [1].Traditional postearthquake evaluation method for building structures is to arrange a large number of structural experts to inspect tens of thousands of damaged structures, check the damage of beams, columns, walls, stories, and roofs, and mark the damaged structures as intact, slightly damaged, medium damage, severe damage, and collapse categories [2]. is procedure is time consuming and laborious, and it is difficult to assess the internal damage which cannot be seen by eyes and the assessment is not objective.For example, after the Kobe earthquake in 1995, postearthquake investigation of 46,000 damaged structures in one disaster-stricken area took 6,000 structural experts spending up to three weeks to complete just a preliminary assessment of the damage structures [3].
is seriously delayed the emergency rescue and reconstruction after the earthquake.In addition, the damaged structures may collapse in the subsequent aftershocks and pose a life threat to the investigating structural experts.On the other hand, vibrationbased structural health monitoring (SHM) provides a new rapid and effective approach for structural health condition determination and seismic damage assessment [4][5][6][7][8][9][10].SHM deploys sensors to acquire the dynamic responses of the prototype structures under real service environments and external loads, extracts damage-sensitive features hidden in the measurements, and provides safety information of the monitored structures in real time or near real time.SHM does not induce additional damage or negative effect to the structural and nonstructural components and does not affect the normal operation and service of the monitored structure.
e assessment of structural damage by SHM is automatic, rapid, and effortless.Hence, SHM provides promising means for assessing the safety condition of structures prone to earthquakes, especially heritage structures such as towers and churches [11][12][13][14].
Damage detection is the core part of SHM, and the first important information needed for damage assessment and seismic retrofit for structures shaken by strong ground motions is to determine whether the structures is damaged and where the locations of damage are.Structural seismic damage is divided into three levels as component-level damage, storylevel damage and structure-level damage.Due to the complexity and huge degrees-of-freedoms (DOFs) of the structure, sensors are only installed on very limited stories, which makes the current seismic damage detection methods for real structures mainly focus on the story-level.e structure excited by strong ground motion will be in the nonlinear stage because the structural material yields under high stress and large deformation.e plastic hinges will exist at the ends of beams and columns, the cracks at components like beams, columns, and walls will open and close, and the forcedisplacement curves will show obvious hysteresis characteristics.
ese kinds of nonlinear behaviour of monitored structure make the damage detection methods based on linear system theory not suitable for detecting structural seismic damage.us, nonlinear indicator function methods, which belong to data-driven approaches, are proposed to detect these kinds of damages by extracting the nonlinear features in the dynamic vibration measurements.Farrar et al. gave a comprehensive review on damage detection using nonlinear system identification and nonlinear indicator functions approaches [15].Due to the nonstationarity of measured structural responses, time-frequency distributions such as Wigner-Ville distribution and time-scale transforms such as wavelet transform have been implemented and widely used to analyse and interpret the nonstationary signals in timefrequency or time-scale domain.Staszewski and Robertson [16], Robertson and Basu [17], and Ceravolo [18] gave comprehensive reviews on using time-frequency and time-scale transforms for SHM.In this paper, references focused on wavelet transform were presented.Robertson et al. used the Holder index derived from the continuous wavelet transform to characterize the singularity of the structure due to nonlinearity and identified damage such as structural looseness, collisions, and cracks [19].Civera et al. proposed a novel approach using bispectral analysis and neural network to identify the locations of breathing cracks in a cantilever beam [20].Liew and Wang implemented wavelet theory to identify crack damage in structures, where discrete wavelet transform was used to decompose the vibrational shape of a simply supported cracked beam, and the crack location was identified by the variation of the imaginary parts of the 8 th and 9 th order wavelet coefficients [21].Loh et al. used the first-order discrete wavelet coefficients, wavelet packet energy, and Holder index to identify the seismic damage of a single-story reinforced concrete structure in a shaking table test [22,23].Vafaei and Adnan used continuous wavelet transform to identify the simulated seismic damage of an airport traffic control tower and found that continuous wavelet transform had stronger noise immunity than discrete wavelet transform and did not need as high sampling frequency as discrete wavelet transform [24].Kumar and Zanotti Fragonara proposed a time-frequency method to identify nonlinear characteristics of a cable element [25].Nagarajaiah and Basu summarized the pros and cons of short time Fourier transform, wavelet transform, and Hilbert-Huang transform application in structural seismic damage detection [26].e features of nonlinear system such as fractal dimension (FD) and chaos were also used as nonlinear indicator function [27][28][29][30].Li et al. proposed a datadriven approach using the time-frequency feature (TFF) for shear-type buildings where the FD of TFF at each DOF for linear system was derived with the help of modal superposition principle, and the FD of TFF for each substructure in nonlinear system was derived with the help of substructure method [31].Carrillo and Avila used FD of breathing cracks to assess the seismic damage of reinforced concrete walls [32].
In this paper, a data-driven approach using the timefrequency feature (TFF) of the interstory displacement is further developed for detecting story-level damage of moment resisting frame (MRF) structure subjected to earthquake.
is nonlinear indicator function method was proposed by the author to detect seismic damage of sheartype buildings.For the integrity and understandability, some important principles and equations are also given in the paper.In Section 2.1, the real part of the Gabor wavelet transform is introduced to stand for the TFF.Section 2.2 uses the FD to quantify the random fractal set of TFF.Section 2.3 deduces the FD property of linear MRF using the definition of TFF and modal superposition technique.Section 2.4 deduces the FD property of damaged MRF with plastic hinges at the ends of beams and columns using the force analogy method, which establishes the relation between interstory inelastic displacements and plastic hinge rotations.Section 3 presents numerical simulations on a 16story 3-bay MRF to demonstrate the validity of the approach, and finally Section 4 gives the conclusions.

Seismic Damage Detection Method of MRF
Using the Time-Frequency Feature (TFF) where x(t) is the time domain signal, b is time, a is scale, and φ * is the conjugate of Gabor wavelet given by where j is the imaginary unit, σ is the frequency bandwidth, and f 0 is the central frequency.e frequency in the timefrequency domain is given by the scale a, sampling frequency f s , and central frequency f 0 as shown in the following equation: As deduced in Reference [31], the TFF is defined as the real part of Gabor wavelet transform in a chosen scale or frequency band given by To keep the linear superposition of TFF and to suppress the interference between different modes, the real part, not the norm of the wavelet coefficients of the signal, is used in this study.Meanwhile, for most damage detection purpose, the scale is chosen within the fundamental frequency of the analysed structure.e TFF of the damaged structure will show much energy in the frequency band below the structural fundamental frequency.
e TFF can be treated a matrix whose columns stand for time and rows for frequency, and it can be visualized as a three-dimensional (3D) surface.

Fractal Dimension (FD) to Quantify Time-Frequency Feature (TFF).
To quantify the calculated TFF of a signal, FD is introduced.FD is an indicator for describing the complexity or irregularity of natural objects, textures, images, and signals [34].Unlike deterministic fractal such as Koch curve or Sierpinski Triangle which has FD at any scale, the TFF is a kind of random fractal set which has FD only in specific scale.e FD of TFF is calculated using the box counting method given by dim where c 1 is the regression coefficient given by where δ is the side length of small box and N δ is the minimal number of δ-length small box to cover completely the 3D surface of TFF.

FD Property of
where X m denotes displacement vector in the translational or master degree-of-freedom (DOF) and X s displacement vector in the rotational or secondary DOF; _ X m , _ X s are the corresponding translational and rotational velocity vectors; € X m , € X s are the corresponding translational and rotational acceleration vectors; M mm , K mm are the translational block matrices of mass, stiffness matrices; K ss is the rotational stiffness matrices; K ms and K sm are the off-diagonal stiffness matrices; α, β are the Rayleigh proportional damping parameters, and € x g is the ground motion acceleration.From Equation ( 7), the rotational or secondary displacement vector can be expressed as Substituting Equation ( 8) into (7), the static condensed motion equation is given by where M e , C e , K e are the effective mass, damping, and stiffness matrices giving by M e � M mm , K e � K mm − K ms K −1 ss K sm , and C e � αM e + βK e .e relative translational displacements of the condensed linear MRF can be given as follows with the help of modal decomposition technique: where is the ith mode shape vector, and q i is the ith generalized modal displacement.Based on the definition of TFF, the following equation can be deduced: Using the property of FD, which is that the FD of a constant times a fractal is the same as the original fractal, the following equation can be deduced: where x mi,j is the partial relative displacement at ith translational DOF contributed by the jth mode.From Equation ( 12), it can be seen that FD of TFF of all relative translational displacements of a MRF contributed by any

Shock and Vibration
given mode is identical to the FD of that generalized mode displacement.

FD Property of Damaged MRF Structure with Beam or
Column Plastic Hinges.When MRF structure is excited by strong ground motion, plastic hinges will be formed at the ends of beams and columns.It is these plastic hinges that make the horizontal displacement of the structure contain inelastic displacements.e force analogy method establishes an explicit expression between the inelastic displacement and the plastic hinge rotation of the beams and columns.e force analogy method is a numerical method for solving dynamic responses of nonlinear structures.It was rst proposed by Lin and developed and re ned by Hart [35].e traditional method to solve the structural dynamic response is to constantly correct the instantaneous sti ness of the structure.However, the force analogy method adjusts the inelastic displacement so that the total restoring force equals the total restoring force of the modi ed sti ness method.e translational displacement of the reduced MRF is given by where X e (t) is the elastic displacement and X ie (t) is the inelastic displacement which is introduced by damage at the ends of beams and columns.e plastic rotation vector of the plastic hinges at the ends of beams and columns is de ned as where θ ie i (t) is the plastic rotation at the ith plastic hinge and n is the total number of plastic hinges.
e nonequilibrium restoring force in the translational DOF generated by the plastic rotations at the ends of beams and columns is given by −K P Θ ie (t) where K P is called member force recovery matrix.
e unbalanced bending moment generated at the beam-column joint is given by −K R Θ ie (t) where K R is called member restoring force matrix.K P and K R are the combination of corresponding matrix of each element.For an element with plastic hinges at both ends shown in Figure 1, the corresponding recovery matrices is given by From Figure 1 and Equation ( 15), it can be seen that only plastic hinges at the ends of columns produce horizontal inelastic displacement, and the plastic hinges at both ends of the beams do not produce horizontal inelastic displacement.
e nonequilibrium restoring force caused by plastic hinges is balanced by the restoring force due to inelastic displacement: Equation ( 16) gives the explicit expression between the inelastic displacement and the plastic hinge rotation of the beams and columns.us, the inelastic displacement at the master DOF of the reduced order frame structure can be deduced as where K e P I −(K −1 ss K sm ) T K P is the member force recovery matrix for the reduced order structure.If denoting a matrix D as then the explicit expression between interstory inelastic displacement in the master DOF and the plastic hinge rotation is given by where T is a sparse transformation matrix from plastic hinge rotation to interstory inelastic displacement.
When the MRF structure is in linear stage and there is no plastic hinge at the ends of beams and columns, then, based on Equation (19), there is no interstory inelastic displacement in the master or translational DOFs, and the TFFs of interstory displacement in the master DOFs have the same FD.When plastic hinges occur at the ends of local beams and columns, due to sparseness of transformation matrix T, plastic hinges only generate inelastic displacements in the vicinity of the stories.Since the inelastic displacement belongs to low frequency signal, the TFFs of the interstory displacements with inelastic displacements are di erent from the TFF of the interstory displacements that do not contain inelastic displacements.Correspondingly, the FD of the TFF is also di erent from that of other stories, so as to determine the story location where exist plastic hinges at the ends of columns or beams.Shock and Vibration

MRF Structural Model.
e 3D nite element model was established by an open source software called the Open System for Earthquake Engineering Simulation (OpenSees) [36].e simulated structural model was a 16-story steel MRF structure as shown in Figure 2, which was a 1 : 8 scale model of a real structure.e height of the story was 0.5 m, and the total height of the structure was 8.0 m. e MRF had three bays along the shaking direction with each span 0.75 m and one bay in the perpendicular direction with span 1.0 m. e size of square column was 50 × 50 × 4 mm, the size of the I-beam was 60 × 4 × 30 × 5 mm, the density of the steel was 7800 kg/m 3 , the modulus of elasticity was 2.06 × 10 5 N/mm 2 , and the additional mass at each story was 300 kg. e beam and column elements were all made of Nonlinear Beam-Column elements.Inelastic ber cross section was adopted for the element section.Axial hardening stress strain (Uniaxial Steel01) was chosen, the ratio of after-yield to preyield sti ness was 0.01, and the Poisson ratio was taken as 0.3.e damping of the entire structure was Rayleigh proportional damping, and the rst-order and third-order modal damping ratios were 0.02.e 3D nite element model ignored the rotational mass and vertical mass.e rigid diaphragm was turned o (RigidDiaphram OFF), and 5 Gaussian integration points were set in nonlinear curvature distribution of the element.Input ground motion used uniform excitation pattern, and the calculation step interval was 0.01 s.
e input ground motions were chosen from PEER's NGA database.e record at El Centro site in October 15, 1979, Imperial Valley earthquake, was chosen as far-eld strong ground motion.
e record code was E12140.e PGA was 143 gal, the peak frequency was 1.76 Hz, and the original sampling interval was 0.05 s.
is record was abbreviated as El Centro.e record at Ridge RT site, in January 18, 1994, Northridge earthquake, was chosen as near-eld strong ground motion.e record code was ORR090.e PGA was 580 gal, the peak frequency was 1.22 Hz, and the original sampling interval was 0.02 s.
is record was abbreviated as Northridge.In order to simulate no damage or di erent damage in speci c stories, the chosen ground motions were scaled up or down, and the yielding stress of beams and columns were carefully set.e simulation included 10 cases such as no damage or linear structure, plastic hinges in single-story beams, plastic hinges in single-story columns, plastic hinges in single-story beams and columns, plastic hinges in multiple story beams, and columns as shown in Table 1.e yielding stress at nonweak elements was set at 300 MPa so that these elements would be not yielded.

Structural Seismic Damage Detection Results.
To detect damage and locate the weak stories, the translational displacements at each translational DOF were extracted from the OpenSees dynamic response analysis.e plastic rotation at the ends of beams and columns were also extracted to demonstrate where the weak stories were truly located.e Gabor wavelet transform program used was WAVELAB, a MATLAB wavelet toolkit provided by David Donoho's group at Stanford University [37].e FD calculation program used was FRA-CLAB, a MATLAB fractal toolkit provided by Canus et al. [38].
e central frequency of the Gabor wavelet was set 2.4 Hz, which was close to the fundamental frequency of the simulated MRF structure, and the bandwidth was set 1.

Linear Structure.
In case 1 when the MRF structure was excited by PGA 28.6 gal El Centro ground motion, the structure was in linear range without yielded element.e rst two modal frequencies were 2.13 Hz and 6.53 Hz, respectively.e contour map of time-frequency distribution, Fast Fourier Transform (FFT) spectra, and the interstory displacement is shown in Figure 3.For simplicity and concision, only the plots at 1 st , 7 th , and 15 th stories were shown here.It could be seen from Figure 3 that the time-frequency distributions at these three stories were quite similar to each other.For the cut-o frequency of 0.1-4.0Hz, the log-log plot of the box counting of the TFF at the three stories is shown in Figure 4. e double logarithmic graph of the TFFs in the scale interval 2 −5 -2 −2 was approximately linear, indicating a fractal feature in this scale, and the corresponding FDs obtained by the least-squares method were 2.049, 2.048, and 2.049, respectively.e TFFs of the interstory displacement of other stories also had fractal features in this scale interval.e corresponding FDs are shown in Figure 5, and they were equal.

Damaged Structure with Plastic Hinges in Single-Story
Beams.In case 2 when the MRF structure was excited by PGA 429 gal El Centro ground motion, the structure was damaged with plastic hinges in the beams at the 1 st story.e rotation of plastic hinges at the ends of columns and beams Shock and Vibration is shown in Figure 6, where columns A-D referred to the left column, left middle column, right middle column, and right column of the frame shown in Figure 2 and beams E-G referred to the left side beam, the centre beam, and the right side beam of the frame shown in Figure 2, and counts referred to the ends of the beams and columns.For example, counts in the column counted as 1 was to indicate the lower end of a column in the rst story and 2 indicated the upper end of a column in the rst story.e beam counted as 5 was to indicate the left end of a beam on the 3 rd story, 6 indicated the right end of a beam on the 3 rd story, and so on so forth.
From Figure 6, it could be seen that the beams at the 1 st story were damaged, and the columns at the 1 st story and beams and columns at other stories did not yield.e contour map of time-frequency distribution, FFTspectra, and the interstory displacement is shown in Figure 7, from which it could be seen that time-frequency distribution at the 1 st story was di erent from those at other stories, especially between 1.2 and 1.5 Hz in the frequency range and 10 and 12 s in the time range.e corresponding FDs in the scale interval 2 −5 − 2 −2 were 2.092, 2.073, and 2.073 as shown in Figure 8. e FDs of TFFs at other stories are shown in Figure 9. e FDs of TFF at

6
Shock and Vibration the 1 st and 2 nd stories were a little di erent from those at other stories, and this is because plastic hinges introduced inelastic displacement at the 1 st story as inferred from Equation (19).

Damaged Structure with Plastic Hinges in Single-Story
Columns.In case 3 when the MRF structure was excited by PGA 429 gal El Centro ground motion, the structure was damaged with plastic hinges in the columns at the 1 st story.e rotation of plastic hinges at the ends of columns and beams is shown in Figure 10, from which it could be seen that the columns at the 1 st story were damaged and the beams at the 1 st story and beams and columns at other stories did not yield.e contour map of time-frequency distribution, FFT spectra, and the interstory displacement is shown in Figure 11, from which it could be seen that timefrequency distribution at the 1 st story was quite di erent from those at other stories, especially between 1.2 and 1.7 Hz in the frequency range and 8 and 12 s in the time range.e corresponding FDs in the scale interval 2 −5 − 2 −2 were 2.134, 2.070, and 2.071 as shown in Figure 12. e FDs of TFFs at other stories are shown in Figure 13.e FD of TFF at the 1 st story was quite di erent from those at other stories, and this could be interpreted by the transformation matrix T given by e coe cients at the rst row were much higher than those at other rows, which meant the plastic hinges in the columns at the 1 st story mainly caused interstory inelastic displacement at the 1 st story, and thus, the FD of TFF at the 1 st story was di erent from others.

Damaged Structure with Plastic Hinges in Single-Story
Beams and Columns.In case 4 when the MRF structure was excited by PGA 429 gal El Centro ground motion, the structure was damaged with plastic hinges in the beams and columns at the 1 st story.e rotation of plastic hinges at the ends of columns and beams is shown in Figure 14, from which it could be seen that the beams and columns at the 1 st story were damaged and the beams and columns at other stories did not yield.e contour map of time-frequency distribution, FFT spectra, and the interstory displacement is shown in Figure 15, from which it can be seen that time-frequency distribution at the 1 st story was quite di erent from those at other stories, especially between 0 and 1.0 Hz in the frequency range and 5 and 20 s in the time range.e corresponding FDs in the scale interval 2 −5 -2 −2 were 2.172, 2.070, and 2.071 as shown in Figure 16.e FD of TFF at other stories is shown in Figure 17. e FDs of TFFs at the 1 st story were quite di erent from those at other stories, and it was because the plastic hinges in the beams and columns at the 1 st story mainly caused interstory inelastic displacement at the 1 st story.

Damaged Structure with Plastic Hinges in Multiple-
Story Beams and Columns.In case 5 when the MRF structure was excited by PGA 715 gal El Centro ground motion, the structure was damaged with plastic hinges in the beams and columns at multiple stories.e rotation of plastic hinges at the ends of columns and beams is shown in Figure 18, from which it could be seen that the columns at Stories 1-6 and beams at Stories 1-10 were damaged.e contour map of time-frequency distribution, FFT spectra, and the interstory displacement is shown in Figure 19, where the plots at Stories 1-7 and 10 and 15 were shown here to demonstrate the behaviour of widespread damages in the frame.From Figure 19, it can be seen that timefrequency distribution at Stories 1-6 were quite di erent from those at other stories, especially between 0 and 1.0 Hz in the frequency range and 5 and 25 s in the time range.e low frequencies generated by plastic hinges at beams and columns were obvious in the lower stories especially in Stories 1-6, and disappeared in the upper stories.e corresponding FDs in the scale interval 2 −5 -2  21. e FDs of TFFs at Stories 1-6 were quite di erent from those at other stories, and it was because the plastic hinges in the beams and columns at Stories 1-6 mainly caused interstory inelastic displacement at Stories 1-6.

Shock and Vibration
To investigate the robustness of the approach to measurement noise, white Gaussian noise with 5%, 10%, 15%, and 20% root mean square (RMS) noise-to-signal were added to the relative displacements at each story.e results for Case 5 are shown here.Time histories, TFFs, and the FFT spectra of the interstory displacement at Story 1   Shock and Vibration noise levels were quite the same as those without noise.When the noise level became 15%, the FDs at lower stories deviated from those values without noise, but it still could be deduced that the lower stories from Stories 1-10 were damaged and the upper stories from Stories 11-16 were not damaged.When the noise level became 20%, the FD along the structure was a zigzag line which may be caused by damage or high-level noise, and the damage could not be determined.

Detection Result for Other Cases Studies under Other Ground Motions.
e results of the FDs of the TFFs of the MRF structure for Cases 6-10 under Northridge ground  To simulate di erent damage level, these ground motions were scaled up or down, and the time interval was set as 0.0035 second.For the cases when the yielding stresses at weak components were chosen as beams and columns 205 MPa, the detection results are shown in Figure 26.e results when the structure was subjected to Chichi, Kobe, and Park eld ground motions were similar as the results when subjected to El Centro or Northridge ground motions, and it could be seen that the lower stories were damaged in these cases.However, the FD of TFFs under Taft ground motion along the structure was a zigzag line, and it could be seen that the structure was damaged but the location of damage could not be determined.According to the de nition of TFF given by Equation (4), the TFF was carefully chosen so that the frequency in the TFF was around or below the fundamental frequency of the analysed structure.If the chosen TFF was contaminated by higher  -12 -10 -  Shock and Vibration order frequency energy, the FD based approach would not give correct detection results.
e rst two modal frequencies of the MRF was 2.13 and 6.53 Hz, and it could be seen from Figure 26 that Taft ground motion had relatively more high frequency energy especially around the second structural modal frequency.In this case, the chosen TFF would have more energy from higher frequency contents and the corresponding FDs would not show the right damage detection result.
e above detection results' comparison under several ground motions implied that the proposed approach should be only used when the structural response was mainly contributed by the fundamental mode.For the structural response with dense mode contribution, time frequency separation should be rst implemented before using the FD of TFF technique to detect damage.

Conclusions
In this paper, an extension of a data-driven approach proposed by the authors for detecting damage of shear-type   Shock and Vibration building structure is re ned and further developed to detect seismic damage of MRF subjected to far-eld or near-eld strong ground motions.e method is one kind of vibration-based approach, which uses the measured displacements to detect damage and locate stories with plastic hinges at the ends of beams and columns.e principle is that the FD of TFF of interstory displacements with inelastic component are di erent from the FD of TFF of the interstory displacements that do not contain inelastic component.e numerical simulation on a 16-story MRF indicates that the method can locate the stories in which there are plastic hinges at the end of beams and columns.e damaged story location accuracy is better when the plastic hinges are at the end of columns than those at the end of  Shock and Vibration beams.As the seismic damage of real MRF structure will introduce plastic hinges at the ends of beams or columns, the method may be utilized for detecting and localizing these kinds of damage.Further study will include experiment validation on a MRF in shaking table test, and testing on insitu MRF under real ground motions will also be conducted.

Figure 1 :
Figure 1: Element with plastic hinges at both ends.

Figure 4 :
Figure 4: Log-log plot of the box counting of the TFF of the structure under PGA 28.6 gal El Centro ground motion (Case 1).(a) Story 1, (b) Story 7, and (c) Story 15.

Figure 5 :
Figure 5: FDs of the TFFs of interstory displacements of the structure under PGA 28.6 gal El Centro ground motion (Case 1).

Figure 6 :(Figure 10 :Figure 11 :
Figure 6: Rotation of plastic hinges at the ends of columns and beams under PGA 429 gal El Centro ground motion (Case 2).(a) Rotation of plastic hinges at the ends of columns.(b) Rotation of plastic hinges at the ends of beams.

Figure 9 :
Figure 9: FDs of the TFFs of interstory displacements of the structure under PGA 429 gal El Centro ground motion (Case 2).

Figure 13 :Figure 12 :
Figure 13: FDs of the TFFs of interstory displacements of the structure under PGA 429 gal El Centro ground motion (Case 3).

Figure 14 :
Figure 14: Rotation of plastic hinges at the ends of columns and beams under PGA 429 gal El Centro ground motion (Case 4).(a) Rotation of plastic hinges at the ends of columns.(b) Rotation of plastic hinges at the ends of beams.

Figure 17 :Figure 18 :Figure 19 Figure 19 :
Figure 17: FDs of the TFFs of interstory displacements of the structure under PGA 429 gal El Centro ground motion (Case 4).

Figure 21 :
Figure 21: FDs of the TFFs of interstory displacements of the structure under PGA 715 gal El Centro ground motion (Case 5).

Figure 26 :
Figure 26: FDs of the TFFs of interstory displacements of the structure under strong ground motions when all the yielding stress of beams and columns as 205 MPa.(a) Chichi, (b) Kobe, (c) Park eld, and (d) Taft.

Table 1 :
Location of nonlinearity and input ground motion.