# Cyclic Loading Test for the Small-Strain Shear Modulus of Saturated Soft Clay and Its Failure Mechanism

Authors: Zhende Zhu; Cong Zhang; Jun Wang; Peng Zhang; Duan Zhu
Journal: Geofluids
 (2021)

---

## Abstract

Small-strain shear modulus,Gmax, is a key evaluation index to study the dynamic characteristics of soil in geotechnical engineering. It is widely adopted to evaluate the stiffness of soft soil in soil dynamic engineering. In this paper, the cyclic triaxial tests and resonance column tests were carried out to explore the variation of Gmax of soft clay with respect to various confining stresses, cyclic shear stress ratios, pore pressures, and effective stress paths. Test results indicated that the effective stress decreased gradually with the increase of the cycle shear stress ratio. The failure points were mainly concentrated in a rectangular area, defined by the normalized effective stress from 0.56 to 0.64 and the normalized shear modulus from 0.72 to 0.78. Additionally, a short pause caused a small increase of 1-2% in Gmax as well as pore pressure. This study demonstrates that Gmax can be effectively used to characterize the failure of saturated soft clay in a more intuitive and convenient way, compared to the commonly used strain failure standards.


---

## Body

## 1. Introduction

Soft clay is found in coastal plains, estuary deltas, lake basins around, mountain valleys, and other areas, most of which belong to the middle or late Fourth Age. Figure1 shows the clay content of China. In general, soft clay has the characteristics of low carrying capacity, shallow water table, high water content, high compression, sensitivity, and large touch degeneration. These characteristics have a very negative impact on the long-term stability and normal operation of the upper structure [1]. That is why it is a subject that people have been interested in for decades [2–6]. The small-strain shear modulus (Gmax) is a key evaluation index in the seismic design and performance of geotechnical disciplines. As a fundamental property of the soil, accurate estimation and determination of Gmax is an important aspect of scientific and rational calculation and evaluation of liquid phase reaction, earthquake-induced deformation, and dynamic and static soil structure interactions [7].Figure 1
Clay content of China [12].For soft clay, soil destruction is a complex and cumulative process defined by strain as a standard of destruction (in general, it is 5 or 10%) [8]. Hyodo et al. [9] proposed the destruction criteria based on the cycle stress ratio (CSR) to explain that most relationship curves of soft clay had a distinct inflection point, and this inflection point corresponds to the strain as a destructive strain. Hyde et al. [10] presented a set of destruction criteria based on the extreme state of the active stress path. That is, the strain value corresponding to the effective stress path reaching the static failure limit state line can be used as the standard of destruction in uniaxial loading. Lei et al. [11] put forward the destruction standard of soft clay, by showing that the strain development of soft clay had a clear inflection point in the coordinate space of ε‐logN, and all the inflections were in a straight line. Thus, there is no uniform standard for the destruction of soft clay under dynamic load. Therefore, we judge the change of the soft clay structure based on the change of Gmax.The undrained triaxial tests were performed by Shirley and Hampton [13] to explore the stress-state dependency of Gmax under a constant void ratio (e). The results indicated that the stress-strain hysteretic loop, resilient modulus, and permanent strain were significantly dependent on the confining stress and the cyclic stress ratio (CSR). Thammathiwat and Chim-Oye [14] employed the cyclic triaxial test to study the interaction between the cyclic strength and the pore pressure of undisturbed Bangkok clay under various confining stresses, frequencies, and stress ratios. Meanwhile, several researchers have developed an uniaxial shear test instrument with multidirectional and looping capabilities to study the dynamic characteristics of soft clay [15–17]. These studies demonstrated that the dynamic shear modulus decreased, and the cyclic strain accumulation increased with an increase in CSR [18]. Besides, as the strain increased, the dependence of the modulus on the direction of the stress path decreased, and shear strain and volume strain also increased [19–21].Many indoor tests have been carried out and widely used to study the stress-strain properties of soil at small shear strain levels/amplitudes (γ<0.01%), for example, resonant column, bender elements, torsional shear tests, and cyclic triaxial tests [22–24]. Darabi et al. [25] conducted the cyclic triaxial tests of variable circumference pressure (VCP) to draw that the coupling of cycle deviatoric and circular confining stress can simulate the combination of cyclic shear stress and cyclic normal stress components. Compared to the resonant column device employed by Aggour and Zhang [26], stress and displacement were seen as control variables that control cyclic loadings. According to the recorded acceleration calculation of the time and orbit of the shear strain, it can be clearly seen that seismic effect is produced by multidirectional shear strain [27, 28]. It was carried out by Bowman and Soga [29] about experimental research on the shear modulus for temporal change effects based on various test materials such as indoor clean sand and Gulf of Mexico soft clay. The results were shown that under constant confining stress, the law of shear modulus growth varied for different material tests, while the shear modulus gradually increased with time.In spite of many studies dealing with the characterization of soft clay, the change law ofGmax is still not clearly known. More importantly, the deformation laws and criteria are not yet perfectly known under cyclic loads. There are numerous factors that affect the determination of Gmax. Moreover, due to sample disturbance and long consolidation time, the effect of time is often ignored, which could result in large errors. As a consequence, the parameters determined by laboratory tests cannot be directly used in engineering practice.To address the above problems, we first prepared the soft clay samples and measured the physical indexes (Section2) and described the test instrument and system parameters in Section 3. Then, we investigated the impact of cyclic loads on Gmax of saturated soft clay. The dynamic characteristics of soft clay were studied under different CSRs. Factors influencing the attenuation of Gmax were examined, and failure criteria of intact saturated soft clay on the basis of Gmax are determined in Section 4. The conclusions drawn from the experimental work and computational analysis done in this research are presented in the last section.


## 2. Sampling and Physical Properties

Soft clay in the test was derived from the Chashan Higher Education Park, Wenzhou University. Soil distribution maps within China and in Zhejiang Province are depicted as Figure2, and the stratigraphic map is illustrated in Figure 3.Figure 2
Soil distribution maps within China and in Zhejiang Province [12].Figure 3
The stratigraphic map of Wenzhou, Zhejiang Province [30].The soil depth is 5-7 m, which is typical for the marine soft clay on the southeast coast of China [31]. In order to obtain a sample of the intact soil with fewer disturbances, the thin-walled pipe cutting method was adopted to decrease the disturbance of intact soft clay. The soil collector was a thin-walled PVC pipe (diameter×height:160mm×300mm). Both ends of the PVC pipe were sealed with wax and maintained in the constant temperature and humidity incubator. To shorten the consolidation drainage time of the samples, round filter paper pieces were attached on the upper and lower ends of the samples, and filter paper was attached on the side and connected to the top filter. After saturation for 24 h, the value of B was set to 0.97 as the saturation completion criterion. The consolidation test was based on the excess pore water pressure dropping to 0 and displacement of less than 60 mm3/h as the standard. The mechanical parameters of soft clay are given in Table 1.Table 1
Mechanical parameters and values [31].
Soil parameterρ (g/m3)ω (%)ωp (%)ωL (%)IpPc (%)Fc (%)Value ranges1.7~1.7542~5026573154.394.6


## 3. Test Apparatus and Parameters

The tests are carried out using a combination of cyclic triaxial test system and piezoelectric ceramic bender element system developed by GDS, UK [31]. The dynamic triaxial test system can be used to do both static and dynamic tests. It includes five important components, as shown in Figure 4(a). The test items and accuracy of the test instrument are presented in Table 2. The piezoelectric ceramic bending element test system includes four key components, which are shown in Figure 4(b).Figure 4
GDS dynamic triaxial test system (a) and bending element test system (b).
(a)(b)Table 2
The test items and accuracy of the dynamic triaxial test system.
No.Test itemsRange valueAccuracy1Axial displacement100 mmMeasurement resolution reaches 0.08μm2Axial stress2 kPaAccuracy can reach 0.1%, that minimum resolution is 0.2 kPa3Confining stress range3 MPaAccuracy reaches 1 kPa4Back pressure range2 MPaAccuracy reaches 1 kPaThe test control and data recording of the dynamic triaxial test system are recorded in real time by the GDSLAB software, which can achieve the conventional triaxial compression test, the stress-controlled dynamic triaxial test, and the strain-controlled dynamic triaxial test. The method of bending element test is a simple and intuitive method with a clear principle. It is noted that the combination device used in this study can improve the accuracy of the wave velocity. Soil is a particle material and has a variety of contact forces-deformation behavior. When a point in the soil is stressed by the release of energy, a small strain will spread in the condition of elastic waves. The velocity of elastic wave depends on the elastic modulus of the soil. Therefore, if it can be measured for the small-strain elastic wave and velocity, the elastic parameters of the soil can be determined.Figure5 is the schematic shear wave velocity test of the soil sample by the bending. During the soil shear wave velocity test, the emission bending and receiving bending must be inserted into both ends of the soil. Then, the voltage pulse of the excitation signal at the specific frequency will be transmitted by the emission bending and be received by receiving bending. The vibrations generated by the shear wave are transformed into electrical signals. The signals are stored, and then, the propagation time (ts) of the shear wave is obtained by the signal comparison. The shear wave velocity (Vs) and shear modulus can be calculated by the distance of propagation (L).Figure 5
Working principle of shear wave velocity tested by bending element [32].


## 4. Impact of Cyclic Loads onGmax of Saturated Soft Clay

### 4.1. Test Method

Dynamic tests are performed using two methods: (1) Stress control: a certain number of constant amplitude cyclic load stress are imposed by the bending element system to testGmax. (2) Strain control: the control variable method is used to apply the dynamic load by controlling the dynamic strain method. Between two adjacent sets of cyclic shear vibration intervals, the shear wave velocity, residual pore pressure, and residual strain of soft clay are measured. The test is set based on short-term vibration interval (1-minute pause). Based on the actual occurrence of earthquakes, the loading frequency is chosen as 0.02 Hz, and 50 test points are recorded per cycle during the test (data is collected with 0.02-second sampling rate). Table 3 provides the details of the loading test method.Table 3
The test scheme design.
No.σm (kPa)f (Hz)CSRCyclic strain (%)1500.020.476—2500.020.551—3500.020.625—41000.020.489—51000.020.516—61000.020.571—71000.020.625—82000.020.362—92000.020.392—102000.020.447—112000.020.517—123000.020.432—133000.020.472—143000.020.569—151000.02—0.5161000.02—1171000.02—2182000.02—0.75192000.02—2202000.02—3Because of its structure, expectedly, intact soft clay has a nonlinear feature between its strength and confining stress. Moreover, dynamic characteristic of soft clay varied significantly under the same confining stress. To analyze the dynamic stress levels and characteristics under different confining stresses, it was expressed by the CSR in the test, i.e.,CSR=qcyc/qf, where qcyc is the imposed dynamics stress and qf is the shear strength under static conditions [31]. The samples are submitted to cyclic loading characterized by a sine wave. The total axial strain generated has stress amplitude of qcyc, and a total axial strain of εa,t is generated under cyclic loading. During the unloading process, part of the axial strain is called the rebound strain and is denoted by εa,r. However, the remaining part cannot be recovered, and it is accumulated with the increase of cycles, which is defined by the cumulative strain εa,p. The total strain is computed as εa,t=εa,r+εa,p. The curves of biaxial cycling obtained from the triaxial test results are depicted as in Figure 6.Figure 6
Curves of biaxial cycling obtained from the triaxial test results: strain development curve with cycle time (a) and a typical hysteresis (b).
(a)(b)


### 4.2. Dynamic Characteristics under Different CSRs

#### 4.2.1. Strain-Pore Pressure Curves

Figure7 depicts the change of axial strain change over time at 200 kPa and 0.392 (CSR), respectively. The tensile strain and compressive strain increased linearly, rapidly, under cyclic load as time increased, and the rebound strain also increased sharply. With the increasing number of cyclic loads, the samples continued to soften until they reached the state of failure. In the meantime, Gmax continuously decreased.Figure 7
Strain and pore pressure curve under 200 kPa confining stress: strain curve (a) and pore pressure (b).
(a)(b)For the confining stress of 100 kPa and 200 kPa, the curves of dynamic strain and different CSRs are illustrated in Figure8. The dynamic strain developed faster and faster with the increasing loading cycles, which indicates an exponential growth relationship. Under the constant confining stress, when the CSR was small, the strain increased very slowly. As the confining stress increased, the CSR gradually decreased, which resulted in an earlier exponential growth. Similarly, when the dynamic strength was large enough, the samples crossed the vibration-shearing stage and directly entered the vibration-destructive stage. In this case, the cyclic flow occurred and the dynamic strain increased approximately linearly.Figure 8
Curves of dynamic strain and different cyclic stress ratios: 100 kPa (a); 200 kPa (b).
(a)(b)Figure9 shows the pore pressure development curves when the confining stress was 100 kPa and 200 kPa, respectively. When the number of the loop loaded was small, the pore pressure increased rapidly and almost linearly. As the cycles increased, the increase in pore pressure slowed down and finally stabilized. The pore pressure increased with increasing CSR at 100 kPa (Figure 9(a)). The inflection points happened around 50 cycles, which was considered the boundary between linear and nonlinear processes. The pore pressure revealed the same growth law at 200 kPa (Figure 9(b)).Figure 9
Pore pressure development curves when confining stress is 100 kPa and 200 kPa.
(a)(b)


#### 4.2.2. Stress-Strain Relationships

When the confining stress is 200 kPa, stress-strain hysteresis loops of four different cyclic stress ratios (0.362, 0.392, 0.447, and 0.512) are shown in Figure10. Some of the hysteresis loops are labeled in Figure 10. Each group is loaded with a stress control mode. Therefore, local stress is fixed when the loop loads.Figure 10
Stress-strain hysteresis loops at different CSRs when confining stress is 200 kPa.One can observe three notable features in Figure10: (1) The axial strain developed slowly when the number of loop loaded was small (no more than ten). The axial strain increased rapidly with increasing cycles. For CSR=0.392, the strain values corresponding to the 10th, 100th, 200th, and 400th laps were 18%, 30%, 42%, and 60%, respectively. Similarly, the area of the hysteresis loop also increased. (2) The strain value increased rapidly with the increase of the CSR. The strain corresponding to the 10th lap increased slightly by 0.432%, 0.619%, 0.776%, and 1.346% in the case of CSR of 0.362, 0.392, 0.447, and 0.512, respectively. (3) In the case of a small cyclic stress ratio such as CSR=0.362, the resulting tensile strain and compressive strain were almost equal to the area of the hysteresis loop. For a larger cyclic stress ratio such as CSR=0.447 or 0.512, the tensile strain developed faster and the compressive strain developed relatively slowly, which was reflected in the shape of the hysteresis loop. The area enclosed by the tensile half-circle was greater than that of the compressed half-circle. It indicated that the tensile capacity was worse than its compressive capacity. Therefore, it was concluded that the tensile strain was more important than the compressive strain for various CSRs and dynamic stresses.


### 4.3. Factors Influencing the Attenuation ofGmax

The attenuation ofGmax is usually caused by joint action of both effective stress and dynamic strain. However, the dominant factors are not clearly known. To clearly understand the main reason for the abrupt change, the failure mechanism of soft clay is experimentally studied using a strain-controlled test. The test is conducted in 100 and 200 kPa to cycle the intact soft clay until damage to the structure of the soil. Figure 11 shows the normalized relationship curves between Gmax and effective stress.Figure 11
The change ofGmax under constant strain.As can be seen in Figure11, Gmax decreased when effective stress decreased under a constant strain control. However, the single strain value at 100 kPa was insignificant (ε=0.5%, 1%). When the effective stress decreased, Gmax gradually stabilized. Soft had not yet reached the full failure state. When ε was 2%, Gmax decreased with respect to the decreasing effective stress. Then, the attenuation amplitude increased until the final failure happened. When the failure strain of the samples did not happen, the modulus decreased as the effective stress decreased during the cyclic loading. Eventually, it stabilized and did not fail. Once the strain reached or exceeded the failure strain, the modulus decreased with respect to the decreasing effective stress. When confining stress was 200 kPa, Gmax and the effective stress revealed the same change law. The results showed that the increase of dynamic strain will weaken the interaction between the soil particles, which will lead to the decay of Gmax. The continuous development of dynamic strain is the main reason for this attenuation [31].


### 4.4. Failure Criteria of Intact Saturated Soft Clay on the basis ofGmax

At 100 kPa, 200 kPa, and 300 kPa, the law of the effect of cyclic loading history onGmax is plotted in Figure 12.Figure 12
Effect of cyclic loading history onGmax under different confining stresses: (a) 100 kPa; (b) 200 kPa; (c) 300 kPa.
(a)(b)(c)As can be seen,Gmax was 10% to 15% greater than the measured modulus after the cyclic load, as it had been shown by Cho [33]. The external force caused the disturbance of the soil structure, and the position and the contact force of the soil particles were changed. Moreover, the soil structure was destroyed when the effective stress decreased. It was difficult to restore the soil structure and strength in a short time, as it had been shown by Wichtmann and Triantafyllidis [34]. Furthermore, at the initial stage of the dynamic test, Gmax was almost constant and only showed a slight linear decline (in the range 22-25%). When the effective stress reached 75 kPa (Figure 12(a)), 120 kPa (Figure 12(b)), and 170 kPa (Figure 12(c)), Gmax decreased rapidly. This happens at a point where the soil is greatly deformed, i.e., the failure point. As can be seen, the failure point was not a fixed value under different CSRs. As the cyclic stress ratio increased, the failure point appeared earlier. The larger the CSR, the greater the effective confining stress corresponding to the failure point.The normalizedGmax/Gmax0 values and the normalized effective stress σm/σm0 values are plotted in Figure 13. Gmax0andσm0 represent the initial small-strain modulus and the initial effective stress, respectively. It shows a good correlation between the normalized Gmax and the normalized effective stress σm. The effective confining stress has a small impact on Gmax of saturated soft clay. Similarly, the rectangular area of failure points can be easily obtained. This area was concentrated in the normalized effective stress σm′/σm0′ from 0.56 to 0.64 and the normalized Gmax/Gmax0 from 0.72 to 0.78. This analysis demonstrates that Gmax can reflect the change and the failure of the soil structure. Figure 14 indicates the change laws between the strain εDA and the numbers of the loop loaded in this test. When confining stress was 100 kPa and 200 kPa, the turning points of various CSRs could be connected by a line. The test results were in agreement with Zhao et al. [31, 35, 36].Figure 13
Determined failure criterion area of intact soft clay based onGmax.Figure 14
Failure criteria on the basis of double-amplitude dynamic strain.By applying the cyclic loading, we found that the turning point of the double-amplitude strain of soft clay structure destruction was not constant, but it was closely related to the load amplitude. The double-amplitude strain value will increase with the increase of dynamic stress amplitude. The greater the dynamic stress amplitude, the earlier the corresponding turning point occurs. Moreover, these turning points are basically distributed on a straight line. The test results were in agreement with References [33, 37–39].From Figure14, we can draw the inference that Gmax can characterize the destruction of soft clay. The turning points of soft clay strain were basically on the same straight line, and the linear equation was basically consistent with the soft clay destruction criterion proposed by Cho [33]. It is represented by the following formula:
(1)εDA0=A×lgNf+B,where εDA0 is the corresponding strain value at the strain turning point and Nf is the vibration corresponding to the turning point.Therefore, we concluded that making use ofGmax as a criterion to assess the failure of soft clay was basically consistent with the commonly used strain failure standards. Compared with the commonly used failure standards, the principle of Gmax was more intuitive and convenient. Gmax can also reveal the characterization of the saturated soft clay. The increased dynamic strain and the decreased effective stress changed the soil structure and resulted in the attenuation of Gmax. A point with a fast change in Gmax will eventually appear as the failure of soil.


## 4.1. Test Method

Dynamic tests are performed using two methods: (1) Stress control: a certain number of constant amplitude cyclic load stress are imposed by the bending element system to testGmax. (2) Strain control: the control variable method is used to apply the dynamic load by controlling the dynamic strain method. Between two adjacent sets of cyclic shear vibration intervals, the shear wave velocity, residual pore pressure, and residual strain of soft clay are measured. The test is set based on short-term vibration interval (1-minute pause). Based on the actual occurrence of earthquakes, the loading frequency is chosen as 0.02 Hz, and 50 test points are recorded per cycle during the test (data is collected with 0.02-second sampling rate). Table 3 provides the details of the loading test method.Table 3
The test scheme design.
No.σm (kPa)f (Hz)CSRCyclic strain (%)1500.020.476—2500.020.551—3500.020.625—41000.020.489—51000.020.516—61000.020.571—71000.020.625—82000.020.362—92000.020.392—102000.020.447—112000.020.517—123000.020.432—133000.020.472—143000.020.569—151000.02—0.5161000.02—1171000.02—2182000.02—0.75192000.02—2202000.02—3Because of its structure, expectedly, intact soft clay has a nonlinear feature between its strength and confining stress. Moreover, dynamic characteristic of soft clay varied significantly under the same confining stress. To analyze the dynamic stress levels and characteristics under different confining stresses, it was expressed by the CSR in the test, i.e.,CSR=qcyc/qf, where qcyc is the imposed dynamics stress and qf is the shear strength under static conditions [31]. The samples are submitted to cyclic loading characterized by a sine wave. The total axial strain generated has stress amplitude of qcyc, and a total axial strain of εa,t is generated under cyclic loading. During the unloading process, part of the axial strain is called the rebound strain and is denoted by εa,r. However, the remaining part cannot be recovered, and it is accumulated with the increase of cycles, which is defined by the cumulative strain εa,p. The total strain is computed as εa,t=εa,r+εa,p. The curves of biaxial cycling obtained from the triaxial test results are depicted as in Figure 6.Figure 6
Curves of biaxial cycling obtained from the triaxial test results: strain development curve with cycle time (a) and a typical hysteresis (b).
(a)(b)


## 4.2. Dynamic Characteristics under Different CSRs

### 4.2.1. Strain-Pore Pressure Curves

Figure7 depicts the change of axial strain change over time at 200 kPa and 0.392 (CSR), respectively. The tensile strain and compressive strain increased linearly, rapidly, under cyclic load as time increased, and the rebound strain also increased sharply. With the increasing number of cyclic loads, the samples continued to soften until they reached the state of failure. In the meantime, Gmax continuously decreased.Figure 7
Strain and pore pressure curve under 200 kPa confining stress: strain curve (a) and pore pressure (b).
(a)(b)For the confining stress of 100 kPa and 200 kPa, the curves of dynamic strain and different CSRs are illustrated in Figure8. The dynamic strain developed faster and faster with the increasing loading cycles, which indicates an exponential growth relationship. Under the constant confining stress, when the CSR was small, the strain increased very slowly. As the confining stress increased, the CSR gradually decreased, which resulted in an earlier exponential growth. Similarly, when the dynamic strength was large enough, the samples crossed the vibration-shearing stage and directly entered the vibration-destructive stage. In this case, the cyclic flow occurred and the dynamic strain increased approximately linearly.Figure 8
Curves of dynamic strain and different cyclic stress ratios: 100 kPa (a); 200 kPa (b).
(a)(b)Figure9 shows the pore pressure development curves when the confining stress was 100 kPa and 200 kPa, respectively. When the number of the loop loaded was small, the pore pressure increased rapidly and almost linearly. As the cycles increased, the increase in pore pressure slowed down and finally stabilized. The pore pressure increased with increasing CSR at 100 kPa (Figure 9(a)). The inflection points happened around 50 cycles, which was considered the boundary between linear and nonlinear processes. The pore pressure revealed the same growth law at 200 kPa (Figure 9(b)).Figure 9
Pore pressure development curves when confining stress is 100 kPa and 200 kPa.
(a)(b)


### 4.2.2. Stress-Strain Relationships

When the confining stress is 200 kPa, stress-strain hysteresis loops of four different cyclic stress ratios (0.362, 0.392, 0.447, and 0.512) are shown in Figure10. Some of the hysteresis loops are labeled in Figure 10. Each group is loaded with a stress control mode. Therefore, local stress is fixed when the loop loads.Figure 10
Stress-strain hysteresis loops at different CSRs when confining stress is 200 kPa.One can observe three notable features in Figure10: (1) The axial strain developed slowly when the number of loop loaded was small (no more than ten). The axial strain increased rapidly with increasing cycles. For CSR=0.392, the strain values corresponding to the 10th, 100th, 200th, and 400th laps were 18%, 30%, 42%, and 60%, respectively. Similarly, the area of the hysteresis loop also increased. (2) The strain value increased rapidly with the increase of the CSR. The strain corresponding to the 10th lap increased slightly by 0.432%, 0.619%, 0.776%, and 1.346% in the case of CSR of 0.362, 0.392, 0.447, and 0.512, respectively. (3) In the case of a small cyclic stress ratio such as CSR=0.362, the resulting tensile strain and compressive strain were almost equal to the area of the hysteresis loop. For a larger cyclic stress ratio such as CSR=0.447 or 0.512, the tensile strain developed faster and the compressive strain developed relatively slowly, which was reflected in the shape of the hysteresis loop. The area enclosed by the tensile half-circle was greater than that of the compressed half-circle. It indicated that the tensile capacity was worse than its compressive capacity. Therefore, it was concluded that the tensile strain was more important than the compressive strain for various CSRs and dynamic stresses.


## 4.2.1. Strain-Pore Pressure Curves

Figure7 depicts the change of axial strain change over time at 200 kPa and 0.392 (CSR), respectively. The tensile strain and compressive strain increased linearly, rapidly, under cyclic load as time increased, and the rebound strain also increased sharply. With the increasing number of cyclic loads, the samples continued to soften until they reached the state of failure. In the meantime, Gmax continuously decreased.Figure 7
Strain and pore pressure curve under 200 kPa confining stress: strain curve (a) and pore pressure (b).
(a)(b)For the confining stress of 100 kPa and 200 kPa, the curves of dynamic strain and different CSRs are illustrated in Figure8. The dynamic strain developed faster and faster with the increasing loading cycles, which indicates an exponential growth relationship. Under the constant confining stress, when the CSR was small, the strain increased very slowly. As the confining stress increased, the CSR gradually decreased, which resulted in an earlier exponential growth. Similarly, when the dynamic strength was large enough, the samples crossed the vibration-shearing stage and directly entered the vibration-destructive stage. In this case, the cyclic flow occurred and the dynamic strain increased approximately linearly.Figure 8
Curves of dynamic strain and different cyclic stress ratios: 100 kPa (a); 200 kPa (b).
(a)(b)Figure9 shows the pore pressure development curves when the confining stress was 100 kPa and 200 kPa, respectively. When the number of the loop loaded was small, the pore pressure increased rapidly and almost linearly. As the cycles increased, the increase in pore pressure slowed down and finally stabilized. The pore pressure increased with increasing CSR at 100 kPa (Figure 9(a)). The inflection points happened around 50 cycles, which was considered the boundary between linear and nonlinear processes. The pore pressure revealed the same growth law at 200 kPa (Figure 9(b)).Figure 9
Pore pressure development curves when confining stress is 100 kPa and 200 kPa.
(a)(b)


## 4.2.2. Stress-Strain Relationships

When the confining stress is 200 kPa, stress-strain hysteresis loops of four different cyclic stress ratios (0.362, 0.392, 0.447, and 0.512) are shown in Figure10. Some of the hysteresis loops are labeled in Figure 10. Each group is loaded with a stress control mode. Therefore, local stress is fixed when the loop loads.Figure 10
Stress-strain hysteresis loops at different CSRs when confining stress is 200 kPa.One can observe three notable features in Figure10: (1) The axial strain developed slowly when the number of loop loaded was small (no more than ten). The axial strain increased rapidly with increasing cycles. For CSR=0.392, the strain values corresponding to the 10th, 100th, 200th, and 400th laps were 18%, 30%, 42%, and 60%, respectively. Similarly, the area of the hysteresis loop also increased. (2) The strain value increased rapidly with the increase of the CSR. The strain corresponding to the 10th lap increased slightly by 0.432%, 0.619%, 0.776%, and 1.346% in the case of CSR of 0.362, 0.392, 0.447, and 0.512, respectively. (3) In the case of a small cyclic stress ratio such as CSR=0.362, the resulting tensile strain and compressive strain were almost equal to the area of the hysteresis loop. For a larger cyclic stress ratio such as CSR=0.447 or 0.512, the tensile strain developed faster and the compressive strain developed relatively slowly, which was reflected in the shape of the hysteresis loop. The area enclosed by the tensile half-circle was greater than that of the compressed half-circle. It indicated that the tensile capacity was worse than its compressive capacity. Therefore, it was concluded that the tensile strain was more important than the compressive strain for various CSRs and dynamic stresses.


## 4.3. Factors Influencing the Attenuation ofGmax

The attenuation ofGmax is usually caused by joint action of both effective stress and dynamic strain. However, the dominant factors are not clearly known. To clearly understand the main reason for the abrupt change, the failure mechanism of soft clay is experimentally studied using a strain-controlled test. The test is conducted in 100 and 200 kPa to cycle the intact soft clay until damage to the structure of the soil. Figure 11 shows the normalized relationship curves between Gmax and effective stress.Figure 11
The change ofGmax under constant strain.As can be seen in Figure11, Gmax decreased when effective stress decreased under a constant strain control. However, the single strain value at 100 kPa was insignificant (ε=0.5%, 1%). When the effective stress decreased, Gmax gradually stabilized. Soft had not yet reached the full failure state. When ε was 2%, Gmax decreased with respect to the decreasing effective stress. Then, the attenuation amplitude increased until the final failure happened. When the failure strain of the samples did not happen, the modulus decreased as the effective stress decreased during the cyclic loading. Eventually, it stabilized and did not fail. Once the strain reached or exceeded the failure strain, the modulus decreased with respect to the decreasing effective stress. When confining stress was 200 kPa, Gmax and the effective stress revealed the same change law. The results showed that the increase of dynamic strain will weaken the interaction between the soil particles, which will lead to the decay of Gmax. The continuous development of dynamic strain is the main reason for this attenuation [31].


## 4.4. Failure Criteria of Intact Saturated Soft Clay on the basis ofGmax

At 100 kPa, 200 kPa, and 300 kPa, the law of the effect of cyclic loading history onGmax is plotted in Figure 12.Figure 12
Effect of cyclic loading history onGmax under different confining stresses: (a) 100 kPa; (b) 200 kPa; (c) 300 kPa.
(a)(b)(c)As can be seen,Gmax was 10% to 15% greater than the measured modulus after the cyclic load, as it had been shown by Cho [33]. The external force caused the disturbance of the soil structure, and the position and the contact force of the soil particles were changed. Moreover, the soil structure was destroyed when the effective stress decreased. It was difficult to restore the soil structure and strength in a short time, as it had been shown by Wichtmann and Triantafyllidis [34]. Furthermore, at the initial stage of the dynamic test, Gmax was almost constant and only showed a slight linear decline (in the range 22-25%). When the effective stress reached 75 kPa (Figure 12(a)), 120 kPa (Figure 12(b)), and 170 kPa (Figure 12(c)), Gmax decreased rapidly. This happens at a point where the soil is greatly deformed, i.e., the failure point. As can be seen, the failure point was not a fixed value under different CSRs. As the cyclic stress ratio increased, the failure point appeared earlier. The larger the CSR, the greater the effective confining stress corresponding to the failure point.The normalizedGmax/Gmax0 values and the normalized effective stress σm/σm0 values are plotted in Figure 13. Gmax0andσm0 represent the initial small-strain modulus and the initial effective stress, respectively. It shows a good correlation between the normalized Gmax and the normalized effective stress σm. The effective confining stress has a small impact on Gmax of saturated soft clay. Similarly, the rectangular area of failure points can be easily obtained. This area was concentrated in the normalized effective stress σm′/σm0′ from 0.56 to 0.64 and the normalized Gmax/Gmax0 from 0.72 to 0.78. This analysis demonstrates that Gmax can reflect the change and the failure of the soil structure. Figure 14 indicates the change laws between the strain εDA and the numbers of the loop loaded in this test. When confining stress was 100 kPa and 200 kPa, the turning points of various CSRs could be connected by a line. The test results were in agreement with Zhao et al. [31, 35, 36].Figure 13
Determined failure criterion area of intact soft clay based onGmax.Figure 14
Failure criteria on the basis of double-amplitude dynamic strain.By applying the cyclic loading, we found that the turning point of the double-amplitude strain of soft clay structure destruction was not constant, but it was closely related to the load amplitude. The double-amplitude strain value will increase with the increase of dynamic stress amplitude. The greater the dynamic stress amplitude, the earlier the corresponding turning point occurs. Moreover, these turning points are basically distributed on a straight line. The test results were in agreement with References [33, 37–39].From Figure14, we can draw the inference that Gmax can characterize the destruction of soft clay. The turning points of soft clay strain were basically on the same straight line, and the linear equation was basically consistent with the soft clay destruction criterion proposed by Cho [33]. It is represented by the following formula:
(1)εDA0=A×lgNf+B,where εDA0 is the corresponding strain value at the strain turning point and Nf is the vibration corresponding to the turning point.Therefore, we concluded that making use ofGmax as a criterion to assess the failure of soft clay was basically consistent with the commonly used strain failure standards. Compared with the commonly used failure standards, the principle of Gmax was more intuitive and convenient. Gmax can also reveal the characterization of the saturated soft clay. The increased dynamic strain and the decreased effective stress changed the soil structure and resulted in the attenuation of Gmax. A point with a fast change in Gmax will eventually appear as the failure of soil.


## 5. Conclusions

In this paper, the dynamic triaxial joint bending element test and the resonance column test were conducted for saturated soft clay in the southeast coast of China to investigate the interaction betweenGmax and soft clay mechanical parameters. The main conclusions obtained from the tests are refined in the following.
(1)
The larger the CSR, the smaller the effective stress reduction, and the earlier the occurrence of exponential growth. The increase in pore pressure is consistent with the trend of strain development. The failure points are not the same for different CSRs. The larger the CSRs, the earlier the occurrence of the failure point(2)
The rectangular area of failure points was drawn in this paper, which it defined by normalized effective stressσm′/σm0′ from 0.56 to 0.64 and the normalized Gmax/Gmax0 from 0.72 to 0.78. Under strain control, the sudden drop of Gmax reflects the broken structure of the saturated soft clay. Therefore, we can use Gmax as an intuitive and easy-to-use parameter to judge the failure of soft clay(3)
A short pause causes a small increase of 1-2% inGmax and pore pressure. The curves of pore pressure and strain were smooth and continuous



---

*Source: 2083682-2021-06-22.xml*