Studyof theLong-TermDeformationCharacteristicsofMunicipal Sludge Solidified Soil under the Coupling Action of Dry-Wet Cycles and Initial Static Deviatoric Stress

A self-developed curing agent is used to solidify the municipal sludge taken from Tianjin. .en, the long-term deformation characteristics of the sludge solidified soil are investigated by means of unconsolidated undrained creep tests with different drywet cycles for considering the influence of climate. .e experimental results show that the attenuation rate of the shear peak strength of municipal sludge solidified soil decreases gradually with the increase of the number of dry-wet cycles, and the strength remains unchanged when the number of dry-wet cycles is greater than 10. .e variation laws under different initial static deviatoric stresses are basically identical. When the applied stress is less than the yield stress of the sludge solidified soil, the duration curves of creep show only attenuated stage, i.e., with very small deformation, and the deformation reaches a constant in a short period of time. When the deviatoric stress reaches the long-term strength of the soil, the instantaneous deformation of the sludge solidified soil becomes large and damage occurs quickly. Under the same deviatoric stress, the creep deformation increases with the increase of the number of dry-wet cycles. When the load applied in each step is of the same magnitude, the higher the initial static deviatoric stress is, the larger the deformation of sludge solidified soil will be. It is found that the stress-strain relationship and the relationship between creep strain and time can be well described by an exponential function and a hyperbolic function, respectively. On this basis, a creep model is proposed for the long-term deformation considering the effect of dry-wet cycle times and initial static deviatoric stress. .e model is further validated by comparing the predictions with the test results under different deviatoric stresses; the good agreement between which shows the potential application of the model to relevant practical engineering.


Introduction
With the rapid development of economy and the acceleration of urbanization, municipal sludge treatment has become a hot issue that a great number of countries cannot avoid and needs to be solved urgently [1,2]. Sludge has high moisture content, poor mechanical properties and contains a large number of pathogens, microorganisms, heavy metals, and other harmful substances, which can cause serious secondary pollution and ecological problems if mishandled [3]. e traditional method of sludge treatment is inefficient and unsafe. e curing agent and skeleton materials are used to solidify municipal sludge, so that it can be recycled. After solidification, the physical and mechanical behaviors of the sludge are significantly improved, while the contents of organic matter and heavy metals are significantly reduced. Moreover, it is easier to be treated and transported than the original sludge, so that the municipal sludge can be transformed into the stable and safe subgrade filler for use to meet the engineering requirements and environmental requirements [4,5].
In recent years, many scholars have made great efforts and carried out a large number of experiments on the solidification of municipal sludge. By using different hybrid cement binders to add the sewage sludge with high organic content, Xin et al. [6] studied the unconfined compressive strength and leachability of heavy metals of sludge solidified soil. But when the cement is used to solidify sludge only, the strength of solidified body increases slowly, which is difficult to meet the requirements of resource utilization or landfill [7]. Qian et al. [8] investigated the potential for utilization of Municipal Solid Waste Incineration (MSWI) ash as solidification binder to treat heavy metals-bearing industrial waste sludge.
e results showed that the matrixes with heavy metals-bearing sludge and MSWI ash have a strong fixing capacity for heavy metals: Zn, Pb, Cu, Ni, and Mn. Cerbo et al. [9] used the fly ash for solidification of heavy metal sludge and added several kinds of cement additives such as sodium sulfate, sodium carbonate, and ethylenediaminetetraacetic acid (EDTA) to the solidified matrix in order to ascertain its physical and chemical characterizations. Many other materials are also used with cement to solidify municipal sludge. In view of the shrinkage and cracking of cement solidified sludge, Wang et al. [10] added different proportions of cement, quicklime, and sand into sludge and found that the strength and stability of solidified soil are greatly improved. When heat treatment, chemical treatment, and physical treatment were performed on already solidified soil, it was found that the structure of the solidified soil was more stable [11]. However, when certain proportion of jarosite and alunite was added to cement solidified sludge soil, it was found that the incorporation of them had a significant reduction in soil strength, while the leaching content of heavy metals was significantly reduced [12]. Ballesteros et al. [13] investigated the solidification/ stabilization of heavy metals in the electroplating sludge with the application of cement binder and admixtures of kaolinite clay, waste latex paint, and calcium chloride. Li et al. [14] found that adding an appropriate amount of waste glass powder was beneficial to the expansion of the sludge, which can make it possible to produce lightweight aggregates with high porosity and proper mechanical property.
On the other hand, the mechanical properties of soil under dry-wet cycles also arouse the interest of many investigators. Using c-ray and X-ray tomography techniques, Pires et al. [15] and Wang et al. [16] found that the microstructures and the pore spatial distribution of soils were significantly altered under the dry-wet cycles, resulting in the changes of its macroscopic mechanical properties. For example, the shear strength and specific surface area of soil decreased significantly when the soil undergoes several drywet cycles [17,18]. In practical engineering, one can always meet the situation that many tensile cracks occur at the surface of mortar block stones covering the terrace soil in several years after the construction and instability and failure of the low-angle low-expansion soil slope. is can lead to the collapses of the soil slope and lots of fractures in the piles of foundation [19].
Likewise, sludge solidified soil, under natural environmental conditions, will be inevitably subjected to the dry-wet cycle caused by environmental and climate change [20,21]. Whether its mechanical properties can meet the engineering requirements needs to be further explored. e aforementioned scholars have made great progress in the preparation of urban sludge solidified soil and understanding its fundamental mechanical properties. However, there are few reports on the mechanical properties of municipal sludge solidified soil under the dry-wet cycle, and there is even a lack of experimental data on the long-term deformation during operation process. In view of this, in this paper, a series of triaxial creep tests were carried out on the sludge solidified soil under unconsolidated and undrained conditions, so as to explore its long-term deformation law under the coupling action of dry-wet cycle and initial static deviator stress. en a model for reflecting the creep characteristics of solidified soil is proposed.
e model can provide valuable information to the prediction of solidified soil deformation under long-term loading.

Preparation of Solidified Earth Samples.
e test adopts marine soft clay taken from Tianjin Binhai as the skeleton material, and its basic physical properties are shown in Table 1. e municipal sludge used in the test was taken from a sludge landfill in Tianjin, which was solid and had a strong and rancid smell. After being retrieved, the sludge was packed into plastic boxes for sealing and preservation under shade, and the temperature was kept at about 20°C. e municipal sludge used in the test was degradation. e basic physical properties of municipal sludge are shown in Table 2.
Due to the high moisture content of municipal sludge and the large amount of refractory organic matter, the direct curing treatment requires a large number of curing materials and has poor curing effect. e sludge is digested with quicklime before solidification, so as to reduce the moisture content and decompose organic matter. e specific content of quicklime added to municipal sludge needs to be tested. By referring to relevant literature, the quicklime content was designed to be 5%, 10%, 15%, 20%, 25%, and 30% of sludge quality. A certain amount of sludge and different amounts of quicklime were mixed and stirred evenly for the purpose of digestion test. After 24 hours' standing, the moisture content and organic matter content of sludge with different quicklime content were measured, respectively. e test results are shown in Figure 1.
It can be seen from Figure 1 that the quicklime has a significant effect on the digestion of sludge. In design range, the moisture content and the organic matter content in sludge declined with the increase of the content of quicklime. At the beginning, the decrease extent was obvious; then the curves became almost stable. For the water content, when the content of quicklime exceeds 15%, the decrease rate becomes stable. For organic matter content, the curve becomes a flat straight line when the lime content exceeds 20%. Considering the requirement of technology and economy, it is optimal when the quicklime content is 20%. e curing agent used in the test is independently developed and this product has been applied for patents [22]. It is mainly composed of the main agent and the auxiliary agent of the curing agent. With the cement as the main agent, it plays a major role in the strength of sludge solidified soil. A small amount of auxiliaries A and B is added to catalyze the increase of strength growth rate and final strength. e toxicity of heavy metal leaching after curing meets the requirements of relevant specifications. e specific preparation scheme of municipal sludge solidified soil specimen is shown in Table 3.

e Manner of Dry-Wet Cycle.
e sludge solidified soil samples prepared were put into the standard curing chamber for curing with 7, 14, 28, and 60 days, respectively. en the basic physical properties of the samples after reaching the curing age and demoulding were tested. Table 4 shows the physical and mechanical properties of the solidified soil samples at the curing age of 28 days. It was found that the growth rate of unconfined compressive strength of the samples at the curing age of 28 days tended to be flat, which indicated that the internal reaction of soil mass was basically completed. erefore, the wetting-drying tests were carried out on the samples at the curing age of 28 days. e specific steps are as follows: first, put the samples with different times of dry-wet cycles into the oven to dry for 12 hours at low temperature (35 ± 3°C). en, take out the samples and leave them at room temperature (20°C) for 1 to 2 hours. Wrap the samples with the plastic wrap and place them on the porous stone. Place the sample, along with the porous stone, into a custom mould (as shown in Figure 2), and fill the mould with water until the porous stone is submerged a little beyond the bottom of the sample. After standing, the samples would absorb water due to the capillary action of pores. When the upper surface of the soil samples is completely wet, the samples would be taken out and placed in an incubator at 20 ± 3°C for 24 hours. Finally, remove the film, clean the samples, and place them at room temperature for 1 hour.
is is a dry-wet cycle process, which is shown in Figure 3. Repeat this operation until the specified number of times.

e Number of Dry-Wet Cycles.
e relation curves between triaxial shear peak strength and the number of drywet cycles n under different initial static deviator stress q are shown in Figure 4.
As can be seen from Figure 4, with the increase of the number of dry-wet cycles, the peak strength of the sludge solidified soil generally showed a trend of significantly decreasing and then tending to be stable, and the variation law was basically identical under different levels of initial loads. It can be also seen that the peak intensity attenuates rapidly with the increase of the number of drywet cycles when n ≤ 5, but the intensity attenuating rate gradually slows down. When n � 5, the attenuating amount accounts for more than 70% of the intensity attenuating amplitude. When n reaches 7 times, the curve flattens out. When n exceeds 10 times, the strength is basically stable and is no longer affected by the dry-wet cycle. is is mainly because the strength of sludge solidified soil mainly comes from the cementing substances produced by a series of hydration reactions between curing agent and soil particles and water. ese cementing substances increase the bonding force between various particles and make the soil have a certain resistance to the dry-wet effects. After the first few times of dry-wet cycles, there are many cracks in the soil due to the dehydrated shrinkage and swelling of water absorption, which leads to a significant decrease in soil strength. However, after a certain number of dry-wet cycles, the intergranular bonding force and the stress generated by dry-wet cycles gradually tend to balance and reach a steady state. So, the cracks no longer continue to expand, and the strength gradually becomes stable. is indicates that the shear resistance of soil decreases with the   increase of the number of dry-wet cycles. When a certain number of dry-wet cycles is reached, the influence of drywet cycles on the triaxial shear strength is no longer significant. According to the above analysis, it is reasonable to choose 0, 1, 3, 5, 7, and 10 times as the dry-wet cycles.

Testing Scheme
In order to research the long-term performance of municipal sludge solidified soil, TSS10 triaxial rheological test instrument was used to conduct creep test under unconsolidated and undrained conditions, so as to simulate the engineering situation with poor drainage conditions in the actual site. Considering that the subgrade soil has certain buried depth and will be subjected to additional static stress of the upper soil, the initial static deviator stress q of 0, 10, 15, and 20 kPa was selected for the study.
At present, there are two loading methods of creep test: graded loading and separate loading. e latter is to test a number of identical soil samples, while the former is to apply various levels of stress step by step on the same soil samples. In order to eliminate the influence of the differences between samples on the test results, the creep tests were carried out by means of graded loading. e creep loading was determined by the conventional triaxial shear test (UU), and the confining pressure σ 3 was chosen to be 20 kPa to simulate the condition of small fill thickness. Let s � (σ 1 − σ 3 ) represent the deviator stress. e creep test loading series is usually 5-8 groups, and then the increment of stress level of each class is calculated. As for the determination of creep stability criteria, there is no clear specification yet. In this paper, if the axial strain ε was less than 0.01 mm within 48 hours,       Advances in Civil Engineering the next stress level test would be conducted. e specific test scheme was shown in Table 5.

Creep Property
According to the available test results, the complete strain and time curves at each level of stress cannot be obtained directly. e curves are step-like and cannot be directly applied to analysis and research. Instead, the superposition numerical method is needed to transform them into curves under separate loading. In this paper, "Chen's method" [23] is used to establish the superposition relationship in the real deformation process.

Creep Testing Curve.
e creep duration curves of the sludge solidified soil under the action of different numbers of dry-wet cycle and initial static bias stress are shown in Figure 5.
It can be seen from Figure 5 that when the deviant stress is less than the yield stress, the creep curves of sludge solidified soil under different numbers of dry-wet cycles all present attenuation type of creep morphology. e overall trend of the duration curves under all steps is similar, which shows that the deformation amount is relatively small and tends to be stable within the deformation time. As time progresses, the strain growth rate slows down and eventually tends to a certain stable value. And the greater the deviator stress is, the longer it takes to reach the stable value. Under the action of the last step, the curve grows approximately linearly, which indicates that the deviator stress level has reached the long-term strength of soil, and the soil has a large instantaneous deformation after being acted on by external forces, resulting in a large deformation and rapid failure at the initial stage.
It can also be seen from Figure 5 that when the load applied in each step is of the same magnitude, the higher the initial static deviatoric stress is, the larger the deformation of sludge solidified soil will be. When the deviator stress levels are the same, the soil deformation at stability increases with the increase of the number of dry-wet cycles. is is because the soil structure disturbance increases with the increase of dry-wet cycles. Its internal connection becomes looser because of the increase of soil structure damage. Consequently, its stability and ability to resist external deformation are weakened. With the increase of stress level, creep deformation increases rapidly when the deviatoric stress exceeds the long-term strength of soil.

Establishment of Creep Model.
According to the creep curves in Figure 5, the strain of the solidified soil shows certain attenuation characteristics. ere are many kinds of creep models.
e models with few and easy-to-get parameters have some advantages. Exponential function, logarithmic function, power function, and hyperbolic function were used to fit the strain-time relation curves, respectively, with the least square method. e results show that the initial strain value is higher and the final settlement value is lower when the exponential function is used for fitting. By using logarithmic function and power function fitting, the strain grows slowly and continuously. It is worth noting that the strain increases rapidly in the later period and deviates from the test curve, while the hyperbolic relation is more practical [24].
According to the above analysis, the stress-strain relationship is the relation between the strain and the deviator stress (σ 1 − σ 3 ), the number of dry-wet cycles n, and the initial static deviator stress q. e stress-strain relation curve is tentatively expressed by the function f ((σ 1 − σ 3 ), n, q), while the strain-time relation curve is expressed by the hyperbolic function. erefore, the creep expression is assumed as shown in the following equation: Command B/A � C. And C is also a function of the number of dry-wet cycles n, initial static bias stress q, and deviator stress (σ 1 − σ 3 ). Transform equation (1) to where F is equal to (1/A) × f((σ 1 − σ 3 ), n, q). It is still a function of the number of dry-wet cycles n, initial static bias stress q, and deviator stress (σ 1 − σ 3 ).

Determination of Creep Model Parameters.
For the municipal sludge solidified soil studied, when the number of dry-wet cycles, initial static deviator stress and stress level are determined, the function F is a constant F 0 , which can be expressed as Transform equation (3) to get On the basis of equation (4), the creep curves can be transformed into the t/ε − t curves, and then linear fitting can be used to obtain the parameter values C 1 and C 2 taking into account the influence of different dry-wet cycles n and different initial static deviator stress q, respectively.
In order to explore the influence of dry-wet cycles on the long-term deformation of sludge solidified soil, the creep curves of initial static bias stress q � 0 kPa were selected as the research object. Figure 6 shows the t/ε − t relationship curves under different numbers of dry-wet cycles.

Advances in Civil Engineering
It can be seen from Figure 5 that the linear fitting of the t/ε − t relation curves get better effects, and the fitting degree R 2 can reach more than 90%. e intercept a 1 and slope b 1 of each fitting curve on the t/ε-axis are shown in Table 6.
As shown in Table 6, under the same initial static deviator stress, when the number of dry-wet cycles is constant, the parameter C 1 does not change much at each stress level.
at is, the magnitude of deviator stress has little influence on the value of parameter C 1 . And the average value of Advances in Civil Engineering parameter C 1 under different dry-wet cycles also changes little, so the final average value of parameter C 1 is set as 0.6478. When t � 96 h, the creep variable ε was taken to calculate the F 0 corresponding to different numbers of dry-wet cycles. Since F is the function of the number of dry-wet cycles n, the initial static deviator stress q, and the deviator stress (σ 1 − σ 3 ), when the initial static deviator stress and the number of dry-wet cycles are constant, F becomes a function of the deviator stress. e exponential function y � ab x is used to describe the relationship between F and the deviator stress level. e results are shown in Figure 7.
As shown in Figure 7, the exponential function can be used to describe the change relationship between F and (σ 1 − σ 3 ), with both R 2 reaching more than 95%. e parameter values of a 1 and b 1 are shown in Table 7. Table 7 shows that parameters a 1 and b 1 change to different degrees with the increase of the number of dry-wet cycles. erefore, the relationship between parameters a 1 and b 1 and the number of dry-wet cycles n is further studied, which is as shown in Figure 8. Figure 8 shows that it is suitable to use power function and linear function, respectively, to describe the relationship between parameters a 1 and b 1 and the number of dry-wet cycles n. And the fitting degree R 2 all can reach more than 96%. e specific expressions are shown in the following equations: In order to discuss the influence of initial static deviatoric stress on the long-term deformation of sludge solidified soil, the creep curves of sludge solidified soil under different initial static deviating stresses at n � 3 were randomly selected as the research object. e creep curves were also transformed into the t/ε − t curves and fitted linearly, as shown in Figure 9.
e parameter values are shown in Table 8.
Similarly, it can be seen from Table 8 that under the same dry-wet cycles, when the initial static deviator stress is  constant, the parameter C 2 does not change much at each stress level. at is, the magnitude of deviator stress has little influence on the value of parameter C 2 . And the average value of parameter C 2 under different initial static deviator stress also changes little, so the final average value of parameter C 2 is set as 0.5716. When t � 96 h, the creep variable ε was taken to calculate F 0 corresponding to different initial static deviator stress. Since F is the function of the number of dry-wet cycles n, the initial static deviator stress q, and the deviator stress (σ 1 − σ 3 ), when the initial static deviator stress and the number of dry-wet cycle are constant, F becomes a function of the deviator stress. e exponential function y � ab x is the same used to describe the relationship between F and the deviator stress level. e results are shown in Figure 10.
As shown in Figure 10, the exponential function can be used to describe the change relationship between F and (σ 1 − σ 3 ), with both R 2 reaching more than 95%. e parameter values of a 2 and b 2 are shown in Table 9. Table 9 shows that parameters a 2 and b 2 change to different degrees with the increase of the initial static deviator stress q. e relationship between parameters a 2 and b 2 and q is as shown in Figure 11. Figure 11 shows that it is suitable to use power function and linear function, respectively, to describe the relationship between parameters a 2 and b 2 and the initial static deviatoric stress q.
In the above analysis process, F is the function of the number of dry-wet cycles n, initial static bias stress q, and deviator stress (σ 1 − σ 3 ). erefore, when considering the influence of the dry-wet cycles, the initial static deviator stress is assumed to be a certain value, and then the relationship between F and the deviator stress level is analyzed. By using the exponential function y � ab x to fit the data values, the functional relationship between the values of a 1 and b 1 and the numbers of dry-wet cycles could be obtained. In the same way, when the influence of initial static deviator stress is taken into account, it assumes that the number of dry-wet cycles is a certain value. By using the exponential function y � ab x to fit the data values, the functional relationship between the values of a 2 and b 2 and the initial static deviator stress, the fitting results show that the number of dry-wet cycles and the initial static deviator stress all have influences on parameters a and b. Just the influence degree of impact varies. So, it can be considered that the parameters a and b are functions of the number of dry-wet cycles and initial static deviator stress, but the idea of single variable method is adopted in the analysis process. In the actual research process, when considering the influence of the dry-wet cycles, the initial static deflection stress is randomly selected. And when considering the influence of initial static deviator stress, the number of dry-wet cycles is also arbitrarily selected. erefore, it can be Advances in Civil Engineering approximately assumed that the functions of parameters a and b on the number of dry-wet cycles and the initial static deflection stress are the same function. On this basis, by the synthesis of equations (5)-(8), functions of a and b on the number of drywet cycles and initial static deviator stress are obtained as shown in the following equations: According to the above analysis results, the influence of deviator stress, dry-wet cycles, and initial static deviator stress on parameter C can be ignored. e parameter C can be approximately taken as the average value of C 1 and C 2 , which is equal to 0.6097.
To sum up, equations (9) and (10) are substituted into equation (2) to obtain the prediction formula of long-term deformation of sludge solidified soil under the action of different number of dry-wet cycles and initial static biased stresses. e results are shown as follows: where q represents the initial static bias stress, n is the number of dry-wet cycles, t represents time, and (σ 1 − σ 3 ) is the deviator stress. e values of other parameter in the equation are shown in Table 10.

Verification of the Creep Model.
In order to verify the accuracy and applicability of the creep prediction formula,          (11).  Advances in Civil Engineering 13

Advances in Civil Engineering
the experimental values under different initial static deviating stresses and different numbers of dry-wet cycles were compared with the creep values calculated by equation (11). Take the initial static bias stress q � 10 kPa and the number of dry-wet cycles n � 7 as an example. e predicted results are shown in Figure 12. Figure 12 shows that the calculated creep value is close to the experimental data. And the overall variation trend of the two curves is basically the same, which indicates that the prediction formula can meet the requirements of strain variation with time.

Conclusions
A series of creep tests were carried out on the sludge solidified soil under the condition of unconsolidated and undrained. en the characteristic of long-term deformation of the soil under coupling action of dry-wet cycle and initial static deviator stress was discussed. e key conclusions are as follows: (1) e attenuation rate of the shear peak strength of municipal sludge solidified soil decreases gradually with the increase of the number of dry-wet cycles, and the strength remains unchanged when the number of dry-wet cycles is greater than 10. e variation laws under different initial static deviatoric stresses are basically identical. After the first few times of dry-wet cycles, there are many cracks in the soil due to the dehydrated shrinkage and swelling of water absorption, which leads to a significant decrease in soil strength. However, after a certain number of dry-wet cycles, the intergranular bonding force and the stress generated by dry-wet cycles gradually tend to balance and reach a steady state. So, the cracks no longer continue to expand, and the strength gradually becomes stable.
(2) When the applied stress is less than the yield stress of the sludge solidified soil, the duration curves of creep show only attenuated stage. e development rate of strain slows down and quickly enters the stable stage. When the deviatoric stress reaches the long-term strength of soil, the instantaneous deformation of the sludge solidified soil becomes large and damage occurs quickly. When the load applied in each step is of the same magnitude, the higher the initial static deviatoric stress is, the larger the deformation will be. Under the identical condition of deviator stress, the creep variable increases with the increase of the number of dry-wet cycles. (3) Based on the hyperbolic strain-time relationship and stress-strain index function relation, the influence of the number of dry-wet cycles and the initial static stress on the deformation was analyzed by using the single-variable method. On this basis, a creep model is proposed for the long-term deformation considering the effect of dry-wet cycle times and initial static deviatoric stress. e model can provide theoretical support for the prediction of deformation of solidified soil under long-term loading.

Data Availability
Some or all data used to support the findings of this study are available from the corresponding author upon request (E-mail: gfxu@whrsm.ac.cn).

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.