Research Article  Open Access
Analysis of the Load Sharing Behaviour and Cushion Failure Mode for a Disconnected Piled Raft
Abstract
Disconnected piled raft (DPR) foundations are characterized by having no structural connection to a raft. The raftpile gap is filled with a granular cushion, which creates a more uniform pressure distribution on the raft and reduces the differential settlement. A series of tests has been performed to investigate the load transfer mechanisms of a DPR. The effects of the thickness and stiffness of the cushion and the pile diameter are presented and discussed. A simple failure mode of the cushion is also put forward according to the cushion slip plane obtained by using the digital image correlation technique that tracks the pixels in the images and generates a displacement field. Then a new theoretical analysis has been presented for rigid pile composite foundation. Through comparative studies, it is found that the present method is validated as reasonable by laboratory tests and is in agreement with several current design methods.
1. Introduction
In situations where a raft foundation alone does not satisfy the design requirements, it may be possible to enhance the performance of a raft by adding piles, producing what is called a piled raft, where the piles are employed as settlement reducers. In the presence of certain structural/geotechnical conditions, the use of a piled raft has been widely adopted as an effective foundation method for designing highrise buildings because of the efficiency of piled rafts in increasing the global bearing capacity and reducing the total and differential settlements, as was first proposed by Burland et al. [1]. As mentioned by Poulos [2], piled rafts embody a possible realization of the socalled creep piling design approach defined by Randolph [3], in which piles work under very high loads. The adoption of piled rafts is also associated with interesting economic and environmental benefits: as reported by Huang et al. [4], this design approach has been widely applied in China, to the extent that designing settlementreducing piles for tall buildings is currently mandatory in Shanghai.
When piles are used as settlement reducers, the objective of the piled raft design is to achieve an economic foundation by minimizing the number of piles while fully utilizing their bearing capacity. Then, high bending moments and concentration of the axial stress in the pile head occur in the piled raft. In seismically active zones, if the piles are structurally connected to the raft, high horizontal shear forces and overturning moments may develop in the pile head due to cyclic or dynamic lateral loads. For example, the foundation system of the RionAntirion Bridge in Greece [5, 6] consists of a gravel bed that limits the shear forces and moments transmitted between the superstructure and the foundation soil. In these cases, the loadcarrying capacity of piles may be governed by their structural capacity rather than by their geotechnical capacity. In addition, constraint reactions between the piles of a raft can be reduced by disconnecting the piles from the raft by inserting a layer of compacted granular material called a cushion [7]. The cushion prohibits loads coming from the superstructure from being directly transmitted to the piles; hence, the disconnected piles function as soil reinforcements rather than as pure structure members [8].
In the latter respect, a disconnected piled raft (DPR) is different from a piled raft (PR) because the piles are disconnected from the raft. Several authors have recently presented the behaviour of a DPR via numerical simulations. Liang et al. [9, 10] showed that a cushion made with sandgravel materials plays an important role in mobilizing the bearing capacity of the subsoil and in modifying the load transfer mechanism of piles. Eslami and Malekshah [11] performed elastoplastic analyses to show that the maximum axial stress along piles occurs at different depths depending on the thickness and stiffness of the cushion. Tradigo et al. [12, 13] performed 3D finite element analyses to show that an increase in cushion thickness will reduce the overall settlement/stiffness.
The behaviour of a DPR has also been investigated by performing model experiments. Nakai et al. [14] performed centrifuge model tests to show that the DPR foundation makes a significant contribution to the dynamic soilstructure interaction. Cao et al. [15] and Wong et al. [16] suggested that piles can be disconnected from a raft to overcome the problem of high stress in the piles and raft and that the gap between the raft and unconnected piles can be filled with a cushion of structural fill material. Fioravante and Giretti [17] and Fioravante [18] performed centrifuge tests on a piled raft model to investigate the load transfer mechanisms between a raft and piles in dry dense sand in cases where the piles are either in direct contact with the raft or separated from the raft by a granular cushion.
In this framework, experimental tests have been performed on piled raft models to explore the influence of a gravel cushion that is interposed between a raft and piles. The insertion of a deformable cushion between a raft and pile heads allows downward soilpile relative displacement, which produces negative skin friction on the upper part of the shaft. Such a mechanism is governed mainly by the stiffness and thickness of the interposed cushion. In the DPR system, piles act as vertical reinforcements of the subsoil and improve its mechanical properties by forming a new composite and stiffer material on which the raft rests [9, 16, 18, 19].
This paper is aimed at investigating the influence of the raftcushionsoilpile interactions on the load transfer mechanisms and cushion failure mode of the piled raft foundation based on experimental tests.
2. Model Tests
A series of model tests are carried out using the test apparatus shown in Figure 1 to study the failure mode of a cushion in a DPR foundation. The test device is mainly composed of a soil box, loading device, and measuring system. The measuring system consists of strain gauges, a load cell, a data acquisition instrument, a linear displacement transducer (LDT), and a digital camera.
2.1. Test Apparatus
Figure 1 presents a schematic diagram of the test apparatus. The soil container size is 1500 mm long, 1000 mm wide, and 1200 mm high. The front of the box is made of tempered glass, while the remaining four sides are made of steel. The model raft is a 300 mm long (), 200 mm wide (), and 50 mm thick () rectangular steel plate characterized by a modulus of elasticity of MPa. The raft can be considered as rigid because the raftsoil stiffness ratio is sufficiently high [20]. Therefore, the settlement of the raft is considered to be uniform in the model test.
2.2. Materials
The experiments were performed using dry fine sand with grains that were 22% finer than 0.075 mm and a relative density of 66%, and the granular cushion was modelled using different types of graded gravel. Figure 3 presents the grainsize distribution of the test soil. The properties of the soil in each layer are given in Table 1.
 
Note. is the effective particle size of soil; is the unit weight of soil; is the internal friction angle; is the cohesion; and is the relative density. 
The piles had various diameters, , of 40 mm, 60 mm, and 80 mm and a length, , of 400 mm. The piles were made from a polyvinyl chloride (PVC) pipe with a compressive modulus of 800 MPa. In addition, full and half piles were used in the tests. In the test interpretation, the piles were considered rigid, and their settlements were assumed to be constant with depth. The external layers of the piles were specially treated via mechanical turning, to ensure that the relative roughness was greater than 0.1, so that the ultimate shear resistance did not depend on pile roughness.
The diameters of the model piles, normalized with respect to particlesize dimension, were greater than 300. This ratio is larger than the lower limit given by Garnier and Konig [21] (1998; ). The main geometrical characteristics of the model and boundary conditions are presented in Figure 1; the ratio between the diameter of the container and the raft side was greater than three, and the distance of the pile tip from the container bottom was always greater than 10 D.
The stresses at the pile head and tip were measured by using two load cells positioned at the top and base of each pile. In addition, the pressure of the subsoil was measured by taking the mean value of four load cells positioned on top of the fine sand, as shown in Figure 2. The pilelongitudinal strains in the central pile as a test pile were measured by using five strain gauges, which were attached to the surface of the pile shaft and covered with epoxy resin. The bondedfoil strain gauges were connected in a Wheatstone bridge and gave an accurate measurement of the axial load.
2.3. Test Cases
A total of 9 tests were conducted to discuss the behaviour of the cushion, as shown in Table 2, in which the main characteristics of the tests are presented, and case 1 was used as a reference test. To find the factors affecting the vertical load transfer mechanism, cases 2 through 7 were focused on the gravel cushion. The behaviour of the cushion was investigated for three types of gravel with different internal friction angles and particle sizes, including graded gravel, coarse gravel, and fine gravel. Case 2, case 3, case 4, and case 5 were used as reference cases to compare the conditions based on the cushion thickness. Case 6 and case 7 were focused on the effect of the cushion stiffness. In this part, case 2, which involved a graded gravel cushion, was compared with case 6 and case 7, which involved coarse gravel and fine gravel. Case 8 and case 9 were focused on the effect of the pile diameter. In this part, case 4, which had a pile diameter of 40 mm, was compared with case 8 and case 9, which had pile diameters of 60 mm and 80 mm, respectively.
 
. is the cushion thickness; is the pile diameter; and is the pile length. 
2.4. Test Procedures
The soil box in the apparatus was filled with fine sand to a height of 1000 mm, and the gravel cushion was laid on the fine sand. The coloured sand was used as a marking line in the soil. The particle size distributions of the soil layers are shown in Figure 3. After completing the gravel cushion, the raft was laid on the cushion, and then a vertical load on the raft was imposed stepwise by using a hydraulic jack. The jacking and loading phases were executed at a constant loading rate with an increment of 400 N, which was measured by the proving ring. For each loading stage, the static load was kept constant for half an hour. Then, the settlement of the raft was obtained by averaging the measurements of the two LDTs that were placed in diametrically opposite positions with respect to the jack, as shown in Figure 1. Meanwhile, the strain on the pile and the soil pressure at the pile head and cushion were collected by the acquisition instrument.
The process of raft sinking was captured using a highresolution camera. Figure 4(a) shows a typical picture of the test, where the analysis grid has been superimposed onto the photography. The images have a resolution of 8.1 pixels per millimetre. Then, using MATLAB’S toolbox for digital image correlation [22], it is possible to track the movement of the soil and to compute the strain field that occurs in the soil mass.
(a)
(b)
Digital image correlation (DIC) is a technique used to calculate displacements from a series of consecutive images. In general terms, the technique involves maximizing a correlation coefficient that is determined from the analysis of a subset of pixels. The analysis grid is specified based on the distances and radii of circular track regions, as shown in Figure 4(b). The optimization process generates a mapping function of marker displacements from one image to the next. The precision of this technique depends on the resolution of the pictures, the analysis grid, and the size of the markers. If the markers are too close, noise is incorporated into the results. On the other hand, if the markers are too far apart, the precision of the estimated displacement is reduced. In this case, the picture was pixels, the centres of markers were 60 pixels apart, and the radius of each marker was 20 pixels. The grid shown in Figure 4(a) has 6000 markers.
3. Test Results and Discussion
The raft load was borne by the pile and fine sand because of the installation of the gravel cushion and then transmitted to the bottom soil layer. During the loading test, each load increment was maintained for approximately half an hour until the deformations of the raft were stable. Therefore, the influence of consolidation on the test results was neglected in the study.
3.1. LoadSettlement Behaviour of DPR
Figure 5(a) presents the loadsettlement behaviour of the disconnected piled raft. The curves of the piled rafts had higher gradients than that of the reference test (case 1) because of the pile being under the cushion. Obviously, the settlements increased as the cushion thickness increased and decreased as the pile diameter increased. Compared to the graded gravel cushion, the settlements increased with fine gravel and decreased with coarse gravel.
(a)
(b)
The settlement performance of a DPR can be quantified by means of the definition of the following settlement efficiency ratio [17]:where and are the settlements of the unpiled raft and the DPR, respectively. Since (and ), the larger the value of , the more efficient the DPR at a given load level.
The dependence of the settlement efficiency on the vertical load is presented in Figure 5(b). The trends in Figure 5(b) show that the settlement efficiency is approximately constant in the experienced range of displacement and that it depends on the thickness and stiffness of the gravel cushion: , 0.61, 0.35, and 0.29 for , 1.0, 1.5, and 2.0, respectively; , 0.69, and 0.16 for the cushions consisting of graded gravel, coarse gravel, and fine gravel, respectively.
3.2. Stress Concentration Ratio
The variations of the pilesoil stress concentration ratio had nearly identical characteristics in all 8 piled raft tests (Figure 6). As shown in Figure 6(a), for different cushion types, the stress concentration ratio versus the vertical load had the same characteristics. With an increase in the vertical load, the stress concentration ratio increased rapidly; it then increased slowly when the vertical load reached 3000 N, except in case 2 and case 6. The larger the ratio, the lower the stress concentration ratio. Case 6 and case 7 showed that the stress concentration ratio decreased with the particle size of the gravel cushion. A lower stress concentration ratio indicates that a higher percentage of the vertical load on the raft is transferred to the soil surrounding the pile through the gravel cushion.
(a)
(b)
Figure 6(b) shows the maximum stress concentration ratio of cases 2–5 and cases 89 versus the ratio; it was observed that the stress concentration ratio decreased as the value of increased. When , the stress concentration ratio decreased gradually and remained within a nearly constant value range of 10–14.
3.3. Axial Forces in Pile
The axial forces in a pile were calculated based on the measured strains of the piles. The deformation of a pile under the loading condition was assumed to be in the elastic state. Thus, the axial force in a pile at a certain depth can be calculated as follows: , where is the average strain of the pile, is the elastic modulus of the pile, and is the crosssectional area of the pile. The elastic modulus and crosssectional area of the pile in our study were 800 MPa and 12.56 cm^{2}, respectively. The average of the measured strains at a certain depth was adopted. There were two strain gauges at certain depths of the pile.
Figure 7 presents the variation of axial forces along the peripheral pile shafts in the tests. The axial forces along the pile shaft were distributed nonlinearly. The distribution patterns of axial forces along the pile shaft were different because of the different cushion thicknesses. As shown in Figure 7(a), the axial force at the top of the pile was the largest, and it decreased with increasing depth of the pile shaft.
(a)
(b)
(c)
(d)
A pile is loaded partially through its head and partially by means of negative skin friction, which acts on its upper shaft. The pilesoil relative displacement is at a maximum at the pile head and tip and is zero at a certain depth with respect to the pile head (neutral plane). As a result, the maximum axial stress of the pile shaft in the downward direction with respect to the pile head was as shown in Figures 7(b)–7(d). The neutral plane was probably at a depth lower than the midpile at the beginning of the load and then moved downwards as the applied load was increased.
3.4. Soil Displacement Nephogram
This test was performed in 20 stages. Figures 89 show the displacements at the end of cases 2–4, cases 67, and case 9 inferred using the DIC technique.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
The horizontal displacement of soil obtained from the DIC technique is shown in Figure 8. To visualize the soil displacement more easily, colour changes are defined in the colour bar. displacement value was the interval for the 6 cases. The unit of the coordinates was mm. The horizontal soil displacement was concentrated at the top of the pile and then developed along both sides, with the pile serving as the axis. In addition, the soil particles were ejected from the pile head.
The common behaviour in the case was a coneshaped area of high density at the pile head, as shown in Figure 9. From the same figure, it was seen that there was a significant change in the density from higher to lower, forming two independent sectors divided by the pile shaft. It was seen that the density at the top of the pile head and that under the raft were greater because the soil was compacted as a result of the settlement of the raft. The density at the edge of sectors decreased because of soil dilatancy as a result of the settlement of the subsoil. The shear plane and dilatancy occurred at the edge of sectors. This resulted in a decrease in the density of soil in the fanshaped area, which was defined by lowdensity areas. The highdensity areas above the pile head represented the soil compaction process.
3.5. Settlement Distribution of Gravel Cushion
The settlement of the gravel cushion was determined based on the soil displacement nephogram, and the displacement curve is shown in Figure 10. The legend of the graph indicates the distance from the bottom of the raft. As seen, the curve shape was similar to a saddle shape, and the settlement of soil on top of the pile head was relatively low, while that beside the pile was much greater, and differential settlement occurred. Figure 10 also shows that the soil displacement outside the raft was negative as a result of the uplift of soil beside the raft.
(a)
(b)
(c)
(d)
(e)
(f)
Figure 10(a) shows the settlement of the gravel cushion for case 2; the pilesoil relative displacement increased as the distance from the bottom of the raft increased. This increase was due to the effect of the compression of the raft and pile. Figures 10(b)–10(f) show the settlement of the gravel cushion for cases 34, cases 67, and case 9. As seen, at a distance of 1.0 cm, no apparent differential settlement occurred. When the distance from the raft was 2.0 cm or more, the pilesoil relative displacement was formed, which increased as the cushion thickness increased.
4. A New Cushion Failure Mode
A fanshaped shear plane occurred at the cushion and a coneshaped area occurred at the pile head according to the results of tests that involved the DIC technique. In addition, the slip plane of the cushion could be determined according to these characteristics. A simple failure mode of the cushion was put forward, as shown in Figure 11.
4.1. Establishment of the Cushion Failure Mode
The cushion failure mode can be divided into three parts. The coneshaped area (area I) at the pile head is the active earth pressure zone. In this area, the cushion moves to the pile top, which is in the passive state, and the pile can be considered as equivalent to a retaining wall. Area II is the transition zone, and the shear zone is formed along the pile top core in this area. The soil of area II slides along the shear zone and compresses the soil in area III. Area III is the passive earth pressure zone. The envelope curve of the three zones forms the slip line, as shown by curves ACDE and ABGF in Figure 11.
The graphical method for determining the cushion failure mode is put forward according to the literature [23]. Draw a right triangle with a proper scale and make the acute angle (∠) equal to the internal friction angle of the gravel. Then, draw a line perpendicular to hypotenuse . Finally, draw three concentric circles from centre with radii of , , and . The smallest, middle, and largest circles are, respectively, plane circle I, vertex circle II, and pole circle III, as shown in Figure 12.
The boundaries of zone I and zone III are determined by using the graphical method. First, draw a chord that is tangent to the plane circle at the point and parallel to the hypotenuse . Afterwards, extend the radius , with the vertex circle and pole circle intersecting at points and , respectively. In addition, retract the radius , with the vertex circle intersecting at point . Line can be analogous to the pile head. Then, the sliding direction of zone I is given by the lines and . In addition, the sliding direction of zone II is given by the lines and .
Line is assumed to be a logarithmic spiral line, as shown in Figure 11. The length of line can be calculated as follows:where is the length of , is the length of , is the angle between and , and is the internal friction angle of the gravel.
Draw the line parallel to to intersect with the gravel surface at the point . The boundary of zone III is given by the lines and . Link points and using a smooth curve, and make line a logarithmic spiral line. Lines and are the radii of the logarithmic spiral line . Therefore, the boundary of zone II is given by lines and .
The shape of the cushion failure mode corresponds to the characteristic circle. In addition, the angles in the failure mode can be calculated based on the characteristic circle; angle is equal to , angle is equal to , and angle is equal to .
4.2. Theoretical Analysis of the Cushion Failure Mode
In developing the present method, the following simplifications were used:(a)The gravel cushion is homogeneous, isotropic, and cohesionless.(b)The gravity of the gravel cushion is ignored.(c)The pile is sufficiently rigid and undergoes insignificant deformation.(d)There is no friction between the pile and the gravel cushion.(e)The gravel cushion undergoes general shear failure.
The cushion failure mode can be simplified to the shape shown in Figure 13 because of the symmetry property of the stress and shape. Area is the active earth pressure zone, while area is the passive earth pressure zone. A formula for determining the stress at the pile head and in the soil between the piles and the pilesoil stress ratio can be derived based on the plastic equilibrium theory.
According to assumptions, the resultant force action point of both the pressure and shear on the logarithmic spiral line is point . Take point as the centre of the moment; then, the moment of the stress on the pile head iswhere is the moment of the stress on the pile head, is the pile diameter, and is the stress on the pile head.
The moment of the active earth pressure iswhere is the moment of the active earth pressure, is the active earth pressure, and is the angle of .
The moment of the passive earth pressure iswhere is the moment of the passive earth pressure and is the active earth pressure.
The moment of the pressure under the raft iswhere is the moment of the pressure under the raft and is the pressure under the raft.
The moment balance equation of point iswhere is the resultant moment of point .
Substituting (3)~(6) into (7), the stress on the pile head is obtained:
The overall equilibrium requires thatwhere is the gravity of the gravel cushion, is the thickness of the cushion, is the earth pressure between piles, and is the area replacement rate of a pile.
Then, the pilesoil stress ratio is given bywhere is the pilesoil stress ratio.
Equations (2)~(10) were deduced under the assumption of general shear failure of the cushion. According to this precondition, the thickness of the cushion should be at least .
5. Comparison Results
5.1. Comparison with Laboratory Experiments
The laboratory experiments described in Section 2 are used to check the validity of the proposed method. The cushion gravel properties and the dimensions of the piles in the experiments are detailed in Tables 1 and 2, respectively. The areas of the pile and raft are 31.4 cm^{2} and 600 cm^{2}, respectively, and the area of the soil between piles is 568.6 cm^{2}, then the value of parameter m is 5.23%.
According to Section 4, the thickness of the cushion should be at least . Therefore, the measured data of case 3, case 4, and case 7 is selected for comparison with the computation results, as shown in Table 3. And the internal friction angles of the gravel of case 3, case 4, and case 7 are 36.7, 36.7, and 33.5 respectively. The value could be calculated by the formula .

When the cushion is filled with graded gravel and coarse gravel, the calculated results of the proposed method and the measured results are very close. The calculated results of case 7 are on the high side because the settlement curve of case 7 undergoes a steep drop, and the cushion failure mode is local shear failure. Therefore, the vertical stress on the pile has the phenomenon of stress focus. In other words, the proposed method is reasonable.
5.2. Comparison with Several Current Design Methods
To verify the proposed method again, the results of another embankment project, calculated using the present method, are also compared with computed values from several existing design methods proposed by Abusharar et al. [24], Low et al. [25], BS8006 [26], and Guido et al. [27]. The geometry of the embankment and the design parameters used in the present case study are obtained from Chen et al. [28] and Abusharar et al. [24]. They are as follows: pile capwidth = 1.13 m; embankment fillheight = 4.52 m; unit weight = 20 kN/m^{3}, and internal friction angle = 30°. A surcharge of 12 kPa is used throughout unless otherwise specified. No partial factors of safety are applied to the design parameters. The results obtained from several current design methods are compared with the results obtained from the present method and summarized in Table 4.

The results show that the methods of Low et al. [25] and BS8006 [26] overpredict the vertical stress on the pile compared to the method of Guido et al. [27] and the present method because the methods of Low et al. [25] and BS8006 [26] consider the reaction of the geosynthetics. Because the present method and the methods of BS8006 [26] and Guido et al. [27] do not consider the support of the subground by the geosynthetics, their predictions of the stress concentration ratio are close. It can be noted that the prediction results of the vertical stress on the pile of the present method are in good agreement with the results of the methods of Guido et al. [27] and Abusharar et al. [24].
6. Conclusions
In this paper, an experimental and theoretical study of the cushion failure mode of a rigid pile foundation was presented based on the results of eight laboratory tests. The focus was on the thickness and gradation of a gravel cushion and the pile diameter.
The experimental procedure consisted of reproducing the cushion failure mode, with a highresolution camera used to track the experience. The displacement field in the soil was inferred using the digital image correlation (DIC) technique, which tracks the pixels in images and generates a displacement field. Moreover, to discuss the load transfer mechanisms of a composite foundation, five strain gauges attached to the pile, five earth pressure cells, and two displacement transducers were instrumented in the model tests. The main conclusions of this study are as follows:
The settlements of the raft increased with the cushion thickness and decreased with the pile diameter. In addition, the pilesoil stress concentration ratio decreased with the ratio and with the particle size of the gravel cushion.
Negative friction resistance of the pile shaft occurred in cases 3–5 because the gravel cushion at the top of each pile had a sufficient thickness. The neutral plane of the pile shaft moved down as the vertical load increased. The neutral plane was probably at a depth lower than the midpile at the beginning of the load and then moved downwards as the applied load was increased.
The soil displacement field was reproduced by using the DIC technique. It was seen that a fanshaped shear plane occurred at the cushion and that a coneshaped area occurred at the pile head. In addition, a shear plane occurred at the edge of sectors because the density at the edge of sectors was decreased due to soil dilatancy as a result of the settlement of the subsoil.
A simple cushion failure mode was put forward according to the slip plane of the cushion. Then, a new theoretical analysis was presented for a rigid pile composite foundation. The results of the proposed methods were verified by laboratory tests results.
A comparison of the amount of vertical stress on a pile predicted by the present method to those of the current design methods showed good agreement with the results of Guido et al. [27] and Abusharar et al. [24]. As for the presented method, the values of the stress concentration ratio predicted by the methods of BS8006 [26] and Guido et al. [27] were close because these methods do not consider the support of the subground by the geosynthetics.
In this study, all the laboratory tests were conducted using rigid piles and a gravel cushion. In future research, experiments with various cushion soils and foundation soils and with flexible piles should be studied to reproduce the cushion failure mode. In addition, further research involving a fullscale test or a centrifuge model test is required and would be useful for accurately investigating the validity of the present method.
Conflicts of Interest
The author declares that he has no conflicts of interest.
Acknowledgments
This study was supported by the National Program on Key Basic Research Project of China (Project no. 2013CB036304), Nature Science Foundation of Jiangsu Province of China (Project no. BK20170509), and Natural Science Foundation for Colleges and Universities in Jiangsu Province (Project no. 15KJB580013).
References
 J. Burland, B. Broms, and V. De Mello, “Behaviour of foundations and structures,” in Proceedings of the 9th International Conference on Soil Mechanics and Foundation Engineering (ICSMFE '77), vol. 2, pp. 495–546, Tokyo, Japan, 1977. View at: Google Scholar
 H. G. Poulos, “Methods of analysis of piled raft foundations—a report prepared on behalf of Technical Committee TC18 on piled foundations,” Tech. Rep., International Society of Soil Mechanics and Geotechnical Engineering, London, UK, 2001. View at: Google Scholar
 M. Randolph, “Design methods for pile groups and piled rafts,” in State of the Art Report: 13th International Conference of Soil Mechanics and Foundation Engineering, New Delhi, 1994. View at: Google Scholar
 M. Huang, F. Liang, and Z. Li, “Recent advances in the analysis of pile foundation in China,” in Advances in Deep Foundations, Y. Kikuchi, J. Otani, M. Kimura, and Y. Morikawa, Eds., pp. 115–124, Taylor & Francis, London, UK, 2007. View at: Google Scholar
 J. Garnier and A. Pecker, “Use of centrifuge tests for the validation of innovative concepts in foundation engineering,” in Proceedings of the 2nd International Conference of Earthquake Geotechnical Engineering, pp. 431–439, 1999. View at: Google Scholar
 A. Pecker, “Design and construction of the rion antirion bridge,” in Geotrans, pp. 216–240, 2004. View at: Google Scholar
 N. Mattsson, A. Menoret, C. Simon, and M. Ray, “Case study of a fullscale load test of a piled raft with an interposed layer for a nuclear storage facility,” Géotechnique, vol. 63, no. 11, pp. 965–976, 2013. View at: Publisher Site  Google Scholar
 H. G. Poulos, “Piled raft foundations: design and applications,” Géotechnique, vol. 51, no. 2, pp. 111–113, 2001. View at: Publisher Site  Google Scholar
 F.Y. Liang, L.Z. Chen, and X.G. Shi, “Numerical analysis of composite piled raft with cushion subjected to vertical load,” Computers and Geotechnics, vol. 30, no. 6, pp. 443–453, 2003. View at: Publisher Site  Google Scholar
 F. Liang, J. Li, and L. Chen, “Optimization of composite piled raft foundation with varied rigidity of cushion,” in Foundation Analysis and Design: Innovative Methods, pp. 29–34, 2006. View at: Publisher Site  Google Scholar
 A. Eslami and S. S. Malekshah, “Analysis of nonconnected piled raft foundations (NCPRF) with cushion by finite element method,” Computational Methods in Civil Engineering, vol. 2, no. 2, 2011. View at: Google Scholar
 F. Tradigo, F. Pisanò, C. Di Prisco, and A. Mussi, “Nonlinear soilstructure interaction in disconnected piled raft foundations,” Computers and Geotechnics, vol. 63, pp. 121–134, 2014. View at: Publisher Site  Google Scholar
 F. Tradigo, F. Pisanò, and C. di Prisco, “On the use of embedded pile elements for the numerical analysis of disconnected piled rafts,” Computers and Geotechnics, vol. 72, pp. 89–99, 2016. View at: Publisher Site  Google Scholar
 S. Nakaia, H. Katoa, R. Ishidaa et al., “Load bearing mechanism of piled raft foundation during earthquake,” in Proceedings of the of 3rd UJNR Workshop on SoilStructure Interaction, pp. 2930, March 2004. View at: Google Scholar
 X. D. Cao, I. H. Wong, and M.F. Chang, “Behavior of model rafts resting on pilereinforced sand,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 130, no. 2, pp. 129–138, 2004. View at: Publisher Site  Google Scholar
 I. H. Wong, M. F. Chang, and X. D. Cao, “Raft foundations with disconnected settlementreducing piles,” in Design Applications of Raft Foundations, J. A. Hensley, Ed., chapter 17, Thomas Telford Ltd, 2000. View at: Google Scholar
 V. Fioravante and D. Giretti, “Contact versus noncontact piled raft foundations,” Canadian Geotechnical Journal, vol. 47, no. 11, pp. 1271–1287, 2010. View at: Publisher Site  Google Scholar
 V. Fioravante, “Load transfer from a raft to a pile with an interposed layer,” Géotechnique, vol. 61, no. 2, pp. 121–132, 2011. View at: Publisher Site  Google Scholar
 M. Huang, F. Liang, and Z. Li, “Recent advances in the analysis of pile foundation in China,” Doklady Akademii Nauk SSSR, vol. 302, no. 2, pp. 480–483, 2007. View at: Google Scholar
 K. Horikoshi and M. F. Randolph, “On the definition of raft—soil stiffness ratio for rectangular rafts,” Géotechnique, vol. 47, no. 5, pp. 1055–1061, 1997. View at: Publisher Site  Google Scholar
 J. Garnier and D. Konig, “Scale effects in piles and nails loading tests in sand,” in Centrifuge '98, Tokyo, Japan, 23–25 September 1998, T. Kimura, O. Kusakabe, and J. Takemura, Eds., pp. 205–210, Balkema, Rotterdam, The Netherlands, 1998. View at: Google Scholar
 E. Christoph, T. Robert, and G. Daniel, Digital Image Correlation and Tracking with Matlab, 2006.
 T. K. Krein, Structural Mechanics of Granular Media, translated by C. Wanjia, China Railway Publishing House, Beijing, China, 1983.
 S. W. Abusharar, J. J. Zheng, B. G. Chen et al., “A simplified method for analysis of a piled embankment reinforced with geosynthetics,” Geotextiles & Geomembranes, vol. 27, no. 1, pp. 39–52, 2009. View at: Google Scholar
 B. K. Low, S. K. Tang, and V. Choa, “Arching in piled embankments,” Journal of Geotechnical Engineering, vol. 120, no. 11, pp. 1917–1938, 1994. View at: Publisher Site  Google Scholar
 British Standards Institution, “Code of practice for strengthened, reinforced soils and other fills,” British Standard BS 8006:1995, 1995. View at: Google Scholar
 V. A. Guido, J. D. Kneuppel, and M. A. Sweeny, “Plate loading tests on geogridreinforced earth slabs,” in Proceedings of Geosynthetics '87, New Orleans, USA, pp. 216–225, IFAI, 1987. View at: Google Scholar
 R. P. C. P. Chen, Y. M. C. M. Chen, J. H. Han, and Z. Z. Xu, “A theoretical solution for pilesupported embankments on soft soils under onedimensional compression,” Canadian Geotechnical Journal, vol. 45, no. 5, pp. 611–623, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Xiaojun Zhu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.