Abstract
As an important support member in the structural system of coastal bridges, the frictional slip and the rubber aging of laminated rubber bearings will affect the service safety of the overall structure in earthquakes. In order to investigate the mechanical properties aging law of the rubber bearings considering frictional slip in the coastal bridges, a frictional slip experiment was carried out on the laminated rubber bearings. Moreover, the influence of rubber aging on the mechanical properties of the bearings with various shape coefficients was analyzed by the finite element method during the 100 years of service life of bridges. The results indicate that (1) the horizontal and vertical stiffness of the bearing increase linearly with the aging time of the rubber. The amplification of the bearing stiffness also grows with the shape coefficient of the bearing. (2) The frictional slip initiation displacement of the bearing grows with the rubber aging time. Furthermore, the larger the shape coefficient of the bearing is, the more the frictional slip initiation displacement of the bearing increases. (3) With the increase of the aging time, the equivalent viscous damping ratio of the bearing continues to increase and more energy is consumed by frictional slip. For the bearing with the shape coefficient of 16.38, the equivalent viscous damping ratio at 100 years of rubber aging time is 1.17 times higher than that of the initial state of the bearing, and 33.21% more energy is consumed through frictional slip. Given that the marine environment accelerates rubber aging and changes the mechanical properties, the effects of the frictional slip and rubber aging properties of the laminated rubber bearings on the seismic dynamic response of bridges should be considered in the seismic design of coastal bridges.
1. Introduction
With the extensive construction of coastal roads in countries around the world, laminated rubber bearings have been widely used in coastal bridges [1]. The laminated rubber bearing is constituted by the evenly distributed rubber layers bonded with steel plates, which is the important loading transmission component connecting the girder with the pier in bridges. Under the coupling effect of the coastal environment and service load, the laminated rubber bearing suffered the problems such as steel plate corrosion and rubber aging [2]. At the same time, the laminated rubber bearing is usually placed directly between the girder and the pier without any anchoring in the small- and medium-span bridges. In earthquakes, the unanchored bearings will effectively release the seismic forces and protect the pier by larger shear deformation and frictional slip [3]. In the 2008 Wenchuan Earthquake, there were only 56 piers in the total 2316 piers of simply supported girder bridges that suffered damage in the VII–XI degrees region and the damage rate was 2.4%. However, the damage rate of the bearing reached 16.6% and the majority is laminated rubber bearings [4].
As the laminated rubber bearing is widely used in high-intensity areas in the world, the mechanical properties of the laminated rubber bearing were extensively studied in earthquakes. McDonald et al. [5, 6] analyzed the effects of vertical forces, support surface roughness and sliding speed on the frictional slip characteristics between the rubber bearing and the supporting member and found that the frictional slip of the bearing has a significant impact on its mechanical properties. Konstantindis et al. [7] investigated the shear performance of the unanchored laminated rubber bearing by the experiment; the results showed that the ultimate shear strain of the bearing can reach 150%–250%, after that the bearing will slip. Steelman et al. [8] studied the frictional slip characteristics of a single anchorage laminated rubber bearing (only the top steel plate is anchored) and concluded that the energy dissipation was the most when the friction coefficient between the bearing and the concrete support was 0.25–0.50. The experiment results provided evidence for the amendment of the bearing provisions in the AASHTO Code. After the 2008 Wenchuan Earthquake, Li et al. [9, 10] and Wang et al. [11] confirmed that the frictional slip of the bearings provided seismic isolation effect but increased the risk of girder falling in the simple support girder bridge by the shaking table experiment.
In addition, Itoh et al. [12] found that salt spray, acid mist, and ultraviolet radiation would influence the mechanical properties of the rubber by the durability study on the rubber materials in the coastal environment. The rubber aging will be accelerated by the complex marine environment, which results in the loss of initial elastic properties and the decrease of the loading transferring performance of the laminated rubber bearing [13–17]. Ab-Malek and Stevenson [18] found the mechanical properties deteriorated after testing the rubber materials immersed in seawater for about 42 years. Brown and Butler [19–21] investigated the rubber material in the marine environment from the perspectives of thermal aging and natural aging, respectively. It was found that the mechanical properties of rubber deteriorated after aging; for example, the tensile strength and elongation at break of rubber were both decreased. Yoshito et al. [22] conducted a durability test on the rubber block and confirmed that the marine environment had a great impact on the mechanical properties of the rubber.
Moreover, based on the rubber aging test, Kato et al. [23] conducted a thorough inspection on the material aging characteristics of the rubber bearings used for about 40 years on the offshore Pelham Bridge and found that the horizontal and vertical stiffness of the bearings increased by about 10% and 6%, respectively. Han et al. [24, 25] determined the basic performance and vertical compression deformation of the rubber bearings after 60 aged years in coastal bridges by conducting experiments on the temperature characteristics of rubber bearings [26]. Tlpakidis [27, 28], Takenaka [29], Gu and ltoh [30, 31], and Chou et al. [32, 33] studied the performance deterioration of rubber bearings from a single marine environmental factor such as thermal aging, ultraviolet radiation, and seawater corrosion, respectively. The results indicated that the rubber aging was accelerated by the complex marine environment, which led to the degradation of the mechanical properties of the rubber bearing and seriously influenced the service time and seismic isolation performance of the rubber bearing.
From the aforementioned literature review, it can be noticed that the frictional slip of the laminated rubber bearing has an important influence on the seismic response of the bridge, and the aging of the rubber material in the coastal environment will further affect the seismic isolation effect of the bearing. The existing research is mainly focused on the mechanical properties of the rubber material and the bearing in the coastal environment. However, little attention has been paid to the effect of the rubber material aging on the mechanical performance of the laminated rubber bearings in coastal bridges considering the frictional slip. Therefore, the effects of the rubber materials aging on the stiffness, frictional slip initiation displacement, energy dissipation capacity, and equivalent viscous damping ratio of the laminated rubber bearing in coastal bridges are analyzed based on the reverse loading test of the laminated rubber bearing and the aging performance theory of the rubber material. Moreover, the influence of the aging of rubber materials on the mechanical performance of the rubber bearings is investigated, which will provide the reference and theory guide for the seismic design of coastal bridges.
2. Mechanical Properties of the Laminated Rubber Bearing after Material Aging
The stiffness of the laminated rubber bearing refers to the ability to resist elastic deformation under loading. According to the multinational specification [34, 35], the vertical stiffness of the bearing is calculated by the following equation [36, 37]:where Kv is the vertical stiffness of the bearing, Ec is the rubber modified elastic compression modulus, A is the actual vertical load bearing area, Eap is the apparent elastic modulus of the rubber, E∞ is the bulk modulus of the rubber, E0 is the elastic compressive modulus of the rubber, k is the correction factor, S1 is the first shape coefficient, and G is the shear modulus of rubber.
The horizontal shear stiffness of the bearing is calculated by the following equation:where Kh is the horizontal stiffness and Tr is the total rubber layer thickness.
Due to the complex molecular properties of the rubber material, Mooney [38] has proposed the large elastic deformation theory of rubber based on the strain energy functions in 1940 to more accurately describe the elastic constitutive behaviors of rubber. The theory assumes that (1) the rubber is isotropic and incompressible (0.48 ≤ μ ≤ 0.5) superelastic body; (2) the rubber follows Hooke’s law in shear deformation and superimposed shear in a plane perpendicular to the unidirectional tensile or compressive axis. Based on the assumption, the rubber strain energy function W is defined here as follows [39–41]:where λ1, λ2, and λ3 are side elongations of the rubber cube and C10 and C01 are the positive definite constants for rubber materials.
From the mathematical point of view, Rivlin [42] has obtained the general expression for the strain invariant considering the general form of the strain energy function, as follows:
Based on the incompressibility of the rubber material and the volume invariance in deformation, the relation I3 = 1 is introduced, so that only two of the three strain invariants are independent, as follows:where I1 and I2 are determined by the stretching ratio of the three main axes. Thus, the strain energy function corresponding to the isotropic incompressible material can be expressed as = (I1, I2, I3). Then, Mooney’s rubber strain energy function is described as follows:
Equation (10) is the Mooney–Rivlin model widely used in the study of the deformation of the elastomer, which reflects the mechanical properties of rubber.
However, on coastal bridges, rubber aging is extremely serious due to the long-term complex and harsh marine environment. Based on the thermal aging experiment of rubber materials, the time-varying laws of the constants C10 and C01 in the Mooney–Rivlin model for rubber materials are obtained as shown in equations (11) and (12) [43, 44]. Derived from equations (11) and (12), the material constants for every 10 aging years are shown in Table 1, and the variation curves are shown in Figure 1.where n is the actual aging time and the unit is year.

Then, the expressions between the mechanical properties and aging time of the laminated rubber bearings are obtained as follows [36]:where Kv (t) is the vertical stiffness corresponding to the aged time t, Kv (0) is the initial vertical stiffness, Kh (t) is the horizontal stiffness corresponding to the aged time t, and Kh (0) is the initial horizontal stiffness.
3. Mechanical Performance Experiments of the Laminated Rubber Bearings
As the laminated rubber bearings are always put on the piers without any connection with the girders in small and medium-span bridges, the frictional slip performance of the bearings was investigated under various axial pressures by the cyclic loading test. The test results are the basis of the investigation of the seismic mechanical properties of the aged laminated rubber bearing in the coastal environment. The loading device is shown in Figure 2 and the test conditions are shown in Table 2.

In order to realistically reflect the working condition of the bearing, the top surface of the bearing is anchored with the loading device through the connection plate under a constant vertical load. The bottom rubber layer of the bearing will slide on the support surface. The deformation of the bearing is shown in Figure 3 during the loading process.

The relationship between the friction force and the loading displacement is shown in Figure 4. It can be found that the friction will convert from static friction to sliding friction with the increase of the displacement. The friction coefficient is inversely proportional to the compressive stress (as shown in Figure 5) and decreases from 0.31 to 0.15 when the vertical pressure grows from 4 MPa to 10 MPa.


To investigate the effect of rubber aging on the mechanical properties, the finite element models of the laminated rubber bearing were established in ABAQUS. The bearing was put on the supporting rigid block, as shown in Figure 6(a).

(a)

(b)
The eight-node linear brick reduced integral three-dimensional solid element was adopted to simulate the steel plate layer of the bearing whose elastic modulus is 2.0E5 MPa and Poisson’s ratio is 0.3. To precisely describe the mechanical properties, the Mooney–Rivlin model was adopted to simulate the rubber material, and the intrinsic constants of the material were taken from Table 1 corresponding with aging time. The rubber layers were meshed by the eight-node linear brick reduced integral hybrid three-dimensional solid element. Additionally, a mesh convergence analysis was carried out to explore the appropriate mesh size. The results suggested that the bearing with a mesh size of about 20 mm can effectively simulate the performance of the studied bearing with an appropriate time-consuming. The tied constraint was used between each steel plate layer and rubber layer.
To model the surface interaction between the bottom rubber layer of the bearing and the top layer of the supporting rigid block, the classical Coulomb friction model and the contact pressure model were designated the tangential and normal directions to the contact surfaces, respectively. Based on previous studies [3], the friction coefficient between the contact surface of the rubber and steel plate was set as 0.3, which is also adopted in the code [45]. All nodes on the top endplate of the bearing were constrained to a reference point at the centroid that was restrained against rotation but allowed to translate vertically and in one direction laterally. Similarly, all nodes on the bottom surface of the supporting rigid block were constrained against translation and rotation in all degrees of freedom. The deformation of the bearing is illustrated in Figure 6(b).
The friction sliding hysteresis curves of the bearing obtained from the finite element model are compared with the test results, as shown in Figure 7 (taken specimen B-1 under a vertical pressure of 6 MPa as an example). It can be found that the simulation results agree well with the experimental results. In addition, the damping ratio, energy dissipation, equivalent stiffness, and friction coefficient are compared with the test results for different equivalent shear deformations after frictional slip of the bearing, as shown in Table 3. At the early loading stage, when the shear deformation is 100%, the bearing does not slip, and the energy dissipation is minor, so there is a numerical error in the damping ratio and energy dissipation, and the deviation between the experimental and numerical results is 27%. However, after the shear deformation beyond 100%, the energy dissipation of the bearing keeps increasing, and the deviations of the damping ratio and the energy dissipation between the experimental and numerical results are within 10%. At this time, the deviation of the equivalent stiffness is only 2.5% at 200%, and the deviation of the friction coefficient is 9.3% at 300%. According to the experimental and numerical results shown in Figure 7 and Table 3, it can be verified that the finite element model could better reflect the mechanical properties of the bearings under the horizontal recycling load.

(a)

(b)

(c)

(d)
4. Influence of Rubber Aging on the Mechanical Properties of Laminated Rubber Bearings
In earthquakes, the frictional slip of the laminated rubber bearing and the aging time of the rubber material will have a significant impact on the seismic performance of the coastal bridge [3, 37]. Therefore, the influence of the rubber aging time on the horizontal and vertical stiffness, frictional slip initiation displacement, and energy dissipation capacity of laminated rubber bearings with four shape coefficients is analyzed, and the bearing parameters are shown in Table 4.
4.1. Effect of Rubber Aging on the Stiffness of Laminated Rubber Bearings with Four Shape Coefficients
Figure 8 shows the relationship between the rubber aging time and the horizontal and vertical stiffness of the bearings with four shape coefficients. It can be found that the horizontal and vertical stiffness of the bearing is increased linearly with the rubber aging time.

(a)

(b)
For the specimen P-1 with a shape coefficient of 16.38, the maximum difference of the horizontal stiffness between the aging 100 years and the initiation reaches 34.3%, and the maximum difference of the vertical stiffness is 23.48%. Moreover, the stiffness growth rate is increasing with the aging time for the larger shape coefficient of the bearing.
4.2. Effect of Aging Time on the Frictional Slip Initiation Displacement of the Bearings
The friction slip initiation displacement of the laminated rubber bearing is an important reference after which the stage will convert from the controlled static friction to the uncontrolled sliding friction in a horizontal earthquake, which has a great impact on the displacement control of the girder. Therefore, in the service life of the bridge, the effect of rubber aging time on the frictional slip initiation displacement of the bearings with four shape coefficients is analyzed when the shear deformation reaches 300%, as shown in Figure 9. It can be found that the initial sliding displacement of the bearing increases with the aging time. At the same time, the larger shape coefficient bearing also has a greater initiation displacement. Therefore, for the bridges with laminated rubber bearings located at the coastal areas, more attention should be paid to controlling the displacement of the girders in earthquakes and preventing the girder from falling. At the same time, for the bearing with a smaller shape coefficient, the initial sliding displacement is also smaller. Namely, for the bearing with a thicker rubber layer, the displacement of the girder is smaller when the bearing enters the sliding state. Therefore, when the high bearings are used in the bridge, more attention should be paid to the seismic damage due to the bearing slip.

4.3. Effect of Aging Time on the Frictional Slip Performance of the Laminated Rubber Bearings
The force and displacement hysteresis curves of the bearings are compared in the initial state, the aging for 50 years, and the aging for 100 years. Taking the shear deformation of 300% as an example (as shown in Figure 10), the horizontal force and displacement hysteresis curve of the bearing is fuller with the rubber aging time increasing. At each rubber aging stage, the larger the shape coefficients of the bearings, namely, the thinner the rubber layer, the fuller the force and displacement hysteresis curve of the bearing under horizontal load, as shown in Figure 11. It is mainly because the thinner the rubber layer is, the greater the horizontal stiffness is, the earlier the bearing enters the frictional slip stage. Therefore, for the bridges in the marine environment, it is necessary to consider the influence of the bearing friction slip on the seismic performance of the bridge to avoid the girder falling when the shape coefficient of the laminated rubber bearing is smaller.


(a)

(b)

(c)
From the relationship between the aging time of rubber material and the energy consumption by the frictional slip of the bearing (as shown in Figure 12), it can be found that the dissipated energy by the friction slip of the bearing will be more with the aging time increasing. For the bearing with a 16.38 shape coefficient, the rubber material aging for 100 years dissipated 33.21% more energy than that of the initial state by frictional slip. When the rubber aging time is 100 years, the bearing with a shape coefficient of 16.38 dissipated 41.87% more energy by frictional slip than the bearing with a shape coefficient of 8.73.

Meanwhile, from the relationship curve between the rubber aging time and the equivalent viscous damping ratio of the bearing as shown in Figure 13, it can be concluded that the equivalent viscous damping ratio of the bearing increases with the aging time. For the bearings with the shape coefficients of 16.38 and 8.73, the equivalent viscous damping ratios of the rubber 100 years aging time are 1.17 times and 1.40 times that of the bearing initial state. Moreover, the shape coefficient has a large effect on the equivalent viscous damping ratio of the bearing. When the rubber aging time reaches 100 years, the equivalent viscous damping ratio of the bearing with the shape coefficient of 16.38 is 1.42 times the bearing with the shape coefficient of 8.73.

5. Conclusions
Earthquake experience shows that the frictional slip of the bearing and the aging of the rubber material will have a significant impact on the seismic mechanical properties of the laminated rubber bearings in coastal bridges. Based on the frictional slip experiment, the effect of the rubber material aging on the mechanical performance of the laminated rubber bearings in coastal bridges was studied, and the main conclusions are as follows:(1)The horizontal and vertical stiffness of the laminated rubber bearing increases linearly with the growth of rubber aging time. For the bearing with a shape coefficient of 16.38, the maximum differences between the horizontal and vertical stiffness of the bearing after 100 years of aging are 34.3% and 23.8% of the initial state, respectively. Moreover, the larger the shape coefficient is, the greater the bearing stiffness increases with the aging time of rubber material.(2)Considering the frictional slip, the initial sliding displacement of the bearing increases with the growth of the aging time. Moreover, the larger the shape coefficient of the bearing is, the quicker the initial displacement of the frictional slip increases.(3)With the growth of aging time, the dissipation energy is also increased by the frictional slip. For the bearing with a shape coefficient of 16.38, the dissipation energy by frictional sliding of the bearing with the 100 years of aging is 33.21% that of the bearing in initial state. When the aging time of rubber materials reaches 100 years, the bearing with the shape coefficient of 16.38 dissipates 41.87% more energy by frictional sliding than that of the bearing with the shape coefficient of 8.73.(4)With the growth of aging time, the equivalent viscous damping ratio of the bearing also increases. For the bearing with the shape coefficients of 16.38 and 8.73, the equivalent viscous damping ratio of the bearing with 100 years aging time is 1.17 times and 1.40 times that in the initial state. Moreover, the larger the shape coefficient of the bearing is, the quicker the equivalent viscous damping ratio increases.(5)As the rubber aging has a significant impact on the seismic response of the laminated rubber bearings, the effect of the rubber aging time and the shape coefficient of the bearings on the service life in the seismic areas with significant marine erosion environment should be considered.
Data Availability
The experimental data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this study.
Acknowledgments
This work was financially supported by the National Natural Science Foundation of China (Grant nos. 51408009 and 51608010), Fundamental Research Funds for the Beijing’s Universities (110052971921/062), and Yuyou Talent Project of North China University of Technology (Grant no. XN012/044).