Model-Free Identification of Nonlinear Restoring Force with Modified Observation Equation

Nonlinearity exists widely in civil engineering structures; for example, the initiation and growth of damage under dynamic loadings is a typical nonlinear process. To date, for the purpose of structural evaluation and a better understanding nonlinear characteristics of complicated structures, a number of parametric and nonparametric methods have been developed for the identification of nonlinear restoring force (NRF). However, due to the highly individualistic nature of nonlinear systems, it would be inefficient to attempt to express the structural NRF in a general parametric form. For many nonparametric techniques, their nonparametric models or approximations may result in undesirable results or oscillations around unsmooth points. In this paper, on the basis of extended Kalman filter (EKF), a model-free NRF identification approach is proposed to circumvent the limitations mentioned above. The NRF to be identified was treated as ‘unknown fictitious input’, and thus, no prior assumptions or approximations for the NRF model were required. With the aid of a projection matrix, a modified version of observation equation was obtained. Based on the principle of EKF, the recursive solution of the proposed approach was analytically derived. The NRFs provided by the nonlinear components were identified by means of least squares estimation (LSE) at each time step. Numerical examples, including building structures equipped with magnetorheological (MR) damper and shape memory alloy (SMA) damper, demonstrated that the proposed approach is capable of satisfactorily identifying NRF without knowledge or intuitive assumptions of any nonlinear model class in advance.


Background
Civil engineering structures are prone to deterioration or damage due to certain factors, such as strong dynamic loadings, material fatigue accumulation and environmental corrosion.To date, a variety of vibration-based damage detection and system identification methods have been developed [1,2].Most of them are applicable when the structure of interest is linear.However, nonlinearity is generic in nature, and linear behavior is an exception [3].The occurrence of a fault in an initially linear structure will, in many cases, result in nonlinear behavior.Therefore, it is highly desirable to develop nonlinear identification algorithms to facilitate the evaluation of the structural performance and condition for the purpose of prevention of potentially catastrophic events.
It is known that typical sources of nonlinearities include geometric nonlinearity, inertia nonlinearity, material nonlinearity and boundary nonlinearity.It is of importance to have an accurate model for the description of the nonlinear behavior of the physical structure prior to parameter estimation.Nonlinear modeling is a challenging step, because nonlinearity may be caused by many different mechanisms and may result in plethora of dynamic phenomena.Without a precise understanding of the nonlinear mechanisms involved, the identification process is bound to fail.A priori knowledge and nonlinearity classification may help to select a reasonably accurate model of the nonlinearity.If little is known about the form of the model before starting the identification process, one may resort to nonparametric approximations or expressions.

Literature Survey
Since the restoring force can act as an indicator of the extent of the nonlinearity and be used to quantitatively evaluate the absorbed energy during vibration, the research on nonlinear restoring force (NRF) identification has been actively conducted over years and the comprehensive literature reviews can be found [3,4].Broadly speaking, the existing NRF identification methods can be divided into two groups, i.e., parametric identification and nonparametric identification.In parametric identification methods, a model structure is specified, whereas the primary quantities obtained from the nonparametric identification methods do not directly specify equations of motion [5].Once the structural model is fixed, the nonlinear identification problem is reduced to parameter estimation as only the coefficients of the model terms remain unspecified.The subject of parametric nonlinear identification is extremely broad, and an extensive literature exists (e.g., [6][7][8][9][10]).It is not intended to provide a comprehensive overview of the past and current approaches for the parametric identification of nonlinearity in this paper.Instead, some relevant research regarding extended Kalman filter (EKF)-based nonlinear parameter identification methods is introduced herein.For example, Yang et al. [11] proposed an EKF-based adaptive tracking technique for the identification of time-vary parameters of linear and nonlinear hysteretic structures.The performance of EKF and Unscented Kalman filter (UKF) for the nonlinear identification was discussed in details by Wu and Smyth [12].Based on the sequential application of extended Kalman estimator and least squares estimation (LSE), Lei et al. [13] proposed an algorithm for the identification of nonlinear structural parameters.On the basis of the orthogonal decomposition of excitation, Ding et al. [14] proposed an EKF-based method for the condition assessment of linear and nonlinear hysteretic structures.By using EKF to estimate time-invariant parameters of the nonlinear material constitutive models, a novel nonlinear finite element model updating framework was proposed by Ebrahimian et al. [15].A two-stage NRF identification method was proposed by Lei et al. [16] by using EKF for the determination of nonlinear locations in the first stage and the estimation of NRFs in the second stage.This method was extended by Liu et al. [17] for identifying strong nonlinearity parameters.With the combination of EKF and dynamic statistical process control, Jin et al. [18] presented a novel real-time damage detection method for identifying the parameters of both linear and nonlinear structures under different damage scenarios.Xiao et al. [19] proposed an adaptive three-stage EKF method to solve state and fault estimation in nonlinear discrete-time system.A dual adaptive filtering approach was proposed by Astroza [20] for nonlinear finite element model updating with the consideration of modeling uncertainty.The performance of different model updating techniques using EKF, UKF and iterated EKF were further compared by Astroza et al. [21] in terms of convergence, accuracy, robustness and computational demand.
The methods mentioned above are all capable of identifying the NRF of a preset nonlinear model.However, accessing prior knowledge of nonlinear characteristics may prove difficult in many circumstances owing to the highly individualistic nature of real nonlinearities.Moreover, it is sometimes difficult to define a specific physically-motivated model to precisely describe the encountered nonlinear phenomena such as the initiation and development of crack.The lack of high-fidelity, physics-based and robust models that could accurately reflect the multidimensional behavior of these nonlinear systems would hinder the application of the computational tools or methods in identifying the complex nonlinearities in the design stage.Thus, nonparametric nonlinear identification solely using input and output data is necessary.Increasing attention has been received in this area and much progress has been made.Early contributions to nonparametric nonlinear identification can be traced back to the 1980s when Masri et al. developed their famous restoring force surface (RFS) methods by using Chebyshev polynomial [22][23][24][25].The data processing and experimental design for the investigation of RFS method were further conducted by Worden [26,27].A two-stage general procedure was presented by Masri et al. [28] for developing data-based, model-free representations of complex nonlinear systems under arbitrary dynamic loads.By using a power series polynomial to represent the NRF, a data-based approach was proposed for the identification of chain-like structure equipped with magneto-rheological (MR) damper [29,30] and shape memory alloy (SMA) damper [31].The generalized two-variable Chebyshev polynomials and their relevant relations were further discussed by Cesarano and Fornaro [32,33].A data-based nonlinear identification approach was proposed by Xu et al. [34] for mass and NRF identification under incomplete loadings.By introducing an equivalent linear structural system and a power series polynomial, Lei et al. [35] proposed an EKF-based model-free NRF identification approach for the identification of the unknown coefficients of the power series polynomials and other structural parameters.Su et al. [36] proposed a two-step approach to identify the nonlinearities of model-free MR damper with partial response data.Based on the Hilbert transform, two novel nonparametric identification approaches were proposed by Yuan et al. [37] for nonlinear piezoelectric mechanical systems.Ondra et al. [38] proposed a zero-crossing method for nonparametric identification of systems with asymmetric NRFs.As the rapid development of computers and their associated computational capabilities, the intelligent black-box methods for nonparametric nonlinear identification using artificial neural networks, fuzzy networks, genetic algorithms and so forth have received increasing attention.The literature review on these natural computing methods with different emphasis can be found [4,39].

Scope
Most of the aforementioned nonparametric identification methods are based on the usage of some specific nonparametric expressions to describe or approximate the nonlinear properties.The key of these methods is to properly determine the values of the coefficients in their nonparametric expressions basing on the input and output data.Such approximations may, however, result in undesirable oscillations or inaccurate results around unsmooth points.In this paper, an EKF-based model-free nonlinear identification approach is proposed for the identification of NRFs and other structural parameters.By using a projection matrix, a modified version of observation equation is obtained.The NRF is treated as 'unknown fictitious input' and identified by means of LSE.Based on the principle of EKF, the recursive solution of the proposed approach is analytically derived.The proposed approach requires input and output data only, without knowledge or intuitive assumptions of any nonlinear model class in advance.
The present research work is organized as follows: Section 2: formulation of the proposed model-free NRF identification procedure; Section 3: validation of the proposed approach via a building structure equipped with MR damper and SMA damper; Section 4: conclusions and discussions of the proposed approach.

Formulas of the Proposed Approach
The general equation of motion of a discrete n-DOFs nonlinear structure under external excitations can be expressed as: in which M is the mass matrix and assumed to be known; C and K are the damping matrix and stiffness matrix, respectively; Appl.Sci.2019, 9, 306 4 of 13 Here, the state vector of the structure is defined as Z(t) = x(t) T , .
x(t) T , θ T T , where θ is the unknown time-invariant structural parameters (i.e., stiffness and damping coefficients).Then, the following state equation can be found: .
where w(t) is model noise vector with zero mean and a covariance matrix Q(t).
It should be noted, besides the unknown structural parameters, the NRF R non k is assumed to be unavailable as well.Since the NRF is treated as 'unknown fictitious inputs' herein, it can be found the prior knowledge or intuitive assumptions of any nonlinear model class of NRF is not required.Let Ẑk|k and Rnon k be the estimations of Z k and R non k , respectively.Using Taylor expansion, Equation (2) can be linearized at time instant t = k × ∆t with the sampling interval of ∆t: where the subscript k denotes the k-th time step: In consideration of .Z(k∆t) = (Z k+1 − Z k )/∆t, one can obtain the following expression with the aid of Equations ( 2) and (3): Generally, the acceleration responses are more easily obtained and reliable as compared with other responses, such as displacement and velocity.Thus, the acceleration is assumed to be measured in this study.The discretized observation equation can be then expressed as: where y k is the measured acceleration vector; the subscript k denotes time step t = k × ∆t; v k is the measurement noise vector assumed to be a Gaussian white noise vector with zero mean and a covariance matrix _ R k Equation ( 6) can be also re-written as: where The re-arrangement of Equation ( 7) leads to the following expression: The closest solution of R non k can be then obtained using LSE as: The error of the aforementioned solution can be calculated as: where I is an identity matrix; D D T D −1 D T is known as the projection matrix.As a limit, the error shown in Equation ( 10) should tend to be zero, leading to: As compared with Equation ( 6), it can be found that multiple unknown issues have been transformed into a single unknown problem.Then, the EKF technique can be used to deal with such single regression problem.
Analogous to the traditional EKF, the priori state estimate Ẑk+1|k can be given as: Then, with the combination of Equation ( 5), the priori estimate error ε k+1|k can be found as: Focusing on the second term on the right-hand-side of Equation ( 13), the following transformation can be obtained: Similarly, h(Z k ) can be linearized as follows: where: Substituting Equations ( 14) and (15) into Equation ( 13) leads to: where: The priori estimate error covariance matrix can be then determined as: Based on Equation (11) and the principle of KF, the posteriori state estimate Ẑk+1|k+1 herein can be calculated as a linear combination of the priori state estimate Ẑk+1|k and a weighted difference between the actual measurements and the corresponding predictions: where G k+1 is the gain matrix at time step t = (k + 1) × ∆t.Based on Equations ( 5), ( 15) and ( 20), the posteriori estimate error ε k+1|k+1 can be calculated as: Appl.Sci.2019, 9, 306 where H k+1|k can be determined using Equation ( 16) when Z k = Ẑk+1|k .Accordingly, the posteriori estimate error covariance matrix can be obtained as: Based on the principle of EKF, the gain matrix can be determined by the following equation: where tr(•) denotes the trace operator.Then, one can obtain: Once the posteriori state estimate Ẑk+1|k+1 is determined through Equation ( 20), the unknown NRF can be determined by Equation ( 9) as follows: It is clear that the recursive solution for the proposed approach includes priori estimation through Equations ( 12) and ( 19), posteriori estimation through Equations ( 20) and ( 22), gain matrix calculation by Equation ( 24), as well as NRF identification by Equation (25).The effectiveness and robustness of the proposed method will be verified via several numerical examples in the following section.

Numerical Application
Here, a six-story shear building structure equipped with different types of nonlinear hysteretic dampers was used to demonstrate the effectiveness of the proposed approach.The parameters of the building structure were assumed to be m i = 300 kg and k i = 1.8 × 10 5 N/m (i = 1, 2, . . ., 6).The Rayleigh damping assumption was adopted with the proportional coefficients for the mass matrix and stiffness matrix being α = 0.2644 and β = 2.578 × 10 −3 , respectively.For all of the examples, the covariance matrices of measurement noise and model noise were chosen to be _ R k = 10 0 × I and Q k = 10 −3 × I, where I is the identity matrix with appropriate dimension.All the measured signals will be superimposed with the corresponding noise with 5% noise-to-signal ratio in root mean square (RMS).

Building Structure Equipped with MR Damper
The building structure equipped with MR damper on the first floor, as shown in Figure 1, was considered as the first example.The modified Dahl model [40], as expressed by Equations ( 26) and ( 27), was employed for the description of NRF provided by MR damper: where k d , c d , f d and f 0 are the coefficients of the model; s and .
s are the relative displacement and relative velocity of the MR damper, respectively.Since the MR damper is located on the first floor, s and .s are equal to the relative displacement and relative velocity of the first floor, respectively.z i is a dimensionless hysteretic coefficient using for the description of Coulomb friction: .
in which σ is the coefficient controlling the shape of the hysteretic loop; sgn(•) denotes signum function.
fictitious input'.No assumption or prior knowledge of the NRF was made.The unknown quantities to be identified include the extended state vector   ,..., , ,..., , ,..., , , and the NRF provided by MR damper.The acceleration responses were measured for the identification.The initial values of the unknown parameters were all chosen to be 50% of the corresponding actual ones.Based on the proposed approach, the unknown structural parameters were identified as shown in Table 1.It can be seen from Table 1 that the maximum error was only 0.42% indicating the proposed approach can accurately identify the structural parameters.Figure 2 gives the convergence process of the parameters as dash lines.Due to the limitation of paper length, only the identified stiffness coefficients of the first and third floor are shown in Figure 2 as examples.Similar results can be found for the remaining parameters.It is clear that the identified parameters can accurately and stably converge to the actual ones.Moreover, the unknown NRF provided by MR damper can be identified at the same time as shown in Figure 3. Obviously, a good agreement between the identified results and the actual ones can be found.The El-Centro seismic excitation with peak ground acceleration (PGA) of 0.34 g was applied to the building structure.The nonlinear structural responses were calculated by using the Runge-Kutta method with a time interval of 0.001 s.It should be noted that after the calculation of structural responses, the NRF provided by MR damper was viewed as 'unknown fictitious input'.
No assumption or prior knowledge of the NRF was made.The unknown quantities to be identified include the extended state vector Z = x 1 , . . ., x 6 , .
x 6 , k 1 , . . ., k 6 , α, β T and the NRF provided by MR damper.The acceleration responses were measured for the identification.The initial values of the unknown parameters were all chosen to be 50% of the corresponding actual ones.
Based on the proposed approach, the unknown structural parameters were identified as shown in Table 1.It can be seen from Table 1 that the maximum error was only 0.42% indicating the proposed approach can accurately identify the structural parameters.Figure 2 gives the convergence process of the parameters as dash lines.Due to the limitation of paper length, only the identified stiffness coefficients of the first and third floor are shown in Figure 2 as examples.Similar results can be found for the remaining parameters.It is clear that the identified parameters can accurately and stably converge to the actual ones.Moreover, the unknown NRF provided by MR damper can be identified at the same time as shown in Figure 3. Obviously, a good agreement between the identified results and the actual ones can be found.

Building Structure Equipped with SMA Damper
In this example, an SMA damper rather than MR damper is employed to investigate the performance of the proposed approach for the identification of NRF around the unsmooth points.The NRF provided by SMA damper in a double-flag form as shown in Figure 4 can be given by the following equation:

Building Structure Equipped with SMA Damper
In this example, an SMA damper rather than MR damper is employed to investigate the performance of the proposed approach for the identification of NRF around the unsmooth points.The NRF provided by SMA damper in a double-flag form as shown in Figure 4 can be given by the following equation: where SMA k are the stiffness provided by SMA in different stages as shown in Figure 4; sgn(•) is the sign function; s is relative displacement of the damper as mentioned before.The definition of sa, sb, sc and sd can be found in Figure 4.It can be seen that the stiffness as well as NRF provided by SMA damper at these points has a sudden change.

Building Structure Equipped with SMA Damper
In this example, an SMA damper rather than MR damper is employed to investigate the performance of the proposed approach for the identification of NRF around the unsmooth points.The NRF provided by SMA damper in a double-flag form as shown in Figure 4 can be given by the following equation: where k 1 SMA and k 2 SMA are the stiffness provided by SMA in different stages as shown in Figure 4; sgn(•) is the sign function; s is relative displacement of the damper as mentioned before.The definition of s a , s b , s c and s d can be found in Figure 4.It can be seen that the stiffness as well as NRF provided by SMA damper at these points has a sudden change.
Here, the SMA damper was assumed to be located on the second floor of the aforementioned six-story building, and a random excitation was applied to the top floor of the structure.The following parameters were selected for SMA damper:   Here, the SMA damper was assumed to be located on the second floor of the aforementioned six-story building, and a random excitation was applied to the top floor of the structure.The following parameters were selected for SMA damper: ,..., , ,..., , ,..., , , and the NRF provided by SMA damper.The initial values of the unknowns were set to be 50% of the corresponding actual ones.By using the proposed approach, the unknown structural parameters could be identified as shown in Table 2. Obviously, the maximum error of the identified parameters was very small.Similarly, only the stiffness of the second and fifth floor is given in Figure 5 as examples.It can be clear that the prompt and stable convergence process can be obtained.By using the proposed approach, the unknown structural parameters could be identified as shown in Table 2. Obviously, the maximum error of the identified parameters was very small.Similarly, only the stiffness of the second and fifth floor is given in Figure 5 as examples.It can be clear that the prompt and stable convergence process can be obtained.Besides the identification of structural parameters including stiffness and damping coefficients, the NRF provided by SMA damper can be simultaneously identified as well.Moreover, to demonstrate the effectiveness of the proposed approach around the unsmooth points, another nonparametric NRF identification method basing on power series polynomial model (PSPM) is employed [24,25].The identified results by means of the proposed approach and PSPM method are plotted in Figure 6.It can be found from Figure 6a that more accurate results can be obtained by using the proposed approach especially around the unsmooth points.Figure 6b gives the comparison of time histories of the identified error by using the proposed approach and PSPM method.Because the initial values of the structural parameters are 50% of the actual ones, the identified NRF error of the proposed approach is larger than that of PSPM method at the beginning.As the parameters are being identified promptly and accurately, the identified error of the proposed approach becomes much smaller than the error obtained from PSPM method.The comparison results show that the proposed approach is capable of satisfactorily identifying the NRF around the unsmooth points without any assumption of the nonlinear model in advance.Besides the identification of structural parameters including stiffness and damping coefficients, the NRF provided by SMA damper can be simultaneously identified as well.Moreover, to demonstrate the effectiveness of the proposed approach around the unsmooth points, another nonparametric NRF identification method basing on power series polynomial model (PSPM) is employed [24,25].The identified results by means of the proposed approach and PSPM method are plotted in Figure 6.It can be found from Figure 6a that more accurate results can be obtained by using the proposed approach especially around the unsmooth points.Figure 6b gives the comparison of time histories of the identified error by using the proposed approach and PSPM method.Because the initial values of the structural parameters are 50% of the actual ones, the identified NRF error of the proposed approach is larger than that of PSPM method at the beginning.As the parameters are being identified promptly and accurately, the identified error of the proposed approach becomes much smaller than the error obtained from PSPM method.The comparison results show that the proposed approach is capable of satisfactorily identifying the NRF around the unsmooth points without any assumption of the nonlinear model in advance.
Besides the identification of structural parameters including stiffness and damping coefficients, the NRF provided by SMA damper can be simultaneously identified as well.Moreover, to demonstrate the effectiveness of the proposed approach around the unsmooth points, another nonparametric NRF identification method basing on power series polynomial model (PSPM) is employed [24,25].The identified results by means of the proposed approach and PSPM method are plotted in Figure 6.It can be found from Figure 6a that more accurate results can be obtained by using the proposed approach especially around the unsmooth points.Figure 6b gives the comparison of time histories of the identified error by using the proposed approach and PSPM method.Because the initial values of the structural parameters are 50% of the actual ones, the identified NRF error of the proposed approach is larger than that of PSPM method at the beginning.As the parameters are being identified promptly and accurately, the identified error of the proposed approach becomes much smaller than the error obtained from PSPM method.The comparison results show that the proposed approach is capable of satisfactorily identifying the NRF around the unsmooth points without any assumption of the nonlinear model in advance.

Building Structure Equipped with SMA and MR Damper
To investigate the performance of the proposed approach for the identification of the structure with multi-type nonlinear characteristics, the six-story building structure with SMA damper and MR damper was considered.All of the parameters for the six-story structure, SMA damper and MR damper were the same as the previous two examples.The SMA damper and MR damper were installed on the first and fourth floor of the building structure, respectively.A random excitation was applied to the top floor, and the structural responses were calculated by the Runge-Kutta

Building Structure Equipped with SMA and MR Damper
To investigate the performance of the proposed approach for the identification of the structure with multi-type nonlinear characteristics, the six-story building structure with SMA damper and MR damper was considered.All of the parameters for the six-story structure, SMA damper and MR damper were the same as the previous two examples.The SMA damper and MR damper were installed on the first and fourth floor of the building structure, respectively.A random excitation was applied to the top floor, and the structural responses were calculated by the Runge-Kutta method with the time interval of 0.001 s.In this example, the unknowns include Z = x 1 , . . ., x 6 , .
x 6 , k 1 , . . ., k 6 , α, β T and the NRFs provided by SMA damper and MR damper.The initial values were chosen to be 50% of the corresponding actual ones.
Based on the proposed approach, the unknown parameters could be identified as shown in Table 3.Clearly, the relative errors were quite small.Figure 7 gives the identified k 1 and k 4 as examples.It can be seen that the identified results can converge to the actual ones fast.Moreover, the NRFs provided by SMA damper and MR damper are identified and plotted in Figure 8 as dashed lines, whereas the theoretical ones are shown as solid lines.It can be observed the identified NRFs are all close to the actual ones.Based on the proposed approach, the unknown parameters could be identified as shown in Table 3.Clearly, the relative errors were quite small.Figure 7 gives the identified k1 and k4 as examples.It can be seen that the identified results can converge to the actual ones fast.Moreover, the NRFs provided by SMA damper and MR damper are identified and plotted in Figure 8   Based on the proposed approach, the unknown parameters could be identified as shown in Table 3.Clearly, the relative errors were quite small.Figure 7 gives the identified k1 and k4 as examples.It can be seen that the identified results can converge to the actual ones fast.Moreover, the NRFs provided by SMA damper and MR damper are identified and plotted in Figure 8

Conclusions
In this paper, a model-free NRF identification approach was proposed basing on the principle of EKF.The NRFs provided by the nonlinear components such as the MR damper and SMA damper were treated as 'unknown fictitious input', and thus, no prior assumptions or approximations for the NRF model were required.A modified version of observation equation was obtained on the basis of the usage of a projection matrix.The recursive solution of the proposed approach was analytically derived.The structural parameters, including the stiffness and damping coefficients, were identified at each time step.The NRFs were also identified by means of LSE at the same time.The effectiveness and robustness of the proposed approach was numerically validated via a six-story building structure equipped with MR damper and/or SMA damper.A further comparison of the proposed approach with the PSPM method was made to demonstrate the reliability of the proposed approach for the identification of NRFs around the unsmooth points.The results show that the proposed approach is capable of satisfactorily identifying the structural parameters and unknown NRFs provided by the nonlinear components.
Although no prior knowledge of NRF model is required, the proposed approach should be performed under the condition that the influence matrix η u associated with NRFs is available.Nonlinearities may often be located easily, at least for simple structures [3].However, for complex structural systems, it is difficult to precisely determine the spatial localization of nonlinearities.Therefore, further research can be conducted towards the model-free NRF identification without the assumption of influence matrix being available.

Figure 1 .
Figure 1.Building structure equipped with magnetorheological (MR) damper.The parameters of the aforementioned model were set to be k d = 25 N/m, c d = 2000 N•s/m, f d = 50 N and σ = 1.0 × 10 3 s/m.The El-Centro seismic excitation with peak ground acceleration (PGA) of 0.34 g was applied to the building structure.The nonlinear structural responses were calculated by using the Runge-Kutta method with a time interval of 0.001 s.It should be noted that after the calculation of structural responses, the NRF provided by MR damper was viewed as 'unknown fictitious input'.No assumption or prior knowledge of the NRF was made.The unknown quantities to be identified
provided by SMA in different stages as shown in Figure4; sgn(•) is the sign function; s is relative displacement of the damper as mentioned before.The definition of sa, sb, sc and sd can be found in Figure4.It can be seen that the stiffness as well as NRF provided by SMA damper at these points has a sudden change.

Figure 3 .
Figure 3.Comparison of the identified MR damper force: (a) Hysteretic loop; (b) Time history of nonlinear restoring force (NRF).
10 −3 m and s b = 8 × 10 −3 m.The acceleration responses were measured with the sampling frequency of 1000 Hz for the identification.Here, the quantities to be identified included the extended state vector Z = x 1 , . . ., x 6 , .x 1 , . . ., .x 6 , k 1 , . . ., k 6 , α, β T and the NRF provided by SMA damper.The initial values of the unknowns were set to be 50% of the corresponding actual ones.
m and sb = 8 × 10 −3 m.The acceleration responses were measured with the sampling frequency of 1000 Hz for the identification.Here, the quantities to be identified included the extended state vector

Figure 6 .
Figure 6.Comparison of SMA damper force: (a) hysteretic loop; (b) time history of the identified error.

Figure 6 .
Figure 6.Comparison of SMA damper force: (a) hysteretic loop; (b) time history of the identified error.

Figure 7 .
Figure 7.The identified stiffness (building with SMA damper and MR damper): (a) first floor; (b) fourth floor.

Figure 8 .
Figure 8.Comparison of NRFs: (a) hysteretic loop of SMA damper; (b) time history of NRF by SMA damper; (c) hysteretic loop of MR damper; (d) time history of NRF by MR damper.

Table 1 .
The identified results of building structure equipped with MR damper.

Table 1 .
The identified results of building structure equipped with MR damper.

Table
The identified results of building structure equipped with SMA damper.

Table 2 .
The identified results of building structure equipped with SMA damper.

Table 3 .
The identified results of structural parameters equipped with SMA and MR damper.