In order to investigate the applicability of dynamic stiffness for bridge cap-pile system, a laboratory test was performed. A numerical model was also built for this type of system. The impact load was applied on the cap top and the dynamic stiffness was analysed. Then, the effect of the effective friction area between pile and soil was also considered. Finally, the dynamic stiffness relationship between the single pile and the cap-pile system was also compared. The results show that the dynamic stiffness is a sensitive index and can well reflect the static characteristics of the pile at the elastic stage. There is a significant positive correlation between the vertical dynamic stiffness index and bearing capacity of the cap-pile system in the similar formation environment. For the cap-pile system with four piles, the dynamic stiffness is about four times as large as the single pile between 10 and 20 Hz.
Fundamental Research Funds for the Central Universities2016JBM034Beijing Jiaotong University1601300441. Introduction
The pile foundation is an important component to bear the upper structures, whose working condition directly affects the safety of the upper buildings or bridges. For existing bridges in service, it is very important to evaluate the capacity of railway bridge piles, especially when the train axial load or train speed needs to increase. Two different types of capacity are needed to distinguish: one is ultimate bearing capacity which is usually controlled by the pile strength and the other is service capacity which is usually controlled by the pile settlement. For most friction piles designed as railway bridge foundations, the service capacity is paid more attentions and is determined by the maximum settlement Sa of the pile head [1]. Static load test (SLT) is a reliable method to evaluate bearing capacity of pile foundations [2]; nevertheless, SLT cannot be employed for testing existing bridge piles in service. Then, the dynamic measurement is a good choice to estimate the service capacity. The transient response method (TRM), also known as the mechanical mobility method, is a dynamic loading method for the pile inspection. It was proposed in the 1970s by Davis and Dumm [3]. The application of TRM in pile foundation test has been carried out by many researches. Rausche et al. [4] compared the pulse echo and TRM and recommended the analyses should be made in both the frequency and time domains. Tham et al. [5] analysed the transient response of single pile under vertical load in layered soil by using the quasilinear time domain boundary element method. By using dynamic test, both low-strain and high-strain testing have been widely used in the pile nondestructive examination [6–9]. The pile dynamic response under vertical loads was also analysed by using various dynamic interaction models [10–13].
TRM analyses both the velocity and force signals in the frequency domain. The velocity spectrum is divided by the force spectrum to determine the mobility or mechanical admittance spectrum [14], which helps provide more information than the traditional pulse echo method to identify defects near the top of the pile [15]. Some key information from the graph, such as the peak/mean mobility ratio, mobility, and damping, is widely used to evaluate the pile integrity and pile length [16–19]. Another important parameter from mobility spectrum is the dynamic stiffness. It is defined as the slope of the low frequency linear portion of the spectrum from the origin to the first peak. The dynamic stiffness Kd can be calculated by(1)Kdf=2πfVf/Ff,where V(f) and F(f) are the velocity and force signals in the frequency domain. The value of Kd is sensitive to the stiffness of the pile shaft under compression. When the frequency approaches to 0 (f→0), the value of the dynamic stiffness approaches to the static stiffness (Kd→Ks). In practice, however, the frequency of the dynamic impulse cannot be 0 Hz. Therefore, a coefficient α is introduced here to describe the ratio between the dynamic and static stiffness: α=Kd/Ks. Then, the pile allowable bearing capacity Q can be estimated by(2)Q=KdSaα,where Sa is a guideline value of the pile settlement. The coefficient α is a bridge between the dynamic test and static capacity estimation; however, it is an empirical value rather than a theoretical one. This capacity estimation method was proposed in China in 1980s [20] and used for pile capacity evaluation [21]. Liu and Ma [22] analysed the parameter sensitivity on single pile and found that a reasonable dynamic stiffness can be used as an alert value for pile capacity. Ma et al. [23] measured the dynamic stiffness for 680 bridge piles and found that an obvious positive correlation exits between the dynamic stiffness and bearing capacity of intact piles. Nevertheless, the existing researches did not consider the influences of pile cap or platform on dynamic stiffness. For a railway bridge with a cap-pile foundation, it was not clear whether the capacity estimation method can still be employed as the existence of the pile cap. Accordingly, in this paper, a laboratory test was performed to learn the dynamic stiffness of the cap-pile system. Then a numerical model was also built for this type of system.
2. Laboratory Test
In order to analyse the relationship between dynamic stiffness and capacity of cap-pile system, a laboratory test was designed. As we mainly care about the basic principle and characteristic of dynamic stiffness, the similarity theory was not considered here. The cap-pile model was made of plaster, and it was then embedded in a sand box. For each foundation, there were four piles under the cap. All the piles were designed as friction piles considering the pile-soil interaction. The length of pile was 0.7 m with a diameter of 0.075 m. The schematic diagram of the model was shown in Figure 1.
The picture and the dimension of cap-pile model (unit: cm).
The dynamic stiffness test was based on the impulse transient response method. The hand hammer was used to impact the cap top. Then, the responses were obtained by the acceleration sensors installed on the top of the cap, and the force curve is obtained by force sensors installed on the hand hammer. The velocity responses can be calculated by integration. The dynamic stiffness can be obtained by (1). In order to ease the wave reflection on the steel plate of the sand box, five polystyrene foam cushions with a thickness of 0.02 m were installed on the bottom and four sides of the sand box (Figure 2).
Laboratory test sand box.
For friction piles, the pile capacity is strongly dependent on the effective friction area around the piles. Therefore, in this test, the effective friction area between piles and soil was designed to decrease in order to model the reduction of capacity. Then, the relation between dynamic stiffness and effective friction area A was analysed.
Suppose the initial embedded depth of piles was H, designed as H=70 cm, and the embedded depth decreased by d=10 cm each time. In Figure 3, it can be found that, with the decrease of the embedded depth of the pile foundation, the value of A was gradually reduced, and accordingly the bearing capacity of the system was also decreased. The dynamic stiffness as a function of frequency was shown in Figure 4.
Sketch of decrease of the effective friction area between piles and soil. The embedded depth of piles was designed as (a) H, (b) H-d, (c) H-2d, and (d) H-3d.
Amplitude frequency curve of dynamic stiffness at low frequency of centre of the cap.
It can be found that, with the decrease of A, dynamic stiffness decreases in all frequencies below 50 Hz. In Figure 5(a), it can clearly show that Kd decrease with A at four different frequencies to be analysed: 5 Hz, 10 Hz, 15 Hz, and 20 Hz. The normalized dynamic stiffness by H=70 cm or A=6597 cm2 was shown in Figure 5(b). At these frequencies, the general attenuation gradients of Kd were similar overall but were different for each A. The normalized dynamic stiffness decayed largest with the value of 35% when the first layer of sand was relieved. When the subsequent soil layers were relieved, the normalized dynamic stiffness decayed gradually with the value of less than 10%. That was the dynamic stiffness largely affected by the top soil.
(a) Dynamic stiffness and (b) normalized dynamic stiffness change with effective friction area.
3. Numerical Analysis3.1. Finite Element Model
In this section, a 3D dynamic finite element (FE) model of cap-pile-soil was established for a bridge foundation (Figure 6). The dimensions of the piles and cap in the calculation are from a tram bridge. There were also four friction piles under the cap, the same as the model in the laboratory test. The length and the diameter of each pile were 15.5 m and 1.5 m, respectively. The dimension of the cap was 6 m × 6 m × 2 m. The centre distance between the four piles along the long side of cap was 3.5 m. The concrete strength class was C25, whose modulus was 2.5 × 1010 Pa, the Poisson’s ratio 0.16, and the bulk density 2300 kg/m3. According to the geological investigation report, the soil was simplified as five layers, and the detail parameters were listed in Table 1, where H was the depth of the soil layer, ρ the density, c the cohesion, φ the internal friction angle, cs and cp the S-wave and P-wave velocity and v the Poisson’s ratio.
Soil parameters.
Soil type
H/m
ρ/(kg/m3)
cs/(m/s)
cp/(m/s)
Poisson’s ratio v
Mucky soil
15.0
1660.2
143.8
285.5
0.33
Silt sand
12.5
1850.0
145.3
288.5
0.33
Fine sand
10
1960.2
190.0
333.6
0.26
Medium coarse sand
7.5
1979.6
240.0
415.7
0.25
Silty clay
5.0
1929.6
181.6
333.9
0.29
FE model of cap-pile-soil system.
The impact load was applied at the centre of the cap and the response at the cap 0.5 m away from the cap centre was analysed. The time history and Fourier spectrum of the dynamic load were shown in Figure 7. The Rayleigh damping assumption was employed. The damping matrix C is defined as a linear superposition of mass matrix M and stiffness matrix K with the coefficients a1 and a2. The two coefficients can be only defined by the damping ratio of soils and the analysis frequencies.
(a) Time history and (b) Fourier spectrum of the impact load.
To avoid the wave reflection at the model boundaries, the artificial viscoelastic boundaries were applied, which were widely used together with FEM [24, 25]. Then, the surface springs were created at four side boundaries and the bottom boundary of the soil layers. The relationship between the spring stiffness and the elastic modulus E can be estimated by an empirical formula:(3)kv=130EAv30-3/4,kh=130EAh30-3/4,where E is the elastic modulus; Av and Ah are the cross-sectional area in the normal and tangential directions at the boundaries.
3.2. Effect of Effective Friction Area
As the capacity of friction piles was strongly related to the effective friction area around the piles. Similar to Section 2, here in the FE model by changing the thickness of the soil around the pile foundation, the effective friction area was changed.
Suppose that the initial embedded depth of piles was H, designed as H=15 m, and the embedded depth decreased by d=2.5 m each time (Figure 8).
Sketch of decrease of the effective friction area between piles and soil. The embedded depth of piles was designed as (a) H, (b) H-d, (c) H-2d (d) H-3d, and (e) H-4d.
H=15 m
H-d=12.5 m
H-2d=10 m
H-3d=7.5m
H-4d=5 m
Figure 9 shows the velocity admittance of the cap response. It can be found that, with the decrease of the effective friction area, the first natural frequency was also decreased. In the low frequency of the admittance curve, the linear regression analysis was carried out (Figure 10). It can be found that when the embedded depth of pile foundation was 15 m, the slope of the linear fitting curve for the low frequency can be determined to be k=6.671×10-10 N/m. Then, the dynamic stiffness of the centre measuring point at the top of the cap can be calculated as Kd=9.419×109 N/m.
The time history and velocity admittance of the cap response.
The linear fitting curve of low frequency curve of admittance curve.
In Figure 11, it can be clearly observed that the normalized dynamic stiffness by H=15 m decayed with the effective friction area. The simulation result was similar to the test result in Section 2.
The normalized dynamic stiffness change with effective friction area.
3.3. Comparison of Single Pile and Cap-Pile System
To analyse the dynamic stiffness relationship between the single pile and the cap-pile system, four FE models were built (M1~M4) (Figure 12). In M1, there was only single pile; in M2, a single pile with a small cap was built; in M3, two piles under the cap were considered; finally, in M4, a cap-pile system with four piles was built. The piles and cap dimensions were listed in Table 2
The pile and cap dimensions of the four FE models.
Pile and cap
M1
M2
M3
M4
Pile
Diameter/m
1.5
1.5
1.5
1.5
Length/m
15.5
15.5
15.5
15.5
Cap
Dimension: length × width × height/(m × m × m)
—
2.5 × 2.5 × 2.0
6.0 × 2.5 × 2.0
6.0 × 6.0 × 2.0
FE models of the cap-pile-soil system.
The impact load was applied at the centre of the cap (or pile in M1) and the responses on the cap head were analysed. The time history and Fourier spectrum of the dynamic load were shown in Figure 7. The velocity admittance and dynamic stiffness curves of the four models were shown in Figure 13. It can be found that, in the frequencies below 30 Hz, models M1 and M2 have similar values of dynamic stiffness. The existence of small cap induces the decrease of the value. Model M4 has the largest dynamic stiffness, which is about four times as large as M1 and M2 between 10 and 20 Hz. Model M3 has the second largest dynamic stiffness, which is about twice as large as M1 and M2.
The (a) velocity admittance and (b) dynamics stiffness of four models.
4. Conclusion
In order to analyse the dynamic stiffness of cap-pile system, both a laboratory test and numerical simulation were performed in this paper. The results show that the dynamic stiffness obtained from the pile mobility curve is a sensitive index and can well reflect the static characteristics of the pile at the elastic stage. In addition, there is a significant positive correlation between the vertical dynamic stiffness index and bearing capacity of the cap-pile system in the similar formation environment. Finally, for the cap-pile system with four piles, the dynamic stiffness is about four times as large as the single pile between 10 and 20 Hz.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors wish to acknowledge the support and motivation provided by the Fundamental Research Funds for the Central Universities (no. 2016JBM034) and the Funding of College Student Scientific Research and Training Projects of Beijing Jiaotong University (no. 160130044).
JozefiakK.ZbiciakA.MaslakowskiM.PiotrowskiT.Numerical modelling and bearing capacity analysis of pile foundationProceedings of the 24th Russian-Polish-Slovak Seminar on Theoretical Foundation of Civil Engineering, TFoCE 2015August 2015Russia35636310.1016/j.proeng.2015.07.1012-s2.0-84947474951BudiG. S.KosasiM.WijayaD. H.Bearing capacity of pile foundations embedded in clays and sands layer predicted using PDA test and static load testProceedings of the 5th Euro Asia Civil Engineering Forum Conference, EACEF 2015September 2015Indonesia4064102-s2.0-8496090980410.1016/j.proeng.2015.11.101DavisA. G.DummC. S.From theory to field experience with the non-destructive vibration testing , Proceeding of Institution of Civil Engineers2, no. 57Proceedings of the Part 21974571593RauscheF.RenkungS.LikinsG. E.Comparison of pulse echo and transient response pile integrity test methods19912127ThamL. G.ChuC. K.LeiZ. X.Analysis of the transient response of vertically loaded single piles by time-domain BEM19961921171362-s2.0-004312686810.1016/0266-352X(95)00064-HNiS.-H.LehmannL.CharngJ.-J.LoK.-F.Low-strain integrity testing of drilled piles with high slenderness ratio2006336-728329310.1016/j.compgeo.2006.08.0012-s2.0-33750682393PaquetJ.BriardM.Nondestructive quality control of concrete piles19984980YanlongL.ShouyiL.YangY.XingT.Temperature stress and surface insulation measures of concrete face slabs during cold wave period20151345015072-s2.0-84949434289LiY.SunY.LiB.XuZ.Penalty function-based method for obtaining a reliability indicator of gravity dam stability201781192510.1016/j.compgeo.2016.07.003WuW. B.WangK. H.ZhangZ. Q.LeoC. J.Soil-pile interaction in the pile vertical vibration considering true three-dimensional wave effect of soil20133717286028762-s2.0-8488782486110.1002/nag.2164RandolphM. F.DeeksA. J.Dynamic and static soil models for axial responseProceedings of the in Proceedings of the 4th International Conference on the Application of stress wave Theory to PilesSeptember 1992The Hague, The Netherlands314MichaelidesO.GazetasG.BouckovalasG.ChrysikouE.Approximate non-linear dynamic axial response of piles199848133532-s2.0-003200740910.1680/geot.1998.48.1.33YesilceY.CatalH. H.Free vibration of piles embedded in soil having different modulus of subgrade reaction20083258899002-s2.0-3834911026810.1016/j.apm.2007.02.015Zbl1130.74032LiangL.BeimJ.Effect of soil resistance on the low strain mobility response of piles using impulse transient response methodProceedings of the in Proceedings of the 8th International Conference on the Application of Stress Wave Theory to PilesSeptember 2008Lisbon, Portugal435441RauscheF.LikinsG.HusseinM. H.Analysis of Post-Installation Dynamic Load Test Data for Capacity Evaluation of Deep FoundationsProceedings of the Symposium Honoring Dr. John H. Schmertmann for His Contributions to Civil Engineering at Research to Practice in Geotechnical Engineering Congress 2008New Orleans, Louisiana, United States31233010.1061/40962(325)8DavisA. G.The nondestructive impulse response test in North America: 1985–2001200336418519310.1016/s0963-8695(02)00065-82-s2.0-0037409463LiaoS. T.RoessetJ. M.Dynamic response of intact piles to impulse loads19972142552752-s2.0-003112628510.1002/(SICI)1096-9853(199704)21:4<255::AID-NAG869>3.0.CO;2-JLoK. F.NiS. H.HuangY. H.Non-destructive test for pile beneath bridge in the time, frequency, and time-frequency domains using transient loading2010621-234936010.1007/s11071-010-9721-7Zbl1359.74318MassoudiN.TefferaW.Non-destructive testing of pile using the low strain integrity methodProceedings of the in Proceedings of the 5th International Conference on Case Histories in Geotechnical EngineeringApril 2004New York, NY, USAXuY. Z.Determination of pile bearing capacity by dynamic test19921417483ZhengC. T.WangY.XuS. J.et al.Pile bearing capacity measurement by steady mechanical impedance method19972885052LiuJ.MaM.Analysis of the dynamic stiffness and bearing capacity for pile foundationsProceedings of the International Conference on Vibroengineering - 2015September 2015chn1341392-s2.0-84954237199MaM.LiuJ.KeZ.GaoY.Bearing Capacity Estimation of Bridge Piles Using the Impulse Transient Response Method201620162-s2.0-8495545961910.1155/2016/41870264187026WanmingZ.ZhenxingH.XiaolinS.Prediction of high-speed train induced ground vibration based on train-track-ground system model2010945455542-s2.0-7875160842410.1007/s11803-010-0036-yXuQ.XiaoZ.LiuT.LouP.SongX.Comparison of 2D and 3D prediction models for environmental vibration induced by underground railway with two types of tracks2015681691832-s2.0-8492855523410.1016/j.compgeo.2015.04.011