Study on Dynamic Response of Inclined Pile Group Foundation under Earthquake Action

: With the development of all types of new special-shaped arch bridges, the deformation of the pile foundation of special-shaped steel arch bridges has a great inﬂuence on the superstructure, and the requirements for it are strict. In addition, the force of the inclined pile group under an earth-quake is complicated. In this paper, based on the Zhuhai Jindao Bridge project, the ﬁnite element numerical model is established in order to study the dynamic response and parameter inﬂuence law of an inclined straight pile group foundation under earthquake action. After a thorough analysis of the results of the simulation, it was found that under the action of an earthquake, the inclination angle of inclined piles has a signiﬁcant inﬂuence on the seismic ability of pile group foundations. When the inclination angle of an inclined pile is 5°, the ability of the inclined pile group foundation to resist deformation can be eﬀectively improved, and the load distribution of the pile top is more reasonable.


Introduction
With the continuous updating of urban infrastructure, special-shaped steel arch bridges with unique shapes have gradually become an important choice for urban architecture.Complex structures put forward new requirements regarding the force of pile foundations.An inclined straight pile group foundation consists of inclined and straight piles, and its working mechanism is different from that of conventional pile foundations.When inclined piles participate in bearing lateral loads in the pile group, the displacement of the pile top and the internal force of the pile body are reduced by varying degrees, making inclined straight piles more suitable for arch bridge foundations [1,2].The dynamic response law of inclined straight pile groups under earthquake action is complex [3].Currently, there is a lack of seismic design provisions for inclined straight pile groups in seismic design specifications.Therefore, it is necessary to conduct research on the dynamic response of inclined straight pile groups.
The bearing capacity and seismic performance of inclined pile groups have been paid a ention to by geotechnical engineering researchers.The research methods of inclined piles mainly include laboratory tests [4][5][6][7][8] and numerical calculations [9][10][11][12][13][14][15].Some scholars also analyzed inclined piles using the p-y curve method.Zhou Dequan et al. [16] set an inclined pile at different angles in the model sand and applied a horizontal load to the top of the pile, so as to obtain the change law of p-y curve characteristics with the inclination angle and pile bo om constraint conditions.Chen Wenlong et al. [17] carried out a vibration table comparison test of inclined and straight piles in saturated sand and dis-cussed the seismic performance of inclined pile foundations from three aspects: superstructure acceleration, pile top shear force and the pile body bending moment.In recent years, many scholars have carried out studies on the seismic performance of inclined pile foundations through numerical simulation [17][18][19][20][21][22].Ref. [23] obtained peak cap acceleration and displacement based on a shaking table test and compared them with FLAC-3D numerical simulation.The results show that FLAC-3D software is able to be er simulate the dynamic response of inclined piles.In this paper, focusing on the stress characteristics of the inclined pile group and relying on a real-life engineering scenario, the finite element software MIDAS Civil was used to establish a three-dimensional model, and the timehistory analysis method was used to analyze the seismic response of inclined pile group foundations.
The special-shaped arch bridge studied in this paper is a separate arch rib structure, and the stress and parameter influence rules during earthquakes are different from those of conventional arch bridges, resulting in a complicated load transfer to the pile foundation.Therefore, based on the stress characteristics of the special-shaped arch bridge studied in this paper, the stress conditions and deformation rules of the inclined straight pile group structure of an arch bridge were analyzed.Relying on a real-life engineering scenario, the finite element software MIDAS Civil was used to establish a three-dimensional model and verify the rationality of the model.On this basis, the seismic response analysis of the inclined straight pile group foundation was carried out using the time-history analysis method.This provides a reference for the design of inclined pile groups for this type of special-shaped arch bridge in practice.
The remainder of this paper is organized as follows.Section 2 presents the details of the engineering scenario.In Section 3, the establishment of the finite element model is explained.The analysis of natural vibration characteristics is shown in Section 4. The computational results are presented and analyzed in Sections 5 and 6, respectively.Finally, Section 7 presents the conclusions.

Project Profile
The Jindao Bridge is located in Jinwan District, Zhuhai.The arch rib is a separate structure.After the completion of the bridge, it will become a beautiful part of the scenery in Zhuhai with its beautiful and unique shape.The bridge spans the main flood drainage channel.The starting point of the bridge is KC0 + 132.00, and the endpoint of the bridge is pile number KC0 + 232.00.An asymmetrically shaped arch bridge structure was adopted.A 3D render of the Golden Island Bridge is shown in Figure 1.

Superstructure Profile
The total width of the bridge is 46.25 m, and the main beam is a steel-concrete composite beam structure.Three box longitudinal beams are set on the main beam.The top and bo om of the side longitudinal beams are the same as the bridge floor, and the transverse slope is 1.5% inward.The width of the side longitudinal beams is 2.0 m, and the height is 1.5 m.The longitudinal beam is provided with a transverse partition every 4 m, and a steel beam is located between the longitudinal beam of the side box and the longitudinal beam of the middle box.The beam spacing is 4.0 m, and the beam is part of a "working" section.The thickness of the upper and lower flanges is 24 mm, the thickness of the web is 12 mm, the height of the end beam is 1.5 m, the thickness of the upper and lower flanges is 30 mm, and the thickness of the web is 20 mm.A small longitudinal beam is arranged in the middle of the span between the beams.The small longitudinal beam consists of an I-section, the beam height is 0.7 m, the thickness of the upper and lower flanges is 16 mm, and the thickness of the web is 12 mm.
In addition, the arch rib is connected with the main beam via a derrick.The derrick comprises 19 × Φs15.24; 27 × Φs15.24; and 37 × Φs15.24 types, whose tensile strength is no less than 1860 MPa.The distance between the derrick and the bridge is 4.0 m.The wind support utilizes a Φ630 × 16 mm seamless steel pipe.
The theoretical span of the No.1 to No.3 arch ribs is 128.125 m, 110.0 m, and 128.125 m; the vector height is 44 m, 40 m, and 44 m; and the vector span ratio is l/2.91,1/2.75, and 1/2.91,respectively.The section outline dimensions of arch ribs are 2.8 m × 2.8 m, and the longitudinal section is equal to this section.The main bridge is a steel arch bridge with a span arrangement of 1 × 100 m.

Pile Foundation Layout Scheme
The arch seat and cap are made of C40 concrete, the cap is rectangular, the size of the #1-3 arch seat cap is 12.4 × 16.1 m, and the thickness is 4.0 m.Nine boreholes and cast-in piles with diameters of 2.2 m and 3 inclined steel pipe piles with diameters of 1.5 m are used.The size of the #4 arch seat cap is 16.8 × 35.9 m, and the thickness is 5.0 m.Twentyfour drilled cast-in piles with a 2.2 m diameter, and fourteen inclined steel pipe piles with a 1.5 m diameter are used.The #4 arch seat cap is close to the abutment, part of the platform is directly set on the #4 arch seat cap, and a 2 cm fault joint is set between the #4 arch seat and the abutment.The structure diagrams of #1~#4 arch are shown in Figures 2 and 3

Establishment of a Finite Element Model
In bridge engineering, when the span of the structure is much larger than the other two dimensions (width and height of the section), it can be approximately simplified to the bar structure in model analysis.A large number of experiments and calculations show that this simplification meets the accuracy requirements of bridge structural analysis [24,25].In this paper, in order to meet the requirements and accuracy of model analysis and to avoid the memory and time required for calculation, the model is properly simplified, which not only reduces the corresponding calculation amount but also speeds up the research progress.
In the modeling of the whole bridge structure, X and Y were defined as the transverse and longitudinal bridge directions, respectively, and Z is defined as the gravity direction.The structure was transformed into a space rod model, which divided the bare arch state and the bridge stage.The suspender was simulated using a tensile element only.The arch rib, wind support, main beam, arch seat, and pile foundation were simulated using beam element, and the bridge panel was simulated using the plate element.The whole bridge was divided into 2250 nodes and 2830 elements.The numerical model of the bridge is shown in Figure 4.In addition, the relevant material properties used in this model are shown in Table 1.

Main Beam Simulation
The whole main beam is composed of three main longitudinal beams, spacer beams, and small and medium-sized longitudinal beams across the beam, all of which are simulated via beam elements.Considering the 1.5% transverse slope of the longitudinal beam, the virtual beam is set on both sides, a rigid connection is adopted between the longitudinal beam and the main beam, and a total of six supports are set at both ends of the middle longitudinal beam and the main longitudinal beam.The constraints of the supports are presented in Table 2.All sections were drawn using CAD and imported into Midas/Civil SPC to calculate the section characteristics, to generate the section of each beam body, and then to assign the section and material to the element.The main beam model is shown in Figure 5.
Here, X is the transverse bridge direction, Y is the longitudinal bridge direction, and Z is the vertical bridge direction.In addition, 1 means that the degree of freedom is relaxed, and 0 means that the degree of freedom is constrained.

Arch Rib and Foundation Simulation
The cross-sections of three types of arch rib segments were drawn in CAD, imported to calculate the section characteristics, and assigned to the section of arch rib.After that, the plane coordinates of the three arch ribs with the steel-concrete junction as the origin were calculated according to the quadratic oblique parabola formula in the design drawing.Since curves cannot be recognized in Midas Civil, multiple distributed arch rib nodes were connected insofar as possible through multiple line segments in CAD to ensure a be er arch rib alignment in the modeling.In addition, corresponding nodes were set at the connection between the boom and the main arch.The arch rib and wind-bracing models are shown in Figure 6.
Arch #4 is located at the intersection of the three arch ribs, and the force is complicated.In order to make the load distribution of pile foundation more reasonable, a bearing cap is set at the large arch seat of the model as a transmission medium, and a bearing cap is set at the arch seat #4.The bearing cap is simulated using the beam element, so that the load of the superstructure can be reasonably distributed to pile foundation.A rigid connection is set at the intersection of the three arch feet projected in the Z direction and the simulated cap of the beam unit.The cap and each single row of piles are rigidly connected with the cap as the main node and the top of the pile as the secondary node, so that the load transmi ed by the cap is relatively evenly distributed to the pile foundation.For the end of the small arch seat, only the single arch transmits the load.Therefore, the pile top displacement can be reasonably distributed through the direct use of the primary-secondary node connection.The arch foot and pile group at the #1-#3 arch seat are rigid and connected through the "arch main pile slave".In addition, in order to consider the influence of the interaction between piles and soil, the soil spring method is used to simulate the effect of soil on the pile, and the stiffness of the equivalent soil spring is calculated through the "m" method.The pile foundation nodes are simulated using the linear node elastic support.The boundary constraints are shown in Table 3.

Load Condition
The special-shaped arch bridge is located in Zhuhai City.Besides the general design load, load considerations should be combined with reality.The reasonable consideration of the load condition is more conducive to the fi ing of the finite element model to reality.The load types include structural dead weight, secondary dead load, crowd load, automobile load, temperature load, the initial tension of the boom, and wind load.In the finite element model, the self-weight of the structure will be automatically calculated.The load grade of vehicles is decided at city level.According to the General Code for the Design of Highway Bridges and Culverts (JTG D60-2004) [26], the crowd load q = 2.5 kN/m 2 , the overall temperature of the whole bridge is calculated according to a temperature rise of 21° and a temperature drop of 29°, and the wind load is considered according to the 100year-old wind standard of the service stage.The initial tension of the boom is shown in Table 4.

Analysis of Natural Vibration Characteristics
The premise of the study of dynamic problems is the calculation of the dynamic characteristics of the structure, including the natural vibration frequency and vibration type of the structure and so on.Each natural vibration mode has a fixed frequency, damping ratio, and shape, and its calculation results are the foundation of seismic analysis.The purpose of modal analysis is to avoid certain external excitations (such as earthquakes, crowds, and wind vibrations) approaching the natural frequency of the structure during structural design, thereby preventing resonance and vibration phenomena in the structure.This allows for a more comprehensive analysis of how the structure responds to different types of dynamic loads, which is beneficial for solving and calculating control parameters in other dynamic analysis.
The principle of modal analysis is based on the undamped free vibration equation.The structure is discretized into a calculation model of finite degrees of freedom, ignoring the influence of structural damping [27].In this case, the damping matrix and the node load matrix are equal to 0, and the undamped free vibration equation of the structure is as follows: where [ ] is the structure mass matrix, [ ] is the structural stiffness matrix, ̈ is the structural acceleration vector, and { } is the structural displacement vector.Hypothetical order: Then, there is: where ω is the vibration frequency, is the time, φ is the phase angle, and Ø is the displacement amplitude vector that does not change with time.Hypothetical order: The Formulas (2), (3), and (4) can be substituted into Formula (1), so that the final arrangement can be obtained: The main purpose of modal analysis is to calculate the eigenvalues and corresponding eigenvectors.The above equation is the characteristic equation of structural vibration.The commonly used methods to calculate eigenvalues in Midas Civil include the eigenvalue vector method and the multiple Ri vector method.In this section, the Lanczos method in the eigenvalue vector method is used for calculation.
According to the Code for Seismic Design of Buildings [28], in order to accurately analyze the seismic action, the modal participation mass of the structure should reach at least 90%.This paper calculates the first 200 vibration modes of the structure, and the cumulative modal participation mass of X, Y, and Z directions reaches 99.5%, 99.6%, and 96.2%, respectively, meeting the requirements of the code.Due to the limited space, the first 10 order frequencies and modes of the structure are selected, and the results are shown in Table 5.

Seismic Response Time History Analysis
The average intensity of an earthquake at the bridge site is 7 degrees, the peak acceleration of ground motion is 0.1 g, and the category of the engineering site is class III.This paper selects a Taft wave suitable for this type of site as the ground motion input.On the basis of the selected seismic wave, the seismic wave is adjusted according to the design of working conditions, and the basic principle of acceleration value adjustment is as follows: the peak value of the actual seismic wave is the same as that of the time-history analysis method (see Formula ( 6)).In this paper, a total of 7 degrees is applied (the peak value of the design acceleration is 0.35 m/s 2 ).The revised time-history curve is shown in Figure 7.
where ́( ) is the adjusted acceleration time-history curve, ( ) is the actual seismic wave before adjustment, is the peak acceleration of the adjusted acceleration curve for time-history analysis, and is the peak acceleration of the acceleration time-history curve used for time-history analysis before adjustment.According to the selected and adjusted seismic waves, the three-way combined input can be used to obtain the top displacement of each single pile in the pile group foundation, and the results are shown in Figure 8.In addition, the maximum displacement of each column of piles is shown in Figure 9. Figures 8a and 9a show the #1, #2, and #3 pile foundations from left to right, respectively.Under the action of earthquake, the displacement of the pile top of the column near the negative direction of Y of each pile group from #1 to #3 changes greatly, and the displacement of pile top reaches the maximum value.The maximum displacements of single pile tops in the #1, #2, and #3 pile groups are 2.556 mm, 2.144 mm, and 2.478 mm, respectively.In addition, the displacement trend of the corresponding single pile along the pile body is roughly the same among the pile groups.At the same time, due to the small depth of the inclined pile embedded in the sandstone, there is still a small amount of horizontal displacement near the bo om of the inclined pile, but the overall displacement trend is still similar to that of the straight pile group.
It can be seen from Figures 8b and 9b that the horizontal displacement of piles on both sides of pile group #4 is slightly larger than the horizontal displacement of piles in the middle two rows.Due to the greater load at #4, the maximum displacement of pile body is greater than from #1 to #3, which is 6.498 mm.The forces of the piles in the back row decrease in turn, but the displacement difference is small, and the maximum difference is 0.043 mm, which can be considered very close.There is still a small amount of horizontal displacement in the bo om part of the inclined pile, but the overall displacement trend is still similar to the straight pile group, and the principle is the same as the analysis of #1~#3 pile group.

Influence of Pile Tilt Angle on Seismic Response
In this paper, the angle between the direction of inclined pile and the direction of gravity is 7°.In order to study the influence of pile inclination on the seismic response of inclined pile group, this paper keeps other parameters unchanged and selects four angles, As can be seen from Figures 10 to 13, the overall displacement of #1~#3 pile group foundation is stable with the change in pile body as the angle of inclined pile gradually increases, but the displacement of pile top gradually decreases with the increase in angle, indicating that the increase in inclined pile angle improves the seismic performance of the #1~#3 pile foundation.In the process of increasing the angle of pile group #4 (assuming that each column along the positive direction of Y successively represents the first, second, third, and fourth columns, the top displacement of piles in the first and second columns gradually increases, while that of the third and fourth columns gradually decreases.When the tilt angle is 5-7°, the top displacement of each single pile in pile group is close, indicating that the force of pile group #4 is more uniform.Furthermore, as can be seen from Figure 14, when the inclination is 5°, the displacement of the pile top is small, which is conducive to the resistance of pile foundation #4 to the deformation caused by earthquakes.

Conclusions
In this paper, the seismic response of inclined straight pile group foundation of this type of special arch bridge was studied and analyzed, the working characteristics of in-clined straight pile group during earthquakes were examined, and the influence of different inclined pile angles on the ability of pile group foundation to resist the deformation caused by earthquakes were discussed.The main conclusions are as follows: 1.Under the action of earthquakes, the deformation resistance of the #1-#3 pile group foundation is gradually enhanced with the increase in inclined pile angle, which improves the seismic performance of the pile foundation to a certain extent.The deformation resistance of the #4 pile group foundation shows a trend of decreasing first and then increasing with the increase in inclined pile angle, and the horizontal deformation resistance is be er when the inclination is 5°. 2. Under the action of earthquakes, when the inclination angle is 5-7°, the pile top load distribution of the #4 pile group foundation is more uniform, which is conducive to the interactions of each single pile.Therefore, the inclination angle of the inclined pile group foundation of the bridge type should be set within a suitable range.

Figure 1 .
Figure 1.A 3D render of the Golden Island Bridge. .

Figure 4 .
Figure 4. Full bridge model of special-shaped steel arch bridge.

Table 1 .
Material properties used in the numerical analysis.

Figure 5 .
Figure 5. Full bridge model of special-shaped steel arch bridge.

Table 2 .
Boundary conditions of the supports.

Table 4 .
Initial pull of the boom.

Table 5 .
First 10 natural vibration frequencies of the structure.