REAL-TIME IDENTIFICATION METHOD FOR NONLINEAR CHARACTERISTIC OF MODEL-FREE IN-SERVICE SEISMIC ISOLATION/VIBRATION REDUCTION DEVICE

Abstract

The present disclosure provides a real-time identification method for a nonlinear characteristic of a model-free in-service seismic isolation/vibration reduction device. The method includes: dividing a structure of a seismic isolation/vibration reduction system into a plurality of substructures, and defining a substructure where the seismic isolation/vibration reduction device is located as a target substructure; and using a General Extended Kalman filter with unknown inputs (GEKF-UI) to identify a linear stiffness and a damping coefficient of the seismic isolation/vibration reduction device and a linear stiffness and a damping coefficient of the target substructure, in the case that an external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state.

Claims

1. A real-time identification method for a nonlinear characteristic of a model-free in-service seismic isolation/vibration reduction device driven based on monitoring data, comprising: dividing a structure of a seismic isolation/vibration reduction system into a plurality of substructures, and defining a substructure where the seismic isolation/vibration reduction device is located as a target substructure; and using a General Extended Kalman filter with unknown inputs (GEKF-UI) to identify a linear stiffness and a damping coefficient of the seismic isolation/vibration reduction device and a linear stiffness and a damping coefficient of the target substructure, in the case that an external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state.

2. The method according to claim 1, wherein the method further comprises: using a General Kalman filter with unknown inputs (GKF-UI) to identify a nonlinear restoring force received by the target substructure by the linear stiffness and the damping coefficient of the seismic isolation/vibration reduction device and the linear stiffness and the damping coefficient of the target substructure, in the case that the external excitation on the seismic isolation/vibration reduction device is large so that the seismic isolation/vibration reduction device generates a nonlinear force.

3. The method according to claim 1, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+(C.sup.r+C){dot over (x)}.sup.r(t)+(K.sup.r+K)x.sup.r(t)=.sup.rf.sup.r(t)+.sup.gf.sup.g(t)=.sup.uf.sup.u(t) wherein, M.sup.r, C.sup.r and K.sup.r are mass, damping and stiffness matrices of the target substructure, respectively; {umlaut over (x)}.sup.r(t), {dot over (x)}.sup.r(t) and x.sup.r(t) are acceleration, velocity and displacement vectors of the target substructure, respectively; f.sup.r(t) is the external excitation received by the target substructure; f.sup.g(t) is a force of an adjacent substructure on the target substructure; .sup.r and .sup.g are position matrices of the external excitation and the force; C and K are additional damping and stiffness provided by the seismic isolation/vibration reduction device to the system structure; .sup.u is a position matrix of the external excitation f.sup.u(t).

4. The method according to claim 3, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state, an extended vector Z=[x.sup.T {dot over (x)}.sup.T .sup.T].sup.T is established, and the GEKF-UI is used to obtain to express the linear stiffness and the damping coefficient of the seismic isolation/vibration reduction device and the linear stiffness and the damping coefficient of the target substructure.

5. The method according to claim 2, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is large so that the seismic isolation/vibration reduction device generates a nonlinear force, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+C.sup.r{dot over (x)}.sup.r(t)+K.sup.rx.sup.r(t)=.sup.rf.sup.r(t)+.sup.gg.sup.r(t).sup.nonf.sup.non(t)=.sup.uf.sup.u(t) wherein C.sup.r and K.sup.r are identified values in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state; the GKF-UI is used to identify f.sup.u(t) by observing a response of a partial structure; the nonlinear force generated by the vibration reduction/seismic isolation device is identified based on the identified f.sup.u(t) and f.sup.u=[f.sup.r(t) g.sup.r(t) f.sup.non(t)].sup.T.

6. A real-time identification device for a nonlinear characteristic of a model-free in-service seismic isolation/vibration reduction device driven based on monitoring data, comprising: a dividing module, for dividing a structure of a seismic isolation/vibration reduction system into a plurality of substructures, and defining a substructure where the seismic isolation/vibration reduction device is located as a target substructure; and a GEKF-UI identification module, for using a GEKF-UI to identify a linear stiffness and a damping coefficient of the seismic isolation/vibration reduction device and a linear stiffness and a damping coefficient of the target substructure based on monitoring data, in the case that an external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state.

7. The device according to claim 6, wherein the device further comprises: a GKF-UI identification module, for using a GKF-UI to identify a nonlinear restoring force received by the target substructure by the linear stiffness and the damping coefficient of the seismic isolation/vibration reduction device and the linear stiffness and the damping coefficient of the target substructure identified by the GEKF-UI identification module, in the case that the external excitation on the seismic isolation/vibration reduction device is large so that the seismic isolation/vibration reduction device generates a nonlinear force.

8. The device according to claim 6, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+(C.sup.r+C){dot over (x)}.sup.r(t)+(K.sup.r+K)x.sup.r(t)=.sup.rf.sup.r(t)+.sup.gf.sup.g(t)=.sup.uf.sup.u(t) wherein, M.sup.r, C.sup.r and K.sup.r are mass, damping and stiffness matrices of the target substructure, respectively; {umlaut over (x)}.sup.r(t), {dot over (x)}.sup.r(t) and x.sup.r (t) are acceleration, velocity and displacement vectors of the target substructure, respectively; f.sup.r(t) is the external excitation received by the target substructure; f.sup.g(t) is a force of an adjacent substructure on the target substructure; .sup.r and .sup.g are position matrices of the external excitation and the force; C and K are additional damping and stiffness provided by the seismic isolation/vibration reduction device to the system structure; .sup.u is a position matrix of the external excitation f.sup.u(t); an extended vector Z=[x.sup.T {dot over (x)}.sup.T .sup.T].sup.T is established, and the GEKF-UI is used to obtain to express the linear stiffness and the damping coefficient of the seismic isolation/vibration reduction device and the linear stiffness and the damping coefficient of the target substructure.

9. The device according to claim 6, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is large so that the seismic isolation/vibration reduction device generates a nonlinear force, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+C.sup.r{dot over (x)}.sup.r(t)+K.sup.rx.sup.r(t)=.sup.rf.sup.r(t)+.sup.gg.sup.r(t).sup.nonf.sup.non(t)=.sup.uf.sup.u(t), wherein C.sup.r and K.sup.r are identified values in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state; the GKF-UI is used to identify f.sup.u(t) by observing a response of a partial structure; the nonlinear force generated by the vibration reduction/seismic isolation device is identified based on the identified f.sup.u(t) and f.sup.u=[f.sup.r(t) g.sup.r(t) f.sup.non(t)].sup.T.

10. The device according to claim 7, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+(C.sup.r+C){dot over (x)}.sup.r(t)+(K.sup.r+K)x.sup.r(t)=.sup.rf.sup.r(t)+.sup.gf.sup.g(t)=.sup.uf.sup.u(t) wherein, M.sup.r, C.sup.r and K.sup.r are mass, damping and stiffness matrices of the target substructure, respectively; {umlaut over (x)}(t), {dot over (x)}.sup.r(t) and x.sup.r(t) are acceleration, velocity and displacement vectors of the target substructure, respectively; f.sup.r(t) is the external excitation received by the target substructure; f.sup.g(t) is a force of an adjacent substructure on the target substructure; .sup.r and .sup.g are position matrices of the external excitation and the force; C and K are additional damping and stiffness provided by the seismic isolation/vibration reduction device to the system structure; .sup.u is a position matrix of the external excitation f.sup.u(t); an extended vector Z=[x.sup.T {dot over (x)}.sup.T .sup.T].sup.T is established, and the GEKF-UI is used to obtain to express the linear stiffness and the damping coefficient of the seismic isolation/vibration reduction device and the linear stiffness and the damping coefficient of the target substructure.

11. The device according to claim 7, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is large so that the seismic isolation/vibration reduction device generates a nonlinear force, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+C.sup.r{dot over (x)}.sup.r(t)+K.sup.rx.sup.r(t)=.sup.rf.sup.r(t)+.sup.gg.sup.r(t).sup.nonf.sup.non(t)=.sup.uf.sup.u(t) wherein C.sup.r and K.sup.r are identified values in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state; the GKF-UI is used to identify f.sup.u(t) by observing a response of a partial structure; the nonlinear force generated by the vibration reduction/seismic isolation device is identified based on the identified f.sup.u(t) and f.sup.u=[f.sup.r(t) g.sup.r(t) f.sup.non(t)].sup.T.

12. The method according to claim 2, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+(C.sup.r+C){dot over (x)}.sup.r(t)+(K.sup.r+K)x.sup.r(t)=.sup.rf.sup.r(t)+.sup.gf.sup.g(t)=.sup.uf.sup.u(t) wherein, M.sup.r, C.sup.r and K.sup.r are mass, damping and stiffness matrices of the target substructure, respectively; {umlaut over (x)}(t), {dot over (x)}.sup.r(t) and x.sup.r (t) are acceleration, velocity and displacement vectors of the target substructure, respectively; f.sup.r(t) is the external excitation received by the target substructure; f.sup.g(t) is a force of an adjacent substructure on the target substructure; .sup.r and .sup.g are position matrices of the external excitation and the force; C and K are additional damping and stiffness provided by the seismic isolation/vibration reduction device to the system structure; .sup.u is a position matrix of the external excitation f.sup.u(t).

13. The method according to claim 12, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state, an extended vector Z=[x.sup.T {dot over (x)}.sup.T .sup.T].sup.T is established, and the GEKF-UI is used to obtain to express the linear stiffness and the damping coefficient of the seismic isolation/vibration reduction device and the linear stiffness and the damping coefficient of the target substructure.

14. The method according to claim 13, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is large so that the seismic isolation/vibration reduction device generates a nonlinear force, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+C.sup.r{dot over (x)}.sup.r(t)+K.sup.rx.sup.r(t)=.sup.rf.sup.r(t)+.sup.gg.sup.r(t).sup.nonf.sup.non(t)=.sup.uf.sup.u(t) wherein C.sup.r and K.sup.r are identified values in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state; the GKF-UI is used to identify f.sup.u(t) by observing a response of a partial structure; the nonlinear force generated by the vibration reduction/seismic isolation device is identified based on the identified f.sup.u(t) and f.sup.u=[f.sup.r(t) g.sup.r(t) f.sup.non(t)].sup.T.

15. The method according to claim 4, wherein in the case that the external excitation on the seismic isolation/vibration reduction device is large so that the seismic isolation/vibration reduction device generates a nonlinear force, a motion equation of the target substructure is:
M.sup.r{umlaut over (x)}.sup.r(t)+C.sup.r{dot over (x)}.sup.r(t)+K.sup.rx.sup.r(t)=.sup.rf.sup.r(t)+.sup.gg.sup.r(t).sup.nonf.sup.non(t)=.sup.uf.sup.u(t) wherein C.sup.r and K.sup.r are identified values in the case that the external excitation on the seismic isolation/vibration reduction device is small so that the seismic isolation/vibration reduction system is in a linear state; the GKF-UI is used to identify f.sup.u(t) by observing a response of a partial structure; the nonlinear force generated by the vibration reduction/seismic isolation device is identified based on the identified f.sup.u(t) and f.sup.u=[f.sup.r(t) g.sup.r(t) f.sup.non(t)].sup.T.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0019] The accompanying drawings illustrate exemplary implementations of the disclosure, and are intended to explain the principles of the disclosure together with the description thereof. The accompanying drawings are included to provide a further understanding of the disclosure, and are included in and constitute part of the specification.

[0020] FIG. 1 shows substructure division of a frame equipped with a seismic isolation device.

[0021] FIG. 2 shows substructure division of a frame equipped with a vibration reduction device.

[0022] FIG. 3 shows a comparison of numerical results of nonlinear force identification of a seismic isolation device.

[0023] FIG. 4 shows a comparison of numerical results of nonlinear force identification of a vibration reduction device.

[0024] FIG. 5 shows a comparison of experimental results of hysteresis loop identification of a nonlinear force of a negative stiffness damper (NSD).

[0025] FIG. 6 shows a comparison of experimental results of nonlinear force identification of a vibration reduction device.

DETAILED DESCRIPTION

[0026] The present invention is described in further detail below with reference to the accompanying drawings and implementations. It should be understood that the specific implementations described herein are merely intended to explain related content, rather than to limit the present disclosure. It should also be noted that, for convenience of description, only the parts related to the present disclosure are shown in the accompany drawings.

[0027] It should be noted that, in case of no conflict, the implementations in the present disclosure and the features in the implementations may be combined with each other. The present disclosure is described in detail below with reference to the accompanying drawings and implementations.

[0028] The present disclosure provides a real-time identification method for a nonlinear characteristic of a model-free in-service seismic isolation/vibration reduction device driven based on monitoring data.

[0029] The method may include the following steps. Step 1: in case an external excitation on the seismic isolation/vibration reduction device is small, an entire seismic isolation/vibration reduction system is in a linear state; at this time, the seismic isolation/vibration reduction device only provides a linear restoring force; a system structure is divided into a plurality of substructures, and a target substructure where the seismic isolation/vibration reduction device is located is defined as a target substrate for analysis; a General Extended Kalman filter with unknown inputs (GEKF-UI) is used to identify a linear stiffness and a damping coefficient of the seismic isolation/vibration reduction device and a linear stiffness and a damping coefficient of the target substructure. Step 2: under a large external excitation, the seismic isolation/vibration reduction device generates a nonlinear force, and the system structure is still in a linear state; the nonlinear restoring force generated by the seismic isolation/vibration reduction device is regarded as an unknown additional virtual force applied to the target substructure; the target substructure is analyzed; the target substructure is subject to the nonlinear restoring force, the external excitation and an interaction force between the substructures. a General Kalman filter with unknown inputs (GKF-UI) is used to identify the nonlinear restoring force, the external excitation and the interaction force between the substructures.

[0030] In this way, the method proposed by the present disclosure only observes a response of a partial structure, and is applicable regardless whether an acceleration response is observed at an action position of the nonlinear force, an action position of the seismic isolation/vibration reduction device and a connection position between the substructures.

[0031] A frame structure is taken as an example. FIG. 1 shows substructure division of a frame equipped with a seismic isolation device, and FIG. 2 shows substructure division of a frame equipped with a vibration reduction device. Due to substructure division, a substructure is separated from an overall structure. For example, an r-th substructure is separated from the overall structure. The substructure is affected by a force from other substructures where the substructure is connected to other substructures. Therefore, the r-th substructure is taken out for parameter identification, and an interaction force between the substructures is unknown. A motion equation of the substructure is:

M.sup.r{umlaut over (x)}.sup.r(t)+C.sup.r{dot over (x)}.sup.r(t)+K.sup.rx.sup.r(t)+.sup.nonf.sup.non(t)=.sup.rf.sup.r(t)+.sup.gf.sup.g(t)(1)

[0032] where, M.sup.r, C.sup.r and K.sup.r are mass, damping and stiffness matrices of the r-th substructure, respectively; {umlaut over (x)}.sup.r(t), {dot over (x)}.sup.r(t) and x.sup.r(t) are acceleration, velocity and displacement vectors of the r-th substructure, respectively; f.sup.non(t) is a restoring force provided by the seismic isolation/vibration reduction device; f.sup.r(t) is an external excitation received by the r-th substructure; f.sup.g(t) is a force of an adjacent substructure on a target substructure (the r-th substructure); .sup.non, .sup.r and .sup.g are corresponding position matrices of the restoring force, the external excitation and the force.

[0033] A seismic isolation/vibration reduction device installed in a large structure will generate a nonlinear force under a large external excitation. For a locally nonlinear structure, substructures including a substructure where the seismic isolation/vibration reduction device is located are regarded as target substructures, and the target sub-structures are separately nonlinearly identified. In this way, the identification efficiency is greatly improved.

[0034] When the external excitation is small, the seismic isolation/vibration reduction device is in a linear state. f.sup.non(t) is a linear restoring force, which means that the seismic isolation/vibration reduction device provides additional stiffness and damping to the structure. The motion equation (1) can be:

M.sup.r{umlaut over (x)}.sup.r(t)+(C.sup.r+C){dot over (x)}.sup.r(t)+(K.sup.r+K)x.sup.r(t)=.sup.rf.sup.r(t)+.sup.gf.sup.g(t)=.sup.uf.sup.u(t)(2)

[0035] where, C and K are the additional stiffness and damping provided by the seismic isolation/vibration reduction device to the structure, respectively.

[0036] An extended state vector

[00001]$H_{k+1k}=\frac{\stackrel{\_}{h}\left(Z\right)}{Z^T}.Math.{}_{Z={\tilde{Z}}_{k+1k}}$

is established:

[00002]$\begin{array}{cc}\stackrel{.}{Z}=\left[\begin{array}{c}\stackrel{.}{x}\\ \stackrel{\mathrm{.Math.}}{x}\\ \stackrel{.}{}\end{array}\right]=\left[\begin{array}{c}\stackrel{.}{x}\\ {\left(M^r\right)}^{-1}.Math.\left(-\left(K^r+K^{}\right)-\left(C^r+C^{}\right).Math.\stackrel{.}{x}+{}^u.Math.f^u\right)\\ 0\end{array}\right]=\stackrel{\_}{g}\left(Z\right)+\stackrel{\_}{B}.Math.f^u& \left(3\right)\end{array}$

[0037] where, g(.Math.) represents a nonlinear function, B=[0 (M.sup.r).sup.1.sup.u 0].sup.T, and represents a structural parameter.

[0038] A supplementary observation equation is as follows:

y=h(Z)+Df.sup.u(4)

[0039] where, h(.Math.) represents a nonlinear function, and D represents a position influence matrix of a known external excitation f.sup.u.

[0040] A first-order Taylor expansion is performed on g(Z) and h(Z).

g(Z)=g({circumflex over (Z)}.sub.k|k)+G.sub.k|k(Z{circumflex over (Z)}.sub.k|k)(5)

h(Z)=h({tilde over (Z)}.sub.k+1|k)+H.sub.k+1|k(Z{tilde over (Z)}.sub.k+1|k)(6)

[0041] where,

[00003]$G_{kk}=\frac{\stackrel{\_}{g}\left(Z\right)}{Z^T}.Math.{}_{Z={\hat{Z}}_{kk}},H_{k+1k}=\frac{\stackrel{\_}{h}\left(Z\right)}{Z^T}.Math.{}_{Z={\tilde{Z}}_{k+1k}},$

{tilde over (Z)}.sub.k|k represents an optimal state estimation at time k, and {tilde over (Z)}.sub.k+1|k represents a state forecast at time k+1.

{tilde over (Z)}.sub.k+1|kA.sub.k{circumflex over (Z)}.sub.k|k+B.sub.k*{circumflex over (f)}.sub.k|k.sup.u+g.sub.k|k(7)

[0042] where, A.sub.k is a state-transition matrix, B.sub.k*=(A.sub.kI) (G.sub.k|kt).sup.1 (Bt) g.sub.k|k=[g({circumflex over (Z)}.sub.k|k)G.sub.k|k{circumflex over (Z)}.sub.k|k]t.

[0043] Assuming that the external excitation f.sup.u adopts first-order hold (FOH dispersion) in a sampling interval, the state equation and the observation equation are discretized as follows:

Z.sub.k+1=A.sub.kZ.sub.k+B.sub.kf.sub.k.sup.u+B.sub.k+1f.sub.k+1.sup.u+g.sub.k|k+w.sub.k(8)

y.sub.k+1=H.sub.k+1|kZ.sub.k+1+h({tilde over (Z)}.sub.k+1|k)H.sub.k+1|k{tilde over (Z)}.sub.k+1|k+Df.sub.k+1.sup.u+v.sub.k+1(9)

[0044] where, Z.sub.k+1 is a state vector at time k+1; Z.sub.k is a state vector at time k, B.sub.k=Bt, B.sub.k+1=Bt; f.sub.k.sup.u is an unknown external excitation at time k; f.sub.k+1.sup.u is an unknown external excitation at time k+1; w.sub.k is a model error, with a mean of 0 and a covariance of Q.sub.k; u.sub.k+1 is an observation vector at time k+1; v.sub.k+1 is a measurement error, with a mean of 0 and a variance of R.sub.k+1.

[0045] An optimal state estimation at time k+1 is:

{circumflex over (Z)}.sub.k+1|k+1=Z.sub.k+1|k+K.sub.k+1(y.sub.k+1H.sub.k+1|kZ.sub.k+1|kD{circumflex over (f)}.sub.k+1|k+1.sup.uh.sub.k+1|k)(10)

where, h.sub.k+1|k=h({tilde over (Z)}.sub.k+1|k)H.sub.k+1|k{tilde over (Z)}.sub.k+1|k(11)

Z.sub.k+1|kA.sub.k{circumflex over (Z)}.sub.k|k+B.sub.k{circumflex over (f)}.sub.k|k.sup.uB.sub.k+1+{circumflex over (f)}.sub.k+1|k+1.sup.u+g.sub.k|k(12)

[0046] where, K.sub.k+1 is a Kalman gain matrix;

[00004]$\begin{array}{cc}{K}_{k+1}={\stackrel{.Math.}{P}}_{k+1k}^Z.Math.{H_{k+1k}^T\left(H_{k+1k}.Math.{\stackrel{\_}{P}}_{k+1k}^Z.Math.H_{k+1k}^T+R_{k+1}\right)}^{-1}& \left(13\right)\\ {\stackrel{.Math.}{P}}_{k+1k}^Z=\left[\begin{array}{cc}{A}_k& B_k\end{array}\right]\left[\begin{array}{cc}{\hat{P}}_{kk}^Z& {\hat{P}}_{kk}^{\mathrm{Zf}}\\ {\hat{P}}_{kk}^{\mathrm{fZ}}& {\hat{P}}_{kk}^f\end{array}\right]\left[\begin{array}{c}{A}_k^T\\ B_k^T\end{array}\right]+Q_k& \left(14\right)\end{array}$

[0047] By minimizing the error vector, an unknown external excitation {circumflex over (f)}.sub.k+1|k+1.sup.u is obtained:

{circumflex over (f)}.sub.k+1|k+1.sup.u=S.sub.k+1[y.sub.k+1H.sub.k+1|k(A.sub.k{circumflex over (Z)}.sub.k|k+B.sub.k{circumflex over (f)}.sub.k|k.sup.ug.sub.k|k)h.sub.k+1|k](15)

where

S.sub.k+1=(H.sub.k+1|kB.sub.k+1+D).sup.TR.sub.k+1.sup.1(IH.sub.k+1|kK.sub.k+1)(H.sub.k+1|kB.sub.k+1D).sup.1(H.sub.k+1B.sub.k+1D).sup.TR.sub.k+1.sup.1(IH.sub.k+1|kK.sub.k+1(16)

[0048] By observing a partial response, the GEKF-UI method identifies the unknown external excitation, the interaction force between the substructures and the extended vector Z=[x.sup.T {dot over (x)}.sup.T .sup.T].sup.T, where represents the linear stiffness and the damping coefficient of the seismic isolation/vibration reduction device and the linear stiffness and the damping coefficient of the target substructure. These physical coefficients identified will be used in step 2.

[0049] The GKF-UI method is similar to the GEKF-UI method, except that the GKF-UI method is used with a known structural parameter. When the external excitation is large, the seismic isolation/vibration reduction device generates a nonlinear force. The nonlinear force provided by the seismic isolation/vibration reduction device is introduced to the right of the equation to form an unknown additional virtual force, and the equation (1) is rewritten as:

M.sup.r{umlaut over (x)}.sup.r(t)+C.sup.r{dot over (x)}.sup.r(t)+K.sup.rx.sup.r(t)=.sup.rf.sup.r(t)+.sup.gg.sup.r(t).sup.nonf.sup.non(t)=.sup.uf.sup.u(t)(17)

[0050] where, C.sup.r and K.sup.r are the values identified by the GEKF-UI in step 1. In order to overcome the limit of observing a response at an excitation point, as in the previous section, assuming that the unknown external excitation f.sup.u adopts FOH in the sampling interval, the motion equation is discretized into the following form:

X.sub.k+1=A.sub.kX.sub.k+B.sub.kf.sub.k.sup.u+G.sub.k+1f.sub.k+1.sup.uw.sub.k(18)

[0051] where, X.sub.k+1 and X.sub.k are state vectors of time t=(k+1) t, kt, respectively; t is a sampling time interval; A.sub.k is a state-transition matrix; B.sub.k and G.sub.k+1 are influence matrices of an unknown force. w.sub.k is a model error, with a mean of 0 and a covariance of Q.sub.k.

[0052] By observing only the partial structure response, the discretized observation equation can be expressed as:

y.sub.k+1=C.sub.k+1X.sub.k+1+H.sub.k+1.sup.uf.sub.k+1.sup.u+b.sub.k+1(19)

[0053] where, y.sub.k+1 is an observation vector; C.sub.k+1 is a state observation matrix; H.sub.k+1.sup.u is an observation matrix of an unknown force f.sub.k+1.sup.u; v.sub.k+1 is an observation noise vector, with a mean of 0 and a covariance of R.sub.k+1.

[0054] A state forecast equation and a state estimation equation are established as follows:

{tilde over (X)}.sub.k+1|k=A.sub.k{circumflex over (X)}.sub.k|k+B.sub.k{circumflex over (f)}.sub.k|k.sup.uG.sub.k+1{circumflex over (f)}.sub.k+1|k+1.sup.u(20)

{circumflex over (X)}.sub.x+1|k+1={tilde over (X)}.sub.k+1|k+K.sub.k+1(y.sub.k+1C.sub.k+1{tilde over (X)}.sub.k+1|kH.sub.k+1.sup.u{circumflex over (f)}.sub.k+1|k+1.sup.u) (21)

[0055] where, K.sub.k+1 is a Kalman gain matrix;

[00005]$\begin{array}{cc}{K}_{k+1}={\stackrel{.Math.}{P}}_{k+1k}^X.Math.{C_{k+1}^T\left(C_{k+1}.Math.{\tilde{P}}_{k+1k}^X.Math.C_{k+1}^T+R_{k+1}\right)}^{-1}& \left(22\right)\\ {\stackrel{.Math.}{P}}_{k+1k}^X=\left[\begin{array}{cc}{A}_k& B_k\end{array}\right]\left[\begin{array}{cc}{\hat{P}}_{kk}^X& {\hat{P}}_{kk}^{X.Math..Math.f}\\ {\hat{P}}_{kk}^{\mathrm{fZ}}& {\hat{P}}_{kk}^f\end{array}\right]\left[\begin{array}{c}{A}_k^T\\ B_k^T\end{array}\right]+Q_k& \left(23\right)\end{array}$

[0056] By minimizing the observation equation, the unknown external excitation {circumflex over (f)}.sub.k+1|k+1.sup.u is obtained:

{circumflex over (f)}.sub.k+1|k+1.sup.u=M.sub.k+1[y.sub.k+1C.sub.k+1(A.sub.k{circumflex over (X)}.sub.k|k+B.sub.k{circumflex over (f)}.sub.k|k.sup.u)](24)

where

M.sub.k+1=[(C.sub.k+1G.sub.k+1+H.sub.k+1.sup.u).sup.TR.sub.k+1.sup.1(IC.sub.k+1K.sub.k+1)(C.sub.k+1G.sub.k+1+K.sub.k+1.sup.u)].sup.1(C.sub.k+1G.sub.k+1+H.sub.k+1.sup.u).sup.TR.sub.k+1.sup.1(IC.sub.k+1K.sub.k+1(25)

[0057] where, I is an identity matrix.

[0058] The GKF-UI method identifies f.sup.u(t) in real time by observing a partial response. f.sup.u=[f.sup.r(t) g.sup.r(t) f.sup.non(t)].sup.T, so that the required nonlinear force generated by the vibration reduction/seismic isolation device, the interaction force between the substructures and the unknown external excitation can be identified.

[0059] Based on the above method of the present disclosure, the technical effects are numerically verified as follows.

[0060] (1) Numerical verification of seismic isolation bearing identification

[0061] A superstructure adopted was a 9-story shear frame, an El-Centro seismic wave was used as a seismic action, and a Bouc-Wen nonlinear model of hysteresis was used to simulate a force-deformation relationship of the seismic isolation bearing. The Bouc-Wen nonlinear model of hysteresis expressed a restoring force of the seismic isolation bearing as:

[00006]$\begin{array}{cc}.Math.R_b^{\mathrm{Bouc}-\mathrm{Wen}}={}_b.Math.k_b.Math.x_b+\left(1-{}_b\right).Math.k_b.Math.z_b^{\mathrm{Bouc}-\mathrm{Wen}}& \left(26\right)\\ {\stackrel{.}{z}}_b^{\mathrm{Bouc}-\mathrm{Wen}}={\stackrel{.}{x}}_b-\left\{{}_b.Math..Math.{\stackrel{.}{x}}_b.Math..Math.{.Math.z_b^{\mathrm{Bouc}-\mathrm{Wen}}.Math.}^{n_b-1}.Math.z_b^{\mathrm{Bouc}-\mathrm{Wen}}+{}_b.Math.{\stackrel{.}{x}}_b.Math.{.Math.z_b^{\mathrm{Bouc}-\mathrm{Wen}}.Math.}^{n_b}\right\}& \left(27\right)\end{array}$

[0062] Structural parameters of a seismic isolation layer: lumped mass of the seismic isolation layer m.sub.b=65 kg initial stiffness k.sub.b=1010.sup.5 N/m, and viscous damping coefficient c.sub.b=800 N.Math.s/m Parameters of the superstructure: element lumped mass m.sub.i=60 kg, element linear stiffness k.sub.i=1.210.sup.5 N/m, and element viscous damping coefficient c.sub.i=1000 N.Math.s/m (i=1,2, . . . , 9) Parameters of the Bouc-Wen nonlinear model of hysteresis: .sub.b=0.1, .sub.b=2000, .sub.b=2000, n.sub.b=1.25.

[0063] The results of the GEKF-UI in step 1 are shown in Table 1:

TABLE-US-00001 TABLE 1 Identification results of linear parameters of structure and damper Stiff- Iden- Iden- ness Actual tified Error Damping Actual tified Error (kN/m) Value Value (%) (N .Math. s/m) Value Value (%) K.sub.b 100 100.00 0.00 c.sub.b 800 809 1.18 k.sub.1 120 120.00 0.00 c.sub.1 1000 985 1.48 k.sub.2 120 120.06 0.05 c.sub.2 1000 1012 1.17 k.sub.3 120 119.92 0.06 c.sub.3 1000 986 1.42 k.sub.4 120 120.08 0.06 c.sub.4 1000 1008 0.79 k.sub.5 120 119.92 0.06 c.sub.5 1000 991 0.92 k.sub.6 120 120.04 0.04 c.sub.6 1000 1005 0.50 k.sub.7 120 119.95 0.04 c.sub.7 1000 992 0.77 k.sub.8 120 120.04 0.03 c.sub.8 1000 1006 0.59 k.sub.9 120 119.81 0.16 c.sub.9 1000 984 1.56

[0064] The results of the GKF-UI in step 2 are shown in FIG. 3. FIG. 3 shows a comparison of numerical results of nonlinear force identification of the seismic isolation device.

[0065] (2) Numerical Verification of Vibration Damper Identification

[0066] A model of a 10-story shear frame was externally excited by a white noise acting on a 9.sup.th story. A vibration damper was installed between 8.sup.th and 9.sup.th stories. A Dahl nonlinear model was used to simulate a force-deformation relationship of the damper, which is expressed as follows:

F.sub.Dahl(t)=k.sup.mrx.sub.i(t)+c.sup.mr{dot over (x)}.sub.i(t)+f.sub.dz.sub.i.sup.Dahl(t)+f.sub.0(28)

[0067] Structural parameters: element lumped mass m.sub.i=60 kg, element linear stiffness k.sub.i=1.210.sup.5 N/m, (i=1,2, . . . , 10); structural damping C=M.sub.0+K which is Rayleigh damping; first and second-order damping ratios =0.3; damping coefficient =0.3002, =0.002257; and mass of the damper m.sup.mr=2 kg Nonlinear model parameters of the damper: k.sup.mr=410.sup.4 N/m, c.sup.mr=100 N.Math.s/m, f.sub.d=200 N, f.sub.0=0, =4000 s/m.

[0068] The results of the GEKF-UI in step 1 are shown in Tables 2 and 3:

TABLE-US-00002 TABLE 2 Identification results of linear parameters of structure Structural Structural Stiffness Actual Identified Error Damping Actual Identified Error (KN/m) Value Value (%) Coefficient Value Value (%) k.sub.1 120 119.33 0.56 0.3002 0.2963 1.30 k.sub.2 120 119.80 0.16 0.002257 0.002340 3.70 k.sub.3 120 120.12 0.10 k.sub.4 120 120.11 0.09 k.sub.5 120 119.90 0.08 k.sub.6 120 120.25 0.21 k.sub.7 120 120.10 0.09 k.sub.8 120 120.04 0.03 k.sub.9 120 124.80 4.00 k.sub.10 120 120.89 0.75

TABLE-US-00003 TABLE 3 Identification results of linear parameters of damper Stiffness of Damping Damper Actual Identified of Damper Actual Identified (kN/m) Value Value Error (%) (N .Math. s/m) Value Value Error (%) k.sup.mr 40 36.32 9.20 c.sup.mr 100 104.75 4.75

[0069] The results of the GKF-UI in step 2 are shown in FIG. 4. FIG. 4 shows a comparison of numerical results of nonlinear force identification of the vibration reduction device.

[0070] Based on the above method of the present disclosure, the technical effects are experimentally verified as follows.

[0071] (1) Experimental Verification of Seismic Isolation Bearing Identification

[0072] Experimental data provided by Professor Satish Nagarajaiah of Rice University in the United States were used to identify a nonlinear force generated by a three-story experimental shear frame equipped with a negative stiffness damper (NSD) on a first story under an earthquake action.

[0073] A known element lumped mass was m.sub.i=8.6 kips (i=1,2,3). A structural initial linear element stiffness was identified by an extended Kalman filter (EKF) as k.sub.1=8.9041 kip/in, k.sub.2=14.0061 kip/in, k.sub.3=17.9544 kip/in and a damping coefficient was identified as =0.5202, =0.0017. The identification results of the linear physical parameters were used to identify a nonlinear force of the NSD.

[0074] A hysteresis loop of the nonlinear force of the NSD was identified by the GKF-UI, as shown in FIG. 5. FIG. 5 shows a comparison of experimental results of hysteresis loop identification of the nonlinear force of the NSD.

[0075] (2) Experimental Verification of Vibration Damper Identification

[0076] A five-story shear frame structure equipped with a magneto-rheological (MR) damper was experimented. An electromagnetic exciter acted on a 3.sup.rd story of the structure, and continuously excited for 6 s, with a sampling frequency of 1,000 Hz.

[0077] The results of the GEKF-UI in step 1 are shown in Table 4:

TABLE-US-00004 TABLE 4 Identification results of linear parameters of structure Structural Structural Stiffness Actual Idntified Error Damping Actual Identified Error (kN/m) Value Value (%) Coefficient Value Value (%) k1 128.43 131.1 2.05 0.14 0.1335 4.6 k2 138.49 139.7 0.92 6.8*10{circumflex over ()} 6.3*10{circumflex over ()} 7.3 6 6 k3 129.34 134.3 3.82 k4 126.72 120.2 5.14 k5 126.96 122.4 3.57

[0078] The results of the GKF-UI in step 2 are shown in FIG. 6. FIG. 6 shows a comparison of experimental results of nonlinear force identification of the vibration reduction device.

[0079] In the description of this specification, the description of the terms one embodiment/implementation, some embodiments/implementations, example, specific example and some examples etc. means that the specific features, structures, materials or characteristics described with reference to the embodiment/implementation or example are included in at least one embodiment/implementation or example of the present application. In this specification, the illustrative expressions of the above terms are not intended to refer to the same embodiment/implementation or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments/implementations or examples. In addition, those skilled in the art may combine different embodiments/implementations or examples described herein or features in different embodiments/implementations or examples without any contradiction.

[0080] Those skilled in the art should understand that the foregoing implementations are merely intended to describe the present disclosure clearly, rather than to limit the scope of the present disclosure. Those skilled in the art may make other changes or modifications based on the foregoing disclosure, but these changes or modifications should fall within the scope of the present disclosure.