# Nonlinear Radial Consolidation Analysis of Soft Soil with Vertical Drains under Cyclic Loadings

Authors: Pyol Kim; Hyong-Sik Kim; Yong-Gun Kim; Chung-Hyok Paek; Song-Nam Oh; Sung-Ryol So; Chang-Son Ok
Journal: Shock and Vibration
 (2020)

---

## Abstract

This paper presents analytical solutions for nonlinear radial consolidation of soft soil with vertical drains under various cyclic loadings. By considering the nonlinear variations of compressibility and permeability expressed by the logarithm relations (e−logσ′ and e−logkh), the governing equation for nonlinear radial consolidation of the soil under equal strain and time-dependent loading is established. The analytical solutions are derived for nonlinear radial consolidation under haversine cyclic loading, trapezoidal cyclic loading, rectangular cyclic loading, and triangular cyclic loading. The presented solution is verified through the degeneration into the existing solutions for nonlinear radial consolidation under constant and ramp loadings, which shows the solution proposed in this paper is more general for nonlinear radial consolidation under time-dependent loading. The nonlinear radial consolidation behavior of the soil with vertical drains subjected to various cyclic loadings is investigated using the solutions developed. The proposed solutions can be effectively utilized in the analysis of nonlinear radial consolidation under various cyclic loadings.


---

## Body

## 1. Introduction

In most coastal and marine areas, the soils beneath many structures are highly compressible and may be subjected to complicated cyclic loadings, such as filling and discharging in silos, tanks, and reservoirs, traffic loads on embankments, highways, railways, or airports, wave and wind actions for offshore structures, pile construction, and machine vibrations. The consolidation behavior of the soft soil under cyclic loadings is very complex and has been studied for many years. In the past decades, a lot of scholars have conducted research on one-dimensional (1D) consolidation of soft soil subjected to cyclic loadings based on either the linear consolidation theory [1–6] or the nonlinear consolidation theory [7–16]. It is worth mentioning that Razouki and Schanz [4] and Razouki et al. [5] studied 1D consolidation of a clay layer subjected to haversine cyclic loading with and without the rest period on the basis of linear assumptions. More recently, Kim et al. [16] derived analytical solutions for 1D nonlinear consolidation of a saturated soil with variable compressibility and permeability under various cyclic loadings. On the contrary, Kim et al. [17] presented analytical solutions for 1D consolidation of unsaturated soils under cyclic loadings.Vertical drains are widely used for improvement of soft soil, and the radial consolidation behavior of the soil with vertical drains under cyclic loadings has also been investigated by many researchers. Indraratna et al. [18] conducted a large-scale triaxial test to examine the effectiveness of a prefabricated vertical drain (PVD) installed in soft clay subjected to cyclic loads, while Ni et al. [19, 20] proposed new constitutive models to predict the radial consolidation of soft soils under cyclic loading. Razouki [21] studied radial consolidation behavior of clay soil with vertical drains under free-strain condition and haversine cyclic loading based on the assumption that the consolidation coefficient in the radial direction is constant during the consolidation process.Although there are many studies on 1D consolidation of the soil under cyclic loadings, research studies on radial consolidation under cyclic loadings are very limited, and any analytical solutions have hardly ever been proposed for nonlinear radial consolidation of the soil with vertical drains subjected to cyclic loadings with the consideration of variable compressibility and permeability. Moreover, the effects of different parameters on nonlinear radial consolidation under cyclic loadings have not been investigated. Therefore, this study is aimed to derive analytical solutions for nonlinear radial consolidation under various cyclic loadings and investigate nonlinear radial consolidation behavior of the soil under various cyclic loadings.This paper is organized as follows. Firstly, the mathematical model and the governing equation for nonlinear radial consolidation of the soil with vertical drains subjected to time-dependent loading are presented in Section2. Secondly, analytical solutions for nonlinear radial consolidation under cyclic loadings are derived in Section 3, and the proposed solutions are verified through the degeneration into special cases for nonlinear radial consolidation in Section 4. Thirdly, the effects of different parameters on nonlinear radial consolidation behavior of the soft soil with vertical drains under various cyclic loadings are investigated in Section 5. Finally, conclusions are given in Section 6.


## 2. Mathematical Model and Governing Equation

The mathematical model for radial consolidation analysis of soft soil with vertical drains is shown in Figure1. Figures 1(a) and 1(b) show the arrangement of vertical drains and the cylindrical unit cell of the soil with a vertical drain, where rw,rs, and re are the radii of the vertical drain, the smear zone, and the influence zone, respectively. The parameters s,ks, and kh are the spacing of vertical drains and the horizontal permeability coefficients of the soil within the smear zone and the influence zone. qt is the uniformly distributed time-dependent loading applied on the top surface of the soil of thickness H. In order to obtain the analytical solution for nonlinear radial consolidation of the soil with vertical drains subjected to time-dependent loading, the main assumptions are summarized as follows:(1)
The soil layer is saturated and homogeneous.(2)
The soil particles and pore water are incompressible.(3)
Only the radial flow of pore water is considered. Darcy’s law is valid for the water flow.(4)
Both the loading and the settlement occur only in the vertical direction, while the soil and the vertical drain have an equal strain at any depth (i.e., the equal strain condition).(5)
The effect of the drain resistance is neglected.(6)
The nonlinear variations of compressibility and permeability during the consolidation process are given by(1)e=e0−Cclogσ′σ0′,(2)e=e0+Cklogkhkh0,where e and e0 are the void ratio at any time and the initial void ratio, respectively; σ′ and σ0′ are the effective stress at any time and the initial effective stress, respectively; kh and kh0 are the horizontal permeability coefficients of the soil within the influence zone at any time and at the initial time, respectively; Cc and Ck are the compression index and the permeability index, respectively.Figure 1
A schematic diagram for radial consolidation analysis. (a) The arrangement of vertical drains; (b) the cylindrical unit cell of the soil with a vertical drain.Based on equations (1) and (2), the following relations are obtained:(3)kh=kh0σ0′σ′Cc/Ck,(4)mv=−11+e0∂e∂σ′=11+e0Ccln101σ′=mv0σ0′σ′,where mv and mv0 are the compressibility coefficients of the soil at any time and at the initial time, respectively.The governing equation for radial consolidation of the soil under the equal strain condition can be written as follows [22]:(5)2rrs2−r2⋅ksγw⋅∂us∂r=∂ε∂t,rw<r<rs,2rre2−r2⋅khγw⋅∂un∂r=∂ε∂t,rs<r<re,where γw is the unit weight of water; us and un are the excess pore water pressure within the smear zone and the influence zone, respectively; ε is the vertical strain of the soil and the vertical drain; r is the radial distance from the center of the vertical drain; t is the time.The average excess pore water pressure for both the smear zone and the influence zone,u¯, can be defined as(6)u¯=1πre2−rw2∫rwrsus2πrdr+∫rsreun2πrdr.The boundary conditions in the radial direction are given by(7)us=0,r=rw,us=un,r=rs,ks∂us∂r=kh∂un∂r,r=rs,∂un∂r=0,r=re.The initial condition for the average excess pore water pressure can be expressed as(8)u¯r,0=q0=u¯0.


## 3. Analytical Solutions for Cyclic Loadings

In this section, analytical solutions are derived for nonlinear radial consolidation of the soil with vertical drains under various cyclic loadings (Figure2), such as haversine cyclic loading, trapezoidal cyclic loading, rectangular cyclic loading, and triangular cyclic loading.Figure 2
Various cyclic loadings. (a) Haversine cyclic loading; (b) trapezoidal cyclic loading; (c) rectangular cyclic loading; (d) triangular cyclic loading.
(a)(b)(c)(d)


### 3.1. Average Excess Pore Water Pressure

The excess pore water pressures within the smear zone and the influence zone can be obtained by integrating (5) and using the boundary conditions (7) as follows:(9)us=γw2ks∂ε∂tre2lnrrw−r2−rw22,rw<r<rs,un=γw2kh∂ε∂tre2lnrrs−r2−rs22+γw2ks∂ε∂tre2lnrsrw−rs2−rw22,rs<r<re.Substituting (9) into (6), the average excess pore water pressure can be determined by(10)u¯=re2ξnγw2kh1+e0∂e∂σ′∂q∂t−∂u¯∂t,where(11)∂ε∂t=11+e0∂e∂t=11+e0∂e∂σ′∂σ′∂t=11+e0∂e∂σ′∂q∂t−∂u¯∂t,and ξn is a parameter that reflects the characteristics of the smear zone and the influence zone of soil with vertical drains derived as follows:(12)ξn=n2n2−1lnn−34+4n2−14n2+khks−1lnm+1−m2n2+m4−14n4,n=rerw,m=rsrw.Defining the time factor in the radial directionTh=ch0t/4re2, in which ch0=kh01+e0/γw∂e/∂σ′t=0 is the initial horizontal consolidation coefficient, and combining with equations (3), (4), and (10), the following expression can be obtained:(13)du¯dTh=−8ξnηu¯+dqdTh,where η=1+q−u¯/σ0′1−Cc/Ck.Since (13) is a nonlinear ordinary differential equation, it is difficult to obtain its exact solution. η varies from 1.0 to 1+qu/σ0′1−Cc/Ck, in which qu is the ultimate value of loading q. Hence, η in (13) can be assumed to be constant and have an average value given by η¯=1/21+1+qu/σ0′1−Cc/Ck [16, 22, 23]. Then, (13) can be linearized to(14)du¯dTh=−8ξnη¯u¯+dqdTh.The general solution to (14) satisfying the initial condition of (8) can be given as follows:(15)u¯=exp−8ξnη¯Thu¯0+∫0ThdqdThexp8ξnη¯τdτ.


#### 3.1.1. Haversine Cyclic Loading

Haversine cyclic loading shown in Figure2(a) can be expressed as follows:(16)qt=qusin2πtt0,where t0 is the period of one loading cycle.Substituting (16) into (15), the complete solution for the average excess pore water pressure under haversine cyclic loading can be obtained as follows:(17)u¯=πquλ2T02+4π2λT0sin2πThT0−2πcos2πThT0−exp−λTh,where λ=8/ξn⋅η¯=4/ξn1+Nσ1−Cc/Ck,Nσ=σ0′+qu/σ′,andT0=ch0t0/4re2.


#### 3.1.2. Trapezoidal Cyclic Loading

Trapezoidal cyclic loading shown in Figure2(b) can be expressed as follows:(18)qt=quαt0t−N−1βt0,N−1βt0≤t≤N−1β+αt0,qu,N−1β+αt0≤t≤N−1β+1−αt0,−quαt0t−N−1βt0−t0,N−1β+1−αt0≤t≤N−1β+1t0,0,N−1β+1t0≤t≤Nβt0,where βt0 is the period of one loading cycle; α and β are the loading parameters corresponding to the rate of loading increment or decrement and the rest period of loading, respectively; N is the cycle number.Similarly, substituting (18) into (17), the average excess pore water pressure for the trapezoidal cyclic loading can be written as(19)u¯=quC1λTcexp−λTh,N−1βt0≤t≤N−1β+αt0,quC2λTcexp−λTh,N−1β+αt0≤t≤N−1β+1−αt0,quC3λTcexp−λTh,N−1β+1−αt0≤t≤N−1β+1t0,quC4λTcexp−λTh,N−1β+1t0≤t≤Nβt0,where(20)C1=DN−1+D1,C2=DN−1+D2,C3=DN−1+D2+D3,C4=DN,Dn=∑N=1nD2+D4,D1=expλTh−expλTb,D2=expλTbexpλTc−1,D3=expλTf−Tc−expλTh,D4=expλTfexp−λTc−1,Tc=ch0αt04re2,Tb=ch0N−1βt04re2,Tf=ch0N−1β+1t04re2.


#### 3.1.3. Rectangular and Triangular Cyclic Loadings

Trapezoidal cyclic loading can be degenerated into rectangular cyclic loading and triangular cyclic loading by changing the value of the loading parameterα. When α=0, trapezoidal cyclic loading reduces into rectangular cyclic loading (Figure 2(c)). On the other hand, trapezoidal cyclic loading reduces to triangular cyclic loading (Figure 2(d)) when α=0.5.Thus, analytical solutions for the average excess pore water pressure under rectangular and triangular cyclic loadings can be obtained substitutingα=0 and α=0.5 into (19), respectively.


### 3.2. Average Degree of Consolidation

Based on the average excess pore water pressures obtained under the above different cyclic loadings, the average degrees of consolidation with respect to excess pore water pressure and settlement,Up and Us, can be expressed, respectively, as follows:(21)Up=q−u¯qu,(22)Us=∫0Hεdz∫0Hεfdz=a∫0Hlogσ′/σ0′dz∫0Hlogσf′/σ0′logdz=lnσ0′+q−u¯/σ0′lnNσ,where εf and σf′ are the final vertical strain and the final effective stress, respectively.


## 3.1. Average Excess Pore Water Pressure

The excess pore water pressures within the smear zone and the influence zone can be obtained by integrating (5) and using the boundary conditions (7) as follows:(9)us=γw2ks∂ε∂tre2lnrrw−r2−rw22,rw<r<rs,un=γw2kh∂ε∂tre2lnrrs−r2−rs22+γw2ks∂ε∂tre2lnrsrw−rs2−rw22,rs<r<re.Substituting (9) into (6), the average excess pore water pressure can be determined by(10)u¯=re2ξnγw2kh1+e0∂e∂σ′∂q∂t−∂u¯∂t,where(11)∂ε∂t=11+e0∂e∂t=11+e0∂e∂σ′∂σ′∂t=11+e0∂e∂σ′∂q∂t−∂u¯∂t,and ξn is a parameter that reflects the characteristics of the smear zone and the influence zone of soil with vertical drains derived as follows:(12)ξn=n2n2−1lnn−34+4n2−14n2+khks−1lnm+1−m2n2+m4−14n4,n=rerw,m=rsrw.Defining the time factor in the radial directionTh=ch0t/4re2, in which ch0=kh01+e0/γw∂e/∂σ′t=0 is the initial horizontal consolidation coefficient, and combining with equations (3), (4), and (10), the following expression can be obtained:(13)du¯dTh=−8ξnηu¯+dqdTh,where η=1+q−u¯/σ0′1−Cc/Ck.Since (13) is a nonlinear ordinary differential equation, it is difficult to obtain its exact solution. η varies from 1.0 to 1+qu/σ0′1−Cc/Ck, in which qu is the ultimate value of loading q. Hence, η in (13) can be assumed to be constant and have an average value given by η¯=1/21+1+qu/σ0′1−Cc/Ck [16, 22, 23]. Then, (13) can be linearized to(14)du¯dTh=−8ξnη¯u¯+dqdTh.The general solution to (14) satisfying the initial condition of (8) can be given as follows:(15)u¯=exp−8ξnη¯Thu¯0+∫0ThdqdThexp8ξnη¯τdτ.


### 3.1.1. Haversine Cyclic Loading

Haversine cyclic loading shown in Figure2(a) can be expressed as follows:(16)qt=qusin2πtt0,where t0 is the period of one loading cycle.Substituting (16) into (15), the complete solution for the average excess pore water pressure under haversine cyclic loading can be obtained as follows:(17)u¯=πquλ2T02+4π2λT0sin2πThT0−2πcos2πThT0−exp−λTh,where λ=8/ξn⋅η¯=4/ξn1+Nσ1−Cc/Ck,Nσ=σ0′+qu/σ′,andT0=ch0t0/4re2.


### 3.1.2. Trapezoidal Cyclic Loading

Trapezoidal cyclic loading shown in Figure2(b) can be expressed as follows:(18)qt=quαt0t−N−1βt0,N−1βt0≤t≤N−1β+αt0,qu,N−1β+αt0≤t≤N−1β+1−αt0,−quαt0t−N−1βt0−t0,N−1β+1−αt0≤t≤N−1β+1t0,0,N−1β+1t0≤t≤Nβt0,where βt0 is the period of one loading cycle; α and β are the loading parameters corresponding to the rate of loading increment or decrement and the rest period of loading, respectively; N is the cycle number.Similarly, substituting (18) into (17), the average excess pore water pressure for the trapezoidal cyclic loading can be written as(19)u¯=quC1λTcexp−λTh,N−1βt0≤t≤N−1β+αt0,quC2λTcexp−λTh,N−1β+αt0≤t≤N−1β+1−αt0,quC3λTcexp−λTh,N−1β+1−αt0≤t≤N−1β+1t0,quC4λTcexp−λTh,N−1β+1t0≤t≤Nβt0,where(20)C1=DN−1+D1,C2=DN−1+D2,C3=DN−1+D2+D3,C4=DN,Dn=∑N=1nD2+D4,D1=expλTh−expλTb,D2=expλTbexpλTc−1,D3=expλTf−Tc−expλTh,D4=expλTfexp−λTc−1,Tc=ch0αt04re2,Tb=ch0N−1βt04re2,Tf=ch0N−1β+1t04re2.


### 3.1.3. Rectangular and Triangular Cyclic Loadings

Trapezoidal cyclic loading can be degenerated into rectangular cyclic loading and triangular cyclic loading by changing the value of the loading parameterα. When α=0, trapezoidal cyclic loading reduces into rectangular cyclic loading (Figure 2(c)). On the other hand, trapezoidal cyclic loading reduces to triangular cyclic loading (Figure 2(d)) when α=0.5.Thus, analytical solutions for the average excess pore water pressure under rectangular and triangular cyclic loadings can be obtained substitutingα=0 and α=0.5 into (19), respectively.


## 3.1.1. Haversine Cyclic Loading

Haversine cyclic loading shown in Figure2(a) can be expressed as follows:(16)qt=qusin2πtt0,where t0 is the period of one loading cycle.Substituting (16) into (15), the complete solution for the average excess pore water pressure under haversine cyclic loading can be obtained as follows:(17)u¯=πquλ2T02+4π2λT0sin2πThT0−2πcos2πThT0−exp−λTh,where λ=8/ξn⋅η¯=4/ξn1+Nσ1−Cc/Ck,Nσ=σ0′+qu/σ′,andT0=ch0t0/4re2.


## 3.1.2. Trapezoidal Cyclic Loading

Trapezoidal cyclic loading shown in Figure2(b) can be expressed as follows:(18)qt=quαt0t−N−1βt0,N−1βt0≤t≤N−1β+αt0,qu,N−1β+αt0≤t≤N−1β+1−αt0,−quαt0t−N−1βt0−t0,N−1β+1−αt0≤t≤N−1β+1t0,0,N−1β+1t0≤t≤Nβt0,where βt0 is the period of one loading cycle; α and β are the loading parameters corresponding to the rate of loading increment or decrement and the rest period of loading, respectively; N is the cycle number.Similarly, substituting (18) into (17), the average excess pore water pressure for the trapezoidal cyclic loading can be written as(19)u¯=quC1λTcexp−λTh,N−1βt0≤t≤N−1β+αt0,quC2λTcexp−λTh,N−1β+αt0≤t≤N−1β+1−αt0,quC3λTcexp−λTh,N−1β+1−αt0≤t≤N−1β+1t0,quC4λTcexp−λTh,N−1β+1t0≤t≤Nβt0,where(20)C1=DN−1+D1,C2=DN−1+D2,C3=DN−1+D2+D3,C4=DN,Dn=∑N=1nD2+D4,D1=expλTh−expλTb,D2=expλTbexpλTc−1,D3=expλTf−Tc−expλTh,D4=expλTfexp−λTc−1,Tc=ch0αt04re2,Tb=ch0N−1βt04re2,Tf=ch0N−1β+1t04re2.


## 3.1.3. Rectangular and Triangular Cyclic Loadings

Trapezoidal cyclic loading can be degenerated into rectangular cyclic loading and triangular cyclic loading by changing the value of the loading parameterα. When α=0, trapezoidal cyclic loading reduces into rectangular cyclic loading (Figure 2(c)). On the other hand, trapezoidal cyclic loading reduces to triangular cyclic loading (Figure 2(d)) when α=0.5.Thus, analytical solutions for the average excess pore water pressure under rectangular and triangular cyclic loadings can be obtained substitutingα=0 and α=0.5 into (19), respectively.


## 3.2. Average Degree of Consolidation

Based on the average excess pore water pressures obtained under the above different cyclic loadings, the average degrees of consolidation with respect to excess pore water pressure and settlement,Up and Us, can be expressed, respectively, as follows:(21)Up=q−u¯qu,(22)Us=∫0Hεdz∫0Hεfdz=a∫0Hlogσ′/σ0′dz∫0Hlogσf′/σ0′logdz=lnσ0′+q−u¯/σ0′lnNσ,where εf and σf′ are the final vertical strain and the final effective stress, respectively.


## 4. Degeneration of the Proposed Solution to Special Cases

It can be easily found that the proposed analytical solution for trapezoidal cyclic loading can be degenerated into the existing solutions for nonlinear radial consolidation under constant loading and ramp loading.


### 4.1. Constant Loading

Whenα=0,β=1, and N=1, trapezoidal cyclic loading reduces to constant loading (Figure 3(a)) expressed by qt=qu. Then, the solution for the average excess pore water pressure can be obtained as follows:(23)u¯=quexp−λTh.Figure 3
Constant loading (a) and ramp loading (b).
(a)(b)Equation (23) is the solution for nonlinear radial consolidation of clays with variable permeability and compressibility under constant loading proposed by Indraratna et al. [22].


### 4.2. Ramp Loading

WhenN=1 and t0 tends to infinity, trapezoidal cyclic loading reduces to ramp loading (Figure 3(b)) given by(24)qt=ttcqu,t≤tc,qu,t≥tc,where tc=αt0.The solution can be derived as(25)u¯=quλTcexpλTh−1,t≤tc,quλTcexpλTh1−exp−λTc,t≥tc.Equation (25) is the solution for nonlinear radial consolidation of a clay layer with variable compressibility and permeability under ramp loading developed by Lekha et al. [23]. As shown in Figure 4, it can be seen that the results from the degenerate solutions and the existing solutions are consistent.Figure 4
Comparison with the existing solutions for constant loading (a) and ramp loading (b).
(a)(b)Through the degeneration mentioned above, it can be seen that the solution developed in this paper is a more general one for nonlinear radial consolidation under time-dependent loading.


## 4.1. Constant Loading

Whenα=0,β=1, and N=1, trapezoidal cyclic loading reduces to constant loading (Figure 3(a)) expressed by qt=qu. Then, the solution for the average excess pore water pressure can be obtained as follows:(23)u¯=quexp−λTh.Figure 3
Constant loading (a) and ramp loading (b).
(a)(b)Equation (23) is the solution for nonlinear radial consolidation of clays with variable permeability and compressibility under constant loading proposed by Indraratna et al. [22].


## 4.2. Ramp Loading

WhenN=1 and t0 tends to infinity, trapezoidal cyclic loading reduces to ramp loading (Figure 3(b)) given by(24)qt=ttcqu,t≤tc,qu,t≥tc,where tc=αt0.The solution can be derived as(25)u¯=quλTcexpλTh−1,t≤tc,quλTcexpλTh1−exp−λTc,t≥tc.Equation (25) is the solution for nonlinear radial consolidation of a clay layer with variable compressibility and permeability under ramp loading developed by Lekha et al. [23]. As shown in Figure 4, it can be seen that the results from the degenerate solutions and the existing solutions are consistent.Figure 4
Comparison with the existing solutions for constant loading (a) and ramp loading (b).
(a)(b)Through the degeneration mentioned above, it can be seen that the solution developed in this paper is a more general one for nonlinear radial consolidation under time-dependent loading.


## 5. Consolidation Behavior under Cyclic Loadings

In this section, the effects of different parameters on nonlinear radial consolidation behavior under different types of cyclic loadings are investigated. Since the parameterCc/Ck mainly changes from 0.5 to 2.0 [24], Cc/Ck=0.5 and Cc/Ck=1.5 are taken for the cases of Cc/Ck<1 and Cc/Ck>1, respectively. The other properties adopted in this study are as follows:(26)n=15,m=4,khks=5,Nσ=2.5,T0=0.1,β=1.5.Figure5 shows the variation of the average degree of consolidation Up and Us with time factor Th at different values of Cc/Ck under different types of cyclic loadings. It can be found that average degree of consolidation Up and Us decreases but fluctuates with the increase in Cc/Ck, suggesting that the dissipation rate of excess pore water pressure and the settlement rate decrease with the fluctuation when Cc/Ck increases. Furthermore, a bigger value of Cc/Ck results in a smaller amplitude of fluctuation in the dissipation rate and the settlement rate. In addition, it can be seen that Us is always bigger than Up in all cases, which means the rate of settlement is faster than that of dissipation of excess pore water pressure during the consolidation process.Figure 5
Nonlinear radial consolidation behavior at different values ofCc/Ck under haversine cyclic loading (a) and trapezoidal cyclic loading (b).
(a)(b)Figure6 presents the variation of the average degree of consolidation Up and Us with time factor Th at different values of Nσ under different types of cyclic loadings. It can be found that a bigger value of Nσ results in a faster dissipation rate of excess pore water pressure in the case of Cc/Ck<1, but a smaller value of Nσ leads to a faster dissipation rate in the case of Cc/Ck>1. Meanwhile, the bigger value of Nσ increases the rate of settlement in all cases of Cc/Ck<1 and Cc/Ck>1. Since Nσ is related to the maximum loading qu, the following conclusion is drawn: the dissipation rate increases when maximum loading qu increases in the case of Cc/Ck<1 but decreases with the increase in maximum loading qu in the case of Cc/Ck>1, whereas the settlement rate increases with the increase in maximum loading qu in all cases of Cc/Ck<1 and Cc/Ck>1.Figure 6
Nonlinear radial consolidation behavior at different values ofNσ under haversine cyclic loading (a) and trapezoidal cyclic loading (b).
(a)(b)Figure7 represents the variation of the average degree of consolidation Up and Us with time factor Th under different types of cyclic loadings with different loading parameter α. It can be seen that the smaller the loading parameter α is, the more quickly the dissipation and the settlement develop in all cases of Cc/Ck<1 and Cc/Ck>1. As the loading parameter α reflects the rate of loading increment or decrement, it can be found that the higher the rate of loading increment or decrement is, the faster the rates of dissipation and settlement are. Moreover, the rates of dissipation and settlement under rectangular cyclic loading (i.e., α=0) are the biggest, and the rates of dissipation and settlement under triangular cyclic loading (i.e., α=0.5) are the smallest.Figure 7
Nonlinear radial consolidation behavior under cyclic loadings with different loading parameterα.Figure8 describes the variation of the average degree of consolidation Up and Us with time factor Th under different cyclic loadings with different loading parameter β. As shown in Figure 8(a), through the results under rectangular cyclic loading (α=0) for different loading parameter β, it can be found that the rates of dissipation and settlement increase with the decrease in the value of loading parameter β in all cases of Cc/Ck<1 and Cc/Ck>1. Moreover, rectangular cyclic loading without the rest period (i.e., β=1) reduces into constant loading, and the rates of dissipation and settlement are the highest and increase monotonically with time in that case. The results under trapezoidal cyclic loading and triangular cyclic loading for different loading parameter β can be observed in Figures 8(b) and 8(c), respectively. Similar to the results in Figure 8(a), the dissipation rate and the settlement rate increase with the decrease in the value of loading parameter β in all cases of Cc/Ck<1 and Cc/Ck>1. As well as, a smaller value of the loading parameter β results in more cycles and a smaller amplitude of fluctuation in the rates of dissipation and settlement. Since the loading parameter β reflects the rest period of loading, it is quite obvious that a longer rest period of loading induces a decrease in the rates of dissipation and settlement, less cycles, and a bigger amplitude of fluctuation.Figure 8
Nonlinear radial consolidation behavior under cyclic loadings with different loading parameterβ. Rectangular cyclic loading (a); trapezoidal cyclic loading (b); triangular cyclic loading (c).
(a)(b)(c)Figure9 outlines the variation of the average degree of consolidation Up and Us with time factor Th under different cyclic loadings with different loading parameter T0. It can be seen that a smaller value of the loading parameter T0 induces a smaller amplitude and more cycles of fluctuation in the rates of dissipation and settlement in all cases of Cc/Ck<1 and Cc/Ck>1. In addition, the loading parameter T0 has a significant influence on the amplitude and the cycles of fluctuation in the dissipation rate and the settlement rate but does not affect the mean value of the rates. The loading parameter T0 is related to the period of loading, and thus it can be found that the longer the period of loading, the bigger the amplitude and less cycles of fluctuation in the dissipation rate and the settlement rate.Figure 9
Nonlinear radial consolidation behavior under cyclic loadings with different loading parameterT0: haversine cyclic loading (a); trapezoidal cyclic loading (b).
(a)(b)It should be noted that unlike constant loading and ramp loading conditions, the dissipation process of the excess pore water pressure does not complete after a long time but rather continues oscillating with a certain amplitude, and the average degrees of consolidation under cyclic loadings do not reach 100% in all cases ofCc/Ck<1 and Cc/Ck>1. It is due to repetitive loading and unloading stages of cyclic loadings, which results in squeezing out of pore water during the loading stage and then absorbing during the unloading stage. Thus, neglecting the effect of cyclic loadings may result in overestimated settlement and cannot consider the continuous volume change.It should also be pointed out that the present solution is only valid for the elastic constitutive model, which is not suitable for the elastic-plastic model of soft soil with creep. However, the analytical solution is first proposed for nonlinear radial consolidation of soft soil with vertical drains subjected to cyclic loadings. Also, the proposed solution is more general for the soil under time-dependent loadings and can be easily used by engineers.


## 6. Conclusion

In this paper, analytical solutions are derived for nonlinear radial consolidation of soft soil with vertical drains subjected to various cyclic loadings, such as haversine, trapezoidal, rectangular, and triangular cyclic loadings. The presented solution is verified through the degeneration into the existing solutions for nonlinear radial consolidation under constant and ramp loadings. The nonlinear radial consolidation behavior under various cyclic loadings is investigated. The results show that under cyclic loadings, both the dissipation rate of excess pore water pressure and the settlement rate slow down with the fluctuation corresponding to cyclic loadings whenCc/Ck increases. The dissipation rate increases when maximum loading qu increases in the case of Cc/Ck<1 but decreases with the increase in the maximum loading qu in the case of Cc/Ck>1, whereas the settlement rate increases with the increase in the maximum loading qu in all cases of Cc/Ck<1 and Cc/Ck>1. The higher the rate of loading increment or decrement is, the faster the dissipation rate and the settlement rate are. A longer rest period of loading induces the decrease in the rates of dissipation and settlement, less cycles, and a bigger amplitude of fluctuation. The longer the period of loading, the bigger the amplitude and less cycles of fluctuation in the rates.



---

*Source: 8810973-2020-06-30.xml*