http://orcid.org/0000-0002-7199-4693
Figures
Abstract
In this paper, a method was developed to assess the vertical response of piles considering crushing characteristics. A load transfer model of a pile is built and the relationship between the crushing characteristics and load transfer function of a pile was analyzed. Based on the analysis results, the transfer equation of load function considering the crushing characteristics was derived using the finite difference method. The ideal side resistance and the load transfer mechanism of a pile can be obtained by the transfer equation of load function. The predicted response of piles with the proposed method is in good agreement with the response observed in some pile tests reported in the literature. Further study was conducted to discuss the effect of the parameters related to the load transfer model considering crushing characteristics on the load-settlement response of a single pile.
Citation: Jiang C, Li Y, He J-l (2019) Analysis of vertically loaded piles considering crushing characteristics of crushed stones. PLoS ONE 14(7): e0219003. https://doi.org/10.1371/journal.pone.0219003
Editor: Jianguo Wang, China University of Mining and Technology, CHINA
Received: March 9, 2019; Accepted: June 10, 2019; Published: July 11, 2019
Copyright: © 2019 Jiang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and Supporting Information files.
Funding: This work is supported by the National Natural Science Foundation of China (Grant No. 51678570;51478479).
Competing interests: The authors have declared that no competing interests exist.
Nomenclature:
qs,
the ultimate side resistance; qb,
the ultimate base resistance; μ,
the axial displacement; μb,
the base axial displacement; τ,
the side resistance; q,
the base resistance; EM,
the Ménard modulus; B,
the pile diameter; Br,
the relative breakage; σ3,
the rolling load; n,
the number of segments of pile; φp,
the peak internal friction angle; c,
the cohesion of crushed stones; σ,
the horizontal stresses; k0,
the lateral pressure coefficient; γ,
the crushed stones bulk density; z,
The depth of calculated point; Pa,
the atmospheric pressure; Ei,
the Young’s modulus; Ai,
the section area of the pile; U,
the circumference of pile; L,
the pile length; h,
the length of each of segment; ,
the derivative of μi; μi,
the displacement of the i-th node; df,
a differential of “percent passing” divided by 100; bp,
the potential for breakage that is significant to soil behavior; a,b,j,k,
the fitting parameters of the particle; bp0,
the original values of bp; bpl,
the original values of bp after loading; μib,
the base displacement of i-th node; μit,
the top displacement of i-th node; αs,
the factor taking into account the influence of the soil type; αb,
the factor taking into account the influence of the rock type; ks,
the initial stiffness of side resistance; kb,
the initial stiffness of base resistance; λs,
the parameter taking into account the influence of the crushed stones type; λb,
the parameter taking into account the influence of the rock type; τ(μ),
the function describing the variation of the mobilized side resistance τ in function of the axial displacement μ; q(μb),
the function describing the variation of the mobilized base resistance q in function of the axial displacement μb
Introduction
Karst caves are inevitably encountered during tunnel construction in western China. In order to strengthen the bearing capacity of rock-socketed pile foundations in the karst cave below the tunnel, the rock-soil aggregate from excavation is backfilled around the rock-socketed pile foundations. However, the effect of backfill is ignored in engineering practice because the design is too conservative. It was found that the pile side resistance from field tests bore most of the total load and increased with the increase of the pile modulus ratio [1]. Therefore, it is of great significance to further study the effect of backfill on the bearing mechanism of the pile foundation.
In past decades, a number of methods have been developed for pile response analysis, included the theoretical load transfer curve method [2–7], the shear displacement method [8–12], finite element method [13–16], the boundary-element method [17–20], the variational approach [21–24] and the field test [25–27], etc. It is worth mentioning that these studies reported that the pile was surrounded by clay, sand or rock-soil aggregate. However, the pile foundation of a tunnel is surrounded by crushed stones in many cases, especially in the case of backfilling. The load—settlement curve of backfill piles in the large karst area is developed by Jiang et al. [28]. Chen et al. [29] proposed a method for the calculation of the load settlement of karst cave backfill pile foundations under low confining pressure. In the process of backfilling, crushed stones will be rolled to make it more compact. However, some research [28–29] failed to consider that crushing is often caused by a rolling load in the backfilling process, and the crushing has great influence on the mechanical properties of crushed stones [30–33]. Therefore, more extensive studies are warranted to assess the vertical response of piles considering crushing characteristics.
In light of this, this study investigates the relationship between the crushing characteristics of crushed stones and the load transfer of the pile foundation, and then analyzes the load transfer mechanism of piles caused by crushing characteristics. The effect of crushing characteristics on the load transfer response of the pile surrounded by crushed stones is also discussed. The results calculated by the proposed method in this paper are compared with the measured results from other research.
Establishment of the pile foundation model in crushed stones backfill
Establishment of the transfer function and assumptions
Generally, researchers study the vertical response of piles using transfer function method, shear displacement method, elastic theory method and numerical method. The primary difference in various methods is the assumption of the pile-soil model. The transfer function method is simple for engineering practice. In fact, the key to the calculation method of rock-socketed pile foundations in crushed stones backfill area (Fig 1) is the selection of the transfer function. A linear model was usually used to calculate the settlement of vertically loaded pile. However, the measured results in some engineering practice show that the load transfer function curves of pile side were generally nonlinear [2, 34, 35]. Subsequently, an exponential model was built to simulate the nonlinear variation of the pile side resistance. The stiffness of the soil around the pile was assumed in the ‘AB1’ model [7] (Fig 2 and Eq (1)).
In a recent study, the complexity associated with characteristics of crushed stones was not thoroughly studied, and the flaw rested with the mobilization of the pile side resistance for crushed stones [28]. It is worthwhile to build a load transfer model considering breakage of crushed stones around the pile. The side resistance and the base resistance is obtained by the integration of Eq (1).
(2)
(3)
where qs, qb, μ and μb are, respectively, the ultimate side resistance, the ultimate base resistance, the axial displacement and the base axial displacement. λs and λb are the parameters taking into account the effect of the crushed stones type and rock type, respectively. t–z curves come from the function τ(μ) which describes the variation of the mobilized side resistance τ in the function of the axial displacement μ. And the function q(μb) describing the variation of the mobilized base resistance q in the function of the base axial displacement μb is obtained.
(4)
(5)
Where B is the pile diameter, EM is the Ménard modulus, αs and αb is the factors taking into account the effect of the soil type. αs is 2.0 and αb is 9.0 in chalk and crushed stones [7, 28]. When μ = 0 and μb = 0, ks and kb are the initial stiffness of the pile side and pile base respectively.
The morphology and the geometrical conditions that are assumed for the pile in backfill area and the detail stresses diagram in the pile side are shown in Fig 3. B is the pile diameter. σ3 is the rolling load and τ is the mobilized side resistance. σc and σr are the horizontal stresses that are applied on the pile in the crushed stones and rock, respectively. Hc and Hr are the depth of the pile in the crushed stones and rock, respectively.
To formulate pile foundation model in backfill area considering the crushing characteristics, in terms of the interaction between the pile and crushed stones, assumptions and statements are made in advance.
- In the same horizontal plane of pile, the settlement of the crushed stones backfill is the same, and the corresponding displacement is the same. There is only vertical displacement for a pile under vertical loading and the radial displacements for a pile is negligible.
- The crushed stones and rock around the pile are considered as a series of independent nonlinear springs, and the pile is modeled as an elastic-plastic material.
- The peak internal friction angle of crushed stones is assumed to be the internal friction angle of the Mohr-Coulomb criterion. Before the crushed stones are backfilled, the rolling load is equivalent to the confining pressure of the triaxial test of compacted crushed stones.
Based on the above assumptions, the relationship between the function τ(μ) of the crushed stones, rock of pile side, and the function q(μb) of rock of the pile base and corresponding displacement is shown in Fig 4.
Relative breakage
Up to now, breakage often takes place when crushed stones are compressed or sheared [36–37]. The breakage (Bg) proposed by Marsal [38] is the most commonly used parameter in particle breakage. Bai [39] proposed B60 to estimate crushing characteristics by measuring the particle size difference before and after loading. The relative breakage Br, proposed by Hardin [30], is stable and can also reflect the size variation of each particle before and after the test. Therefore, the relative breakage Br of crushed stones is used to measure the crushing characteristics.
(8)
Where bp represents the potential for breakage that is significant to soil behavior; df is a differential of “percent passing” divided by 100; bp0 is the original values of bp; bpl is the original values of bp after loading.
Crushed stones compaction is warranted to improve the pile bearing capacity. The effect of rolling loads on the crushed stones breakage depends on the types of rolling load. Based on the crushing test of equal-pressure consolidation particles under different confining pressure by large Triaxial tester, the relationship between breakage Br and the rolling load σ3 was obtained [31]. The plots show a good linear relationship. By fitting the formula of the linear relationship, the relationship between the relative breakage Br and the rolling load σ3 are given:
(9)
Where a, b are fitting parameters, which are obtained by fitting the test results of
, and Pa is atmospheric pressure to convert relative breakage Br into dimensionless quantity.
Based on the fitting experimental data, the relationship between the peak internal friction angle φp and the relative breakage Br was obtained by Indraratna [32]:
(10)
Where j and k are test parameters of the particle. The fitting result is shown in Table 1 [31]:
Ultimate side resistance
The Eqs 9 and 10 show that the rolling load and the relative breakage affect the internal friction angle of crushed stones. The internal friction angle is related to the ultimate side resistance of crushed stones. Thus, it is necessary to analyze the effect of relative breakage on the ultimate side resistance. Triaxial tests on rock particles were carried out to determine the strength characteristics of the particle under different horizontal stress [33]. The cohesion, internal friction angles, Mohr’s circle and the shear strength envelope of the samples under different horizontal stress are in good agreement with the Mohr-Coulomb strength calculation formula, as shown in Fig 5.
(11)
Where c, σ and φp are the cohesion of crushed stones, the horizontal stress and the peak internal friction angle of crushed stones, respectively.
The shear plane occurs on the vertical plane and the horizontal stresses on the pile side varies with depth. Therefore, the horizontal stress can be expressed as Eq (12).
(12)
Where k0, γand z are the lateral pressure coefficient, the crushed stones bulk density and calculated point depth, respectively; qsc and qsr are the ultimate side resistance at the horizontal stress σc and the ultimate side resistance at the horizontal stress σr, respectively.
Analysis of vertical load transfer
Governing equation and solution
As shown in Fig 6, the whole length L of pile is divided into n segments by the finite difference method. Each segment is h = L/n. The node of pile top and pile base are denoted as 0 and n, respectively. The i-th segment is subjected to the side resistance τ along the pile. The force of the upper segment and the lower segment are represented by Pit and Pib, respectively. Based on the force equilibrium of the circular pile segment located at depth z, Eq (13) is obtained:
(13)
Where Ei, U, and Ai are the Young’s modulus, the circumference and the segment area of the pile, respectively. Eq (14) is obtained by substituting the load transfer function Eq (2) and the ultimate side resistance expression Eq (11) into Eq (13).
In this model, the boundary conditions between each segment are as follows:
(15)
Where b and t correspond to the base and top of the section, respectively; i is the i-th section; the displacement of the i-th node is μi;
is the derivative of μi.
According to the base axial displacement μb, Eq (15) can be transformed into the two forms of boundary conditions as follows.
(2) μb >0
(17)
Where s0 is the settlement of the pile top; P is the load applied on the pile top.
A program was written to compute the Q-S curves, side resistance and axial force according to boundary conditions. The program flowchart is shown in Fig 7.
Comparison with other research
To validate the proposed model for the analysis of pile response, the proposed method in this paper is applied in this section to some pile load tests, which are well documented in the literature. The computed pile responses are compared to the pile test results. The pile tests considered are the tests reported by Carrubba [40] (Case 1) and Radhakrishnan [41] (Case 2).
Case1. The first pile tests were performed to analyze the side resistance of large-diameter rock-socketed piles [40]. The pile diameter is 1.2 m; the pile length is 18.5 m; the pile depth in the granular soil and rock are 11m and 7.5 m, respectively; the Young’s modulus of the pile Eb is 31500MPa. It is difficult to determine the fitting parameters of breakage. Therefore, the value of qs is used to verify the example. Based on the measured pile load-settlement curves, the back-calculated analysis and numerical method were used to determine the pile side resistance [42]. The parameters of the exponential load transfer function are summarized in Table 2. The pile responses calculated by the proposed method are compared to the pile test results [40] and to the pile responses computed by the method from Xiao [42], which is shown in Fig 8.
As illustrated in Fig 8, the methods in this paper and from Xiao [42] consider the nonlinear mobilization of side resistance. The results computed by the method from Xiao [42] are in good agreement with the measured results before the pile settlement of 9 mm while the results from the proposed method in this paper compare well with the measured results during the loading tests.
Case 2. In the pile tests conducted by Radhakrishnan [41], the pile TP1 was driven into highly weathered and fragmented siltstone. The pile diameter is 0.81m and the length of the pile is 11.5m. The unit shaft friction/shaft movement relationship with buried depth of 2m, 4m,8m and 10m is fitted by the parameters of the first, second, third and fourth layers, respectively. The parameters of the pile base are obtained by combing the figure of the variation of load versus depth with the load-settlement curve. The parameters of the exponential load transfer function are given in Table 3.
The load distribution and side resistance versus depth relationships are shown in Fig 9 for different vertical load. The predicted results are in good agreement with field tests in Fig 9A. However, the predicted side resistance is inconsistent with the measured result when Q≥6000KN. The small discrepancy may be attributable to the decreased shear displacement of the interface which inhibits the mobilization of side resistance of the pile and then induces that the EM value of predicted result is higher than that of the measured result. As shown in Fig 10, a discrepancy from the comparison may be caused by the difference of side resistance between measured results and predicted results. Therefore, the realistic Q-S curve is below the predicted curve at the same settlement.
(a) Load distribution versus depth relationship; (b) side resistance versus depth relationship.
Parametric study
The parameters of the crushing characteristics of crushed stones include lateral pressure coefficient, rolling load, relative breakage and the lithology of crushed stones. The effect of crushing characteristics on pile bearing capacity is analyzed with piles of 1.2m diameter, 14m depth in the crushing stones and 4.5m depth in the rock. The parameters for crushed stones are αs = 2.0, EM = 120MPa, c = 30KPa and γ = 20.7KN/m3. The parameters for the rock are qs = 135KPa,αs = 2.0 and EM = 230MPa. The parameters for the pile base are qb = 5000 KPa, αb = 9.0, EM = 330MPa and the Young’s modulus Eb = 31500MPa.
Effect of lateral pressure coefficient
The lateral pressure coefficient determines the mobilization of the horizontal stress around the pile. Figs 11–14 show the Q-S curves, τ-μ curves, τ/qs-z curves and P/P0-z curves of the pile in backfill area with different lateral pressure coefficient. The lithology of sample 2 is limestone and the rolling load is σ3 = 1500KPa.
As shown in Fig 11, with the increase of the lateral pressure coefficient, the ultimate bearing capacity of the pile becomes larger and the pile settlement decreases at the same load. When the ultimate bearing capacity of the pile is reached, the lateral pressure coefficient increases by 0.1 and the ultimate bearing capacity increases by 9.22% and 10.67%. Fig 12 shows that the increase of lateral pressure coefficient directly increases the ultimate side resistance. The lateral pressure coefficient increases by 0.1 and the ultimate side resistance increases by 33.33% and 33.33%. And the displacement μ increases when the ultimate side resistance is reached. In Figs 13 and 14, when the lateral pressure coefficient increases, the τ/qs ratio and the P/P0 ratio decreases with the increase of pile depth z. When the τ/qs ratio of the embedded rock section is stable, the τ/qs ratio of the backfill section continues to increase. When the vertical load reaches 9000KPa, the lateral pressure coefficient increases by 0.1 and the τ/qs ratio decreases by 22.98% and 15.31%, and the P/P0 ratio decreases by 5.07% and 3.12%.
Effect of the rolling load and the relative breakage
The relative breakage Br of sample 2 was 5.0%, 12.5%, and 20.0%, respectively. The rolling load was 500KPa, 1500KPa, and 2500KPa. The lateral pressure coefficient is k0 = 0.4. Fig 15A and 15B show the Q-S curves of piles in crushed stones backfill area with the relative breakage of the crushed stones and the rolling load change, respectively.
(a) Effect of the relative breakage on Q-S curves; (b) Effect of the rolling load on Q-S curves.
It can be seen that the ultimate bearing capacity of the pile decreases with the increase of the rolling load and relative breakage in Fig 15. This trend becomes inconspicuous with the increase of the rolling load and the relative breakage. The decreases ratio by 10.63%, 4.8% in relative breakage, respectively, and the decreases ratio by 7.26%, 1.84% in the rolling load, respectively. This is because the increase in the relative breakage of the crushed stones reduces the pile side resistance. Fig 16 shows that the relative breakage increases by 7.5% and the ultimate side resistance decreases by 33.6% and 16.26%, and rolling load increases by 1000KPa and the ultimate side resistance decreases by 21.62% and 6.18%. As shown in Fig 17, the relative breakage increases and the growth rate becomes slow with the increase of the rolling load. Pile side resistance decreases as the relative breakage of crushed stones increases. In addition, the relative breakage of sample 1 and sample 2 are very close when the rolling load reaches 2500KPa.
(a) Effect of the relative breakage on τ-μ curves; (b) Effect of the rolling load on τ-μ curves.
Figs 18 and 19 show that the side resistance increases gradually as the relative breakage and the rolling load increase. However, the effect of the relative breakage and the rolling load on the mobilization of the side resistance is gradually insignificant as the relative breakage and the rolling load increase. At the vertical load of 9000KN, the τ/qs ratio increases initially and then decreases as the depth increases. Because the axial displacement decreases and the ultimate shear strength increases with the depth, the maximum side resistance is reached at a certain depth. Generally, the increase of the relative breakage and the rolling load promotes the mobilization of the side resistance, which makes the axial load shorn off more easily.
(a) Effect of the relative breakage on τ/qs-z curves; (b) Effect of the rolling load on τ/qs-z curves.
(a) Effect of the relative breakage on P/P0-z curves; (b) Effect of the rolling load on P/P0-z curves.
Effect of crushed stones lithology
Figs 20–23 show the Q-S curves, τ-μ curves, τ/qs-z curves, and P/P0-z curves of the pile in backfill area with different crushed stones lithology, respectively. The lithologies for eight samples are quartz sandstone, Limestone et.al. The lateral pressure coefficient is 0.4 and the rolling load is 1500KPa.
It is obvious that there is little difference between the eight samples and the ultimate load of the pile is about 10600KN in Fig 20. As shown in Fig 21, the ultimate side resistance of the eight samples ranged from 90KPa to 105KPa and the ratio of the length of the interval to the maximum ultimate side resistance is 0.143. It illustrates that the effect of different samples on the ultimate side resistance is very small. Because the relative breakage and the peak internal friction angle curve in the literature [31] have little change in lithology, the ultimate side resistance is less affected by lithology. From Figs 22 and 23, it can be seen that the τ/qs ratio and the P/P0 ratio are not changed by the lithology of the crushed stones.
Conclusion
A load transfer model for the rock-socketed pile considering the relationship between the crushing characteristics and the nonlinear load transfer function was built in the large karst cave. Using the finite difference method, the transfer equation of load function was derived. In addition, the effect of the key parameters including the lateral pressure coefficient, the relative breakage, the rolling load and the crushed stones lithology on the bearing mechanism of the rock-socketed pile in the crushed stones backfill was analyzed. Main conclusions can be summarized as follows.
- A load transfer function model was built to reflect the load-settlement behavior of the rock-socketed pile in the crushed stones backfill area. The proposed model was verified by the comparison of two cases.
- The lateral pressure coefficient has a significant effect on the ultimate bearing capacity and the mobilization of side resistance of pile. The ultimate side resistance and the ultimate bearing capacity gradually increase as the lateral pressure coefficient increases. The side resistance increases initially and then decreases as the depth increases at high load levels.
- The increase in relative breakage and the rolling load increases the ultimate bearing capacity of the pile and the mobilization of side resistance, this trend becomes not obvious with the increase of the rolling load and the relative breakage. At the vertical load of 9000KN, the τ/qs ratio increases initially and then decreases as the depth increases. The crushed stones lithology has little effect on the ultimate bearing capacity (about 10600KN) and the ultimate side resistance (ranging from 90KPa to 105KPa) by the comparison of eight samples. Thus, to improve the pile bearing capacity, the relative breakage of crushed stones should be reduced.
- By analyzing the Q-S curves, τ-μcurves, τ/qs-z curves, and P/P0-z curves, it can be concluded that the side resistance of the pile in the crushed stones backfill area cannot be ignored in term of avoiding the waste of engineering materials. If the test parameters of the crushed stones and the rolling load are obtained, the relative breakage and the bearing capacity of pile can be estimated. The present study proposes a method to assess the vertically loaded pile response considering crushing characteristics. The analysis can be further extended to account for characteristics of crushed stones and to study on combined loads.
References
- 1. Rosenbeng P, Journeaux NL. Friction and end bearing tests on bedrock for high capacity socket design. Can Geotech J 1976;13(3):324–333.
- 2. Kraft L. M, Ray R. P, Kagawa T. Theoretical t-z curves. J Geotech Eng 1981;107(11):1543–61.
- 3. Liu J, Xiao HB, Tang J, Li QS. Analysis of load-transfer of single pile in layered soil. Comput Geotech 2004;31(2):127–135.
- 4. Zhang Q, Zhang Z, He J. A simplified approach for settlement analysis of single pile and pile groups considering interaction between identical piles in multilayered soils. Comput Geotech 2010;37(7–8):969–976.
- 5. Li J, Tan Y, Liang F. A modified analysis method for the nonlinear load transfer behaviour of axially loaded piles. KSCE J Civ Eng 2012;16(3):325–333.
- 6. Zhang Q, Li L, Chen Y. Analysis of Compression Pile Response Using a Softening Model, a Hyperbolic Model of Skin Friction, and a Bilinear Model of End Resistance. J Eng Mech 2014;140(1):102–111.
- 7. Abchir Z, Burlon S, Frank R, Habert J, Legrand S. t–z curves for piles from pressuremeter test results. Géotechnique 2016;66(2):137–148.
- 8. Randolph M. F, Wroth CP. An analysis of the vertical deformation of pile groups. Geotechnique 1979;29(4):423–439.
- 9. Lee CY. Discrete layer analysis of axially loaded piles and pile groups. Comput Geotech 1991;11(4):295–313.
- 10. Guo WD, Randolph MF. An effcient approach for settlement prediction of pile groups. Geotechnique 1999;49(2): 161–179.
- 11. Xie H., Jiang C., He J., & Han H. Analytical Solution for the Steady-State Karst Water Inflow into a Tunnel. Geofluids, 2019, 9.
- 12. Lee KM, Xiao ZR. A simplified nonlinear approach for pile group settlement analysis in multilayered soils. Can Geotech J 2001;38(5):1063–1080.
- 13. Sheng D, Eigenbrod KD, Wriggers P. Finite element analysis of pile installation using large-slip frictional contact. Comput Geotech 2005;32(1):17–26.
- 14. Ma D., Duan H.Y., Liu J.F., Li X.B., Zhou Z.L. The role of gangue on the mitigation of mining-induced hazards and environmental pollution: An experimental investigation. Science of the Total Environment, 2019, 664, 436–448. pmid:30754011
- 15. Said I, De Gennaro V, Frank R. Axisymmetric finite element analysis of pile loading tests. Comput Geotech 2009;36(1–2):6–19.
- 16. Comodromos EM, Papadopoulou MC, Rentzeperis IK. Pile foundation analysis and design using experimental data and 3-D numerical analysis. Comput Geotech 2009;36(5):819–836.
- 17. Mandolini A, Viggiani C. Settlement of piled foundations. Géotechnique 1997;47(4):791–816.
- 18. Ma D., Cai X., Li Q., Duan H.Y. In-situ and numerical investigation of groundwater inrush hazard from grouted karst collapse pillar in longwall mining. Water, 2018, 10(9), 1187.
- 19. Xu KJ, Poulos HG. General elastic analysis of piles and pile groups. Int J Numer Anal Met 2000;24(15):1109–1138.
- 20. Ai ZY, Han J. Boundary element analysis of axially loaded piles embedded in a multi-layered soil. Comput Geotech 2009;36(3):427–434.
- 21. Shen WY, Chow YK, Yong KY. A variational approach for vertical deformation analysis of pile group. Int J Numer Anal Met 1997;21(11):741–752.
- 22. Ma D., Rezania M., Yu H.S., Bai H.B. Variations of hydraulic properties of granular sandstones during water inrush: Effect of small particle migration, Engineering Geology, 2017, 217: 61–70
- 23. Shen WY, Chow YK, Yong KY. Variational solution for vertically loaded pile groups in an elastic half-space. Geotechnique 1999;49(2):199–213.
- 24. Shen WY, Chow YK, Yong KY. A variational approach for the analysis of pile group ± pile cap interaction. Geotechnique 2000;50(4):349–357.
- 25. O'Neill MW, Lecture TTKT. Side Resistance In Piles and Drilled Shafts. J Geotech Geoenviron 2001;127(1):3–16.
- 26. DONG Yun. Experimental study on intensity character of rock-soil aggregate Mixture. Rock and Soil Mechanics 2007;28(6):1269–1274.
- 27. Zhang Q, Zhang Z. Complete load transfer behavior of base-grouted bored piles. J Cent South Univ 2012;19(7):2037–2046.
- 28. JIANG Wu-jun, YAN Ding-yuan, WANG Ming-ming. Reasearch on load-settlement regular pattern and calculation of foundation pile in karst cave. J Geotech Eng 2017;39(S2):67–70.
- 29. CHEN Hua-wei, CHEN Zhao, JIANG Chong, LI Yu. Influences of ultimate frictional resistance at low confining pressure on bearing characteristics of backfill piles. J Geotech Eng 2018;40(S2):87–90.
- 30. Bobby O. Hardin. Crushing of soil particles. Journal of Geotechnical Engineering, American Society of Civil Engineers 1985;111(10):1177–1192.
- 31. LIU Han-long, QIN Hong-yu, GAO Yu-feng, ZHOU Yun-dong. Experimental study on particle breakage of rockfill and coarse aggregates. Soil Mechanics 2005;26(4):562–566.
- 32. Indraratna B., Ionescu D., Christie H. D. Shear Behavior of Railway Ballast Based on Large-Scale Triaxial Tests. Journal of Geotechnical and Geoenvironmental Engineering 1998;124(5):439–449.
- 33. QIN HLYG Hong-yu. Research on strength and deformation behavior of coarse aggregates based on large-scale Triaxial tests. Soil Mechanics 2004;25(10):1575–1580.
- 34. XIAO Zhao-ran, WANG Lu-min, LI Jian-wen. Non-linear behavior of a single pile. J Geotech Eng 2002; (05):640–644.
- 35. Jiang C, Li Y, Liu L, Lin H. Nonlinear Analysis of Flexible Pile near Undrained Clay Slope under Lateral Loading. ADV CIV ENG 2018;1–13.
- 36. Zhang Sheng, Tong C.-X, Xi Li, Sheng Daichao. A new method for studying the evolution of particle breakage. Géotechnique, 2015; 65: 911–922.
- 37. Tong C, Burton GJ, Zhang S, et al. A simple particle-size distribution model for granular materials[J]. Canadian Geotechnical Journal, 2018; 55(2):246–257.
- 38. Marsal RJ. Mechanical properties of rockfill in embankment-dam engineering. textbook. eds. r.c. hirschfeld and s.j. poulos. john wiley and sons inc. pub. ny,1973,92p. International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts 1967;12(4):67.
- 39. Shu-tian BAI, Yi-hao CUI. Mechanics characteristic of rockfill. Journal of Hydroelectricity 1997(03):22–31.
- 40. Carrubba P. Skin friction on large-diameter piles socketed into rock. Can Geotech J 1997;34(2):230–240.
- 41. Radhakrishnan R, Leung CF. Load transfer behavior of rock-socketed piles. J Geotech Eng 1989;115(6):755–768.
- 42. Xiao HB, Luo QZ, Tang J, Li QS. Prediction of load–settlement relationship for large‐diameter piles. The Structural Design of Tall Buildings 2002;11(4):285–293.