Ying Wu · Xuecheng Bian · Chong Cheng ·Jianqun Jiang
Abstract An efficient computational approach based on substructure methodology is proposed to analyze the viaduct-pile foundation-soil dynamic interaction under train loads. The train-viaduct subsystem is solved using the dynamic stiffness integration method, and its accuracy is verified by the existing analytical solution for a moving vehicle on a simply supported beam. For the pile foundation-soil subsystem, the geometric and material properties of piles and soils are assumed to be invariable along the azimuth direction. By introducing the equivalent stiffness of grouped piles, the governing equations of pile foundationsoil interaction are simplified based on Fourier decomposition method, so the three-dimensional problem is decomposed into several two-dimensional axisymmetric finite element models. The pile foundation-soil interaction model is verified by field measurements due to shaker loading at pile foundation top. In addition, these two substructures are coupled with the displacement compatibility condition at interface of pier bottom and pile foundation top. Finally,the proposed train-viaduct-pile foundation-soil interaction model was validated by field tests. The results show that the proposed model can predict vibrations of pile foundation and soil accurately, thereby providing a basis for the prediction of pile-soil foundation settlement. The frequency spectra of the vibration in Beijing-Tianjin high-speed railway demonstrated that the main frequencies of the pier top and ground surface are below 100 and 30 Hz, respectively.
Keywords Train moving load · Viaduct · Pile foundation · Soil · Vibration · Substructure method
Pile-supported viaducts are widely used in high-speed railways across soft soil deposited areas or densely populated urban areas. In China, viaducts are used in more than 70% of newly built high-speed railways. When trains pass through urban areas, the ground vibrations induced by trains running over viaducts transmit to the surrounding buildings due to the interaction between the soil and the structure, leading to serious noise interference and surrounding ground settlement, which brings detrimental effects to people’s work and lives. There are many studies on the environmental noise problem; for example, in Refs. [1-3], the researchers proposed some theoretical models to predict the environmental noise generated by trains passing through viaducts and discussed noise reduction measures. However, less attention was paid to the prediction and evaluation of settlement caused by ground vibrations due to trains running on viaducts. To this end, it is necessary to build a train-viaduct-pile foundation-soil interaction model to predict effectively the vibrations of pile foundation and soil.
To date, much work has been done to study train-bridge(viaduct) dynamic interaction for the purposes of safe train operation and railway structure design. Frýba [4] derived an analytical solution for the beam structure subjected to a moving load. Yang et al. [5] simplified train as a series of moving load to investigate parameters that control the dynamic bridge responses, such as bridge’s span length,car length, and train speed. Zeng et al. [6] proposed a train-track-bridge coupled dynamics modeling method which greatly reduces the matrix dimension and computer performance requirements. Gou et al. [7] studied the dynamic responses of a long-span cable-stayed railway bridge with steel-concrete hybrid girders in different train load scenarios through in situ testing. Cheng et al. [8-10]proposed a bridge-track-vehicle element for analyzing vibrations of railway bridges under a moving train, which shows that bridge vibration has a significant effect on the dynamic response of railway track.
The foregoing studies on train-bridge (viaduct) dynamics provide theoretical basics for understanding train-viaduct dynamic interaction and vibration generation. So far, traininduced ground vibrations have been widely studied, but most of the published works focus on vibrations induced by trains running on ground surface [11-14], and very limited literature is available on ground vibrations induced by train running on viaduct [15-18]. Zhang et al. [15] proposed a semi-analytical model to investigate structure and groundborne vibrations caused by trains running over viaducts,taken into consideration the vibration-isolating effects of the multi-layered elastic supports. Shi et al. [16] presented a theoretical evaluation on vibrations of saturated ground induced by high-speed train operating on an elevated bridge and identified two characteristic frequencies of the superstructure that dominates the ground displacement response.In addition, Takemiya and Bian [17] proposed a substructure method to study the ground vibration induced by Shinkansen train and validated it with field measurements. Fujikake[18] proposed a vehicle-bridge-ground dynamic interaction model to investigate the environment vibration along Shinkansen railway and analyzed vibration characteristics of the bridge structure and alongside soil.
In the literature, most analyses on bridge-pile-soil interaction under the excitation of trains fall into two categories:(1) simplified analytical solution and (2) fully three-dimensional numerical model. The simplified analytical solution can get the vibration results of a bridge and its nearby ground at a very low computational cost and get insight into structure vibration and propagation in soil. But this kind of models usually imposed some assumptions on bridge structure, pile foundation, and soil, so they can be described by simple analytical formulations. The fully three-dimensional finite models can well describe the details of bridge, pile foundation, and ground soil and incorporate more realistic material attributes, but computations of these numerical models are very time-consuming, which prevents its application in engineering.
In this paper, the coupled viaduct-pile foundation-soil dynamic interaction system is divided into two subsystems based on the dynamic substructure method. The train-viaduct subsystem is solved by dynamic stiffness integration method, and the pile foundation-soil subsystem is solved using Fourier series decomposition in the azimuth direction. The two subsystems are coupled with the displacement compatibility condition at the interface of viaduct pier bottom and pile foundation top. After the models of the two subsystems are verified, they are integrated to analyze the coupled viaduct-pile foundation-soil interaction under train moving load. Finally, field test validation was used to demonstrate the applicability of the proposed analysis model in predicting high-speed train induced structure and ground vibrations.
The finite element method (FEM) can deal with the complex geometry of viaduct and pile foundation effectively, but the full three-dimensional finite element model of viaduct-pile foundation-soil with dynamic excitation of train loadings can be very computationally expensive since the ground is infinite in nature. The whole viaduct-pile foundation-soil system with train moving load is shown in Fig. 1. It can be divided into two subsystems based on dynamic substructure method to reduce the computation complexity. One is the train-viaduct subsystem, and the other is the pile foundation-soil subsystem. By assuming that the interface of viaduct pier and pile foundation is in close contact, these two subsystems are coupled with the displacement continuity condition at interface. All of these solutions take place in the frequency domain, since some important parameters controlling structure and soil vibrations are more accurately and easily identified in the frequency domain, such as the impedance functions and damping. The time histories of viaduct and ground vibrations are then obtained by the inverse Fourier transform from the counterpart solution.
Several steps are required for solving the whole dynamic system based on the substructure method. The methodology flow chart is shown in Fig. 2. First, the impedance function of pile foundation is derived from the pile foundation-soil interaction, and the transfer function of vibration transmission from pile foundation to surrounding soil is determined. The impedance function of pile foundation is integrated into viaduct governing equations at all pier bottoms. After solving the governing equation of viaduct under train moving load, the loadings acting on pile foundation top can be retrieved and are used in computation of ground vibration based on the transfer function derived in the first step.
Fig. 1 Train-viaduct-pile foundation-soil interaction model
Fig. 2 Flow chart of the methodology
Fig. 3 Train-beam coupling model
In this study, the train car is modeled as a multi-rigid body and its motions are represented by 4 degrees of freedom (DOFs) in the vertical plane, as shown in Fig. 3.Although such simplification can affect the response of ground vibration in high frequency domain, it has negligible influence on the subsequent calculation of pile foundation-soil settlement. The dynamic wheel-rail force and fastener force due to track irregularity mostly affect track dynamic responses in high-frequency band and decays faster, so it has little influence on the pile settlement and the dynamic response of far field. In order to develop a simplified model to predict the pile foundation-soil settlement, the track irregularity is not considered here. The motions of the train car include vertical displacements of the car body and two wheels,and the rotation of the car body. The suspensions of the car body are simplified as primary springs with constantsk1andk2and dashpots with damping coefficientsc1andc2, respectively. The masses of car body and wheels are represented bymcandmwand their vertical movements byycandyw, respectively. The torsion angle and wheelbase are denoted byφcandlc, respectively. The beam stiffness and moment inertia of the car body are denoted bykbandIc, respectively.
The motion of train car can be described by the following governing equation:
whereKv,Cv, andMvare stiffness matrix, damping matrix and mass matrix of car body;yvis displacement vector; andfvis force vector. Their detailed expressions are given as below:
in whichKv,Cv, andMvare symmetric matrices.
In this study, the train load on the viaduct is obtained by solving Eq. (1); then, it is substituted into the viaduct’s governing equation (i.e., Eq. (2)) as the equivalent node load in the frequency domain, for which one can refer to the author’s previous work [17]. The vertical displacement at wheel-viaduct contact position is approximated by linear interpolation of the neighboring two nodes on viaduct.
The model of a moving vehicle on a simply supported beam is used to verify the proposed train-viaduct dynamic interaction analysis model, as shown in Fig. 3, since there is an analytical solution available in Ref. [4]. Parameters of the train car and beam are given in Table 1. As shown in Figs. 4 and 5, the developed model shows good agreement with Frýba’s analytical solution [4], both for the vibrations of beam center and car body center.
The viaduct deck and piers are modeled by three-dimensional beam elements. The track, including rails, fastener,and underlying slab layer, is condensed into the stiffness ofviaduct deck. The supporting of pile foundation to viaduct is represented by frequency-dependent impedance functions.
Table 1 Parameters of train car and beam
Fig. 4 Vibration at beam center
The governing equation of the viaduct with multiple pier supports is formulated in frequency as follows:
whereM,C, andKare the mass, damping, and stiffness matrices, respectively;Utis the displacement matrix;Ftis the force vector;ωdenotes angular frequency and i the imaginary unit.
By separating the displacement and load vectors at superstructure nodes (denoted by subscript ‘s’) and at interface nodes at pier bottoms (denoted by subscript ‘i’), Eq. (2) is rewritten as
Fig. 5 Vibration at car body center
where the matrixes with subscripts ‘‘ss’’ and ‘‘ii’’ denote the self-matrices of the train-viaduct subsystem and pile foundation-soil subsystem, respectively, and the ones with subscripts ‘‘si’’ and ‘‘is’’ denote the interaction matrices between subsystems.
Here,Ftis only the train load on the viaduct, and it is not a stationary process. The train loadFtcan be obtained by solving Eq. (1), and then, it is substituted into the viaduct’s governing equation (i.e., Eq. (2)) as equivalent node load in the frequency domain. Thus, the global dynamic equations of viaduct structure are expressed as
whereKfis the impedance function of pile foundation.
The governing equation of train car movements is integrated into Eq. (4) according to the train-beam interaction model proposed in Sect. 3.1.
Pile foundation is used to support viaduct widely, but the pile foundation-soil interaction has been paid less attention to in previous works. However, soil-structure interaction significantly affects the dynamic response of viaduct under high-speed train loadings [17]. Besides, study on ground vibrations caused by train running over viaduct is another important issue from environmental amenity considerations. Therefore, it is indispensable for soil-pile interaction analysis to determine the vibration propagation from pile foundation to adjacent ground soil. In this regard,semi-analytical solution with axisymmetric FEM is effective to study the dynamic problem of pile foundation with the ring-like arrangement of pile group [19, 20]. Here, we extend this method to study pile-foundation dynamic interaction under loading transferred from viaduct due to train moving loads, as shown in Fig. 6,whereFx,Fz, andMydenote the horizontal force, vertical force,and bending moment, respectively. One of the basic assumptions is that the geometric and material properties are invariable in the azimuth direction. Then, for the axisymmetric bodies subject to non-axisymmetric loadings, the Fourier expansion is applied for both the displacement and load components. This solution has advantages in saving the memory and computational efforts,compared with the full three-dimensional models.
Fig. 6 Pile foundation-soil interaction model
In this paper, the pile foundations and soil are modeled by three-dimensional axisymmetric finite elements with Fourier expansion about the azimuth. A transmitting boundary based on the thin layer element method [20] is introduced into the analysis to model the infinitely extending layered soil.
For the dynamic problem of pile foundation-soil interaction under train loading on viaduct, displacements and stresses of the pile foundation and surrounding soil do not vary dramatically around the circumference. In this case, the Fourier expansion of the displacement along the circumference is performed and motions of the pile foundation and surrounding soil are approximated expressed by the first several series items.
The ground motions are governed by the Navier-Stokes equations as follows:
whereUandPare the soil displacement vector and force vector, respectively;ρis the soil density;λcandμcare the Lame constants of the soil including a small part of complex number to consider the ground damping effect and can be expressed as
whereβis the internal damping of soil and the superscript“c” denotes complex number.
14. Slipped from the palace: The princess in this story is by no means passive such as many of her better known fairy tale counterparts. She attempts to protect herself as much as she can from the threat of immorality79. When she fails, she flees the palace and the threat, abandoning her birthright and the luxury of her life.
The displacement vectorUand force vectorPcan be expressed in cylindrical coordinate system (u,r,θ) using Fourier expansion, as follows:
The superscripts ‘s’ and ‘a’ stand for the symmetric and anti-symmetric terms, respectively. In case of the smoothly deformed configuration, then= 0 and 1 terms would be enough to describe the pile and soil responses. Namely, then= 0 symmetric harmonic represents a vertical motion in theZdirection as well as a symmetric dilation, while then= 0 anti-symmetric harmonic does a torsional motion about theZaxis. Then= 1 symmetric harmonic relates to a coupled motion of translation along theXaxis (θ=0o) and rotation about theYaxis (θ=90o), and then= 1 anti-symmetric harmonic does a coupled motion of translation along theYaxis and rotation about theXaxis. Since the symmetric mode is independent to the anti-symmetric mode, these two motions can be decoupled in the governing equations.bnis body force;tnis surface traction;sis arc length;B1n,B2n,D1,andD2are matrices defined as
A pile foundation mostly comprises a group of piles with a footing on top of them. Therefore, the equilibrium condition between the discrete piles and the continuum soils in this axisymmetric finite element model should be considered carefully. The behavior of piles on the same radius from the center is also approximated by the limited Fourier harmonics ofn= 0 and 1 along the circumferential direction connecting the relevant piles axis. This assumption is reasonable so long as the rings of piles maintain the original axisymmetric configuration during the motion at any depth. In the analysis of the interaction between discrete ring piles and the surrounding continuum soils, the pile element and the neighboring soil element are assumed to have the same translation displacements at the connection nodes. The piles are modeled by the twonode three-dimensional beam elements, and their stiffness is integrated into the soil stiffness.
Field tests of pile foundations subjected to vertical and horizontal dynamic loadings were conducted in Taiwan high-speed railway to verify the proposed analysis method for pile foundation-soil dynamic interaction. The profile of the pile foundation and its finite element model are shown in Figs. 7 and 8, respectively. The piles are modeled into two rings of beam element, and parameters of the model in Fig. 8 are taken from the literature [21]. Because of the high rigidity of the footing structure, it is regarded as a rigid block in the analysis. The motions of the footing are transferred to the soil nodes and pile nodes at the connecting interface with a master-slave relation. In the numerical computation, the amplitude of the dynamic loading is kept as 10 kN to accord with the output force from the shaker in field tests.
Fig. 7 Pile foundation of Taiwan high-speed railway viaduct (in mm): a front view of pile foundation; b top view of section A-A; c top view of section B-B
Fig. 8 Axisymmetric finite element model of pile foundation and site condition
Fig. 9 Comparison of simulated and recorded ground surface velocity responses due to the horizontal loading (prefix ‘F’ stands for computation results and ‘T’ for field test data): a horizontal velocity responses; b vertical velocity responses
Fig. 10 Comparison of simulated and recorded ground surface velocity responses due to the vertical loading (prefix ‘F’ stands for computation results and ‘T’ for field test data): a horizontal velocity responses; b vertical velocity responses
The field test results [21] by the National Earthquake Engineering Center of Taiwan for the pile foundation subjected to shaker loading in Fig. 7 are used to validate the three-dimensional axisymmetric finite element model. In the shaker tests, the output force is proportional to the eccentric mass and the square of the test frequency, such that the frequency of input force varies from 2 to 11 Hz, to guarantee that the output force is within the design capacity of the exciter. Either vertical or horizontal loadings was applied at the top of the pile foundation. Since the harmonic exerting force by shakers generally depends on the driving frequency,the measured data were normalized with respect to a fixed maximum excitation force of 10 kN for both the horizontal and vertical directions. The velocity response of ground surface was recorded at the positions of 0, 10, 20, 30, and 50 m from the center of pile foundation. Test results are presented in Figs. 9 and 10 for horizontal and vertical shaker loadings, respectively. The horizontal and vertical velocity of the computed results under horizontal loads and vertical loads is also presented in the same figures for comparison purposes.Generally, the simulation results show reasonable agreement with the field measurements, which proves the reliability of the proposed computation procedure to evaluate ground responses under pile foundation excitations.
The Beijing-Tianjin intercity express railway is one of the newly built high-speed railways in China with an operation speed up to 350 km/h. About 87% of this railway uses viaducts with span lengths of 32 m, which are supported by pile foundation with pile length up to 50 m. Field tests of ground vibrations due to trains running on a viaduct wereperformed. The parameters of pile foundation and site conditions as shown in Fig. 11 are taken from the literature [22,23]. Observation points at pier top and ground surface in a line perpendicular to the viaduct were selected to record vibrations induced by trains running at different speeds. Distance of these observation points from the pier foot was 10,20, 30, and 50 m, respectively.
Fig. 11 Test site in Beijing-Tianjin intercity express railway
Table 2 Train parameters
Based on the proposed viaduct-pile foundation-soil dynamic analysis model with train moving loads, the computed vibration with time at pier top and on ground surface 20 m from pier foot is plotted in Figs. 12 and 13 for train speeds of 150 and 300 km/h, respectively; field test data at the same locations are plotted in the same figures for comparison. Thecomputed vibrations at these two positions show reasonable agreement with the test data despite some discrepancies. The peak responses due to train wheel axle movement can be easily identified in these results. Both computation results and field test data show that train with a higher speed causes stronger vibrations at the viaduct and nearby ground.
Table 3 Viaduct parameters
The frequency spectra of the vertical velocities of the pier top and ground surface are shown in Figs. 14 and 15 for train speeds of 150 and 300 km/h, respectively. The computed vibrations at these two positions show reasonable agreement with the test data despite some discrepancy. The difference in the range of 40-80 Hz arises from site conditions. The site conditions in FEM are ideal, while the actual site conditions are complex due to the spatial nonuniformity of soil and the variability of soil damping in each layer. The soil stiffness and damping are variable because the soil shear strain is significantly related to high frequency, which affects the response of ground vibration at high frequency. Both computation results and field test data show that the main frequencies of the pier top are below 100 Hz while the main frequencies of the ground surface are below 30 Hz. It means that the higher-frequency characteristics decay faster than lowers. Figure 16 shows that train speed has noticeable effect on the vibration at pier top and nearby soil, but this effect will diminish with the distance from pier at ground surface.
Fig. 12 Vertical responses of the sensors at train speed of 150 km/h: a pier top; b ground surface 20 m from pier foot
Fig. 13 Vertical responses of the sensors at train speed of 300 km/h: a pier top; b ground surface 20 m from pier foot
Fig. 14 Frequency spectra of the vibration at train speed of 150 km/h: a pier top; b ground surface 20 m from pier foot
Fig. 15 Frequency spectra of the vibration at train speed of 300 km/h: a pier top; b ground surface 20 m from pier foot
Fig. 16 Attenuation of the vertical vibration at pier top and ground surface with the distance from pier foot
The Shanghai-Hangzhou high-speed railway is also one of the newly built high-speed railways in China with an operation speed up to 350 km/h. About 90% of this railway uses viaducts which are supported by pile foundation.Field tests of bridge vibrations due to trains running on a viaduct with a bridge span of 24 m were performed.Observation points at bridge midspan and endspan as shown in Fig. 17 were selected to record the vibrations induced by trains running at different speeds. The points at pier top and on ground surface 10 m from the pier foot were also observed.
Fig. 17 Test site in Shanghai-Hangzhou high-speed railway
Based on the proposed viaduct-pile foundation-soil dynamic analysis model with train moving loads, the computed vibration with time at midspan and endspan is plotted in Figs. 18 and 19 for train speeds of 80 and 220 km/h,respectively; field test data at the same locations are also plotted in the same figures for comparison. The computed vibrations at these two positions show reasonable agreement with the test data despite some discrepancy. The peak responses due to train wheel axle movement can be easily identified in these results.
The frequency spectra of vertical velocities at bridge midspan and endspan for train speeds of 80 and 220 km/h are shown in Figs. 20 and 21, respectively. The computed vibrations at these two positions show reasonable agreement with the test data despite some discrepancy.
Vibration reaction at bridge structure due to train running on viaduct is presented in Fig. 22 for two train speeds of 80 and 220 km/h. Both computation results and field test data show a consistent trend of vibration attenuation.
Fig. 18 Vertical response of the sensors at train speed of 80 km/h: a midspan; b endspan
Fig. 19 Vertical response of the sensors at train speed of 220 km/h: a midspan; b endspan
Fig. 20 Frequency spectra of the vibration at train speed of 80 km/h: a midspan; b endspan
Fig. 21 Frequency spectra of the vibration at train speed of 220 km/h: a midspan; b endspan
Fig. 22 Attenuation of the vertical vibration at bridge and ground surface
An efficient train-viaduct-pile foundation-soil interaction analysis model has been proposed herein to study the ground vibration induced by train running on viaduct. Dynamic substructure method is applied to divide the whole system into two subsystems. The dynamic interaction of the train-viaduct subsystem and the pile foundation-soil subsystem is solved by dynamic stiffness integration method and Fourier series decomposition method, respectively. The three-dimensional problem of pile-soil interaction is decomposed into several twodimensional axisymmetric finite element models according to the axisymmetric or anti-axisymmetric modes of loads. The two subsystems are coupled with the displacement compatibility condition at interface of viaduct pier bottom and pile foundation via the impedance function at pile foundation top. The proposed analysis model has been applied to evaluate the ground vibrations induced by train moving load on a viaduct in Beijing-Tianjin intercity express railway, and the computed results are compared with field measurements both at pier top and on ground surface. The model has also been applied to evaluate the viaduct vibrations induced by train moving load in Shanghai-Hangzhou high-speed railway, and the computed results are compared with field measurements both at midspan and at endspan. The results of the case studies have demonstrated the applicability and accuracy of the proposed analysis model in viaduct-ground interaction analysis with train moving loads.
AcknowledgementsThis work was supported by the National Natural Science Foundation of China (Nos. 52125803, 51988101 and 52008369).
Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.
扩展阅读文章
推荐阅读文章
推荐内容
老骥秘书网 https://www.round-online.com
Copyright © 2002-2018 . 老骥秘书网 版权所有