Next Article in Journal
Characterization of Volume Gratings Based on Distributed Dielectric Constant Model Using Mueller Matrix Ellipsometry
Next Article in Special Issue
Study on Seismic Performance and Damage Analysis of Steel Plate Shear Wall with Partially Encased Composite (PEC) Columns
Previous Article in Journal
Sorption of Eu (III) onto Nano-Sized H-Titanates of Different Structures
Previous Article in Special Issue
Hysteretic Behavior of Steel Reinforced Concrete Columns Based on Damage Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Output-Only Parameters Identification of Earthquake-Excited Building Structures with Least Squares and Input Modification Process

1
Heilongjiang Key Laboratory of Disaster Prevention, Mitigation and Protection Engineering, College of Civil Engineering, Northeast Petroleum University, Daqing 163318, China
2
College of Civil Engineering, Hunan University, Changsha 410082, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(4), 696; https://doi.org/10.3390/app9040696
Submission received: 19 January 2019 / Revised: 12 February 2019 / Accepted: 14 February 2019 / Published: 18 February 2019
(This article belongs to the Special Issue Structural Damage Detection and Health Monitoring)

Abstract

:
Damage detection and system identification with output-only information is an important but challenging task for ensuring the safety and functionality of civil structures during their service life. In this paper, a relatively simple and efficient iteration identification method consisting of the least squares estimation (LSE) technique and an input modification process is proposed for the simultaneous identification of structural parameters and the unknown ground motion. The spatial distribution characteristics of ground acceleration on earthquake-excited building structures are considered as additional information for parameters identification in each iterative step. First, the unknown input is estimated using the measured responses and the initial guesses of the structural parameters. The estimated input is then modified on the basis of the property of its spatial distribution. This modified input is further employed for providing the updated estimation of structural parameters. The iterative procedure would continue until the preset convergence criterion is satisfied. The accuracy of the proposed approach is numerically validated via a shear building model under the El Centro earthquake. The effects of signal noise, the number of sample points, and the initial guesses of structural parameters are discussed. The results show that the proposed approach can satisfactorily identify the structural parameters and unknown earthquakes.

1. Introduction

1.1. Background

Many civil structures in service may deteriorate or degrade due to many factors, such as environmental corrosion, long-term fatigue, strong winds, severe earthquakes, and so forth. The failure of these deficient structures could be catastrophic, not only in terms of loss of life and the economy, but also of subsequent social and psychological impacts [1]. It is therefore imperative to identify the parameters of the structure, detect the damage, and guide the maintenance to ensure structural safety and functionality. Moreover, since the loadings applied to the structures are of importance for structural evaluation and future life prediction, the development of loading identification techniques is also highly required.

1.2. Literature Survey

To date, a number of vibration-based system identification and damage detection methods have been developed, and review articles with different emphases can be also found [2,3,4]. In general, structural parameters identification can be performed in three different paradigms: (i) time domain (e.g., [5,6,7,8,9,10]), (ii) frequency domain (e.g., [11,12,13,14]), and (iii) time–frequency domain (e.g., [15,16,17]).
To identify structural parameters and their variations caused by damages, time domain analyses have been studied extensively. In particular, research on system identification based on least-squares estimation (LSE) has been carried out for a relatively long time, and much progress has been made. Caravani et al. [18] first proposed a time domain identification method by using recursive LSE in 1970s. By using the recursive LSE-based algorithm to upgrade the diagonal elements of the adaptation gain matrix at each time step, Lin et al. [19] proposed an on-line parametric identification algorithm for the identification of non-linear hysteretic structures. By using LSE for the coefficient estimation, Masri et al. [20] developed a system identification approach for non-linear restoring force identification. Based on the LSE approach, Yang and Lin [21] proposed an adaptive tracking technique to identify and track the system parameter variations due to damages. Yang et al. [22] further developed a sequential non-linear LSE approach for the on-line identification of structural parameters by employing the proposed adaptive tracking technique. Gui et al. [23] presented an adaptive system identification approach with the aid of the least mean square algorithm and least mean fourth algorithm. Kazemi and Arefi [24] proposed a fast iterative recursive LSE algorithm for identifying non-linear systems with output colored noise.
Although the aforementioned LSE-based approaches are capable of accurately identifying structural parameters, the external excitations in these approaches are all assumed to be available. However, in many cases, the excitations are unmeasurable or difficult to be precisely measured. Thus, the development of LSE-based algorithms for parameters identification under unknown loadings is required.
Based on recursive LSE and the extended Kalman filter (EKF), Wang and Haldar [25] proposed a system identification approach for the estimation of structural parameters under unknown excitations. Using the Taylor series approximation to transform a set of non-linear equations to linear equations, Ling and Haldar [26] developed a modified iterative LSE approach for the element-level system identification in consideration of Rayleigh damping under unknown input. Based on the LSE technique and partially known input information, Chen and Li [27] presented a time domain method for the simultaneous identification of the structural parameters and the time history of the unknown external excitations. A recursive least-squares estimation with unknown inputs (RLSE-UI) approach was proposed by Yang et al. [28] for the identification of the structural parameters and the unmeasured excitations. To reduce the number of sensors required, Yang and Huang [29] further proposed a sequential non-linear LSE method for damage identification of structures with unknown inputs and unknown outputs. With the aid of a weighted positive definite matrix and a learning coefficient, Xu et al. [30] proposed a weighted adaptive LSE approach for the identification of structural parameters and partially unknown loadings. This method was further extended for the identification of large-scale structures [31] and non-linear structures [32]. By sequentially using the EKF for the parameters identification of the concerned structure and LSE for its unknown external excitations, Lei et al. [33] proposed an approach for estimating structural parameters with limited input and output measurements. Based on the fluctuating properties of wind, Lei et al. [34] further extended their approach to identify the parameters of the high-rise wind-excited buildings. Based on recursive LSE with adaptive multiple forgetting factors, Askari et al. [35] proposed an online tracking technique for estimating the abrupt changes in structural parameters with unknown inputs.

1.3. Scope

The above-mentioned methods can all simultaneously identify the parameters of the structure and the excitation applied to it. The analytical solutions of most of these methods and their derivation, however, are rather complicated. In this paper, an iterative identification procedure based on the LSE technique and an input modification process is proposed for the simultaneous identification of the parameters of building structures and the unknown ground motion. The spatial distribution characteristics of the ground acceleration in earthquake-excited structures are used through the input modification process in each iterative step. The proposed iterative identification procedure can be described as follows. First, based on the equation of motion, the unknown ground motion is estimated using the measured responses and initial guesses of the structural parameters. The estimated input is then modified with the aid of the spatial distribution characteristics of the ground acceleration. This modified input is further employed to provide an updated estimation of the structural parameters. The iterative procedure would continue until the structural parameters satisfy the preset convergence criterion. The feasibility and accuracy of the proposed approach is numerically verified via a shear building model under the El Centro earthquake.
The framework of the paper is presented as follows: Section 2: description of the proposed iterative identification method; Section 3: investigation of the effect of some factors including the noise, the number of sample points, and the initial values of structural parameters; Section 4: concluding remarks.

2. Proposed Iterative Identification Method

The governing equation of motion of an n degree-of-freedom (DOF) linear structure under earthquake can be expressed as:
M x ¨ ( t ) + C x ˙ ( t ) + K x ( t ) = M 1 x ¨ g ( t )
where M, C, and K represent the n × n mass, damping, and stiffness matrices of the structure, respectively; x ¨ ( t ) , x ˙ ( t ) and x ( t ) denote the structural acceleration, velocity, and displacement responses, respectively; 1 denotes an n × 1 matrix with all the elements being equal to 1; and x ¨ g ( t ) is the time history of ground acceleration. In this study, the Rayleigh damping assumption is adopted:
C = αM + βK
where α and β are mass-proportional and stiffness-proportional damping coefficients, respectively. Substituting Equation (2) into Equation (1) results in
M x ¨ ( t ) + ( α M + β K ) x ˙ ( t ) + K x ( t ) = M 1 x ¨ g ( t )
Here, the mass matrix M is assumed to be known, and then Equation (3) can be rearranged as:
α M x ˙ ( t ) + β K x ˙ ( t ) + K x ( t ) = M 1 x ¨ g ( t ) M x ¨ ( t )
Notably, since both β and K are the parameters to be identified, their coupled products would introduce a source of non-linearity into the identification procedure when the LSE technique is used. To handle this situation, Equation (4) can be rewritten as:
α M x ˙ ( t ) + Q x ˙ ( t ) + K x ( t ) = M 1 x ¨ g ( t ) M x ¨ ( t )
where Q can be viewed as a transformation matrix (Q = βK) with the same dimension of K.
By rearranging Equation (5), the identification equation can be expressed as follows:
H ( t ) θ = P ( t )
in which H(t) is the response matrix with the dimension of ns × h (ns = n × s); h and s are defined as the number of unknown coefficients and sample points to be used, respectively; P(t) is ns × 1 inertia force vector; θ is the h × 1 (h = 2 × n + 1) unknown coefficients vector defined as:
θ = [ k 1 , , k n , α , q 1 , , q n ] T
where ki (i = 1, …, n) denotes the stiffness of the i-th floor; qi = β × ki (i = 1, …, n) is the element-wise coefficient in Q matrix.
The unknown parameters can then be estimated by LSE as follows:
θ = [ H ( t ) T H ( t ) ] 1 H ( t ) T P ( t )
It follows from Equation (8) that the input should be known for the identification. However, in some practical situations, for example an earthquake-excited building, the measurements of loadings are sometimes unavailable. Therefore, the parameters cannot be directly identified according to Equation (8) in this case. Additional information or supplementary equations are required for the identification.
Based on the spatial distribution characteristics of the ground motion in earthquake-excited structures, an input modification process is proposed in this study. A time-domain iterative identification procedure is then developed for the simultaneous identification of structural parameters and unknown ground motion. The details are given below:
Step 1: Assign the initial values of the unknown parameters named as θ ˜ .
Step 2: Estimate ground acceleration by using θ ˜ and structural response measurements. Based on the Equation (5), the ground acceleration can be estimated as:
X ¨ ˜ g ( t ) = M 1 ( M x ¨ ( t ) + α ˜ M x ˙ ( t ) + Q ˜ x ˙ ( t ) + K ˜ x ( t ) )
where α ˜ , Q ˜ and K ˜ are the estimated values of α, Q and K, respectively. They can be determined with the aid of θ ˜ and structural topology. Notably, the estimated ground acceleration X ¨ ˜ g obtained from Equation (9) is an n × s matrix in the form of
X ¨ ˜ g ( t ) = [ x ¨ ˜ g 1 ( t 1 ) x ¨ ˜ g 1 ( t 2 ) x ¨ ˜ g 1 ( t s ) x ¨ ˜ g n ( t 1 ) x ¨ ˜ g n ( t 2 ) x ¨ ˜ g n ( t s ) ]
where x ¨ ˜ g i ( t j ) (i = 1, …, n; j = 1, …, s) denotes the ground acceleration distributed on the i-th DOF at the time step of tj.
Step 3: Modify the estimated ground acceleration and form the corresponding inertia force vector P ^ ( t ) . For earthquake-excited building structures, the ground acceleration should be equally distributed on each floor of the building structure. It can be concluded that, if the unknown parameters are not stably converged to the actual ones, the estimated ground acceleration values shown in Equations (9,10) for different DOFs would not be close to each other. Based on this distribution characteristic, the following input modification process is presented in this study:
x ¨ ^ g ( t j ) = 1 n i = 1 n x ¨ ˜ g i ( t j ) , ( j = 1 ,   2 ,   , s )
where x ¨ ^ g stands for the modified ground acceleration. Obviously, the mean values of the estimated ground acceleration are used as the modified ones, and the corresponding inertia force vector can be determined as:
P ^ ( t ) = [ P ^ ( t 1 ) T P ^ ( t 2 ) T P ^ ( t s ) T ] T
where
P ^ ( t j ) = M 1 x ¨ ^ g ( t j ) M x ¨ ( t j ) ,   ( j = 1 ,   2 ,   , s )
Step 4: Form the response matrix H(t), and obtain the updated parameters θ ^ according to Equation (8).
θ ^ = [ H ( t ) T H ( t ) ] 1 H ( t ) T P ^ ( t )
Step 5: Calculate the error between θ ^ and θ ˜ . Let θ ^ be θ ˜ , and repeat Steps 2–5 until the following convergence criterion is satisfied:
| θ ^ l θ ˜ l θ ˜ l | ε   ( l = 1 ,   2 ,   ,   h )
where θ ^ l and θ ˜ l are the l-th element of θ ^ and θ ˜ , respectively; ε is a preset value.
When the iteration procedure is complete, the ground acceleration can be directly found from Equation (11). Moreover, since both ki and qi are identified, the damping coefficient β can be determined as follows:
β ^ = 1 n i = 1 n q ^ i k ^ i
where q ^ i and k ^ i are the identified values of qi and ki, respectively.
The flowchart of the proposed approach is shown in Figure 1. As can be seen from Figure 1, the LSE technique is employed for parameter identification and the input modification process is used to force the estimated ground acceleration to comply with its spatial distribution along the structure. The key point of this approach is the input modification procedure as described in step 3, which converts the spatial distribution of ground motion into the additional information for parameters identification. It can also be observed that there are no complicated formulas and their derivations involved, implying the proposed approach is quite simple and can be easily realized. The accuracy of the proposed approach for the identification of structural parameters and unknown ground motion is numerically validated via a shear building model under the El Centro earthquake. More details can be found in the following section.

3. Numerical Investigation

3.1. General Description

A five-story shearing building model as shown in Figure 2 is employed in this study for the investigation of the effectiveness of the proposed iterative identification procedure. The structural parameters are set as mi = 500 kg and ki = 6.0×105 N/m (i = 1, 2, …, 5). The Rayleigh damping coefficients are set as α = 0.6 and β = 8.0 × 10−5 providing approximately 3% damping ratio for the first mode of vibration. To simulate the structural damage, the stiffness of the first floor is assumed to be 80% of its designed value in the healthy state (i.e., k1 = 4.8 × 105 N/m).
In this numerical example, the building structure is assumed to be subject to El-Centro earthquake with a peak ground acceleration (PGA) of 0.34 g. The corresponding structural responses are calculated using the Newmark method with a time interval of 0.02 s. Based on Equation (6), the response matrix H(t) used for the identification in this study can be expressed as
H ( t ) = [ H 1 , 1 H 2 , 1 H 3 , 1 H 1 , k H 2 , k H 3 , k H 1 , s H 2 , s H 3 , s ]
where
H 1 , k = [ x 1 , k x 1 , k x 2 , k x 2 , k x 1 , k x 2 , k x 3 , k x 3 , k x 2 , k x 3 , k x 4 , k x 4 , k x 3 , k x 4 , k x 5 , k x 5 , k x 4 , k ] ; H 2 , k = [ m 1 x ˙ 1 , k m 2 x ˙ 2 , k m 3 x ˙ 3 , k m 4 x ˙ 4 , k m 5 x ˙ 5 , k ] ; H 3 , k = [ x ˙ 1 , k x ˙ 1 , k x ˙ 2 , k x ˙ 2 , k x ˙ 1 , k x ˙ 2 , k x ˙ 3 , k x ˙ 3 , k x ˙ 2 , k x ˙ 3 , k x ˙ 4 , k x ˙ 4 , k x ˙ 3 , k x ˙ 4 , k x ˙ 5 , k x ˙ 5 , k x ˙ 4 , k ] ; x i , k and x ˙ i , k (i = 1, …, 5; k = 1, …, s) denote the displacement and velocity of the i-th floor at the k-th sampling point, respectively; the subscript ‘s’ denotes the number of sample points as defined before.
The convergence threshold ε in Equation (15) is set to be 1%. The effects of signal noise, the number of sample points, and the initial values of the structural parameters are discussed in detail in the following subsections.

3.2. The Effect of Signal Noise

Here, the structural responses are first assumed to be noise-free. In practical situations, the measured signals are inevitably contaminated by noise. In this example, the noise level is set to be 5% and 10%, which means the ratios of the root mean square (RMS) of the noise to the RMS of the corresponding responses are 5% and 10%, respectively. All of the measured time duration from 0 s to 50 s with 2500 sample points is used for identification. The initial guess of the stiffness is set to be 3.0 × 105 N/m for each floor. The damping coefficients α and β are assumed to be 1.0 and 1.5 × 10−4, respectively. The element-wise coefficient qi is therefore initially equal to 45. By using the proposed approach, the identification results are given in Table 1.
It can be observed from Table 1 that the stiffness can be precisely identified in the noise-free case. Although the identification errors increase with increasing noise intensity, it is still clear that the proposed approach is capable of satisfactorily identifying the structural parameters, including stiffness and damping coefficient α. The maximum identification errors for the stiffness in the case of 5% noise and 10% noise are 0.10% and 0.41%, respectively. The errors for identified damping coefficients are relatively large, but even with 10% noise the maximum error is only 4.33%. The comparison of the effect of noise intensity on the iteration process is given in Figure 3. Only the identified stiffness of the 5th floor (k5) is plotted in Figure 3 as an example. Similar results can be found for the remaining parameters. It can be observed from Figure 3 that the larger the noise intensity involved, the more iterations are required. For the cases of noise-free, 5% noise, and 10% noise, the number of iterations is 147, 172, and 248, respectively.
To further show the convergence process during the iteration, the identified stiffness in the case of 10% noise is given in Figure 4 as an example. Obviously, although the iteration is conducted for 248 times due to the relatively stringent tolerance defined before (i.e., ε = 1%), the identified stiffness parameters are stably converged to actual ones, the number of iterations being approximately 50. On the basis of the identified stiffness and the element-wise coefficient qi in matrix Q, the damping coefficient β can be calculated according to Equation (16) as 7.556 × 10−5, 6.737 × 10−5, and 4.103 × 10−5 for the cases of noise-free, 5% noise, and 10% noise, respectively. Since the value of β is extremely small and difficult to assess, these identified results for β, to some extent, are acceptable. Since only the standard LSE technique and the input modification process are involved, the computational time for the proposed approach is very short. In this example, only 10 seconds are required for the identification.
Besides the identification of structural parameters, the unknown ground motion can be estimated by the proposed approach as well. Similarly, taking the case of 10% noise as an example, the comparison of the identified ground acceleration with the actual one is given in Figure 5a. The dashed line represents the identified ground acceleration, whereas the solid line is its actual value. For clarity of comparison, the time segment from 4 s to 6 s is plotted in Figure 5b. Clearly, the proposed approach can identify the ground acceleration very well, even when the signals are contaminated by 10% noise. Based on these identified parameters and ground motion, the structural responses can also be calculated using the equation of motion. Figure 6 gives the comparison of the actual structural responses with the estimated ones. Due to the limitation of the length of this paper, only the responses on the top floor are shown as examples. It is obvious that the estimated responses are in good agreement with the actual ones.

3.3. The Effect of the Number of Sample Points

In this section, the effect of the number of sample points is discussed. Similarly, the initial values of the stiffness and damping coefficients α and β are also assumed to be 3.0 × 105 N/m and 1.0 and 1.5 × 10−4, respectively. The element-wise coefficient qi is thus initially equal to 45. For practical consideration, noise of 10% intensity as mentioned before is introduced to the measured signals. Herein, the time durations of 6 s, 12 s, and 18 s with 300, 600, and 900 sample points, respectively, are considered. The identification results for these three cases are given in Table 2. It can be seen in Table 2 that as the number of sample points increases, the identification error decreases. By comparing the identified results using 600 sample points with those using 900 sample points, the results are improved slightly when even more sample points are being used. This means that the effect of the number of sample points on the identification results becomes small when sufficient sample points are used. To further explore the effect on the convergence procedure, the comparison of the identified k5 for the cases of 300, 600, and 900 sample points is shown in Figure 7 as an example. It can be seen that the identified parameters can be stably converged to the actual ones. Moreover, the iterative process is conducted 415, 312, and 291 times for the cases of 300, 600, and 900 sample points, respectively, which means that more iterations are required if fewer sample points are employed.
Moreover, by using the proposed approach, the ground acceleration can be identified as well. Taking the case of 900 sample points as an example, Figure 8 gives the identified excitation as a dashed line, whereas the actual one is plotted as a solid line. It is clear that the identified ground motion is close to the corresponding actual one.

3.4. The Effect of Initial Values of Structural Parameters

The influence of the initial values of structural parameters on the convergence is investigated in this section. The structural responses are assumed to be contaminated by 10% noise. The time durations from 2 s to 20 s with 900 sample points are used for the identification. The initial values of the structural parameters are set to be 5 times, 10 times, and 15 times larger than the actual ones, respectively. The identification results are given in Table 3. It can be observed that the maximum identification error for the stiffness is only 1.08%, even though the initial guesses are 15 times larger than the actual ones. To further investigate the effect of the initial values on the iteration procedure, the comparison of the identified stiffness with different initial values is given in Figure 9. Similarly, only the identified k5 is plotted in Figure 9 as an example. Obviously, the identified values are stably and accurately converged to the actual ones for the aforementioned three cases. Similar results can be obtained for the remaining parameters.
Similarly, the unknown ground acceleration can be also identified. The comparison results for the case of 15 times initial values are plotted in Figure 10. It is clear that satisfactory identification results can be obtained.

4. Conclusions

In this paper, a relatively simple and efficient iteration identification procedure is proposed for the simultaneous identification of the parameters of building structure and the unknown ground motion. There are no complicated formulas and their derivations involved, which implies the proposed iterative identification procedure can be easily understood and implemented. The standard LSE technique is employed for the estimation of unknown parameters, and an input modification process is developed to ensure that the parameters to be identified would be converged to actual ones. The feasibility and accuracy of the proposed approach is verified via numerical investigation of a shear building model under the El Centro earthquake. The effects of signal noise, the number of sample points, and the initial guesses of structural parameters are discussed. It can be concluded from the numerical results that the proposed approach provides an alternative way for simultaneously identifying the parameters of building structures and the unknown ground motion.
It should be noted that the proposed approach is performed under the condition of the structural responses all being available. Moreover, the effect of the bias caused by a faulty sensor or environmental changes is not considered in this study. Further research can be conducted with the consideration of limited output information and measurement bias.

Author Contributions

J.J. wrote the original draft and participated in paper revision; J.H. presented the framework of the proposed approach and reviewed the manuscript; M.Y. and L.J. wrote the MATLAB main function code; Z.T., Y.L., and H.S. gave the numerical examples for the validation of the proposed approach.

Funding

This research was funded by the Petro China Innovation Foundation (2016D-5007-0608), the Northeast Petroleum University Graduate Innovation Research Project (JYCX_CX05_2018), the China Petroleum and Chemical Industry Federation Science and Technology Guidance Project (2017-11-5), and the Special Foundation of The Education Department of Heilongjiang Province (No.2016YSFX-02). The support from the National Natural Science Foundation of China (No. 51708198) is also greatly appreciated.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Xu, Y.L.; He, J. Smart Civil Structures; CRC Press: Boca Raton, FL, USA, 2017; pp. 3–30. [Google Scholar]
  2. Farrar, C.R.; Worden, K. An Introduction to Structural Health Monitoring. Philos. Trans. R. Soc. A 2007, 365, 303–315. [Google Scholar] [CrossRef] [PubMed]
  3. Ou, J.P.; Li, H. Structural Health Monitoring in Mainland China: Review and future trends. Struct. Control Health Monit. 2010, 9, 219–231. [Google Scholar]
  4. Palacz, M. Spectral Methods for Modelling of Wave Propagation in Structures in Terms of Damage Detection-A Review. Appl. Sci. 2018, 8, 1124. [Google Scholar] [CrossRef]
  5. Ling, Y.; Li, D.; Song, G.; Zhan, C. Damage Detection of Shear Buildings through Structural Mass-stiffness Distribution. Smart Struct. Syst. 2017, 19, 11–20. [Google Scholar] [CrossRef]
  6. Pappalardo, C.M.; Guida, D. System Identification Algorithm for Computing the Modal Parameters of Linear Mechanical Systems. Machines 2018, 6, 12. [Google Scholar] [CrossRef]
  7. Zhang, P.; Huo, L.; Song, G. Impact Fatigue of Viscoelastic Materials Subjected to Pounding. Appl. Sci. 2018, 8, 117. [Google Scholar] [CrossRef]
  8. He, J.; Zhang, X.; Qi, M.; Xu, B. Model-Free Identification of Nonlinear Restoring Force with Modified Observation Equation. Appl. Sci. 2019, 9, 306. [Google Scholar] [CrossRef]
  9. Ji, J.; Xu, Z.C.; Jiang, L.Q. Nonlinear Buckling Analysis of H-type Honeycombed Composite Column with Rectangular Concrete-filled Steel Tube Flanges. Int. J. Steel Struct. 2018, 18, 1153–1166. [Google Scholar] [CrossRef]
  10. Ji, J.; Yu, D.Y.; Jiang, L.Q.; Xu, Z.C. Effect of Post-Fire Curing on the Compressive Properties of Fire-Damaged Ultra-High Toughness Cementitious Composites. J. Test. Eval. 2019, 47, 140–152. [Google Scholar] [CrossRef]
  11. Hamey, C.S.; Lestari, W.; Qiao, P.Z.; Song, G. Experimental Damage Identification of Carbon/Epoxy Composite Beams Using Curvature Mode Shapes. Int. J. Struct. Health Monit. 2004, 3, 333–353. [Google Scholar] [CrossRef]
  12. Di Massa, G.; Russo, R.; Strano, S.; Terzo, M. System structure identification and adaptive control of a seismic isolator test rig. Mech. Syst. Signal Proc. 2013, 40, 736–753. [Google Scholar] [CrossRef]
  13. Shabbir, F.; Khan, M.I.; Ahmad, N.; Tahir, M.F.; Ejaz, N.; Hussain, J. Structural Damage Detection with Different Objective Functions in Noisy Conditions Using an Evolutionary Algorithm. Appl. Sci. 2017, 7, 1245. [Google Scholar] [CrossRef]
  14. Moreno-Gomez, A.; Amezquita-Sanchez, J.P.; Valtierra-Rodriguez, M.; Perez-Ramirez, C.A.; Dominguez-Gonzalez, A.; Chavez-Alegria, O. EMD-Shannon Entropy-Based Methodology to Detect Incipient Damages in a Truss Structure. Appl. Sci. 2018, 8, 2068. [Google Scholar] [CrossRef]
  15. Nagarajaiah, S.; Basu, B. Output Only Modal Identification and Structural Damage Detection Using Time Frequency & Wavelet Techniques. Earthq. Eng. Eng. Vib. 2009, 8, 583–605. [Google Scholar]
  16. Yao, R.; Pakzad, S.N. Time and Frequency Domain Regression-based Stiffness Estimation and Damage Identification. Struct. Control Health Monit. 2014, 21, 356–380. [Google Scholar] [CrossRef]
  17. Perez-Ramirez, C.A.; Jaen-Cuellar, A.Y.; Valtierra-Rodriguez, M.; Dominguez-Gonzalez, A.; Osornio-Rios, R.A.; Romero-Troncoso, R.D.J.; Amezquita-Sanchez, J.P. A Two-Step Strategy for System Identification of Civil Structures for Structural Health Monitoring Using Wavelet Transform and Genetic Algorithms. Appl. Sci. 2017, 7, 111. [Google Scholar] [CrossRef]
  18. Caravani, P.; Watson, M.L.; Thomson, W.T. Recursive Least-squares Time Domain Identification of Structural Parameters. J. Appl. Mech. 1977, 44, 135–140. [Google Scholar] [CrossRef]
  19. Lin, J.W.; Betti, R.; Smyth, A.W.; Longman, R.W. On-line Identification of Non-linear Hysteretic Structural Systems Using a Variable Trace Approach. Earthq. Eng. Struct. Dyn. 2001, 30, 1279–1303. [Google Scholar] [CrossRef]
  20. Masri, S.F.; Caffery, J.P.; Caughey, T.K.; Smyth, A.W.; Chassiakos, A.G. A General Data-based Approach for Developing Reduced-order Models of Non-linear MDOF Systems. Nonlinear Dyn. 2005, 39, 95–112. [Google Scholar] [CrossRef]
  21. Yang, J.N.; Lin, S. Identification of Parametric Variations of Structures Based on Least Square Estimation and Adaptive Tracking Technique. J. Eng. Mech. 2005, 131, 290–298. [Google Scholar] [CrossRef]
  22. Yang, J.N.; Huang, H.W.; Lin, S. Sequential Non-linear Least-square Estimation for Damage Identification of Structures. Int. J. Non-Linear Mech. 2006, 41, 124–140. [Google Scholar] [CrossRef]
  23. Gui, G.; Peng, W.; Adachi, F. Adaptive System Identification Using Robust LMS/F Algorithm. Int. J. Commun. Syst. 2014, 27, 2956–2963. [Google Scholar] [CrossRef]
  24. Kazemi, M.; Arefi, M.M. A Fast Iterative Recursive Least Squares Algorithm for Wiener Model Identification of Highly Nonlinear Systems. ISA Trans. 2017, 67, 382–388. [Google Scholar] [CrossRef]
  25. Wang, D.; Haldar, A. System Identification with Limited Observations and Without Input. J. Eng. Mech. 1997, 123, 504–511. [Google Scholar] [CrossRef]
  26. Ling, X.L.; Haldar, A. Element Level System Identification with Unknown Input with Rayleigh Damping. J. Eng. Mech. 2004, 130, 877–885. [Google Scholar] [CrossRef]
  27. Chen, J.; Li, J. Simultaneous Identification of Structural Parameters and Input Time History from Output-only Measurements. Comput. Mech. 2004, 33, 365–374. [Google Scholar] [CrossRef]
  28. Yang, J.N.; Pan, S.; Lin, S. Least-squares Estimation with Unknown Excitations for Damage Identification of Structures. J. Eng. Mech. 2007, 133, 12–21. [Google Scholar] [CrossRef]
  29. Yang, J.N.; Huang, H.W. Sequential Non-linear Least-square Estimation for Damage Identification of Structures with Unknown Inputs and Unknown Outputs. Int. J. Non-Linear Mech. 2007, 42, 789–801. [Google Scholar] [CrossRef]
  30. Xu, B.; He, J.; Rovekamp, R.; Dyke, S.J. Structural Parameters and Dynamic Loading Identification from Incomplete Measurements: Approach and validation. Mech. Syst. Signal Proc. 2012, 28, 244–257. [Google Scholar] [CrossRef]
  31. Xu, B.; He, J. Substructural parameters and dynamic loading identification with limited observations. Smart Struct. Syst. 2015, 15, 169–189. [Google Scholar] [CrossRef]
  32. Xu, B.; He, J.; Masri, S.F. Data-based Model-free Hysteretic Restoring Force and Mass Identification for Dynamic Systems. Comput.-Aided Civil Infrastruct. Eng. 2015, 30, 2–18. [Google Scholar] [CrossRef]
  33. Lei, Y.; Jiang, Y.; Xu, Z. Structural Damage Detection with Limited Input and Output Measurement Signals. Mech. Syst. Signal Proc. 2012, 28, 229–243. [Google Scholar] [CrossRef]
  34. Lei, Y.; Wang, Q.; Liu, L. Structural Parameter Identification of Wind-excited High-rise Buildings with Limited Output Measurements. Appl. Math. Inform. Sci. 2016, 10, 235–241. [Google Scholar] [CrossRef]
  35. Askari, M.; Li, J.; Samali, B. A Compact Self-adaptive Recursive Least Square Approach for Real-time Structural Identification with Unknown Inputs. Adv. Struct. Eng. 2016, 19, 1118–1129. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed iterative identification procedure.
Figure 1. Flowchart of the proposed iterative identification procedure.
Applsci 09 00696 g001
Figure 2. Five-story building model in the numerical example.
Figure 2. Five-story building model in the numerical example.
Applsci 09 00696 g002
Figure 3. Comparison of the effect of noise intensity on the iteration process: (a) noise-free; (b) 5% noise; (c) 10% noise.
Figure 3. Comparison of the effect of noise intensity on the iteration process: (a) noise-free; (b) 5% noise; (c) 10% noise.
Applsci 09 00696 g003
Figure 4. The identified structural parameters (10% noise): (a) k1; (b) k2; (c) k3; (d) k4; (e) k5
Figure 4. The identified structural parameters (10% noise): (a) k1; (b) k2; (c) k3; (d) k4; (e) k5
Applsci 09 00696 g004
Figure 5. Identified results of the ground acceleration (10% noise): (a) the whole time history; (b) the time segment from 4 s to 6 s.
Figure 5. Identified results of the ground acceleration (10% noise): (a) the whole time history; (b) the time segment from 4 s to 6 s.
Applsci 09 00696 g005
Figure 6. Comparison of the structural responses on the top floor (10% noise): (a) acceleration; (b) velocity; (c) displacement.
Figure 6. Comparison of the structural responses on the top floor (10% noise): (a) acceleration; (b) velocity; (c) displacement.
Applsci 09 00696 g006
Figure 7. Comparison of the identified k5 using different numbers of sample points: (a) using 300 sample points; (b) using 600 sample points; (c) using 900 sample points.
Figure 7. Comparison of the identified k5 using different numbers of sample points: (a) using 300 sample points; (b) using 600 sample points; (c) using 900 sample points.
Applsci 09 00696 g007
Figure 8. Identified ground acceleration using 900 sample points.
Figure 8. Identified ground acceleration using 900 sample points.
Applsci 09 00696 g008
Figure 9. Comparison of the identified k5 with different initial values: (a) 5 times; (b) 10 times; (c) 15 times.
Figure 9. Comparison of the identified k5 with different initial values: (a) 5 times; (b) 10 times; (c) 15 times.
Applsci 09 00696 g009
Figure 10. Identified results of the ground acceleration (15 times).
Figure 10. Identified results of the ground acceleration (15 times).
Applsci 09 00696 g010
Table 1. Identified structural parameters under different noise levels.
Table 1. Identified structural parameters under different noise levels.
Structural ParametersNoise-Free5% Noise10% Noise
IdentifiedError (%)IdentifiedError (%)IdentifiedError (%)
k14.800 × 1050.004.795 × 1050.104.781 × 1050.41
k26.000 × 1050.005.995 × 1050.085.979 × 1050.35
k36.000 × 1050.005.994 × 1050.105.979 × 1050.35
k46.000 × 1050.005.995 × 1050.085.981 × 1050.32
k56.000 × 1050.005.996 × 1050.075.983 × 1050.28
α0.6152.500.6172.830.6264.33
Note: The error of 0.00% in this table means it is less than 10−4%.
Table 2. Identified structural parameters with different numbers of sample points.
Table 2. Identified structural parameters with different numbers of sample points.
Structural Parameters300 Sample Points600 Sample Points900 Sample Points
IdentifiedError (%)IdentifiedError (%)IdentifiedError (%)
k14.715 × 1051.774.756 × 1050.924.768 × 1050.67
k25.907 × 1051.555.952 × 1050.805.964 × 1050.60
k35.912 × 1051.475.953 × 1050.785.966 × 1050.57
k45.914 × 1051.435.954 × 1050.775.967 × 1050.55
k55.932 × 1051.135.963 × 1050.625.974 × 1050.43
α0.53411.000.5587.000.5714.83
Table 3. The identified parameters with different initial values.
Table 3. The identified parameters with different initial values.
Structural Parameters5 Times10 Times15 Times
IdentifiedError (%)IdentifiedError (%)IdentifiedError (%)
k14.758 × 1050.884.751 × 1051.024.748 × 1051.08
k25.954 × 1050.665.948 × 1050.875.944 × 1050.93
k35.956 × 1050.735.950 × 1050.835.946 × 1050.90
k45.957 × 1050.715.951 × 1050.815.947 × 1050.88
k55.964 × 1050.605.960 × 1050.675.954 × 1050.77
α0.6386.330.6427.000.6488.00

Share and Cite

MDPI and ACS Style

Ji, J.; Yang, M.; Jiang, L.; He, J.; Teng, Z.; Liu, Y.; Song, H. Output-Only Parameters Identification of Earthquake-Excited Building Structures with Least Squares and Input Modification Process. Appl. Sci. 2019, 9, 696. https://doi.org/10.3390/app9040696

AMA Style

Ji J, Yang M, Jiang L, He J, Teng Z, Liu Y, Song H. Output-Only Parameters Identification of Earthquake-Excited Building Structures with Least Squares and Input Modification Process. Applied Sciences. 2019; 9(4):696. https://doi.org/10.3390/app9040696

Chicago/Turabian Style

Ji, Jing, Maomao Yang, Liangqin Jiang, Jia He, Zhenchao Teng, Yingchun Liu, and Huayu Song. 2019. "Output-Only Parameters Identification of Earthquake-Excited Building Structures with Least Squares and Input Modification Process" Applied Sciences 9, no. 4: 696. https://doi.org/10.3390/app9040696

APA Style

Ji, J., Yang, M., Jiang, L., He, J., Teng, Z., Liu, Y., & Song, H. (2019). Output-Only Parameters Identification of Earthquake-Excited Building Structures with Least Squares and Input Modification Process. Applied Sciences, 9(4), 696. https://doi.org/10.3390/app9040696

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