#### Abstract

Based on the axisymmetric large-strain consolidation (ALSC) model with the void ratio as the variable under equal strain condition, difference schemes of model’s equation, initial condition, and boundary condition were given. Taking phosphatic clay in Florida as a research object, the consolidation behaviors of soil with high water content by axisymmetric large-strain theory and one-dimensional large-strain theory were analyzed. The effect of different kinds of consolidation theories and self-weight stress on an average degree of consolidation was evaluated. The development of the void ratio and excess pore water pressure along the soil layer was clarified. The results show that the theoretical value of Terzaghi’s consolidation degree is always less than that of ALSC (, the average degree of consolidation defined by strain)-vertical drainage in the consolidation process. Terzaghi’s solution overestimates the dissipation rate of excess pore water pressure during the earlier consolidation period but underestimates it during the later consolidation period. The degree of consolidation calculated by Hansbo develops faster than ALSC (, the average degree of consolidation defined by stress)-radial drainage, but slower than ALSC ()-radial drainage. In the ALSC model, is always been faster than . The effect of self-weight on the consolidation degree of axisymmetric large-strain consolidation theory is relatively small (maximum error is less than 16%), while it can accelerate the consolidation rate of soil in one-dimensional large-strain consolidation theory largely. When only the vertical drainage occurs, the consolidation rate in the middle of the soil is obviously lagging the upper and lower parts, while the radial drainage can reduce the void ratio and the excess pore water pressure along the soil layer uniformly and more rapidly.

#### 1. Introduction

In dredging projects, such as port and dock construction projects, the treatment of large-area soft soils is often encountered. Due to the low permeability of soft soils, the consolidation rate under self-weight and applied loads is slow, and prefabricate vertical drains (such as plastic drainage boards) are often used in the site works to speed up the consolidation process of soil. Consolidation with the vertical drainage body in soft soil can be simplified as a single well consolidation problem. Axisymmetric consolidation theory proposed by Barron [1] is widely accepted, and the degree of soil’s consolidation can be predicted by analytical or numerical methods. The consolidation properties of soil can be revealed by this theory. Over the years, based on Barron’s axial symmetry consolidation theory, the in-depth study of soil consolidation behavior was carried out from many aspects, such as well resistance and smearing [2–5], loading changes over time [6–10], nonlinear compression and penetration parameters [11, 12, 13], negative pressure loading [14, 15, 16, 17], and multi-layer [18, 19]. It provides theoretical support for design and calculation of the consolidation method with vertical drains in practical projects.

However, the accuracy of prediction and analysis of soil consolidation behavior is not only related to the parameters used in the model and calculation methods but also essentially related to the consolidation model itself. When using the above-mentioned theoretical models for consolidation analysis of supersoft foundations, especially those with high water ocean soft soils or dredged soils that are in a flow state, there is a big difference between the theoretical and the field values [20]. This is mainly because these theoretical models are all based on the assumption of small strain, without considering the effect of large geometric strains.

Fox et al. [21] used a series of small-strain calculations to perform large-strain radial consolidation analysis by displacement-corrected coordinates. However, the variation of material nonlinearity with geometric large deformation was not considered. Indraratna et al. [22] analyzed the large-strain consolidation behavior considering radial nonlinear flow and non-Darcy’s law under vacuum loading without considering vertical seepage. Jiang et al. [20] coupled geometric nonlinearity and material nonlinearity, considering vertical and radial simultaneous seepage, and expanded Gibson’s one-dimensional large-strain consolidation theory to axisymmetric large-strain consolidation theory. More recently, Zhang et al. [23], Hu et al. [24], and Sun et al. [25] studied the sand-drained large-strain consolidation model with double-layer foundations, creep, and partially penetrated sand wells, respectively. Although the smearing effect controls the consolidation rate of soil [3], its effect on consolidation properties cannot be quantitatively described by these theories.

It should be pointed out that the above-mentioned theory of large-strain consolidation is under free-strain conditions (equal stress consolidation). Some scholars believe that under axisymmetric consolidation, the soil is closer to the state of free strain [7]. From the perspective of practical engineering, some scholars have pointed out that equal settlement is easier to achieve in the project because the horizontal spacing of vertical drainage bodies is the same [21]. The soil consolidation near the sand well is fast, which will cause the contact stress redistribution [26]. It is undeniable that results from the equal strain solutions and the free-strain solutions are very close [1, 3, 7, 27]. The theoretical solution based on the equal strain hypothesis is relatively simple. It can also easily reflect the effect of soil compression characteristics on consolidation in smeared and undisturbed zones. At present, researchers generally assume that the consolidation problem of sand wells is an equal vertical strain problem [26]. In order to facilitate the comparative analysis of large- and small-strain consolidation theories, Cao et al. [28, 29] established the axisymmetric large-strain consolidation model in positive loading pressure (named ALSC) and negative vacuum pressure (named NALSC) with the void ratio as the variable under equal strain condition and verified the validity of these models.

The axisymmetric large-strain consolidation theory breaks through the small-strain theoretical framework and can accurately analyze the settlement characteristics of soils. However, the focus of researchers is mainly on verifying the theory through numerical or analytical solutions or on the development of consolidation rate [20, 23, 24, 25, 28, 29]. Very few studies investigated the difference between axisymmetric large-strain and small-strain theories considering self-weight stress. Also, the effect of gravity stress and radial drainage on the large-strain consolidation characteristics was not adequately investigated, especially for under-consolidated dredged soils with high water content.

In this study, based on the ALSC model, difference schemes of the consolidation equation, initial conditions, and boundary conditions were included. Taking the phosphatic clay in Florida as a research object [30], vertical drainage bodies were arranged based on scenario C (quiescent consolidation and surcharge loading of a pond having a uniform initial void ratio) to investigate the consolidation behavior of soils with high water content. The difference between the calculated values of large- and small-strain under axial symmetry and one-dimensional consolidation theory was analyzed. The influence of self-weight on the average degree of consolidation was considered. The variation of the void ratio and excess pore water pressure along soil height and consolidation time under different drainage modes was studied.

#### 2. ALSC Model

Figure 1 shows the schematic representation of the case of a circular soil cylinder where a vertical drain is surrounded by a smeared zone and undisturbed soil. The soil has an initial layer thickness , and it is freely draining at the top and impermeable at the bottom. The center of the cylinder contains a vertical drain well of radius surrounded by a zone of remolded soil by radius . The smear zone is surrounded by undisturbed soil with a radius of influence . Let and represent the horizontal permeability of the smear zone and undisturbed zone, respectively, and represent the vertical permeability of soil.

The convective coordinate measured downwards in the direction of gravity is the same with one-dimensional large-strain consolidation. The Euler polar coordinates () are used for radical flows, where is measured away from the drain (Figure 2). and represent the vertical pore water inflow and outflow of the unit cell in the soil at a unit time, respectively. and represent the radial pore water inflow and outflow of the unit cell in the soil at a unit time, respectively.

Coupled flow equation, continuous equation of the saturated soil, stress balance equation, effective stress principle, etc. [28] established the ALSC model with the void ratio as variable under equal strain station. The model accounts for smear effect, soil self-weights, radial and vertical flows, and variable permeability and compressibility during the consolidation progress. Well resistance is not considered. The consolidation equation is expressed as shown in the following equation:where is the specific gravity of the solid particles, is the unit weight of water, is the loading pressure; is the ratio , is the ratio , is the effective vertical stress, is the void ratio, is the spatial vertical coordinates in the Lagrange coordinate system, and is the diameter of influence. *µ* is expressed as follows:

Equation (1) is based on the following assumptions:(1)The soil is completely homogeneous and water-saturated.(2)Solid particles and water are incompressible, and the deformation of the soil is completely caused by the discharge of pore water.(3)The vertical and radial flows obey Darcy’s law, and the coefficient of permeability changes with the void ratio during the consolidation.(4)All compressive strains within the soil occur in a vertical direction. The soil particles do not move along the radial and tangential directions, and no creep is considered.(5)Horizontal sections remain horizontal during the consolidation.(6)All vertical loads are applied instantaneously, and the load distribution is uniform over the whole cylindrical area.

The above assumptions remove the limits of small-strain hypothesis of soil. The vertical nonlinear compression characteristics, nonlinear permeability for radial and vertical flows, and skeleton deformation generated by vertical movement of soil particles are considered in model ALSC. Therefore, geometric nonlinearity and material nonlinearity are brought into equal strain consolidation for large-area loading condition, which agrees well with the settlement law with vertical drain in real engineering problems. Equation (1) can be rewritten into the one-dimensional large-strain consolidation equation as given by Gibson et al. [31, 32] by ignoring the horizontal radial flow (i.e., ).

If the assumption of small strain is considered, and there is no vertical flow (i.e. ), Equation (1) can be rewritten into the axisymmetric consolidation control equation as given by Hansbo [3].

#### 3. Difference Schemes of the ALSC Model

As it can be seen, the theoretical model of formula (1) is a highly complex nonlinear equation, and it is difficult to solve it analytically. The finite difference method can be used for calculation and analysis. For convenience, Equation (1) can be written aswhere,

Each of Equations (4)–(7) has its own unique parameters and perform different consolidation behavior. and constitute the gravity item, is relevant to the radial drainage, and represents the coefficient of one-dimensional large-strain consolidation.

Equation (3) is a variable coefficient convection-diffusion equation, , , , and all have upper and lower limits, such as , , , and .

Let represent the space step and time step, respectively, and . The difference scheme is like the convection-diffusion equation with constant coefficients. According to the upwind difference scheme, Equation (3) can be written in the following format, and the discrete schematic diagram of space and time is shown in Figure 3:where

Difference equations can better reflect the material nonlinear characteristics of soil, which is much closer to the actual situation than special assumptions of soil permeability and compressibility for obtaining analytical solutions easily. During the calculation, good results for any given form of will be obtained by substituting , , , , into , , , .

#### 4. Difference Schemes for Initial and Boundary Conditions

Assume that the initial thickness of the high water-saturated homogeneous soil is (thickness in the spatial vertical coordinates in the Lagrange coordinate system, ). The initial void ratio at any position of soil for time zero is the same. Initial condition is described as

The difference scheme is given as

The average void ratio can be obtained from the relationship of effective vertical stress with the void ratio. The top of the soil is freely draining, and the excess pore water pressure dissipated instantaneously. According to the compression properties of soil, the final void ratio at the soil surface is given by

The difference scheme for the upper boundary condition is given as

If the soil at the bottom is impermeable during the consolidation, the bottom boundary condition is determined as

The difference scheme for the bottom boundary condition is given as

#### 5. Case Analysis

According to the ALSC model and the difference structure, the axisymmetric large-strain consolidation calculation program is compiled. Townsend and McVay [30] have described the predictions of ponds’ (solids) elevation histories and one-year pore water pressure and void ratio profiles for four different wasted clay disposal scenarios (Scenario A: quiescent consolidation, uniform initial void ratio; Scenario B: stage filling, nonuniform initial void ratio; Scenario C: quiescent consolidation and surcharge loading of a pond having a uniform initial void ratio; and Scenario D: two-layer quiescent consolidation, sand/clay surcharge, and nonuniform initial void ratio). Station for Scenario C is closer to the assumptions (1) and (6) of the ALSC model, so Scenario C is selected for analysis. The nonlinear parameters, i.e., Equations (16) and (17), are used for prediction by nine different large-strain consolidation models. However, the predictions are varying from program to program:

Figure 4 shows the quiescent consolidation with the top freely draining and impermeable at the bottom. The sludge is water-saturated and homogeneous. A 7.2 m deep waste pond with a uniform initial void ratio of 14.8 is capped with a 9.48 kPa surcharge. The specific gravity of the solid particles is 2.82. To illustrate the difference between large- and small-strain consolidations with vertical drains, the PVD is quincunx arranged with the spacing of 1.0 m based on Scenario C. The width () ∗ thickness () of the PVDs is 100 mm ∗ 4 mm. So, the equivalent radius of the vertical drain (36.3 mm) can be calculated according to the formula proposed by Hansbo [33]. The radius of influence is set to 525 mm. The large-strain material nonlinear parameters change with the void ratio seen in Equations (16) and (17), but the parameters of small-strain consolidation are the constant values. According to Equations (16) and (17), compressibility shall be obtained at each level of the loading, and the vertical coefficient of consolidation can be calculated from its definition. Assuming the soil is isotropy, we found the relationship . The smearing is not considered during the simulation progress. The specific calculation parameters are shown in Table 1.

In the calculation program of the ALSC model, the soil is divided into 100 layers along the depth, and time step length is 0.1 day. Fifteen different kinds of simulations were performed, such as consolidation with or without self-weight, and only radial or vertical flows occurs. The whole computation time is 10000 days except for scenario “only vertical drainage without self-weight”, where the time is 30000 days. The calculation parameters were obtained from [30]. The unit conversion is as follows: 1 ft = 304.8 mm and 1 psf = 47.88 Pa.

In small-strain calculations, the time step is 1.0 day, the computation time of models “Hansbo” (which states that the axisymmetric small-strain consolidation theory of only radial flow is considered) and “Hansbo + Terzaghi” (which refers to that the vertical consolidation degree is calculated by Terzaghi’s theory, the radial consolidation degree is calculated by Hansbo’s theory, and then the average degree of consolidation is calculated according to Carrillo’s theorem) is 10000 days. However, the computation time of model “Terzaghi” (which means computing by the Terzaghi classical one-dimensional small-strain consolidation theory—vertical flows only) needs 70000 days to achieve stablity.

#### 6. Differences of Average Degree of Consolidation between Large-Strain and Small-Strain Theories

Figure 5 shows the average consolidation degree versus time for large- and small-strain consolidation theories. “Double drainage” represents the simultaneous occurrence of vertical and radial flows, “radial drainage” indicates that only radial flow occurs (, when calculating), and “vertical drainage” means that only vertical flows occur (, when computing). In the calculation process, self-weight stress is considered, i.e., . At the same time, to reflect the difference between the settlement rate and the dissipating rate of excess pore water, the stress consolidation degree (defined by stress) and the strain consolidation degree (defined by strain) were calculated separately.(1)The difference between one-dimensional consolidations. During the early stage of consolidation (within 1000 days), the degree of consolidation for Terzaghi’s one-dimensional small-strain consolidation theory is faster than that for ALSC ()-vertical drainage. This is because the consolidation coefficient remains unchanged in Terzaghi theory. After 1000 days, the development of the Terzaghi consolidation degree is obviously lagging and is always less than that of ALSC ()-vertical drainage in the same consolidation time. This is because the large-strain theory considers the reduction in drainage distance during the consolidation process [34]. Terzaghi’s solution overestimates the dissipation rate of the excess pore water pressure in the early consolidation stage and underestimates it during the later period. We can see from the consolidation stability time, the Terzaghi solution clearly overestimates the settling time of the soil. Take the degree of consolidation of 80% as an example. Terzaghi’s theory needs 12853 days, while the time required for ALSC ()-vertical drainage and ALSC ()-vertical drainage are 876 days and 3247 days, and the consolidation time is 15 and 4 times those of the ALSC theory, respectively.(2)The difference of radial consolidation. The development of consolidation by Hansbo is always faster than ALSC ()-radial drainage consolidation. This is because Hansbo’s radial permeability coefficient is constant, while the ALSC ()-radial permeability coefficient decreases with time. When the degree of consolidation is less than 90%, Hansbo’s solution is smaller than ALSC ()-radial drainage value; after 90% consolidation degree, the two are closer. This is mainly because ALSC () considers the vertical geometrical nonlinearities, and settlement changes faster before 90% consolidation. In the consolidation stability time, taking the degree of consolidation of 80% as an example, Hansbo’s small-strain consolidation time requires 47 days, while the time required for and calculated by the ALSC model considering only radial flows are 27 days and 234 days, respectively. The large-strain consolidation time is 1/1.7 times and 5 times that of the small-strain theory, respectively.(3)The degree of consolidation for ALSC-double drainage in this case is like that for ALSC-radial drainage. However, the degree of consolidation for ALSC-vertical drainage is much smaller than that for ALSC-radial drainage. It shows that when laying vertical drainage boards for working condition C, drainage path is shortened significantly, and radial soil consolidation (flows) occurs mainly in the soil. “Hansbo” solution is consistent with “Hansbo+Terzaghi” solution, and “Terzaghi” solution is much smaller than “Hansbo” solution, which also illustrates this point.(4)In small-strain consolidation theory, the degrees of consolidation defined by stress and by strain are the same due to the assumption of microdeformation. In the ALSC theory, material nonlinearity and geometric nonlinearity are considered. In the same consolidation time, development of has always been faster than [20]. Thus, the dissipation rate of excess pore water pressure always lags the settlement deformation rate. If the settlement is used as the guidance in the construction process, it will overestimate the shear strength of the soil.

#### 7. Effect of Self-Weight on Large-Strain Consolidation Behavior under the ALSC Model

To discuss the effect of self-weight on large-strain consolidation behavior, let (self-weight is not considered) and (self-weight is considered). Figure 6 shows the effect of self-weight on the degree of consolidation of double drainage and radial flows only. It can be seen that, in the case of and , development of the average consolidation degree with double drainage is basically similar to that with radial flows only. Therefore, the time-course curve of coupled radial and vertical flows is not drawn in Figure 7 in order to clearly reflect the effect of self-weight on the degree of consolidation.

Figure 7 shows that when only radial flows occur, compared with the consolidation rate when , the consolidation rate of when is faster. The error between the two is 16% on the first day and is less than 10% on the fourth day. The consolidation curve basically coincided after 100 days. The curves, comparing the situation where self-weight is considered versus the situation where it is not, basically coincided.

For large-strain consolidation where only vertical flows occur, the degree of consolidation of and with considering self-weight is faster than those without considering it. Therefore, for this kind of unconsolidated soil with high water content, the self-weight stress has little effect on the radial consolidation degree. However, self-weight has some effects on the degree of vertical consolidation. If the self-weight of the soil is neglected, the settlement rate of the soil and the dissipation rate of excess pore water pressure are underestimated [35, 36].

#### 8. Change Law of Void Ratio and Excess Pore Water Pressure with Soil Height and Consolidation Time

The variations of the void ratio and excess pore water pressure along the soil height are shown in Figures 8 and 9. As the consolidation time increases, the void ratio gradually decreases. When only vertical flows occur (Figure 8), the change rate of the void ratio varies greatly at different depths. Under the effect of external load and self-weight, the consolidation rate in the middle of the soil lags than that of the upper and lower parts [35]. When only radial flows occur, the void ratio is uniform along the height. In the same consolidation time, the reduction rate of the void ratio for radial flows only is faster than that of the vertical flows only, indicating that the arrangement of vertical drainage body can shorten the drainage path and accelerate soil consolidation rate.

From the change law of excess pore water pressure, it can be seen that, after the arrangement of the vertical drainage body, although the excess pore pressure of the soil at different depths is evenly dissipated, excess pore pressure along the depth is linearly increasing which is affected by the self-weight when only radial flows occur. However, distribution of the excess pore water pressure along the depth shows nonlinearity when only vertical flows occur. Like the change law of the void ratio, the dissipation rate of the excess pore pressure in the middle is slower than those in the upper and lower parts. From this point of view, radial flows significantly accelerate the dissipating rate of excess pore water pressure. For example, when the consolidation time is 361 days, the excess pore water pressure at the bottom of the soil when only radial flows occur is 1.5 kPa but when only vertical flows occur is approximately 15 kPa (10 times larger than the values of the previous case).

Table 2 summarizes the pertinent data and statistics of nine research institutes [30] and model ALSC (vertical flows only) for Scenario C. The final average thickness of soil and other consolidation parameters after one year of nine research institutes are close to those of the model ALSC, which confirms that it is a reasonable model and good precision for soil or dredged sludge at high water content.

#### 9. Conclusions

In this paper, based on the ALSC model under equal strain conditions (established by Cao et al. [28]), the consolidation behavior of soils with high water content considering self-weight was investigated, and the following conclusions are obtained:(1)The theoretical values obtained from Terzaghi’s one-dimensional small-strain consolidation degree are always less than those obtained from ALSC ()-vertical drainage. Terzaghi’s solution overestimates the dissipation rate of excess pore water pressure in the earlier consolidation stage and underestimates it during the later consolidation stage. Terzaghi’s solution overestimates the consolidation time of soil.(2)In the axisymmetric large-strain system and the small-strain theory system, the vertical drainage body can shorten the drainage path, and the consolidation of soils with high water content is dominated by radial flows. At the same consolidation time, the average degree of consolidation obtained from Hansbo’s theory is smaller than the values of ALSC()-radial drainage, but it is higher than values obtained from ALSC()-radial drainage. In the ALSC model, development of is always faster than .(3)For the unconsolidated soil with high water content, the effect of self-weight stress on axisymmetric large-strain radial consolidation is small, but on one-dimensional large-strain vertical consolidation, it is large.(4)In one-dimensional large-strain consolidation, the consolidation rate in the middle of the soil lags than those in the upper and lower parts under the combined action of external load and self-weight. The vertical drainage body can evenly dissipate the excess pore water pressure and accelerate the consolidation rate. Affected by the self-weight, the void ratio and the excess pore water pressure are distributed linearly along the soil. The calculation results based on the ALSC model are close to the average statistics of the nine institutions, which verifies the accuracy and applicability of the ALSC model and the difference solution method.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This research was supported by the Natural National Science Foundation of China (Grant no. 51608312) and Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents (Grant no. 2016RCJJ021).