Research Article  Open Access
Han Liu, Ping Tan, Fulin Zhou, "Stiffness Analysis of Rectangular Isolators Reinforced by Engineering Plastics", Shock and Vibration, vol. 2019, Article ID 3074834, 13 pages, 2019. https://doi.org/10.1155/2019/3074834
Stiffness Analysis of Rectangular Isolators Reinforced by Engineering Plastics
Abstract
A novel costeffective isolator reinforced by engineering plastics has been designed and manufactured for seismic protection for lowrise buildings in less developed areas. The reinforcement is flexible in tension, which is similar to fiberreinforced isolators. However, available solutions for fiberreinforced isolators are not applicable, because the Poisson effect of engineering plastics cannot be neglected, which is done for fiber reinforcement. In this paper, analytical solutions for compression and bending stiffness for rectangular isolators reinforced by engineering plastics are proposed, with both the Poisson effect of the reinforcement and the effect of rubber compressibility taken into consideration. Then, the simplified solutions are also derived, which can greatly improve calculation efficiency. To validate the solutions, finite element analysis is conducted on a set of isolators with different reinforcement stiffnesses. The results show the superiority of the proposed solutions to the previous solutions for fiberreinforced isolators. A series of experimental tests of the isolators are also carried out to verify the solutions. Both the analytical and the simplified solutions match well with the experimental results.
1. Introduction
Recent earthquakes in China such as Wenchuan and Yushu earthquakes caused numerous loss of lives, which has emphasized that a greater lose may happen when earthquakes occur in rural areas and less developed areas. Ge et al. [1] investigated the buildings in these areas to find that most of them are properly designed masonry structures with insufficient seismic capacity. About 80% land of China is located in seismically active zones with an intensity magnitude of six or above, and 58% building area locates in rural areas, according to the statistics by Cao et al. [2]. These buildings need to be protected to mitigate potential risk due to severe earthquakes.
Zhou [3] systematically introduced the method for seismic protection, and one of the most effective techniques is seismic isolation technology. However, conventional isolators are expensive, large, and heavy. Kelly [4] developed fiberreinforced isolators employing fibers to replace steel plates in laminated rubber isolators, which could reduce both the weight and the cost. In these isolators, the connecting plates are removed, which leads to an unbonded boundary condition. Kelly [5] and Kelly and Takhirov [6] made experimental attempts to prove that such kind of isolators could be produced and could provide enough capacity and proper stiffness. Unfortunately, the fibers have no flexure rigidity to restrain the rollover deformation of the isolator. ToopchiNezhad et al. [7] observed the rollover deformation of rectangular fiberreinforced isolators. When it happened, the horizontal stiffness decreased apparently. Kelly and Calabrese [8] derived that the displacement of this kind of isolators was limited to 1.67 times the height of the isolator. Angeli et al. [9] tested 17 pairs of isolators and proposed the deformed pattern of the rolloff deformation. The horizontal deformation is limited as long as the free surfaces contact with the loading plane. Pauletta et al. [10] observed an unstable rollout deformation of which the forcedisplacement curves had negative slopes. As a promotion, Tan et al. [11] proposed an innovative isolator which employed engineering plastics as reinforcement. Xu et al. [12] conducted systematical experiments of the isolators, employing two different kinds of plastic material. It was proved that engineering plastics have much higher outofplane stiffness than fibers, leading to higher stiffness and better ultimate properties of the isolator. As the mechanic assumptions of engineering plastics are different from fiber layers, available stiffness formulas of fiberreinforced isolators are not appropriate to calculate this kind of isolators.
One of the earliest theoretical analyses of bonded rubber was conducted by Rocard [13] using energy method. Solutions of strip isolators and circular isolators were derived by Gent and Lindley [14] with the incompressibility assumption of the rubber. This solution was extended by Gent and Meinecke [15] to compression and bending stiffness of rectangle isolators. Koh and Lim [16] considered rubber compressibility for rectangular isolators and found that the incompressible assumption overestimated the stiffness of the isolators considerably. Tsai and Pai [17] proposed simplified solutions for circular, square, and stripshaped isolators based on analytical solutions which considered rubber compressibility. All the above solutions were based on the assumption that the reinforcements were rigid.
Fiber reinforcement is flexible in tension, so the problem becomes much more complex. A widely acceptable method to simplify the problem is the pressure solution method introduced by Kelly [18]. The normal stresses are assumed to be equal, and some of the shear stress components are assumed to be zero. Based on the method, several analytical solutions for flexible reinforcement isolators have been developed. Kelly and Takhirov [6] derived stiffness of striptype isolators considering reinforcement extensibility, but the solution was exact only for plane problems. Tsai and Kelly [19] studied the compression and bending modulus of rectangular isolators based on the assumption of rubber incompressibility. Kelly and Calabrese [20] solved the stiffness solutions, considering both reinforcement extension and rubber compressibility for circular isolators. Angeli et al. [9] further derived the solutions for rectangular isolators, but the Poisson effect of the reinforcement is neglected. They considered the fibers laid out in two orthogonal directions and simply overlapped, and the fibers did not sustain shear stress. Pinarbasi and Mengi [21] solved the stress and modulus solutions of fiberreinforced isolators without the assumptions of pressure solution method. But the equations were complex, and only the solutions for stripshaped isolators were achieved.
In engineering practice, even though the value of Poisson’s ratio of rubber may be very close to 0.5, the effect of compressibility is considerable to isolator stiffness and should not be neglected. For engineering plastics, the shear stress and Poisson effect cannot be neglected either. In this paper, the solutions for compression and bending stiffness of rectangular isolators reinforced by engineering plastics considering both rubber compressibility and reinforcement Poisson effect is presented. The derived solutions are expressed in terms of Fourier series which are complex and to the disadvantage of popularizing isolation technology. Thus, the weighted residual approach is introduced to solve the partial differential equation of the stress, and a semianalytical solution in a much simpler form is achieved. To validate the solutions, experimental studies are conducted, and finite element analysis is also carried out. The proposed solutions are proved to have enough accuracy.
2. FiberReinforced Plastic Plate Isolator
Fiberreinforced plastic plate is usually employed as the reinforcement in the novel rectangular isolator, which is a common kind of engineering plastic and belongs to composite material. It was made up of unsaturated polyester matrix and glass fiber cloth in previous research [11]. The isolator was fabricated into a prism, as shown in Figure 1. Its dimension is decided by the demand of its mechanical properties as well as its installation convenience. The cross section of the isolator is illustrated in Figure 2.
The isolator is much cheaper and lighter than conventional steelreinforced isolators of the same size, and the processing of plastic is much easier than steel. The isolator promotes the displacement ability of fiberreinforced isolators. Unlike fiberreinforced isolators, the reinforcement of this isolator is composite material made up of plastic matrix and embedded glass fiber. So, the outofplane stiffness of the reinforcement is much higher than fiber sheets. As a result, the rollover phenomenon is limited and the vertical surface of the isolator does not contact with the loading plane even at very large horizontal deformation. It is to be noted that the stresses in two orthogonal directions in this isolator cannot be decoupled as they are in fiberreinforced isolators because of the plastic matrix. So, the analytical deviations should consider the Poisson effect of the reinforcement.
3. Compression Analysis of the Isolator
3.1. Equilibrium in Rubber Layer
A fiberreinforced plastic plate isolator is constituted of multiple layers. Assuming the stress components of each layer along the height of the isolator is uniform, a layer of rubber and two layers of plastic plates are extracted as an analytical unit. Considering the analytical unit subjected to a vertical compression load P, the coordinate system and its deformation under P are plotted in Figure 3. The thickness of the rubber layer is t, and the thickness of the plastic plate layer is t_{f}. The side length is 2b along the xaxis and 2a along the yaxis. The total compression of the rubber is . The aspect ratio is defined by .
The displacements of the rubber are assumed to obey the kinematic assumption as followed by Kelly and Calabrese [20]:where and are due to the lateral bulging of the rubber in the x and y directions, respectively, and and are due to the extension of the plastic plate.
The pressure solution method [18] is used here as follows:where p is the internal pressure. The shear stress components τ_{xz} and τ_{yz} are of order , and τ_{xy} is of order , where l is a characteristic length in the xy plane and equals 2b in this problem. Thus,
According to the above equations and the constitutive equations of linear elastic material, there are relationships in rubber derived as follows [8]:
The compressibility of rubber layer is considered bywhere K is the bulk modulus of the rubber. Substituting equations (1)–(3) into equation (10) and integrating through the thickness of the layer lead towhere .
3.2. Equilibrium in Plastic Plate Layer
In a single layer of plastic plate, several layers of glass fiber cloth are embedded, of which the directions are centrosymmetric. The base material of unsaturated polyester is isotropic. Thus, the plastic plate is assumed to be an isotropic material with a modulus of E_{f}. An analytical unit of the size dx × dy × t_{f} is extracted, and its stress state is shown in Figure 4, where F_{x} and F_{y} represent the normal forces of thickness t_{f} per unit length in x and y directions, and F_{xy} represents the inplane shear force. These internal forces are balanced by the shear stresses of rubber, τ_{xz} and τ_{yz}, which are transmitted from adjacent layers.
The equilibrium equations in x and y directions, respectively, are
The constitutive equations considering Poisson ratio of the plates are as follows:where ν is the Poisson ratio of the plate. When ν equals zero, equations (14)–(16) become the same as the ones derived by Angeli et al. [9] for fiberreinforced elastomeric bearings.
3.3. Equilibrium Equation of Pressure
The inplane shear strain component of the rubber is as follows:
Substituting equations (1) and (2) into equation (17), and then integrating it through the thickness, and then multiplying the result by the shear modulus , the following equation is obtained:
Equation (18) means is of the same order as or higher. From the pressure method assumption, as shown in equation (5), can be assumed to be zero. With this assumption, equation (16) is simplified and we have the following equation:
Substituting equations (6)–(9) and (19) into equations (12) and (13), the equilibrium equations of the plate are expressed as follows:
It is to be noted that the rubber and plastic plate layers are free in the xy plane. Thus, there are boundary conditions of ; ; ; and .
Integrating equations (20) and (21) with respect to x and y, respectively, and then substituting boundary conditions into them, the following equations are obtained:
Combining equations (14), (15), (22), and (23) leads to
Differentiating equations (8) and (9) with respect to x and y, respectively, and adding them up, and substituting the sum and equation (24) into equation (11) leads towhere is the Laplace operator (i.e., ) and
The Poisson ratio of the reinforcement and the bulk modulus of rubber are both included in the equilibrium equation (25), which is derived for the first time for fiberreinforced plastic plate rectangular isolators.
3.4. Solution of Compression Modulus
The method of separation of variables is used here combining with Fourier transformation to solve equation (25). It is observed that equation (25) is in the same form as the equilibrium function derived by Angeli et al. [9]. The difference is the expression of parameter . Thus, the analytical solution of pressure can be solved and expressed in the same form:
Considering the vertical normal stress of the plate, there is a relationship as follows:where P is the vertical resultant load. Combined with , the compression modulus is solved and normalized as follows:
The normalized compression modulus versus α_{1}b and βb for a square case r = 1 is plotted in Figure 5(a). For given and , the compression modulus decreases obviously as either α_{1}b or βb increases. For clarification, a parameter α is defined as . The variation of normalized compression modulus with Poisson ratio ν for different values of βb is plotted in Figure 5(b). The value of at ν = 0 means the modulus ignoring the effect of the Poisson ratio, which corresponds to the solutions by Angeli et al. [9]. It can be observed that, for smaller value of βb, is more sensitive to the change of ν. Figure 5(c) shows as a function of βb with ν = 0.3 for different values of αb. For both rigid reinforcement with αb = 0 and very flexible reinforcement with αb = 3, is very sensitive to βb. As βb is a function of K, it is obvious that the compressibility of rubber is not negligible. The curves of versus r for different values of βb are plotted in Figure 5(d), which reveals that the normalized compression modulus of square type is the largest among rectangle isolators.
(a)
(b)
(c)
(d)
3.5. Simplified Solution of Compression Modulus
To simplify the solution of the square isolator, a numerical method of weighted residuals (MWR) introduced by Xu [22] is used to solve equation (25). After several tests, a polynomial is selected as a trial function of pressure p for r near 1 and automatically satisfies boundary conditions as follows:
Replacing with in equation (25), the lefthand side of the equation does not equal to zero in the whole domain, and the residual is defined by as follows:
Using collocation method to eliminate the residuals bywhere represents a point in the definition domain, and N is the number of such points, which can be determined by the number of undetermined coefficients in above equations. The selection of is very important for the solution accuracy. After tests, five points at (0, 0), (±b, 0), and (0, ±a) are selected which can lead to three equations due to symmetry. In order to solve simplified modulus E_{cs}, an additional equation is given by assuming that satisfies equation (28). Then, C_{i} and E_{cs} can be solved as follows:where is a dimensionless parameter defined by the following equation:
For a special case of a square isolator with r equals 1, equation (36) reduces to
is plotted versus α_{1}b and βb in Figure 6(a). It has similar shape as the analytical solution plotted in Figure 5(a). It is to be noted that the simplified solution has enough accuracy only when the polynomial assumption of pressure p in equation (30) approaches its real distribution. E_{cs}/E_{c} is plotted versus γ in Figure 6(b). The error of simplified solution increases as γ increases. It is observed that when γ is smaller than 2.60, which covers most of the square fiberreinforced plastic plate isolators, the error of simplified compression modulus is smaller than 10%. A more accurate simplified solution can be obtained by improving the order of the trial function of p in equation (30). For isolators with other aspect ratios, a simplified solution can also be obtained by this method on the basis of a proper selection of the trial function of p.
(a)
(b)
3.6. Comparison with Finite Element Analysis
A special case of square cross section with typical geometry parameters is considered. Six models with different elastic moduli of the reinforcement E_{f} are established within the environment of ABAQUS. All models have a side length of 230 mm and a total thickness of 11 mm. The thickness of the reinforcement layer and the rubber layer is t_{f} = 3 mm and t = 5 mm, respectively. A constant value of Poisson’s ratio of 0.3 is used for the reinforcement material. The rubber had a shear modulus of = 0.46 MPa and bulk modulus of = 2040 MPa. E_{f} varies from 0.5 GPa to 20 GPa, and the resulting γ distribute in a relatively large range from 1.33 to 3.88, as shown in Table 1.

The Ogden model was used to describe the hyperelastic behavior of rubber. Reinforcement material was modeled by linear elastic isotropic material. The outofplane translational degrees of freedom and the rotational degrees of freedom of the reinforcement layers were restrained by reference points to keep the layer plain. The inplane translational degrees of freedom were released to allow the inplain deformation. The meshing is shown in Figure 7. A vertical pressure of 5 MPa is applied. The results of compression modulus are compared with the solutions proposed in this paper, the solutions for incompressiblerubber isolators by Tsai and Kelly [19] and the solutions for fiberreinforced isolators by Angeli et al. [9], as shown in Figure 8.
It is observed that the proposed analytical and simplified solutions always present better evaluation of than other solutions, except for model P_{6}. When α_{1}b is very small, the stiffness of the reinforcement is high, and the proposed solutions tend to underestimate the constraint effect of the reinforcement. However, E_{f} of P_{6} is as high as steel, which is almost impossible for a plastic plate employed in this isolator. In other points, the analytical curve matches well with FEA results. The error of the simplified solution becomes obvious when α_{1}b is large, but it is still smaller than previous solutions and can be optimized as discussed in Section 3.5. The proposed solutions coincide with Angeli solution at α_{1}b = 0, because the reinforcement is strictly rigid and the effect of Poisson ratio vanishes. When α_{1}b is larger than 0, Angeli solution always underestimates the modulus, showing that the effect of the Poisson ratio cannot be neglected. It can also be seen that Tsai solution always overestimates the compression modulus, showing that the effect of the compressibility of rubber need to be included.
3.7. Compression Stiffness
For an isolator with n layers, the effective compression modulus of the isolator is [23]where is the modulus of exterior layers. When n is big, is assumed to be equal to E_{c} [23]. The vertical stiffness can be calculated using equation (29):
The simplified form is calculated using equation (36):
4. Bending Analysis of the Isolator
Applying moment M about yaxis to the analytical unit, the resultant deformation is illustrated in Figure 9, in which the outofplane deformation of the reinforcement is neglected to simplify the analysis.
4.1. Equilibrium Equation of Pressure
The displacement functions of rubber in the xy plane are the same as the compression condition and are expressed by equations (1) and (2), respectively. In z direction, it is expressed by
The equilibrium of stress in rubber taking into account of rubber compressibility in a pure bending condition is derived by Angeli et al. [9], and there are two equations:
In addition, the expression of is the same as equation (9). Although the load of the bending condition is different from the compression condition, the equilibrium and boundary conditions of stresses in the reinforcement layer remain the same. Thus, equation (24) is also tenable in the bending problem. Substituting equations (9), (24), and (44) in (11), the following equation is obtained:
This equilibrium equation of p under pure bending is derived for the first time for fiberreinforced plastic plate rectangular isolators.
4.2. Solution of Bending Modulus
Equation (45) is in the same form as the equilibrium equation derived by Angeli et al. [9] except the expression of α_{1}. A similar procedure of solving this equation is conducted employing the method of separation of variables. The following equation is obtained:
The applied moment M can be expressed as follows:
Combining with , where I is the moment of inertia about the yaxis, the bending modulus is solved and normalized as follows:
The normalized bending modulus versus α_{1}b and βb for a square case (r = 1) is plotted in Figure 10(a). For given and , the bending modulus decreases obviously as either α_{1}b or βb increases. A parameter analysis is conducted of Poisson ratio ν, βb, and aspect ratio r, as plotted in Figures 10(b)–10(d). It can be observed that varies nearly linearly with the variation of Poisson ratio ν and is more sensitive to βb than to ν. The normalized bending modulus of square type is the largest among rectangle isolators of different aspect ratios.
(a)
(b)
(c)
(d)
4.3. Simplified Solution of Bending Modulus
The MWR method is also used here to achieve a simplified solution of bending modulus. After analyzing the shape of the function of E_{b} in equation (48) and conducting several tests of trial functions of pressure p, a polynomial function as shown in equation (49) is selected as a balance point between accuracy and simplification:
This function automatically satisfies the boundary condition that the stresses on the boundary equal zero. Substituting it into equation (45) and eliminating the residual by collocation method at , , and (0,0), D_{i} can be solved and finally we have the following equation:
For a square case, it can be further simplified as follows:
is plotted versus α_{1}b and βb in Figure 11(a). It has similar shape as the analytical solution plotted in Figure 10(a). E_{bs}/E_{b} is plotted versus γ in Figure 11(b). The error of simplified solution increases as γ increases. It is observed that, when γ is smaller than 3.4, the error of simplified bending modulus is smaller than 10%, which is 30% larger than the suitable range of E_{cs} and can cover most of the common square isolators. A more accurate simplified solution can be obtained by improving the order of the trial function of p in equation (49). For isolators with other aspect ratios, a simplified solution can also be obtained by this method on the basis of a proper selection of the trial function of p.
(a)
(b)
4.4. Comparison with Finite Element Analysis
The models established in section 3.6 are employed, and a rotation of 0.005 rad is applied to the reinforcement layer to analyze the bending modulus. The FEA results are normalized and compared with the proposed solutions, Tsai solution and Angeli solution, as shown in Figure 12.
The validity of the proposed solutions is confirmed by the curves. Both the analytical and simplified solutions are obviously closer to the FEA results. As E_{f} is commonly larger than 0.5 GPa, the solutions exhibit very good approximation of the modulus at the five points from the left in Figure 12. Similar to the conclusions for the compression case, Angeli solution always underestimates the bending modulus and Tsai solution always overestimates it, showing that the effects of rubber compressibility and reinforcement Poisson ratio are also important for bending modulus.
4.5. Horizontal Stiffness
The effect of the bending modulus of a single layer on the horizontal stiffness of an isolator is considered based on the theory of Haringx [24]. A simplified formula has enough accuracy to calculate the initial effective horizontal stiffness [25]:where is the bucking load of the isolator, which has a simplified form as follows:where E_{b} is calculated by the proposed solutions, and I is the inertia moment of the cross section.
The analytical solution of is obtained as follows:
The simplified solution for is
The analytical and simplified solutions for horizontal stiffness can be obtained by substituting equations (54) and (55) into equation (52), respectively.
5. Experimental Verification
5.1. Test Specimen
Eight representative specimens were tested to verify the proposed solutions. The specimens have been designed based on commonly used materials and dimensions of fiberreinforced plastic plate isolators. The detailed parameters are summarized in Table 2. Each specimen is described by aspect ratio r, elastic modulus of reinforcement E_{f}, number of rubber layers n, side length 2b, thickness of each rubber layer t, thickness of each reinforcement layer t_{f}, and shear modulus of rubber . The dimensionless parameter γ is reported in the last column.

5.2. Setup and Test Procedure
All the bearings were tested on the compressionshear testing machine in the Earthquake Engineering Research and Test Center of Guangzhou University, China. The vertical actuator can impose a maximum compression force of 10,000 kN, and the horizontal actuator is capable of a maximum force of ±500 kN.
Standard cyclic vertical and horizontal tests as shown in Figure 13 were conducted according to the method prescribed in the Chinese standard [26]. The isolators were subjected to a cyclic loading for four cycles in each test, and the data of the third cycle were used for calculation. The vertical tests were conducted under a design pressure of 5 MPa for each specimen, and the horizontal tests were loaded to a shear strain of 100% to achieve equivalent stiffness.
(a)
(b)
5.3. Results and Discussion
5.3.1. Vertical Tests
Samples of the results of vertical tests of r = 1 and r = 2 are shown in Figures 14(a) and 14(b). Specimen B3 has smaller shape factor S than B5 but shows obvious larger vertical stiffness in Figure 14(a), which proves the increase of γ results in reduction in compression modulus as derived by equations (29) and (36). The comparison between B7 and B8 shows that the increase of results in increase of compression stiffness, despite a resultant increase of γ, which indicates has a larger effect.
(a)
(b)
5.3.2. Horizontal Tests
A sample of the results of horizontal tests of B7 is shown in Figures 15(a) and 15(b). It can be observed that there is a big difference between the stiffness in x and y directions. This phenomenon cannot be explained by traditional stiffness formula of . But it is reasonable in equation (48) in which different aspect ratio r leads to different bending modulus. The damping ratio of the isolator in x direction is 4.0% and in y direction is 3.0%.
(a)
(b)
5.3.3. Discussion
The test results and analytical and simplified solutions are listed in Table 3 for r = 1, including pressure σ_{z}, experimental vertical and horizontal stiffness K_{v.exp} and K_{H.exp}, analytical vertical and horizontal stiffness K_{v.ana} and K_{H.ana}, and simplified vertical and horizontal stiffness K_{v.sim} and K_{H.sim}. Test results and analytical solutions are listed in Table 4 for r = 2. The vertical stiffness was tested under design pressure of each specimen, and the relativity between vertical stiffness and pressure is not discussed.


As can be seen from Tables 3 and 4, results by the proposed analytical and simplified solutions agree well with the experimental results for both cases. The relative differences of analytical solutions for vertical and horizontal stiffness are less than 8.7% and 11.7%, and the relative difference of simplified solutions for both vertical and horizontal stiffness is less than 11.7%.
6. Conclusion
The analytical solutions for compression and bending modulus of rectangular fiberreinforced plastic plate isolators are proposed, and simplified solutions are also derived for a special case of square shape. The solutions show advantages to the existing solutions and agree well with experimental tests.
Based on the analytical analysis, finite element analysis, and experimental results, the following conclusions can be made:(1)The effects of the Poisson ratio of reinforcement and compressibility of rubber cannot be neglected for both compression and bending modulus. The proposed solutions of compression and bending modulus can be used to calculate the vertical and horizontal stiffness of this type of isolators.(2)Simplified solutions for compression and bending modulus are derived by a method of weighted residuals. The accuracy is controlled by the selection of trial function of pressure p. For square isolators, the error of simplified solution is less than 10% for compression modulus when dimensionless parameter γ is less than 2.6 and for bending modulus when γ is less than 3.4. The applicable range covers most types of fiberreinforced plastic plate isolators.(3)Finite element analysis results show that neglecting the compressibility of rubber overestimates both compression and bending modulus and neglecting the effect of the Poisson ratio underestimates both the moduli. The proposed solutions include both of the effects and have obvious advantage over the solutions by previous papers.(4)The validity of the proposed solutions is verified by a series of experimental tests. The error of analytical solution for vertical stiffness is less than 8.7% and for horizontal stiffness is less than 11.7%. The error of simplified solution for vertical stiffness is less than 11.7% and for horizontal stiffness is less than 8.0%.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by National Key R&D Program of China under Grant no. 2017YFC0703600, the Program for Changjiang Scholars under Grant no. IRT13057, in part by program for Innovation Team of Department of Education of Guangdong Province under Grant no. 2016KCXTD016, in part by the Program for Ram City Scholar under Grant no. 1201541630, and Guangdong Special Program under Grant no. 2014TX01C141.
References
 X. L. Ge, L. X. Zhu, and W. Yu, “Overview of researches on disaster prevention of buildings in village and town,” Journal of Natural Disasters, vol. 23, no. 2, pp. 8–14, 2014, in Chinese. View at: Google Scholar
 W. L. Cao, Z. Y. Zhou, Q. Wang et al., “Experimental study on base vibration isolation and antiseismic masonry structure in rural areas by shaking table test,” Journal of Vibration and Shock, vol. 30, no. 11, pp. 209–213, 2011, in Chinese. View at: Google Scholar
 F. L. Zhou, Vibration Control of Engineering Structures, Seismological Press, Beijing, China, 1997, in Chinese.
 J. M. Kelly, “Analysis of fiberreinforced elastomeric isolators,” Journal of Seismology and Earthquake Engineering, vol. 2, no. 1, pp. 19–34, 1999. View at: Google Scholar
 J. M. Kelly, “Seismic isolation systems for developing countries,” Earthquake Spectra, vol. 18, no. 3, pp. 385–406, 2002. View at: Publisher Site  Google Scholar
 J. M. Kelly and S. M. Takhirov, “Analytical and experimental study of fiberreinforced strip isolators,” PEER, Berkeley, CA, USA, 2002, PEER Report, vol. 11. View at: Google Scholar
 H. ToopchiNezhad, M. J. Tait, and R. G. Drysdale, “Testing and modeling of square carbon fiber‐reinforced elastomeric seismic isolators,” Structural Control and Health Monitoring, vol. 15, no. 6, pp. 876–900, 2010. View at: Publisher Site  Google Scholar
 J. M. Kelly and A. Calabrese, “Mechanics of fiber reinforced bearings,” PEER, Berkeley, CA, USA, 2012, PEER Report, vol. 2. View at: Google Scholar
 P. Angeli, G. Russo, and A. Paschini, “Carbon fiberreinforced rectangular isolators with compressible elastomer: analytical solution for compression and bending,” International Journal of Solids and Structures, vol. 50, no. 2223, pp. 3519–3527, 2013. View at: Publisher Site  Google Scholar
 M. Pauletta, A. Cortesia, and G. Russo, “Rollout instability of small size fiberreinforced elastomeric isolators in unbonded applications,” Engineering Structures, vol. 102, pp. 358–368, 2015. View at: Publisher Site  Google Scholar
 P. Tan, K. Xu, B. Wang, C. Chang, H. Liu, and F. Zhou, “Development and performance evaluation of an innovative lowcost seismic isolator,” Science China Technological Sciences, vol. 57, no. 10, pp. 2050–2061, 2014. View at: Publisher Site  Google Scholar
 K. Xu, P. Tan, B. Wang et al., “Experimental comparative study of both types of novel simple isolators,” Earthquake Engineering & Engineering Dynamics, vol. 34, no. 6, pp. 11–18, 2014, in Chinese. View at: Google Scholar
 Y. Rocard, “Note sur le calcul des propriétès élastiques des supports en caoutchouc adhérent,” Journal de Physique et le Radium, vol. 8, no. 5, pp. 197–203, 1937. View at: Publisher Site  Google Scholar
 A. N. Gent and P. B. Lindley, “Compression of bonded rubber blocks,” Proceedings of the Institution of Mechanical Engineers, vol. 173, no. 3, pp. 111–122, 1959. View at: Publisher Site  Google Scholar
 A. N. Gent and E. A. Meinecke, “Compression, bending, and shear of bonded rubber blocks,” Polymer Engineering and Science, vol. 10, no. 1, pp. 48–53, 2010. View at: Publisher Site  Google Scholar
 C. G. Koh and H. L. Lim, “Analytical solution for compression stiffness of bonded rectangular layers,” International Journal of Solids and Structures, vol. 38, no. 3, pp. 445–455, 2001. View at: Publisher Site  Google Scholar
 H.C. Tsai and W.J. Pai, “Simplified stiffness formulae for elastic layers bonded between rigid plates,” Engineering Structures, vol. 25, no. 11, pp. 1443–1454, 2003. View at: Publisher Site  Google Scholar
 J. M. Kelly, “Behavior of multilayered bearings under compression and bending,” in EarthquakeResistant Design with Rubber, pp. 131–159, Springer, London, UK, 1997. View at: Publisher Site  Google Scholar
 H.C. Tsai and J. M. Kelly, “Stiffness analysis of fiberreinforced rectangular seismic isolators,” Journal of Engineering Mechanics, vol. 128, no. 4, pp. 462–470, 2002. View at: Publisher Site  Google Scholar
 J. M. Kelly and A. Calabrese, “Analysis of fiberreinforced elastomeric isolators including stretching of reinforcement and compressibility of elastomer,” Ingegneria Sismica, vol. 30, no. 3, pp. 5–16, 2013. View at: Google Scholar
 S. Pinarbasi and Y. Mengi, “Elastic layers bonded to flexible reinforcements,” International Journal of Solids and Structures, vol. 45, no. 34, pp. 794–820, 2008. View at: Publisher Site  Google Scholar
 C. D. Xu, “Weighted residuals method for solid mechanics,” Mechanics in Engineering, vol. 2, no. 4, pp. 12–20, 1980, in Chinese. View at: Google Scholar
 H. C. Tsai, “Compression stiffness of circular bearings of laminated elastic material interleaving with flexible reinforcements,” International Journal of Solids & Structures, vol. 43, no. 1112, pp. 3484–3497, 2006. View at: Publisher Site  Google Scholar
 J. A. Haringx, “On highly compressible helical springs, and rubber rods, and their application for vibrationfree mountings,” Philips Research Reports, vol. 4, pp. 49–80, 1949. View at: Google Scholar
 Architectural Institute of Japan, Recommendation for the Design of Base Isolated Building, Marozen Corporation, Tokyo, Japan, 1993.
 Standardization Administration of the People’s Republic of China, GB/T 20688.12007: Rubber Bearings—Part 1: Seismic Protection Isolators Test Methods, SAC, Beijing, China, 2007.
Copyright
Copyright © 2019 Han Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.