Study on seismic response characteristics of the interlayer isolation structure

Zhang, Shangrong; Liu, Fanglan; Hu, Yuchen; Li, Shihao; Zhu, Lei

doi:10.21595/jve.2021.22003

Journal of Vibroengineering

Browse Journal

Submit article

Published: 30 August 2021

Check for updates

Study on seismic response characteristics of the interlayer isolation structure

Shangrong Zhang¹

Fanglan Liu²

Yuchen Hu³

Shihao Li⁴

Lei Zhu⁵

^{1, 2, 3, 4, 5}School of Civil and Hydraulic Engineering, Ningxia University, Yinchuan, 750021, China

²Northwest Engineering Corporation Limited, Powerchina, Xi’an, 710065, China

Corresponding Author:

Fanglan Liu

Cite the article Download PDF

Downloads 787

WoS Core Citations 4

CrossRef Citations 3

Abstract

Based on Timoshenko beam theory and considering the layer bending deformation and shear deformation of the structure, a semi analytical solution method (SAM) is proposed for the interlayer isolation structure system. The equivalent mechanical model of cantilever beam is established by using the distributed parameter system, on which the dynamic partial differential equation is derived; the boundary condition and continuity condition are used to solve the equation to obtain the natural vibration characteristics; the orthogonal conditions are deduced by Betti’s law and the seismic response of the structure can be solved by the mode superposition method. A 10 story isolated structure is designed as well as analysed by the semi analytical method (SAM) and the ETABS finite element simulation (FESE) separately. The results show that: the existence of the isolation layer may reduce the response of the superstructure and amplify the response of the substructure, and the amplification can be solved by increasing the mass, bending stiffness of the superstructure and shear stiffness of the isolation layer. In regular interlayer isolation structure, the optimal frequency ratio of the superstructure and substructure increases with the rising of isolation layer position. The results also prove the correctness and accuracy of the SAM, which enriches the analysis theory of interlayer isolation structure and provides reference for the accurate analysis of structural response in engineering.

Study on seismic response characteristics of the interlayer isolation structure

Highlights

Using the distributed parameter system, a semi analytical solution method (SAM) is proposed for the interlayer isolation structure.
Through the comparison with the solution results of finite element method for the interlayer isolation structure, the SAM is proved to be correct and effective.
The existence of the isolation layer may reduce the response of the superstructure and amplify the response of the substructure, and the amplification can be solved by increasing the mass, bending stiffness of the superstructure and shear stiffness of the isolation layer.

1. Introduction

In recent years, researchers continue to improve the theoretical system of interlayer isolation structure on the basis of engineering practice [1-14]. At present, the particle system model is mainly used in the dynamic analysis of interlayer isolation structure. Zhou Fulin et al. [15] used two particle simplified model to study the mechanical characteristics of interlayer isolation structure, and obtained the change of isolation mechanism when the location of isolation layer is different; Cao Jialiang [16], He Dongdong [17] also used two particle model to study the model parameter sensitivity and model parameter optimization of interlayer isolation structure. Based on the multi-particle layer shear model, using the random vibration analysis method, Qi Ai [18] obtains the relationship between the isolation effect of the interlayer isolation structure and the damping of the isolation layer and the frequency ratio of each substructure. However, the actual interlayer isolation structure is a continuum. Its mass, stiffness and damping are continuously distributed in space, which makes the structure have infinite degrees of freedom. Therefore, for the equivalence of the continuous structure, the distributed parameter system model is closer to the reality than the particle system model. When using the theory of parameter distribution system to simplify the interlayer isolation structure to beam model for dynamic analysis, according to different assumptions, Euler Bernoulli beam theory or Timoshenko beam theory can be used [19]. Compared with Euler Bernoulli beam theory, Timoshenko beam theory takes shear action into account, and has higher accuracy and wider application range of slenderness ratio [20]. When the position of the isolation layer rises gradually, the deformation of the substructure and superstructure of the interlayer isolation structure changes from the bending type to the bending shear type, and the influence of the shear deformation on the higher-order mode is a source of the modal error that cannot be ignored [21]. Therefore, it is more accurate and practical to adopt Timoshenko beam theory to establish the distributed parameter system model of interlayer isolation structure.

For the distributed parameter system, there are few researches on the interlayer isolation structure, but it has been mature in other scientific fields. Li Yuanchun [22, 23] calculated the kinetic energy, potential energy and virtual work of the whole system for the flexible loading system of manipulator, and gave the distributed parameter dynamic model reflecting the essence of the system by using Hamiltonian principle. Based on the distributed parameter system, Chen Yangyang et al. [24, 25] analyzed the single pier particle system of the high pier beam bridge, and obtained the analytical formula with higher accuracy and wider application range for studying the factors affecting the mode mass participation coefficient of the system; Also based on the vibration theory of distributed parameter beam, Du Yongfeng [26, 27] deduced and verified the semi analytical method of seismic response of the isolation system of series and parallel electrical equipment support, which provided the thinking and reference for the study of parameter model analysis of building structure such as interlayer isolation; Song Xiao et al. [28] established the distributed parameter model of the interlayer isolation structure, and used the virtual excitation method to solve the seismic response. It was found that the participation of the isolation layer adjusted the quality participation coefficient of each vibration mode of the structure, and then affected the isolation effect.

In this paper, the Timoshenko beam theory is used to establish the equivalent mechanical model of structural distribution parameters. Then the dynamic equation is deduced, the natural frequency and mode are calculated, and the mode decomposition method is used to solve the response of the interlayer isolation structure. Finally, for a 10 story isolated structure, the correctness of the model and the structure solution method is verified by the finite element method (FEM) and the ETABS finite element software simulation (FESE). In addition, the influence of the parameters of substructures on the isolation effect is also studied. There are many symbols in this paper. In order to make it more convenient to read, the important symbols and meanings are listed in Table 1. If there is special case, it will be explained in the text.

Table 1Symbols and the corresponding meanings

Symbol	Meaning	Symbol	Meaning
$m$	Mass of the beam	$ω$	Natural frequency of the structure
$E I$	Bending stiffness of the beam	$ϕ$	Mode function of the structure
$G A$	Shear stiffness of the beam	${\ddot{u}}_{g}$	Acceleration of ground motion
$K$	stiffness	$C$	The damping coefficients of horizontal velocity
$J$	Moment of inertia of the beam	$C_{S}$	The damping coefficients of strain velocity
$l$	Length of the beam	$C_{J}$	The damping coefficients of rotational velocity
$ξ$	Damping ratio
Subscript	Meaning	superscript	Meaning
1	Substructure	$\bar{}$	Per unit length
2	superstructure
iso	Isolation layer

2. Establishment of dynamic analysis model

2.1. Basic assumptions

Based on the theory of distributed parameter system, the mass, stiffness and damping of the substructure, the isolation layer and the superstructure of interlayer isolation structure are regarded as the uniformly distributed parameters,

It is also assumed that:

(1) The rigidity of the floor is infinite in the plane;

(2) The axial deformation of each component is ignored;

(3) The torsional deformation of the structure is ignored.

2.2. Establishment of analysis model

Based on Timoshenko beam theory and considering its shear deformation and bending deformation, the interlayer isolation structure is simplified as a cantilever beam, and the mass, stiffness and damping of each substructure are uniformly distributed along the axis. Its equivalent mechanical model is shown in Fig. 1. Among them, ${\bar{m}}_{2}$ , $E_{2} I_{2}$ , $G_{2} A_{2}$ , ${\bar{m}}_{1}$ , $E_{1} I_{1}$ , $G_{1} A_{1}$ are the equivalent line mass, equivalent bending stiffness and equivalent shear stiffness of the superstructure and substructure respectively; $m_{i s o}$ , $K_{i s o}$ are the equivalent mass and equivalent shear stiffness of the isolation layer respectively. When the interlayer isolation structure is determined, all the above parameters can be determined.

Fig. 1Simplified equivalent mechanical model of interlayer isolation structure

Fig. 2Stress analysis of micro element of the beam

In the case of linear elasticity, for the equivalent mechanical model of interlayer isolation structure as shown in Fig. 1, take the micro segment $d x$ for stress analysis, as shown in Fig. 2, and the free vibration equation of the structure can be obtained as Eq. (1), which has the same form as Timoshenko Equation [19] of the beam with constant cross section:

1

E I \frac{\partial^{4} y (x, t)}{\partial x^{4}} = (\frac{E I \bar{m}}{G A} + \bar{J}) \frac{\partial^{4} y (x, t)}{\partial t^{2} \partial x^{2}} - \frac{\bar{J} \bar{m}}{G A} \frac{\partial^{4} y (x, t)}{\partial t^{4}} - \bar{m} \frac{\partial^{2} y (x, t)}{\partial t^{2}},

where, $A$ is the cross-sectional area of the beam, $l$ is the length of the beam, $E I$ and $G A$ are the bending stiffness and shear stiffness of the beam respectively, $\bar{m}$ are the line mass of the beam, and $\bar{J}$ is the moment of inertia of the beam to the neutral axis per unit length, $x$ is the length of the beam and $y$ is the vibration displacement of the beam. Neglect the influence of the moment of inertia of the section according to the existing research [29], which is means $\bar{J} =$ 0. Then Eq. (1) can be simplified and written as follows:

2

y^{''''} = \frac{\bar{m}}{G A} {\ddot{y}}^{''} - \frac{\bar{m}}{E I} \ddot{y},

where $y = (x, t)$ is the displacement function; $\dot{y}$ is the first partial derivative of displacement $y$ to time $t$ and $y'$ is the first partial derivative of displacement $y$ to coordinate $x$ .

3. Solution of structural natural vibration characteristics

The local coordinate system $x_{i} o_{i} y_{i}$ ( $i =$ 1, 2) and the overall coordinate system $x o y$ of the superstructure and substructure are established respectively, as shown in Fig. 1. The solution of Eq. (2) is expressed in the form of separation of variables:

3

y (x, t) = ϕ (x) s i n ω t,

where, $ω$ is the natural frequency of the structure and $ϕ (x)$ is the mode function of the structure.

Substituting Eq. (3) into Eq. (2), we can get:

4

ϕ'''' + \frac{\bar{m}}{G A} ω^{2} ϕ'' - \frac{\bar{m}}{E I} ω^{2} ϕ = 0 .

The parameter $a^{2} = \bar{m} / G A$ and $b^{2} = \bar{m} / E I$ are introduced, so that:

5

ϕ'''' + a^{2} ω^{2} ϕ'' - b^{2} ω^{2} ϕ = 0 .

The general solution of Eq. (5) is:

6

ϕ (x) = A_{1} s i n (S_{1} x) + A_{2} c o s (S_{1} x) + A_{3} s i n h (S_{2} x) + A_{4} c o s h (S_{2} x),

where, $A_{1}$ , $A_{2}$ , $A_{3}$ , $A_{4}$ are coefficients that determine the vibration amplitude of the structure and $S_{1} = \frac{a ω}{\sqrt{2}} \sqrt{\sqrt{1 + \frac{4 b^{2}}{a^{4} ω^{2}}} + 1}$ , $S_{2} = \frac{a ω}{\sqrt{2}} \sqrt{\sqrt{1 + \frac{4 b^{2}}{a^{4} ω^{2}}} - 1}$ .

In the local coordinate system, combined with Eq. (6) the boundary and continuity conditions of each substructure of the interlayer isolation structure are given respectively.

At the bottom of the structure, where $x_{1} =$ 0, the displacement $y$ and the angle $θ$ are both zero. The specific boundary conditions can be expressed as follows:

7

ϕ_{1} (0) = 0,

8

{ϕ_{1}}^{'''} (0) + (\frac{{b_{1}}^{2}}{{a_{1}}^{2}} + {a_{1}}^{2} {ω_{1}}^{2}) {ϕ_{1}}^{'} (0) = 0 .

At the top of the structure, where $x_{2} = l_{2}$ , the shear force $Q$ and the bending moment $M$ are both zero. The specific boundary conditions can be expressed as follows:

9

{ϕ_{2}}^{''} (l_{2}) + {a_{2}}^{2} {ω_{2}}^{2} ϕ_{2} (l_{2}) = 0,

10

{ϕ_{2}}^{'''} (l_{2}) + {a_{2}}^{2} {ω_{2}}^{2} {ϕ_{2}}^{'} (l_{2}) = 0 .

At the isolation layer, where $x_{1} = l_{1}$ and $x_{2} =$ 0, the displacement $y$ , the angle $θ$ , the shear force $Q$ and the bending moment $M$ are continuous. The specific boundary conditions can be expressed as follows:

11

ϕ_{1} (l_{1}) - \frac{E_{1} I_{1} [{ϕ_{1}}^{'''} (l_{1}) + {a_{1}}^{2} {ω_{1}}^{2} {ϕ_{1}}^{'} (l_{1})]}{K_{i s o}} = ϕ_{2} (0),

12

E_{1} I_{1} [{ϕ_{1}}^{'''} (l_{1}) + {a_{1}}^{2} {ω_{1}}^{2} {ϕ_{1}}^{'} (l_{1})] = E_{2} I_{2} [{ϕ_{2}}^{'''} (0) + {a_{2}}^{2} {ω_{2}}^{2} ϕ_{2}'' (0)] - m_{i s o} {ω_{2}}^{2} ϕ_{2} (0),

13

E_{1} I_{1} [{ϕ_{1}}^{''} (l_{1}) + {a_{1}}^{2} {ω_{1}}^{2} ϕ_{1} (l_{1})] = E_{2} I_{2} [{ϕ_{2}}^{''} (0) + {a_{2}}^{2} {ω_{2}}^{2} ϕ_{2} (0)] - J_{i s o} {ω_{2}}^{2} ({ϕ_{2}}^{'} (0) + \frac{{a_{2}}^{4}}{{b_{2}}^{2}} {ω_{2}}^{2} {ϕ_{2}}^{'} (0) + \frac{{a_{2}}^{2}}{{b_{2}}^{2}} {ϕ_{2}}^{'''} (0)),

14

{ϕ_{1}}^{'''} (l_{1}) + (\frac{{a_{1}}^{4}}{{b_{1}}^{2}} + {a_{1}}^{2} {ω_{1}}^{2}) {ϕ_{1}}^{'} (l_{1}) = {ϕ_{2}}^{'''} (0) + (\frac{{a_{2}}^{4}}{{b_{2}}^{2}} + {a_{2}}^{2} {ω_{2}}^{2}) {ϕ_{2}}^{'} (0) + \frac{1}{G_{2} A_{2} K_{i s o}} [{ϕ_{2}}^{''} (0) + {a_{2}}^{2} {ω_{2}}^{2} ϕ_{2} (0)] .

In Eq. (12-14), $ϕ_{1}$ and $ϕ_{2}$ are the modal functions of the substructure and superstructure respectively, and the solutions are as follows:

15

ϕ_{1} (x) = A_{11} s i n (S_{11} x) + A_{12} c o s (S_{11} x) + A_{13} s i n h (S_{12} x) + A_{14} c o s h (S_{12} x),

16

ϕ_{2} (x) = A_{21} s i n (S_{21} x) + A_{22} c o s (S_{21} x) + A_{23} s i n h (S_{22} x) + A_{24} c o s h (S_{22} x),

where:

S_{11} = \frac{a_{1} ω_{1}}{\sqrt{2}} \sqrt{\sqrt{1 + \frac{4 {b_{1}}^{2}}{{a_{1}}^{4} {ω_{1}}^{2}}} + 1}, S_{12} = \frac{a_{1} ω_{1}}{\sqrt{2}} \sqrt{\sqrt{1 + \frac{4 {b_{1}}^{2}}{{a_{1}}^{4} {ω_{1}}^{2}}} - 1},

S_{21} = \frac{a_{2} ω_{2}}{\sqrt{2}} \sqrt{\sqrt{1 + \frac{4 {b_{2}}^{2}}{{a_{2}}^{4} {ω_{2}}^{2}}} + 1}, S_{22} = \frac{a_{2} ω_{2}}{\sqrt{2}} \sqrt{\sqrt{1 + \frac{4 {b_{2}}^{2}}{{a_{2}}^{4} {ω_{2}}^{2}}} - 1} .

By substituting Eq. (15-16) into Eq. (7-14) respectively, eight equations about coefficient $A_{i j}$ ( $i =$ 1, 2; $j =$ 1, 2, 3, 4) are obtained. The above 8 equations s are combined into a system of equations and can be written in matrix form:

17

[D] [A] = [0],

where, $[A] = {[\begin{array}{l} A_{11} & A_{12} & A_{13} & A_{14} \begin{array}{l} A_{21} & A_{22} & A_{23} & A_{24} \end{array} \end{array}]}_{1 \times 8}^{T}$ , [ $D$ ] is determined coefficient matrix of equations about unknown [ $A$ ]. In order to make Eq. (17) have solutions, it is required that $|D| =$ 0, and the natural frequency of the structure can be obtained by solving it; if the frequency $ω$ is substituted into Eq. (17), the coefficient $A_{i j}$ ( $i =$ 1, 2; $j =$ 1, 2, 3, 4) can be obtained, that is, the mode function $ϕ (x)$ corresponding to $ω$ can be obtained, thus the mode of the structure can be obtained.

4. Seismic response solution

Betti’s law is applied to the above distributed parameter system model, that is, the work done by the inertial force of the $m$ -th mode on the $n$ -th mode displacement is equal to the work done by the inertial force of the $n$ -th mode on the $m$ -th mode displacement. It is easy to obtain the first orthogonal condition of interlayer isolation structure:

18

\begin{array}{l} \int_{0}^{l_{1}} ϕ_{1 m} \bar{m} ϕ_{1 n} d x + \int_{0}^{l_{2}} ϕ_{1 m} \bar{m} ϕ_{1 n} d x + ϕ_{2 m} (0) m_{i s o} ϕ_{2 n} (0) \\ + [ϕ_{2 m}' (0) + \frac{a^{4}}{b^{2}} {ω_{m}}^{2} ϕ_{2 m}' (0) + \frac{a^{2}}{b^{2}} ϕ_{2 m}''' (0)] \cdot J_{i s o} \end{array}

\cdot [ϕ_{2 n}' (0) + \frac{a^{4}}{b^{2}} {ω_{n}}^{2} ϕ_{2 n}' (0) + \frac{a^{2}}{b^{2}} ϕ_{2 n}''' (0)] = 0

The response of the $n$ -th order modes of the superstructure and substructure are respectively (the $m$ -th order modes are similar):

19

y_{1 n} (x, t) = ϕ_{1 n} (x) Y_{n} s i n ω_{n} t,

20

y_{2 n} (x, t) = ϕ_{2 n} (x) Y_{n} s i n ω_{n} t,

Eq. (19-20) are substituted into the equations of motion of the superstructure and substructure and Eq. (11-14), respectively:

21

E I (ϕ_{, 1 n}'''' (x) + a^{2} {ω_{n}}^{2} ϕ_{1 n}'' (x)) = {ω_{n}}^{2} \bar{m} ϕ_{1 n} (x)

22

E I (ϕ_{2 n}'''' (x) + a^{2} {ω_{n}}^{2} ϕ_{2 n}'' (x)) = {ω_{n}}^{2} \bar{m} ϕ_{2 n} (x),

23

E I (ϕ_{2 n}''' (0) + a^{2} {ω_{n}}^{2} ϕ_{2 n}' (0)) - E I (ϕ_{1 n}^{'''} (l_{1}) + a^{2} {ω_{n}}^{2} ϕ_{1 n}^{'} (l_{1})) = m_{i s o} {ω_{n}}^{2} ϕ_{2 n} (0),

24

E I ({ϕ_{2 n}}^{''} (0) + a^{2} {ω_{n}}^{2} ϕ_{2 n} (0)) - E I ({ϕ_{1 n}}^{''} (l_{1}) + a^{2} {ω_{n}}^{2} ϕ_{1 n} (l_{1})) = - J_{i s o} {ω_{n}}^{2} ({ϕ_{2 n}}^{'} (0) + \frac{a^{4}}{b^{2}} {ω_{n}}^{2} {ϕ_{2 n}}^{'} (0) + \frac{a^{2}}{b^{2}} {ϕ_{2 n}}^{'''} (0)) .

Operate Eq. (20-23) as follows: Eq. (20)× $ϕ_{1 m} (x)$ is integrated on [0, $l_{1}$ ]; Eq. (21)× $ϕ_{2 m} (x)$ is integrated on [0, $l_{2}$ ]; Eq. (22)× $ϕ_{2 m} (0)$ ; Eq. (23)× $(ϕ_{2 m}' (0) + \frac{a^{4}}{b^{2}} {ω_{m}}^{2} ϕ_{2 m}' (0) + \frac{a^{2}}{b^{2}} ϕ_{2 m}''' (0))$ . The above four equations are added into one, and the first orthogonal condition is substituted into it, then the second orthogonal condition of interlayer isolation structure is obtained:

25

\int_{0}^{l_{1}} ϕ_{1 m} (x) (ϕ_{1 n}'''' (x) + a^{2} {ω_{n}}^{2} ϕ_{1 n}'' (x)) d x

+ \int_{0}^{l_{2}} ϕ_{2 m} (x) (ϕ_{2 n}'''' (x) + a^{2} {ω_{n}}^{2} ϕ_{2 n}'' (x)) d x

+ ϕ_{2 m} (0) [(ϕ_{2 n}''' (0) + a^{2} {ω_{n}}^{2} ϕ_{2 n}' (0)) - (ϕ_{1 n}''' (l_{1}) + a^{2} {ω_{n}}^{2} ϕ_{1 n}' (l_{1}))]

+ (ϕ_{2 m}' (0) + \frac{a^{4}}{b^{2}} {ω_{m}}^{2} ϕ_{2 m}' (0) + \frac{a^{2}}{b^{2}} ϕ_{2 m}''' (0))

\times [(ϕ_{2 n}'' (0) + a^{2} {ω_{n}}^{2} ϕ_{2 n} (0)) - (ϕ_{1 n}'' (l_{1}) + a^{2} {ω_{n}}^{2} ϕ_{1 n} (l_{1}))] = 0 .

Under the action of earthquake, the equations of motion of the damping system of the superstructure and substructure and the isolation layer are as follows:

26

E_{1} I_{1} \frac{\partial^{4} y_{1} (x, t)}{\partial x^{4}} + C_{S 1} I_{1} \frac{\partial^{4}}{\partial x^{4}} (\frac{\partial y_{1} (x, t)}{\partial t}) + {\bar{m}}_{1} \frac{\partial^{2} y_{1} (x, t)}{\partial t^{2}} + C_{1} \frac{\partial y_{1} (x, t)}{\partial t} - \frac{E_{1} I_{1} {\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial^{4} y_{1} (x, t)}{\partial t^{2} \partial x^{2}} - C_{S 1} I_{1} \frac{{\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial}{\partial t} (\frac{\partial^{4} y_{1} (x, t)}{\partial t^{2} \partial x^{2}}) = - {\bar{m}}_{1} {\ddot{u}}_{g},

27

E_{2} I_{2} \frac{\partial^{4} y_{2} (x, t)}{\partial x^{4}} + C_{S 2} I_{2} \frac{\partial^{4}}{\partial x^{4}} (\frac{\partial y_{2} (x, t)}{\partial t}) + {\bar{m}}_{2} \frac{\partial^{2} y_{2} (x, t)}{\partial t^{2}} + C_{2} \frac{\partial y_{2} (x, t)}{\partial t} - \frac{E_{2} I_{2} {\bar{m}}_{2}}{G_{2} A_{2}} \frac{\partial^{4} y_{2} (x, t)}{\partial t^{2} \partial x^{2}} - C_{S 2} I_{2} \frac{{\bar{m}}_{2}}{G_{2} A_{2}} \frac{\partial}{\partial t} (\frac{\partial^{4} y_{2} (x, t)}{\partial t^{2} \partial x^{2}}) = - {\bar{m}}_{2} {\ddot{u}}_{g},

28

m_{i s o} \frac{\partial^{2} y_{2} (0, t)}{\partial t^{2}} + C_{i s o} \frac{\partial y_{2} (0, t)}{\partial t} - E_{1} I_{1} (\frac{\partial^{3} y_{1} (l_{1}, t)}{\partial x^{3}} - \frac{{\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial^{3} y_{1} (l_{1}, t)}{\partial x \partial t^{2}})

- C_{S 1} (l_{1}) I_{1} \frac{\partial}{\partial t} (\frac{\partial^{3} y_{1} (l_{1}, t)}{\partial x^{3}} - \frac{{\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial^{3} y_{1} (l_{1}, t)}{\partial x \partial t^{2}}) + E_{2} I_{2} (\frac{\partial^{3} y_{2} (0, t)}{\partial x^{3}} - \frac{{\bar{m}}_{2}}{G_{2} A_{2}} \frac{\partial^{3} y_{2} (0, t)}{\partial x \partial t^{2}})

+ C_{S 2} (0) I_{2} \frac{\partial}{\partial t} (\frac{\partial^{3} y_{2} (0, t)}{\partial x^{3}} - \frac{{\bar{m}}_{2}}{G_{2} A_{2}} \frac{\partial^{3} y_{2} (0, t)}{\partial x \partial t^{2}}) = - m_{i s o} {\ddot{u}}_{g},

29

J_{i s o} \frac{\partial^{2}}{\partial t^{2}} (\frac{\partial y_{2} (0, t)}{\partial t} - \frac{E_{2} I_{2} {\bar{m}}_{2}}{(G_{2} A_{2})^{2}} \frac{\partial^{3} y_{2} (0, t)}{\partial x \partial t^{2}} + \frac{E_{2} I_{2}}{G_{2} A_{2}} \frac{\partial^{3} y_{2} (0, t)}{\partial x^{3}})

+ C_{J i s o} \frac{\partial}{\partial t} (\frac{\partial y_{2} (0, t)}{\partial t} - \frac{E_{2} I_{2} {\bar{m}}_{2}}{(G_{2} A_{2})^{2}} \frac{\partial^{3} y_{2} (0, t)}{\partial x \partial t^{2}} + \frac{E_{2} I_{2}}{G_{2} A_{2}} \frac{\partial^{3} y_{2} (0, t)}{\partial x^{3}})

+ E_{1} I_{1} (\frac{\partial^{2} y_{1} (l_{1}, t)}{\partial x^{2}} - \frac{{\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial^{2} y_{1} (l_{1}, t)}{\partial t^{2}}) + C_{S 1} (l_{1}) I_{1} \frac{\partial}{\partial t} (\frac{\partial^{2} y_{1} (l_{1}, t)}{\partial x^{2}} - \frac{{\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial^{2} y_{1} (l_{1}, t)}{\partial t^{2}})

- E_{1} I_{1} (\frac{\partial^{2} y_{2} (0, t)}{\partial x^{2}} - \frac{{\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial^{2} y_{2} (0, t)}{\partial t^{2}})

- C_{S 2} (0) I_{1} \frac{\partial}{\partial t} (\frac{\partial^{2} y_{2} (0, t)}{\partial x^{2}} - \frac{{\bar{m}}_{1}}{G_{1} A_{1}} \frac{\partial^{2} y_{2} (0, t)}{\partial t^{2}}) = 0,

where, ${\ddot{u}}_{g}$ is the acceleration of ground motion, $C$ and $C_{S}$ are the damping coefficients of horizontal velocity and strain velocity respectively. If the superstructure and substructure have the same damping characteristics, then $C_{1} = C_{2} = C$ , $C_{S 1} = C_{S 2} = C_{S}$ . $C_{i s o}$ and $C_{J i s o}$ are damping coefficients of horizontal velocity and rotational velocity of isolation layer respectively. Assuming that the structural damping is expressed as Rayleigh damping [29], that is to say, the damping effect is directly proportional to the structural mass and stiffness, then: $C = α \bar{m}$ , $C_{S} = β E$ , $C_{i s o} = α m_{i s o}$ , $C_{J} = α J_{i s o} .$

Use the following formula to change the geometric coordinates into regular coordinates:

30

y (x, t) = \sum_{i = 1}^{\infty} ϕ_{i} (x) Y_{i} (t) .

Eq. (30) is substituted into Eq. (26-29), Because the amplitude $Y_{i} (t)$ satisfies the motion equations of the superstructure and substructure and the isolation layer with damping effect at the same time, the four formulas can be added. Using the first and second orthogonal conditions, the sum formula after treatment is simplified:

31

{\ddot{Y}}_{m} (t) + (α + β ω_{m}^{2}) {\dot{Y}}_{m} (t) + ω_{m}^{2} Y_{m} (t) = \frac{P_{m}}{M_{m}} {\ddot{u}}_{g},

where:

M_{m} = \int_{0}^{l_{1}} {\bar{m}}_{1} {ϕ^{2}}_{1 m} (x) d x + \int_{0}^{l_{2}} {\bar{m}}_{2} {ϕ^{2}}_{2 m} (x) d x + m_{i s o} ϕ_{_{2 m}}^{2} (0)

+ J_{i s o} {[ϕ_{2 m}^{'} (0) + \frac{E_{2} I_{2} {\bar{m}}_{2}}{(G_{2} A_{2})^{2}} ω_{m}^{2} ϕ_{2 m}^{'} (0) + \frac{E_{2} I_{2}}{G_{2} A_{2}} ϕ_{2 m}^{'''} (0)]}^{2},

P_{m} = - (\int_{0}^{l_{1}} {\bar{m}}_{1} ϕ_{1 m} (x) d x + \int_{0}^{l_{2}} {\bar{m}}_{2} ϕ_{2 m} (x) d x + m_{i s o} ϕ_{2 m} (0)) .

By introducing the damping ratio of the $m$ -th mode $ξ_{m} = \frac{α}{2 ω_{m}} + \frac{β ω_{m}}{2}$ , we can get:

32

{\ddot{Y}}_{m} (t) + 2 ξ_{m} ω_{m} {\dot{Y}}_{m} (t) + ω_{m}^{2} Y_{m} (t) = \frac{P_{m}}{M_{m}} {\ddot{u}}_{g} .

The above formula can be solved by selecting appropriate seismic excitation. Then, the solution results are substituted into Eq. (29), and the dynamic response of the superstructure and substructure of the interlayer isolation structure can be obtained by the mode superposition method. The dynamic response of the whole system is obtained by simply assembling the results of the superstructure and substructure.

5. Structural calculation example and model establishment

In order to verify the correctness of the distributed parameter system model and the corresponding solution method based on Timoshenko beam theory, a 10-story interlayer isolation structure is designed and taken as an example, as shown in Fig. 3. The semi analytical method (SAM), and the ETABS finite element software simulation (FESE) are separately used to solve the seismic response of the example. In the structural design of this paper, the seismic fortification intensity is 8 degrees, the site category is three, and the design earthquake group is the second group.

Fig. 3Model design of the interlayer isolation structure

a) Three dimensional diagram of the model

b) Component distribution of the model

The non-isolation floor of the structure is designed to have 10 floors in total, the bottom floor is 3.6 m high and the rest each floor is 3.0 m high; the floor slab is 120 mm thick, the column section is 600 mm×600 mm, the beam section is 650 mm×300 mm and 550 mm×300 mm respectively; the concrete strength grade of each component is C30. The isolation layer is located between the 3-4 floors of the structure, and the floor of the isolation layer is 180 mm thick. According to the relevant provisions of code for seismic design of buildings (GB50011-2010) [30], 20 LRB600 and LRB500 lead-rubber bearings are selected. The finite element model is shown in Fig. 3.

When the story isolation structure with 10 stories and one isolation story is designed in ETABS software, the floor mass, stiffness, damping ratio and other data of the structure can be derived through the software, and the corresponding linear mass and stiffness can be obtained by dividing the concentrated mass stiffness by the floor height. The distribution mass of the superstructure and substructure $\bar{m}$ is 7.54×10⁴ kg/m, the equivalent bending stiffness $E I$ is 1.05×10¹¹ N·m², and the equivalent shear stiffness $G A$ is 9.0×10¹⁰ N. The total mass of the isolation layer $m_{i s o}$ is 2.15×10⁵ kg, the total shear stiffness $K_{i s o}$ is 2.37×10⁷ N/m, and the equivalent moment of inertia of isolation layer $J_{i s o}$ is 12196 kg·m². The damping ratio of the superstructure and substructure is 0.05 and the isolation layer’s is 0.15. The seismic excitation is recorded by Hector Mine (1999) seismic station (referred to as Hector Mine wave). The peak amplitude of acceleration of seismic wave is adjusted to 0.7 m/s², as shown in Fig. 4.

Fig. 4The seismic wave of Hector Mine

6. Results and discussion

6.1. Structural natural vibration characteristics and response analysis

6.1.1. Mode analysis

Fig. 5 shows the structural vibration mode. It can be observed that there is a sudden change in the isolation layer in the vibration mode of each stage, and vibration modes of the rest part of the structure change evenly. That is to say, the deformation of the structure is mainly concentrated in the isolation layer, which conforms to the vibration law of the interlayer isolation structure. The main period of the isolated structure is 2.529 s, which is 2.31 times of the main period of the corresponding non-isolated structure of 1.097 s, indicating that the design of interlayer isolation effectively extends the period of the structure and has a good isolation effect.

6.1.2. Periodic comparative analysis

The semi-analytical method was used to calculate the structural natural vibration period, and the results were compared with the ETABS finite element simulation results.

Fig. 5First three mode shapes

a) 1st mode shape

b) 2nd mode shape

c) 3rd mode shape

Table 2Period comparison of interlayer isolation structure

Mode $T$ (s)	$T$ (s)
Mode $T$ (s)	FESE	SAM
1	2.529	2.5284
2	0.526	0.6826

It is proved that the method of considering the shear deformation and using the distributed parameter theory to establish the equivalent model of the interlayer isolation structure, and using the semi-analytical method in this paper to solve the structural natural vibration characteristics is correct and effective. The comparison of the first two natural vibration periods of the structure is shown in Table 2.

It is found that the second-order periodic error obtained by semi-analytical method is large. The main reason is that the structure is equivalent to a cantilever beam at the beginning of model establishment and the stiffness of the beam model is estimated by the stiffness of the combined columns, which results in the smaller equivalent bending stiffness and shear stiffness of the structure and makes the whole period of each order larger. In addition, the smaller damping equivalence of the model for the isolation layer will also make the periodic calculation result larger.

6.1.3. Comparison of top layer’s responses

The acceleration and displacement response of the top layer of the structure by SAM and FESE. Results are compared and analyzed, as shown in Fig. 6. It can be found that the change trends of acceleration and displacement time history curve are the same, and the curves are in good agreement as well. Through simple analysis, the absolute error of top acceleration peak value is 16.26 %, and the absolute error of displacement peak value is 8.49 %. The results further prove that the model and solution method are effective and have good accuracy.

The response of the top layer of the structure calculated by SAM is slightly enlarged in 5-8 s, and the displacement response is advanced after 18 s. The reason is that when SAM is used to solve the structural period, the second-order period has more errors than the first-order period, which is larger than the finite element simulation results, making the top layer response of the structure enlarged. When using the FESE method to calculate the structural response, more modes are considered in order to be closer to the actual situation, but the SAM does not consider higher-order modes due to the calculation efficiency. Thus the three times history curves are slightly different. But on the whole, the acceleration and displacement responses of the top layer obtained by this model are in good agreement with the finite element results.

Fig. 6Response comparison of the top layer

a) Acceleration response comparison

b) Displacement response comparison

6.1.4. Comparison of isolation layer’s responses

The acceleration and displacement response of the isolation layer of the structure by SAM and FESE are compared and analyzed as shown in Fig. 7. It can be found that the change trend of acceleration and displacement time history curve of the isolation layer is consistent, which is in good agreement.

After 12.5 s, the acceleration and displacement response of the isolation layer by SAM are amplified. In order to explain this, the frequency spectrum of the Hector Mine seismic wave is analyzed in Fig. 4. At 12.5 s, the seismic acceleration reaches the peak, and the components of the seismic wave with a period of 0.5 s account for a large proportion. The second-order period of the finite element simulation structure is 0.526 s, the period is close, the structure resonance occurs, the response increases, and the acceleration response increase is greater than the speed response increase. The second-order period calculated by the semi analytical method is 0.6826 s. Because of the error, the resonance phenomenon is not obvious. Compared with the response of the top layer of the structure, there is no response amplification phenomenon, which is considered to be due to the filtering effect of the isolation layer. The large deformation of the isolation layer results in the absorption of seismic energy, which leads to the reduction of the low frequency wave components and the reduction of the structural response. Therefore, the structure can have the function of isolation. Of course, the filtering effect of the isolation layer needs further verification.

Fig. 7Response comparison of the isolation layer

a) Acceleration response comparison

b) Displacement response comparison

Based on the above, it can be concluded that the distributed parameter system model of the interlayer isolation structure considering both shear deformation and bending deformation is effective, which can accurately reflect the natural vibration characteristics of the interlayer isolation structure, and can obtain more accurate structural time history response. By using the distributed parameter system and the algorithm in this paper, the response at any position of the structure can be obtained, which provides a more detailed reference for the engineering analysis of the structure, and enriches the theoretical system of the interlayer isolation structure.

6.2. Parameter analysis

6.2.1. Single parameter analysis

For the convenience of expression, order:

33

θ_{i} = \frac{X'}{X}, (i = X), μ_{j} = \frac{Y'}{Y}, (j = Y),

where, $X$ is the parameter before adjustment, and the value can be ${\bar{m}}_{2}$ , $E_{2} I_{2}$ , $G_{2} A_{2}$ , ${\bar{m}}_{1}$ , $E_{1} I_{1}$ , $G_{1} A_{1}$ , $m_{i s o}$ , $K_{i s o}$ ; $X'$ is the corresponding adjusted parameter; $θ$ is the adjusted proportion. $Y'$ is the response of the interlayer isolation structure corresponding to $X'$ , and the values include the peak acceleration response of the superstructure $A_{2 m a x}$ , the peak displacement response of the superstructure $D_{2 m a x}$ , the peak acceleration response of the substructure $A_{1 m a x}$ , the peak displacement response of the substructure $D_{1 m a x}$ , and the deformation of the isolation layer $D_{i s o}$ ; $Y$ is the response of the corresponding non-isolated structure; $μ$ is the corresponding response ratio of isolated structure and non-isolated structure.

(1) Parameters of the superstructure.

The distributed mass, stiffness and damping ratio of superstructure are adjusted respectively to study the response of the superstructure and substructure and the deformation of the isolation layer. The results are shown in Fig. 8.

Fig. 8Influence curves of the superstructure parameters on structure response

a) Distributed mass

b) Bending stiffness

c) Damping ratio

In Fig. 8(a), the adjustment ratio of the superstructure mass is 0.1-15. It is found that when the mass of the superstructure of the interlayer isolation structure is smaller, the ratio $μ$ is greater than 1, which means the structure has obvious response amplification phenomenon. And with the increase of the superstructure mass, the response ratio decreases and the isolation effect is enhanced. In Fig. 8(b), with the increase of the bending stiffness of the superstructure, the acceleration response peak ratios of the superstructure and the substructure increases and decreases, and the amplification phenomenon of the substructure response disappears. In Fig. 8(c), increasing the damping of the superstructure increases the displacement response peak and the acceleration response peak of the substructure, which may cause the response amplification of the superstructure and substructure at the same time. In a word, the displacement response of the superstructure is most sensitive to the change of the mass of the superstructure, and the amplification of the response of the substructure can be eliminated by adjusting the bending stiffness of the superstructure.

In addition, the deformation of the isolation layer can be reduced by increasing the damping ratio of the superstructure. In practical engineering application, the isolation layer is the weak layer of the interlayer isolation structure, and reducing its deformation will reduce the overturning moment of the superstructure, reduce the overall $p$ - $Δ$ effect of the structure, and make the structure safer. Thus the deformation of the isolation layer should be controlled.

(2) Parameters of the substructure.

Fig. 9Influence curves of the substructure parameters on structure response

a) Distributed mass

b) Bending stiffness

c) Damping ratio

The distributed mass, stiffness and damping ratio of the substructure are adjusted respectively and the results are shown in Fig. 9. It can be concluded that the acceleration response of the substructure is more sensitive to the change of the mass of the substructure. If the value is too large, the acceleration response of the substructure will be amplified. With the increase of the bending stiffness of the substructure, the acceleration response of the superstructure increases continuously. Increasing the damping ratio of the substructure will significantly reduce the acceleration response of the substructure, but increase the response of the superstructure. Therefore, increasing the damping ratio of the substructure is one of the effective methods to avoid the acceleration amplification of the substructure, and there is an optimal $θ_{d a m p 1}$ to make the best isolation effect.

(3) Parameters of the isolation layer.

The mass, stiffness and damping ratio of the isolation layer are adjusted respectively to study the response of the superstructure and substructure and the deformation of the isolation layer. The results are shown in Fig. 10. The conclusions are as follows: too large mass of isolation layer will enlarge the acceleration response of substructure, but has little effect on the deformation of isolation layer. The larger the shear stiffness of the isolation layer is, the smaller the deformation of the isolation layer is. However, the larger the stiffness of the isolation layer is, the larger the response of the substructure will be. In Fig. 10, increasing the damping ratio of the isolation layer can effectively reduce the response of the structure, enhance the isolation effect, and effectively reduce the deformation of the isolation layer as well, which make the structure safer. It explains that in the practical application of interlayer isolation structure, damping devices are often added to the isolation layer to increase the damping ratio of the isolation layer, so as to ensure the safety of the structure and optimize the isolation effect.

Fig. 10Influence curves of the isolation layer parameters on structure response

a) Distributed mass

b) Bending stiffness

c) Damping ratio

6.2.2. Multi parameter analysis

According to the previous analysis, the response of interlayer isolation structure is more sensitive to the change of mass and bending stiffness of substructures. The mass and stiffness of substructure are not isolated for a building whose shape or function is determined. For reinforced concrete frame structure, $E I$ is linear with $m$ when the shape is determined. Therefore, this section will study the seismic response of the structure with the mass and stiffness changing at the same time. Make:

34

λ_{i} = \frac{Z'}{Z}, (i = 1,2, \dots, 9), μ_{j} = \frac{Y'}{Y}, (j = Y),

where $Z$ is the total mass of the structure, $Z^{'}$ is the mass of the substructure after adjustment, $λ$ is the adjustment proportion, $μ$ is the corresponding response ratio of the isolated structure and non-isolated structure. The isolation layer is located between $i$ -th and ( $i + 1$ )-th layers of the structure.

With the parameter $λ$ varies from 0 to 1, the distributed parameter system model based on the FEM is used to calculate the structure responses when the isolation layer is located in different positions, and the corresponding comparison with the corresponding non-isolation structure is made to draw the response peak ratio curve and the isolation layer deformation curve.

(1) Lower isolation layer ( $i =$ 1, 2, 3).

Fig. 11Influence curves of the superstructure and the substructure mass on structural response with lower isolation layer

a) Distributed mass

b) Bending stiffness

c) Damping ratio

When $i =$ 1, the $λ_{1}$ - $μ$ curves and the $λ_{1}$ - $D_{i s o}$ are as Fig. 11, and the natural frequency ratio of the superstructure and substructure $R$ is 0.0178. With the increase of the ratio of substructure mass to total structure mass $λ_{1}$ , The $λ_{1}$ - $μ (A_{2 m a x})$ , $λ_{1}$ - $μ (D_{2 m a x})$ , and $λ_{1}$ - $μ (D_{1 m a x})$ curves first rise and then fall, presenting “^” shape, while the $λ_{1}$ - $μ (A_{1 m a x})$ curve first falls and then rises, presenting “v” shape. And the deformation curve of the isolation layer fluctuating rises and then falls. Therefore, there is a value of $λ_{1}$ to make the isolation effect optimal and ensure the safety of the structure. When $i =$ 2, 3, the curves change similarly.

(2) Middle isolation layer ( $i =$ 4, 5, 6).

Fig. 12Influence curves of the superstructure and the substructure mass on structural response with middle isolation layer

a) Distributed mass

b) Bending stiffness

c) Damping ratio

When $i =$ 4, the $λ_{4}$ - $μ$ curves and the $λ_{4}$ - $D_{i s o}$ are as Fig. 12, and the natural frequency ratio of the superstructure and substructure is 0.4900. With the increase of $λ_{4}$ , the $λ_{4}$ - $μ (D_{2 m a x})$ and $λ_{4}$ - $μ (D_{1 m a x})$ curves first rise and then fall, presenting “^” shape. But there is no obvious regular pattern of curve $λ_{4}$ - $μ (A_{2 m a x})$ and $λ_{4}$ - $μ (A_{1 m a x})$ . The $λ_{1}$ - $D_{i s o}$ curve still presents “^” shape broadly. Similarly, there is an optimal $λ_{4}$ and $λ_{4}$ , $λ_{5}$ as well.

(3) Upper isolation layer ( $i =$ 7, 8, 9).

Fig. 13Influence curves of the superstructure and the substructure mass on structural response with upper isolation layer

a) Distributed mass

b) Bending stiffness

c) Damping ratio

When the isolation layer is respectively placed on the 7th, 8th and 9th floor, the $λ$ - $μ$ curves and the $λ$ - $D_{i s o}$ are as Fig. 13. With the increase of the ratio of substructure mass to total structure mass $λ$ , the $λ$ - $μ (A_{1 m a x})$ curve first rises and then falls, presenting “^” shape, while the $λ$ - $μ (A_{2 m a x})$ , $λ$ - $μ (D_{1 m a x})$ and $λ$ - $μ (D_{2 m a x})$ , curves first fall and then rise, presenting “v” shape. In the upper isolated structures, by changing the mass of the substructure and the superstructure, the acceleration response of the substructure is always enlarged compared with that of the non-isolated structure, so that the substructure is in a more dangerous state. The $λ$ - $D_{i s o}$ curves still present “^” shape broadly.

It can be found that when the isolation layer is located in different positions of the structure, the response of the superstructure and substructure varies with the change of mass ratio $λ$ , and there will be an optimal $λ_{d}$ to make the isolation effect the best, as shown in Table 3. And the $λ_{d}$ – $i$ curve is as Fig. 14.

Table 3Comparison of optimal isolation effect when the isolation located differently

Isolation layer’s location		$R$	$λ_{d}$	$μ$				$D_{i s o}$ (mm)
Isolation layer’s location		$R$	$λ_{d}$	$D_{2 m a x}$	$A_{2 m a x}$	$D_{1 m a x}$	$A_{1 m a x}$	$D_{i s o}$ (mm)
Lower isolation layer	1	0.0178	0.15	0.6742	0.4606	0.4606	0.3807	21.40
	2	0.0756	0.20	0.7635	0.4475	0.5467	0.7635	20.34
	3	0.2090	0.28	0.6984	0.3995	0.4452	0.6984	21.72
Middle isolation layer	4	0.4900	0.46	0.5734	0.5734	0.3397	1.152	23.97
	5	1.0816	0.64	0.6457	0.6457	0.3575	0.9758	30.84
	6	2.4025	0.73	0.7099	0.4840	0.4840	1.902	37.59
Upper isolation layer	7	5.7600	0.82	1.0770	0.3713	0.5814	1.0770	37.27
	8	16.8100	0.82	0.8642	0.4783	0.4573	0.8642	33.33
	9	84.6400	0.85	0.8629	0.5535	0.4880	0.8629	30.40

The curve can be fitted as $λ_{d} =$ 0.1 $i$ +0.049. The fitting residual norm is 0.19442, which shows that there is a linear relationship between the optimal mass ratio and the location of isolation layer in regular isolated structure. This means that in the design of irregular structure, such as large chassis tower structure, we can refer to the optimal mass ratio of the corresponding regular structure to adjust the structural mass and story height, so as to make the isolation effect of the isolation structure better.

Fig. 14Curve of optimal mass ratio with changing of the position of isolation layer

When $λ$ is equal to $λ_{d}$ in different $i$ , the corresponding peak response ratio of the superstructure and substructure and the deformation of the isolation layer can be obtained, and draw them into curves as Fig. 15. In the form of lower isolation, the response ratio $μ$ is less than 1, and the deformation of isolation layer is small, so the isolation effect is the best; in the form of middle isolation, the acceleration response of the substructure is amplified, and the deformation of isolation layer is increased; in the form of upper isolation, the response ratio $μ$ close to 1 and the deformation of isolation layer is more than 30 mm. Therefore, for the isolation effect, the lower isolation is better than the upper isolation, and the middle isolation should be avoided as far as possible, and the upper isolation should be designed to avoid excessive deformation of the isolation layer.

Fig. 15Structural responses at different optimal mass ratio

a) The responses of the substructures

b) The deformation of the isolation layer

7. Conclusions

1) In this paper, based on certain assumptions, a distributed parameter system model considering both bending deformation and shear deformation is established, and its dynamic equation is derived. And the orthogonal conditions of mode shape are derived by Betti’s law, which is used to obtain the equivalent vibration equation in single mode. Finally, the calculation method of the response of the interlayer isolation structure is obtained by the mode superposition method.

2) Comparing the results of SAM and FESE, the error of the peak value of acceleration and displacement response of the natural vibration period and the top layer and isolation layer of the structure is less than 10 %, which can prove the correctness of the model and calculation method.

3) The existence of the isolation layer can effectively reduce the response of the superstructure, and the isolation effect is good; but the response of the substructure is amplified. With the increase of the mass and stiffness of the superstructure, the mass of the substructure and the damping of the isolation layer can effectively reduce the amplification. In regular interlayer isolation structure, the lower isolation form can get better isolation effect and the optimal frequency ratio of the superstructure and substructure increases with the rising of isolation layer position.

4) In order to further reduce the error, the future research will be based on this paper to optimize the parameters, simplify the calculation process, and use this method to study the impact of the isolation layer of the interlayer isolation structure on the superstructure and substructure, so that the method can be more effectively applied to research and practical engineering.

References

Yadao Qiuyan et al., “Precautions for construction of interlayer isolation structure,” (in Japanese), Construction, Vol. 390, pp. 72–77, 1998.

Search CrossRef
Murakami K., Ozaki H., and Teramoto T., “Design and analysis of a building with the middle-story isolation structural system,” in Proceedings of 12th WCEE, 2000.

Search CrossRef
Toshiyuki Sueoka, Shingo Torii, and Yasuhiro Tsuneki, “The application of response control design using middle-story isolation system to high-rise building,” in Proeeedings of 13th WCEE Vallcouver, Vol. 3457, Jan. 2004.

Search CrossRef
R. Villaverde, “Roof isolation system to reduce the seismic response of buildings: a preliminary assessment,” Earthquake Spectra, Vol. 14, No. 3, pp. 521–532, Aug. 1998, https://doi.org/10.1193/1.1586013

Publisher
R. Villaverde, “Aseismic roof isolation system: feasibility study with 13-story building,” Journal of Structural Engineering, Vol. 128, No. 2, pp. 188–196, Feb. 2002, https://doi.org/10.1061/(asce)0733-9445(2002)128:2(188)

Publisher
Roberto Villaverde and Gilberto Mosqueda, “Aseismic roof isolation system: analytic and shake table studies,” Earthquake Engineering and Structural Dynamics, Vol. 28, No. 3, pp. 217–234, Feb. 1999, https://doi.org/10.1002/(sici)1096-9845(199903)28:3

Publisher
Xu Zhonggen, Zhou Fulin, and Qiu Xiangquan, “Structural dynamic analysis of Shantou Museum,” (in Chinese), World Earthquake Engineering, Vol. 12, No. 2, pp. 33–36, 1996.

Search CrossRef
Jiang Chunqin, “Study on seismic behavior of isolated buildings with intermediate floors,” (in Chinese), National Taiwan University, Taiwan, 2007.

Search CrossRef
Fu-Lin Z. et al., “Research and application of seismic isolation system for building structures,” Journal of Architecture and Civil Engineering, Vol. 2, pp. 1–8, 2006.

Search CrossRef
K. Tian, W. Liu, D. Feng, and Q. Yang, “Dynamic characteristic analysis and shaking table test for a curved surface isolated structure,” Engineering Structures, Vol. 203, p. 109847, Jan. 2020, https://doi.org/10.1016/j.engstruct.2019.109847

Publisher
A. Di Matteo, C. Masnata, and A. Pirrotta, “Simplified analytical solution for the optimal design of tuned mass damper inerter for base isolated structures,” Mechanical Systems and Signal Processing, Vol. 134, p. 106337, Dec. 2019, https://doi.org/10.1016/j.ymssp.2019.106337

Publisher
Yin Chuanyin et al., “Design comparison of the seismic design of base-isolated RC frames using two Chinese codes,” (in Chinese), Engineering Mechanics, Vol. 36, No. 9, pp. 197–204, 2019.

Search CrossRef
D. Fang, Y. Du, C. Liu, and Y. Deng, “Analysis for displacement response characteristics of complex multi-layer base-isolated structure under near-field ground motion,” (in Chinese), Journal of Southwest Jiaotong University, Vol. 55, No. 1, Feb. 2020, https://doi.org/10.3969/j.issn.0258-2724.20181038

Publisher
Dang Yu and Tian Hongtu., “Analysis of calculation accuracy and improvement of real mode superposition response spectrum method of isolated structures,” (in Chinese), Building Structure, Vol. 49, No. 16, pp. 120–126, 2019.

Search CrossRef
Zhou Fulin, Zhang Ying, and Tan Ping, “Theoretical study on story isolation system,” (in Chinese), Journal of Civil Engineering, Vol. 42, No. 8, pp. 1–8, 2009.

Search CrossRef
Cao Jialiang, Shi Weixing, and Liu Wenguang, “Sensitivity of dynamic response and characteristics of isolated structure to calculation model parameters,” (in Chinese), Engineering Mechanics, Vol. 28, No. 8, pp. 168–176, 2011.

Publisher
He Dongdong, Miao Qisong, and Li Wenfeng, “Parameter optimization design of high position story isolated structure subjected to earthquake,” (in Chinese), Journal of Vibration Engineering, Vol. 26, No. 6, pp. 895–900, 2013.

Search CrossRef
Qi Ai, Lin Yunteng, and Zheng Guochen, “Research on working mechanism of interlayer isolation structure,” (in Chinese), Earthquake Engineering and Engineering Vibration, Vol. 2006, No. 4, pp. 239–243, 2006.

Search CrossRef
Wang Guangyuan, Vibration of Building Structure. (in Chinese), Beijing: Science Press, 1978.

Search CrossRef
Li Luxian et al., “Study on physical essence of the Timoshenko beam theory and its relation with the Euler beam theory,” (in Chinese), Chinese Science: Technology Science, Vol. 48, No. 4, pp. 360–368, 2018.

Publisher
Liu Liangkun et al., “The seismic response analysis of frame-shear wall structures based on Laplace transformation,” (in Chinese), Journal of Vibration Engineering, Vol. 30, No. 5, pp. 790–797, 2017.

Search CrossRef
Y.-C. Li, S.-Y. Liu, and Z.-G. Tang, “Two-time scale control of two manipulators handling flexible payload based on distributed parameter,” (in Chinese), Kongzhi yu Juece/Control and Decision, Vol. 27, No. 5, pp. 686–690+696, May 2012.

Search CrossRef
X. Cao and Y. Li, “Variable structure force control of flexible manipulator based on distributed parameter model,” (in Chinese), Control and Decision, Vol. 2008, No. 9, pp. 1035–1039, 2008.

Search CrossRef
Chen Yangyang et al., “A theoretical study on the horizontal dominant vibration mode of a girder bridge with tall piers,” (in Chinese), Vibration and Impact, Vol. 35, No. 12, pp. 13–21, 2016.

Search CrossRef
Chen Yangyang et al., “Semi-analytical method for time-history responses of single column and mass systems used to seismic isolated regular bridges,” (in Chinese), Earthquake Engineering and Engineering Vibration, Vol. 33, No. 3, pp. 183–191, 2013.

Search CrossRef
Du Yongfeng, Liu Yanhui, and Li Hui, “Semi-analytical solution of seismic response of serially isolated equipment and supporter,” (in Chinese), Journal of Mechanics, Vol. 41, No. 3, pp. 440–448, 2009.

Search CrossRef
Du Yongfeng, Liu Yanhui, and Li Hui, “Response of serial system of equipment and supporter under earthquake based on semi-analytical method,” (in Chinese), Journal of Computational Mechanics, Vol. 27, No. 2, pp. 225–231, 2010.

Search CrossRef
Song Xiao et al., “Analysis of parameters and aseismic properties of inter-story isolation system,” (in Chinese), Earthquake Engineering and Engineering Vibration, Vol. 38, No. 5, pp. 41–49, 2018.

Search CrossRef
R. Kraft and J. Peng Jin, Structural dynamics: Translated by Wang Guangyuan. (in Chinese), Beijing: Higher Education Press, 2006.

Search CrossRef
“Code for seismic design of buildings: GB50011-2010,” (in Chinese), China Architecture and Building Press, Beijing, 2010.

Search CrossRef

Cited by

Seismic response of a mid-story isolated stilted structure in mountainous areas

Jianhua Li | Dewen Liu | Ahad Javanmardi

(2024)

Study on the shock-absorbing effect of a new staggered story isolated structure under the long-period earthquake motion

(2023)

Broadband seismic isolation effect of three-dimensional seismic metamaterials in a semi-infinite foundation: modelled by fast multipole indirect boundary element method

Zhongxian Liu | Mingkai Zhang | Lei Huang | Hai Zhang

(2023)

About this article

Received

19 April 2021

Accepted

06 August 2021

Published

30 August 2021

SUBJECTS

Seismic engineering and applications

DOI

https://doi.org/10.21595/jve.2021.22003

Keywords

interlayer isolation structure

distributed parameter system

Timoshenko beam theory

semi analytical method

seismic response

Acknowledgements

The research is funded by the Project of Outstanding Youth Program of Natural Science Foundation of Ningxia [Grant No. 2021AAC05003] and Young Scholars in the West [Grant No. XAB2017AW09] and Natural Science Foundation of Ningxia [Grant No. 2019AAC03041] and The National Natural Science Foundation of China [Grant No. 51608283].

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous article in issue Previous Next article in issue Next

Research article

2020 06 30

Dynamic reliability sensitivity analysis for inter-story isolation structure under stochastic excitations

Shangrong Zhang, Fanglan Liu, Jianqiang Xiao, Yuchen Hu, Hairui Deng

Research article

2020 02 15

Performance-based seismic isolation design using the theory of spatially concave friction distribution

Shanshan Li, Biao Wei, Chaobin Li, Weihao Wang, Yunji Fu

Research article

2018 10 19

Mean square responses of a viscoelastic Timoshenko cantilever beam with different damping mechanisms

Qingzhao Zhou, David He, Zhongcheng Wu, Yaping Zhao

Research article

2018 09 30

Analytical solution for modal analysis of Euler-Bernoulli and Timoshenko beam with an arbitrary varying cross-section

Fatemeh Sohani, H. R. Eipakchi

S. Zhang, F. Liu, Y. Hu, S. Li, and L. Zhu, “Study on seismic response characteristics of the interlayer isolation structure,” Journal of Vibroengineering, Vol. 23, No. 8, pp. 1765–1784, Aug. 2021, https://doi.org/10.21595/jve.2021.22003

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2021.22003
UR  - https://doi.org/10.21595/jve.2021.22003
TI  - Study on seismic response characteristics of the interlayer isolation structure
T2  - Journal of Vibroengineering
AU  - Liu, Fanglan
AU  - Zhang, Shangrong
AU  - Hu, Yuchen
AU  - Li, Shihao
AU  - Zhu, Lei
PY  - 2021
DA  - 2021/08/30
PB  - JVE International Ltd.
SP  - 1765-1784
IS  - 8
VL  - 23
SN  - 1392-8716
SN  - 2538-8460
ER  - 

Copy Ris

Copied to clipboard!

@article{Liu_2021,
	doi = {10.21595/jve.2021.22003},
	url = {https://doi.org/10.21595/jve.2021.22003},
	year = 2021,
	month = {aug},
	publisher = {{JVE} International Ltd.},
	volume = {23},
	number = {8},
	pages = {1765--1784},
	author = {Fanglan Liu and Shangrong Zhang and Yuchen Hu and Shihao Li and Lei Zhu},
	title = {Study on seismic response characteristics of the interlayer isolation structure},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]F. Liu, S. Zhang, Y. Hu, S. Li, and L. Zhu, “Study on seismic response characteristics of the interlayer isolation structure,” Journal of Vibroengineering, vol. 23, no. 8, pp. 1765–1784, Aug. 2021, doi: 10.21595/jve.2021.22003.

Copy IEEE

Copied to clipboard!

Liu, Fanglan, Shangrong Zhang, Yuchen Hu, Shihao Li, and Lei Zhu. “Study on Seismic Response Characteristics of the Interlayer Isolation Structure.” Journal of Vibroengineering 23, no. 8 (August 30, 2021): 1765–84. https://doi.org/10.21595/jve.2021.22003.

Copy Chicago

Copied to clipboard!

Study on seismic response characteristics of the interlayer isolation structure

Abstract

Highlights

1. Introduction

2. Establishment of dynamic analysis model

2.1. Basic assumptions

2.2. Establishment of analysis model

3. Solution of structural natural vibration characteristics

4. Seismic response solution

5. Structural calculation example and model establishment

6. Results and discussion

6.1. Structural natural vibration characteristics and response analysis

6.1.1. Mode analysis

6.1.2. Periodic comparative analysis

6.1.3. Comparison of top layer’s responses

6.1.4. Comparison of isolation layer’s responses

6.2. Parameter analysis

6.2.1. Single parameter analysis

6.2.2. Multi parameter analysis

7. Conclusions

References

Cited by

About this article

Related Articles