Numerical Analysis of the Preloading Consolidation of the Soft Foundation with a Plastic Drainage Plate

The permeability coeﬃcient and void ratio are related to the stress state of the soil, which manifests as the temporal and spatial variability of consolidation. Based on the modiﬁed Cambridge model, this paper quotes the functional relationship between the soil permeability coeﬃcient and void ratio and soil stress obtained by Tylor’s research and then uses ﬁnite-element software ABAQUS and its secondary development platform to establish a plastic drainage plate surcharge preloading numerical model that considers the inﬂuence of the void ratio and permeability coeﬃcient with stress changes. The results show that considering the variability of the permeability coeﬃcient in the consolidation process of the soft foundation has a great inﬂuence on the change of the consolidation state but has a small eﬀect on the ﬁnal consolidation state. For soft foundation surcharge preloading reinforcement treatment, if the strength of the soft foundation is low but the surcharge is large, the geometric nonlinearity and material nonlinearity should be considered. Using the theory of large deformation consolidation to calculate the consolidation of soil is closer to the actual situation.


Introduction
Soft clay deposited by marine, lacustrine, and river facies is widely distributed along the coast and inland in China. is kind of soil is characterized by high compressibility, low strength, high water content, and poor water permeability. Under the effect of load, considerable settlement and settlement differences occur, which affects the construction quality and safety of the project and the normal use of the building.
With the development of economic construction in recent years, a large number of construction projects need to be built in soft soil and supersoft soil areas, and newer and higher requirements have been put forward for soft soil foundation treatment technology. At present, the main reinforcement treatment methods for soft soil foundations generally include the replacement method, the drainage consolidation method, the lime pile method, the highpressure jet grouting method, the deep mixing method, the sandstone pile method, reinforcement method, and so on [1][2][3][4][5][6][7]. Drainage plate preloading is a very common application in the drainage consolidation method. As an important member of foundation reinforcement technology, it is mainly used to treat silt or silt soil, flush fill, and saturated cohesive soil. e drainage plate preloading method is an effective soft foundation reinforcement method developed on the basis of the traditional sand well drainage preloading method. It is mainly used to treat foundations such as silt or silty soil, flush fill soil, and saturated clay soil. e basic principle of this method is to increase the drainage performance of the soil, improve the drainage conditions, and accelerate the drainage process of the soil by preloading, so as to promote the compression and consolidation of the soft foundation and improve its bearing capacity and stability.
As the basic theory of analysis and calculation in soft foundation reinforcement engineering, consolidation theory is the basis for soft foundation consolidation calculation. Terzaghi was the first to establish the basic differential equation of one-dimensional consolidation based on the principle of effective stress and the continuum equation of seepage flow. Because of its simple form and obtainable analytical solutions under certain initial and boundary conditions, it has been widely used in engineering. However, Terzaghi's consolidation theory is based on the assumption of constant compressibility and the permeability coefficient, which is obviously inconsistent with the actual situation that the permeability coefficient changes nonlinearly along the depth of the soil and changes with the change of the consolidation stress state. In addition, an important assumption in this theory is that the deformation of the soil is based on the theory of small deformation, which often results in large deviations in the calculation of the consolidation of soft clay foundations or foundations under heavy loads. In order to reflect the actual situation, Gibson et al. [8] assumed that the large-strain consolidation coefficient in the governing equation is linearly related to the void ratio, and the gravity coefficient is a constant. e finite-difference method was used to study the consolidation of thin soil layers. Poskitt [9] elaborated on the Gibson theory and considered the linear changes of soil compressibility and permeability during the consolidation process, using the perturbation method to obtain an analytical solution to the one-dimensional largestrain consolidation equation. Gibson et al. [10] assumed that both the large strain area coefficient and the gravitational coefficient are constant during the consolidation process, solved the equation and compared the calculation result with the solution. e research results show that the calculation result of the consolidation degree of the largestrain consolidation theory is safe. Davis and Raymond [11] obtained an analytical solution under constant load based on the assumption that the permeability coefficient changes linearly with depth, and the initial effective stress is constant. On the basis of Gibson's one-dimensional large-strain theory and Barron's free strain theory, Geng and Yu [12] present a governing equation which considers the radial and vertical flows, the weight of the soil, time-dependent surcharge, and vacuum preloading, along with variable hydraulic conductivity and compressibility during consolidation. Based on the analytical solution of one-dimensional linear elastic consolidation of layered foundations, Hui [13] used a semianalytical method to comprehensively analyze the properties of soil nonlinearity, rheology, soil weight, and variable load. Hong et al. [14] used a combination of explicit and implicit differences to solve the equation on the basis of the theory and studied the influence of factors such as initial conditions and boundary conditions on the calculation results. Geng et al. [15] presented an analytical solution that reflected the difference between the vacuum preloading systems under time-dependent surcharge preloading, and the results were in accordance with the laboratory and field results. Xie et al. [16] established a one-dimensional large-strain consolidation control equation that can consider the influence of load changes, soil weight, and other factors under Lagrange coordinates, using excess pore pressure as the control variable. Xie Xinyu [17] derived the one-dimensional large-strain surface bonding equation based on the mixture theory and believes that the consolidation theory based on the mixture theory is by far the most common consolidation theory. Applying the group transformation method, the explicit or implicit solution of the one-dimensional large-strain consolidation equation under certain specific conditions is obtained. In order to investigate the reinforcement effect of the SLR under vacuum-surcharge preloading, Wang et al. [18] conducted a series of laboratory tests for combination of vacuum pressure and surcharge load with different loading rates. Liu et al. [19] proposed an analytical solution that considered the effect of vacuum pressure loss, radial and vertical flow, smear effect, void ratio-dependent compressibility and permeability, and time-dependent surcharge preloading.
Relying on actual engineering, this article uses finiteelement software ABAQUS to perform a numerical simulation study on the soft foundation treatment of a plastic drainage plate with surcharge pressure considering the variability of the permeability coefficient and large deformation consolidation. e numerical calculation results are compared with the classical Terzaghi calculation results and the actual measured data to obtain the influence law of considering the variability of seepage and large deformation on the numerical calculation of soft soil consolidation and predict its development trend. In addition, the results are used to analyze the consolidation settlement of the soft foundation reinforced by the preloading of the drainage plate by the numerical simulation method's practicality and effectiveness in forecasting.

Project Conditions.
is project is the construction of a water system and landscape engineering on the west side of Fenghuang Avenue on Lingshan Island in Guangdong Province. ere are 5 rivers planned in this area, with a total length of 3868.1 m. Among them, the two existing rivers in the planned Zong Er Chong and Heng Er Chong are going to be rebuilt, with a length of 1584.2 m. e main task of the project is to reconstruct the internal water system on the west side of Lingshan Island to form a river network and ensure the connection of the water system. e main construction tasks are to complete the planned soil excavation, foundation treatment, ecological slope protection construction, and so on. In this simulation analysis, the planned longitudinal second swell was selected as the numerical simulation object for the preloaded soft foundation treatment, and the cross section BM650C with complete and reliable monitoring data was selected as the numerical modeling sample section on the swell.

Soil Properties.
e Zong Er Chong was rebuilt on the basis of the original Dong Wei Chong. According to the site geological survey report, the distribution of the foundation soil layer can be obtained as (1) silty clay layer, (2) silty sand layer, (3) sand layer, (4) silty clay layer, (5) sandy clay layer, and (6) bedrock layer. e simplified soil section is shown in Figure 1.
e physicomechanical parameters of soils are summarized in Table 1. 2 Advances in Civil Engineering

Foundation Reinforcement Plan.
e width of the Zong Er Chong is 43 m, and the soft foundation is reinforced. After cleaning the basement, the elevation of the bottom of the creek is 2.5 m. After the bottom of the creek is filled with sand to an elevation of 4.5 m, a 50 cm medium-coarse sand cushion is laid; the permeability coefficient is greater than

Advances in Civil Engineering
0.01 cm/s, and the drainage plate is preloaded. After the drainage plate is completed, a layer of 300 g/m 2 nonwoven geotextile is laid and then stacked and preloaded in layers. e stacking preloading backfilling shall be filled in layers, each with a thickness of 50 cm. e stacking load should be loaded evenly in layers to avoid excessive local stacking loads, and the loading time of the lower layer should be determined according to the settlement data. During the filling process, the settlement and stability of each layer are monitored. An ordinary B-type plastic drainage plate is used as the vertical drainage channel, arranged in a regular triangle, with a distance of 1.2 m. e plastic drainage plate has a cross-sectional width of 100 ± 2 mm, a thickness of 4.0 mm, longitudinal water flow, and a membrane permeability coefficient of ≥5 × 10 −4 m/s. e plastic drainage plate is required to pass through the silty soft soil layer and the weak underlying layer into the sand layer about 1.0 m.

Secondary Development of the Constitutive Model.
In this paper, the modified Cambridge constitutive model is used for the soft foundation. e modified Cambridge model in ABAQUS defaults so that the void ratio of the soil and the permeability coefficient remain unchanged. In order to make up for the deficiencies of the built-in model, user subroutines are written to realize the change of the void ratio and permeability coefficient with the consolidation stress. A large number of experiments and practices have proved that a semilogarithmic linear correlation between the permeability coefficient and void ratio empirical formula proposed by Taylor [20] in 1948 can well describe the relationship between the permeability coefficient and void ratio of soft soil. e formula is as follows: where e 0 is the initial void ratio, k 0 is the initial permeability coefficient, and C k is the slope of the e-logk curve, under normal circumstances, C k � 0.5e 0 . According to the formula, combined with the initial conditions, the relationship between k and e can be obtained, as shown in Figure 2.

Choice of Material Constitutive
(1) Constitutive model of plastic drainage plate: the plastic drainage plate is a typical composite material, and it is difficult to determine its constitutive relationship. For the sake of simplification, the plastic drainage plate is simplified into a continuous medium in engineering, and the linear elastic model is used to analyze it.
(2) Constitutive model of the dredger fill, sand cushion, and landfill: both the D-P model and the M-C model are suitable for granular materials. In order to keep the yield surface convex on the plane in the D-P model, a friction angle of less than 22°is required. e friction angle of materials such as the dredger fill and landfill is often greater than 22°, so the M-C model is used. (3) Foundation structure constitutive model: a lot of engineering practices and theories prove that the MCC model is suitable for normal consolidated clay and underconsolidated clay. erefore, the MCC constitutive model with nonlinear and hardening characteristics in this numerical simulation is selected.

Material Parameters.
According to the survey report and indoor geotechnical test and convert, the parameters of the drainage plate and the soils are shown in Tables 2 and 3.

Loading
Steps. According to the actual construction plan of the project, the loading process is divided into three main stages: sand filling (including sand cushion), layered filling, and preloading. ere is a certain time interval between each stage and each loading process, and it can be regarded as an isokinetic process. e loading process can be simplified as follows: 1 to 3 days are to fill the first layer of sand with a thickness of 1 m. After an interval of 2 days, the second layer of sand shall be blow-filled in 6 to 8 days with a thickness of 1 m. After an interval of 2 days, a 0.5 m thick sand cushion is laid out in 11 to 12 days. After an interval of 109 days, the landfill will start. e total thickness of the landfill is 2 m, and the landfill is divided 4 times. e first landfill thickness is about 72 cm, the total time is 3 days, and the preloading is 100 days. e thickness of the second filling is about 17 cm, which takes 2 days, and the precompression is 12 days. e thickness of the third filling is about 23 cm, which takes 2 days, and the preloading is 5 days. e thickness of the fourth filling is about 88 cm, which takes 3 days, and the preloading is 243 days. e loading curve is shown in Figure 3.  foundation model stacking area is 40 m, the horizontal influence area is 80 m, the calculating depth is taken to 47 m below the ground surface and the groundwater level is taken to the ground surface. From top to bottom, they are landfill, sand cushion, sand-filled layer, and different soil layers down. According to the grid division, the calculation model has a total of 6236 calculation nodes and 6077 calculation units. e soil element, drainage plate element, dredger fill element, and sand cushion element adopt the four-node plane pore flow element CPE4P, which is calculated by coupling of seepage and consolidation. e landfill element adopts the four-node plane strain element CPE4 without calculating its pore pressure. e simplified model is shown in Figure 4

Results and Analysis
According to the construction situation, the soil stress and strain, the change of excess pore pressure, and the displacement change during the preloading process are numerically simulated. e numerical calculation results and the monitoring values are compared and analyzed.

Analysis of Surface Settlement Results.
e project monitors the surface settlement of the left and right landfills at three locations. ey are BM650A-ZEC, BM650B-ZEC, and BM650C-ZEC (left pile, middle pile, and right pile), which are referred to as measuring points A, B, and C, respectively. Take Point B as an example. From the monitoring values in Figure 5, it can be seen that the settlement changes stepwise with the loading process, and the final settlement result is between 2.23 and 2.27 m.
From Figure 6, it can be clearly seen that the settlement value obtained based on the constant k is obviously greater than the monitoring value, reaching 2.62 m. e settlement change rate is relatively larger in each construction step. e settlement value calculated based on the nonlinear change of k with the void ratio is consistent with the monitoring values. Although there is a significant difference in the early stage of loading (within 70 days), the difference in the middle and late stages gradually decreases and tends to be stable. e main reason for this difference is that the construction and drainage conditions in the early stage of construction and the initial stage of the installation of the drainage plate are better than those of numerical simulation, which means the settlement speed obtained from the monitoring value is faster. However, when the plastic drainage plate works, the drainage conditions of the two are similar; the settlement value will gradually approach. e settlement curve calculated based on the large-strain theory is between the constant k and the change k curve, which means that the numerical simulation of the fluid-structure interaction based on largestrain theory is better than conventional consolidation theory.
From Figure 7, it can be seen that in the initial stage of construction, the relative error rate between the numerical calculation and the monitoring value is very large. As the construction progresses, the relative error rate under the condition of k changes gradually decreases and tends to be stable and remains at about 1%. Based on the constant k, the error rate also tends to decrease in the later period, but it remains at a high value and reaches about 18%. Figure 8 shows the time-history change curve of the layered settlement value measured by the 1#-6# magnetic ring of the layered settlement pipe on the left. From the monitoring results, the settlement law of the soft foundation at each depth corresponds to the construction loading process. In addition, the   Advances in Civil Engineering settlement deformation value measured from the 1# to 6# magnetic rings decreases successively, and the deformation amplitude also decreases. Due to stress diffusion, the additional stress generated by the upper load decreases as the depth increases. Figure 9 shows the layered settlement monitoring values and numerical calculation results of the 1#-3# magnetic ring with obvious settlement changes. e gap between the monitoring results and numerical values is large in the early stage of construction. As the construction progresses, the trend of change gradually approaches. is result is similar to the aforementioned surface settlement change law; the numerical result based on the nonlinear change of the permeability coefficient k with the void ratio is very close to the settlement result obtained by monitoring, which is obviously better than the result of constant permeability coefficient.

Analysis of Horizontal Displacement
Results. It can be seen from Figure 10 that there is a significant difference between the numerical calculation values and monitoring    e numerical calculation result changes significantly during the loading stage, while the monitoring result changes relatively smoothly. is difference is caused by complex field conditions, and numerical calculations are difficult to accurately simulate the actual situation. In addition, the displacement increases rapidly during loading, while the displacement decreases slowly during static preloading. Moreover, without considering the variability of the permeability coefficient, the horizontal displacement decreases more significantly. e soil is compressed during the consolidation process of the soft foundation. Under the effect of the drainage plate, this compression deformation is faster, so the pores of the soil become smaller and smaller, and the permeability worsens. Low permeability will hinder the drainage and consolidation of soil, which will slow down the deformation. Without considering the variability of the permeability coefficient, such obstacles will not occur. e soil will continue to compress and deform, resulting in an increase in the vertical compression and a decrease in the horizontal displacement.
is also explains that at the same time node, the horizontal displacement calculated without considering the permeability variability is smaller and more obvious than considering the permeability coefficient variability. Figure 11 shows the horizontal displacement curves of the soft foundation at different depths during each construction process, which are calculated based on different theoretical calculations. Figure 11   1# magnetic ring-measured data 2# magnetic ring-measured data 3# magnetic ring-measured data 4# magnetic ring-measured data 5# magnetic ring-measured data 6# magnetic ring-measured data   Horizontal displacement (cm) Advances in Civil Engineering 7 nonlinear variability of the permeability coefficient. Figure 11(b) shows the variation curve of the soft foundation in each construction stage when the permeability coefficient is constant. From Figures 11(a) and 11(b), the horizontal displacement from top to bottom generally shows a trend of large and small, which is more obvious before the construction stage of the plastic drainage plate. After entering the preloading stage, the horizontal displacement increases and then decreases in depth. is phenomenon is caused by the formation of stress bubbles in the soft foundation when the soft soil is subjected to a large upper load, and the soft soil is caused by the coordination of deformation. During the preloading stage of the soft foundation, its horizontal displacement will gradually decrease. During the loading stage, the surrounding soil is squeezed to cause horizontal displacement. e excess pore water pressure in the soft foundation is then quickly dissipated under the action of the plastic drainage plate, and the soft foundation undergoes a large vertical settlement. e surrounding soft soil tends to move closer to the center, which causes the horizontal displacement of the surrounding soft foundation to be reversed. erefore, compared to the previous loading process, the surrounding horizontal displacement is smaller.

Analysis of Excess Pore Water Pressure Results.
e distribution and dissipation of excess pore water pressure is an important indicator, reflecting the development of the consolidation process of soft foundations. A total of six pore pressure gauges are buried in the selected section, and they are buried symmetrically at different depths in the middle of the left and right landfills. ey are named A, B, C and D, E, F, and the pore pressure gauges are buried at 6, 11, and 18 m. Take the left as an example. e monitoring point is located within the depth of the plastic drainage plate. From Figure 12, the pore pressure changes significantly, showing the characteristics of rapid increase and dissipation. When considering the permeability variability of soft soil, the measured excess pore pressure at the same part of the soil is larger than the excess pore pressure without considering the permeability variability, and the dissipation rate of pore pressure is slower.
is gap is especially obvious in shallow soils. It decreases as the depth increases and eventually becomes consistent. As the upper soil compresses and deforms significantly under the load, the permeability coefficient of the soil will be significantly reduced, resulting in a slower dissipation rate of excess pore pressure than that obtained without considering the change of the permeability coefficient. As the depth increases, the compression of the soil layer gradually decreases, and the gap between the two gradually decreases.
From Figure 13, in the shallow and medium soils directly affected by the drainage plate, the excess pore pressure changes significantly, which is characterized by a rapid rise in the loading stage and a rapid decline in the precompression stage, especially in the filling construction stage and the large heap load. e excess pore pressure value corresponding to the construction stage of Landfills 1 and 4 is much larger, which exceeds the upper loading value. It is especially obvious in the sand filling construction stage, where there is no plastic drainage plate. is is called the Mandel-Cray effect [21,22].   Figure 11: Horizontal displacement changes in depth.

Conclusion
is article focuses on the nonlinear change of the permeability coefficient of soft soil along the depth. e consolidation theory of soft soil is developed by the ABAQUS secondary development platform. On this basis, combined with an engineering example of soft foundation reinforcement treatment in Lingshan Island rivers, finite-element calculations are carried out on the basis of field measuring data. Compared with the measuring data, the influence law of soft soil consolidation characteristics on the calculation of soft foundation reinforcement treatment is studied. e conclusions are as follows: (1) e choice of a constitutive model of soft soil consolidation has a great influence on the calculation results of consolidation, and it is more appropriate to use the modified Cambridge model for normally consolidated soft clay. (2) Considering the variability of the permeability coefficient during the consolidation process of soft foundations has a great influence on the change of the consolidation state, and the most direct factor is to affect the speed of consolidation. e consolidation speed, considering the permeability variability, is significantly lower than the consolidation speed without considering the permeability variability but has little effect on the final consolidation state. (3) For the preloading and precompression reinforcement of soft foundations, if the strength of the soft foundation is low and the surcharge is large, the geometric nonlinearity and material nonlinearity should be considered, and the large deformation consolidation theory is used to calculate the consolidation problem of the soil. (4) e calculated results by considering the permeability variability of soft foundations are more consistent with actual project monitoring values, which can better guide the selection of construction schemes and better predict and evaluate the development trend of soft foundation surcharge preloading reinforcement treatment.

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

Conflicts of Interest
e authors declare that they have no conflicts of interest. Excess pore water pressure (kPa) 6m-K is constant 18m-K is constant 11m-K is variable 6 m-measured data 11m-K is constant 6m-K is variable 18m-K is variable Load (2kPa) Figure 12: Comparison in soil excess pore water pressure at monitoring locations. Advances in Civil Engineering 9