The Inﬂuence of Freeze-Thaw Cycles on the Mechanical Properties and Parameters of the Duncan-Chang Constitutive Model of Remolded Saline Soil in Nong’an County, Jilin Province, Northeastern China

: In seasonally frozen areas, physical and mechanical soil properties change dynamically under the e ﬀ ect of the freeze-thaw cycles (FTCs), which is a problem that cannot be ignored in geotechnical engineering. In order to study the e ﬀ ect of the FTC on the strength and the mechanism of deformation and failure of saline soil, this paper took Nong’an saline soil as the research object. In total, 105 groups of remolded samples with di ﬀ erent salt contents (S) after FTCs were examined in unconsolidated-undrained tests. On the basis of the experiment results, the inﬂuence of FTCs on the mechanical properties of Nong’an saline soil was analyzed. The failure principal stress di ﬀ erence ( σ 1 − σ 3 ) f and cohesion (c) were both decreased with FTC. This occurred especially rapidly after the ﬁrst cycle and became stable between 30 and 60 cycles. The internal friction angle ϕ increased at ﬁrst and then decreased. According to experimental data, a modiﬁed Duncan-Chang model was established. Compared with the experiment results, this model was reasonable to simulate the stress-strain relationship of Nong’an saline soil. Furthermore, the empirical formulas of Duncan-Chang model parameters were obtained by regression analysis. This provides a theoretical basis for saline–soil foundation and subgrade engineering in seasonal frozen areas.


Introduction
In seasonal frozen-soil regions, the freeze-thaw cycle (FTC) is an important factor affecting soil's physical and mechanical properties [1]. As the temperature decreases, water in soil turns into ice, volume expansion increases pores, and the freezing process leads to aggregate separation and the breaking of soil particles. Changes in soil structure during the freezing process do not completely recover as the temperatures rises [2]. Under the effect of repeated FTCs, the soil structure is repeatedly adjusted, so its physical and mechanical properties are changed, such as porosity [3,4], severity [5], hydrologic properties [6], permeability [7,8], and strength and deformation [9]. Thus, soil mechanical properties are key considerations in construction projects in cold areas.
Because soil-structure changes caused by FTCs have a significant influence on soil properties, there are many previous studies on shear-strength changes of soils after FTCs [5,[9][10][11]. In one, the expansive soil remained stable after the first cycle, but soil volume increased, and c and ϕ decreased as

Study Area and Soil Properties
The soil of the remolded samples was collected at the county of Nong'an in the western province of Jilin. According to the curve of depth and S, 40 cm was usually the turning point and the highest point [18,19]. Therefore, samples in this paper were made from a 40 cm layer soil. Through laboratory tests, the maximum dry density of this layer was found to be 1.63g·cm −3 , and the optimal moisture content was 21.26%. The moisture content of all samples was therefore set as 22%, and the compaction degree was 90%. Zhang et al. (2016) [20] studied the physicochemical characteristics and water-salt transport regularity of saline soil in Nong'an of the province of Jilin. In 2013, soluble soil salt from the 0-40 cm layer in two sampling points of Nong'an varied from 0.1% to 0.8%. In 2014, the soluble salt of the soil above 1 m depth varied from 0.2% to 3.0%. In this experiment, the total soluble salt contents of samples were set as 0%, 0.5%, 1%, 2%, and 3%, respectively. Soluble salt is the mixture of NaHCO 3 and Na 2 SO 4 , and NaHCO 3 :Na 2 SO 4 = 1:1.
Nong'an meteorological data from 1981 to 2010 shows that temperatures are below 0 • C every year from November to March of the following year, and the soil-freezing period lasts for five months. Temperatures of the remaining months are above 0 • C, which is the soil-melting period. The monthly mean temperatures are in the range of −16.3 to 23.1 • C; the extreme minimal temperature is −37.2 • C, and the extreme maximal temperature is as high as 37.8 • C. So Nong'an in western Jilin is a typical seasonal frozen-soil area. Freezing greatly influences soil properties [21], and it is not conducive to the stability of geotechnical engineering builds, such as foundations, highways, railways, bridges, and tunnels, or slope engineering [19,[22][23][24][25]. Soil is dispersed saline soil [26], and an FTC causes salt expansion [22] leading to pavement deformation, cracking, and bulging, and road boiling.

Sample Preparation
In order to prepare samples of the triaxial shear test with different S, desalination was the first step. We put the distilled water and soil collected from the field in a plastic bucket and stirred well to ensure the soil was soaked and dissolved in the distilled water. After soaking for 24 h, HNO 3, AgNO 3 , and a small amount of hydrochloric acid and BaCl 2 were dropped into the 2 tubes with filter liquor, Appl. Sci. 2019, 9, 4941 3 of 18 respectively. If white precipitate was generated, water was changed repeatedly in the bucket until the white precipitate was no longer generated. Then, we placed the desalted soil on trays and let it naturally dry.
According to the standard for the soil-test method (GB/T50123-1999) [27], the crushed dry soil under a 2 mm diameter sieve was used to make the remolded samples. First, NaHCO 3 and Na 2 SO 4 solutions were mixed with dry soil. Then, the wet soil was sealed in fresh-keeping bags and placed in a wet cylinder for 24 h so moisture was evenly distributed in the soil. Finally, the soil sample was compacted into a 100 mm high and 50 mm in diameter cylinder with the stratified-compaction method.

Freeze-Thaw Cycle Test
The freeze-thaw cycle test was carried out in the laboratory, and the comprehensive simulation platform of rock and soil freeze-thaw cycle test under an ultracold environment was adopted. Freezing temperature was set to −20 • C. The wrapped triaxial samples were frozen in a closed freezing chamber for 12 h. Then, temperature was set to 20 • C for 12 h for samples to be completely thawed. FTCs in this study were set as 0, 1, 5, 10, 30, 60, and 120. When a sample completed the set FTCs, the soil column was trimmed into a cylinder with a diameter of 39.1 mm and a height of 80 mm.

Triaxial Shear Test
At present, most models are based on consolidated drainage triaxial tests. In fact, shallow surface soil is often in an unsaturated state, and clay soil has low permeability. Therefore, in the compaction process, the soil consolidated and crept at the moment of loading, and generally stayed in an unconsolidated and undrained state. So, the stress-strain relationship and strength of saline soil after freeze-thaw cycles being studied through an unconsolidated and undrained triaxial shear test was suitable in this paper. The triaxial shear test adopted a strain-control 3-axis instrument. Whenever axial deformation changed by 0.2 mm, axial deformation and dynamometer readings were recorded. If total axial strain was greater than 3%, every time axial deformation changed by 0.5 mm, readings were recorded. When total axial strain was 15%-20% or the sample was destroyed, the test was completed. As shown in the Table 1, the moisture content of the samples was 22%, compactness was 90%, and the salinity was set to 0%, 0.5%, 1%, 2%, or 3%, FTCs were set to 0, 1, 5, 10, 30, 60, and 120, and confining pressure (σ 3 ) was set to 100, 200, or 300 kPa. A total of 105 groups of unconsolidated and undrained shear tests were carried out.

Experiment Results and Analysis
The relational curves between principal-stress deviation (σ 1 − σ 3 ) and axial strain ε 1 of samples with 0.5% salt content are shown in Figure 1. Because the laws of stress-strain curve affected by FTCs were similar under different S, the stress-strain curve with 0.5% salt content is illustrated here.
As shown in Figure 1, the principal stress deviation (σ 1 − σ 3 ) had nonlinear growth with the increase of axial strain ε 1 . When σ 3 was the same, the (σ 1 − σ 3 ) -ε 1 curve-variation tendency of the samples with different FTC was similar. All (σ 1 − σ 3 ) -ε 1 curves were of the strain-hardening type, and the shape of curves was hyperbolic. The stress-strain curve could be divided into three stages: the quasi-elastic stage, the elastic-plastic deformation stage, and the stable stage.
In the repeated process of freezing and thawing, salt expansion, pore water migration, particle breaking and rearranging, and fracture development, gradually changed the force between soil particles, affecting the macromechanical behavior of saline soil [2, 28,29]. The initial yield point of the samples decreased with more FTCs (Figure 1). The samples that underwent more FTCs were easier to damage, and their strength was lower.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 18 samples with different FTC was similar. All (σ1 − σ3) -ε1 curves were of the strain-hardening type, and the shape of curves was hyperbolic. The stress-strain curve could be divided into three stages: the quasi-elastic stage, the elastic-plastic deformation stage, and the stable stage.
In the repeated process of freezing and thawing, salt expansion, pore water migration, particle breaking and rearranging, and fracture development, gradually changed the force between soil particles, affecting the macromechanical behavior of saline soil [2, 28,29]. The initial yield point of the samples decreased with more FTCs (Figure 1). The samples that underwent more FTCs were easier to damage, and their strength was lower. For strain-hardening soil, there is no peak on the stress-strain curves, so principal-stress deviation when the strain value is 15% is usually defined as the shear strength of Nong'an saline soil. It was recorded as (σ1 − σ3)f and called failure principal-stress difference, which is summarized in Table 2. (σ1 − σ3)f of soil varied with FTC change, as shown in Figure 2. The following exponential function was used for curve fitting:  For strain-hardening soil, there is no peak on the stress-strain curves, so principal-stress deviation when the strain value is 15% is usually defined as the shear strength of Nong'an saline soil. It was recorded as (σ 1 − σ 3 ) f and called failure principal-stress difference, which is summarized in Table 2. (σ 1 − σ 3 ) f of soil varied with FTC change, as shown in Figure 2. The following exponential function was used for curve fitting: where x is the FTC number that the samples underwent, and A 1 , B 1 , C 1, and D 1 are the fitting parameters. Fitting results are shown in Figure 2, and parameter values are shown in Table 3.  The (σ 1 − σ 3 ) f of the samples with different S decreased rapidly within 10 cycles, especially about 31.22%-36.42% after the first cycle. At this stage, the structure of Nong'an saline soil was damaged. Then, it slowly declined between 10 and 60 cycles, and it declined faster after samples underwent 120 cycles (24.85%-57.13%). During the first freezing, water in the soil became ice and the volume expanded, which made small particles agglomerate and the aggregates separated from each other. Moreover, pore volume in the soil increased. When the ice melted, aggregate changes and the pore volume in the soil could not recover completely during the thawing process. After the first cycle, the structure of Nong'an saline soil was seriously damaged and its strength decreased rapidly. After 5 and 10 cycles, the micro and small pores in the soil gradually became medium and large pores, large particles or aggregates were separated, and small particles were reunited; finally, the soil entered into a new dynamic and stable equilibrium state [10,30]. Therefore, the decreasing rate of (σ 1 − σ 3 ) f was obviously reduced when samples underwent 30 and 60 cycles. However, as the number of FTCs increased, pores in the soil gradually spread into cracks, so soil strength had an accelerated decline after 120 cycles.
The exponential function can reflect the relationship between (σ 1 − σ 3 ) f and FTC. Correlation coefficient R 2 was above 0.93, so this function is suitable to explain the relationship between (σ 1 − σ 3 ) f and FTC.
(σ 1 − σ 3 ) f increased with confining pressure, which indicated that the latter could inhibit pore-volume changes caused by axial pressure. So, confining pressure can slightly enhance the strength of the samples and have a certain inhibitory effect on deformation and failure [11,31].
Appl. Sci. 2019, 9, x FOR PEER REVIEW 6 of 18 and 10 cycles, the micro and small pores in the soil gradually became medium and large pores, large particles or aggregates were separated, and small particles were reunited; finally, the soil entered into a new dynamic and stable equilibrium state [10,30]. Therefore, the decreasing rate of (σ1 − σ3)f was obviously reduced when samples underwent 30 and 60 cycles. However, as the number of FTCs increased, pores in the soil gradually spread into cracks, so soil strength had an accelerated decline after 120 cycles. The exponential function can reflect the relationship between (σ1 − σ3)f and FTC. Correlation coefficient R 2 was above 0.93, so this function is suitable to explain the relationship between (σ1 − σ3)f and FTC.
(σ1 − σ3)f increased with confining pressure, which indicated that the latter could inhibit porevolume changes caused by axial pressure. So, confining pressure can slightly enhance the strength of the samples and have a certain inhibitory effect on deformation and failure [11,31]. c and the φ are two important indices for evaluating soil's shear strength that are mainly affected by particle-size composition, arrangement, and cementation [32][33][34]. Recurrent FTC and changes in S change these factors, thus affecting soil c and φ. According to the experiment data, c and φ of the sample could be obtained, as shown in Table 4. As shown in Figure 3, c and φ were related to FTCs.  c and the ϕ are two important indices for evaluating soil's shear strength that are mainly affected by particle-size composition, arrangement, and cementation [32][33][34]. Recurrent FTC and changes in S change these factors, thus affecting soil c and ϕ. According to the experiment data, c and ϕ of the sample could be obtained, as shown in Table 4. As shown in Figure 3, c and ϕ were related to FTCs. The exponential function (Equation (1)) could explain the relationship between c and FTC, while the Fourier function could explain the relationship between ϕ and FTC. The two functions' parameters are shown in Table 5. where x is the number of FTCs the samples underwent, and A2, B2, C2, and D2 are the fitting parameters. According to Figure 3, the exponential function can reflect well the variation regulation of c with FTC. The overall trend of c decreased with FTC number, and c decreased considerably when FTCs varied from 0 to 10, especially when samples underwent the first cycles. After that, decline was slow and c was basically stable at around 30 kPa between 10 and 60 cycles; φ increased when below 30 cycles, and then decreased rapidly after 60 cycles. The Fourier function reflects the law of φ with FTC very well.

Determination of Duncan-Chang Model Parameters
As shown in Figure 1, relational curves between principal-stress deviation (σ1 − σ3) and axial strain ε1 of samples are hyperbolic. Therefore, the Duncan-Chang model was suitably applied to analyze the relationship between sample stress and strain [13]. The relationship between model parameters and FTC is discussed below.
Kondner et al. (1963) [13] analyzed the stress-strain curve of soil morphology on the basis of a large number of triaxial shear-test data, and suggested that the hyperbolic function was suitable for simulating the stress-strain relationship of soils: where a and b are parameters that could be determined by conventional triaxial shear test.
Since dσ2 = dσ3 = 0 in the triaxial test, tangent elastic modulus Et could be obtained by taking the derivative of Equation (3) and eliminating axial strain ε1: The Fourier function is as follows: where x is the number of FTCs the samples underwent, and A 2 , B 2 , C 2, and D 2 are the fitting parameters. According to Figure 3, the exponential function can reflect well the variation regulation of c with FTC. The overall trend of c decreased with FTC number, and c decreased considerably when FTCs varied from 0 to 10, especially when samples underwent the first cycles. After that, decline was slow and c was basically stable at around 30 kPa between 10 and 60 cycles; ϕ increased when below 30 cycles, and then decreased rapidly after 60 cycles. The Fourier function reflects the law of ϕ with FTC very well.

Determination of Duncan-Chang Model Parameters
As shown in Figure 1, relational curves between principal-stress deviation (σ 1 − σ 3 ) and axial strain ε 1 of samples are hyperbolic. Therefore, the Duncan-Chang model was suitably applied to analyze the relationship between sample stress and strain [13]. The relationship between model parameters and FTC is discussed below.
Kondner et al. (1963) [13] analyzed the stress-strain curve of soil morphology on the basis of a large number of triaxial shear-test data, and suggested that the hyperbolic function was suitable for simulating the stress-strain relationship of soils: where a and b are parameters that could be determined by conventional triaxial shear test.
Since dσ 2 = dσ 3 = 0 in the triaxial test, tangent elastic modulus E t could be obtained by taking the derivative of Equation (3) and eliminating axial strain ε 1 : where E i is the initial elastic modulus, and (σ 1 − σ 3 ) ult is the ultimate principal-stress deviation. Ultimate principal-stress deviation is difficult to be directly obtained through the test, but failure principal-stress difference is easy to determine. Therefore, in order to obtain the ultimate principal-stress deviation, the failure ratio was defined as follows: According to Mohr-Coulomb strength theory: We transformed the coordinate system of Equation (3) to get the following equation: According to Equation (7), triaxial test data are transformed into another coordinate, and curves of ε 1 /(σ 1 − σ 3 ) and ε 1 were drawn in Figure 4. It can be seen from Figure 4 that there was a first-order linear correlation between stress and strain in the ε 1 /(σ 1 − σ 3 ) − ε 1 coordinate. The stress-strain curve of the sample conformed to the hyperbolic function. The y-intercept of the fitting line was parameter a, and the slope was parameter b, as shown in Table 6. The reciprocal of parameters a and b are E i and (σ 1 − σ 3 ) ult , respectively. Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 18    In Figure 5, under different confining pressures, the variation tendencies of initial tangent modulus E i and FTC were similar, decreasing rapidly below 30 cycles, and then decreasing slowly after 30 cycles. The regression relation of E i to FTC could be fitted by the exponential function, and fitting parameters are shown in Table 7. In Figure 5, under different confining pressures, the variation tendencies of initial tangent modulus Ei and FTC were similar, decreasing rapidly below 30 cycles, and then decreasing slowly after 30 cycles. The regression relation of Ei to FTC could be fitted by the exponential function, and fitting parameters are shown in Table 7.  The failure principal-stress deviator (σ1 − σ3)f and limiting principal-stress deviator (σ1 − σ3)ult were determined. According to Equation (5), fracture ratio Rf could be obtained. As shown in Figure  6, Rf increased after the first cycle, and then decreased between five and 60 cycles. After 120 cycles, Rf obviously increased. The following rational function was used for curve fitting: = (A 3 * 3 + B 3 * 2 + C 3 * + D)/( + E 3 ), (8) where x is the FTC number that the samples underwent, and A3, B3, C3, D3, and E3 are fitting parameters. The fitting result is shown in Figure 6, and parameter values are shown in Table 8.  The failure principal-stress deviator (σ 1 − σ 3 ) f and limiting principal-stress deviator (σ 1 − σ 3 ) ult were determined. According to Equation (5), fracture ratio R f could be obtained. As shown in Figure 6, R f increased after the first cycle, and then decreased between five and 60 cycles. After 120 cycles, R f obviously increased. The following rational function was used for curve fitting: where x is the FTC number that the samples underwent, and A 3 , B 3 , C 3 , D 3 , and E 3 are fitting parameters. The fitting result is shown in Figure 6, and parameter values are shown in Table 8.  According to Janbu's suggestion, initial tangent modulus Ei is related to σ3, and its empirical relationship is as follows: where K and N are material constant and Pa = 101.4 kPa is the standard atmospheric pressure. We drew the diagram of log10(σ3⁄Pa) and log10 (Ei/Pa), as shown in Figure 7. According to the relationship between log10(σ3/Pa) and log10(Ei/Pa), the y-intercept and slope of the fitting line were log10(K) and N, respectively, so material constants K and N were obtained, as shown in Table 9.
With the increase of FTC, the change rule of material constant K rapidly decreased below 30 cycles, and then slowly decreased after 60 cycles. The regulation of N increased below 10 cycles, and then decreased between 10 and 60 cycles; it then increased after 120 cycles. The regression relation of K-FTC and N-FTC could be fitted by the exponential function; parameters of this function are shown in Table 10. And the fitting result is shown in Figure 8.   According to Janbu's suggestion, initial tangent modulus E i is related to σ 3 , and its empirical relationship is as follows:

Freez-thaw Cycles
where K and N are material constant and P a = 101.4 kPa is the standard atmospheric pressure. We drew the diagram of log 10 (σ 3 /P a ) and log 10 (E i /P a ), as shown in Figure 7. According to the relationship between log 10 (σ 3 /P a ) and log 10 (E i /P a ), the y-intercept and slope of the fitting line were log 10 (K) and N, respectively, so material constants K and N were obtained, as shown in Table 9.   With the increase of FTC, the change rule of material constant K rapidly decreased below 30 cycles, and then slowly decreased after 60 cycles. The regulation of N increased below 10 cycles, and then decreased between 10 and 60 cycles; it then increased after 120 cycles. The regression relation of K-FTC and N-FTC could be fitted by the exponential function; parameters of this function are shown in Table 10. And the fitting result is shown in Figure 8.   Stress level SL is defined as the ratio of actual principal-stress difference (σ1 − σ3) to failure principal-stress difference (σ1 − σ3)f: By substituting Equations (5), (6), (9), and (10) into Equation (4), we can get the following formula: Stress level S L is defined as the ratio of actual principal-stress difference (σ 1 − σ 3 ) to failure principal-stress difference (σ 1 − σ 3 ) f : By substituting Equations (5), (6), (9) and (10) into Equation (4), we can get the following formula: Simplified: In conclusion, the experience formula of the five material constants contained in the tangent deformation modulus of the Duncan-Chang model were established, so the variation regular of the tangent deformation modulus with the FTC could be determined by unconsolidated undrained tests.
In the triaxial unconsolidated, undrained shear tests, no volumetric strain occurred in the samples. So, material constants G, F, and D that calculate tangential Poisson's ratio in the E-µmodel (initial tangent modulus -tangent Poisson's ratio model) could not be obtained, nor could the K b and m of the E-B model (initial tangent modulus -bulk modulus model) be obtained. According to the generalized Hooke's law, the value of µ t was 0.5 when the sample strain was 0. In order to ensure the stability of numerical calculation, Poisson's ratio is usually 0.49 [35].

Model Validation
By using Equations (5), (6) and (9), we could calculate hyperbolic model parameters a and b. According to Equation (6), the calculated results could be drawn in a stress-strain diagram, and we compared them with the experimental stress-strain curves. As shown in Figure 9, stress-strain curves that were determined by the modified Duncan-Chang model were, considering the FTC, basically close to the experiment stress-strain curves; the simulation effect was better in the stage when the strain was small and there was nearly elastic deformation. When the strain increased to 4%, there was a difference between some simulated and experimental curves. On the whole, it was reasonable to use this model to simulate the stress-strain curve of Nong'an saline soil affected by FTCs under the condition of an unconsolidated, undrained triaxial test.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 15 of 18 In conclusion, the experience formula of the five material constants contained in the tangent deformation modulus of the Duncan-Chang model were established, so the variation regular of the tangent deformation modulus with the FTC could be determined by unconsolidated undrained tests.
In the triaxial unconsolidated, undrained shear tests, no volumetric strain occurred in the samples. So, material constants G, F, and D that calculate tangential Poisson's ratio in the E-μmodel (initial tangent modulus -tangent Poisson's ratio model) could not be obtained, nor could the Kb and m of the E-B model (initial tangent modulus -bulk modulus model) be obtained. According to the generalized Hooke's law, the value of μt was 0.5 when the sample strain was 0. In order to ensure the stability of numerical calculation, Poisson's ratio is usually 0.49 [35].

Model Validation
By using Equations (5), (6), and (9), we could calculate hyperbolic model parameters a and b. According to Equation (6), the calculated results could be drawn in a stress-strain diagram, and we compared them with the experimental stress-strain curves. As shown in Figure 9, stress-strain curves that were determined by the modified Duncan-Chang model were, considering the FTC, basically close to the experiment stress-strain curves; the simulation effect was better in the stage when the strain was small and there was nearly elastic deformation. When the strain increased to 4%, there was a difference between some simulated and experimental curves. On the whole, it was reasonable to use this model to simulate the stress-strain curve of Nong'an saline soil affected by FTCs under the condition of an unconsolidated, undrained triaxial test.

Conclusions
In this paper, Nong'an saline soil from Western Jilin was taken as the research object. Through unconsolidated, undrained triaxial shear tests, the relationship between the mechanical properties of Nong'an saline soil and freeze-thaw cycles was studied, and a modified Duncan-Zhang constitutive model was established.
The stress-strain curves of the remolded saline soil were strain-hardening. There was linear elastic deformation in the initial loading stage, and it entered into the elastic-plastic deformation stage when strain reached 2%. With the increase of FTCs, (σ1 − σ3)f and c both decreased, rapidly decreasing after the first cycle and becoming stable between 30 and 60 cycles. φ increased at first, and then decreased. Using the exponential function could explain the relationship of (σ1 − σ3)f and c to FTC, and φ to FTC could be fitted by the Fourier function. Fitting results showed that the two functions reflected their variation rule well.
On the basis of Duncan-Chang model theory and regression analysis of the experimental data, a modified Duncan-Chang model was established for Nong'an saline soil, and the relationships of parameters Ei, K, N, and Rf to FTCs were analyzed. The analysis results indicated that Ei and K decreased rapidly at first and then became stable. Rf and N increased at first and then decreased. After 120 cycles, Rf and N all obviously increased. The relationships of Ei, K, and N to FTC could be explained by the exponential function, and the relationship of Rf to FTC could be explained by the rational function.
The rationality of the model was verified by comparing the simulation and experiment curves. It was found that, in the low-strain elastic-deformation stage, the model simulation results were consistent with the experiment results. When the strain was higher than 4%, there was little difference between simulated and experimental curves. Therefore, the Duncan-Chang model could well reflect the stress-strain relationship of the remolded Nong'an saline soil as a whole. Fitting curve 1 Fitting curve 5 Fitting curve 10 Fitting curve 30 Fitting curve 60 Fitting curve 120 Fitting curve

Conclusions
In this paper, Nong'an saline soil from Western Jilin was taken as the research object. Through unconsolidated, undrained triaxial shear tests, the relationship between the mechanical properties of Nong'an saline soil and freeze-thaw cycles was studied, and a modified Duncan-Zhang constitutive model was established.
The stress-strain curves of the remolded saline soil were strain-hardening. There was linear elastic deformation in the initial loading stage, and it entered into the elastic-plastic deformation stage when strain reached 2%. With the increase of FTCs, (σ 1 − σ 3 ) f and c both decreased, rapidly decreasing after the first cycle and becoming stable between 30 and 60 cycles. ϕ increased at first, and then decreased. Using the exponential function could explain the relationship of (σ 1 − σ 3 ) f and c to FTC, and ϕ to FTC could be fitted by the Fourier function. Fitting results showed that the two functions reflected their variation rule well.
On the basis of Duncan-Chang model theory and regression analysis of the experimental data, a modified Duncan-Chang model was established for Nong'an saline soil, and the relationships of parameters E i , K, N, and R f to FTCs were analyzed. The analysis results indicated that E i and K decreased rapidly at first and then became stable. R f and N increased at first and then decreased. After 120 cycles, R f and N all obviously increased. The relationships of E i , K, and N to FTC could be explained by the exponential function, and the relationship of R f to FTC could be explained by the rational function.
The rationality of the model was verified by comparing the simulation and experiment curves. It was found that, in the low-strain elastic-deformation stage, the model simulation results were consistent with the experiment results. When the strain was higher than 4%, there was little difference between simulated and experimental curves. Therefore, the Duncan-Chang model could well reflect the stress-strain relationship of the remolded Nong'an saline soil as a whole.