Soil squeezing effect and formation disturbance caused by tunnel excavation can be simulated by cylindrical cavity expansion due to the comparability between tunneling and cavity expansion. Although most of the existing theoretical derivation is based on simple constitutive model of soil foundation, not only the relation between principal stress components was simplified in the solution process, but also the stress history, initial stress anisotropy, and stress-induced anisotropy of structural soil were neglected. The mechanical characteristics of soil are closely related to its stress history, so there is a gap between the above research and the actual engineering conditions. A three-dimensional elastoplastic solution of cylindrical cavity expansion is obtained based on the theory of critical state soil mechanics and engineering characteristics of shield tunneling. In order to fully consider the influence of initial anisotropy and induced anisotropy on the mechanical behavior of soils, the soil elastoplastic constitutive relation of cavity expansion is described in the course of K0-based modified Cam-clay (K0-MCC) model after soil yielding. An equation with equal number of variables is obtained under the elastic-plastic boundary condition based on the Lagrange multiplier method. By solving the extreme value of the original function, the analytical solution of radial, tangential, and vertical effective stresses distribution around the circular tunnel excavation is obtained. In addition, changes of elastic deformation area and plastic deformation area for soil during the shield excavation have been analyzed. Calculation results are compared with the numerical solutions which usually consider isotropic soil behavior as the basic assumption. In this paper, a constitutive model which is more consistent with the actual mechanical behavior of the soil and the construction process of the shield tunnel is considered. Therefore, the numerical solutions are more realistic and suitable for the shield excavation analysis and can provide theoretical guidance required for design of shield tunneling.
National Natural Science Foundation of China411722381. Introduction
Theoretical analysis of cylindrical cavity expansion theory in geotechnical engineering has been widely used in pile sinking, shield tunneling, static cone penetration test (CPT), pressuremeter test (PMT), and mixing pile construction disturbance problems [1–6]. Hill (1950) first proposed the spherical cavity expansion method and derived a general solution of stresses and displacements in Tresca material [7]. Vesic (1972) obtained the classical solutions considering the plastic zone volume change of ideal elastic-plastic expansion problem based on Mohr-Coulomb model [8]. Cao et al. (2001) presented closed form solutions for both spherical and cylindrical cavities by assuming the ultimate deviator stress distributions in the plastic zone and taking the in situ stress conditions into consideration [9]. Alonso et al. (2003) obtained the self-similar solution for the circular tunnel in strain-softening rock masses [10]. Wheeler et al. (2003) and Nakano et al. (2005) adjusted yield surface equation on the basis of K0-based modified Cam-clay, considering K0 consolidation caused by induced anisotropy and its evolution law in the loading process [11, 12]. Yin and Hicher (2008) identified parameters of modified Cam-clay model with viscosity of soil from the cavity expansion [13]. The problem of cavity expansion and cavity contraction has attracted much attention in geotechnical problems with application to the bearing capacity of deep foundations, interpretation of pressuremeter tests, breakout resistance of anchors, pile driving, wellbore instability, underground excavation, and blasting fracturing by explosives.
However, based on the above-mentioned researches, the conventional cylinder expansion theory assumes that initial stress is isotropic (static earth pressure coefficient K0=1.0). Due to the sedimentation history, mineralogical characteristics, and other several factors, the initial stress of soil layer is usually anisotropic. For example, the vertical stress and horizontal stress of the cylinder bore are not the same in the horizontal tunneling construction of the underground pipe laying, tunnel engineering. As a result, the traditional theory of cylindrical expansion generally used in practical engineering fails in explaining some practical phenomenon. Also, the traditional cylindrical expansion theory assumes that the cylindrical boundary condition is controlled by displacement, which presents around the cylinder in the form of symmetrical distribution as a function of distance. The traditional control theory can be applied in geotechnical problems analysis such as cylindrical pile and static cone penetration test. But in the practical engineering of shield tunnel, boundary condition of the cylinder expansion is controlled by stress and displacement of the cylinder, which was unable to satisfy the condition of axial symmetry. Hence, it may not be reasonable to explain such geotechnical problems considering traditional theory of cylindrical expansion.
In this paper, a three-dimensional elastoplastic solution of cylindrical cavity expansion is obtained based on the theory of critical state soil mechanics and engineering characteristics of shield tunneling to calculate the undrained cylindrical cavity expansion considered anisotropic clay of K0. This research work not only contributes to the theoretical values in field of geotechnical engineering but also has an important practical engineering significance.
2. Definition of Soil-Compacting Effect and Mechanical Model2.1. Definition of Soil-Compacting Effect
Shield tunneling is a typical three-dimensional problem. Soil deformation is closely related to the relative position of the shield machine. In actual shield construction, earth pressure and support pressure ahead of excavation face are not completely balanced [15]. When support pressure is greater than passive earth pressure, the soil is extruded and then soil-compacting effect presents. And the free surface may collapse if support pressure is less than active earth pressure. This paper mainly aims to study the soil-compacting effect produced when the support pressure is larger than passive earth pressure. Soil-compacting effect of shield tunneling is shown in Figure 1.
Schematic diagram of compacting soil by shield.
Excavation stage
Driving stage
Effect of initial stress field
(a) Along with the shield approach, the soil, which is squeezed by the excavation face, will be pushed around the shield. And then, the soil is output by the screw conveyor to form a void. The surrounding soil is excavated and unloaded, the stress is released to cause soil expansion, and the front soil moves to the void (Figure 1(a)).
(b) As the shield machine moves forward, the soil will be squeezed around the shield and excavation face. In the plane perpendicular to the axis of the tunnel, the soil is assume subjected to uniform radial extrusion force, which generates an equal amount of outward radial movement of the tunnel (Figure 1(b)).
(c) The extruded soil moves outwards of the shield machine and is also affected by the friction force arise from the driving process of shield machine. Radial deformation along the outer side of the tunnel finally formed under the influence of initial stress field (Figure 1(c)).
It can be seen that the squeezing effect of shield construction makes the surface subsidence or uplift and layered soil movement, soil stress, water content, pore water pressure, elastic modulus, Poisson’s ratio, strength, and bearing capacity of the physical and mechanical properties are likely to change. It is feasible to simulate the soil-compacting effect of shield tunneling by using cylindrical cavity expansion theory.
2.2. Mechanical Model
Figure 2 shows a cylindrical cavity expansion in an infinite K0 consolidated saturated clay with an in-plane initial horizontal pressure stress σh=K0σv and vertical pressure stress σv at infinity, as well as the hydrostatic pressure u0. Note that cylindrical polar coordinate system is adopted for present analysis, and the occurring compressive stresses and strains are taken as positive. As average internal pressure gradually increases from σh to σa, the soil adjacent to the cavity wall yields first. Further increase in average internal cavity wall pressure will lead to current cavity radius to the elastic-plastic (EP) boundary rp with a further increased internal cavity pressure with the continuation of the shield excavation. The symbol rp0 represents the initial position of the soil when the soil becomes plastic state.
Expansion of a cylindrical cavity in K0 consolidated anisotropic clays.
Under undrained condition, relationship between current radius rx, initial radius rx0 of a material particle, and current radius a, and initial radius a0 of the cylindrical cavity is as follows:(1)rx2-rx02=a2-a02.
At any stage of the shield excavation, any soil element within the surrounding soil mass satisfies the following equilibrium equation in both elastic and plastic regions (effective stress form): (2)∂σr′∂r+∂u∂r+σr′-σθ′r=0or alternatively in the total stress form(3)∂σr∂r+σr-σθr=0,where σr′ and σθ′ are effective radial and tangential stresses, respectively; σr and σθ are total radial and tangential stresses, respectively; u is pore water pressure.
2.3. Constitutive Relations
The K0 consolidated clay is assumed to be linearly elastic and infinitesimal deformation until the onset of yield. Thus, according to the Hook’s Law, the elastic stress-strain relationship can be expressed in the increment form as follows:(4)dεije=1+νEdσij-νEdσmmδij,where E and G are the elastic modulus and the shear modulus, respectively, which are defined in the K0-MCC model as follows [6]:(5)E=2G1+νG=31-2νυp′21+νκ,where κ is the slope of swelling line, υ is the specific volume, and ν is the Poisson’s ratio.
After yielding, elastoplastic behavior of soil can be described by large deformation theory and K0-MCC model, where a relative stress ratio η∗=η-η0 is adopted instead of the stress ratio η=q/p′ as used in the MCC model to consider the effect of initially stress-induced anisotropy on mechanical behavior of K0 consolidated clay. Yield function f for associative plasticity also serves as a plastic potential function g of the K0-MCC model can be expressed as follows [6, 11, 12]:(6)f=q-η0p′2+M2p′p′-pc′=0,where η0=31-K0/2K0+1 is stress ratio under K0 consolidated condition, M=6sinφ′/(3-sinφ′) is slope of critical state line, φ′ is effective friction angle, pc′ is effective yield pressure, and stress parameters p′ and q are defined respectively as follows:(7)p′=13σii′(8)q=32σij′-p′δijσij′-p′δij,where δij is Kronecker’s delta.
After taking strain hardening rule into consideration, yield function or plastic potential function can be given as follows:(9)f=λ-κ1+e0lnp′p0′+λ-κ1+e0ln1+q-η0p′2M2p′2-εvp=0,where λ and κ are slopes of consolidation and swelling lines, respectively; e0 is initial void ratio of soil for p′=p0′; εvp is plastic volumetric strain which is used as a hardening parameter.
As Figure 3 shows, plastic potential surface of K0-MCC model in p′-q space is a rotated ellipse, and degree of rotation portrays the extent of anisotropy. When shearing starts from isotropic stress η∗=η=q/p′ because η0=0, the K0-MCC model becomes the MCC model. The yield surface rotation can also be found in more detail in Yin et al. (2011, 2015) [16, 17].
Yield curves of anisotropic elastoplastic model.
For elastoplastic constitutive model, the total strain increments dεij are decomposed into elastic and plastic strain increments; i.e., (10)dεij=dεije+dεijp,where elastic strain increments dεijp can be still calculated using (3), and plastic strain increments dεijp in the K0-MCC model can be derived by taking account the associated plastic flow rule as follows:(11)dεijp=Ω∂f∂σij′,where the scalar multiplier Ω is given as follows:(12)Ω=-∂f/∂p′dp′+∂f/∂qdq∂f/∂εvp∂f/∂σii′,where(13)∂f∂p′=λ-κ1+e01p′M2p′2-q2-η02p′2M2p′2+q-η0p′2(14)∂f∂q=λ-κ1+e02q-η0p′M2p′2+q-η0p′2(15)∂p′∂σij′=δij3(16)∂q∂σij′=3σij′-p′δij2q(17)∂f∂εvp=-1(18)∂f∂σij′=∂f∂p′∂p′∂σij′+∂f∂q∂q∂σij′.
Substituting (13)–(16) into (18) yields the following equation:(19)∂f∂σij′=cpM2p′2-q2-η02p′2M2p′2+q-η0p′2δij3p′+3σij′-p′δijM2p′2+q-η0p′2q-η0p′q,where cp=λ-κ/1+e0.
Substituting (18) and (13)–(18) into (11), the scalar multiplier Ω can be obtained as follows:(20)Ω=dp′+2η-η0M2-η2-η02dq.
Then, detailed relationship between the plastic strain increment dεijp and the stress increment dσij can be expressed as follows:(21)dεijp=Ωcpp′2M2+η-η02M2-η2-η02p′δij3+3σi-p′δij1-η0/η.
3. Exact Numerical 3D Solution3.1. Elastic Zone
Based on the small-strain theory and (2) and (3), stress and displacement in the spherical coordinates of the elastic region can be obtained as follows [9, 14, 18]:(22)σr=σh+σp-σhrpr2(23)σθ=σh-σp-σhrpr2(24)σz=σv(25)Ur=σp-σh2G0rp2r(26)Δu=0,where σp is the total radial stress in the phase plane at EP boundary. It is worth noticing that there is no excess pore water pressure in plastic region, because effective mean stress is constant under undrained condition (dυ=31-2vdp′/E=0) and total mean stress, p, also keeps constant in elastic zone.
3.2. Elastoplastic Zone
For the problem of shield excavation, soil position around the tunnel cavity can be described as r, θ, and z in cylindrical polar coordinate system. Correspondingly, there are only variables of principal stress and principal strain without the deviator stress and deviator strain induced during shield excavation process. Therefore, elastoplastic constitutive relation of the K0-MCC model can be simplified to a much simple matrix equation:(27)dεrdεθdεz=1E+bar2-vE+baraθ-vE+baraz-vE+baraθ1E+baθ2-vE+baθaz-vE+bazar-vE+bazaθ1E+baz2dσr′dσθ′dσz′,where (28)ar=p′M2-η2-η023+3σr′-p′1-η0ηaθ=p′M2-η2-η023+3σθ′-p′1-η0ηaz=p′M2-η2-η023+3σz′-p′1-η0ηb=cpp′3M2+η-η02M2-η2-η02.
Taking the inverse calculation to (27) gives(29)dσr′dσθ′dσz′=1ξc11c12c13c21c22c23c31c32c33dεrdεθdεz,where(30)c11=1E21-v2+Eaθ2b+2Evaθazb+Eaz2bc12=1E2-Earaθ+vazb+v1+v-Eaθazb+Eaz2bc13=1E2-Earvaθ+azb+v1+v-Eaθazb+Eaθ2bc22=1E21-v2+Ear2b+2Evarazb+Eaz2bc23=1E2v+v2+Evar2b-Eaθazb-Evaraθ+azbc33=1E21-v2+Ear2b+2Evaraθb+Eaθ2bc21=c12c23=c32(31)ξ=-1+vE3-1+v+2v2+Eb-1+var2+aθ2+az2-2Evbaθaz+aθar+aθaz.
For undrained cylindrical cavity expansion, dεv=dεz; thus, according to the large deformation theory, radial and tangential strains can be given in incremental form as follows:(32)dεr=-dεθ=drr.
Substituting (32) into (29), constitutive matrix can be reduced to a set of first-order ordinary differential equations as follows:(33)dσr′dr-b11-b12ξr=0(34)dσθ′dr-b12-b22ξr=0(35)dσz′dr-b31-b32ξr=0which can be solved by Lagrangian analysis method as an initial value problem with r starting at rxp, provided that initial values of σr′rxp, σθ′rxp, and σz′rxp are given. Here rxp denotes the position of a specific particle which comes into plastic state instantly. All of the initial values required for solving the differential equations can be derived from the EP boundary conditions, which are further analyzed as follows.
3.3. Elastic-Plastic Boundary Conditions
As shown in (22), (23), (24), and (26), the change of mean effective stress, Δp′, is equal to zero in elastic region, which is still valid at the EP boundary; i.e.,(36)pp′=p0′,where pp′ denotes mean effective stress at the EP boundary.
Substituting (36) into (6), the deviator stress at the EP boundary, qp, can be obtained as follows:(37)qp=η0p0′+Mp0′OCR-1,where OCR is overconsolidation ratio, defined as pc0′/p0′, and pressure, pc0′, is the maximum mean preconsolidation stress.
According to (7), the deviator stress at the EP boundary can also be written as(38)qp=12σrp′-σθp′2+σzp′-σθp′2+σrp′-σzp′2,where σrp′, σθp′, and σzp′ are effective radial, tangential, and vertical stresses, respectively, at the EP boundary.
From (22)–(24) and (36), one obtains(39)σrp′+σθp′=σr0′+σθ0′=2σh′(40)σzp′=σz0′=σv0′,where σr0′, σθ0′, and σz0′ are initial effective radial, tangential, and vertical stresses, respectively.
For the K0 consolidated clay, relationship between in situ horizontal effective stress, σr0′ (or σθ0′), and in situ vertical stress, σz0′, can be given as follows:(41)σr0′=σθ0′=K0σz0′,where symbol K0 denotes the anisotropic degree of clay.
Substituting (39), (40), and (41) into (38), effective radial and tangential stresses at the EP boundary can be derived as (42)σrp′=σr0′+13qp2-1-K02σz0′2σθp′=σr0′-13qp2-1-K02σz0′2.
The position, rxp, where the specific particle becomes plasticity, can be derived by taking the displacement at EP boundary, Urp, and the undrained condition into consideration. According to (25), one can obtain(43)Urp=rxp-rx0=σrp′-σr0′2G0rxp.
Substituting (1) into (43) gives [14](44)rxp=2G0a2G0-σrp′-σr0′rxa2+a0a2-1.
Furthermore, current location of the EP boundary, rp, which is required for full analysis of stress distributions around the cavity, can be obtained by equating both rxp and rx to rp [18]:(45)rpa=a0/a2-1σrp′-σr0′/2G0-12-1.
So far, (40), (42), and (44) provide initial conditions for solving the first-order ordinary differential equations (33), (34), and (35).
After solving these first-order ordinary differential equations by numerical method, excess pore water pressure at an arbitrary location in the plastic zone, Δurx, can be computed by integrating (2) from the EP boundary rp to the point rx: (46)Δurx=σrp′-σrx′+∫rxrpσr′-σθ′rdr.
4. Closed Form 3D Analytical Solution
In the front section (refer to Figure 3), the semianalytical solution of expansion of the cylindrical hole has been presented with the limitation of high computation cost. In order to reduce the computational cost and achieve the purpose of practical application, this paper deduces a closed 3D total stress solution for cylindrical cavity undrained expansion in K0 consolidated clay.
Under the undrained condition, the total volume of the soil is equal to zero. That is,(47)κdp′υp′+λ-κdpc′υpc′=0.
Taking EP boundary conditions into consideration in (47), the following equation can be found: (48)p~c′=p~′OCRp~′p~0′-κ/λ-κ.
According to (37) and (48), the relationship between average effective stress and the plastic zone can be presented as follows:(49)q=η0p′+Mp′OCRp′p0′-1/Λ-11/2,where Λ is the plastic volumetric strain ratio, defined as (1-κ/λ).
For the K0-MCC model, critical state condition can be expressed as follows:(50)η=M2+η02
Substituting (50) into (49), ultimate effective mean stress, pf′, and ultimate deviator stress, qf, can be obtained as follows:(51)pf′=p0′M2OCRM2+η02-η02+M2Λ(52)qf=p0′M2+η02M2OCRM2+η02-η02+M2Λ.
On the other hand, the closed form 3D solutions in elastic region have been already derived, as shown in (22)–(26). Only the solutions in plastic region are derived in this study. Total radial stresses at an arbitrary location in plastic region, σrx, can be obtained by integrating (4) from the EP boundary, rp, to the point, rx,(53)σrx=σrp-23∫rprxqrdr,where σrp is the total radial stress at the EP boundary. The deviator stress, q, is related to the location rx, so that (53) cannot be integrated directly. Therefore, in order to obtain the closed form 3D solution, the following two assumptions have to be made at first.
(a) Relationship between Deviator Stress and Radial Distance in the Plastic Region. The deviator stress, q, in the plastic zone is assumed to be equal to the ultimate deviator stress, qf, which is similar to q.
(b) Relationship between and among Total Radial, Tangential, and Vertical Stresses in the Plastic Region. The vertical total stress σz is assumed to be the mean of radial stress σr and tangential stress σθ after soil yields.
Although these assumptions seem to be bold and strong, the approximate analytical results match the exact numerical results fairly well in the plastic region, which will be compared and discussed in detail later in the paper.
According to assumption (a), substituting (25) into (45) and ignoring higher order terms of qf2-1-K02σv02/G, position of the EP boundary, which is related to the current cavity radius, can be written as (54)rpa2=3Gqf2-1-K02σv021-a0a2.
Then, taking assumption (a) into consideration, radial stress in the plastic zone can be derived from (53) and can be written as follows:(55)σrx=σh0+qf2-1-K02σv0′231+ln3Gqf2-1-K02σv0′2a2-a02rx2.
According to the previous two assumptions, substituting (55) into (8) and taking (52) into consideration, total tangential and vertical stresses in the plastic region can be obtained:(56)σθ=σh0+qf2-1-K02σv0′23ln3Gqf2-1-K02σv0′2a2-a02rx2-1σz=σh0+qf2-1-K02σv0′23ln3Gqf2-1-K02σv0′2a2-a02rx2.
Furthermore, excess pore water pressure in plastic region can be derived by effective stress principle:(57)Δurx=K0-1p0′1+2K0+qf2-1-K02σv0′23ln3Gqf2-1-K02σv0′2a2-a02rx2+p0′1-M2OCRM2+η02+M2Λ.
Equating rx to a, pressure-expansion relationship and excess pore water pressure at cavity wall can be derived from (55) and (57), respectively, as follows:(58)σra=σh0+qf2-1-K02σv0′23ln3Gqf2-1-K02σv0′21-a02a2+1(59)Δua=K0-1p0′1+2K0+qf2-1-K02σv0′23ln3Gqf2-1-K02σv0′21-a02a2+p0′1-M2OCRM2+η02+M2Λ.
From (58), as the cavity continues to expand to infinity (a/a0→∞), the limit cavity wall pressure, σrult, and excess pore water pressure, Δuault, for undrained cavity expansion in K0 consolidated anisotropic clay can be obtained:(60)σrult=σh0+qf2-1-K02σv0′23ln3Gqf2-1-K02σv0′2+1(61)Δuault=K0-1p0′1+2K0+qf2-1-K02σv0′23ln3Gqf2-1-K02σv0′21-a02a2+p0′1-M2OCRM2+η02+M2Λ.
It is interesting to note here that, taking qf2-1-K02σv0′2 as q-f, the present solutions are similar to previous MCC model based solutions obtained by Cao et al. [9]. Especially for isotropic clay (K0=1), the closed form 3D solutions presented in this paper reduce to the solutions obtained by Cao et al. [9], which demonstrates that the solution obtained by Cao et al. [9] is just a special case of the presented solution. Actually, for nature clay, the horizontal in situ stress is different from the vertical one, and the soil is initial anisotropy (i.e., K0≠1). Note that the term of qf2-1-K02σv0′2 in the presented solution reflects the effect of initially stress-induced anisotropy on shield excavation. Therefore, presented solutions can yield more reasonable results for cylindrical cavity expansion in anisotropy clay. In addition, the term of K0-1p0′/1+2K0 in (57), (59), and (61) implies that initial anisotropy also has notable effects on excess pore water pressure during shield excavation process.
5. Calculation Result Analysis and Discussion
In order to verify the accuracy and suitability of this method considering the influencing factors of initial anisotropy and induced anisotropy, calculation results of present method and Chen and Abousleiman’s solution [18] using MCC constitutive model results were, respectively, validated in this paper.
Physical and mechanical parameters of soil used in this paper are taken from Chen and Abousleiman [18]. Table 1 shows the properties of three typical K0 consolidated clays.
Physical and mechanical parameters of soil [14].
M=1.2, λ=0.15, κ=0.03, v=0.278
K0
OCR
σv0 (kPa)
G0 (kPa)
υ0
u0 (kPa)
0.625
1
160
4348
2.09
100
1
3
120
4113
1.97
100
2
10
72
3756
1.8
100
5.1. Pressure-Expansion Relationship
As shown in Figure 4, comparison of pressure-expansion relationship around shield tunnel is calculated by the proposed 3D closed form solution and Chen’s solution. All calculated stresses have been normalized. Effective radial stress and vertical stress have been overestimated, while the tangential stress is noted to be underestimated in the plastic zone based on Chen’s solution apart from K0=1. Therefore, the anisotropy induced by initial stress has a significant influence on the stress distribution in the plastic zone. If the tunnel is in the isotropic stress conditions, that is to say K0=1, K0-MCC constitutive model will degenerate into a MCC constitutive model; in this case, the results of the two calculations are in complete agreement.
Comparison of pressure-expansion relationships around shield tunnel calculated by 3D closed form solution and Chen’s solution.
σz′– expansion relationship
σr′– expansion relationship
σθ′– expansion relationship
Δu– expansion relationship
It is also interesting to note in Figure 4 that vertical stress is always equal to the average of radial and tangential stresses in critical state zone for all cases, which can be proved exactly as follows. Combining (33)–(35) gives (62)dσr′+σθ′-2σz′dr=-31-2var-aθσr′+σθ′-2σz′bEξr1-η0η.
Since the stress components keep constant in the critical region, therefore, the left side of (62) is equal to zero; thereby (62) can be reformulated as(63)σz′=12σr′+σθ′.
Equation (63) indicates that the second assumption for approximate closed form 3D solution holds in the critical region.
It can also be seen from Figure 4(d) that excess pore water pressure is underestimated near the cavity wall but overestimated as the radial distance increases. Besides, the excess pore water pressure keeps constant in critical region, and thus the second assumption for the approximate closed form 3D solution still holds well in the critical region. When the shield is compacted to a certain degree, excess pore water pressure in the soil squeezing effect zone is gradually stabilized, and the same level of pressure can be achieved, regardless of the value of K0.
5.2. Shield Excavation Process Analysis
According to cylindrical cavity expansion theory, effect of soil compaction caused by shield tunneling is mainly revealed in dynamic change of the normalized excess pore water pressure Δua/p0′, the soil pressure inside the cavity σa/p0′, and the cavity radius rp/a. Figure 5 shows variations of normalized cavity pressure, soil pressure, and excess pore water pressure with cavity radius during shield excavation process. Results computed by the three different typical values of K0 (0.625, 1, and 2) [14], as well as the results from the MCC model based solutions, are presented in the same figures for comparison. Similar to previous analysis, all calculation results have been normalized.
Evolution of cavity pressure and excess pore water pressure with cavity radius during shield excavation process.
K0=0.625
K0=1
K0=2
It can be seen from Figure 5 that excess pore water pressure and internal cavity pressure are affected by the initial stress-induced anisotropy significantly when a/a0>2. In the process of gradual expansion of inner wall of the tunnel, the result of MCC constitutive model is smaller than that of the K0-MCC constitutive model. Stress values at the inner wall of the tunnel and excess pore water pressure around the tunnel increase dramatically when a/a0<2. Both parameters gradually achieve a stable level as the shield continues to advance. Compared to Figures 5(a), 5(b), and 5(c), it is clear that, as the K0 value increases, internal cavity pressure caused by shield tunneling is effectively suppressed, while the impact of excess pore water pressure and plastic radius is not very obvious. This is because, at the initial stage of pore expansion, negative pore pressure first occurs at the hole wall. Then, with the increase of pore size, the excess pore pressure will turn to a positive value from negative to a critical state. Furthermore, when the surrounding soil is in isotropic initial stress state condition (K0=1), the K0-MCC constitutive model is reduced to the MCC constitutive model, and the results are consistent with numerical results. However, effect of anisotropy degree on discrepancy of excess pore water pressure calculated from the two solutions is not obvious, which implies that the OCR also has a pronounced effect on excess pore water pressure at cavity wall during cavity expansion process.
Evolution of normalized plastic radius, rp/a, with normalized instant cavity radius, a/a0, for three different values of K0 and OCR, together with the MCC model based results is also plotted in Figure 5 for comparison. It is seen that, for normally consolidated clay with K0=0.625, plastic radius during whole cavity expansion process is significantly affected by initially stress-induced anisotropy. However, for heavily overconsolidated clay, discrepancy between presented solution and MCC model based solution is not apparent, because OCR also plays a major role in plastic radius for heavily consolidated clay during cavity expansion process.
5.3. Effective Stress Path in p′-q Plane
Figure 6 shows effective stress path (ESP) in p′-q plane for a soil particle located at cavity wall in three typical K0 consolidated clays, where two stress parameters p′ and q have been normalized by p0′. Note that points O, P, and F in Figure 6 denote the in situ stress point, yield stress point, and failure stress point, respectively. It is also clear that initial yield loci are rotational ellipses for anisotropic clays (K0≠1), and degree of rotation reflects degree of anisotropic. For isotropic clay (K0=1), K0 line overlaps with p′ axis, and yield locus is symmetric to p′ axis, and thus K0-MCC model reduces to MCC model. For anisotropic clays, yield locus and CSL of the K0-MCC model are higher than that of MCC model. Therefore, yield and failure stresses during cavity expansion process are underestimated by previous MCC model based solutions.
ESP in p′-q plane for a soil particle located at cavity wall during shield excavation process.
OCR = 1, K0=0.625
OCR = 3, K0=1
OCR = 10, K0=2
For normally (OCR=1) consolidated clay with K0=0.625 in Figure 6, yielding occurs immediately after cavity expansion, so that in situ stress point O overlaps the yield stress point P. The corresponding ESP starts form point O which is located in K0 line and moves up to the left until ESP reaches critical state line (CSL) at point F. For lightly (OCR=3) and heavily overconsolidated clay (OCR=10), ESP firstly moves vertically and soil keeps elastic until the ESP approaches initial yield locus at point P. Then soil yields, and the ESP turns up to the right until reaching CSL at point F.
5.4. Stress Distribution of Shield Tunneling under Different K0 Conditions
Figure 7 shows stress distribution of shield tunneling under different K0 conditions calculated by the approximate closed form 3D solution, where normalized vertical total stress σz/p0′, mean of radial stress σr/p0′, tangential stress σθ/p0′, and excess pore water pressure Δua/p0′ are plotted against normalized radial distance r/a0. It is seen when r/a<2, the rise of K0 leads to a significant change of stress distribution, which will reach a relatively stable state after EP boundary presents. Besides the value of K0 has no effect on the final stress except the vertical total stress σz/p0′. Stress values under different K0 values finally reach the same magnitude near the cavity wall. It also indirectly verifies that the two basic assumptions of the closed form 3D solution is feasible and reasonable. Therefore, the closed form 3D solution is exact enough to interpret and solve general geotechnical problems.
Stress distribution of shield tunneling under different K0 conditions calculated by the approximate closed form 3D solution.
σz
σr
σθ
Δu
6. Conclusions
This paper aims to study the soil-compacting effect produced by shield tunneling by a precise semianalytical solution. Particular emphasis is given to develop an approximate closed 3D solution of cylindrical cavity undrained expansion in saturated anisotropic clay based on the K0-based modified Cam-clay model. Furthermore, efforts have been made to examine (a) the effects of the initial anisotropic and initially stress-induced anisotropy on shield excavation and (b) stress distribution of shield tunneling under different K0 conditions calculated by the approximate closed form 3D solution. Some of key observation and findings from the study are as follows.
(a) Both numerical solution and closed from 3D solution presented in this study improve the understanding of conventional solutions by considering anisotropic properties of natural clay and reducing MCC model based solutions when K0=1. Therefore, presented solutions can yield a more realistic stress field induced by cylindrical cavity expansion in natural clay.
(b) Vertical total stress σz, radial stress σr, and tangential stress σθ become intermediate principal stress and major and minor principal stress, respectively, when soil came into critical state, and vertical total stress σz is equal to the mean of radial stress σr and tangential stress σθ in critical state region around the cavity for all cases.
(c) Initially stress-induced anisotropy has significant influence on stress field and plastic radius around the cavity, especially in normally consolidated clay. Excess pore water pressure is not only affected by the initially stress-induced anisotropy but also by the initial anisotropy of clay. The stress state of the soil around the borehole wall changes greatly due to the expansion of the column hole in the K0 consolidated clay, and the stress around the hole decays logarithmically with the increase of the radial distance and tends to the stress state at the time of departure.
(d) The initial stress anisotropy and stress history have obvious influence on the expansion process of cylindrical cavity, especially for the ultimate expansion pressure of cylindrical cavity expansion process. Stress field caused by shield excavation is overestimated, whereas plastic radius during cavity expansion process is underestimated by MCC constitutive model based on solutions in anisotropic clays.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (no. 41172238). China Railway 20 Bureau Group Co., Ltd., has also given a lot of practice guidance in the study of the paper.
DurbanD.Large strain solution for pressurized elasto/plastic tubes197946122823010.1115/1.3424511Zbl0408.730382-s2.0-0018446495De Sousa CoutinhoA. G. F.Radial expansion of cylindrical cavities in sandy soils: application to pressuremeter tests19902767377482-s2.0-002557082010.1139/t90-087LadanyiB.In-situ determination of undrained stress-strain behavior of sensitive clays with the pressuremeter19729331331910.1139/t72-034PalmerA. C.Undrained plane-strain expansion of a cylindrical cavity in clay: A simple interpretation of the pressuremeter test19722234514572-s2.0-000086582310.1680/geot.1972.22.3.451LeeY.-K.PietruszczakS.A new numerical procedure for elasto-plastic analysis of a circular opening excavated in a strain-softening rock mass20082355885992-s2.0-4304910612510.1016/j.tust.2007.11.002ChangM. F.TehC. I.CaoL. F.Undrained cavity expansion in modified Cam clay II: Application to the interpretation of the piezocone test200151433535010.1680/geot.51.4.335.393902-s2.0-0035351952HillR.1950London, UKOxford University PressMR1650059VesicA. C.Expansion of cavities in infinite soil mass19729826590CaoL. F.TehC. I.ChangM. F.Undrained cavity expansion in modified Cam clay I: theoretical analysis200151432333410.1680/geot.51.4.323.393952-s2.0-0035351936AlonsoE.AlejanoL. R.VarasF.Fdez-MańinG.Carranza-TorresC.Ground response curves for rock masses exhibiting strain-softening behaviour20032713115311852-s2.0-024263660310.1002/nag.315Zbl1098.74039WheelerS. J.NäätänenA.KarstunenM.LojanderM.An anisotropic elastoplastic model for soft clays20034024034182-s2.0-003875758110.1139/t02-119NakanoM.NakaiK.NodaT.AsaokaA.Simulation of shear and one-dimensional compression behavior of naturally deposited clays by super/subloading Yield Surface Cam-clay model20054511411512-s2.0-34249706081YinZ.-Y.HicherP.-Y.Identifying parameters controlling soil delayed behaviour from laboratory and in situ pressuremeter testing20083212151515352-s2.0-5044909288510.1002/nag.684Zbl1273.74327ChenS. L.AbousleimanY. N.Exact undrained elasto-plastic solution for cylindrical cavity expansion in modified cam clay soil201262544745610.1680/geot.11.p.0272-s2.0-84861393165MarshallA. M.Tunnel-pile interaction analysis using cavity expansion methods201213810123712462-s2.0-8487951920410.1061/(ASCE)GT.1943-5606.0000709YinZ.-Y.KarstunenM.ChangC. S.KoskinenM.LojanderM.Modeling time-dependent behavior of soft sensitive clay201113711110311132-s2.0-8375519657710.1061/(ASCE)GT.1943-5606.0000527YinZ.-Y.YinJ.-H.HuangH.-W.Rate-Dependent and Long-Term Yield Stress and Strength of Soft Wenzhou Marine Clay: Experiments and Modeling201533179912-s2.0-8490728655310.1080/1064119X.2013.797060ChenS. L.AbousleimanY. N.Exact drained solution for cylindrical cavity expansion in modified cam clay soil201363651051710.1680/geot.11.p.0882-s2.0-84878091239