Abstract
In order to grasp the characteristics and mechanism of saline soil deformation under the coupling action of cooling and dynamic load, three types of subgrade fillings from the Qarhan-Golmud Highway in Qinghai-Tibet Plateau were selected as experimental soil samples for the experimental study. Firstly, the freezing temperature experiment was carried out on SS (sandy silt) and HS (high-sulfate silty clay). The test results showed that the freezing temperature of SS is -0.32°C, while that of HC (high-chloride silty clay) and HS will not freeze at -20°C, due to the presence of salt. Secondly, the three soil samples were subjected to deformation characteristic test under coupling action of cooling and dynamic load, respectively, and the time history curve of temperature gradient change, the time history curve of the change rate of deformation, and deformation rate were summarized. Finally, the model of deformation rate vs. time and the model of change rate of deformation vs. deformation rate under the coupling action of cooling and dynamic load are proposed. The test results found that (1) the cooling rate of the temperature gradient curve of the three soil samples showed a rapid cooling rate in the early stage, and it tends to stabilize in the later stage. The distance of the 0°C line from the top gradually decreases, which is affected by the freezing temperature and the salt content. (2) Affected by the freezing temperature and salt type, SS exhibits frost heave, and HS and HC appear to settlement. The final deformation is 1.0%, -0.73%, and -1.10%, respectively. (3) The model of deformation rate vs. time and the model of change rate of deformation vs. deformation rate under the coupling action of cooling and dynamic load were proposed and verified, which are helpful for the evaluation of engineering stability on saline soil subgrade fillings.
1. Introduction
In recent years, with the rapid development of China’s economy, traffic construction occupies an increasingly important position. At present, motor vehicle ownership continues to rise, which also makes more and more scholars begin to pay attention to the study of the impact of deformation of roadbed pavement by traffic load. The number of projects carried out on traffic construction is rapidly increasing, and the scale of the projects is getting bigger and wider, so more and more geotechnical problems are manifested. These geotechnical engineering problems are mainly caused by two kinds of reasons: the deformation caused by the weak foundation resistance or its own foundation bearing capacity is low. Therefore, it is important to study the upper load of foundation and the deformation and strength of geotechnical engineering.
Akoto and Singh [1] determined the elastic and viscoelastic parameters of lime-improved red clay through cyclic triaxial tests and analyzed the effect of dynamic load on the modulus of resilience. Xu and Zhou [2] studied the shape of saturated soft clay under undrained cyclic loading, measured the pore water pressure of different cycles under cyclic loading, and established a characteristic suitable for saturated soft clay in Shanghai. Pritz [3] conducted dynamic triaxial liquefaction tests on undisturbed loess and remolded loess in the eastern United States under constant-amplitude cyclic loading and obtained the same liquefaction standards and failure standards of saturated loess as Prakash. Zhang [4] improved the expansive soil, using quicklime as the improver. Both the plain soil and the improved soil have been tested with static load and dynamic load, and a more comprehensive understanding of the engineering properties of lime-improved expansive soil has been obtained. Li and Zhou [5] conducted a dynamic study on fly ash type subgrade fillers and proposed that the vibration frequency, initial static stress level, and water content are all influencing factors of the dynamic strength of fly ash. Zhang et al. [6] conducted freeze-thaw cycle experiments on sulfate saline soils with 1.0% salt content and studied soil deformation, and the experimental results showed that salt crystallization was the main influencing factor of soil deformation. Wang et al. [7] conducted freezing experiments on chloride saline soils with salinity in the range of 0.1%-1% and concluded that it is not considered constant whether the soil system is in a frozen state, which should depend on soil texture, moisture, salinity, etc. Darrow et al. [8] conducted frost swelling tests on five natural fine-grained soils under a load of 200 kPa, and the experimental results showed that soils with colloidal organic components removed became less susceptible to frost. From the above studies, it can be found that most of the current studies on deformation characteristics of subgrade fillers and related load experiments are for nonsaline soils under frozen conditions or saline soils with low salt content, and there are few experimental studies on deformation of highly chlorinated saline soils under nonfrozen conditions.
Zhang et al. [9] conducted top freezing experiments for saline soils with 0.6% salt content. It was found that the temperature in the frozen zone changed faster than that in the unfrozen zone, and a model was developed to describe the coupled hydro-thermal-salt-mechanical (HTSM) process. Tian et al. [10] conducted deformation tests on a railroad roadbed fill in Inner Mongolia under the coupling effect of load and temperature and found that both static and dynamic loads had a certain inhibitory effect on soil frost swelling, and the frost swelling rate of this fine-grained soil gradually decreases with increasing external load values. Zhou [11] and Ye and Lai [12] studied soft clay and Hangzhou Bay silty soil, respectively, and found that the greater the vibration frequency, the smaller the dynamic strain generated by the sample and the higher the dynamic strength. Zhou [13] also used the stress-controlled cyclic triaxial test and found that the soil damage value is related to the cyclic stress ratio. At the same time, the test also found that the cyclic stress value is the main factor affecting the generation and growth of cyclic pore water pressure. Matesic and Vucetic [14] studied the effect of strain rate on the dynamic secant modulus under the condition of small strain for three types of clays and three types of sands. Khosla and Singh [15] conducted a series of cyclic loading tests on air-dry Ottawa sand and studied the influence of the number of loadings on the strain. As the number of cyclic loading increases, the sand will appear in three strain states. Sarkar and Mishra [16] established a saline soil erosion sensitivity model using binary logistic regression algorithm and artificial neural network. Zhang et al. [17] used BP neural network to influence the effect of expansion load on the suppression effect of salt expansion of 6 different types of output saline soil, including fine-grained soil gravel sulfate saline soil and fine-grained soil gravel sulfite saline soil. Zhao et al. [18] took saline soil as the research object; based on the British GDS dynamic triaxial test system, the freezing triaxial test was carried out under the temperature change state, and its dynamic stress-dynamic strain and dynamic shear modulus and other changing laws were studied under different temperature, frequency, and confining pressure. Zhang et al. [19, 20] investigated the effect of hydrothermal behavior on the frost swelling of saturated chalky clay soils. Four different pressures were applied to the top of each soil sample. Experimental results showed that the recharge water was the main component of the frost swelling. The increased applied pressure can limit the water migration and reduce the frost swelling during soil freezing and proposed that the frost swelling behavior of saturated chalk soil during unidirectional freezing is related to many factors, such as soil properties, water content, soil sample size, and environmental conditions. Previous studies have mostly focused on the dynamic load mechanical properties of nonfrozen soil under low-temperature environments and the factors affecting soil deformation. Few studies have focused on the deformation of nonfrozen saline soils with different salt types under coupled action of the cooling process and the dynamic loads.
Saline soils exist in most parts of China, including the Songnen Plain, Huang-Huai-Hai Plain, the coastal areas, the Hetao area of Inner Mongolia, and the inland areas of Northwest China [21, 22]. Especially in the western part of China, there are a large number of low-grade roads and rural roads, which are affected by climate, technology, and economic factors, and some of them have frost swelling and slurry diseases, which seriously affect the transportation business of China and even affect the pace of development of western China. A total number of geotechnical properties of subgrade soil should be investigated. According to the geotechnical properties and chemical analysis of soil samples from more than 200 sampling points along the Qarhan-Golmud Highway, the saline soil in this area is mostly low-plasticity fine clay (CL), with a salt content ranging from 0.2% to 48.5%. Chloride ion content is mainly concentrated at about 3% to 25%, and the sulfate ion content is basically stable around 0.3% to 2% in the whole process [23]. The salt content and the salt type directly affect the mechanical behavior of saline soil; the two typical high-chlorine saline soil and high-sulfur saline soil are selected as the research objects. Through the experimental study on the deformation of representative soil samples under the coupling effect of cooling and dynamic load, the deformation laws and adaptability of saline soil of different salt types are summarized, which are sufficient reference for the selection of subgrade soil in similar areas.
2. Engineering Background
The Qarhan Salt Lake saline soil area is located in the northeastern part of the Qinghai-Tibet Plateau in the Qaidam Basin. It is an inland basin at the northwestern tip of Qinghai Province, a component of the Central Asian desert, and belongs to the cold zone of the Qinghai-Tibet Plateau in the natural zoning of the highway. The annual precipitation in this area is only 20-143 mm, but the annual evaporation is as high as 2000-3000 mm, and the geographical location is inside the northwestern continent, surrounded by high mountains, making it difficult for warm and humid airflows to enter and being affected by high-temperature sunlight throughout the year. This is a typical continental desert climate; there are a large number of saline soils with more types and complex characteristics. Figure 1 is a map of salt crystals along the roadside of the Qarhan-Golmud Highway, where salt crystals are often seen on the ground surface.

The starting point of the road from Garsu to Golmud is located in Garsu, and the endpoint is located in Golmud City, with a total length of 34.7 km, and various types of saline soils are widely distributed. The major difference between saline soils and general soils is that saline soils contain a certain amount of soluble salts, which are completely dissolved in water when the water content is high, and the supersaturated salts will exist between soil particles in the form of crystals when the water content is low, playing a certain skeletal role [24, 25]. Therefore, saline soils are hard when the water content is low, and there is a certain cementation between the soil particles, possessing a certain bearing capacity. When the external conditions change (rainfall, irrigation, etc.), the soluble salt in the soil dissolves completely, and the soil structure is damaged to the extent that the strength decreases, and it is very likely that salt swelling, slurry, sinking and corrosion, and other disaster problems will occur, which are the main types of road diseases in saline soil areas in China [26, 27].
3. Experimental Design of Soil Column under Dynamic Load
3.1. Physical Parameters of Soil Samples
For the sandy chalk soil, the roadbed filler from the typical frost swelling and slurry section of National Highway 109 Rubber Mountain in Qinghai Province was used, while the samples of two saline soils were taken from along the Qarhan-Golmud Highway in Qinghai Province and divided into two groups according to the content of chloride salt or sulfate. The soil samples were subjected to basic physical properties and soluble salt chemical tests. The basic parameters of the soil samples are shown in Table 1. The contents of different salts in the two saline soil test soil samples are shown in Table 2. The soil samples were numbered. SS (sandy silt) is a sandy chalky soil, the salts in soil sample HC (high-chloride low-sulfate silty clay) in saline soils are mainly chloride salts with high total salt content, while soil sample HS (high-sulfate low-chloride silty clay) has a relatively high sulfate content. Figure 2 shows the particle size distribution of the three soil samples. According to [28], soil samples HC and HS are low-plasticity fine clay (CL), and soil sample SS is low-liquid limit sandy silt (SM).

3.2. Dynamic Load Test
The test equipment is shown in Figure 3. The specific test steps are as follows: (1)In this test, a plexiglass cylinder with a diameter of 150 mm and a wall thickness of 10 mm was used. A thin layer of petroleum jelly was evenly applied to the cylinder wall, and a plastic bag with an open bottom was pasted along the cylinder wall to ensure water inhalation during the test. Then, the plexiglass cylinder and its matching bottom plate were combined and the plastic bag was made to stick to the top of the plate, and then, the required amount of soil 670.57 g was calculated according to the compaction degree of 0.96, and the soil was divided into five layers and put into the plexiglass cylinder, and each layer was compacted according to the requirements. The control height of each layer is 20 mm. After compaction of one layer, the hair should be scraped and then the next layer should be compacted. A piece of filter paper is placed at the top and bottom end of the specimen to act as a water permeable and water impermeable(2)After the soil sample was prepared, the plate was mounted at the top and the temperature sensor was inserted into the temperature hole. From the top to the bottom cylinder, each inserted a temperature sensor 2 cm. After the specimen was made, the outer layer was wrapped with a sponge to provide insulation and left to stand for 24 hours so that the moisture in the soil sample was evenly distributed(3)Keep the temperature of the walk-in cryostat chamber at 5°C and adjust the sample temperature of the top and bottom plates of the cryostat to 5°C. The temperature of the soil sample is determined by an internal sensor until the temperature in the soil is 5°C from top to bottom(4)The temperature of the cryostat on the top plate is adjusted to -10°C, and the temperature of the bottom plate of the cryostat remains unchanged, allowing one-way cooling; at the same time, static load and dynamic load () are applied to the top of the soil sample. The pressure sine wave is adjusted to the dynamic load and its frequency is set to 1 Hz, and the load curve is shown in Figure 4. During the experiment; the temperature changes and soil deformation processes of the soil samples were recorded. The experiment ends when the internal temperature stability and deformation of the soil sample are stabilized

(a)

(b)

The height of the subgrade in this area is about 1.5 m~1.8 m. The maximum density of subgrade filler is 1960 kg/m3 and is taken; it is calculated that the applicable static load for this experiment is 30 kPa. The dynamic load condition in the test is based on [29]; a 20 t vehicle traveling at 60 km/h has a peak dynamic stress of 34 kPa on the surface of the subgrade. Combined with the existing test equipment and minimum speed limit of 80 km/h on Qarhan-Golmud Highway, the dynamic load strength is set to 60 kPa in this test.
4. Results and Discussion
4.1. Freezing Temperature Test Results
It is verified according to the fitting equation of the freezing temperature of saline soil given by Chen et al. [30] where is the optimal water content, Cl- is the chloride ion content, and SO42- is the sulfate ion content.
Figure 5 is the cooling and freezing time history curve of the test soil sample. It can be known from the test results that the freezing temperature of the SS is -0.32°C, and the freezing temperature of HS should be below -20°C.

According to Equation (1), the greater the chloride ion content, the lower the freezing temperature of soil. The content of the major negative ions in the soil sample was brought into Equation (1), the theoretical freezing temperature of the HS is -23°C, and the salt content of HC is higher than that of HS, so the freezing temperature of HC should be much less than -23°C.
4.2. Temperature Gradient Test Results
In the process of simulating the temperature gradient cooling of the field roadbed, all three types of soil samples were subjected to dynamic load tests during the temperature gradient cooling process at the optimum moisture content and 96% compaction. The height of the test soil sample is 10 cm, the top is the cold end (target temperature -10°C), and the bottom of the soil sample is the warm end (target temperature 5°C); under the effect of this temperature gradient, the soil sample characteristics, temperature gradient, and dynamic loading test conditions are shown in Table 3.
As the experiment progresses, the internal temperature of the soil will gradually change. The time course diagrams of temperature variation and temperature gradient diagrams of the three soil samples are shown in Figures 6–8.



Figure 6 shows the change process of the temperature field during the freezing process of sandy silt soil samples. During the first period of cooling, the temperature changes rapidly, but as time goes by, the temperature tends to stabilize. It was also found that the freezing temperature (-0.32°C) line is at a distance of about 8.6 cm from the top and the 0°C line is at a distance of about 8.9 cm from the top when the temperature tends to stabilize.
Figures 7 and 8 show the temperature time history change curves of HC and HS soil samples during freezing process, and their temperature time history changes are similar. It can be seen from the temperature time history chart that the temperature changes very quickly during the first period of cooling, but as time goes by, the temperature tends to stabilize. However, since the HC and HS soils will not freeze during the temperature gradient cooling process, the water in the soil migrates to the cold end at a slower rate than the SS soil, and the migration distance is also longer than that of the SS soil, so it shows the cooling process of SS is also slower than that of SS soil. After the temperature is stabilized, the distance between the 0°C line of HC soil and the top was 6.3 cm, and the distance between the 0°C line of HS soil and the top was 7.6 cm.
Therefore, the three types of soil samples under the same temperature gradient cooling effect all show a faster cooling rate in the early stage, and the temperature gradient change tends to be stable after about 40 hours of the test. However, due to the influence of the freezing temperature and water content of the soil, the temperature is lowered. The velocity and isotherm position are quite different.
4.3. Test Results of Soil Deformation under the Coupling Action of Cooling and Dynamic Load
The change rate of deformation, that is () to the original height of the soil (). According to Figure 9, it is the ratio of the amount of change in the height of the soil column to the initial height of the soil column during a certain test time. HC and HS do not freeze under the effect of this cooling process; therefore, the deformation of these two soil samples is caused by salt type and dynamic load.

Deformation rate, that is , change of soil sample height per unit time reflects the change of soil deformation gradient.
The soil deformation under the dynamic load of SS (sandy silt) can be seen from Figure 10. In the rapid freezing zone and transition zone, the gradient of the soil sample’s frost heave curve is relatively large. The amount of swelling is getting larger and larger, but the gradient of the curve gradually decreases, and finally, the amount of frost heaving of the soil tends to stabilize at about 1%. The frost heave rate of SS is shown in Figure 11. It can be seen from the figure that the frost heave rate of the soil sample increases first and then decreases with time. In the first few hours, the soil sample is in a state of rapid freezing, and the frost heaving rate gradually increases and reaches the maximum value. After entering the freezing transition zone and freezing-like stable zone, the water in the unfrozen zone in the soil sample migrates to the frozen zone, and the change rate of deformation gradually decreases until it finally stabilizes.


The curve of soil deformation rate versus time for HC under dynamic load can be seen in Figure 12. Under dynamic load, in the initial stage of cooling loading, the deformation is larger, and it is mainly due to the compression deformation of soil caused by static load and dynamic load; into the late stage of loading, its deformation rate tends to stabilize, and the final deformation rate is stabilized at -1.10%; the main reason for the decrease in soil deformation rate is the compression of soil voids prompted by dynamic load, in addition to the small salt shrinkage contraction deformation of the HC saline soil itself due to the decrease in temperature [31]. Figure 13 shows the deformation rate graph of HC, and it can be seen that the salt expansion rate varies greatly at the beginning of the test; however, it tends to level off for a period of time and then tends to decrease until finally leveling off. This is because at the beginning of the test, the soil is not highly compacted, and there are still more soil particle voids. After that, as the test proceeded, the soil temperature changed and salt shrinkage occurred in the saline soil, after which it tended to level off.


The relationship curve between soil deformation and time of the HS test soil sample under dynamic load can be seen in Figure 13. Under the action of dynamic load, the soil deformation rate is -0.73%, showing settlement. The initial stage of deformation is also due to the compression deformation of the soil caused by static load and dynamic loading. However, because HS soil samples contain more sulfates, during the cooling process, although the soil does not freeze, when the temperature drops to 15°C, the Na2SO4 in the soil begins to crystallize into Na2SO4·10H2O, the volume of the salt crystals. The growth rate reached 320% [32]. The salt swelling force in the soil can be balanced in a part of the dynamic load and finally lead to a settlement less than the HC soil sample. Figure 14 and 15 show the change rate of deformation and deformation rate for HS, which has a similar trend to HC. In the early stage of loading, the deformation rate increases rapidly with time and then enters a stable state. In the same way, at the initial stage of soil deformation, the originally existing soil particle interstices are compacted due to the 96% compaction. As the temperature decreases, salt crystallization begins to occur in the soil, which increases the volume of the soil. Finally, after the crystallization reaches equilibrium, the internal deformation of the soil tends to be smooth and no major changes occur.


5. Model Analysis
5.1. The Change Rate of Deformation and Time Model
(1)The change rate of deformation and time model of sandy silt
Figure 16 shows the deformation rate and time model of SS. The red curve is the fitted curve, and the white origin is the experimental data. The equation can be fitted according to the test results:

In the equation, is the deformation rate, is the time, and , , and are the test parameters, fitting to get , , , and . (2)The change rate of deformation and time model of saline soil
Figures 17 and 18 are, respectively, the deformation rate and time model diagrams of HC and HS. The red curve is the fitted curve, and the white origin is the experimental data. The equation can be fitted according to the test results:


The fitting parameters of HC and HS are shown in Table 4.
5.2. The Change Rate of Deformation and Deformation Rate Model
(1)The change rate of deformation and deformation rate model of sandy silt
Figure 19 shows the relationship between the freezing rate and frost heave rate of sandy silt. It can be seen that when the top cooling temperature is constant, the freezing rate of the soil sample decreases with the increase of the frost heave rate.

According to the test results, the following equation can be fitted: where is the strain of SS soil (%); is the freezing rate of SS soil (mm/h); and are the experimental parameters related to load, temperature, and soil properties. The experimental results fit the value of to 29.38409 and the value of to -0.16089.
It can also be seen from the above figure that the greater the freezing rate, the smaller the frost heave rate; the smaller the freezing rate, the greater the frost heave rate, and the relationship between the two is inversely proportional. This shows that when the freezing rate is small, the freezing front moves slowly, and the increase in freezing depth is small, but at this time, the frost heave rate of the soil sample is very large, resulting in the formation of ice lenses. (2)The change rate of deformation and deformation rate model of saline soil
Figure 20 shows the relationship between the salt swelling rate and instantaneous strain of HC. It can be seen that the instantaneous strain of the soil sample decreases with the increase of the salt expansion rate, and the equation can be fitted according to the test results: where is the strain of HC soil (%); is the frost heave rate of HC soil (mm/h). The experimental results fit the value of to 0.02832 and the value of to -12.49289.

Figure 21 shows the relationship between HS salt swelling rate and instantaneous strain. It can be seen that the instantaneous strain of the soil sample decreases with the increase of the salt expansion rate and finally tends to be flat. The equation can be fitted according to the test results: where is the strain of HS soil (%); is the frost heave rate (mm/h). The numerical results fit the value of to 0.0324 and the value of to -8.59993.

The fitting parameter values of the three experimental soils are shown in Table 5.
The distribution of salts in the soil sample is more inhomogeneous in high-salinity saline soils, and the sensitivity and dispersion of their phase change and water-salt migration under cooling conditions are more obvious than those of nonsaline soils, making the model fitting effect of high-salinity saline soils (HC and HS) lower than that of SS.
6. Conclusions
From the above research on the deformation characteristics of saline soil and sandy silt under dynamic load, the following conclusions can be drawn: (1)Under the same temperature gradient cooling effect, the three types of soil samples all show a faster cooling rate in the early stage, and the temperature gradient change tends to be stable after about 40 hours of the test. Due to factors such as the freezing temperature and water content of the soil, the cooling rate and the location of the isotherm are quite different, the freezing temperature (-0.32°C) line of SS is near 8.6 cm from the top, and the 0°C line is near 8.9 cm from the top; the 0°C line of HC is at a distance of 6.3 cm from the top; the 0°C line of HS is at a distance of 7.6 cm from the top. Affected by test result of freezing temperature and salt content, HC and HS will not freeze during the one-way cooling test process. The water in the soil migrates to the cold end at a slower rate than SS, and the migration distance is also longer than that of SS, so its cooling process is also slower than that of SS(2)Affected by freezing temperature and salt type, SS, HC, and HS show different deformation characteristics. SS showed frost heave of 1.0%, and HC and HS occurred settlement of -0.73% and -1.10%, respectively. In addition, salt expansion effect is caused by the presence of high sodium sulfate in HS and the action of cooling, which will balance part of settlement due to the action of dynamic load, so the deformation of HS is smaller than that of HC(3)According to the deformation characteristics of SS, HC, and HS, the model of deformation rate vs. time and the model of change rate of deformation vs. deformation rate under coupling action of dynamic load and one-way cooling are given and verified
Data Availability
The [DATA TYPE] data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The work herein is funded by the National Natural Science Foundation of China (No. 41501062) and the Longyuan Youth Innovation and Entrepreneurship Talent (Team) Project of Gansu Province. This support is much appreciated.