1. Introduction
Corals grow widely in tropical oceans between the Tropic of Cancer and the Tropic of Capricorn. When coral groups die, their skeletons and shells accumulate together, forming inanimate marine rock and sediment over a long period of geological processes. Due to the different settling speeds of sediment particles with different sizes, various sediment layers are formed under the seabed. Such layers can be divided into a coral sand layer and a coral mud layer, according to the general particle size. There are several-meter-thick coral mud layers on the seabed of the South China Sea. It is of great significance to carry out one-dimensional loading, unloading, and reloading characteristic analyses of coral mud for reclamation projects and island and reef construction in the South China Sea.
The formation principle of coral mud and coral sand is the same, and both are light yellow, have extremely high calcium carbonate contents, and their particles have marked edges and corners. At present, there have been many studies on the static and dynamic characteristics of coral sand in the South China Sea [
1,
2]. However, few studies have been carried out on coral mud. In fact, although coral mud has the particle-size characteristics of clay, it has an extremely high permeability and many properties similar to those of coral sand [
3]. Compared with coral sand, coral mud has cohesion and a high level of long-term deformation. Therefore, the study of coral mud is not only of practical engineering significance but also a link for a unified model of clay and sand.
The deformation of sediment under constant stress has an obvious nonlinear relationship with time, and many studies have been carried out on this nonlinear deformation. Terzaghi’s contribution to the principle of effective stress and consolidation theory [
4,
5] marked the beginning of modern soil mechanics. It is generally believed that the nonlinear deformation of a saturated sediment under the upper load consists of two parts; namely, a primary consolidation deformation due to the increase in effective stress and a secondary consolidation deformation due to the constant effective stress [
6]. The increasing effective stress stage is usually called the primary consolidation stage, and the stage after this, where the effective stress is constant but the deformation continues to develop, is called the secondary consolidation stage. Through a large number of tests, scholars have found that for clay, the secondary consolidation deformation under the action of constant effective stress accounts for a high proportion of the total deformation, which cannot be ignored. Some tests have even proven that this part of deformation lasts for decades and continues to develop [
7,
8]. Ladd et al. [
9] first asked whether there is secondary consolidation deformation in the primary consolidation stage when the excess pore-water pressure dissipates. There are two different assumptions about the secondary consolidation deformation of sediment in answer to this question; namely, hypothesis A and hypothesis B. Hypothesis A assumes that the secondary consolidation deformation only occurs after the primary consolidation stage, while hypothesis B considers that the secondary consolidation deformation occurs during the entire consolidation process. A large number of studies have proven that hypothesis B is more consistent with actual sediment deformation.
The isotach concept was first introduced by Sǔklje [
10], which is used to describe the effective rate of the compressibility of clay. Based on the isotach concept, Crawford [
11] carried out a long-term one-dimensional
consolidation test for normally consolidated clay and determined that the logarithmic curves of the void ratio and vertical stress for different loading times are a group of approximately parallel curves. Bjerrum [
12] proposed that the deformation process can be divided into instantaneous deformation and delayed deformation. Yin et al. [
13] proposed the EVP model, by introducing the equivalent-time concept. Yao et al. [
14] proposed the concept of an instant normal compression line on the basis of existing theories and experiments and introduced the conversion time into the yield equation of the unified hardening model, to establish the elasto–visco–plastic constitutive model under complex conditions. Feng et al. [
15] produced the EVP model, considering unloading expansion in the long-term time-dependent behavior of Hong Kong marine deposits (HKMD), under the loading stage and unloading stage in a multistage loading oedometer test. In addition, with the development of computer technology, some scholars have also used the grey theory method to predict the long-term one-dimensional deformation of sediment, such as Weng et al. [
16]. Due to the complex deformation characteristics of sediment, it is difficult to make an analogy with the deformation of other well-known materials. Some scholars have chosen to use basic element bodies in combination, to better simulate the deformation characteristics of sediment, such as in the Maxwell model [
17] and Merchant model [
18].
In engineering practice, it is found that some sediment layers have similar overconsolidation characteristics, which is called a quasi-overconsolidated sediment, although there is no unloading process in geological history. Bjerrum [
19] stated that, due to the aging effect, quasi-overconsolidated sediment has certain overconsolidation characteristics, and its shear strength is also higher than that of normal consolidated sediments. Most current long-term nonlinear deformation models only consider the time–strain relationship under a certain constant stress. Some overconsolidation models only consider the historical maximum vertical stress and cannot comprehensively consider all stresses borne by the sediment in its history and the impact of the duration of each level of stress on the deformation. In this article, based on the deformation of coral mud in one-dimensional oedometer tests, some elements of the Burgers model are modified. Furthermore, based on the modified Burgers model, a modified Burgers model considering stress history is obtained by considering all stress histories and loading durations.
3. Modified Burgers Model
3.1. Burgers Model
The strain of sediment after loading is generally divided into instantaneous strain and delayed strain. According to whether the strain can be recovered, the strain can be divided into elastic strain and plastic strain. Therefore, the total volumetric strain can be written as follows:
This is an example of Equation (1):
where
is the instantaneous elastic strain,
is the instantaneous plastic strain,
is the delayed elastic strain, and
is the delayed plastic strain.
To accurately reflect the deformation characteristics of coral mud after loading, the Burgers model was selected on a theoretical basis in this study, and a sketch is shown in
Figure 4. The Burgers model is composed of a Maxwell body and a Kelvin body in series. The model has four parameters, and each parameter has a clear physical meaning. The spring of the Maxwell body represents the instantaneous elastic strain after loading, while the other three elements jointly reflect the delayed strain.
The function of the Burgers model is
where
is the instantaneous strain of coral mud after compression,
;
is the nonrecoverable delayed plastic strain,
; and
is the recoverable delayed elastic strain,
. It is obvious that the Burgers model ignores the instantaneous plastic strain
. That is, the three factors of the Burgers model represent the instantaneous elastic strain, delayed plastic strain, and delayed elastic strain. Therefore, the Burgers model can generally reflect the nonlinear deformation characteristics of sediment after loading.
3.2. Shortcomings of the Burgers Model
In combination with the deformation characteristics of coral mud under loading in
Figure 3, it is not difficult to find problems with the Burgers model. At the moment of loading, the instantaneous strain of the sediment includes not only the instantaneous elastic deformation but also the instantaneous plastic strain. To more accurately describe the instantaneous strain of coral mud during unloading, it is necessary to modify the spring of the Maxwell body in the Burgers model. Taking the derivation of Equation (
2) yields
After loading, is always greater than 0; that is, the Burgers model maintains the feature of a monotonic increase, and the rate of strain increase finally tends to . This indicates that the Burgers model implies that the strain rate of sediment will tend to deform endlessly and never stop, which is obviously inconsistent with the actual situation.
When the Burgers model is under incomplete unloading conditions, it is assumed that the vertical stress at time
is
unloaded to
. The unloading equation is obtained as follows:
At the moment of incomplete unloading, the instantaneous elastic strain partially recovers; the rate of delayed plastic strain decreases, but still develops; and the delayed elastic strain partially recovers after unloading. Superimposing the strain changes of these three parts, the incomplete unloading curve is as shown in
Figure 5. The strain first decreases and finally tends to increase slowly and never stops, which is also inconsistent with the actual situation.
At the moment of incomplete unloading, i.e.,
in Equation (
4), the unloading equation is as follows:
Taking the derivation of Equation (
5) yields:
When unloading completely, the instantaneous elastic strain can be completely recovered in an instant; the delayed plastic strain is maintained at and does not decrease with unloading; and the delayed elastic strain will recover slowly. Therefore, when t tends to infinity, the final strain also tends to .
3.3. Modified Method
The instantaneous strain related to includes the instantaneous elastic strain and instantaneous plastic strain, i.e., . It is difficult to distinguish the two parts of the strain only from the loading test; therefore, the unloading test is needed for analysis. As the Burgers model only considers the instantaneous elastic strain, it is necessary to modify the spring of the Maxwell body; during the loading process, the modified Burgers model still retains the form of a spring; during unloading, only the instantaneous elastic strain can be recovered; thus, . Therefore, is no longer equal before and after continuous unloading and reloading.
Long-term indoor tests [
8] and practical engineering have proven that the deformation of sediment may “never stop”, but it can be clearly observed that this speed will tend to zero indefinitely. Therefore, it is unreasonable that the delayed plastic strain rate in the Burgers model tends to a fixed value; that is, the dashpot element of the Maxwell body is unreasonable.
Most component models are modified by adding more components. However, every additional element means that at least one parameter is added to the model equation. For a four-element model such as the Burgers model, the equation of the model is sufficiently complex. Adding elements does not improve the problem of infinite development and also increases the complexity and difficulty of the solution. Therefore, only the element properties are modified. In combination with the strain characteristics shown in
Figure 3, this paper considers modifying the dashpot into a nonlinear dashpot; that is, the viscosity coefficient of the dashpot is no longer constant but increases with time and tends to infinity, which is also consistent with the real strain of the sediment under loading. The modified Burgers model expression is as follows:
For the first level analysis of T1, as shown in
Figure 6, the predicted value of the modified Burgers model tends to be a straight line, which is more reasonable than the original Burgers model, which develops rapidly and constantly destroying.
3.4. Model Stability Verification
The fitting parameters may vary with the increase or decrease in the data points involved in the model calculation. To verify the stability of the modified Burgers model, the data of the first stage of T1 with different time lengths were substituted into the modified Burgers model for fitting, and the changes in the model parameters obtained are shown in
Figure 7.
With the increase in time involved in the model calculation, is relatively stable, and decrease slightly, and greatly increases. Namely, the instantaneous strain and delayed plastic strain are relatively stable. Although and jointly reflect the delayed elastic strain, this part of the strain finally approaches ; that is, it is only related to , while determines the speed tending to the maximum strain.
Although the time of substitution into the model is different, only the parameter in the modified Burgers model is significantly affected; that is, it only affects the speed at which the sediment tends to the maximum strain. The model can still accurately grasp the total strain of the sediment. The modified Burgers model has good stability.
4. Modified Burgers Model Considering the Stress History
At present, when analyzing a compression test of sediment under different loads, most methods analyze the time–strain relationship under each level of load. These analysis methods are fragmented and cannot take into account the influence of the historical loading, unloading, and reloading processes of the sediment [
20].
Figure 3 shows that the time–strain relationship under each load level is related to all previous load levels; therefore, this study extends the modified Burgers model.
As shown in
Figure 8, when the sediment is subjected to the first level load,
, the time–strain relationship is:
where
,
,
, and
are the first level model parameters, which are obtained by fitting the time–strain relationship under the first level load.
At time
, the load increases from the first level load
to the second level load
. At this time, the second level load
is divided into
and
. That is, the total strain is the sum of the strain produced by
in the whole continuous process and the strain produced by
from
. Then, the time–strain relationship is
where
,
,
, and
are obtained by fitting Equation (
8), and
,
,
, and
are the second-level model parameters, which are obtained by fitting the time–strain relationship under the second-level load.
Then, the time–strain relationship for the
i-th level load is
The time–strain relationship thus obtained can reflect the entire stress history that the sediment has been subjected to and the duration of each level of loading. In the fitting process of each level, only four parameters need to be solved, which ensures the model is not too complex.
The experimental data of T1 were substituted into the calculation, as shown in
Figure 9. The modified Burgers model considering the stress history can well represent the time–strain relationship of sediment under different loads.