Journal of Applied Mathematics

Volume 2013 (2013), Article ID 306849, 10 pages

http://dx.doi.org/10.1155/2013/306849

Research Article

Analytical Solutions of Spherical Cavity Expansion Near a Slope due to Pile Installation

Jingpei Li,1,2 Yaguo Zhang,1,2 Haibing Chen,1,2 and Fayun Liang1,2

¹Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Tongji University, Shanghai 200092, China

²Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China

Received 5 June 2013; Accepted 18 August 2013

Academic Editor: Ga Zhang

Copyright © 2013 Jingpei Li et al. 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.

Abstract.

Based on the hypothesis that the penetration of a single pile can be simulated by a series of spherical cavity expansions, this paper presents an analytical solution of cavity expansion near the sloping ground. Compared with the cavity expansion in the half-space, the sloping free boundary has been taken into account as well as the horizontal free boundary. The sloping and horizontal free surfaces are considered by the introduction of a virtual image technique, the harmonic function, and the Boussinesq solution. The results show that the sloping free boundary and the variation of the inclination angle have pronounced influences on the distribution of the stress and displacement induced by the spherical cavity expansion. The present solution provides a simplified and realistic theoretical method to predict the soil behaviors around the spherical cavity near the sloping ground. The approach can also be used for the determination of the inclination angle of the slope according to the maximum permissible displacement.

1. Introduction

There are many situations where foundations need to be located on the top of a slope, such as the piled bridge abutment adjacent to a slope crest. Hence, pile installation in sloping ground has attracted wide concerns [1–3]. In contrast to the horizontal ground surface case, the boundary effects of a slope should be considered for situations of piles embedded adjacent to the slope. The existing boundary not only affects the bearing capacity of piles, but also adds the risk of slope failure [4]. For instance, a riverbank dike, located along the Bailianjing River in Shanghai, was damaged by pile driving in soft clay during the construction of a newly elevated dike [5].

Following the early suggestion [6, 7], solutions of the limit pressures of spherical and cylindrical cavities are used to predict the end bearing and shaft capacities of piles [8, 9], as well as the stress fields and lateral displacements of the surrounding subsoil induced by installation of a pile [10–13]. However, the solutions of a cavity expansion in an infinite medium do not satisfy the stress conditions at the free surface during the pile installation. Sagaseta et al. [14] and Sagaseta [15] considered the problem as strain controlled and obtained strains by using only the incompressibility condition. The presence of the top free surface was considered by means of a virtual image technique and some results for the elastic half-space. Besides, Keer et al. [16] derived a solution for the expansion of spherical cavity in a half-space by using the image source method [17] and the concept of cavity expansion source [18]. These methods can be well used to analyze the boundary effects of the free surface of the half-space, but they are not directly applicable to cavity expansion near a slope. Compared with cavity expansion in a semi-infinite half-space, the slope surface should be taken into account.

In this paper, the expansion caused by pile tip is simulated as a spherical cavity expansion. Theoretical solutions for the expansion of a single spherical cavity near slope are derived by using the virtual image approach. Meanwhile, the correction stress functions and the Boussinesq solutions are introduced to consider the effects of both horizontal ground surface and slope surface in this analysis.

2. Basic Theories and Geometry of the Problem

The concept of the cavity expansion source was first used by Hopkins [18]. The model shown in Figure 1 is a cavity under a uniform pressure on the internal spherical surface with radius .

Figure 1: Cavity expansion source induced by the uniform pressure on the internal spherical surface.

The fields of stress and displacement induced by the pressure are as follows:

where the is the distance between the center and the calculation point, the is the shear modulus of the soil, and the , , and refer to radial, hoop, and tangential stresses, respectively. Due to the spherical symmetry, the shear stresses are equal to zero, that is, .

The expansion caused by the pile tip is usually simulated as a spherical cavity expansion [19, 20]. Hence, attention will be paid to the solutions of a single cavity expansion adjacent to a slope. As shown in Figure 2, the depth of the cavity below the horizontal ground surface is denoted by , and the angle between the slope and vertical direction is denoted by β. Unlike cavity expansion in an infinite medium, the analytical solutions of cavity problem near a slope are currently only possible in elastic materials. Accordingly, the soil is assumed to be an isotropic, homogeneous, and linear elastic material, and only small strains occur during the process of the cavity expansion. For simplicity, the gravitational stresses are ignored.

Figure 2: Expansion of a spherical cavity near the sloping ground.

3. The Solution Method

3.1. The Expansion of a Spherical Cavity and Its Image in the Half-Space

The theoretical solution for the expansion of a single spherical cavity in a half-space has been presented by Keer et al. [16], where the free surface of the half-space is horizontal (). Similarly, the solutions can also be derived when the free surface is vertical (). Taking the vertical surface as the plane of symmetry, another virtual spherical cavity is put at the image point which is shown in Figure 3. The coordinate of the calculation point is , and is the horizontal distance between the center of spherical cavity and the vertical free surface. is the distance from the spherical cavity to point , and is the distance from the image to point . and are the angles from direction to and , respectively.

Figure 3: A cavity expansion source and its image.

The stress and displacement components of the cavity expansion in a cylindrical coordinate system can be written as

By substituting (1) to (3) into (4) to (7) and using the principle of superposition, the stress and displacement of the spherical cavity and its image in the cylindrical coordinate are as follows:

where and .

The horizontal ground and the slope are both free surfaces, on which the normal and shear stresses are zero. Because of the symmetry, the source and its image produce zero shear stress () and a nonzero normal stress () on the free surface . Thus, the presence of the normal stress violates the free surface boundary condition

where .

3.2. The Expansion of a Spherical Cavity and Its Image Near a Slope

The boundary effects of both ground surface and slope should be investigated during the analysis of the cavity expansion near the sloping ground. The problem turns to be more difficult due to the increase of the slope boundary. Similarly, the virtual image technique is employed to consider the boundary effects of the slope. Taking the slope as the plane of symmetry, the image source is set at the image point as shown in Figure 4.

(a)

(b)

(c)

Figure 4: The expansion of the single spherical cavity near a slope. (a) The image source is located above the ground surface. (b) The image source is located on the ground surface. (c) The image source is located below the ground surface.

The solutions of the two cavities expansion (the actual spherical cavity and its virtual image) can be obtained by the principle of superposition. As a result, the stress and displacement components in the cylindrical coordinate system are shown to be

where the cavity depth below the ground surface is denoted by , the distance from the source to the calculation point is , and the image to the point is . According to the relationship of geometry shown in Figure 4, the expressions for and are

With the increasing of depth of the spherical cavity, the image cavity gradually moves from above ground (Figure 4(a)) to the ground below (Figure 4(c)). Accordingly, the expressions of and will be changed with the depth of cavities. Specifically, when the virtual image cavity is just on the ground surface as shown in Figure 4(b), the can be simplified to be

where and are defined as the horizontal ( direction) and vertical ( direction) distances between the source and its image, respectively.

3.3. The Correction of the Stresses on the Horizontal Ground

According to the models shown in Figure 4, the actual expansion cavity and its image do produce not only nonzero normal stress but also shear stress on the horizontal ground surface , which can be shown to be as follows:

where and .

In order to satisfy the condition of the free horizontal boundary (), the different correction functions are introduced to deal with the normal stress (see (25)) and shear stress (see (26)) on the boundary. Based on the theory of superposition, the stresses on the ground surface can be divided into two parts:

(i)only normal stress on the horizontal ground surface: , ;(ii)only shear stress on the horizontal ground surface: , .

Using the axially symmetric stress function of Kassir and Sih [21] in the first part (i), the corresponding stress and displacement solutions are written as follows:

in which

where is a constant. The horizontal ground surface is treated as the plane of symmetry. denotes the distance between the point and the symmetrical position of the actual cavity . Likewise, is the distance from the point to the symmetrical position of the image cavity . As a result, their expressions can be written as follows:

is a harmonic function, that is, . The corresponding stresses on the free surface are:

where , and .

Thus, the normal stress on the horizontal ground surface could be eliminated by the stress function with accurate value of . According to the following equation:

it can be obtained that

Substituting the above expression of back into (33), one has

and then the stress and displacement components can be obtained by substituting (39) into (27) to (32).

For the second part (ii), there is only shear stress on the horizontal ground surface. In the same way, the stress harmonic function is given to eliminate the shear on the boundary, and the solutions of the stress and displacement are as follows:

Similarly, the expression of stress function can be written as follows:

where is a constant. The stress function satisfies the equilibrium , and the corresponding stress components on the ground surface are

Accordingly, it can be obtained that .

The final expression for the stress function is then as follows:

Substituting (48) into (40) to (45), the stress and displacement components can be obtained.

3.4. The Correction of the Stresses on the Slope Surface

Considering a virtual source, positive mirror image of the actual cavity with respect to a slope surface will produce the same normal stresses and opposite shear stresses as the actual cavity, the shear stress is eliminated, and the normal stress increases doubly, as shown in Figure 5.

Figure 5: The stress distribution of the slope surface.

Using the method of coordinate transformation, the normal and shear stresses can be obtained as follows:

where ,  , and .

In order to satisfy the free surface boundary condition as much as possible, the Boussinesq solution has been introduced to correct the normal stress. Stress is applied on the surface of the slope, which is equal to the normal stress in value but opposite in direction.

As shown in Figure 6, is the intersection line of the slope and the horizontal plane. A small element with an area of is taken out of the slope lor′ (i.e., or′ in Figure 5) for analysis. Further, the force exerted on the small element is equal to q′ρdθdρ. Using the Boussinesq solutions, the stress and displacement components of soil under the action of the force (q′ρdθdρ) can be derived as follows:

in which

Figure 6: The stress analysis of calculation point by the force on the small element.

The directions of stress and displacement induced by the force on the small element are not always in accordance with coordinate axes (Figure 6). By means of coordinate transformation, the solutions of stresses and displacements in the coordinate system r′oz′ (Figure 5) can be obtained. After that, integrating results with respect to the whole sloping ground surface can lead to the following:

Hence, the stress and displacement components in the coordinate system roz can be derived as follows:

where , , and + .

With (53), the stress and displacement fields induced by the stresses q′ on the slope can be derived. Using a virtual source of the actual cavity at the image point, the lateral deformations of soil around spherical cavity were predicted by Rao et al. [22]. Meanwhile, the Cerruti solutions are used to eliminate the shear stresses produced by the expansion of both the real and the imaginary spherical cavities in an infinite space. However, their results are inappropriate because of the incorrect using of the Cerruti solutions. It is known that a distributed stress is clearly not a point force in the elementary sense. Hence, the stress should be integrated with respect to the slope surface when it is substituted into the Cerruti solutions.

Accordingly, when considering the effects of both horizontal and sloping free boundaries, the final results of the expansion of a single spherical cavity near a slope (Figure 2) can be obtained by superposition of all the parts stresses and displacements:

4. Discussion of the Solutions

The presence of the horizontal and sloping free surfaces is considered in this paper. Consequently, the present solutions have more extensive applications compared with solutions of Keer et al. [16]. According to Section 3, the solutions for this problem can be derived in four steps (cavity expansion in an infinite medium, cavity and its image expansion in an infinite medium, and the corrections of stresses on horizontal surface and sloping surface). The range of the angle β between the slope and the vertical plane is .

When , the slope turns to be vertical plane, as shown in Figure 3. Meanwhile, it can be derived from (23) and (24) that the horizontal distance between the actual cavity and its image is , while the vertical distance is . As a result, (15) to (20) convert into (8) to (13) in sequence.

When , there is only a free surface in the horizontal direction. In this case, substituting β into (23) and (24), and can be obtained. Accordingly, (15) converts into equations that were proposed by Keer et al. [16].

When , there are horizontal and sloping free boundaries as described in this paper. Therefore, the solutions of expansion of spherical cavity in half-space with horizontal and vertical free boundaries or only a horizontal free boundary are the particular cases of present solutions.

As stated above, (15) to (20) can degenerate to the existing solutions for the extreme cases of horizontal ground and vertical slope, which demonstrates the correctness of solutions derived by the first two steps. Steps 3 and 4 involve the correction of stresses on free surface, which is based on an understanding that existing stresses can be offset by the stresses with the same magnitude and opposite direction. In fact, effects of Steps 3 and 4 are further demonstrated by the analysis cases below.

5. Results and Parameters Analysis

In order to consider all the stress components together, the variation of the Mises stress during the expansion of the cavity is analyzed. The parameters in the example analyses are  kPa,  kPa, .25 m, and (incompressible undrained clay). As the distance increases or the sloping free boundary is approached, the Mises stress decreases very rapidly in the range , but then it decreases it decreases more slowly with further increase of the distance , as shown in Figure 7. In order to give a further discussion of the influence from the inclination angle of the slope, the analysis covers with different angles: , , , , and . The corresponding Mises stresses are shown in Figure 7, respectively. Evidently, the angle β has a pronounced influence on the distributions of the Mises stress. At the same point near the spherical cavity, the Mises stress increases with the increase of the angle β. For example, the Mises stress rises from  kPa to  kPa, again at the instant of , while the β increases from 15° to 75°. This is consistent with the results for spherical cavity expansion in half-space reported by Keer et al. [16]. The function in Figure 7 represented the Mises stress: . Here, ( is the second principle invariance).

Figure 7: The Mises stress varies with the distance .

If the free surfaces are not taken into account, normal and shear stresses will appear at the place where there ought to be the free surface, so violating the imposed boundary condition of a free surface [13, 15, 16, 23]. As described in Section 3.4, the shear stress on the sloping ground can be eliminated by using the virtual image technique. In order to cancel the normal stress on sloping surface as much as possible, the Boussinesq solutions are introduced. Another group of normal and shear stresses is produced again when the Boussinesq solutions are applied to correct stresses on the sloping surface, although the original stresses on horizontal free surface have been removed by introduction of harmonic functions (see (39) and (48)). That is because the Boussinesq solutions are aimed at half-space problems.

Figures 8(a) and 8(b) show the normal and shear stresses that vary along the horizontal free surface, respectively, where is the distance from origin of the coordinate on the horizontal free surface to the opposite direction of the -axis. The stresses induced by the cavity expansion in an infinite space show large variations, which are determined by the distance from the center to the calculation point. Compared with solutions in an infinite space, both normal and shear stresses on the horizontal free surface corrected by the introduction of the harmonic functions are in close proximity to zero (Figure 8), which demonstrates the validity of the present method. For example, the maximum of the normal stress calculated by present method falls from −0.95 kPa to −0.029 kPa, and the shear stress falls from 0.14 kPa to 0.047 kPa in case that . With the increase of the distance , the normal and shear stresses on the horizontal free surface decrease gradually. This is because the influence from the stresses on sloping surface weakens when the distance increases.

(a)   Normal stress

(b) Shear stress

Figure 8: Stress along the horizontal free surface.

For the sloping free surface, although stresses have been offset by means of virtual image technique, both normal and shear stresses on the sloping surface are produced again when the harmonic functions (see (39) and (48)) are used to correct the stresses on horizontal free surface. Figures 9(a) and 9(b) show, respectively, the variations of the normal and shear stresses vary along the sloping free surface in cases and , where is the distance from origin of the coordinate to an arbitrary point on the sloping free surface and is the distance between the slope and the center of the spherical cavity. Compared to solutions in an infinite space, the stresses on the sloping free surface have been partly corrected by the introduction of the correction functions. Figures 9(a) and 9(b) show that, after being corrected, both the normal and shear stresses still have larger values in the range and then decrease with further increase of . The larger distance of , the smaller the stresses on the slope surface will be produced. The condition of zero stress at the horizontal and sloping free surface is not strictly satisfied. However, this approximate approach is purposeful in reality.

(a)  Normal stress

(b) Shear stress

Figure 9: Stress along the sloping free surface.

From Figures 8 and 9, the stresses, particularly the shear stresses, on the two free surfaces are close to zero after correction. Thus, the stresses will decrease until zero if the correction processes are iterated continually. It is obvious that the stresses on the free surface decline slowly with the iteration increase, so further trivial corrections are not carried out here.

The displacement of soil induced by the cavity expansion is discussed in this section. Figure 10 illustrates that the slope and its inclination angles have pronounced influence on the distributions of the displacements induced by a cavity expansion source. The displacement decreases with the increase of the distance . Similarly, the displacement decreases rapidly in the range , but then more slowly with further increase of . Meanwhile, the displacement induced by a cavity expansion source approaching a free surface, in general, is larger than that in infinite medium (no boundary effect). The findings are in accordance with the results proposed by Keer et al. [16] and Chai et al. [20]. Here, ( and are the displacements in the and directions, resp.).

Figure 10: The displacement varies with distance.

Using the theory of the spherical cavity expansion in an infinite space implies that there is an infinitely thick “soil wall” existing on the side of the slope, which will restrain the lateral displacement induced by the spherical expansion. According to the Figure 2, the slope tends to be flat with the increase of the angle β. As a result, the displacement at the same point in the soil decreases due to the increase of the thickness of “restriction.” For instance, when °, the displacement is 2.4 times of the displacement induced by cavity with no boundary effect at the instant of . Thus, the presence and inclination angle of the sloping free boundary have a great influence on the displacement due to cavity expansion.

6. Conclusions

Analytical solutions of the cavity expansion near the sloping ground were proposed based on the understanding that expansion caused by pile tip can be simulated as a spherical cavity expansion. Both the horizontal and sloping free surfaces are taken into account by using of a virtual image technique, harmonic functions, and the Boussinesq solutions, and the solutions will convert into the solutions reported by Keer et al. [16] when the sloping ground turn to the horizontal direction.

The results show that the presence and inclination angle of sloping free boundary have a considerable influence on the distributions of the stress and displacement fields induced by the spherical cavity expansion. As the distance from the cavity increases or when the boundary is approached, the Mises stress decreases. With the increase of the angle β between the slope and the vertical plane, the slope tends to flat and the displacement at the same point in the soil decreases with the increase of the “lateral restriction.” Likewise, the displacement increases with the decrease of the angle β. Therefore, the existence of the slope increases the risk of the slope failure due to the pile installation.

Acknowledgment

The work reported herein was supported by the National Natural Science Foundation of China (Grant no. 41272288). The above financial support is gratefully acknowledged.

References

K. Georgiadis and M. Georgiadis, “Undrained lateral pile response in sloping ground,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 136, no. 11, pp. 1489–1500, 2010.

C. W. W. Ng, L. M. Zhang, and K. K. S. Ho, “Influence of laterally loaded sleeved piles and pile groups on slope stability,” Canadian Geotechnical Journal, vol. 38, no. 3, pp. 553–566, 2001.

L. M. Zhang, C. W. W. Ng, and C. J. Lee, “Effects of slope and sleeving on the behavior of laterally loaded piles,” Soils and Foundations, vol. 44, no. 4, pp. 99–108, 2004.

C. Sagaseta and A. J. Whittle, “Prediction of ground movements due to pile driving in clay,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 127, no. 1, pp. 55–66, 2001.

S. L. Shen, J. Han, H. H. Zhu, and Z. S. Hong, “Evaluation of a dike damaged by pile driving in soft clay,” Journal of Performance of Constructed Facilities, vol. 19, no. 4, pp. 300–307, 2005.

R. F. Bishop, R. Hill, and N. F. Mott, “The theory of indentation and hardness tests,” Proceedings of the Physical Society, vol. 57, no. 3, pp. 147–159, 1945.

R. Hill, The Mathematical Theory of Plasticity, Oxford University Press, Oxford, UK, 1950.

H. S. Yu and G. T. Houlsby, “Finite cavity expansion in dilatant soils: loading analysis,” Geotechnique, vol. 41, no. 2, pp. 173–183, 1991.

M. F. Randolph, J. Dolwin, and R. Beck, “Design of driven piles in sand,” Geotechnique, vol. 44, no. 3, pp. 427–448, 1994.

M. F. Randolph, J. P. Carter, and C. P. Wroth, “Driven piles in clay—the effects of installation and subsequent consolidation,” Geotechnique, vol. 29, no. 4, pp. 361–393, 1979.

J. M. Pestana, C. E. Hunt, and J. D. Bray, “Soil deformation and excess pore pressure field around a closed-ended pile,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 128, no. 1, pp. 1–12, 2002.

S. L. Chen and N. Y. Abousleiman, “Exact undrained elasto-plastic solution for cylindrical cavity expansion in modified Cam Clay soil,” Geotechnique, vol. 62, no. 5, pp. 447–456, 2012.

J. C. Chai, N. Miura, and H. Koga, “Lateral displacement of ground caused by soil-cement column installation,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 131, no. 5, pp. 623–632, 2005.

C. Sagaseta, G. T. Houlsby, J. Norbury, and A. A. Wheeler, “Quasi-static undrained expansion of a cylindrical cavity in clay in the Presence of shaft friction and anisotropic initial stresses,” Report, Department of Engineering Science, University of Oxford, 1984.

C. Sagaseta, “Analysis of undrained soil deformation due to ground loss,” Geotechnique, vol. 37, no. 3, pp. 301–320, 1987.

L. M. Keer, Y. Xu, and V. K. Luk, “Boundary effects in penetration or perforation,” Journal of Applied Mechanics, vol. 65, no. 2, pp. 489–496, 1998.

R. D. Mindlin and D. H. Cheng, “Nuclei of strain in the semi-infinite solid,” Journal of Applied Physics, vol. 21, pp. 926–930, 1950.

H. G. Hopkins, “Dynamic expansion of spherical cavities in metals,” in Progress in Solid Mechanics, I. N. Sneddon and R. Hill, Eds., vol. 1, chapter 3, pp. 83–164, North-Holland, Amsterdam, The Netherlands, 1960.

A. S. Vesic, “Expansion of cavities in infinite soil mass,” Journal of the soil Mechanics and Foundations Division, vol. 98, no. 3, pp. 265–290, 1972.

J. C. Chai, J. P. Carter, N. Miura, and H. H. Zhu, “Improved prediction of lateral deformations due to installation of soil-cement columns,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 135, no. 12, pp. 1836–1845, 2009.

M. K. Kassir and G. C. Sih, Three-Dimensional Crack Problems, Noordhoff International Publishing, Dordrecht, The Netherlands, 1975.

P. P. Rao, J. P. Li, and Y. Liu, “Analysis of soil squeezing effect for spherical cavity expansion with inclined nonaxisymmetric displacement boundary,” Rock and Soil Mechanics, vol. 32, no. 9, pp. 2681–2687, 2011 (Chinese).

H. S. Yu, Cavity Expansion Method in Geomechanics, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2000.