# Rural Masonry Isolating Structure and LRB Experiment: Seismic Resistance Properties and Statistical Data Curve Analyses

Authors: Jiayao Zheng
Journal: Advances in Materials Science and Engineering
 (2020)

---

## Abstract

The low masonry structure is the most commonly applied building type in rural China. It is possible to install small-diameter, low cost, and easily constructed laminated rubber bearing (LRB) components. Isolation technology has broad application prospects in rural buildings. We developed a small-diameter LRB in this study wherein the isolation layer is set above the floor for easy installation and replacement. We built and tested 4 walls to observe the effects of different LRB thicknesses; we assessed test respective phenomena and seismic parameters accordingly. We ran another test on five small-diameter LRB components with varying horizontal stiffness, different forms of shear strain-equivalent horizontal stiffness, and postyield stiffness while changing the fitting formula for the second shape coefficient to give small-diameter LRB design providing gist.


---

## Body

## 1. Introduction

The isolation structure system serves to protect structural and nonstructural components in the building in which it is installed [1]. There is an urgent demand for new research on low masonry structure seismic building technology [2]. Rural buildings rarely contain laminated rubber bearing (LRB) isolation technologies [3]. Currently available seismic isolation techniques yet merit improvement to meet today’s demand for low masonry structures [4]. Previous researchers have added sliding isolated bearings [5] to rural masonry structures, but these bearings are limited to relatively few specific types of structures.In this study, we conducted an experiment consisting of four 1/2-ratio pseudodynamic seismic tests on masonry wall structures. We used one piece of a masonry wall with no isolation bearing (SRB) and three with different LRB specifications, completing 4 groups of masonry building bearing wall models to the seismic performance test. The equivalent horizontal stiffness, postyield stiffness, and damping ratio [6] of the isolation layer and bearing were determined throughout the experiment to gather workable data for subsequent mechanical analyses. We propose a method for the design of seismic isolation bearings accordingly [7]. The masonry wall (Figure 1(a)) was numbered “MLQT-1.” The other three walls using LRB all have a diameter of 120 mm but heights of 50 mm, 70 mm, and 90 mm, respectively (Figure 1(b)), in addition to the 90 mm LRB wall. These specimens are numbered MLQT-2, MLQT-3, and MLQT-4, respectively.Figure 1
Masonry wall specimens: (a) MLQT-1; (b) MLQT-4.
(a)(b)The bearing wall test loading device is shown in Figure2. According to the actual loading state of the wall, a loading device was utilized to impose a vertical uniform load as a voltage regulator was used to maintain stable axial force during the test [8]. A horizontal maximum load was applied to the top of the wall to simulate an El Centro wave peak acceleration. An electric hydraulic servo actuator [9] was used to apply different peaks of acceleration waves until the failure of the LRB or wall, respectively.Figure 2
Bearing wall experiment loading device: (a) photo; (b) drawing.
(a)(b)We also conducted an experiment to test various LRB mechanical parameters based on the “isolated rubber bears experimental method” [10]. Under these standards, the shear strain was 50%, 100%, and 250% corresponding to maximum displacement to calculate the LRB equivalent horizontal stiffness and damping ratio [11]. Based on the wall size of typical existing rural masonry structures, the design of experimental LRB is under 190 mm, has 12 rubber layers each 1.8～3.5 mm thick, and has a lead core of approximately 12 mm in diameter. We used five isolated bars all with an external diameter of 120 mm as lead core laminated rubber supports. The LRB is composed of a cover steel plate, thin steel plate, layer rubber, and lead core; all are cylindrical in shape. The LRB structure and specifications are listed in Table 1. Each type of LRB was tested in two replications.Table 1
Experimental bar specifications.
Bar modelRubber shear elasticity modulus (MPa)Effective diameter (mm)Number of layersRubber layer thickness (mm)Thin steel plate thickness (mm)Thin steel plate numberLead core diameter (mm)Cover steel plate thickness (mm)Bear overall height (mm)Second shape coefficient = dia/total thickness of the rubber layerGZ10.39110121.8211121063.65.1GZ20.39110122.52111310723.7GZ30.39110103.5291410733.1GZ40.39110121.9211101064.84.8GZ50.55110121.8211121063.65.1We conducted the horizontal stiffness experiment on these LRB specimens at the Guangzhou University Engineering Seismic Research Center on a 1000 t electro-hydraulic servo press and shear tester with the horizontal dynamic loading equipment shown in Figure3. The experiment horizontal loading was controlled by maximum displacement [12]. The vertical loading was 30 kN as the same type of two LRB specimens was inputting 0.1 Hz and 0.5 Hz sine waves separately with different values of peak acceleration in the horizontal direction. The maximum displacement changed as the shearing strain was set to 50%, 100%, and 250% [13]. Loading was recycled four times and the third hysteresis loop was calculated [14], where the LRB number end-term represents the sine wave frequency as GZ1-0.1 GZ1-0.5.Figure 3
Testing LRB and equipment: (a) dynamic loading equipment; (b) LRB.
(a)(b)We conducted five small-diameter LRB horizontal stiffness tests to observe the effects of different LRB horizontal shear strains and second shape coefficients [15] with different forms of equivalent horizontal stiffness, postyield stiffness, and hysteretic behavior-related mechanical properties [16]. We wrote a fitting formula for changes in shear strains with equivalent horizontal stiffness, postyield stiffness, and second shape coefficient. We also analyzed the relationship between the structural mode damping coefficient and horizontal displacement; this formed a correlation curve.The small-diameter LRB mechanical property parameters for rural dwellings provided a reference for our LRB test specimen design. The horizontal isolation layer increases the fundamental horizontal period of the structure significantly; this layer is what isolates the structural system from strong seismic disturbances in the ground. Structural deformation is concentrated in the isolation layer [17], so the LRB should be designed according to predetermined seismic aims and applicable targets to function properly.


## 2. Working Principles of Masonry Wall Experiment and LRB Design

### 2.1. Simulating El Centro Seismic Waves

We adopted pseudodynamic loading conditions for the purposes of this test. A horizontal simulating was applied to the top of the wall [18] with an El Centro seismic wave. The vertical load was applied to simulate the actual loading state of the wall through the vertical device, the dynamic load acceleration, the force of inertia, and the velocity damping force. The displacement of the elastic restoring force can be expressed as follows:(1)RIFd,v,a+RSSd,v,a=F,where RIF is the isolation layer counterforce, RSS is the superstructure counterforce, F is the external load, and d,v,anda are displacement, velocity, and acceleration, respectively. The constant equation of motion [19] is(2)Mai+Cvi+Kdi=F.The following assumption [20] was imposed for the i step velocity and acceleration:(3)vi=di+1−di−12Δt,(4)ai=di+1−2di+di−1Δt2,where M,C,andK are the mass, damping, and stiffness matrix, respectively. The i + 1 step displacement expression can be written by integrating equations (2)–(4).


### 2.2. LRB Design

#### 2.2.1. Shinnev et al.’s Method

The superstructure was fitted here with the LRB and considered as a rigid body. Displacement on one side moves along a similar trapezoid (first mode) as the oscillating displacement increases in equal proportion at the same height. The double bilinear isolation system can be simplified using this model [21]. The effective nature period of vibration can be determined via(5)TB=2πMKB,(6)ξB=ξb+ξh,(7)Tb=CbTB4πM,(8)ξh−2πAh4SbXb,where ξb is the effective velocity damping coefficient, KB is the secant stiffness of the LRB, Ah is the hysteretic curve area, and Sband Xb represent the effective restoring force and displacement.The equivalent seismic response in the bilinear [22] isolation system can be derived from the response spectrum of the earthquake, per the given seismic design:(9)Xb≈CFSDTB,ξBSb≈Qy+Kb2Xb−Xy,where Qy is the yield force and Kb2 is the second stiffness.


#### 2.2.2. AMER UBC Section Method

Assuming that seismic shear force [23] is discontinuous in the isolation system. Its isolated interface base shear formula is(10)Vb=kmaxD1.5,and the superstructure base shear formula is(11)VB=kmaxDRw1,D=10ZNS1T1B,T1=2πWKming,where Kmax is the maximum effective stiffness, RW1 is a ductility coefficient of the superstructure, and Z is the zone coefficient. B is a coefficient associated with material damping, N is a coefficient associated with the structure and the active fracture distance, S1 is the ground influence coefficient, W is the total weight of the structure above isolation interface, and Kmin is the isolated system minimum effective stiffness.The superstructure seismic force distribution is rectangular along its height when the base is isolated:(12)Fi=VBWi∑i=1nWi,Keff=F+−F−Δ+−Δ−,where Wi is the floor variate, n is the number of layers, Keff is the equivalent stiffness, and F−and Δ− are the minimum force/displacement in one cycle.


#### 2.2.3. Zhou Method

The superstructure is equivalent to a particle, as is the isolated layer [24]. The equivalent mass of the superstructure can be obtained when the two systems’ kinetic energy and base shear are considered to be equal.(13)Ms,eq=∑i=1Nmlxli2∑i=1Nmlxl2i.The equivalent stiffness is(14)Ks,eq=∑i=1NMlxli∑mlxl2i2∑i=1Nklxl2i,where N is the number of layers, mi is the layer I mass, ki is the layer I stiffness, and x1i is a fixed base fundamental mode.The period is(15)Tl≈Ts,l2+Tb2,T2=Ts,l1+μ1+Ts,l2Tb2,Ts,l=2πMs,eqKs,eq.The first model participation coefficient is(16)xsxb=1+Ts12Tb2γ1=1+μ1+Ts,121+μ1+2Ts,12/Tb2,and the second model participation coefficient [25] is(17)xbxs≈−μ+1+μTs,12Tb2γ2=Ts,12Tb21μ1+2Ts,12/Tb2,where μ=Ms.eq/Mb is the mass ratio, Tb=2πMb+Ms,eq/Kb is the hypothetical superstructure rigid body isolated structure period, Mb is the mass of the isolated layer, and Kb is the stiffness of the isolated layer.


#### 2.2.4. Proposed Method

We calculated the horizontal equivalent stiffness and equivalent viscous damping ratio in the isolated layer using an oscillation equation’s plural damping theory. The natural vibration period of a rural masonry structure is short, so the appropriate horizontal damping coefficient is 1.2 [25]. A simplified calculation gets the isolation layer horizontal stiffness and equivalent damping ration:(18)Kh=nKj,ςeq=ςj,where Kh is the isolated layer’s horizontal equivalent stiffness, Kj is each signal’s LRB horizontal equivalent stiffness, ςeq is the isolated layer’s equivalent viscous damping ratio, ςj is the LRB’s equivalent viscous damping ratio, and n is the number of bearings.The horizontal direction seismic reduction coefficient is expressed as(19)β=1.2η2TgT1γ,where β is the horizontal seismic reduction coefficient, η2 is the damping adjustment coefficient, γ is the curve descending section attenuation coefficient, Tg is the Eigen period, and T1 is the isolated structure essential nature period of vibration.(20)η2=1+0.05−ξeq0.08+1.6ξeq,γ=0.9+0.05−ξeq0.3+6ξeq,T1=2πGKhg.The LRB horizontal displacement under strong seismic conditions can be calculated as follows:(21)μi≤μi,where μi is the LRB’s horizontal displacement under strong seismic conditions; μi is the limit of displacement when the LRB effective diameter is multiplied by less than 0.55 and the inner rubber total thickness is multiplied threefold [26].(22)μe=λsαhGKh,μi=1.15μe,where μe is the isolation layer centroid of horizontal displacement in a “rare occurrence” earthquake scenario and λs is the near field coefficient.


### 2.3. LRB Mechanical Model

#### 2.3.1. Equivalent Linear Model

We used a linear stiffness and damping equivalent calculation method to assess LRB mechanical properties. The linear stiffness is secant in this case:(23)F=kμ.We analyzed the LRB damping force separately according to a unidirectional damping resilience-displacement relation:(24)F=cμ˙,(25)ξ=Wd4πWe=Wd2πκμ2,where ξ is the equivalent stiction damping ratio, Wd is the equivalent damping energy consumption area, We is the elastic resilience work area, k is the horizontal equivalent linear stiffness, F is the LRB’s horizontal resilience, and μ is the LRB’s horizontal displacement.Plugging the formulas (23) and (25) to (24), we get the following formula:(26)c=4πPξTg,where c is the equivalent stiction damping ratio, P is the LRB’s vertical loads, T is the isolated structure’s first natural period of vibration, and g is the gravitational acceleration.The linear elastic unit [27] comprised of stiffness (Kx, Ky, Kr), combination stiffness center (xK, yK) related to the isolation layer centroid (xB, yB), and stiffness eccentricity [28] is exB=xK−xB,eyB=yK−yB.(27)Fx=∑ikixμxB−eyBμrB,Fy=∑ikiyμyB−exBμrB,T=∑ikixyi2+kiyxi2μrB+∑ikiyexBeyB−∑ikixeyBexB,where kix,kiy,xi2,andyi2 represent the i-th LRB’s x- and y-direction shear stiffness and coordinates, respectively.(28)Fx=Cxμ˙xB−eyBμ˙rB,Fy=Cyμ˙yB−exBμ˙rB,T=Crμ˙rB+CyexBμ˙yB−CxeyBμ˙xB.


#### 2.3.2. Double Linear Model

There are three main types of the double linear model: the ideal elastic-plastic model, linear strengthen elastic-plastic model, and negative stiffness property elastic-plastic model [29]. We used the linear strengthen elastic-plastic model to unify force-displacement under analysis in this study.(29)Fb=Keμb,μb<μy,Fb=αKeμb+1−αFy,μb>μy,where α is the ratio of postyield stiffness to preyield stiffness, Fy is the yield force, μy is the yield displacement, Ke is the preyield stiffness, Fb is LRB resilience, and μb is LRB displacement.We took the double linear model property LRB for equivalent linear model simulation. The equivalent stiffness and equivalent stiction damping ratio in this case is(30)Keff=Fbμb,ξeff=2π1−αμb−μyμbFyFb.We assessed the LRB equivalent stiffness and equivalent stiction damping ratio under different horizontal displacement conditions. The double linear model force parameter conversion formula is(31)Qy=Fb1μb2−Fb2μb1μb2−μb1,Kp=Fb2−Fb1μb2−μb1,μy=Fb2μb2ξb2μb1−Fb1μb1ξb1μb2Fb2μb2ξb2−Fb1μb1ξb1,Fy=Fb1Fb2μb2ξb2−μb1ξb1Fb2μb2ξb2−Fb1μb1ξb1,Ke=Fb1Fb2μb2ξb2−μb1ξb1Fb2μb2ξb2μb1−Fb1μb1ξb1μb2,where Fb1,ub1,Kb1,ξb1,Fb2,ub2,Kb2,ξb2 reveal various displacement levels of LRB equivalent stiffness and equivalent stiction damping ratios. ξb1−Fb1,μb1,ξb2−Fb2,μb2 represent the equivalent stiction damping ratio.The LRB surface was modeled here using a Prager mobile induration model [30]. We hypothesized that ax,ay are the x, y-direction ratio of postyield shear stiffness and preyield shear stiffness, and that Fxy,Yx,Fyy,andYy are the x, y-direction uniaxial load yielding strength and displacement. The yield surface function is(32)S=Fx−Fx0Fxyn+Fy−Fy0Fyyn−1=0,where S is the yield function, Fx0,Fy0 are the x, y-direction loading yield surface center of inertia, and n is the loading yield surface curve exponent.The yield center mobile increases as follows:(33)dFx0dFY0=1Fx−Fx0/Fxy2+Fy−Fy0/Fyy2Fx−Fx0Fxy2Fx−Fx0Fy−Fy0Fyy2Fx−Fx0Fy−Fy0Fxy2Fy−Fy0Fyy2dFxdFy=Fx−Fx0Fxy2Fx−Fx0Fy−Fy0Fyy2Fx−Fx0Fy−Fy0Fxy2Fy−Fy0Fyy2dFxdFy,and the circle yield surface [31] is expressed as(34)duxduy=KexKey−1+1αx/1−αxKexFx−Fx0/Fxy2+αy/1−αyKeyFy−Fy0/Fyy2Fx−Fx0Fxy22Fx−Fx0Fxy2Fy−Fy0Fyy2Fy−Fy0Fyy2Fx−Fx0Fxy2Fy−Fy0Fyy22dFxdFy,Fx−Fx0Fxy2Kexdμx+Fy−Fy0Fyy2Keydμy>0,loading,=0,neutral,<0,unload.


#### 2.3.3. Triple Linear Model

A triple linear model can be adopted upon extensive deformation in the structure (i.e., stiffness hardening). The resilience-displacement relationship [32] in this case is(35)F=K1μ,μ≤D2,K1μ+K2−K1μ−D1,D1<μ≤D2,K1μ+K2−K1μ−D1+K3−K2μ−D2,μ>D2,where F is resilience, K1 is the shear stiffness when μ≤D1,K2 is the shear stiffness when D1<μ≤D2, K3 is the shear stiffness when μ>D2, and μ is the displacement.The stiction damping model [33] is(36)FD=∑i=1nF0i+Ciμ˙ρisgnμ˙,where F0i is 0 velocity damping, Ci is the damping coefficient, and ρi is a damping exponential. Coulomb damping ρi = 0, linear ratio damping ρi=1, and nonlinear damping ρi varies.The force-displacement relation [34] is(37)FD=Cisgnμ˙FD=Ciμ˙sgnμ˙FD=Ciμ˙2sgnμ˙.Our model is a combination of double and triple linear models, wherein we consider five LRB equivalent damping ratios with a total damping ratio of 10～20% and a high damping ratio of 15～20% shown in Table2. We observed a 250% hysteretic curve shape shown in Figure 4(b) that emerged on both sides of the curve and degenerated as displacement resilience increased.Table 2
LRB horizontal stiffness calculation results.
Bearing numberGZ1-0.5GZ2-0.5GZ3-0.5GZ4-0.5GZ5-0.5Equivalent horizontal stiffness (N/mm)278.47177.12156.13211.73328.51Equivalent damping ratio (%)10.617.2519.2419.0712.52Postyield stiffness (N/mm)236140205156270Yield stress9451160139014201280Figure 4
Supporter third hysteretic data: (a) 100% shear strain hysteretic loop; (b) 250% shear strain hysteretic loop.
(a)(b)


## 2.1. Simulating El Centro Seismic Waves

We adopted pseudodynamic loading conditions for the purposes of this test. A horizontal simulating was applied to the top of the wall [18] with an El Centro seismic wave. The vertical load was applied to simulate the actual loading state of the wall through the vertical device, the dynamic load acceleration, the force of inertia, and the velocity damping force. The displacement of the elastic restoring force can be expressed as follows:(1)RIFd,v,a+RSSd,v,a=F,where RIF is the isolation layer counterforce, RSS is the superstructure counterforce, F is the external load, and d,v,anda are displacement, velocity, and acceleration, respectively. The constant equation of motion [19] is(2)Mai+Cvi+Kdi=F.The following assumption [20] was imposed for the i step velocity and acceleration:(3)vi=di+1−di−12Δt,(4)ai=di+1−2di+di−1Δt2,where M,C,andK are the mass, damping, and stiffness matrix, respectively. The i + 1 step displacement expression can be written by integrating equations (2)–(4).


## 2.2. LRB Design

### 2.2.1. Shinnev et al.’s Method

The superstructure was fitted here with the LRB and considered as a rigid body. Displacement on one side moves along a similar trapezoid (first mode) as the oscillating displacement increases in equal proportion at the same height. The double bilinear isolation system can be simplified using this model [21]. The effective nature period of vibration can be determined via(5)TB=2πMKB,(6)ξB=ξb+ξh,(7)Tb=CbTB4πM,(8)ξh−2πAh4SbXb,where ξb is the effective velocity damping coefficient, KB is the secant stiffness of the LRB, Ah is the hysteretic curve area, and Sband Xb represent the effective restoring force and displacement.The equivalent seismic response in the bilinear [22] isolation system can be derived from the response spectrum of the earthquake, per the given seismic design:(9)Xb≈CFSDTB,ξBSb≈Qy+Kb2Xb−Xy,where Qy is the yield force and Kb2 is the second stiffness.


### 2.2.2. AMER UBC Section Method

Assuming that seismic shear force [23] is discontinuous in the isolation system. Its isolated interface base shear formula is(10)Vb=kmaxD1.5,and the superstructure base shear formula is(11)VB=kmaxDRw1,D=10ZNS1T1B,T1=2πWKming,where Kmax is the maximum effective stiffness, RW1 is a ductility coefficient of the superstructure, and Z is the zone coefficient. B is a coefficient associated with material damping, N is a coefficient associated with the structure and the active fracture distance, S1 is the ground influence coefficient, W is the total weight of the structure above isolation interface, and Kmin is the isolated system minimum effective stiffness.The superstructure seismic force distribution is rectangular along its height when the base is isolated:(12)Fi=VBWi∑i=1nWi,Keff=F+−F−Δ+−Δ−,where Wi is the floor variate, n is the number of layers, Keff is the equivalent stiffness, and F−and Δ− are the minimum force/displacement in one cycle.


### 2.2.3. Zhou Method

The superstructure is equivalent to a particle, as is the isolated layer [24]. The equivalent mass of the superstructure can be obtained when the two systems’ kinetic energy and base shear are considered to be equal.(13)Ms,eq=∑i=1Nmlxli2∑i=1Nmlxl2i.The equivalent stiffness is(14)Ks,eq=∑i=1NMlxli∑mlxl2i2∑i=1Nklxl2i,where N is the number of layers, mi is the layer I mass, ki is the layer I stiffness, and x1i is a fixed base fundamental mode.The period is(15)Tl≈Ts,l2+Tb2,T2=Ts,l1+μ1+Ts,l2Tb2,Ts,l=2πMs,eqKs,eq.The first model participation coefficient is(16)xsxb=1+Ts12Tb2γ1=1+μ1+Ts,121+μ1+2Ts,12/Tb2,and the second model participation coefficient [25] is(17)xbxs≈−μ+1+μTs,12Tb2γ2=Ts,12Tb21μ1+2Ts,12/Tb2,where μ=Ms.eq/Mb is the mass ratio, Tb=2πMb+Ms,eq/Kb is the hypothetical superstructure rigid body isolated structure period, Mb is the mass of the isolated layer, and Kb is the stiffness of the isolated layer.


### 2.2.4. Proposed Method

We calculated the horizontal equivalent stiffness and equivalent viscous damping ratio in the isolated layer using an oscillation equation’s plural damping theory. The natural vibration period of a rural masonry structure is short, so the appropriate horizontal damping coefficient is 1.2 [25]. A simplified calculation gets the isolation layer horizontal stiffness and equivalent damping ration:(18)Kh=nKj,ςeq=ςj,where Kh is the isolated layer’s horizontal equivalent stiffness, Kj is each signal’s LRB horizontal equivalent stiffness, ςeq is the isolated layer’s equivalent viscous damping ratio, ςj is the LRB’s equivalent viscous damping ratio, and n is the number of bearings.The horizontal direction seismic reduction coefficient is expressed as(19)β=1.2η2TgT1γ,where β is the horizontal seismic reduction coefficient, η2 is the damping adjustment coefficient, γ is the curve descending section attenuation coefficient, Tg is the Eigen period, and T1 is the isolated structure essential nature period of vibration.(20)η2=1+0.05−ξeq0.08+1.6ξeq,γ=0.9+0.05−ξeq0.3+6ξeq,T1=2πGKhg.The LRB horizontal displacement under strong seismic conditions can be calculated as follows:(21)μi≤μi,where μi is the LRB’s horizontal displacement under strong seismic conditions; μi is the limit of displacement when the LRB effective diameter is multiplied by less than 0.55 and the inner rubber total thickness is multiplied threefold [26].(22)μe=λsαhGKh,μi=1.15μe,where μe is the isolation layer centroid of horizontal displacement in a “rare occurrence” earthquake scenario and λs is the near field coefficient.


## 2.2.1. Shinnev et al.’s Method

The superstructure was fitted here with the LRB and considered as a rigid body. Displacement on one side moves along a similar trapezoid (first mode) as the oscillating displacement increases in equal proportion at the same height. The double bilinear isolation system can be simplified using this model [21]. The effective nature period of vibration can be determined via(5)TB=2πMKB,(6)ξB=ξb+ξh,(7)Tb=CbTB4πM,(8)ξh−2πAh4SbXb,where ξb is the effective velocity damping coefficient, KB is the secant stiffness of the LRB, Ah is the hysteretic curve area, and Sband Xb represent the effective restoring force and displacement.The equivalent seismic response in the bilinear [22] isolation system can be derived from the response spectrum of the earthquake, per the given seismic design:(9)Xb≈CFSDTB,ξBSb≈Qy+Kb2Xb−Xy,where Qy is the yield force and Kb2 is the second stiffness.


## 2.2.2. AMER UBC Section Method

Assuming that seismic shear force [23] is discontinuous in the isolation system. Its isolated interface base shear formula is(10)Vb=kmaxD1.5,and the superstructure base shear formula is(11)VB=kmaxDRw1,D=10ZNS1T1B,T1=2πWKming,where Kmax is the maximum effective stiffness, RW1 is a ductility coefficient of the superstructure, and Z is the zone coefficient. B is a coefficient associated with material damping, N is a coefficient associated with the structure and the active fracture distance, S1 is the ground influence coefficient, W is the total weight of the structure above isolation interface, and Kmin is the isolated system minimum effective stiffness.The superstructure seismic force distribution is rectangular along its height when the base is isolated:(12)Fi=VBWi∑i=1nWi,Keff=F+−F−Δ+−Δ−,where Wi is the floor variate, n is the number of layers, Keff is the equivalent stiffness, and F−and Δ− are the minimum force/displacement in one cycle.


## 2.2.3. Zhou Method

The superstructure is equivalent to a particle, as is the isolated layer [24]. The equivalent mass of the superstructure can be obtained when the two systems’ kinetic energy and base shear are considered to be equal.(13)Ms,eq=∑i=1Nmlxli2∑i=1Nmlxl2i.The equivalent stiffness is(14)Ks,eq=∑i=1NMlxli∑mlxl2i2∑i=1Nklxl2i,where N is the number of layers, mi is the layer I mass, ki is the layer I stiffness, and x1i is a fixed base fundamental mode.The period is(15)Tl≈Ts,l2+Tb2,T2=Ts,l1+μ1+Ts,l2Tb2,Ts,l=2πMs,eqKs,eq.The first model participation coefficient is(16)xsxb=1+Ts12Tb2γ1=1+μ1+Ts,121+μ1+2Ts,12/Tb2,and the second model participation coefficient [25] is(17)xbxs≈−μ+1+μTs,12Tb2γ2=Ts,12Tb21μ1+2Ts,12/Tb2,where μ=Ms.eq/Mb is the mass ratio, Tb=2πMb+Ms,eq/Kb is the hypothetical superstructure rigid body isolated structure period, Mb is the mass of the isolated layer, and Kb is the stiffness of the isolated layer.


## 2.2.4. Proposed Method

We calculated the horizontal equivalent stiffness and equivalent viscous damping ratio in the isolated layer using an oscillation equation’s plural damping theory. The natural vibration period of a rural masonry structure is short, so the appropriate horizontal damping coefficient is 1.2 [25]. A simplified calculation gets the isolation layer horizontal stiffness and equivalent damping ration:(18)Kh=nKj,ςeq=ςj,where Kh is the isolated layer’s horizontal equivalent stiffness, Kj is each signal’s LRB horizontal equivalent stiffness, ςeq is the isolated layer’s equivalent viscous damping ratio, ςj is the LRB’s equivalent viscous damping ratio, and n is the number of bearings.The horizontal direction seismic reduction coefficient is expressed as(19)β=1.2η2TgT1γ,where β is the horizontal seismic reduction coefficient, η2 is the damping adjustment coefficient, γ is the curve descending section attenuation coefficient, Tg is the Eigen period, and T1 is the isolated structure essential nature period of vibration.(20)η2=1+0.05−ξeq0.08+1.6ξeq,γ=0.9+0.05−ξeq0.3+6ξeq,T1=2πGKhg.The LRB horizontal displacement under strong seismic conditions can be calculated as follows:(21)μi≤μi,where μi is the LRB’s horizontal displacement under strong seismic conditions; μi is the limit of displacement when the LRB effective diameter is multiplied by less than 0.55 and the inner rubber total thickness is multiplied threefold [26].(22)μe=λsαhGKh,μi=1.15μe,where μe is the isolation layer centroid of horizontal displacement in a “rare occurrence” earthquake scenario and λs is the near field coefficient.


## 2.3. LRB Mechanical Model

### 2.3.1. Equivalent Linear Model

We used a linear stiffness and damping equivalent calculation method to assess LRB mechanical properties. The linear stiffness is secant in this case:(23)F=kμ.We analyzed the LRB damping force separately according to a unidirectional damping resilience-displacement relation:(24)F=cμ˙,(25)ξ=Wd4πWe=Wd2πκμ2,where ξ is the equivalent stiction damping ratio, Wd is the equivalent damping energy consumption area, We is the elastic resilience work area, k is the horizontal equivalent linear stiffness, F is the LRB’s horizontal resilience, and μ is the LRB’s horizontal displacement.Plugging the formulas (23) and (25) to (24), we get the following formula:(26)c=4πPξTg,where c is the equivalent stiction damping ratio, P is the LRB’s vertical loads, T is the isolated structure’s first natural period of vibration, and g is the gravitational acceleration.The linear elastic unit [27] comprised of stiffness (Kx, Ky, Kr), combination stiffness center (xK, yK) related to the isolation layer centroid (xB, yB), and stiffness eccentricity [28] is exB=xK−xB,eyB=yK−yB.(27)Fx=∑ikixμxB−eyBμrB,Fy=∑ikiyμyB−exBμrB,T=∑ikixyi2+kiyxi2μrB+∑ikiyexBeyB−∑ikixeyBexB,where kix,kiy,xi2,andyi2 represent the i-th LRB’s x- and y-direction shear stiffness and coordinates, respectively.(28)Fx=Cxμ˙xB−eyBμ˙rB,Fy=Cyμ˙yB−exBμ˙rB,T=Crμ˙rB+CyexBμ˙yB−CxeyBμ˙xB.


### 2.3.2. Double Linear Model

There are three main types of the double linear model: the ideal elastic-plastic model, linear strengthen elastic-plastic model, and negative stiffness property elastic-plastic model [29]. We used the linear strengthen elastic-plastic model to unify force-displacement under analysis in this study.(29)Fb=Keμb,μb<μy,Fb=αKeμb+1−αFy,μb>μy,where α is the ratio of postyield stiffness to preyield stiffness, Fy is the yield force, μy is the yield displacement, Ke is the preyield stiffness, Fb is LRB resilience, and μb is LRB displacement.We took the double linear model property LRB for equivalent linear model simulation. The equivalent stiffness and equivalent stiction damping ratio in this case is(30)Keff=Fbμb,ξeff=2π1−αμb−μyμbFyFb.We assessed the LRB equivalent stiffness and equivalent stiction damping ratio under different horizontal displacement conditions. The double linear model force parameter conversion formula is(31)Qy=Fb1μb2−Fb2μb1μb2−μb1,Kp=Fb2−Fb1μb2−μb1,μy=Fb2μb2ξb2μb1−Fb1μb1ξb1μb2Fb2μb2ξb2−Fb1μb1ξb1,Fy=Fb1Fb2μb2ξb2−μb1ξb1Fb2μb2ξb2−Fb1μb1ξb1,Ke=Fb1Fb2μb2ξb2−μb1ξb1Fb2μb2ξb2μb1−Fb1μb1ξb1μb2,where Fb1,ub1,Kb1,ξb1,Fb2,ub2,Kb2,ξb2 reveal various displacement levels of LRB equivalent stiffness and equivalent stiction damping ratios. ξb1−Fb1,μb1,ξb2−Fb2,μb2 represent the equivalent stiction damping ratio.The LRB surface was modeled here using a Prager mobile induration model [30]. We hypothesized that ax,ay are the x, y-direction ratio of postyield shear stiffness and preyield shear stiffness, and that Fxy,Yx,Fyy,andYy are the x, y-direction uniaxial load yielding strength and displacement. The yield surface function is(32)S=Fx−Fx0Fxyn+Fy−Fy0Fyyn−1=0,where S is the yield function, Fx0,Fy0 are the x, y-direction loading yield surface center of inertia, and n is the loading yield surface curve exponent.The yield center mobile increases as follows:(33)dFx0dFY0=1Fx−Fx0/Fxy2+Fy−Fy0/Fyy2Fx−Fx0Fxy2Fx−Fx0Fy−Fy0Fyy2Fx−Fx0Fy−Fy0Fxy2Fy−Fy0Fyy2dFxdFy=Fx−Fx0Fxy2Fx−Fx0Fy−Fy0Fyy2Fx−Fx0Fy−Fy0Fxy2Fy−Fy0Fyy2dFxdFy,and the circle yield surface [31] is expressed as(34)duxduy=KexKey−1+1αx/1−αxKexFx−Fx0/Fxy2+αy/1−αyKeyFy−Fy0/Fyy2Fx−Fx0Fxy22Fx−Fx0Fxy2Fy−Fy0Fyy2Fy−Fy0Fyy2Fx−Fx0Fxy2Fy−Fy0Fyy22dFxdFy,Fx−Fx0Fxy2Kexdμx+Fy−Fy0Fyy2Keydμy>0,loading,=0,neutral,<0,unload.


### 2.3.3. Triple Linear Model

A triple linear model can be adopted upon extensive deformation in the structure (i.e., stiffness hardening). The resilience-displacement relationship [32] in this case is(35)F=K1μ,μ≤D2,K1μ+K2−K1μ−D1,D1<μ≤D2,K1μ+K2−K1μ−D1+K3−K2μ−D2,μ>D2,where F is resilience, K1 is the shear stiffness when μ≤D1,K2 is the shear stiffness when D1<μ≤D2, K3 is the shear stiffness when μ>D2, and μ is the displacement.The stiction damping model [33] is(36)FD=∑i=1nF0i+Ciμ˙ρisgnμ˙,where F0i is 0 velocity damping, Ci is the damping coefficient, and ρi is a damping exponential. Coulomb damping ρi = 0, linear ratio damping ρi=1, and nonlinear damping ρi varies.The force-displacement relation [34] is(37)FD=Cisgnμ˙FD=Ciμ˙sgnμ˙FD=Ciμ˙2sgnμ˙.Our model is a combination of double and triple linear models, wherein we consider five LRB equivalent damping ratios with a total damping ratio of 10～20% and a high damping ratio of 15～20% shown in Table2. We observed a 250% hysteretic curve shape shown in Figure 4(b) that emerged on both sides of the curve and degenerated as displacement resilience increased.Table 2
LRB horizontal stiffness calculation results.
Bearing numberGZ1-0.5GZ2-0.5GZ3-0.5GZ4-0.5GZ5-0.5Equivalent horizontal stiffness (N/mm)278.47177.12156.13211.73328.51Equivalent damping ratio (%)10.617.2519.2419.0712.52Postyield stiffness (N/mm)236140205156270Yield stress9451160139014201280Figure 4
Supporter third hysteretic data: (a) 100% shear strain hysteretic loop; (b) 250% shear strain hysteretic loop.
(a)(b)


## 2.3.1. Equivalent Linear Model

We used a linear stiffness and damping equivalent calculation method to assess LRB mechanical properties. The linear stiffness is secant in this case:(23)F=kμ.We analyzed the LRB damping force separately according to a unidirectional damping resilience-displacement relation:(24)F=cμ˙,(25)ξ=Wd4πWe=Wd2πκμ2,where ξ is the equivalent stiction damping ratio, Wd is the equivalent damping energy consumption area, We is the elastic resilience work area, k is the horizontal equivalent linear stiffness, F is the LRB’s horizontal resilience, and μ is the LRB’s horizontal displacement.Plugging the formulas (23) and (25) to (24), we get the following formula:(26)c=4πPξTg,where c is the equivalent stiction damping ratio, P is the LRB’s vertical loads, T is the isolated structure’s first natural period of vibration, and g is the gravitational acceleration.The linear elastic unit [27] comprised of stiffness (Kx, Ky, Kr), combination stiffness center (xK, yK) related to the isolation layer centroid (xB, yB), and stiffness eccentricity [28] is exB=xK−xB,eyB=yK−yB.(27)Fx=∑ikixμxB−eyBμrB,Fy=∑ikiyμyB−exBμrB,T=∑ikixyi2+kiyxi2μrB+∑ikiyexBeyB−∑ikixeyBexB,where kix,kiy,xi2,andyi2 represent the i-th LRB’s x- and y-direction shear stiffness and coordinates, respectively.(28)Fx=Cxμ˙xB−eyBμ˙rB,Fy=Cyμ˙yB−exBμ˙rB,T=Crμ˙rB+CyexBμ˙yB−CxeyBμ˙xB.


## 2.3.2. Double Linear Model

There are three main types of the double linear model: the ideal elastic-plastic model, linear strengthen elastic-plastic model, and negative stiffness property elastic-plastic model [29]. We used the linear strengthen elastic-plastic model to unify force-displacement under analysis in this study.(29)Fb=Keμb,μb<μy,Fb=αKeμb+1−αFy,μb>μy,where α is the ratio of postyield stiffness to preyield stiffness, Fy is the yield force, μy is the yield displacement, Ke is the preyield stiffness, Fb is LRB resilience, and μb is LRB displacement.We took the double linear model property LRB for equivalent linear model simulation. The equivalent stiffness and equivalent stiction damping ratio in this case is(30)Keff=Fbμb,ξeff=2π1−αμb−μyμbFyFb.We assessed the LRB equivalent stiffness and equivalent stiction damping ratio under different horizontal displacement conditions. The double linear model force parameter conversion formula is(31)Qy=Fb1μb2−Fb2μb1μb2−μb1,Kp=Fb2−Fb1μb2−μb1,μy=Fb2μb2ξb2μb1−Fb1μb1ξb1μb2Fb2μb2ξb2−Fb1μb1ξb1,Fy=Fb1Fb2μb2ξb2−μb1ξb1Fb2μb2ξb2−Fb1μb1ξb1,Ke=Fb1Fb2μb2ξb2−μb1ξb1Fb2μb2ξb2μb1−Fb1μb1ξb1μb2,where Fb1,ub1,Kb1,ξb1,Fb2,ub2,Kb2,ξb2 reveal various displacement levels of LRB equivalent stiffness and equivalent stiction damping ratios. ξb1−Fb1,μb1,ξb2−Fb2,μb2 represent the equivalent stiction damping ratio.The LRB surface was modeled here using a Prager mobile induration model [30]. We hypothesized that ax,ay are the x, y-direction ratio of postyield shear stiffness and preyield shear stiffness, and that Fxy,Yx,Fyy,andYy are the x, y-direction uniaxial load yielding strength and displacement. The yield surface function is(32)S=Fx−Fx0Fxyn+Fy−Fy0Fyyn−1=0,where S is the yield function, Fx0,Fy0 are the x, y-direction loading yield surface center of inertia, and n is the loading yield surface curve exponent.The yield center mobile increases as follows:(33)dFx0dFY0=1Fx−Fx0/Fxy2+Fy−Fy0/Fyy2Fx−Fx0Fxy2Fx−Fx0Fy−Fy0Fyy2Fx−Fx0Fy−Fy0Fxy2Fy−Fy0Fyy2dFxdFy=Fx−Fx0Fxy2Fx−Fx0Fy−Fy0Fyy2Fx−Fx0Fy−Fy0Fxy2Fy−Fy0Fyy2dFxdFy,and the circle yield surface [31] is expressed as(34)duxduy=KexKey−1+1αx/1−αxKexFx−Fx0/Fxy2+αy/1−αyKeyFy−Fy0/Fyy2Fx−Fx0Fxy22Fx−Fx0Fxy2Fy−Fy0Fyy2Fy−Fy0Fyy2Fx−Fx0Fxy2Fy−Fy0Fyy22dFxdFy,Fx−Fx0Fxy2Kexdμx+Fy−Fy0Fyy2Keydμy>0,loading,=0,neutral,<0,unload.


## 2.3.3. Triple Linear Model

A triple linear model can be adopted upon extensive deformation in the structure (i.e., stiffness hardening). The resilience-displacement relationship [32] in this case is(35)F=K1μ,μ≤D2,K1μ+K2−K1μ−D1,D1<μ≤D2,K1μ+K2−K1μ−D1+K3−K2μ−D2,μ>D2,where F is resilience, K1 is the shear stiffness when μ≤D1,K2 is the shear stiffness when D1<μ≤D2, K3 is the shear stiffness when μ>D2, and μ is the displacement.The stiction damping model [33] is(36)FD=∑i=1nF0i+Ciμ˙ρisgnμ˙,where F0i is 0 velocity damping, Ci is the damping coefficient, and ρi is a damping exponential. Coulomb damping ρi = 0, linear ratio damping ρi=1, and nonlinear damping ρi varies.The force-displacement relation [34] is(37)FD=Cisgnμ˙FD=Ciμ˙sgnμ˙FD=Ciμ˙2sgnμ˙.Our model is a combination of double and triple linear models, wherein we consider five LRB equivalent damping ratios with a total damping ratio of 10～20% and a high damping ratio of 15～20% shown in Table2. We observed a 250% hysteretic curve shape shown in Figure 4(b) that emerged on both sides of the curve and degenerated as displacement resilience increased.Table 2
LRB horizontal stiffness calculation results.
Bearing numberGZ1-0.5GZ2-0.5GZ3-0.5GZ4-0.5GZ5-0.5Equivalent horizontal stiffness (N/mm)278.47177.12156.13211.73328.51Equivalent damping ratio (%)10.617.2519.2419.0712.52Postyield stiffness (N/mm)236140205156270Yield stress9451160139014201280Figure 4
Supporter third hysteretic data: (a) 100% shear strain hysteretic loop; (b) 250% shear strain hysteretic loop.
(a)(b)


## 3. Wall Experiment and Seismic Behavior

### 3.1. Relevant Phenomena

Our specimen MLQT-1 is an uninstalled LRB wall. The corresponding relation between earthquake intensity and peak acceleration of seismic waves was imposed under 70, 100, 300, and 400 gal El Centro wave excitation models. We found no obvious cracks until reaching displacement of 9.6 mm and horizontal force of +137.879 kN, at which point the corresponding peak acceleration of seismic waves was 300 gal. We continued reverse loading until finding oblique cracks in the wall as shown in Figure5; these shear cracks developed in the wall during pushing and pulling in a typical “X” shape [35], at which point the wall presented shear failure.Figure 5
MLQT-1 collapse state.MLQT-2 was equipped with a 50 mm LRB masonry specimen and input with a damping ratio of 0.15 and 70, 100, 300, 400, and 510 gal El Centro wave excitation models [36]. We calculated the LRB’s ultimate displacement according to the upper and lower beam relative displacement. When the acceleration peak reached 300 gal, the maximum displacement was 13.6 mm (exceeding LRB allowable limit of displacement 10.8 mm). The middle of the LRB was broken when the acceleration peak reached 510 gal; the wall was not damaged until the test was stopped, as shown in Figure 6.Figure 6
MLQT-2: (a) LRB limit of displacement state; (b) LRB broken state.
(a)(b)As shown in Figure6, MLQT-2 does not meet the design standard though the LRB and upper structure were not broken in the test. LRB thickness variations, in this case, were −1 mm, −2 mm, and −1 mm, which indicate that the wall was not compressed by the vertical load.Specimen MLQT-3 was equipped with a 57 mm LRB masonry specimen designed to be resistant to seismic waves of 400 gal peak acceleration. Once the peak acceleration reached 400 gal, the maximum displacement was 14.2 mm, i.e., below the allowable limit of 22.5 mm. The contact between the LRB and steel plate separated only slightly, which did not affect the continuation of the experiment. The maximum displacement was observed upon reaching 620 gal, at 26.4 mm, which does exceed the allowable limit. The specimen was considered to be broken at this time. The LRB internal steel plate was broken and fully exposed upon reaching 800 gal (Figure7), though the wall did not crack by the end of the experiment.Figure 7
MLQT-3: (a) LRB limit of displacement state; (b) LRB broken state.
(a)(b)MLQT-4 was equipped with a 90 mm LRB masonry specimen and designed to resist the peak acceleration of an 800 gal seismic wave [37]. The maximum displacement reached 39.1 mm at the peak acceleration of 800 gal, which is below the allowable limit of 60.5 mm; thus, this specimen met the predesign specifications. Upon reaching 1,240 gal, the maximum displacement was 65.7 mm (exceeding the allowable limit of 60.5 mm) and the specimen was considered to be broken (Figure 8). Upon reaching 1,400 gal, the maximum displacement of the model structure was 76.9 mm (push), the maximum base shear was 17 kN (push), the LRB was broken, and the wall showed no cracks.Figure 8
MLQT-4: (a) LRB limit of displacement state; (b) LRB broken state.
(a)(b)As shown in Figure8, the MLQT-4 LRB thickness variations after the vertical loads were applied were −1 mm, −1 mm, and −2 mm. The maximum values of the additional loads were, respectively, 187.3 kN and 184.5 kN. This difference is in the allowable range and the wall was not compressed.


### 3.2. Data Operating of Maximum Displacement and Shear

We found that MLQT-1 280–400 gal showed a maximum displacement difference of −3 mm by −4 minus −1 mm, −7 minus −4 mm; the value changed suddenly and dramatically by the end of the test and diagonal cracks appeared through the wall. MLQT-2 reached a maximum displacement exceeding the allowable limit at 300 gal, as listed in Table3.Table 3
MLQT-1\MLQT-2 displacement.
Specimen numberPeak acceleration (gal)Maximum displacement (mm)MLQT-1140−1.879200−4.280280−7.2514008.779MLQT-22009.99030013.61040016.65751021.497The MLQT-3 allowable displacement is 22.5 mm [38], which was exceeded at 620 gal. The MLQT-4 allowable displacement is 60.05 mm, which was exceeded at 1,240 gal (Table 4).Table 4
MLQT-3\MLQT-4 displacement.
Specimen numberPeak acceleration (gal)Maximum displacementMLQT-340014.231 mm (allowable disp 22.50 mm)26.421 mm620800MLQT-4100050.943 mm (allowable disp 60.05 mm)65.780 mm12401400MLQT-1, from 35 to 70 gal, applied acceleration increase of 100%, and then max base shear increased by 130%. For an increased acceleration of 40%, from 70 to 100 gal, the max base shear increased by 34% (Table5). From 100 to 140 gal, acceleration increased by 40%, and max base shear increased by 40%. From 140 to 200 gal, acceleration increased by 40%, and max base shear increased by 70%. From 200 to 280 gal, acceleration increased by 40%, and max base shear increased by 16.8%. From 280 to 400 gal, acceleration increased by 40%, and max base shear increased by 13.7%. As shown in Table 6, the numerical percentages varied each time the peak acceleration increased by 40% [39]; the corresponding base shear percentage increased, on the whole, until the 200–280 gal range at which point it appears to decrease. In effect, the wall was in an elastic state before this acceleration and then entered the elastic-plastic stage until it broke. These numerical values are consistent with the experimental phenomena we observed.Table 5
Maximum shear with peak acceleration (MLQT-1).
Specimen numberPeak acceleration (gal)Max base shear (kN)MLQT-13513.4970−31.149100−41.830140−58.136200−98.837280115.429400−131.231Table 6
Maximum shear with peak acceleration (MLQT-2).
Specimen numberPeak acceleration (gal)Max base shearMLQT-27011.67710016.00620034.85530055.54640072.75351098.816MLQT-2, at 70 to 100 gal, increased in applied acceleration by 40%, and max base shear increased by 37.04%. From 100 to 200 gal, acceleration increased by 50%, and max shear increased by 117.40%. From 200 to 300 gal, acceleration increased by 50%, and max shear increased by 59.36%. From 300 to 400 gal, acceleration increased by 30%, and max shear increased by 30.94%. From 400 to 510 gal, acceleration increased by 26%, and max shear increased by 35.82%. The numerical percentage variations are listed in Table6. We found that the LRB did not work before 200 gal and did work from 200 to 510 gal, and that its performance degraded from 400 to 510 gal.MLQT-3, from 100 to 200 gal, increased in applied acceleration by 50% [40] and max base shear increased by 40%. From 200 to 280 gal, acceleration increased by 40%, and max shear increased by 20%. From 280 to 400 gal, acceleration increased by 140%, and max shear increased by 36%. From 400 to 620 gal, acceleration increased by 50%, and max shear increased by 70%. From 620 to 800 gal, acceleration increased by 30%, and max shear increased by 10%. The numerical percentage variations are listed in Table 7. It appears that the LRB performed effective isolation from 100 to 400 gal and 200 to 400 gal, but this performance was unstable after 400 gal.Table 7
Maximum shear with peak acceleration (MLQT-3).
Specimen numberPeak acceleration (gal)Max base shearMLQT-31009.79120014.09728017.66040024.07362040.82580044.778MLQT-4, from 70 to 140 gal, increased in applied acceleration by 100% and in max base shear by 50%. From 140 to 200 gal, acceleration increased by 40%, and max base shear increased by 30%. From 200 to 280 gal, acceleration increased by 40%, and max base shear increased by 75%. From 280 to 400 gal, acceleration increased by 40%, and max base shear increased by 10%. From 620 to 800 gal, acceleration increased by 30%, and max base shear increased by 20%. From 800 to 1,000 gal, acceleration increased by 25%, and max base shear increased by 10%. From 1,000 to 1,024 gal, acceleration increased by 20%, and max base shear increased by 10%. From 1,240 to 1,400 gal, acceleration increased by 13%, and max base shear increased by 3%. The numerical percentage variations are listed in Table8. The LRB was either not working or had unstable performance before 280 gal; its performance was stable from 620 to 1,024 gal. The LRB showed optimal isolation performance once it started working and the specimen was on the brink of destruction.Table 8
Maximum shear with peak acceleration (MLQT-4).
Specimen numberPeak acceleration (gal)Max base shear (kN)MLQT-4702.0451403.0192004.0002807.0004007.90962011.20680013.825100015.207124016.556140017.021


### 3.3. Seismic Behavior

#### 3.3.1. Hysteretic Curve

We drew a hysteretic curve under the action of various earthquake intensities according to the test data presented above [41]. We used the MLQT-3 MLQT-4 corresponding to peak acceleration under simulated earthquake action as an example. The hysteretic curve characteristics of the bearing walls with different shape coefficients S1 and S2 [42] were compared to assess seismic absorption and isolation performance of two types of LRB wall. The hysteretic curve and hysteretic loop area of MLQT-4 under 800 gal are shown in Figure 9.Figure 9
MLQT-4 (800 gal) hysteretic curve: (a) hysteresis curve of 90 mm LRB wall; (b) maximum hysteresis loop area.
(a)(b)Upon reaching the yield load, MLQT-1 displacement growth accelerated though the hysteretic curve continued to present an inverse S shape. The load applied in the opposite direction decreased the wall’s stiffness due to the slip effect and yet-unclosed fractures [43]. Prior to fracture closure, the specimen had relatively low stiffness with cracked reinforcement, an uncracked block, and stressed concrete. Once the damaged wall came under pressure during the loading process, the deformation velocity decreased and the wall stiffness was enhanced.When the LRB started working in the MLQT-2 specimen, the area of the hysteresis ring was approximately full [44]. The wall’s carrying capacity remained stable after the LRB reached its maximum displacement. After the specimen reached the maximum load, the capacity of the wall decreased slowly, the displacement growth began to accelerate, the yield load of the hysteresis curve formed an inverse S shape, and the central part of the hysteretic curve grew fusiform.The MLQT-3 hysteretic curve presented an inverse S trend in the negative direction (Figure10) [45], which indicates that the shear deformation of the structure influenced the stiffness at the beginning of the load application process. This may be due to microcracks created inside the structure or slight relative displacement occurring at the junction. As the maximum acceleration reached 620 gal, the displacement continued to increase alongside a certain amount of plastic deformation. This formed a clear inverse S shape (Figure 10(d)) as the stiffness continued to degenerate.Figure 10
Specimens MLQT-3 hysteretic curve: (a) 100 gal; (b) 280 gal; (c) 400 gal; (d) 620 gal.
(a)(b)(c)(d)As per the MLQT-4 hysteretic curve, as shown in Figure11, the maximum lateral displacement of the model was within 6 mm when peak acceleration was 100 gal. Inputting a peak acceleration of 400 gal did not produce a sharp inverse S in the curve. The maximum displacement is closer to four times greater than the 100 gal maximum displacement in the hysteretic curve, which suggests that the stiffness of the specimen degraded or the LRB deformed horizontally in the test [46]. From 620 gal to 1,000 gal, the shape of the hysteretic curve did not markedly change as the work of the LRB was stable [47]. The hysteresis ring area increased when the peak acceleration reached 1,024 gal; the LRB was damaged at this time and the sample reached its limit of seismic resistance.Figure 11
MLQT-4 hysteretic curve: (a) 70 gal; (b) 100 gal; (c) 200 gal; (d) 280 gal; (e) 400 gal; (f) 620 gal; (g) 800 gal; (h) 1000 gal; (i) 1240 gal.
(a)(b)(c)(d)(e)(f)(g)(h)(i)


#### 3.3.2. Equivalent Horizontal Stiffness

Equivalent horizontal stiffness [48] is defined as a given component’s resistance to deformation under pressure. It obtains per unit displacement of the force in the elastic range. The equivalent horizontal stiffness in the isolation layer can be calculated as follows:(38)Kh=Q1−Q2X1−X2,where Kh is the isolation layer’s equivalent horizontal stiffness, Q1 is the maximal shear force, Q2 is the minimum shear force, X1 is the maximum displacement, and X2 is the minimum displacement.The relationship between isolation layers with each LRB equivalent horizontal stiffness is expressed as(39)Kh=∑Kj,where Kj is each LRB’s equivalent horizontal stiffness.We determined each LRB’s equivalent horizontal stiffness from MLQT-2, MLQT-3, and MLQT-4 for comparison against the design values. As shown in Figure12, until the LRB was broken and the test was stopped, the LRB equivalent horizontal stiffness in MLQT-2 and MLQT-3 was higher than the LRB 250% equivalent shear strain stiffness. The equivalent horizontal stiffness of MLQT-2-4 did not exceed the maximum seismic resistance design value at any point [49]. Theoretically, stiffness should decrease linearly along the curve as the seismic wave strength increases. Our curve, however, shows an increase under a certain internal seismic strength in MLQT-2 and -3; growth points emerge along the seismic intensity generally corresponding to the LRB seismic resistance design. Abnormal changes in stiffness curves did not affect the LRB in achieving the expected target despite the theoretical inconsistencies.Figure 12
Equivalent horizontal stiffness: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)We noted in our experiment that the lower connection plate of the LRB tended to separate from the rubber layer. This reduced the area of the connection between the rubber bottom and steel plate, causing a slight increase in vertical stress and slightly improving LRB horizontal stiffness. In MLQT-2 and MLQT-3, we found that the LRB components were broken due to the separation of the steel plate and rubber layer. The corresponding peak acceleration postyield stiffness and equivalent viscous damping ratio curve are further evidence of this phenomenon and further indicate that damage to the bearing affected the sample’s mechanical properties. The LRB must be repaired properly after an “intense” earthquake occurs.


#### 3.3.3. Postyield Stiffness

We obtained hysteretic curves with different levels of simulated seismic intensity according to the experimental data. We then analyzed the LRB equivalent horizontal stiffness, postyield stiffness, and equivalent viscous damping ratio mechanical properties [50] to compare them against those measured in the LRB factory to test the applicability of the proposed method.The postyield stiffness of the structural isolation layer is calculated as follows:(40)Kd=Qd1−Q1/X1+Qd2−Q2/X22,where Kd is the postyield stiffness, Qd1 is the hysteretic curve positive direction intersecting with the shear axis, and Qd2 is the hysteretic curve negative direction intersecting with the shear axis.We calculated the LRB postyield stiffness for MLQT-2-4 for comparison with the design values as shown in Figure13.Figure 13
Postyield stiffness: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)The MLQT-4 postyield stiffness curve appears to be consistent with the relevant theory, as the elastic deformation range of the LRB is large. After the LRB yields, the stiffness rapidly increases where the contact area between the rubber layer and the connecting plate is stripped away. Compared to other LRB components, this is more obvious in the lower LRB (which has less rubber).


#### 3.3.4. Equivalent Viscous Damping Ratio

The equivalent viscous damping ratio is a key parameter determining the equivalent period based on displacement in a seismic design [51]. The seismic isolation layer equivalent viscous damping ratio is calculated as follows:(41)ζ=2ΔWπKhX1−X22,where ζ is the structural equivalent viscous damping ratio, ΔW is the envelope area of the hysteretic curve, and Kh is the structural isolation layer equivalent horizontal stiffness.The equivalent viscous damping ratio of each LRB can be calculated as(42)ζ=∑KjζjKh,where ζ is the structural isolation layer equivalent viscous damping ratio, ζj represents each LRB’s equivalent viscous damping ratio, Kh is the structural isolation layer equivalent horizontal stiffness, and Kj represents each LRB’s equivalent horizontal stiffness in the experiment.Each LRB’s equivalent viscous damping ratio was obtained from MLQT-2-4 and compared with the 250% LRB shear deformation horizontal property as shown in Figure14.Figure 14
Equivalent viscous damping ratio: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)As shown in Figure14, the equivalent viscous damping ratios of MLQT-2-4 were not lower than the 250% horizontal deformation of each LRB at any point until the specimens were broken. The MLQT-2 curve fluctuated irregularly as earthquake intensity increased, which indicates that the elastic deformation range of the LRB is small. This may be related to its relatively low rubber content [52]. The MLQT-3 curve is more in line with the theoretical rationality of this design. The MLQT-4 curve begins to increase at 620 cm/s2, which is equivalent to 10 degrees of earthquake protection. The LRB was damaged under higher seismic intensity; it appears that a higher body coefficient renders higher seismic resistance.


## 3.1. Relevant Phenomena

Our specimen MLQT-1 is an uninstalled LRB wall. The corresponding relation between earthquake intensity and peak acceleration of seismic waves was imposed under 70, 100, 300, and 400 gal El Centro wave excitation models. We found no obvious cracks until reaching displacement of 9.6 mm and horizontal force of +137.879 kN, at which point the corresponding peak acceleration of seismic waves was 300 gal. We continued reverse loading until finding oblique cracks in the wall as shown in Figure5; these shear cracks developed in the wall during pushing and pulling in a typical “X” shape [35], at which point the wall presented shear failure.Figure 5
MLQT-1 collapse state.MLQT-2 was equipped with a 50 mm LRB masonry specimen and input with a damping ratio of 0.15 and 70, 100, 300, 400, and 510 gal El Centro wave excitation models [36]. We calculated the LRB’s ultimate displacement according to the upper and lower beam relative displacement. When the acceleration peak reached 300 gal, the maximum displacement was 13.6 mm (exceeding LRB allowable limit of displacement 10.8 mm). The middle of the LRB was broken when the acceleration peak reached 510 gal; the wall was not damaged until the test was stopped, as shown in Figure 6.Figure 6
MLQT-2: (a) LRB limit of displacement state; (b) LRB broken state.
(a)(b)As shown in Figure6, MLQT-2 does not meet the design standard though the LRB and upper structure were not broken in the test. LRB thickness variations, in this case, were −1 mm, −2 mm, and −1 mm, which indicate that the wall was not compressed by the vertical load.Specimen MLQT-3 was equipped with a 57 mm LRB masonry specimen designed to be resistant to seismic waves of 400 gal peak acceleration. Once the peak acceleration reached 400 gal, the maximum displacement was 14.2 mm, i.e., below the allowable limit of 22.5 mm. The contact between the LRB and steel plate separated only slightly, which did not affect the continuation of the experiment. The maximum displacement was observed upon reaching 620 gal, at 26.4 mm, which does exceed the allowable limit. The specimen was considered to be broken at this time. The LRB internal steel plate was broken and fully exposed upon reaching 800 gal (Figure7), though the wall did not crack by the end of the experiment.Figure 7
MLQT-3: (a) LRB limit of displacement state; (b) LRB broken state.
(a)(b)MLQT-4 was equipped with a 90 mm LRB masonry specimen and designed to resist the peak acceleration of an 800 gal seismic wave [37]. The maximum displacement reached 39.1 mm at the peak acceleration of 800 gal, which is below the allowable limit of 60.5 mm; thus, this specimen met the predesign specifications. Upon reaching 1,240 gal, the maximum displacement was 65.7 mm (exceeding the allowable limit of 60.5 mm) and the specimen was considered to be broken (Figure 8). Upon reaching 1,400 gal, the maximum displacement of the model structure was 76.9 mm (push), the maximum base shear was 17 kN (push), the LRB was broken, and the wall showed no cracks.Figure 8
MLQT-4: (a) LRB limit of displacement state; (b) LRB broken state.
(a)(b)As shown in Figure8, the MLQT-4 LRB thickness variations after the vertical loads were applied were −1 mm, −1 mm, and −2 mm. The maximum values of the additional loads were, respectively, 187.3 kN and 184.5 kN. This difference is in the allowable range and the wall was not compressed.


## 3.2. Data Operating of Maximum Displacement and Shear

We found that MLQT-1 280–400 gal showed a maximum displacement difference of −3 mm by −4 minus −1 mm, −7 minus −4 mm; the value changed suddenly and dramatically by the end of the test and diagonal cracks appeared through the wall. MLQT-2 reached a maximum displacement exceeding the allowable limit at 300 gal, as listed in Table3.Table 3
MLQT-1\MLQT-2 displacement.
Specimen numberPeak acceleration (gal)Maximum displacement (mm)MLQT-1140−1.879200−4.280280−7.2514008.779MLQT-22009.99030013.61040016.65751021.497The MLQT-3 allowable displacement is 22.5 mm [38], which was exceeded at 620 gal. The MLQT-4 allowable displacement is 60.05 mm, which was exceeded at 1,240 gal (Table 4).Table 4
MLQT-3\MLQT-4 displacement.
Specimen numberPeak acceleration (gal)Maximum displacementMLQT-340014.231 mm (allowable disp 22.50 mm)26.421 mm620800MLQT-4100050.943 mm (allowable disp 60.05 mm)65.780 mm12401400MLQT-1, from 35 to 70 gal, applied acceleration increase of 100%, and then max base shear increased by 130%. For an increased acceleration of 40%, from 70 to 100 gal, the max base shear increased by 34% (Table5). From 100 to 140 gal, acceleration increased by 40%, and max base shear increased by 40%. From 140 to 200 gal, acceleration increased by 40%, and max base shear increased by 70%. From 200 to 280 gal, acceleration increased by 40%, and max base shear increased by 16.8%. From 280 to 400 gal, acceleration increased by 40%, and max base shear increased by 13.7%. As shown in Table 6, the numerical percentages varied each time the peak acceleration increased by 40% [39]; the corresponding base shear percentage increased, on the whole, until the 200–280 gal range at which point it appears to decrease. In effect, the wall was in an elastic state before this acceleration and then entered the elastic-plastic stage until it broke. These numerical values are consistent with the experimental phenomena we observed.Table 5
Maximum shear with peak acceleration (MLQT-1).
Specimen numberPeak acceleration (gal)Max base shear (kN)MLQT-13513.4970−31.149100−41.830140−58.136200−98.837280115.429400−131.231Table 6
Maximum shear with peak acceleration (MLQT-2).
Specimen numberPeak acceleration (gal)Max base shearMLQT-27011.67710016.00620034.85530055.54640072.75351098.816MLQT-2, at 70 to 100 gal, increased in applied acceleration by 40%, and max base shear increased by 37.04%. From 100 to 200 gal, acceleration increased by 50%, and max shear increased by 117.40%. From 200 to 300 gal, acceleration increased by 50%, and max shear increased by 59.36%. From 300 to 400 gal, acceleration increased by 30%, and max shear increased by 30.94%. From 400 to 510 gal, acceleration increased by 26%, and max shear increased by 35.82%. The numerical percentage variations are listed in Table6. We found that the LRB did not work before 200 gal and did work from 200 to 510 gal, and that its performance degraded from 400 to 510 gal.MLQT-3, from 100 to 200 gal, increased in applied acceleration by 50% [40] and max base shear increased by 40%. From 200 to 280 gal, acceleration increased by 40%, and max shear increased by 20%. From 280 to 400 gal, acceleration increased by 140%, and max shear increased by 36%. From 400 to 620 gal, acceleration increased by 50%, and max shear increased by 70%. From 620 to 800 gal, acceleration increased by 30%, and max shear increased by 10%. The numerical percentage variations are listed in Table 7. It appears that the LRB performed effective isolation from 100 to 400 gal and 200 to 400 gal, but this performance was unstable after 400 gal.Table 7
Maximum shear with peak acceleration (MLQT-3).
Specimen numberPeak acceleration (gal)Max base shearMLQT-31009.79120014.09728017.66040024.07362040.82580044.778MLQT-4, from 70 to 140 gal, increased in applied acceleration by 100% and in max base shear by 50%. From 140 to 200 gal, acceleration increased by 40%, and max base shear increased by 30%. From 200 to 280 gal, acceleration increased by 40%, and max base shear increased by 75%. From 280 to 400 gal, acceleration increased by 40%, and max base shear increased by 10%. From 620 to 800 gal, acceleration increased by 30%, and max base shear increased by 20%. From 800 to 1,000 gal, acceleration increased by 25%, and max base shear increased by 10%. From 1,000 to 1,024 gal, acceleration increased by 20%, and max base shear increased by 10%. From 1,240 to 1,400 gal, acceleration increased by 13%, and max base shear increased by 3%. The numerical percentage variations are listed in Table8. The LRB was either not working or had unstable performance before 280 gal; its performance was stable from 620 to 1,024 gal. The LRB showed optimal isolation performance once it started working and the specimen was on the brink of destruction.Table 8
Maximum shear with peak acceleration (MLQT-4).
Specimen numberPeak acceleration (gal)Max base shear (kN)MLQT-4702.0451403.0192004.0002807.0004007.90962011.20680013.825100015.207124016.556140017.021


## 3.3. Seismic Behavior

### 3.3.1. Hysteretic Curve

We drew a hysteretic curve under the action of various earthquake intensities according to the test data presented above [41]. We used the MLQT-3 MLQT-4 corresponding to peak acceleration under simulated earthquake action as an example. The hysteretic curve characteristics of the bearing walls with different shape coefficients S1 and S2 [42] were compared to assess seismic absorption and isolation performance of two types of LRB wall. The hysteretic curve and hysteretic loop area of MLQT-4 under 800 gal are shown in Figure 9.Figure 9
MLQT-4 (800 gal) hysteretic curve: (a) hysteresis curve of 90 mm LRB wall; (b) maximum hysteresis loop area.
(a)(b)Upon reaching the yield load, MLQT-1 displacement growth accelerated though the hysteretic curve continued to present an inverse S shape. The load applied in the opposite direction decreased the wall’s stiffness due to the slip effect and yet-unclosed fractures [43]. Prior to fracture closure, the specimen had relatively low stiffness with cracked reinforcement, an uncracked block, and stressed concrete. Once the damaged wall came under pressure during the loading process, the deformation velocity decreased and the wall stiffness was enhanced.When the LRB started working in the MLQT-2 specimen, the area of the hysteresis ring was approximately full [44]. The wall’s carrying capacity remained stable after the LRB reached its maximum displacement. After the specimen reached the maximum load, the capacity of the wall decreased slowly, the displacement growth began to accelerate, the yield load of the hysteresis curve formed an inverse S shape, and the central part of the hysteretic curve grew fusiform.The MLQT-3 hysteretic curve presented an inverse S trend in the negative direction (Figure10) [45], which indicates that the shear deformation of the structure influenced the stiffness at the beginning of the load application process. This may be due to microcracks created inside the structure or slight relative displacement occurring at the junction. As the maximum acceleration reached 620 gal, the displacement continued to increase alongside a certain amount of plastic deformation. This formed a clear inverse S shape (Figure 10(d)) as the stiffness continued to degenerate.Figure 10
Specimens MLQT-3 hysteretic curve: (a) 100 gal; (b) 280 gal; (c) 400 gal; (d) 620 gal.
(a)(b)(c)(d)As per the MLQT-4 hysteretic curve, as shown in Figure11, the maximum lateral displacement of the model was within 6 mm when peak acceleration was 100 gal. Inputting a peak acceleration of 400 gal did not produce a sharp inverse S in the curve. The maximum displacement is closer to four times greater than the 100 gal maximum displacement in the hysteretic curve, which suggests that the stiffness of the specimen degraded or the LRB deformed horizontally in the test [46]. From 620 gal to 1,000 gal, the shape of the hysteretic curve did not markedly change as the work of the LRB was stable [47]. The hysteresis ring area increased when the peak acceleration reached 1,024 gal; the LRB was damaged at this time and the sample reached its limit of seismic resistance.Figure 11
MLQT-4 hysteretic curve: (a) 70 gal; (b) 100 gal; (c) 200 gal; (d) 280 gal; (e) 400 gal; (f) 620 gal; (g) 800 gal; (h) 1000 gal; (i) 1240 gal.
(a)(b)(c)(d)(e)(f)(g)(h)(i)


### 3.3.2. Equivalent Horizontal Stiffness

Equivalent horizontal stiffness [48] is defined as a given component’s resistance to deformation under pressure. It obtains per unit displacement of the force in the elastic range. The equivalent horizontal stiffness in the isolation layer can be calculated as follows:(38)Kh=Q1−Q2X1−X2,where Kh is the isolation layer’s equivalent horizontal stiffness, Q1 is the maximal shear force, Q2 is the minimum shear force, X1 is the maximum displacement, and X2 is the minimum displacement.The relationship between isolation layers with each LRB equivalent horizontal stiffness is expressed as(39)Kh=∑Kj,where Kj is each LRB’s equivalent horizontal stiffness.We determined each LRB’s equivalent horizontal stiffness from MLQT-2, MLQT-3, and MLQT-4 for comparison against the design values. As shown in Figure12, until the LRB was broken and the test was stopped, the LRB equivalent horizontal stiffness in MLQT-2 and MLQT-3 was higher than the LRB 250% equivalent shear strain stiffness. The equivalent horizontal stiffness of MLQT-2-4 did not exceed the maximum seismic resistance design value at any point [49]. Theoretically, stiffness should decrease linearly along the curve as the seismic wave strength increases. Our curve, however, shows an increase under a certain internal seismic strength in MLQT-2 and -3; growth points emerge along the seismic intensity generally corresponding to the LRB seismic resistance design. Abnormal changes in stiffness curves did not affect the LRB in achieving the expected target despite the theoretical inconsistencies.Figure 12
Equivalent horizontal stiffness: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)We noted in our experiment that the lower connection plate of the LRB tended to separate from the rubber layer. This reduced the area of the connection between the rubber bottom and steel plate, causing a slight increase in vertical stress and slightly improving LRB horizontal stiffness. In MLQT-2 and MLQT-3, we found that the LRB components were broken due to the separation of the steel plate and rubber layer. The corresponding peak acceleration postyield stiffness and equivalent viscous damping ratio curve are further evidence of this phenomenon and further indicate that damage to the bearing affected the sample’s mechanical properties. The LRB must be repaired properly after an “intense” earthquake occurs.


### 3.3.3. Postyield Stiffness

We obtained hysteretic curves with different levels of simulated seismic intensity according to the experimental data. We then analyzed the LRB equivalent horizontal stiffness, postyield stiffness, and equivalent viscous damping ratio mechanical properties [50] to compare them against those measured in the LRB factory to test the applicability of the proposed method.The postyield stiffness of the structural isolation layer is calculated as follows:(40)Kd=Qd1−Q1/X1+Qd2−Q2/X22,where Kd is the postyield stiffness, Qd1 is the hysteretic curve positive direction intersecting with the shear axis, and Qd2 is the hysteretic curve negative direction intersecting with the shear axis.We calculated the LRB postyield stiffness for MLQT-2-4 for comparison with the design values as shown in Figure13.Figure 13
Postyield stiffness: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)The MLQT-4 postyield stiffness curve appears to be consistent with the relevant theory, as the elastic deformation range of the LRB is large. After the LRB yields, the stiffness rapidly increases where the contact area between the rubber layer and the connecting plate is stripped away. Compared to other LRB components, this is more obvious in the lower LRB (which has less rubber).


### 3.3.4. Equivalent Viscous Damping Ratio

The equivalent viscous damping ratio is a key parameter determining the equivalent period based on displacement in a seismic design [51]. The seismic isolation layer equivalent viscous damping ratio is calculated as follows:(41)ζ=2ΔWπKhX1−X22,where ζ is the structural equivalent viscous damping ratio, ΔW is the envelope area of the hysteretic curve, and Kh is the structural isolation layer equivalent horizontal stiffness.The equivalent viscous damping ratio of each LRB can be calculated as(42)ζ=∑KjζjKh,where ζ is the structural isolation layer equivalent viscous damping ratio, ζj represents each LRB’s equivalent viscous damping ratio, Kh is the structural isolation layer equivalent horizontal stiffness, and Kj represents each LRB’s equivalent horizontal stiffness in the experiment.Each LRB’s equivalent viscous damping ratio was obtained from MLQT-2-4 and compared with the 250% LRB shear deformation horizontal property as shown in Figure14.Figure 14
Equivalent viscous damping ratio: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)As shown in Figure14, the equivalent viscous damping ratios of MLQT-2-4 were not lower than the 250% horizontal deformation of each LRB at any point until the specimens were broken. The MLQT-2 curve fluctuated irregularly as earthquake intensity increased, which indicates that the elastic deformation range of the LRB is small. This may be related to its relatively low rubber content [52]. The MLQT-3 curve is more in line with the theoretical rationality of this design. The MLQT-4 curve begins to increase at 620 cm/s2, which is equivalent to 10 degrees of earthquake protection. The LRB was damaged under higher seismic intensity; it appears that a higher body coefficient renders higher seismic resistance.


## 3.3.1. Hysteretic Curve

We drew a hysteretic curve under the action of various earthquake intensities according to the test data presented above [41]. We used the MLQT-3 MLQT-4 corresponding to peak acceleration under simulated earthquake action as an example. The hysteretic curve characteristics of the bearing walls with different shape coefficients S1 and S2 [42] were compared to assess seismic absorption and isolation performance of two types of LRB wall. The hysteretic curve and hysteretic loop area of MLQT-4 under 800 gal are shown in Figure 9.Figure 9
MLQT-4 (800 gal) hysteretic curve: (a) hysteresis curve of 90 mm LRB wall; (b) maximum hysteresis loop area.
(a)(b)Upon reaching the yield load, MLQT-1 displacement growth accelerated though the hysteretic curve continued to present an inverse S shape. The load applied in the opposite direction decreased the wall’s stiffness due to the slip effect and yet-unclosed fractures [43]. Prior to fracture closure, the specimen had relatively low stiffness with cracked reinforcement, an uncracked block, and stressed concrete. Once the damaged wall came under pressure during the loading process, the deformation velocity decreased and the wall stiffness was enhanced.When the LRB started working in the MLQT-2 specimen, the area of the hysteresis ring was approximately full [44]. The wall’s carrying capacity remained stable after the LRB reached its maximum displacement. After the specimen reached the maximum load, the capacity of the wall decreased slowly, the displacement growth began to accelerate, the yield load of the hysteresis curve formed an inverse S shape, and the central part of the hysteretic curve grew fusiform.The MLQT-3 hysteretic curve presented an inverse S trend in the negative direction (Figure10) [45], which indicates that the shear deformation of the structure influenced the stiffness at the beginning of the load application process. This may be due to microcracks created inside the structure or slight relative displacement occurring at the junction. As the maximum acceleration reached 620 gal, the displacement continued to increase alongside a certain amount of plastic deformation. This formed a clear inverse S shape (Figure 10(d)) as the stiffness continued to degenerate.Figure 10
Specimens MLQT-3 hysteretic curve: (a) 100 gal; (b) 280 gal; (c) 400 gal; (d) 620 gal.
(a)(b)(c)(d)As per the MLQT-4 hysteretic curve, as shown in Figure11, the maximum lateral displacement of the model was within 6 mm when peak acceleration was 100 gal. Inputting a peak acceleration of 400 gal did not produce a sharp inverse S in the curve. The maximum displacement is closer to four times greater than the 100 gal maximum displacement in the hysteretic curve, which suggests that the stiffness of the specimen degraded or the LRB deformed horizontally in the test [46]. From 620 gal to 1,000 gal, the shape of the hysteretic curve did not markedly change as the work of the LRB was stable [47]. The hysteresis ring area increased when the peak acceleration reached 1,024 gal; the LRB was damaged at this time and the sample reached its limit of seismic resistance.Figure 11
MLQT-4 hysteretic curve: (a) 70 gal; (b) 100 gal; (c) 200 gal; (d) 280 gal; (e) 400 gal; (f) 620 gal; (g) 800 gal; (h) 1000 gal; (i) 1240 gal.
(a)(b)(c)(d)(e)(f)(g)(h)(i)


## 3.3.2. Equivalent Horizontal Stiffness

Equivalent horizontal stiffness [48] is defined as a given component’s resistance to deformation under pressure. It obtains per unit displacement of the force in the elastic range. The equivalent horizontal stiffness in the isolation layer can be calculated as follows:(38)Kh=Q1−Q2X1−X2,where Kh is the isolation layer’s equivalent horizontal stiffness, Q1 is the maximal shear force, Q2 is the minimum shear force, X1 is the maximum displacement, and X2 is the minimum displacement.The relationship between isolation layers with each LRB equivalent horizontal stiffness is expressed as(39)Kh=∑Kj,where Kj is each LRB’s equivalent horizontal stiffness.We determined each LRB’s equivalent horizontal stiffness from MLQT-2, MLQT-3, and MLQT-4 for comparison against the design values. As shown in Figure12, until the LRB was broken and the test was stopped, the LRB equivalent horizontal stiffness in MLQT-2 and MLQT-3 was higher than the LRB 250% equivalent shear strain stiffness. The equivalent horizontal stiffness of MLQT-2-4 did not exceed the maximum seismic resistance design value at any point [49]. Theoretically, stiffness should decrease linearly along the curve as the seismic wave strength increases. Our curve, however, shows an increase under a certain internal seismic strength in MLQT-2 and -3; growth points emerge along the seismic intensity generally corresponding to the LRB seismic resistance design. Abnormal changes in stiffness curves did not affect the LRB in achieving the expected target despite the theoretical inconsistencies.Figure 12
Equivalent horizontal stiffness: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)We noted in our experiment that the lower connection plate of the LRB tended to separate from the rubber layer. This reduced the area of the connection between the rubber bottom and steel plate, causing a slight increase in vertical stress and slightly improving LRB horizontal stiffness. In MLQT-2 and MLQT-3, we found that the LRB components were broken due to the separation of the steel plate and rubber layer. The corresponding peak acceleration postyield stiffness and equivalent viscous damping ratio curve are further evidence of this phenomenon and further indicate that damage to the bearing affected the sample’s mechanical properties. The LRB must be repaired properly after an “intense” earthquake occurs.


## 3.3.3. Postyield Stiffness

We obtained hysteretic curves with different levels of simulated seismic intensity according to the experimental data. We then analyzed the LRB equivalent horizontal stiffness, postyield stiffness, and equivalent viscous damping ratio mechanical properties [50] to compare them against those measured in the LRB factory to test the applicability of the proposed method.The postyield stiffness of the structural isolation layer is calculated as follows:(40)Kd=Qd1−Q1/X1+Qd2−Q2/X22,where Kd is the postyield stiffness, Qd1 is the hysteretic curve positive direction intersecting with the shear axis, and Qd2 is the hysteretic curve negative direction intersecting with the shear axis.We calculated the LRB postyield stiffness for MLQT-2-4 for comparison with the design values as shown in Figure13.Figure 13
Postyield stiffness: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)The MLQT-4 postyield stiffness curve appears to be consistent with the relevant theory, as the elastic deformation range of the LRB is large. After the LRB yields, the stiffness rapidly increases where the contact area between the rubber layer and the connecting plate is stripped away. Compared to other LRB components, this is more obvious in the lower LRB (which has less rubber).


## 3.3.4. Equivalent Viscous Damping Ratio

The equivalent viscous damping ratio is a key parameter determining the equivalent period based on displacement in a seismic design [51]. The seismic isolation layer equivalent viscous damping ratio is calculated as follows:(41)ζ=2ΔWπKhX1−X22,where ζ is the structural equivalent viscous damping ratio, ΔW is the envelope area of the hysteretic curve, and Kh is the structural isolation layer equivalent horizontal stiffness.The equivalent viscous damping ratio of each LRB can be calculated as(42)ζ=∑KjζjKh,where ζ is the structural isolation layer equivalent viscous damping ratio, ζj represents each LRB’s equivalent viscous damping ratio, Kh is the structural isolation layer equivalent horizontal stiffness, and Kj represents each LRB’s equivalent horizontal stiffness in the experiment.Each LRB’s equivalent viscous damping ratio was obtained from MLQT-2-4 and compared with the 250% LRB shear deformation horizontal property as shown in Figure14.Figure 14
Equivalent viscous damping ratio: (a) MLQT-2; (b) MLQT-3; (c) MLQT-4.
(a)(b)(c)As shown in Figure14, the equivalent viscous damping ratios of MLQT-2-4 were not lower than the 250% horizontal deformation of each LRB at any point until the specimens were broken. The MLQT-2 curve fluctuated irregularly as earthquake intensity increased, which indicates that the elastic deformation range of the LRB is small. This may be related to its relatively low rubber content [52]. The MLQT-3 curve is more in line with the theoretical rationality of this design. The MLQT-4 curve begins to increase at 620 cm/s2, which is equivalent to 10 degrees of earthquake protection. The LRB was damaged under higher seismic intensity; it appears that a higher body coefficient renders higher seismic resistance.


## 4. LRB Experiment and Fitting Formula

### 4.1. Experiment Method

We conducted a horizontal stiffness experiment on five types of small LRB to explore various horizontal shearing strains and second shape coefficients as they affect the equivalent horizontal stiffness, yield stiffness, and hysteretic properties of the structure. We gave each shearing strain fitting formula an equivalent horizontal stiffness and yield stiffness while varying the second shape coefficient [2]. We analyzed the LRB horizontal stiffness related to the structure model seismic coefficient and horizontal displacement to create a curve-fitting formula. We used a small-diameter LRB to reflect the seismic effects on existing rural low masonry structures, as our goal was to gather results that may serve as a workable reference for bearing wall designs.Residential buildings generally contain circle rubber bearings, the vertical stiffness and bearing capacity of which are related to the single rubber effective pressure area and free sides of the superficial area (first shape coefficient). Their stability is related to the inner rubber layer diameter divided by the total thickness of the inner rubber (second shape coefficient).The LRB has the same steady vertical bearing capacity as other bearing components, in addition to favorable isolation, self-resetting, damping energy dissipation properties, and a large damping ratio (±25%). The LRB resilience model, which is commonly utilized in urban sky-rises, is a double linear model that has a much larger natural rubber bearing and high damping rubber bearing damping ratio; it is ineffective for rural dwelling isolation components. We used technical methods in this study to decrease the LRB damping ratio and observe the isolation effects in small- and medium-intensity earthquakes. We focused on the horizontal property design indicators of postyield stiffness, yield force, equivalent stiffness, and equivalent damping ratio.We adopted inverse calculation methods for our LRB designs and built the seismic response spectrum method into our parameter calculations [10]. The LRB lead core rod has about 7 MP yield stress with certain plastic deformation and repetitive loading properties. The original properties are restored upon the end of the deformation period. The LRB must be designed with sufficient horizontal stiffness to suit the opposite deformation between the foundation and new construction; the LRB can then absorb sufficient earthquake energy to remit seismic resistance to the structure.(43)We≫E,We=12Qmaxδmax,E=12Mve2.The maximum shear coefficient can be calculated as follows:(44)α=QmaxMg,Qmax=KHδmax,and the mass, stiffness, and period relation is formulated here as(45)We=Mg2T2α28π2,where M is the total mass, g is the gravitational acceleration, KH is the horizontal stiffness coefficient, and T is the natural period.


### 4.2. Horizontal Stiffness

#### 4.2.1. Horizontal Stiffness Calculations

We obtained an LRB horizontal displacement-load-related curve from our experimental data adopting a sine wave of 0.5 Hz and shear strain of 100% third hysteresis. We calculated the experimental LRB equivalent horizontal stiffnessKh, equivalent damping ratio heq, and postyield stiffness Kd [6] as follows:(46)Kh=Q1−Q2X1−X2,(47)heq=2ΔWπKhX1−X22,(48)Kd=12Q1−Qd1X1+Q2−Qd2X2,where Q1 is the maximum shear force, Q2 is the minimum shear force, X1 is the maximum displacement, X2 is the minimum displacement, Qd1 and Qd2 are the hysteresis curve positive and negative direction with the shear force coordinate as the crossing point, and ΔW is the hysteresis curve envelope size. The calculation results are listed in Table 2.


#### 4.2.2. Effects on Shear Strain

We calculated the equivalent horizontal stiffness related to shear strain from the experimental LRB hysteresis curve as shown in Figure15. The equivalent horizontal stiffness decreased as the shear strain increased in a nonlinear relation—a higher shear strain gave the LRB a better ability to withstand horizontal stiffness degeneration. The LRB equivalent horizontal stiffness affected shear strain when it was relatively low (e.g., γ = 50%) and was stable when the shear strain was larger (e.g., γ = 100%; γ = 250%).Figure 15
Relative equivalent horizontal stiffness with shear strain.Figure16 shows our experimental LRB postyield stiffness related to horizontal shear strain. In comparison to the effects on equivalent horizontal stiffness, the shear strain does not act strongly upon the postyield degenerate ability of the structure. The experimental LRB postyield decreased as the shear strain increased. When the shear strain was γ = 100–γ = 250%, the experimental LRB postyield change was confined to 15%.Figure 16
Relative postyield stiffness with shear strain.


#### 4.2.3. Effects on Second Shape Coefficient

The second shape coefficient (S2) is determined by dividing the respective experimental LRB layers’ rubber diameters [7]. The fifth LRB in our setup has a different rubber shear elasticity modulus than the other four LRBs, so we only compared the first four LRBs’ equivalent horizontal stiffness related to the S2 curve (Figure 17). We found that the equivalent horizontal stiffness increased, as did the shear strain and increasing rate, as S2 increased.Figure 17
Relative equivalent horizontal stiffness with the second shape coefficient.Formulas (49)–(51) show where shear strains of 50%, 100%, and 250% change alongside the equivalent horizontal stiffness Kh with S2. The results of these formulas provide a workable reference for small-diameter LRB designs.(49)Kh50%=39.1S22−267.4S2+678.9S2≤4.8,(50)Kh100%=34.9S22−236.0S2+560.3S2≤4.8,(51)Kh250%=52.2S22−374.7S2+777.5S2≤4.8.As shown in Figure17, when S2 is 4.8, the equivalent horizontal stiffness-second shape coefficient curve slope mean value is 35.89; its maximum value does not exceed 63. When S2 exceeds 4.8, the equivalent horizontal stiffness slope increases more steeply and the slope mean value is 226.22. The upper formula listed above is more accurate than the others in terms of the stable equivalent horizontal stiffness when S2 is 4.8.We also calculated the experimental LRB postyield stiffnessKd with the S2 relative curve, as shown in Figure 18.Figure 18
Postyield stiffness related to the second shape coefficient.Under the same vertical pressure action, the postyield stiffnessKd increased as S2 increased. At a certain S2, the Kd decreased and trended toward a constant value. The two relations can be described linearly using formulas (52)–(54). A higher S2 can raise the LRB postyield stiffness within a certain range, but this effect is not obvious when S2 exceeds 4.8. When selecting the LRB for use in rural masonry structures, restricting the second shape coefficient to 4.8 can ensure a safe structure while minimizing cost and reducing material waste.(52)Kd50%=21.1S2+29.7S2≤4.8,(53)Kd100%=69.2S2+112.6S2≤4.8,(54)Kd250%=48.3S2+4.5S2≤4.8.


### 4.3. Analysis Progress

We next drew a pressure displacement-pressure curve to explore the vertical contraction stiffness after three loading cycles as shown in Table9 [8]. We found that the supporter vertical stiffness discreteness increased over this curve, which indicates that the component’s vertical pressure property is not stable.Table 9
Pressure stiffness Kv calculation results.
Bearing1#2#3#4#5#Number vertical1-11-22-12-23-13-24-14-25-15-2StiffnessKv (kN/mm)47.9157.2242.0442.0429.8629.0138.8736.7962.4271.03Figure4 shows the shear strain hysteretic loop at a frequency of 0.5 Hz.As shown in Figure4, the 1#–4# hysteretic loops are rhomboid and have straight lines along their sides. Component 5# was designed with a different rubber than the other four; its hysteretic loop is slightly bent in a positive direction. The 1# bearing hysteretic curve shear strain was 100% as shown in the third loop. The 0.1 Hz hysteretic loop areas for 1–5# are 115.98, 178.71, 222.59, 137.52, and 139.77, respectively. The 0.5 Hz hysteretic loop areas for 1–5# are 88.01, 175.2, 216.46, 134.46, and 122.7, respectively. The 0.1 Hz hysteretic loop area is larger than the others [53]. We drew its curve again for comparison as shown in Figure 19 to explore the relatively stable changes in the LRB load-displacement curve at higher loading frequency.Figure 19
1# bearing third hysteretic loop comparison: (a) 0.5 Hz loading frequency; (b) 0.1 Hz loading frequency.
(a)(b)We separately input a 0.5 Hz sine wave to the five experimental LRB components and imposed maximum displacement loading of 50%, 100%, and 250%. There was no apparent difference among the LRBs in this case, so we redrew the LRB GZD-0.5 hysteretic curve for comparison as shown in Figure20.Figure 20
GZ5-0.5 hysteretic curves with varying shear strain.When the shear strain was 50% or 100%, the LRB hysteretic curves were rhomboid and every loop could be divided into four stages according to the loading process. We observed linear regularity among the horizontal loading and displacement indicators in this case. When the shear strain was 250%, the hysteretic curve formed an inverse S in a positive direction. The LRB stiffness degenerated to a certain extent as the inner rubber of the bearing slid along the steel plate. The area of the hysteretic curve increased as the shear strain increased, as the supporter degraded the energy resistance.


### 4.4. Effect of Supporter Stiffness on Isolation

#### 4.4.1. Model Design

The LRB isolation effects can be balanced according to the horizontal seismic coefficient and horizontal displacement. Seismic fortification intensity is generally 6–8 degrees (0.20 g); in rural areas, the horizontal seismic coefficient should be less than 0.4. The isolation structure’s horizontal displacement should be less than 0.55 times the bearing’s effective diameter or 3-fold greater than the inner rubber’s total thickness [54]. For the purposes of our experiment, we set the LRB horizontal displacement to an effective diameter of 0.55 multiplied by 60.5 mm.We used a typical rural masonry dwelling to design our finite element calculation model and shaking table test model. The model is a two-layer single-bay brick masonry house with a reduced-scale (1 : 2) size of 2.25 m × 1.95 m × 3.3 m and a wall thickness of 120 mm. The four corners of the “house” were given four identical LRBs. The door was a one-layer 450 mm wide opening with sliding supporters (LRBs) on either side; the upper part of the slider was connected to the structure while the lower part slid freely along the pedestal plane.We used finite element software SAP2000 to construct the isolated structure model. LRBs 1–4# were, respectively, imposed as model supporters and El Centro waves were input separately inX and Y directions. The peak acceleration was 8.5 degrees and adjusted to 300 gal. We used the model to calculate the structure’s horizontal seismic coefficient and maximum horizontal displacement.


#### 4.4.2. Rural Building Isolated Parameter Estimation

We analyzed the above model’s horizontal seismic coefficient and maximum horizontal displacement to determine the experimental LRB postyield stiffness and equivalent horizontal stiffness, as shown in Figures21–24.Figure 21
Seismic coefficient related to postyield stiffness.Figure 22
Horizontal displacement related to postyield stiffness.Figure 23
Seismic coefficient related to equivalent horizontal stiffness.Figure 24
Horizontal displacement related to equivalent horizontal stiffness.The model seismic coefficient showed an increasing and nearly linear trend with LRB postyield addition. The coefficients were similar inX and Y directions when the supporter postyield stiffness was relatively small; the difference between them increased as this stiffness increased.Formulas (46) and (47) are the X- and Y-direction seismic coefficient (βx, βy)-postyield stiffness (Kd) relations, respectively.(55)βx=1.6×10−3Kd+0.122,βy=8.6×10−4Kd+0.178.The model horizontal displacement decreased as the supporter postyield stiffness increased in a nearly linear relation [55]. The differences in X- and Y-direction horizontal displacement were within 3 mm, which suggests that postyield stiffness does not affect the disparity between positive and negative displacement values. We established formula (56) as a reference for selecting the largest and safest possible displacement value according to the horizontal displacement (u)-postyield stiffness (Kd) curve:(56)u=−0.147Kd+84.709.As shown in Figures21–24, the model seismic coefficient and horizontal displacement change regularly with equivalent horizontal stiffness. Both indicators showed similar postyield relations, but neither are linear. We fitted the curve to write a formula for the XY direction seismic coefficient βxβy-equivalent horizontal stiffness Kh relation (formulas (57) and (58)). The horizontal displacement μ-equivalent horizontal stiffness Kh relative curve, which contains a quadratic function, is arguably more accurate; the X-direction is fitted in this case by formula (59).(57)βx=8.1×10−4Kh+0.147,(58)βy=1.51×10−3Kh+0.065,(59)μ=9.1×10−4Kh2−0.583kh+139.87.LRB postyield stiffness should be 149.7–167.8 N/mm and equivalent horizontal stiffness should be 193.9–218.65 N/mm when the horizontal displacement is less than 60.5 mm and the seismic coefficient is below 0.4 requirements, so we selected the 4# LRB accordingly to conduct our shaking table test.


#### 4.4.3. Seismic Coefficient

The isolated masonry structure horizontal seismic coefficient can be determined by multistorey masonry calculation. A smaller seismic coefficient indicates stronger seismic resistance but also generally reflects greater horizontal displacement in the LRB. The horizontal displacement thus must be restricted within a certain value; in this case, the isolation structure’s horizontal seismic coefficient should be less than 0.4.We constructed a test bearing based on actual rural dwelling design and material requirements and then subjected it to a shaking table test. After inputting an 8.5 degree (300 gal peak acceleration) El Centro wave, we observed the deformation concentration in the isolation layer; the upper structure was altered but did not present any apparent cracks. When the seismic wave peak acceleration reached 800 gal, the structure remained intact and its maximum horizontal displacement value was permissible. To this effect, the proposed LRB is safe and feasible.In an earthquake, the isolation layer presents horizontal displacement which creates a displacement in the upper structure. The seismic structure’s top storey may show greater displacement under certain earthquake conditions even when the isolation layer’s displacement is relatively small. The support will be destroyed when the isolation layer’s horizontal displacement exceeds the LRB limit; in this case, the entire upper structure may fail. We assert that support displacement should not exceed a minimum effective diameter of 0.55- or 3-fold greater than the inner rubber’s total thickness. The effective LRB diameters in this experiment were 110 mm, in which case the horizontal displacement limit is 60.5 mm.We simulated an earthquake wave acting on this model in SAP2000. Again, five LRBs were modeled with 8.5 fortification intensity and a 0.3 g structure base displacement average value. As shown in Table10, only the 4# LRB satisfied both base displacement less than 60.5 mm and seismic coefficient less than 0.4.Table 10
LRB base displacement and seismic coefficient results.
Supporter numberBase displacement (mm)Seismic coefficient1#50.4120.5072#72.2610.2653#64.60.3824#56.110.3735#44.2270.544


## 4.1. Experiment Method

We conducted a horizontal stiffness experiment on five types of small LRB to explore various horizontal shearing strains and second shape coefficients as they affect the equivalent horizontal stiffness, yield stiffness, and hysteretic properties of the structure. We gave each shearing strain fitting formula an equivalent horizontal stiffness and yield stiffness while varying the second shape coefficient [2]. We analyzed the LRB horizontal stiffness related to the structure model seismic coefficient and horizontal displacement to create a curve-fitting formula. We used a small-diameter LRB to reflect the seismic effects on existing rural low masonry structures, as our goal was to gather results that may serve as a workable reference for bearing wall designs.Residential buildings generally contain circle rubber bearings, the vertical stiffness and bearing capacity of which are related to the single rubber effective pressure area and free sides of the superficial area (first shape coefficient). Their stability is related to the inner rubber layer diameter divided by the total thickness of the inner rubber (second shape coefficient).The LRB has the same steady vertical bearing capacity as other bearing components, in addition to favorable isolation, self-resetting, damping energy dissipation properties, and a large damping ratio (±25%). The LRB resilience model, which is commonly utilized in urban sky-rises, is a double linear model that has a much larger natural rubber bearing and high damping rubber bearing damping ratio; it is ineffective for rural dwelling isolation components. We used technical methods in this study to decrease the LRB damping ratio and observe the isolation effects in small- and medium-intensity earthquakes. We focused on the horizontal property design indicators of postyield stiffness, yield force, equivalent stiffness, and equivalent damping ratio.We adopted inverse calculation methods for our LRB designs and built the seismic response spectrum method into our parameter calculations [10]. The LRB lead core rod has about 7 MP yield stress with certain plastic deformation and repetitive loading properties. The original properties are restored upon the end of the deformation period. The LRB must be designed with sufficient horizontal stiffness to suit the opposite deformation between the foundation and new construction; the LRB can then absorb sufficient earthquake energy to remit seismic resistance to the structure.(43)We≫E,We=12Qmaxδmax,E=12Mve2.The maximum shear coefficient can be calculated as follows:(44)α=QmaxMg,Qmax=KHδmax,and the mass, stiffness, and period relation is formulated here as(45)We=Mg2T2α28π2,where M is the total mass, g is the gravitational acceleration, KH is the horizontal stiffness coefficient, and T is the natural period.


## 4.2. Horizontal Stiffness

### 4.2.1. Horizontal Stiffness Calculations

We obtained an LRB horizontal displacement-load-related curve from our experimental data adopting a sine wave of 0.5 Hz and shear strain of 100% third hysteresis. We calculated the experimental LRB equivalent horizontal stiffnessKh, equivalent damping ratio heq, and postyield stiffness Kd [6] as follows:(46)Kh=Q1−Q2X1−X2,(47)heq=2ΔWπKhX1−X22,(48)Kd=12Q1−Qd1X1+Q2−Qd2X2,where Q1 is the maximum shear force, Q2 is the minimum shear force, X1 is the maximum displacement, X2 is the minimum displacement, Qd1 and Qd2 are the hysteresis curve positive and negative direction with the shear force coordinate as the crossing point, and ΔW is the hysteresis curve envelope size. The calculation results are listed in Table 2.


### 4.2.2. Effects on Shear Strain

We calculated the equivalent horizontal stiffness related to shear strain from the experimental LRB hysteresis curve as shown in Figure15. The equivalent horizontal stiffness decreased as the shear strain increased in a nonlinear relation—a higher shear strain gave the LRB a better ability to withstand horizontal stiffness degeneration. The LRB equivalent horizontal stiffness affected shear strain when it was relatively low (e.g., γ = 50%) and was stable when the shear strain was larger (e.g., γ = 100%; γ = 250%).Figure 15
Relative equivalent horizontal stiffness with shear strain.Figure16 shows our experimental LRB postyield stiffness related to horizontal shear strain. In comparison to the effects on equivalent horizontal stiffness, the shear strain does not act strongly upon the postyield degenerate ability of the structure. The experimental LRB postyield decreased as the shear strain increased. When the shear strain was γ = 100–γ = 250%, the experimental LRB postyield change was confined to 15%.Figure 16
Relative postyield stiffness with shear strain.


### 4.2.3. Effects on Second Shape Coefficient

The second shape coefficient (S2) is determined by dividing the respective experimental LRB layers’ rubber diameters [7]. The fifth LRB in our setup has a different rubber shear elasticity modulus than the other four LRBs, so we only compared the first four LRBs’ equivalent horizontal stiffness related to the S2 curve (Figure 17). We found that the equivalent horizontal stiffness increased, as did the shear strain and increasing rate, as S2 increased.Figure 17
Relative equivalent horizontal stiffness with the second shape coefficient.Formulas (49)–(51) show where shear strains of 50%, 100%, and 250% change alongside the equivalent horizontal stiffness Kh with S2. The results of these formulas provide a workable reference for small-diameter LRB designs.(49)Kh50%=39.1S22−267.4S2+678.9S2≤4.8,(50)Kh100%=34.9S22−236.0S2+560.3S2≤4.8,(51)Kh250%=52.2S22−374.7S2+777.5S2≤4.8.As shown in Figure17, when S2 is 4.8, the equivalent horizontal stiffness-second shape coefficient curve slope mean value is 35.89; its maximum value does not exceed 63. When S2 exceeds 4.8, the equivalent horizontal stiffness slope increases more steeply and the slope mean value is 226.22. The upper formula listed above is more accurate than the others in terms of the stable equivalent horizontal stiffness when S2 is 4.8.We also calculated the experimental LRB postyield stiffnessKd with the S2 relative curve, as shown in Figure 18.Figure 18
Postyield stiffness related to the second shape coefficient.Under the same vertical pressure action, the postyield stiffnessKd increased as S2 increased. At a certain S2, the Kd decreased and trended toward a constant value. The two relations can be described linearly using formulas (52)–(54). A higher S2 can raise the LRB postyield stiffness within a certain range, but this effect is not obvious when S2 exceeds 4.8. When selecting the LRB for use in rural masonry structures, restricting the second shape coefficient to 4.8 can ensure a safe structure while minimizing cost and reducing material waste.(52)Kd50%=21.1S2+29.7S2≤4.8,(53)Kd100%=69.2S2+112.6S2≤4.8,(54)Kd250%=48.3S2+4.5S2≤4.8.


## 4.2.1. Horizontal Stiffness Calculations

We obtained an LRB horizontal displacement-load-related curve from our experimental data adopting a sine wave of 0.5 Hz and shear strain of 100% third hysteresis. We calculated the experimental LRB equivalent horizontal stiffnessKh, equivalent damping ratio heq, and postyield stiffness Kd [6] as follows:(46)Kh=Q1−Q2X1−X2,(47)heq=2ΔWπKhX1−X22,(48)Kd=12Q1−Qd1X1+Q2−Qd2X2,where Q1 is the maximum shear force, Q2 is the minimum shear force, X1 is the maximum displacement, X2 is the minimum displacement, Qd1 and Qd2 are the hysteresis curve positive and negative direction with the shear force coordinate as the crossing point, and ΔW is the hysteresis curve envelope size. The calculation results are listed in Table 2.


## 4.2.2. Effects on Shear Strain

We calculated the equivalent horizontal stiffness related to shear strain from the experimental LRB hysteresis curve as shown in Figure15. The equivalent horizontal stiffness decreased as the shear strain increased in a nonlinear relation—a higher shear strain gave the LRB a better ability to withstand horizontal stiffness degeneration. The LRB equivalent horizontal stiffness affected shear strain when it was relatively low (e.g., γ = 50%) and was stable when the shear strain was larger (e.g., γ = 100%; γ = 250%).Figure 15
Relative equivalent horizontal stiffness with shear strain.Figure16 shows our experimental LRB postyield stiffness related to horizontal shear strain. In comparison to the effects on equivalent horizontal stiffness, the shear strain does not act strongly upon the postyield degenerate ability of the structure. The experimental LRB postyield decreased as the shear strain increased. When the shear strain was γ = 100–γ = 250%, the experimental LRB postyield change was confined to 15%.Figure 16
Relative postyield stiffness with shear strain.


## 4.2.3. Effects on Second Shape Coefficient

The second shape coefficient (S2) is determined by dividing the respective experimental LRB layers’ rubber diameters [7]. The fifth LRB in our setup has a different rubber shear elasticity modulus than the other four LRBs, so we only compared the first four LRBs’ equivalent horizontal stiffness related to the S2 curve (Figure 17). We found that the equivalent horizontal stiffness increased, as did the shear strain and increasing rate, as S2 increased.Figure 17
Relative equivalent horizontal stiffness with the second shape coefficient.Formulas (49)–(51) show where shear strains of 50%, 100%, and 250% change alongside the equivalent horizontal stiffness Kh with S2. The results of these formulas provide a workable reference for small-diameter LRB designs.(49)Kh50%=39.1S22−267.4S2+678.9S2≤4.8,(50)Kh100%=34.9S22−236.0S2+560.3S2≤4.8,(51)Kh250%=52.2S22−374.7S2+777.5S2≤4.8.As shown in Figure17, when S2 is 4.8, the equivalent horizontal stiffness-second shape coefficient curve slope mean value is 35.89; its maximum value does not exceed 63. When S2 exceeds 4.8, the equivalent horizontal stiffness slope increases more steeply and the slope mean value is 226.22. The upper formula listed above is more accurate than the others in terms of the stable equivalent horizontal stiffness when S2 is 4.8.We also calculated the experimental LRB postyield stiffnessKd with the S2 relative curve, as shown in Figure 18.Figure 18
Postyield stiffness related to the second shape coefficient.Under the same vertical pressure action, the postyield stiffnessKd increased as S2 increased. At a certain S2, the Kd decreased and trended toward a constant value. The two relations can be described linearly using formulas (52)–(54). A higher S2 can raise the LRB postyield stiffness within a certain range, but this effect is not obvious when S2 exceeds 4.8. When selecting the LRB for use in rural masonry structures, restricting the second shape coefficient to 4.8 can ensure a safe structure while minimizing cost and reducing material waste.(52)Kd50%=21.1S2+29.7S2≤4.8,(53)Kd100%=69.2S2+112.6S2≤4.8,(54)Kd250%=48.3S2+4.5S2≤4.8.


## 4.3. Analysis Progress

We next drew a pressure displacement-pressure curve to explore the vertical contraction stiffness after three loading cycles as shown in Table9 [8]. We found that the supporter vertical stiffness discreteness increased over this curve, which indicates that the component’s vertical pressure property is not stable.Table 9
Pressure stiffness Kv calculation results.
Bearing1#2#3#4#5#Number vertical1-11-22-12-23-13-24-14-25-15-2StiffnessKv (kN/mm)47.9157.2242.0442.0429.8629.0138.8736.7962.4271.03Figure4 shows the shear strain hysteretic loop at a frequency of 0.5 Hz.As shown in Figure4, the 1#–4# hysteretic loops are rhomboid and have straight lines along their sides. Component 5# was designed with a different rubber than the other four; its hysteretic loop is slightly bent in a positive direction. The 1# bearing hysteretic curve shear strain was 100% as shown in the third loop. The 0.1 Hz hysteretic loop areas for 1–5# are 115.98, 178.71, 222.59, 137.52, and 139.77, respectively. The 0.5 Hz hysteretic loop areas for 1–5# are 88.01, 175.2, 216.46, 134.46, and 122.7, respectively. The 0.1 Hz hysteretic loop area is larger than the others [53]. We drew its curve again for comparison as shown in Figure 19 to explore the relatively stable changes in the LRB load-displacement curve at higher loading frequency.Figure 19
1# bearing third hysteretic loop comparison: (a) 0.5 Hz loading frequency; (b) 0.1 Hz loading frequency.
(a)(b)We separately input a 0.5 Hz sine wave to the five experimental LRB components and imposed maximum displacement loading of 50%, 100%, and 250%. There was no apparent difference among the LRBs in this case, so we redrew the LRB GZD-0.5 hysteretic curve for comparison as shown in Figure20.Figure 20
GZ5-0.5 hysteretic curves with varying shear strain.When the shear strain was 50% or 100%, the LRB hysteretic curves were rhomboid and every loop could be divided into four stages according to the loading process. We observed linear regularity among the horizontal loading and displacement indicators in this case. When the shear strain was 250%, the hysteretic curve formed an inverse S in a positive direction. The LRB stiffness degenerated to a certain extent as the inner rubber of the bearing slid along the steel plate. The area of the hysteretic curve increased as the shear strain increased, as the supporter degraded the energy resistance.


## 4.4. Effect of Supporter Stiffness on Isolation

### 4.4.1. Model Design

The LRB isolation effects can be balanced according to the horizontal seismic coefficient and horizontal displacement. Seismic fortification intensity is generally 6–8 degrees (0.20 g); in rural areas, the horizontal seismic coefficient should be less than 0.4. The isolation structure’s horizontal displacement should be less than 0.55 times the bearing’s effective diameter or 3-fold greater than the inner rubber’s total thickness [54]. For the purposes of our experiment, we set the LRB horizontal displacement to an effective diameter of 0.55 multiplied by 60.5 mm.We used a typical rural masonry dwelling to design our finite element calculation model and shaking table test model. The model is a two-layer single-bay brick masonry house with a reduced-scale (1 : 2) size of 2.25 m × 1.95 m × 3.3 m and a wall thickness of 120 mm. The four corners of the “house” were given four identical LRBs. The door was a one-layer 450 mm wide opening with sliding supporters (LRBs) on either side; the upper part of the slider was connected to the structure while the lower part slid freely along the pedestal plane.We used finite element software SAP2000 to construct the isolated structure model. LRBs 1–4# were, respectively, imposed as model supporters and El Centro waves were input separately inX and Y directions. The peak acceleration was 8.5 degrees and adjusted to 300 gal. We used the model to calculate the structure’s horizontal seismic coefficient and maximum horizontal displacement.


### 4.4.2. Rural Building Isolated Parameter Estimation

We analyzed the above model’s horizontal seismic coefficient and maximum horizontal displacement to determine the experimental LRB postyield stiffness and equivalent horizontal stiffness, as shown in Figures21–24.Figure 21
Seismic coefficient related to postyield stiffness.Figure 22
Horizontal displacement related to postyield stiffness.Figure 23
Seismic coefficient related to equivalent horizontal stiffness.Figure 24
Horizontal displacement related to equivalent horizontal stiffness.The model seismic coefficient showed an increasing and nearly linear trend with LRB postyield addition. The coefficients were similar inX and Y directions when the supporter postyield stiffness was relatively small; the difference between them increased as this stiffness increased.Formulas (46) and (47) are the X- and Y-direction seismic coefficient (βx, βy)-postyield stiffness (Kd) relations, respectively.(55)βx=1.6×10−3Kd+0.122,βy=8.6×10−4Kd+0.178.The model horizontal displacement decreased as the supporter postyield stiffness increased in a nearly linear relation [55]. The differences in X- and Y-direction horizontal displacement were within 3 mm, which suggests that postyield stiffness does not affect the disparity between positive and negative displacement values. We established formula (56) as a reference for selecting the largest and safest possible displacement value according to the horizontal displacement (u)-postyield stiffness (Kd) curve:(56)u=−0.147Kd+84.709.As shown in Figures21–24, the model seismic coefficient and horizontal displacement change regularly with equivalent horizontal stiffness. Both indicators showed similar postyield relations, but neither are linear. We fitted the curve to write a formula for the XY direction seismic coefficient βxβy-equivalent horizontal stiffness Kh relation (formulas (57) and (58)). The horizontal displacement μ-equivalent horizontal stiffness Kh relative curve, which contains a quadratic function, is arguably more accurate; the X-direction is fitted in this case by formula (59).(57)βx=8.1×10−4Kh+0.147,(58)βy=1.51×10−3Kh+0.065,(59)μ=9.1×10−4Kh2−0.583kh+139.87.LRB postyield stiffness should be 149.7–167.8 N/mm and equivalent horizontal stiffness should be 193.9–218.65 N/mm when the horizontal displacement is less than 60.5 mm and the seismic coefficient is below 0.4 requirements, so we selected the 4# LRB accordingly to conduct our shaking table test.


### 4.4.3. Seismic Coefficient

The isolated masonry structure horizontal seismic coefficient can be determined by multistorey masonry calculation. A smaller seismic coefficient indicates stronger seismic resistance but also generally reflects greater horizontal displacement in the LRB. The horizontal displacement thus must be restricted within a certain value; in this case, the isolation structure’s horizontal seismic coefficient should be less than 0.4.We constructed a test bearing based on actual rural dwelling design and material requirements and then subjected it to a shaking table test. After inputting an 8.5 degree (300 gal peak acceleration) El Centro wave, we observed the deformation concentration in the isolation layer; the upper structure was altered but did not present any apparent cracks. When the seismic wave peak acceleration reached 800 gal, the structure remained intact and its maximum horizontal displacement value was permissible. To this effect, the proposed LRB is safe and feasible.In an earthquake, the isolation layer presents horizontal displacement which creates a displacement in the upper structure. The seismic structure’s top storey may show greater displacement under certain earthquake conditions even when the isolation layer’s displacement is relatively small. The support will be destroyed when the isolation layer’s horizontal displacement exceeds the LRB limit; in this case, the entire upper structure may fail. We assert that support displacement should not exceed a minimum effective diameter of 0.55- or 3-fold greater than the inner rubber’s total thickness. The effective LRB diameters in this experiment were 110 mm, in which case the horizontal displacement limit is 60.5 mm.We simulated an earthquake wave acting on this model in SAP2000. Again, five LRBs were modeled with 8.5 fortification intensity and a 0.3 g structure base displacement average value. As shown in Table10, only the 4# LRB satisfied both base displacement less than 60.5 mm and seismic coefficient less than 0.4.Table 10
LRB base displacement and seismic coefficient results.
Supporter numberBase displacement (mm)Seismic coefficient1#50.4120.5072#72.2610.2653#64.60.3824#56.110.3735#44.2270.544


## 4.4.1. Model Design

The LRB isolation effects can be balanced according to the horizontal seismic coefficient and horizontal displacement. Seismic fortification intensity is generally 6–8 degrees (0.20 g); in rural areas, the horizontal seismic coefficient should be less than 0.4. The isolation structure’s horizontal displacement should be less than 0.55 times the bearing’s effective diameter or 3-fold greater than the inner rubber’s total thickness [54]. For the purposes of our experiment, we set the LRB horizontal displacement to an effective diameter of 0.55 multiplied by 60.5 mm.We used a typical rural masonry dwelling to design our finite element calculation model and shaking table test model. The model is a two-layer single-bay brick masonry house with a reduced-scale (1 : 2) size of 2.25 m × 1.95 m × 3.3 m and a wall thickness of 120 mm. The four corners of the “house” were given four identical LRBs. The door was a one-layer 450 mm wide opening with sliding supporters (LRBs) on either side; the upper part of the slider was connected to the structure while the lower part slid freely along the pedestal plane.We used finite element software SAP2000 to construct the isolated structure model. LRBs 1–4# were, respectively, imposed as model supporters and El Centro waves were input separately inX and Y directions. The peak acceleration was 8.5 degrees and adjusted to 300 gal. We used the model to calculate the structure’s horizontal seismic coefficient and maximum horizontal displacement.


## 4.4.2. Rural Building Isolated Parameter Estimation

We analyzed the above model’s horizontal seismic coefficient and maximum horizontal displacement to determine the experimental LRB postyield stiffness and equivalent horizontal stiffness, as shown in Figures21–24.Figure 21
Seismic coefficient related to postyield stiffness.Figure 22
Horizontal displacement related to postyield stiffness.Figure 23
Seismic coefficient related to equivalent horizontal stiffness.Figure 24
Horizontal displacement related to equivalent horizontal stiffness.The model seismic coefficient showed an increasing and nearly linear trend with LRB postyield addition. The coefficients were similar inX and Y directions when the supporter postyield stiffness was relatively small; the difference between them increased as this stiffness increased.Formulas (46) and (47) are the X- and Y-direction seismic coefficient (βx, βy)-postyield stiffness (Kd) relations, respectively.(55)βx=1.6×10−3Kd+0.122,βy=8.6×10−4Kd+0.178.The model horizontal displacement decreased as the supporter postyield stiffness increased in a nearly linear relation [55]. The differences in X- and Y-direction horizontal displacement were within 3 mm, which suggests that postyield stiffness does not affect the disparity between positive and negative displacement values. We established formula (56) as a reference for selecting the largest and safest possible displacement value according to the horizontal displacement (u)-postyield stiffness (Kd) curve:(56)u=−0.147Kd+84.709.As shown in Figures21–24, the model seismic coefficient and horizontal displacement change regularly with equivalent horizontal stiffness. Both indicators showed similar postyield relations, but neither are linear. We fitted the curve to write a formula for the XY direction seismic coefficient βxβy-equivalent horizontal stiffness Kh relation (formulas (57) and (58)). The horizontal displacement μ-equivalent horizontal stiffness Kh relative curve, which contains a quadratic function, is arguably more accurate; the X-direction is fitted in this case by formula (59).(57)βx=8.1×10−4Kh+0.147,(58)βy=1.51×10−3Kh+0.065,(59)μ=9.1×10−4Kh2−0.583kh+139.87.LRB postyield stiffness should be 149.7–167.8 N/mm and equivalent horizontal stiffness should be 193.9–218.65 N/mm when the horizontal displacement is less than 60.5 mm and the seismic coefficient is below 0.4 requirements, so we selected the 4# LRB accordingly to conduct our shaking table test.


## 4.4.3. Seismic Coefficient

The isolated masonry structure horizontal seismic coefficient can be determined by multistorey masonry calculation. A smaller seismic coefficient indicates stronger seismic resistance but also generally reflects greater horizontal displacement in the LRB. The horizontal displacement thus must be restricted within a certain value; in this case, the isolation structure’s horizontal seismic coefficient should be less than 0.4.We constructed a test bearing based on actual rural dwelling design and material requirements and then subjected it to a shaking table test. After inputting an 8.5 degree (300 gal peak acceleration) El Centro wave, we observed the deformation concentration in the isolation layer; the upper structure was altered but did not present any apparent cracks. When the seismic wave peak acceleration reached 800 gal, the structure remained intact and its maximum horizontal displacement value was permissible. To this effect, the proposed LRB is safe and feasible.In an earthquake, the isolation layer presents horizontal displacement which creates a displacement in the upper structure. The seismic structure’s top storey may show greater displacement under certain earthquake conditions even when the isolation layer’s displacement is relatively small. The support will be destroyed when the isolation layer’s horizontal displacement exceeds the LRB limit; in this case, the entire upper structure may fail. We assert that support displacement should not exceed a minimum effective diameter of 0.55- or 3-fold greater than the inner rubber’s total thickness. The effective LRB diameters in this experiment were 110 mm, in which case the horizontal displacement limit is 60.5 mm.We simulated an earthquake wave acting on this model in SAP2000. Again, five LRBs were modeled with 8.5 fortification intensity and a 0.3 g structure base displacement average value. As shown in Table10, only the 4# LRB satisfied both base displacement less than 60.5 mm and seismic coefficient less than 0.4.Table 10
LRB base displacement and seismic coefficient results.
Supporter numberBase displacement (mm)Seismic coefficient1#50.4120.5072#72.2610.2653#64.60.3824#56.110.3735#44.2270.544


## 5. Conclusions

We conducted an experiment to test the seismic performance and mechanical failure of several LRB structural specimens in this study, each with different parameters and excited by seismic waves of varying intensity. Our conclusions can be summarized as follows.(1)
The 50 mm thick LRB specimen MLQT-2 was ineffective as its rubber layer is too thin (3 × 1.2 mm). The 57 mm thick LRB specimen MLQT-3 LRB exceeds the limit displacement for strong-intensity seismic waves. LRB specimen MLQT-4, which is 90 mm thick, can withstand middle-intensity seismic waves. The LRB exceeds the limit displacement for strong-intensity seismic waves in this case, but the LRB and structure were unbroken by the end of the test (though any stronger seismic waves would cause the LRB to break)(2)
MLQT-4 was shown to withstand stronger peak acceleration seismic waves than other specimens due to its effective masonry performance and mortar strength. Changes in LRB thickness changes do not appear to affect bearing capacity but do change the position of the centroid. We recommend that the superstructure be designed according to extant standards and that the sliding bearing isolator support be placed in the appropriate position according to the stresses acting on the structure(3)
We wrote a seismic coefficient fitting formula that reflects the regular changes in small-diameter LRB horizontal displacement according to our shaking table test results. This formula may serve as a workable reference for practical LRB designs(4)
The structure’s horizontal seismic coefficient and horizontal displacement have a fittable relation with the LRB’s horizontal stiffness (provided that said stiffness falls within a reasonable range in the shaking table test)(5)
Various shear strains produce different forms of LRB equivalent horizontal stiffness and postyield indicators according to the second shape coefficient. When this coefficient is less than 4.8, the LRB equivalent horizontal stiffness changes in a stable manner and the curve of the two indicators can be accurately fitted. Restricting the second shape coefficient below 4.8 can assist in reducing the cost and material waste of the structure’s fabrication(6)
We analyzed the LRB horizontal stiffness with the structure model seismic coefficient and horizontal displacement relation to establish a curve-fitting formula. We input our shaking table test results to find that a small-diameter LRB remits ideal seismic resistance in typical, rural low masonry structures



---

*Source: 1345920-2020-12-29.xml*