Abstract

The cross-sea bridges play an important role to promote the development of regional economy. These bridges located in earthquake-prone areas may be subjected to severe earthquakes during their lifetime. Group pile foundations have been widely used in cross-sea bridges due to their structural efficiency, ease of construction, and low cost. This paper investigates the seismic performance of bridge pile foundation based on the seismic fragility analysis. Based on the analysis platform OpenSees, the three-dimensional finite model of the bridge pile foundation is developed, where the pile-water interaction is replaced by the added mass method, nonlinear p-y, t-z, and q-z elements are used to simulate pile-soil interaction, and the displacement of the surface ground motion due to seismic excitations is applied on all spring supports. The seismic fragility curves of the bridge pile foundation are generated by using the earthquake records recommended by FEMA P695 as input motions. The curvature ductility based fragility curves are obtained using seismic responses for different peak ground accelerations. The effects of pile-water interaction, soil conditions, and different types of ground motions on the bridge pier fragilities are studied and discussed. Seismic fragility of the pier-group pile system shows that Sec C (the bottom section of the pier) is the most vulnerable section in the example fluid-structure-soil interaction (FSSI) system for all four damage LSs. The seismic responses of Sec E (a pile section located at the interface of the soil layer and water layer) are much lower than other sections. The parameter analysis shows that pile-water interaction has slight influence (less than 5%) on the fragility curves of the bridge pier. For the bridge group pile foundations considering the fluid-pile-soil interaction, PNF may induce larger seismic response than far-field (FF) and no-pulse near field (NNF). The bridge pile foundation in stiff soil is most vulnerable to seismic damage than soft condition.

1. Introduction

With the development of economy and advancement of technology, more and more cross-sea bridges have been constructed around the world in recent decades. The cross-sea bridges located in earthquake-prone areas, such as the eastern coast of China or the western coast of the USA, may be subjected to earthquake action. The strong ground motion can destruct the cross-sea bridges, which will cause great economic loss and confusion of society. Therefore, seismic design of cross-sea bridges has become an urgent research topic. In recent years, pile group foundation has been widely used in cross-sea bridges or brides in deep water due to the structural efficiency, ease of construction, and low cost [1], where most of these foundations have their piles partially or totally submerged in deep water. It should be noted that the dynamic behavior of structures in deep water requires special considerations due to the interaction of water with structure [24]. The earthquake-induced hydrodynamic force on the structure may increase the structural responses. The experiment conducted by Wei et al. [5] also indicated that fluid-structure interaction had significant effects on the dynamic response of bridge pile foundations submerged in water. In addition to the fluid-structure interaction (FSI), the soil-structure interaction (SSI) can play an important role in the seismic responses of bridges [610]. Consequently, it is important to evaluate the seismic responses of cross-sea bridges considering SSI. This study mainly focuses on the seismic responses of bridge pile group foundation considering fluid-structure-soil interaction (FSSI).

The earthquake-induced hydrodynamic forces on a single cylinder [1113] and the seismic responses of cylinder surrounded by water [1416] have been investigated widely in recent decades. The results indicated that effect of free surface waves on the earthquake-induced hydrodynamic force can be neglected and effect of water compressibility is quite important when the dimensionless frequency is greater than unity. Wang et al. [17] developed an accurate time-domain model to replace the compressible water in the water-cylinder interaction analysis. It should be noted that the incompressible water can be replaced by an added mass [3]. Recently, a series of simplified formulas to evaluate the added mass of a single cylinder are presented by assuming the cylinder rigid, such as the circular cylinder [18, 19], the elliptical cylinder [20, 21], the round-ended and rectangular cylinders [22], and the circular tapered cylinders [23]. In addition, the simple added mass models for a flexible cylinder vibrating in water were also proposed by Han and Xu [24], Yang and Li [25], and Wang et al. [26]. A few studies were carried out to investigate the hydrodynamic interaction of multiple cylinders subjected to earthquake loading [27, 28]. The results indicated that the effects of interaction of multiple cylinders with water on the seismic responses of the cylinder increase as the relative distance decreases. Wang et al. [29] developed an accurate finite element model to evaluate the earthquake-induced hydrodynamic forces on multiple cylinders. The method presented by Wang et al. [29] is used to calculate the coefficients of the pile group in the present study.

It is well known that the SSI can have significant influence on the seismic responses of structures. Generally, two methods including the integrated finite element model [3032] and the substructure model [3335] are always adopted to model the SSI on the dynamic responses of offshore structures subjected to earthquake action. The integrated finite element model is inconvenient to use the model in structural seismic design because of its long time costs in modeling and calculation. The substructure models are more applicable to be used in structural seismic design than the above model. In the substructure model, the pile-soil interaction can be replaced by linear springs [35] or nonlinear springs based on p-y curves [36], where all of the soil springs are attached to the pile. There are mainly two different methods for applying ground motion when the pile-soil interaction replaced by soil springs. One method is uniform seismic excitations [37], where the surface ground motion is applied to all soil spring boundaries from the bedrock and to the seabed. The other method is nonuniform seismic excitations [38], where different motion obtained by one-dimensional free-field analysis is applied at each support of soil springs. The research conducted by [39] indicated that the offshore structures will be safe by using the uniform seismic excitations. Therefore, the nonlinear springs based on p-y curves and the uniform seismic excitations are adopted in the present study.

However, a few studies have been carried out to analyze the seismic response of cross-sea bridges considering FSSI. The main objective of this study is to obtain fragility curves for a typical bridge pile group foundation considering FSSI. Fragility curves are an appropriate tool to evaluate the seismic vulnerability of structures and to estimate the probability of exceeding a certain damage level for a specific seismic ground shaking intensity [40]. In this paper, fragility curves for the bridge pile group foundation are estimated by an analytical approach, specifically by conducting an incremental dynamic analysis. The effects of water, soil, and different kinds of ground motions on the seismic behavior of the bridge pile group foundation are quantitatively analyzed.

1.1. Numerical Modeling

Based on a sea-crossing bridge in the southeast coastal areas of China, a bridge pier on high-cap pile group foundations is selected as an example. The bridge pier comprises a concrete hollow pier, a pile cap, and 15 piles. Figure 1 shows the elevation view and component cross sections of the bridge pier. A three-dimensional finite element model of the pier-pile system considering FSSI is developed using the open-source software platform OpenSees [41], as shown in Figure 2. The detailed description of the pier, pile, pile-soil interaction, and pile-water interaction is provided herein.

1.2. Bridge and Pile

The heights of the pier, pile cap, and bridge pile foundation are 38 m, 5 m, and 80 m, respectively. The upper 40 m of pile group is submerged in seawater, and the lower 40 m is embedded in subsea soils. Figure 1 shows the cross sections ranged from A-A to E-E and locations of reinforcing bars of the pier and the pile foundations. Poisson’s ratio of steel and concrete is 0.3 and 0.2, respectively. The rebar’s diameter of the pier and the piles are 28 mm and 36 mm, respectively.

Figure 3 shows the material behavior of steel of HRB335 and HRB400 and the concrete of C40 and C30. The C40 concrete is used for the pier column and pile cap, while C30 concrete is used for pile group foundations The standard compressive strengths of the concrete are 26.8 and 20.1 MPa for the pier column and the piles, respectively. The yield strengths of longitudinal bars used for the pier column and pile foundations are 335 MPa and 400 MPa, respectively, while the yield strengths of stirrups are 335 MPa. The torsional and moments of inertial in the two in-plane directions are 57283.3 m4, 10953.6 m4, and 46329.7 m4 for the pile cap, respectively. The section area and elastic modulus of concrete for the pile cap are 519.9 m2 and 32.5 GPa, respectively. The moments of inertia of the pile cap cross section about the x-axis is 10953.6 m4, while about the y-axis is 46329.7 m4.

The bridge pier is modeled with the displacement-based beam-column element, which is based on the displacement formulation, and considers the spread of plasticity along the element. Based on the fiber model in the OpenSees program, the hysteretic behavior of bridge column specimens is investigated by Han et al. [42]. The results indicated that the fiber model agrees well with the experimental results. In the present study, Concrete 01 from the material library in OpenSees proposed by Scott et al. [43] and Reinforcing Steel from the material library in OpenSees proposed by Chang and Mander [44] are selected to simulate the concrete behavior and the reinforcement behavior. Figure 2 shows the finite element model for the connections between the piles and surrounding soil layers. For each pile, the bottom 40 m embedded in the subsea soil layers consists of 20 pile elements and 21 pile nodes and the upper 40 m submerged in seawater consists of 20 pile elements and 20 pile nodes.

1.3. Pile-Soil Interaction

The soil is assumed to be soft clay. The nonlinear p-y, t-z and q-z elements developed by Boulanger et al. [45, 46] are used to describe the interaction of pile with soil. Figure 4 shows the characteristics of the nonlinear p-y spring components, which is conceptualized as consisting of elastic, plastic, and gap components in series. The gap component consists of a nonlinear closure spring in parallel with a nonlinear drag spring. In addition, a dashpot in parallel with the elastic component is used to model the radiation damping.

The constitutive behavior of the abovementioned springs can be defined by uniaxial materials PySimple1, TzSimple1, and QzSimple1 in OpenSees, where parameters pult, y50, and Cd are required for PySimple1, tult and z50 for TzSimple1, and qult and z50 for QzSimple1. Distributed nonlinear zero-length elements with p-y, t-z and q-z materials are created at the pile nodes. The properties of each p-y, t-z, and q-z element are derived based on the properties of each layer at the corresponding depth. It should be noted that pult and y50 proposed by [47] and tult, qult, and z50 based on the API [48] recommendations are adopted in the present study. This study uses Cu = 40 kPa,  = 0.01, J = 0.375, and γ = 10 kN/m3, where Cu, , J, and γ denote the undrained shear strength, the strain corresponding to one-half the maximum stress on laboratory undrained compression tests of undisturbed soil, an empirical dimensionless constant, and the submerged unit weight.

The spring nodes and the fixed nodes with three translational degrees of freedom are created, which are initially fixed in all degrees of freedom. The spring nodes of the zero-length elements share equal degrees of freedom in all translational directions with the pile nodes. The nodes on the other end are subjected to the ground motion of the seismic loading. The uniform seismic excitations are adopted in the present study. Therefore, the displacement of the surface ground motion is applied on all spring supports.

1.4. Pile-Water Interaction

The fluid damping effects on structural response include the viscous damping due to the fluid viscosity and the radiation damping due to the water compressible and free surface waves. Neglecting the fluid viscosity, many studies were performed to study the earthquake-induced hydrodynamic pressures on cylinders [3, 11, 15, 17, 21, 24, 26, 29]. The study conducted by Liaw and Chopra [3] indicated that water compressibility and free surface waves can be neglected in solving the earthquake-induced hydrodynamic pressure on circular cylinders. Then, the pile-water interaction can be conveniently modeled as the product of an added mass of water and the acceleration of the pile. In this study, the added mass is used to represent the lateral resistance of the pile-water interaction. The simplified formula of the added mass for a single pile can be expressed as [21]where , , h is the water depth, r0 is the radius of the pile, is the water density, and z is the distance above seabed.

Recently, Wang et al. [29] proposed an accurate and efficient finite element model to calculate the earthquake-induced hydrodynamic forces on multiple circular cylinders. This numerical model is used to calculate the coefficients of the pile group for the added mass in the present study. The coefficients of the pile group for these 15 piles with different water depth are illustrated in Table 1.

1.5. FEM and Validation

Based on the OpenSees platform, a 3D finite element analysis model of group pile bridge towers was established. The schematic diagram of FEM is shown in Figure 5(a). In the model, the soil-pile interaction is simulated by the nonlinear zero-length elements. The pile-water interaction is considered by added mass shown in Figure 5(a).

For the FEM model comparison, the FEM of bridge group pile foundations is established by SAP2000. The pier and piles are simulated by the linear frame element. Except for the sections of piers, the section definition of the SAP2000 model is completely equal to the OpenSees model. For the pier section, the displacement-based beam-column element with fiber is used in the OpenSees model. In contrast, the pier in SAP2000 is simulated by linear frame elements. Therefore, the cross section parameters of pier section in the SAP2000 model are equivalently converted according to the OpenSees model. The FEM based on SAP2000 is shown in Figure 5(b).

Both vibration periods of the first 3 modes and seismic responses on critical sections of the OpenSees model and SAP2000 model are compared in Table 2. It can be seen that both the vibration periods and seismic responses agree well with each model. Therefore, the validated FEM model based on the OpenSees platform can be used in the following studies.

2. Seismic Fragility Analysis of the System

2.1. Methodology of Component Fragility Approach

The seismic fragility function, which defines the probability of physical damage as a function of ground motion intensity, is widely used to develop vulnerability information for bridges [4952].

The seismic fragility curves are the relation between the ground motion intensity measure (IM) and the seismic demand (D) of structures and modeled as lognormal distribution [53]:where C is structural capacity; SC is median estimate of the capacity; SD and are the median value and dispersion of the structural demand as a function of an IM, respectively; and Φ[⋅] denotes the cumulative standard normal distribution function.

Generally, the relation between SD and IM can be characterized by a powerful log-linear model [53], defined aswhere a and b are the power-law model coefficient parameters. And the dispersion of the demand can be estimated bywhere Di and N are the ith realization of structural seismic demand and the number of conducted nonlinear time history analyses (NHTA), respectively.

Substituting equation (5) into equation (4), the fragility function can be rewritten asin which and represent the median and dispersion of the generated fragility curve.

2.2. Incremental Dynamic Analysis

Incremental dynamic analysis (IDA) is a computer analysis technique that is widely performed to assess the seismic performance of structures subjected to earthquake loads. IDA is used to generate the analytical fragility curves in the present study. For carrying out the IDA, the PGA is utilized as the intensity measure (IM) of each record, and PGA of all ground motion records is scaled from 0.1 g to 1.5 g with the increment of 0.1 g.

The longitudinal curvature ductility (μ) of the pier and pile sections is assumed as engineering demand parameters (EDPs). The curvature ductility limit states are also evaluated using the relative curvature ductility ratio of the pier, which is defined aswhere is the relative pier top curvature obtained from dynamic analysis and is the relative pier curvature when the longitudinal reinforcement reaches its first yield. All EDPs for the highway bridge column are presented in Table 3.

In the structural response analyses, the curvature ductility of the pier or piles is recorded and the maximum response is employed as the damage index (DI) of the corresponding critical sections. The curvature ductility capacity, namely, μC in equation (4), can be obtained based on expert judgments, experimental data, or numerical simulations. In this study, the value of μC is determined by referencing the capacity-based limit state (LS) models defined in previous studies [5557]. The curvature ductility of different sections is graduated into four damage LSs, namely, slight, moderate, extensive, and collapse. Table 3 summarizes the DI values corresponding to the four LSs of each section. In this paper, the index values provided by Choi et al. [54] in the table are used for seismic vulnerability analysis.

2.3. Input Ground Motions

In this study, the far-field suite of 22 records introduced in the FEMA P695 (2009) is considered to investigate the seismic fragility of the FSSI system, which occurred between 1971 and 1999. Among them, 8 events were from California earthquakes and 6 events were from five different foreign countries. For the far-field record set, event magnitudes range from M6.5 to M7.6 with an average magnitude of M7.0. Figure 6 represents the average spectrum of the far-field record set from FEMA and different acceleration response spectra with 5% damping ratio.

2.4. Seismic Fragility of the Pier-Group Pile System

Figure 7 shows the computed demand values of different sections, i.e., Sec C and Sec D, with respect to the ground motion IM (PGA) in the logarithm coordinates. The linear regression results are also illustrated. It can be observed that the relationships between the various structural demands and the PGA value fit well with the corresponding log-linear models. Using the parameters obtained from the regression analyses in equations (2) and (3), the median and dispersion values for the seismic fragilities of various bridge components are calculated.

Table 4 summarizes the corresponding fragility median and dispersion values of different critical sections of the pier and piles. The generated component fragility curves under the four damage LSs are plotted in Figure 8. As shown in Figure 8, Sec C is the most vulnerable section in the example FSSI system for all of the four damage LSs, followed by Sec B and Sec A. The fragilities of Sec E are the lowest. This is owing to the fact that the bottom section of the pier suffers very large moment during earthquake excitations.

The seismic responses of Sec E are much lower since the stiffness and damping provided by surrounding soil are very limited. It can be observed from Figures 8(c) and 8(d) that the seismic fragilities of the critical sections other than Sec C are very low under the extensive and collapse damage, and the corresponding fragility median values are very high (as listed in Table 4), which means that the damage probabilities of these sections are negligible and the seismic performance of the FSSI system under higher damage LSs is controlled by Sec C.

3. Parameter Analysis

3.1. Effect of Pile-Water Interaction

To analyze the effect of pile-water interaction on different pier systems, there are two models including Model A and Model B are discussed in this section. Model A considers the pile-soil interaction, and the detail model is introduced in Section 2. Model B is assuming the soil rigid. In Model B, the piles are fixed on the interface section between soil layer and water.

Table 5 shows the fundamental period (T1) for Model A and Model B, which is 3.51 s and 2.25 s. It indicates that assuming the fixed boundary condition will obviously decrease the fundamental period by 36%. For the different model, the fundamental periods of foundation when considering FSI and without FSI are shown in Table 5. It can be seen that pile-water interaction can increase the fundamental period of the bridge pile foundation and it has more influence for Model A. When water depth is 40m, the fundamental period is increased about 4% for Model A while it is increased about 1% for Model B.

The effects of pile-water interaction on the seismic responses of the bridge pier are firstly investigated, where four earthquake records shown in Table 6 are used. Table 7 shows the maximum displacement on the pier head versus water depth. Table 8 shows the maximum curvature at Sec C versus water depth. It can be seen that the seismic responses of the bridge pier in the case of Model B are increased by the pile-water interaction and this influence increases as water depth increasing. However, for Model A, the seismic responses of the bridge pier may be increased or decreased by the pile-water interaction, which is relevant to the characteristics of the earthquake records, water depth, and dynamic properties of the bridge pier significantly changed by pile-soil interaction and pile-water interaction, as shown in Table 5. It should be noted that pile-water interaction can increase the seismic responses of the bridge pier in most cases. When water depth is 40m, the seismic responses of the bridge pier can increase about 10% in some cases.

Figures 9 and 10 show the fragility curves of Model A and Model B at Sec C for different water depth. It can be seen that pile-water interaction has slight influence on the fragility curves of the bridge pier in the present study and this influence is less than 5%. Tables 9 and 10 show the probability ratio of fragility curves between considering FSI and without FSI. It can be seen that the value of fragility curves in the case of Model B is increased by the pile-water interaction and the probability ratio increases as water depth increases. However, pile-water interaction makes the value of fragility curves in the case of Model A decrease about 2% when water depth is 40m. It also can be seen that the effects of pile-water interaction on fragility curves decrease as earthquake intensity increases. It should be noted that the numerical results conducted by Wang et al. [17] indicated that the FSI has significant effect on the seismic responses of the cylinder when the water depth is much larger than half height of the cylinder, while the effect of FSI on the seismic responses of the cylinder can be neglected when the water depth is less than half height of the cylinder. This study also indicates that the effect of FSI on the fragility curves can be neglected when the water depth is less than half height of the bridge pier.

3.2. Effect of Different Types of Ground Motions

The previous analysis only discussed the seismic vulnerability of the bridge pile foundation system subjected to far-field (FF) earthquakes. In fact, cross-sea piers may also suffer from near-field earthquakes as the widespread existence of submarine faults. Therefore, this section will discuss the seismic vulnerability of the bridge pile foundation system under different types of earthquakes. This paper still selects 28 recommended near-field (NF) records from the US specification FEMA P695 as input motions. All NF motions are recorded in the earthquakes occurred between 1976 and 2002. Among them, 7 records are from United States earthquakes (six in California) and another 7 records are from five different foreign countries. Event magnitudes range from M6.5 to M7.9 with an average magnitude of M7.0. As judged by wavelet analysis classification of the records [59], recommended NF motions are divided into two sets of 14 NF motions: pulse near fault (PNF) and no-pulse near fault. The acceleration response spectra of each motion of PNF set and NNF set (PGA = 1.0 g) are shown in Figures 11(a) and 11(b). The mean acceleration response spectrum comparison is shown in Figure 12. For these three sets of motions, the FF has larger value in short-period range (T < 1 s). However, PNF has larger value than FF and NNF when the period larger than 1 s.

The fragility curves of Sec C and Sec D subject to NF ground motions in four LSs are shown in Figure 13. The damage probabilities of Sec C are always larger than Sec D when subjected to the NF ground motions, which is similar to the results of FF ground motions. It means that Sec C is always the most critical section for the FSSI system both for the FF and NF excitations. On the other hand, the fragility curves suffered to different type motions are compared. It is obvious that the damage probabilities suffer to PNF ground motions are significantly larger than NNF and FF for both Sec C and Sec D. And those of NNF and FF ground motions are closer. Taking excitation of PGA=1.0g as an example, the damage probability comparisons are shown in Figure 14. For Sec C at the slight damage LS, the damage probabilities under PNF are about 127% than the FF motions, while those under NNF are 87% than FF motions. This phenomenon can be explained from the acceleration response spectrum curve in Figure 14. For the fundamental period of the FSSI system (T1= 3.6 s), the response of the PNF ground motion is significantly larger than NNF and FF.

3.3. Effect of Different Soil Conditions

In order to examine the effects of plies at different soil conditions on the seismic fragilities of the bridge pile foundation system, the three different soil conditions including stiff, moderate stiff, and soft are discussed in this section. The value of the undrained shear strength Cu is suggested by [60]. The values of and J are determined according to [47]. The parameters of three soil conditions are shown in Table 11.

The seismic fragility of the bridge pile foundation system under different soil conditions is calculated. The fragility curves of Sec C in four LSs are compared in Figure 15. We can observe that the fragilities in stiff soil are the highest in all LSs. It means that the FSSI in stiff soil is the most vulnerable to earthquake damage than other soil conditions. The fragilities in moderate soil are slightly less than stiff soil, while that in soft soil is obviously less than others. Taking the PGA = 1.0 g as an example, the damage probabilities for stiff are 0.89 in the slight damage LS, which are 1.06 times and 1.19 times comparing to the moderate and soft soil site. For the other three damage LSs, the fragilities of bridge pile foundation at stiff soil are also slightly larger than moderate soil. And both are obviously larger than soft soil. The results seem can be explained from the changing of the system fundamental period. When the piles are surrounded by different soil types, the dynamic characteristics of the bridge piles systems are various, whose fundamental periods are 3.14 s, 3.39 s, and 3.95 s, respectively. It means the bridge pile system in the stiff site condition leads to fundamental period decrease of 8% for moderate soil and 21% for soft soil. In the result, the bridge pile foundation in stiff soil is most vulnerable to earthquake damage than other soil conditions.

4. Conclusions

A three-dimensional finite model is developed to reliably assess the seismic fragilities of the bridge pile foundation system considering FSSI, where effects of hydrodynamic added mass, soil conditions, and different types of ground motions on the fragilities are comprehensively taken into account. Based on the numerical examples, the obtained results can be drawn as follows:(1)Pile-water interaction can increase the fundamental period of the bridge pile foundation, and it has more influence when pile-soil interaction is considered. In general, FSI can increase the seismic responses and fragilities of the bridge pier, and this influence increases as water depth increases and decreases as earthquake intensity increases. The effect of FSI on the fragility curves can be neglected when water depth is less than half height of the bridge pier.(2)Compared to the FF ground motions, the PNF motions yield higher seismic responses and fragilities for the example bridge pile foundation (up to 27%), and the NNF ground motions do not always lead to higher fragilities. In general, the influence of PNF motions should be seriously considered in seismic fragility analysis of the sea-crossing bridge pile foundation considering fluid-pile-soil interaction.(3)The bridge pile foundation in stiff soil is most vulnerable to earthquake damage than other soil conditions. The damage probabilities for stiff case are 1.06 times and 1.19 times comparing to the moderate and soft soil site.

It should be noted that the influence of vertical excitations is not considered in this paper. As the notable influence on the long-span whole bridge, it will be carefully considered in the subsequent studies on the seismic response of the whole offshore bridge.

Data Availability

No data were used to support this study.

Disclosure

The results and conclusions presented are of the authors and do not necessarily reflect the view of the sponsors.

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 the Postdoctoral Science Foundation of China (2018M631292), the Postdoctoral Science Foundation of Beijing (2018-ZZ-032), Natural Science Foundation of Fujian Province (2019J01233), and National Natural Science Foundation of China (51708010). The support is gratefully acknowledged.