Undrained Elastoplastic Solution for Cylindrical Cavity Expansion in Structured Cam Clay Soil Considering the Destructuration Effects

: Soil structure has signiﬁcant inﬂuences on the mechanical behaviors of natural soils, although it is rarely considered in previous cavity expansion analyses. This paper presents an undrained elastoplastic solution for cylindrical cavity expansion in structured soils, considering the destructuration effects. Firstly, a structural ratio was deﬁned to denote the degree of the initial structure, and the Structured Cam Clay (SCC) model was employed to describe the subsequent stress-induced destructuration, including the structure degradation and crushing. Secondly, combined with the large strain theory, the considered problem was formulated as a system of ﬁrst-order differential equations, which can be solved in a simpliﬁed procedure with the introduced auxiliary variable. Finally, the signiﬁcance and efﬁciency of the present solution was demonstrated by comparing with the previous solutions, and parametric studies were also conducted to investigate the effects of soil structure and destructuration on the cavity expansion process. The results show that the soil structure has pronounced effects on the mechanical behavior of structured soils around the cavity. For structured soils, a cavity pressure that is larger than the corresponding reconstituted soils when the cavity expands to the same radius is required, and the effective stresses ﬁrst increase to a peak value before decreasing rapidly with soil structure degradation and crushing. The same ﬁnal critical state is reached for soils with different degrees of the initial structure, which indicates that the soil structure is completely destroyed during the cavity expansion. With the increase of the destructuring index, the soil structure was destroyed more rapidly, and the stress release during the plastic deformation became more signiﬁcant. Moreover, the present solution was applied in the jacking of a casing during the sand compact pile installation and in situ self-boring pressuremeter (SBPM) tests, which indicates that the present solution provides an effective theoretical tool for predicting the behavior of natural structured soils around the cavity. Contributions: Conceptualization, Y.Z. and T.L.; methodology, Z.Z.; software, Z.Z. and S.X.; validation, Z.Z.; data curation, Z.Z.; writing—original draft preparation, Z.Z.; writing—review and editing, Y.Z. and S.X.;


Introduction
Cavity expansion theory provides a useful, accurate, and simple tool for modelling many complex geotechnical problems [1], such as the interpretation of the results of pressuremeter and cone penetration tests [2][3][4], the estimation of stress responses around installed piles [5][6][7], as well as the stability and deformation analyses for wellbore drilling and tunnel excavation [8][9][10][11]. Therefore, it has received wide attention and made great progress over the past few decades.
In previous studies, a variety of constitutive models were employed to develop analytical/semi-analytical cavity expansion solutions [12]. The early solutions were mostly based on the Mohr-Coulomb criterion [2,13,14] and the Drucker-Prager criterion [15], in which the soil was assumed to be the elastic-perfectly plastic materials. Consequently, represents the destructuration of soils. For isotropic compression virgin yield (p > p yi ), the additional voids ratio ∆e of structured soils can be given by [37] where p is the current mean effective stress, p yi is the initial yield stress of structured soils under isotropic compression, and b is the isotropic compression destructuring index. ∆e i is the value of ∆e at p = p yi , and c is a stable value of ∆e when p is sufficiently large.
Appl. Sci. 2022, 12, x 3 of 25 naturally structured soils e is larger than that of the corresponding reconstituted soils e* at the same stress state [39], the isotropic compression behavior of structured soils in the SCC model is idealized, as shown in Figure 1, and the decrease of Δe (the difference between e and e*) represents the destructuration of soils. For isotropic compression virgin yield (p′ > p′yi), the additional voids ratio Δe of structured soils can be given by [37] ( ) e c c p (1) where p′ is the current mean effective stress, p′yi is the initial yield stress of structured soils under isotropic compression, and b is the isotropic compression destructuring index. Δei is the value of Δe at p′ = p′yi, and c is a stable value of Δe when p′ is sufficiently large. As a Cam Clay model, it is assumed that the yield surface of structured soils in the p′-q plane is elliptical, as shown in Figure 2, which is larger than that of the corresponding reconstituted soils, and can be formulated as [37] ( ) where M, named the critical state stress ratio, represents the slope of the critical state line (CSL), q is the deviator stress that follows its standard three-dimensional definition in this study, and p′s is the isotropic yield stress denoting the size of the yield surface. As a Cam Clay model, it is assumed that the yield surface of structured soils in the p -q plane is elliptical, as shown in Figure 2, which is larger than that of the corresponding reconstituted soils, and can be formulated as [37] f = q 2 − M 2 p p s − p = 0 (2) where M, named the critical state stress ratio, represents the slope of the critical state line (CSL), q is the deviator stress that follows its standard three-dimensional definition in this study, and p s is the isotropic yield stress denoting the size of the yield surface.
Appl. Sci. 2022, 12, x 3 of 25 naturally structured soils e is larger than that of the corresponding reconstituted soils e* at the same stress state [39], the isotropic compression behavior of structured soils in the SCC model is idealized, as shown in Figure 1, and the decrease of Δe (the difference between e and e*) represents the destructuration of soils. For isotropic compression virgin yield (p′ > p′yi), the additional voids ratio Δe of structured soils can be given by [37] ( ) e c c p (1) where p′ is the current mean effective stress, p′yi is the initial yield stress of structured soils under isotropic compression, and b is the isotropic compression destructuring index. Δei is the value of Δe at p′ = p′yi, and c is a stable value of Δe when p′ is sufficiently large. As a Cam Clay model, it is assumed that the yield surface of structured soils in the p′-q plane is elliptical, as shown in Figure 2, which is larger than that of the corresponding reconstituted soils, and can be formulated as [37] ( ) where M, named the critical state stress ratio, represents the slope of the critical state line (CSL), q is the deviator stress that follows its standard three-dimensional definition in this study, and p′s is the isotropic yield stress denoting the size of the yield surface. The initial yield surface of structured soils enlarges with the increase of structure, as shown in Figure 2. The equivalent yield surface for structured soils, i.e., the dashed line in Figure 2, is defined as a yield surface for the soils remolded with the same voids ratio and stress state [37]. The size of the equivalent yield surface is denoted by p e , whose value can be determined by the states (p and e) of structured soils, as follows: where e * IC is the voids ratio for reconstituted soils when p = 1 kPa, as shown in Figure 1. λ and κ denote the slopes of the compression and loading-reloading lines of the reconstituted soils, respectively.
When the stress state reaches the structural yield surface with dp s = 0, virgin yielding occurs and is accompanied by the destructuration, which is an irreversible process. Considering the contribution of shearing on the destructuration of soils, the plastic volumetric strain increment can be described as where η = q/p is the stress ratio, γ is the shearing destructuring index, and < > is the Macaulay bracket. Additionally, it was observed that the flow rule, the strain increment ratio dε p d /dε p v , is also significantly affected by the soil structure [40]. The non-associated flow rules, defining the yield surface as inconsistent with the plastic potential surface, are employed in the SCC model to incorporate the influences of structure. For virgin yielding with η < M, hardening occurs, and the flow rule is given by where ω is a material constant describing the flow rule. As suggested by Carter and Liu [37], the value of ω is normally set to 1 when the experimental data is insufficient. For virgin yielding with η > M, softening occurs. During this process, the yield surface shrinks with the destructuration of soils. As a result, it is necessary to change the sign of the plastic deviatoric strain associated with destructuration to ensure the strain increment direction is always outwards of the structured yield surface [35,37]. In this situation, the flow rule can be expressed as Structured soils that undergo virgin yielding may reach a critical state (η = M), but the structure is not completely destroyed (∆e = 0), which means that the plastic deformation will develop continuously with the further loading [38]. This stage is termed the crushing of soil structure and was first introduced into the SCC model by Carter and Liu [37]. During the structure crushing, the stress state always varies along the CSL, which follows Appl. Sci. 2022, 12, 440

of 24
Meanwhile, the additional voids ratio ∆e gradually diminishes with the structure crushing, and tends to be c when the stress state reaches a final critical state. During this process, the flow rule is given by where sign( ) denotes the signum function.

Problem Description
The undrained expansion of a cylindrical cavity with an initial radius a 0 in an infinite structured soil is schematically shown in Figure 3, where σ h and σ v denote the initial horizontal and vertical total stress, respectively, and u 0 is the initial pore pressure. As the internal pressure σ a increases gradually from its initial value σ h , the cavity expands from a 0 to a larger cavity radius a. Meanwhile, a typical soil element moves outward along the radial direction from the initial position r x0 to the current position r x . During this process, the soil element first undergoes elastic deformation until the initial yielding occurs, and then experiences plastic hardening or softening with structure degradation. The soil structure begins to be crushed if the soil element with an incompletly destroyed structure may attain a critical state before reaching the current position, r x . In Figure 3, r xp and r xc denote the positions of the soil element at the precise moment when it enters the plastic state and when its structure begins to be crushed. Appl. Sci. 2022, 12, x 5 of 25 Meanwhile, the additional voids ratio Δe gradually diminishes with the structure crushing, and tends to be c when the stress state reaches a final critical state. During this process, the flow rule is given by where sign( ) denotes the signum function.

Problem Description
The undrained expansion of a cylindrical cavity with an initial radius a0 in an infinite structured soil is schematically shown in Figure 3, where σh and σv denote the initial horizontal and vertical total stress, respectively, and u0 is the initial pore pressure. As the internal pressure σa increases gradually from its initial value σh, the cavity expands from a0 to a larger cavity radius a. Meanwhile, a typical soil element moves outward along the radial direction from the initial position rx0 to the current position rx. During this process, the soil element first undergoes elastic deformation until the initial yielding occurs, and then experiences plastic hardening or softening with structure degradation. The soil structure begins to be crushed if the soil element with an incompletly destroyed structure may attain a critical state before reaching the current position, rx. In Figure 3, rxp and rxc denote the positions of the soil element at the precise moment when it enters the plastic state and when its structure begins to be crushed. Cylindrical coordinates (r, θ, z) are considered because of the symmetry of the cylindrical cavity expansion problem. Considering a soil element located at the radial position r with unit height in z direction, the radial and tangential total stresses, σr and σθ, satisfy the force equilibrium equation as follows: Cylindrical coordinates (r, θ, z) are considered because of the symmetry of the cylindrical cavity expansion problem. Considering a soil element located at the radial position r with unit height in z direction, the radial and tangential total stresses, σ r and σ θ , satisfy the force equilibrium equation as follows: According to the Terzaghi effective stress principle, Equation (9) can be written as where σ r and σ θ are the effective components corresponding to σ r and σ θ , respectively, and u is the pore pressure at position r.

Elastic Analysis
The structured soils undergo elastic deformation before the stress state reaches the initial yield surface. The elastic stress-strain relationship obeys Hooke's law, and can be written in the following incremental form where E is the Young's modulus, v is the effective Poisson's ratio, and δ ij is the Kronecker delta. In the SCC model, v is treated as a constant, and E can be expressed in terms of v and shear modulus G, namely E = 2G(1 + v), where G is related to the mean effective stress p and the voids ratio e is Combining the elastic stress-strain relationship (Equation (11)) and the effective equilibrium equation (Equation (10)), three effective stress components, σ r , σ θ , σ z , the radial displacement U r , and the pore pressure u can be explicitly expressed as [1,22]: where σ rp is the specific value of σ r at the elastoplastic boundary r p .

Elastoplastic Analysis
For the natural soils, the plastic deformation is usually accompanied by the elastic deformation; therefore, the components of the total strain increment dε i consists of the elastic strain increment dε e i and plastic strain increment dε p i , i.e., where dε e i can be calculated using Equation (11), and dε p i is controlled by the yield function, the hardening (or softening) law, and the plastic flow rules of the SCC model.
In order to derive the components of plastic strain increment dε p i , the consistency condition, i.e., df = 0, is applied to the yield function in Equation (2), and yields where Therefore, the plastic volumetric strain increment dε p v is given as where g denotes the plastic potential function, Λ is the plastic multiplier, and During the hardening (η < M) and softening (η > M) processes, the plastic volumetric strain can be described using Equation (4), and yields ∂p s ∂ε During the crushing process (η = M), the effective stress states satisfy Equation (7). Therefore, ∂p s ∂ε The total volumetric strain increment dε v is zero under the undrained condition, which leads to the elastic strain and plastic strain, which can be compensated for each other, namely, dε v p = −dε e v . According to Equation (11), the term ∂p /∂ε p v in Equation (25) can be obtained as follows: ∂p ∂ε By combining Equation (22) and the non-associated flow rules Equations (5), (6), or (8), the plastic deviator strain increment dε p d can be expressed in the following unified form: where β = ∂ε p d /∂ε p v represents the plastic flow rule, which follows the different forms in the processes of hardening (β = β 1 ), softening (β = β 2 ), and crushing (β = β 3 ).
Combining Equations (22) and (27), the components of the plastic strain increments can be obtained as follows: Rewriting Equation (28) in matrix form, the following can be obtained: where Substituting Equations (11) and (29) into Equation (18), the relationship between total strain components and effective stress components can be expressed as follows: By taking the inverse form of Equation (32), the components of stress increment can be expressed in terms of the components of strain increment as where Under the undrained and plane-strain assumptions, the volumetric strain increment dε v and vertical strain increment dε z are zero, which leads to dε r = −dε θ , and therefore Equation (33) becomes According to the large strain theory, the tangential strain increment dε θ in Equation (35a-c) can be expressed in the logarithmic strain form as Substituting Equation (36) into Equation (35a-c), it can be seen that three stress components, σ r , σ θ , and σ z , are functions of radial distance r. Note that the expressions of C ij and Γ in Equation (35a-c) is inconsistent in the hardening/softening and crushing processes. For the hardening/softening process of the structured soils, the variations of the stress components can be calculated by taking the values of σ r (r xp ), σ θ (r xp ), and σ z (r xp ) as the initial values. If the soil element attains the critical state but with ∆e = c, the crushing of the structure will occur at a subsequent stage, which needs to be further solved by taking the values of σ r (r xc ), σ θ (r xc ), and σ z (r xc ) as the initial values. However, it is complex to determine r xc because it depends not only on the current radial location r, but also on the nonlinear displacement U r during the hardening/softening process.
To simplify the solution procedure, the auxiliary variable ξ, defined by Chen and Abousleiman [23] as U r /r, is skillfully introduced, where U r = r − r 0 is the radial displacement for a specific soil element. This auxiliary variable takes the incremental form as follows: Combining Equations (36) and (37), it can be seen that the tangential strain increment dε θ is merely related to ξ; therefore, Equation (35a-c) can be written as Equation (38a-c) formulate the problem as a system of first-order differential equations with respect to a single variable ξ, which are valid for any soil element around the cavity. The values of σ r , σ θ , and σ z at the instant of ξ = ξ p and ξ = ξ c are required as the initial values of the hardening/softening and crushing processes, where ξ p and ξ c correspond to the values of u r /r at the positions r xp and r xc , respectively. The values of σ r (ξ c ), σ θ (ξ c ), and σ z (ξ c ) as the initial values of the crushing process can be determined from the hardening/softening process, owing to the continuity of deformation; therefore, the values of σ r (ξ p ), σ θ (ξ p ), and σ z (ξ p ) are sufficient to solve the stress variation for a specific soil element. Meanwhile, because the cavity expansion is geometrically self-similar, the soil elements at any position r follow the same stress path [41]. Consequently, only the stress variation of the soil element at the cavity wall that undergoes the longest stress path needs to be calculated. For other soil elements around the cavity, the stress state only depends on ξ, which can be found from the experienced stress state of the soil element at the cavity wall by seeking the same ξ.

Initial Values at the Elastoplastic Boundary
To determine the initial values σ r (ξ p ), σ θ (ξ p ), and σ z (ξ p ), prior knowledge of the initial stress state and the degree of over-consolidation of structured soils is needed. As shown in Figure 4, (p 0 , q 0 ) and (p 0 , q p ) denote the initial stress point and initial yield stress point, respectively. The reference point is the size of a virtual yield surface passing through the point (p 0 , q 0 ). The isotropic over-consolidation ratio of the reconstituted soils, R, is defined as by Chen and Abousleiman [22], which is different from the conventional over-consolidation ratio, OCR = σ v,max /σ v [6]. For structured soils, a higher stress level is required to reach the initial yield surface than that for the corresponding reconstituted soils due to p s > p C , which means that the structured soils exhibit partial 'false' overconsolidated behavior [34]. To denote the degree of structure, the initial structural ratio, R s , is defined as follows: soils due to S C p p ′ ′ > , which means that the structured soils exhibit partial 'false' over-consolidated behavior [34]. To denote the degree of structure, the initial structural ratio, Rs, is defined as follows: It can be seen that Rs is greater than 1 for structured soils. When the soil is in a reconstituted state, S p′ equals C p′ ; therefore, Rs = 1. The stronger the initial structure, the greater Rs for the soils.
It should be mentioned that the stress components at the elastoplastic boundary rp (corresponding to ξ = ξp) satisfy both the elastic solutions given by Equations (13)- (15) and the yield function given by Equation (2). Combining this with the definition of Rs, σ′r(ξp), σ′θ(ξp), and σ′z(ξp) can be easily obtained following the method of Chen and Abousleiman [22]: It can be seen that R s is greater than 1 for structured soils. When the soil is in a reconstituted state, p s equals p C ; therefore, R s = 1. The stronger the initial structure, the greater R s for the soils.
It should be mentioned that the stress components at the elastoplastic boundary r p (corresponding to ξ = ξ p ) satisfy both the elastic solutions given by Equations (13)- (15) and the yield function given by Equation (2). Combining this with the definition of R s , σ r (ξ p ), σ θ (ξ p ), and σ z (ξ p ) can be easily obtained following the method of Chen and Abousleiman [22]: The initial values determined from Equations (40)- (44) are sufficient to solve the system of Equation (38a-c). Note that these equations are related to the auxiliary variable ξ; therefore, the obtained stress components vary with the ξ instead of the radial position r.

Determination of the Radial Position and Excess Pore Pressure
For the completeness of the present solution, the relationship between ξ and r needs to be established. For the undrained and plane-strain conditions, the incompressibility of the soil skeleton and pore fluid renders the current position of any soil element r as related to its initial position r 0 , its initial cavity radius a 0 , and its current cavity radius a, which can be formulated as Combining Equation (45) with ξ = (r − r 0 )/r, the relationship between ξ and r can be obtained as: In particular, the position of the elastoplastic boundary r p can be obtained by setting ξ in Equation (46) as ξ p , and then the explicit relationship between r p and the cavity radii (a and a 0 ) can be written as Once the effective stress variations are determined, the pore pressure u at any radial position r x can be calculated from the equilibrium equation. By integrating Equation (10) from r p to r x , u is expressed as Therefore, the excess pore pressure ∆u(r x ) generated around the cavity can be easily obtained as u(r x ) − u 0 .

Results and Discussion
Following the solution procedure shown in Figure 5, the present solution was implemented in the numerical software MATLAB 2018b. Firstly, the present solution was compared with the previous solutions both in reconstituted and structured soils [22,36]. Subsequently, parametric analyses were performed to investigate the cavity expansion responses for the soils with different structure degrees R s and isotropic over-consolidation ratios R. Finally, the influences of the structural parameter on the cavity expansion responses were analyzed.
The present solution was performed in Boston Blue Clay because it has been widely used in the analyses of cavity expansion [22][23][24]. The basic Cam Clay parameters (M, λ, κ, v) and initial state parameters (R, σ h , σ v , e 0 , G 0 , and u 0 ) of Boston Blue Clay are listed in Tables 1 and 2, respectively. Because of the absence of experimental data in the literature, the values of the structural parameters (b, c, ω, γ) were estimated in a reasonable range. In addition, three different initial structural ratios, R s = 1.0, 1.2, and 1.5, were adopted, where the special case of R s = 1.0, was chosen to highlight the difference between structured and reconstituted soils. Table 1. Soil property parameters used in calculation.

Basic Cam Clay Parameters
Structural Parameters Table 2. Soil state parameters used in calculation.  The present solution was performed in Boston Blue Clay because it has been widely used in the analyses of cavity expansion [22][23][24]. The basic Cam Clay parameters (M, λ, κ, v) and initial state parameters (R, σh, σv, e0, G0, and u0) of Boston Blue Clay are listed in Tables 1 and 2, respectively. Because of the absence of experimental data in the literature, the values of the structural parameters (b, c, ω, γ) were estimated in a reasonable range. In addition, three different initial structural ratios, Rs = 1.0, 1.2, and 1.5, were adopted, where the special case of Rs = 1.0, was chosen to highlight the difference between structured and reconstituted soils.

Basic Cam Clay Parameters Structural Parameters
where n is as large as possible, e.g. n≥200 Solving the system of differential equations (38a)-(38c) Until η=M, and then setting β=β3 Solving the system of differential equations (38a)-(38c) Until ξ=ξp, and then convert ξ into the radial position r Calculating the pore pressure from equation (48) Yes No Calculating the stress components in the elastic zone

Comparisons between the Present Solution and Previous Solution
Theoretically, the existing undrained solution based on the MCC model [22] is a special case of the present solution based on the SCC model. This is because the SCC model can evolve to the MCC model by setting the relevant structural parameters as sufficiently small values, i.e., ∆e→c, c→0, b→0, γ→0, and ω→0. In this situation, the initial yield surface of structured soils, as shown in Figure 4, tends to be the same as that of the corresponding reconstituted soils; therefore, the initial structural ratio R s tends to be 1.
The present solution was implemented in the case of R s →1, c→0, b→0, γ→0, and ω→0, and then compared with the existing undrained solution in reconstituted soils of Chen and Abousleiman [22], as shown in Figure 6. Note that only the initial state parameters of normally consolidated soils (R = 1) in Table 2 were adopted, and the stresses in Figure 6a were normalized by the undrained shear strength s u (=54.3 kPa) for convenience of comparison. It can be seen that whether the distributions of the stress components along the radial distance, or the effective stress path for a soil element at the cavity wall, the results of present solution for structured soils are highly consistent with those of Chen and Abousleiman [22] for reconstituted soils. Thus, the solution proposed in this study can be regarded as a unified solution for undrained cavity expansion in structured and reconstituted soils.
in Figure 6a were normalized by the undrained shear strength su (=54.3 kPa) for convenience of comparison. It can be seen that whether the distributions of the stress components along the radial distance, or the effective stress path for a soil element at the cavity wall, the results of present solution for structured soils are highly consistent with those of Chen and Abousleiman [22] for reconstituted soils. Thus, the solution proposed in this study can be regarded as a unified solution for undrained cavity expansion in structured and reconstituted soils. The present solution was further compared with the undrained solution in structured soils proposed by Li et al. [36] to show the significances of destructuration. For the sake of convenience, the case of R = 3 in Table 2, which was calculated in [36], was selected, and all of the stresses were normalized by the undrained shear strength su (=115.0 kPa). Note that the same initial values at the elastoplastic boundary were adopted to ensure the soils containing the same initial structure. Figure 7 showed the comparisons of distributions of the mean effective stress p′, deviator stress q, and excess pore pressure Δu. It can be seen that soil structure has significant influences on the stress distributions around the cavity. Although the influences of the structure were considered in the solution of Li et al. [36], the values of p′ and q tends to be constant with the decrease of radial distance r/a in the plastic zone. This illustrated that only the initial structure was accounted for in their solution, and that the structure degradation was neglected, which resulted in the inability to capture the softening behavior caused by the destructuration. In contrast, the present solution could reflect the effects of structure degradation and crushing. Therefore, the p′ and q calculated in the present solution could finally reach the corresponding values in reconstituted soils at the cavity wall r/a = 1 after experiencing The present solution was further compared with the undrained solution in structured soils proposed by Li et al. [36] to show the significances of destructuration. For the sake of convenience, the case of R = 3 in Table 2, which was calculated in [36], was selected, and all of the stresses were normalized by the undrained shear strength s u (=115.0 kPa). Note that the same initial values at the elastoplastic boundary were adopted to ensure the soils containing the same initial structure. Figure 7 showed the comparisons of distributions of the mean effective stress p , deviator stress q, and excess pore pressure ∆u. It can be seen that soil structure has significant influences on the stress distributions around the cavity. Although the influences of the structure were considered in the solution of Li et al. [36], the values of p and q tends to be constant with the decrease of radial distance r/a in the plastic zone. This illustrated that only the initial structure was accounted for in their solution, and that the structure degradation was neglected, which resulted in the inability to capture the softening behavior caused by the destructuration. In contrast, the present solution could reflect the effects of structure degradation and crushing. Therefore, the p and q calculated in the present solution could finally reach the corresponding values in reconstituted soils at the cavity wall r/a = 1 after experiencing strain softening, indicating that the structure of the soils at the cavity wall was destroyed completely.

Internal Cavity Pressure and Excess Pore Pressure at the Cavity Wall
In this section, the characteristics of the cavity expansion curve (i.e., variations of the total and effective internal cavity pressure, σa and σ′a) and the variations of excess pore pressure at the cavity wall were analyzed. Note that the results of Rs = 1.0 are

Internal Cavity Pressure and Excess Pore Pressure at the Cavity Wall
In this section, the characteristics of the cavity expansion curve (i.e., variations of the total and effective internal cavity pressure, σ a and σ a ) and the variations of excess pore pressure at the cavity wall were analyzed. Note that the results of R s = 1.0 are equivalent to that of Chen and Abousleiman [22] (i.e., the results in the MCC model), which would be taken as a benchmark to highlight the influences of the initial structure and subsequent destructuration of soils on these variations.
As shown in Figures 8 and 9, both the total internal cavity pressure σ a and the excess pore pressure ∆u increase sharply at the initial stage of cavity expansion, but then become slowly with an increase in a/a 0 . When a/a 0 is greater than 5.0, both σ a and ∆u reach constant values. The comparisons between these results for soils with different R s indicate that the structure has significant influences on the variations of σ a and ∆u. Specifically, the greater the R s , the larger σ a is required and the larger ∆u will be generated when the cavity expands to the same a/a 0 . In addition, for heavily over-consolidated soils (R = 10), a negative excess pore pressure is generated at the instant of cavity expansion, but it dissipates immediately with the further expansion, as shown in Figure 9.

Internal Cavity Pressure and Excess Pore Pressure at the Cavity Wall
In this section, the characteristics of the cavity expansion curve (i.e., variations of the total and effective internal cavity pressure, σa and σ′a) and the variations of excess pore pressure at the cavity wall were analyzed. Note that the results of Rs = 1.0 are equivalent to that of Chen and Abousleiman [22] (i.e., the results in the MCC model), which would be taken as a benchmark to highlight the influences of the initial structure and subsequent destructuration of soils on these variations.
As shown in Figures 8 and 9, both the total internal cavity pressure σa and the excess pore pressure Δu increase sharply at the initial stage of cavity expansion, but then become slowly with an increase in a/a0. When a/a0 is greater than 5.0, both σa and Δu reach constant values. The comparisons between these results for soils with different Rs indicate that the structure has significant influences on the variations of σa and Δu. Specifically, the greater the Rs, the larger σa is required and the larger Δu will be generated when the cavity expands to the same a/a0. In addition, for heavily over-consolidated soils (R = 10), a negative excess pore pressure is generated at the instant of cavity expansion, but it dissipates immediately with the further expansion, as shown in Figure 9.  The variations in the effective internal cavity pressure σ′a with a/a0 are plotted in Figure 10. For Rs = 1.0, the variation trends of σ′a along a/a0 are consistent with those of the total internal cavity pressure σa. For Rs = 1.2 and 1.5, σ′a first increases to the peak values, and then decreases gradually with a further increase in a/a0. Finally, σ′a reaches the same value as that of the corresponding reconstituted soils (Rs = 1.0), although the peak value of σ′a increases with Rs. The softening behavior of structured soils can be attributed to the release of accumulated stress due to the degradation and crushing of the structure. The same effective stress state is reached for the structured and the corre- The variations in the effective internal cavity pressure σ a with a/a 0 are plotted in Figure 10. For R s = 1.0, the variation trends of σ a along a/a 0 are consistent with those of the total internal cavity pressure σ a . For R s = 1.2 and 1.5, σ a first increases to the peak values, and then decreases gradually with a further increase in a/a 0 . Finally, σ a reaches the same value as that of the corresponding reconstituted soils (R s = 1.0), although the peak value of σ a increases with R s . The softening behavior of structured soils can be attributed to the release of accumulated stress due to the degradation and crushing of the structure. The same effective stress state is reached for the structured and the corresponding reconstituted soils, which indicates that the structure of soils at the cavity wall was completely destroyed. In fact, σ a can be calculated from the measured values of σ a and ∆u, which can be used to estimate the degree of destructuration for structured soils, as mentioned by Sivasithamparam and Castro [31]. The variations in the effective internal cavity pressure σ′a with a/a0 are plotted in Figure 10. For Rs = 1.0, the variation trends of σ′a along a/a0 are consistent with those of the total internal cavity pressure σa. For Rs = 1.2 and 1.5, σ′a first increases to the peak values, and then decreases gradually with a further increase in a/a0. Finally, σ′a reaches the same value as that of the corresponding reconstituted soils (Rs = 1.0), although the peak value of σ′a increases with Rs. The softening behavior of structured soils can be attributed to the release of accumulated stress due to the degradation and crushing of the structure. The same effective stress state is reached for the structured and the corresponding reconstituted soils, which indicates that the structure of soils at the cavity wall was completely destroyed. In fact, σ′a can be calculated from the measured values of σa and Δu, which can be used to estimate the degree of destructuration for structured soils, as mentioned by Sivasithamparam and Castro [31].

Distributions of Stress and Excess Pore Pressure around the Cavity
The stress distributions along the radial distance, including the stress components σ′r, σ′θ, σ′z, mean effective stress p′, and deviator stress q, as well as the excess pore pressure Δu, are plotted in Figures 11 and 12. These stress distributions at the instant of a/a0 = Figure 10. Relationship between the effective internal cavity pressure σ a and the cavity radius a/a 0 .

Distributions of Stress and Excess Pore Pressure around the Cavity
The stress distributions along the radial distance, including the stress components σ r , σ θ , σ z , mean effective stress p , and deviator stress q, as well as the excess pore pressure ∆u, are plotted in Figures 11 and 12. These stress distributions at the instant of a/a 0 = 2 for R s = 1.2 and 1.5 were also compared with those for R s = 1.0 (i.e., the solution in [22]) to investigate the influences of the initial structure and the subsequent stress-induced destructuation. Note that, in these figures, the radial distance r is normalized by the current cavity radius a, and then plotted in logarithmic coordinates to highlight the stress features near the cavity.
From Figure 11, it can be seen that the distributions of σ r , σ θ , and σ z for R s = 1.0, 1.2, and 1.5 are consistent in the far field but have significant differences in the near field around the cavity. According to the stress state, three zones, i.e., the elastic zone, plastic zone, and critical state zone, can be generally defined around the cavity. For reconstituted soils (R s = 1.0), σ r , σ θ , and σ z vary monotonously in the plastic zone with a decrease in r/a, and tend to be constant in the critical state zone. For structured soils (R s = 1.2, 1.5), σ r first increases to a maximum value corresponding to the peak strength of structured soils, and then decreases gradually with a decrease in r/a until the value is the same as that for R s = 1.0. It can also be seen that the larger the R s that is adopted, the larger the magnitudes of σ θ and σ z that are obtained in the plastic zone. However, the final critical state values are the same for structured and reconstituted soils, irrespective of the initial structure.
The influences of soil structure on the development of the stress-state zones around the expanding cavity can also be seen in Figure 11. For normal consolidated reconstituted soils (R = 1 with R s = 1.0), there is no elastic zone, as shown in Figure 11a. However, for structured soils (R s = 1.2, 1.5), the elastic zone appears because of the partial 'false' over-consolidated behavior of structured soils. For all cases, it can be seen that, as R s increases, the position of the elastoplastic boundary gradually shifts to the left, the ranges of the plastic zone and critical state zone become narrow, and the critical state zone even disappears in some cases. Figure 12 shows the distributions of the mean effective stress p , the deviator stress q, as well as the excess pore pressure ∆u around the cavity. It can be seen that, for reconstituted soils (R s = 1.0), p and q remain almost unchanged in the vicinity of the cavity wall, which indicates that the soils have reached the critical state. As r/a increases, p and q may increase or decrease in the plastic zone, which actually depends on the values of R. For structured soils (R s = 1.2 and 1.5), q increases rapidly with a decrease in r/a, and then decreases to a stable value with the destructuration of soils after reaching the peak value. The variations in p can be divided into two patterns, for normally consolidated and lightly over-consolidated soils (Figure 12a,b); p remains a constant in the elastic zone and decreases monotonously with the decrease in r/a during the plastic deformation. For moderately and heavily overconsolidated soils (Figure 12c,d), p increases to a peak value with the accumulation of the plastic strain and then decreases to the critical state value, which indicates that the soil undergoes plastic softening deformation after yielding.
Generally, when examining the stresses from the far field to the cavity wall, it is observed that the stresses for the soils with different R s are consistent in the elastic zone, and that they then separate from each other approximately in the plastic zone, and finally come together again near the cavity wall. This is because the initial stresses for soils with a certain R and different R s are the same, while the structure affects the stress responses of the soils and enhances the differences between structured and reconstituted soils. After the peak strength, the differences between the stress states of the soils with different R s decrease and, finally, tend to be zero with the degradation and crushing of the structure.
2 for Rs = 1.2 and 1.5 were also compared with those for Rs = 1.0 (i.e., the solution in [22]) to investigate the influences of the initial structure and the subsequent stress-induced destructuation. Note that, in these figures, the radial distance r is normalized by the current cavity radius a, and then plotted in logarithmic coordinates to highlight the stress features near the cavity.
From Figure 11, it can be seen that the distributions of σ′r, σ′θ, and σ′z for Rs = 1.0, 1.2, and 1.5 are consistent in the far field but have significant differences in the near field around the cavity. According to the stress state, three zones, i.e., the elastic zone, plastic zone, and critical state zone, can be generally defined around the cavity. For reconstituted soils (Rs = 1.0), σ′r, σ′θ, and σ′z vary monotonously in the plastic zone with a decrease in r/a, and tend to be constant in the critical state zone. For structured soils (Rs = 1.2, 1.5), σ′r first increases to a maximum value corresponding to the peak strength of structured soils, and then decreases gradually with a decrease in r/a until the value is the same as that for Rs = 1.0. It can also be seen that the larger the Rs that is adopted, the larger the magnitudes of σ′θ and σ′z that are obtained in the plastic zone. However, the final critical state values are the same for structured and reconstituted soils, irrespective of the initial structure.
The influences of soil structure on the development of the stress-state zones around the expanding cavity can also be seen in Figure 11. For normal consolidated reconstituted soils (R = 1 with Rs = 1.0), there is no elastic zone, as shown in Figure 11a. However, for structured soils (Rs = 1.2, 1.5), the elastic zone appears because of the partial 'false' over-consolidated behavior of structured soils. For all cases, it can be seen that, as Rs increases, the position of the elastoplastic boundary gradually shifts to the left, the ranges of the plastic zone and critical state zone become narrow, and the critical state zone even disappears in some cases. It can also be seen that the variation trends of excess pore pressures Δu with r/a are similar for Rs = 1.0, 1.2, and 1.5, i.e., Δu remains zero in the elastic zone, and increases sharply and nearly linearly with the decrease in the logarithm of r/a in the plastic and critical zones. Meanwhile, the gradient of Δu versus lnr/a becomes larger for soils with larger Rs. For heavily over-consolidated soils (R = 10), negative excess pore pressures Δu are generated around the cavity, which illustrates that shear dilation occurs during cavity expansion.

Stress Path for a Soil Element around the Cavity
For the soil element located at any position around the cavity, the same effective stress path (ESP) in p′-q plane is followed; therefore, only the ESP for a soil element at the cavity wall is analyzed. The ESPs for R = 1, 1.2, 3, and 10 are plotted in Figure 13, respectively. In these figures, points A, B, C, and D denote four stress variation points, i.e., the initial stress point (p0, q0), the initial yield stress point, the crushing beginning point, and the final critical state point, respectively. The B, B′, and B″ were used to distinguish the initial yield stresses for Rs = 1.0, 1.2, and 1.5, respectively. Note that the comparisons of ESPs between the present solution and Chen and Abousleiman's solution [22] were shown by comparing the case of Rs = 1.0 with that of Rs = 1.2 and 1.5.
From Figure 13, it is clear that the initial yield surface and yield stress increase as Rs increases. During the elastic deformation, the ESPs follow the vertical line until the stress state reaches the initial yield stress point B, B′, or B″, except for the ESP for R = 1 with Rs = 1.0. Thereafter, the ESPs turn left or right with the development of plastic deformation, which depends on whether the soils experience hardening (the stress ratio at B, B′, or B″ is smaller than M) or softening (the stress ratio at B, B′, or B″ is larger than M).
For Rs = 1.0, the ESPs go monotonously to the left or right after the initial yielding, until the final critical state point D is reached. For Rs = 1.2 and 1.5, the ESPs turn left and reach the point C′ or C″ when hardening occurs, and then travel along the CSL, owing to the structure crushing, until reaching the final critical state point D, as shown in Figure 11a,b. In contrast, the ESPs turn right and move toward the critical line when softening occurs, and then change their directions and travel downwards until they reach point D, as shown in Figure 13c,d. That the same point D is reached for the soils with different Rs indicates that the final critical state is irrespective of the initial structure, although it has significant influences on the ESPs, which can also be found in the shearing test results of structured soils [37]. In addition, it can be seen from Figure 13b that because of the 'false' over-consolidated behavior of structured soils, the slightly over-consolidated soils with R = 1.2 exhibit partial characteristics of the heavily It can also be seen that the variation trends of excess pore pressures ∆u with r/a are similar for R s = 1.0, 1.2, and 1.5, i.e., ∆u remains zero in the elastic zone, and increases sharply and nearly linearly with the decrease in the logarithm of r/a in the plastic and critical zones. Meanwhile, the gradient of ∆u versus lnr/a becomes larger for soils with larger R s . For heavily over-consolidated soils (R = 10), negative excess pore pressures ∆u are generated around the cavity, which illustrates that shear dilation occurs during cavity expansion.

Stress Path for a Soil Element around the Cavity
For the soil element located at any position around the cavity, the same effective stress path (ESP) in p -q plane is followed; therefore, only the ESP for a soil element at the cavity wall is analyzed. The ESPs for R = 1, 1.2, 3, and 10 are plotted in Figure 13, respectively. In these figures, points A, B, C, and D denote four stress variation points, i.e., the initial stress point (p 0 , q 0 ), the initial yield stress point, the crushing beginning point, and the final critical state point, respectively. The B, B , and B were used to distinguish the initial yield stresses for R s = 1.0, 1.2, and 1.5, respectively. Note that the comparisons of ESPs between the present solution and Chen and Abousleiman's solution [22] were shown by comparing the case of R s = 1.0 with that of R s = 1.2 and 1.5.

Influences of Structural Parameters
In this section, the influences of structural parameters (b, γ, c, ω) on the cavity expansion responses were further investigated. Herein, the soil property parameters in Table 1 and the state parameters in the case of R = 1 were chosen, and the calculation results were displayed in terms of effective stress path (ESP) in the p′-q plane.
The isotropic compression destructuring index b controls the decrease of the additional voids ratio Δe of structured soils, which represents the destructuration caused by the isotropic compression. As suggested by Carter and Liu [35], the values of b depend on soil structure and are generally b ≥ 1 for soft clay and b < 1 for stiff clays. Herein, b = 0.25, 1, and 30 were adopted to analyze the influences of b on the ESPs. The calculated results for b = 0.25, 1, and 30 were shown in Figure 14. As b increases, the ESP experiences a lower stress level, indicating that the stress release becomes more significant. From Figure 13, it is clear that the initial yield surface and yield stress increase as R s increases. During the elastic deformation, the ESPs follow the vertical line until the stress state reaches the initial yield stress point B, B , or B , except for the ESP for R = 1 with R s = 1.0. Thereafter, the ESPs turn left or right with the development of plastic deformation, which depends on whether the soils experience hardening (the stress ratio at B, B , or B is smaller than M) or softening (the stress ratio at B, B , or B is larger than M).
For R s = 1.0, the ESPs go monotonously to the left or right after the initial yielding, until the final critical state point D is reached. For R s = 1.2 and 1.5, the ESPs turn left and reach the point C or C when hardening occurs, and then travel along the CSL, owing to the structure crushing, until reaching the final critical state point D, as shown in Figure 11a,b. In contrast, the ESPs turn right and move toward the critical line when softening occurs, and then change their directions and travel downwards until they reach point D, as shown in Figure 13c,d. That the same point D is reached for the soils with different R s indicates that the final critical state is irrespective of the initial structure, although it has significant influences on the ESPs, which can also be found in the shearing test results of structured soils [37]. In addition, it can be seen from Figure 13b that because of the 'false' overconsolidated behavior of structured soils, the slightly over-consolidated soils with R = 1.2 exhibit partial characteristics of the heavily over-consolidated soils when the R s increases to 1.5, namely, the initial yield stress point B was located at the left side of CSL, as shown in Figure 13b.

Influences of Structural Parameters
In this section, the influences of structural parameters (b, γ, c, ω) on the cavity expansion responses were further investigated. Herein, the soil property parameters in Table 1 and the state parameters in the case of R = 1 were chosen, and the calculation results were displayed in terms of effective stress path (ESP) in the p -q plane.
The isotropic compression destructuring index b controls the decrease of the additional voids ratio ∆e of structured soils, which represents the destructuration caused by the isotropic compression. As suggested by Carter and Liu [35], the values of b depend on soil structure and are generally b ≥ 1 for soft clay and b < 1 for stiff clays. Herein, b = 0.25, 1, and 30 were adopted to analyze the influences of b on the ESPs. The calculated results for b = 0.25, 1, and 30 were shown in Figure 14. As b increases, the ESP experiences a lower stress level, indicating that the stress release becomes more significant. The shearing destructuring index γ reflects the influences of shear stress on the destructuration of soils. For the several natural and artificial structured soils reported in Carter and Liu [37], the values of γ are in the range of 0.5-2.1. Herein, γ = 0.5, 1, and 2 are adopted to analyze its influences on the cavity expansion responses. The calculation results are displayed in Figure 15, from which the same variation trends as Figure 14 can be observed. Obviously, γ together with b make prominent contributions to the destructuration of soils. The structural parameter c is defined as the stable value of the additional voids ratio Δe when p′ is sufficiently large during the isotropic compression. Actually, the soil structure is almost completely removed under extremely large pressure, which means that the value of c is very small. Taking the natural Eastern Osaka Clay reported in The shearing destructuring index γ reflects the influences of shear stress on the destructuration of soils. For the several natural and artificial structured soils reported in Carter and Liu [37], the values of γ are in the range of 0.5-2.1. Herein, γ = 0.5, 1, and 2 are adopted to analyze its influences on the cavity expansion responses. The calculation results are displayed in Figure 15, from which the same variation trends as Figure 14 can be observed. Obviously, γ together with b make prominent contributions to the destructuration of soils. The shearing destructuring index γ reflects the influences of shear stress on the destructuration of soils. For the several natural and artificial structured soils reported in Carter and Liu [37], the values of γ are in the range of 0.5-2.1. Herein, γ = 0.5, 1, and 2 are adopted to analyze its influences on the cavity expansion responses. The calculation results are displayed in Figure 15, from which the same variation trends as Figure 14 can be observed. Obviously, γ together with b make prominent contributions to the destructuration of soils. The structural parameter c is defined as the stable value of the additional voids ratio Δe when p′ is sufficiently large during the isotropic compression. Actually, the soil structure is almost completely removed under extremely large pressure, which means that the value of c is very small. Taking the natural Eastern Osaka Clay reported in The structural parameter c is defined as the stable value of the additional voids ratio ∆e when p is sufficiently large during the isotropic compression. Actually, the soil structure is almost completely removed under extremely large pressure, which means that the value of c is very small. Taking the natural Eastern Osaka Clay reported in Adachi et al. [42] as an example, Carter and Liu [37] concluded that the value of c for this clay can be assumed as zero. The parametric study shows that little differences occurred between various c values.
The structural parameter ω describes the flow rule, which is normally taken as 1 when the experimental data is insufficient, as suggested by Carter and Liu [37]. Actually, the values of ω also have little effect on the results.

Applications of Proposed Solution in Geotechnical Problems
In this section, two practical geotechnical problems are analyzed using the present solution. The first one is the jacking of a casing during the sand compaction pile installation, which is a typical penetration problem and is commonly modelled as the expansion of a cylindrical cavity. The second one is the in situ SBPM tests, in which the cylindrical probe is self-bored into soils to minimize the disturbance, and then the rubber membrane covering the probe expands the drilling.

Estimation of Stress Variations Caused by the Casing Installation
Yi et al. [43] conducted a centrifuge model test to investigate the installation of a sand compaction pile. Before the 'compacted' sand column is formed, a cylindrical casing is jacked into the reconstituted Malaysian Kaolin clay, in which the water level is 0.25-m (prototype scale) higher than the ground. In centrifuge modelling, the total stress and pore pressure transducers were inserted at 4.6 times the casing radius and at prototype depths of 1.5 m, 2.75 m, and 4 m. The measured total radial stress and pore pressure at the instant of the casing tip's arrival at these specific depths was extracted.
The present solution was used to predict the variations of total radial stress and pore pressure caused by the casing installation. The parameters of Malaysian Kaolin clay are taken from Purwana et al. [44] and Yi et al. [43]. The basic Modified Cam Clay parameters are M = 0.9, λ = 0.244, κ = 0.053, e * IC = 2.35, and v = 0.33. The lateral stress coefficient K 0 and effective unit weight γ are 0.6 and 6 kN/m 3 , respectively. The initial shear modulus G 0 is taken as 20p 0 , as in Yi et al. [43]. The Malaysian Kaolin clay used in the centrifuge model test is reconstituted; therefore, it is reasonable to take the structural parameters (b, γ, c, ω) as sufficiently small. Because there was no cavity existing before the casing was jacked, the value of a/a 0 should be set large enough (larger than 5, see Figure 8) so that the cavity pressure is stable.
Taking the above parameters as input data, the distribution of effective radial stress and pore pressure can be obtained. According to the location of the transducers (at the radial distance r/a = 4.6), the predicted radial effective stress and excess pore pressure was extracted from the stress distribution. The predicted stresses were compared with the corresponding measured values, as shown in Figure 16. It can be seen that the stress variations caused by the casing installation can be well predicted by the present solution.
Taking the above parameters as input data, the distribution of effective radial stress and pore pressure can be obtained. According to the location of the transducers (at the radial distance r/a = 4.6), the predicted radial effective stress and excess pore pressure was extracted from the stress distribution. The predicted stresses were compared with the corresponding measured values, as shown in Figure 16. It can be seen that the stress variations caused by the casing installation can be well predicted by the present solution. Figure 16. Calculated against measured total radial stress and pore pressure during the casing installation.

Interpretation of Pressuremeter Tests
The in situ SBPM test results reported by Rouainia et al. [45] were selected to further illustrate the applicability of the present solution in structured soils. These SBPM tests were performed at depths of 14 m and 20 m in London Clay, and the relationships between the cavity strain and pressure on the inner borehole wall were well documented.
The relevant parameters of London Clay can be found from Gasparre and Coop [46] and Rouainia et al. [45], i.e., M = 0.87, λ = 0.097, κ = 0.003, v = 0.3. The bulk unit weight below the water table (5.6-m depth at the test site) is 20 kN/m 3 . The over-consolidation ratio at the depths 14 m and 20 m are 3.0 and 4.0, respectively. According to the isotropic compression test results, the values of e * IC is estimated to be 1.4. The lateral stress coefficients for normally consolidated and over-consolidated conditions are estimated using [47]: where ϕ is the effective friction angle. For London Clay, ϕ varies in the range of 20-25 • [48], and its average value is adopted here. As suggested by Liu and Carter [35], the value of b is smaller than 1 for stiff clays, and therefore, 0.8 is adopted for the stiff London Clay in the calculation. The values of c, γ, and ω refer to Carter and Liu [37], and are taken as 0.005, 1.0, and 1.5. Because the initial structure is enhanced with depth, R s = 3.0 and 3.2 are adopted for 14 m and 20 m, respectively.
The variations of the pressure on the inner borehole wall with the cavity strain were predicted using the present solution. Adopting the aforementioned parameters, the relationships between the total internal cavity pressure and cavity radial can be calculated using the present solution, as shown in Figure 8. It is known that the pressure applied by the probe to the inner borehole wall is equivalent to the total internal cavity pressure, and the cavity strain is (a − a 0 )/a 0 . Thus, the results of SBPM tests can be predicated from the cavity expansion curves by converting the cavity radius into the cavity strain. As shown in Figure 17, the general trend of SBPM test response curves are well captured, and the present solution predicted a pressure on the inner borehole wall very close to the measured data from the SBPM tests.
using the present solution, as shown in Figure 8. It is known that the pressure applied by the probe to the inner borehole wall is equivalent to the total internal cavity pressure, and the cavity strain is (a − a0)/a0. Thus, the results of SBPM tests can be predicated from the cavity expansion curves by converting the cavity radius into the cavity strain. As shown in Figure 17, the general trend of SBPM test response curves are well captured, and the present solution predicted a pressure on the inner borehole wall very close to the measured data from the SBPM tests.

Conclusions
For the problem of cylindrical cavity expansion in naturally structured soils, an undrained elastoplastic solution was presented by using the SCC model of Carter and Liu [37]. Although the flow rules in the SCC model follow the different forms during the hardening, softening, and crushing processes, a unified system of differential equations was established for the considered problem, which can be solved in a simplified procedure by introducing an auxiliary variable. The present solution was compared with the previous solutions in reconstituted and structured soils to show the validity and significances of this study. Subsequently, parametric studies were conducted to explore the characteristics of cavity expansion responses in the structured soils, and then the influences of structural parameters were investigated. Finally, the casing installation and in situ SBPM tests were analyzed using the present solution. The following conclusions were drawn: (1) Because the initial structure and destructuration effects are sufficiently considered, the present solution is capable of capturing the typical mechanical responses of naturally structured soils around the cavity. With the increase of the initial structure, the yielding stress increases and the 'false' over-consolidated behavior is strengthened. Moreover, as a result of the destructuration, the soils with different initial structures will reach the same critical state when the structure is completely destroyed; (2) The soil structure has significant influences on the cavity responses. As the cavity expands, the effective internal cavity pressure and the deviatoric stress decreases rapidly once the peak values are reached, which can be attributed to the release of accumulated stress due to the degradation and crushing of the structure; (3) The present solution can evolve to the existing solution in reconstituted soils when the structural parameters are sufficiently small. Therefore, the present solution can be regarded as a unified solution for cavity expansion in structured and reconstituted soils; (4) Compared with the previous solutions, the present solution takes the destructuration into account and captures the softening behavior of natural soils well; (5) The simulations of the casing installation and in situ self-boring pressuremeter tests indicate that the present solution provides an effective theoretical tool for the analyses of practical geotechnical problems involving cylindrical cavity expansion.