M20
M20
M20
Goa, India
th
1 Introduction
For modelling piled foundation, the Plaxis 3D Foundation Program provides embedded pile model, in which the pile is assumed as slender beam elements which are virtually connected to the soil by means of the skin and foot interfaces. Since these elements may have arbitrary inclination and cross the soil elements at any arbitrary position, special-purposed interface elements. The interaction between the pile and the soil at the skin is modelled by means of line-to-volume interface elements and the interaction at the base by means of point-to volume interface elements (Septanika 2005a and 2005b), i.e. in addition to the embedded beam approach as developed by Sadek and Shahrour (2004). Further, the current embedded pile approach considers: (i) different types of skin traction/slippage model (constant/linear, multi-linear and layer-dependent), and (ii) foot slippage model. In the foot slippage model, the foot resistance corresponds to the maximum force that can be sustained by the foot interface during compression. For single pile, validation studies for both compression case and tensile tests case have previously been performed in which the numerical results according to embedded pile are compared with the field test data (Engin et al. 2007, Septanika et al. 2007). To make this paper self-contained, the background of the embedded pile model is shortly recapitulated. Next, the results of the single pile studies are also shortly discussed. As an extension to the previous studies, this paper considers pile group behaviour by applying embedded piles onto simplified cases. First, the results of embedded pile for the idealized problem as described in e.g. Poulos (2001) are compared to various methods. Then a qualitative study of the piled raft foundation is presented, considering the group effect in case the cap is detached from the soil and in case the cap touches the soil. In case the cap detaches the soil, the spacing effect is first examined, excluding the cap contribution. In case the cap touches the soil, complex interaction between the raft-pile-soil will be involved. Besides the soil-pile interaction and the pile spacing, the raft-to-soil interaction and the interaction between the cap and the pile groups will also be involved. Finally, the application of embedded piles as soil reinforcements is also presented. This is only to illustrate the range of applications into which the current embedded pile model can be applied. (Note that soil reinforcements are modelled by simply excluding the foot resistance, while in standard piled foundation applications both the skin and the foot resistances are considered).
3231
2 Background
skin
with
rel
= (
(1)
where t is the traction increments at the integration points, T is the material stiffness matrix of the pseudo skin interface and urel = (up us) represents the relative displacement vector between the soil and the pile. The element stiffness matrix Kskin representing the pile-soil interaction at the mantle has been derived based on the following internal virtual work consideration (Septanika 2005a). To include slippage at the pile-soil contact in 3D Foundation Program, one may use the skin-traction model according to: (a) constant/linear traction-depth model, (b) multi-linear model, and (c) layer-dependent model which relates the allowable traction to the strength of the adjacent soil layer.
skin
The interaction at the foot is modeled by a special-purposed spring element to represent the foot stiffness against the relative movements at the foot. This spring connects the pile foot to the soil in the vicinity of the foot. The force acting on this spring Ffoot vector is determined as
F foot = D foot
foot rel
with
foot rel
= (
(2)
where Ffoot is the force increment, Dfoot represents the material stiffness matrix of the pseudo spring element at
foot the foot, rel represents the relative displacement vector between the soil and the pile at the foot. For the maximum foot resistance representing the base failure (due to penetration or pulled-out) at the pile foot, the following simplified criterion has been utilized
3232
foot foot foot (compression) and Faxial Faxial Fmax = 0 (tension) foot Faxial
(3)
is the axial component of the force at the pile foot. In case of reinforcement the base resistance is where set to zero. Further, the mesh-dependent behaviour has been improved by the application of an elastic zone approach for the soil region within the pile. This approach appears sufficient for reducing the undesirable meshdependent effects (Engin et al., 2007).
400
Load (kN)
300
200
100
Pile Capacity - 5.3 m long pile 5.3 m Test Pile Embedded Pile - 5.3 m long pile Pile Capacity - 3.3 m long pile 3.3 m Test Pile Embedded Pile - 3.3 m long pile
0 10 20 30 40
Settlement (mm)
(a)
(b)
Figure 2. (a) Load-displacement curve of the Alzey bridge pile load test together with embedded pile results and (After Engin et al., 2007) (b) Load-displacement curve of the South Surra pile load test together with embedded pile results. (Engin, 2007)
4 Pile Group
3233
considered as schematized in Figure 3a. The raft foundation is supported by nine piles. Bearing capacity of raft is 0.3 MPa, load capacity of each pile is 0.873 MN (compression) and 0.786 MN (tension). The pile diameter is 0.5 m and the length is 10 m. The pile and raft properties are Ep = Er = 30 GPa and p = r = 0.2.
(a)
(b)
Figure 3. (a) Idealized problem and (b) FE-model in PLAXIS 3D Foundation Program. The results of analyses in 3D Foundation Program (Figure 3b) have been compared to the results from six different methods as reported in Poulos (2001), according to: Poulos and Davis (1980) Randolph (1994) Strip on spring analysis using GASP (Poulos, 1991) Plate on spring approach using GARP (Poulos, 1994) Finite element and boundary element method of Ta and Small (1996) Finite element and boundary element method of Sinha (1996) Figure 4 show the comparison results of the computed characteristics of behaviour of a raft supported by nine piles. The applied load of 12 MN exceeds the ultimate capacity of the piles alone, and hence the global behaviour will be highly non-linear. Despite some differences between various methods, the computed characteristic results are somewhat similar
Moment (MNm/ m)
1,2 1,0 0,8 0,6
Plaxis3D Fnd
R andolph
15,0
(a)
(b)
Figure 4. (a) Average settlements and (b) Bending moments according to various methods.
3234
S trip (GAS P)
FE+ BE S inha
Plate (GAS P)
20,0
FE Ta & S mall
FE+ BE S inha
detached from the soil), the total capacity for all spacing is according to the capacity of a single pile times the number of piles. It is to be noted that for s/d=2 a very large settlement at failure is found. This is probably due to overlapping of the pile influence regions (causing stress bulbs). The soil between the piles may be locked, and hence, the skin capacity is mainly be contributed by the corner piles. However, the traction profiles at failure for different spacing appear uniform as shown in Figure 8. This indicates that for the present case the skin traction is fully mobilized at each pile when the raft foundation fails.
P > Pult
10 m
20 m soil = pile = 20 kN/m3 c = 50 kPa = 0 D=0.5 m (pile diameter) Fskin = c..D = ~ 78.5 kN/m assumed single pile skin capacity Ftip = 9.cu.. D2/4 = ~ 88.4 kN ..assumed single pile tip resistance
30 m
s/d=2.0
s/d=4.0
s/d=6.0
s/d=8.0
Raft only
Figure 6. Patterns used in the analyses for the investigation of pile spacing effect.
3235
Load (kN)
1.00
2.00
3.00
4.00
Settlement (cm)
Figure 8. Skin traction (shear force per unit length) profiles at embedded piles for various spacing.
3236
50000
Cap Thickness, d = 0.75 m
45000 40000
50000
Allowable Pile Group Capacities
35000
40000
9 x Single Pile - Maximum Allowable Capacity Footing without pile - Minimum Allowable Capacity
Load (kN)
30000 25000
s/d=2.0
33196 30000
s/d=4.0 s/d=6.0 s/d=8.0 Single Pile with Small Cap 9 x SP Big Cap without piles 9 x SP without cap
20000
19983
10000
20704
21192
21847
16903
10
10 15 Settlement (cm)
20
25
s/d Ratio
(a)
(b)
Figure 9. (a) Force-displacement curves for different spacing and (b) Allowable loads in case of allowable settlement is 5 cm.
(a)
(b)
Figure 10. (a) FE-model of reinforced slope and (b) Computed safety factor based on the phi-c method.
3237
5 Concluding Remarks
This paper shortly describes the performance of the current embedded pile model available in 3D Foundation Program. The accuracy of single pile model has been validated, by considering the pile compression tests in Frankfurt and the pile tension tests in Kuwait. For both cases, the results are reasonably in agreement with the field test results. The numerical performance of the current embedded pile implementation has been verified by considering the idealized problem and comparing the results to the various methods. Further, to verify the accuracy of the skin and foot interaction models in case of pile groups, the piled raft has been considered in which the cap is detached from the soil. Full mobilization of the skin tractions and foot resistance illustrates the accuracy of the present interaction models. In the piled raft case with the cap touching the soil, the embedded pile group appears capable of catching the spacing effects on the resulting allowable capacity. For all cases the piled raft capacities are higher that the raft capacity alone. For this particular case, the allowable capacity seems to increase with spacing and reaching the upper limit at spacing nearly eight times pile diameter. Finally, to illustrate the range of applicability of the current embedded pile model, an idealized reinforced slope problem is also presented.
6 Acknowledgements
The authors are grateful to Dr. Paul Bonnier (Plaxis BV) for his supports during the code implementation.
7 References
El-Mossallamy Y. and Franke E. 1997. Numerical Modelling to Simulate the Behaviour of Piled Raft Foundations, August 1997, Darmstadt (Germany). El-Mossallamy Y. 1999. Load-settlement behaviour of large diameter bored piles in over-consolidated clay. Proc. 7th Int. Symp.on Numerical Models in Geotechnical Engineering - NUMOG VII, Graz, 1-3 September 1999, 443-450. Balkema. Rotterdam (The Netherlands). Engin H.K., Septanika E.G. and Brinkgreve R.B.J. 2007. Improved embedded beam elements for the modelling of piles. Proc. 10th Int. Symp. on Numerical Models in Geotechnical Engineering NUMOG X, Rhodes (Greece). April 2007. Engin, H.K. 2007. Report on tension pile testing using embedded piles. Plaxis internal report. Delft (The Netherlands) Ismael N.F., Al-Sanad H.A. and Al-Otaibi F. 1994. Tension tests on bored piles in cemented desert sands. Canadian Geotechnical Journal, 31(4), 597-603. Poulos H.G. and Davis, E.H. 1980. Pile foundation analysis and design. New York (USA). Wiley. Poulos H.G. 1991. In Computer Methods and Advances in Geomechanics (eds Beer et al.). 183-191. Balkema. Rotterdam (The Netherlands). Poulos H.G. 1993. Piled rafts in swelling or consolidating soils. J. Geotechical Div., ASCE, 119(2), 374-380. Poulos H.G. 1994. An approximate numerical analysis of pole-raft interaction. Int. J. NAM Geomech. 18, 73-92. Poulos H.G. 2001. Piled raft foundations: design and applications. Geotechnique, 51(2), 95-113. Randolph M.F. 1994. Design method for pile groups and piled rafts: state-of-the-art report. Proc. 13th Int. Conf. Soil Mech.Found. Engng, New Delhi (India), 5, 61-82. Sadek M. and Shahrour I. 2004. A three dimensional embedded beam element for reinforced geomaterials. Int. J. for Numerical and Analytical Methods in Geomechanics 28, 931946. Sinha J. 1996. Analysis of piles and piled rafts in swelling and shrinking soils. PhD Thesis, Univ. of Sydney, Australia. Septanika E. G. 2005a. A finite element description of the embedded pile model. Plaxis internal report. Delft (The Netherlands). Septanika E. G. 2005b. Validation testing embedded pile in Plaxis 3D Foundation. Plaxis internal report. Delft (The Netherlands). Septanika E.G., Bonnier P.G., Brinkgreve R.B.J. and Bakker, K.J. 2007a. An efficient 3D modelling of (multi) pile-soil interaction. Proc. of 3rd Int. Geomechanics Conference, 67-76, Nessebar (Bulgaria). Septanika E.G., Bonnier P.G., Bakker, K.J. and Brinkgreve R.B.J. 2007b. 3D soil reinforcement modelling by means of embedded pile. Proc. of 3rd Int. Conference of IS-Kyushu, November 2007, Fukuoka (Japan). Skempton A.W. 1951. The Bearing Capacity of Clays. Proceeding of Building Research Congress, 1, 180-189. Ta L.D. and Small J.C. 1996. Analysis of piled raft systems in layered soils. Int. J. NAM Geomech. 2, 57-72.
3238