Effects of Drying–Wetting Cycle and Fines Content on Hysteresis and Dynamic Properties of Granite Residual Soil under Cyclic Loading

: In southern China, granite residual soil (GRS) is widely used as road base material. Thus, it is important to study the effects of hot and rainy climates and cyclic loads generated by trains on the dynamic properties of GRS. In this work, by means of dynamic triaxial tests, the effects of the number of drying and wetting (D–W) cycles, ﬁnes content and number of load cycles on the hysteresis curve, dynamic shear modulus G d and damping ratio λ of GRS are systematically investigated. The experimental results illustrate the changes in the morphology of the hysteresis curve and dynamic parameters with the numbers of load and D–W cycles, as well as the ﬁnes content. Namely, the area S , center offset d , and residual strain ε sp of the hysteresis curve decrease with the increase of load cycle number, increase with the growth of ﬁnes content, and ﬁrst decrease with the increase of D–W cycle number, then slowly increase to stabilized values. However, the major axis gradient k exhibits exactly the opposite relationships. Meanwhile, the dynamic shear modulus G d increases with the growing load cycles and decreases with the addition of ﬁnes content, and the damping ratio λ shows the opposite behavior. It is also shown that G d and λ vary linearly with respect to log N , where N is the number of D–W cycles. The dynamic properties of GRS are mostly affected by the number of load cycles, which is followed by the ﬁnes content and then the number of D–W cycles. The inﬂuence of the latter two factors on the dynamic properties of GRS may be primarily due to contact form changes between soil particles.


Introduction
Granite residual soil (GRS), a product of the in situ weathering of igneous parent rocks, is widely distributed in the humid tropical regions of southern China [1,2], especially in the Fujian province. GRS covers approximately three-fifths of areas in many cities such as Xiamen, Fuzhou, Quanzhou, etc. Since the essential mineral grains of GRS are quartz and kaolinite, while the secondary mineral grains are feldspar, chlorite and other ferric minerals, GRS is characterized by high mica content, loose structure, water sensitivity, and low cohesion [3][4][5]. GRS can satisfy most engineering conditions under dry conditions, but its mechanical and compressive properties change drastically after exposure to water. However, the degradation of the mechanical properties for GRS in the railway subgrade is strongly influenced by the rise of moisture content, which causes an interlayer of coarse and fine particles due to the release of contacts between the soil particles subjected to cyclic loading [6,7]. Owing to the subgrade interlayer, the reduction of the shear strength and compressibility of GRS results in severe damage to the roadbed, such as uneven settlement, track irregularity and mud pumping of the railway subgrade. so do the dynamic parameters. Chen et al. [36] investigated the effect of G d and λ of GRS under different confining pressure and stress amplitude and found that G d decreases with the increase of stress amplitude and λ inversely. In addition, increasing the surrounding pressure can suppress the decay of G d with dynamic stress. The results of Yin et al. [37] on the effect of confining pressure and dynamic stress amplitude on G d showed the same trends as those of Chen et al. [36]. Furthermore, Yin et al. [37] also found that the greater the water content in the soil, the faster G d decays.
The development of a dynamic constitutive model to quantitatively describe the dynamic deformation characteristics of soil under cyclic loading is a key issue [38]. The common soil models can be divided into three categories. The first type is the Masing model [39], which assumes that the hysteresis curve is obtained by doubling the skeleton curve, and defines the "upper skeleton curve criterion" and the "upper great circle criterion". The second type is the equivalent linear model, which assumes that the soil is a viscoelastic material, and introduces E d and λ as dynamic model parameters into the dynamic stress-strain model, for example, the Hardin-Drnevich model [40] and Ramberg-Osgood model [41], etc. The third type is the Iwan model [42], which describes the dynamic stress-strain relationship of the soil by combining different elastoplastic units either in series or in parallel. Except for the above three categories, there are also empirical models of E d with the number of load cycles, often denoted by N, confining pressure and deviator stress dynamic intrinsic relationship proposed by scholars. Wu et al. [43] and Xu et al. [33] established the empirical models for E d versus N in the study of soft clay soils, both in the form of E d = alnN + b. In addition, Xu et al. [33] correlated the fitted parameters a and b with the stress amplitude, confining pressure and loading frequency on this basis. Mu et al. [44] divided the stability zone, damage zone, and critical zone based on the law of accumulated deformation of undisturbed red clays under different confining pressures and stress amplitudes. Based upon this, the concept of equivalent dynamic elastic modulus E d-equal was proposed, and an empirical model was established. The model is E d-equal = alnN + b (stability zone a > 0, damage zone a < 0). Although this type of model is not directly associated with strain, it can also better reflect the dynamic property changes and can be associated with the studied factors, which is easy for practical application.
Previous works focused on the effect of confining pressure, loading frequency, or stress amplitude on the hysteresis curve and dynamic parameters for peaty soils, soft soils or GRS. In the meantime, the empirical models of E d with the number of load cycles are proposed according to the influencing factors for specific soils. It is noted that the D-W cycles and fine particles have a large influence on the dynamic characteristics of GRS, and the process of D-W cycles will affect the change of fines content. However, few studies describe in detail the dynamic properties and explain the dynamic mechanism of GRS under the D-W cycles and fines content.
In view of the phenomenon of pulping and mudging under the train cyclic dynamic load in a hot and rainy climate, this work focuses on the dynamic characteristics of GRS under cyclic loading by D-W cycles and fines content through hysteresis curves and further reveals the mechanism underlying the effect of fine particles and D-W cycles on the dynamic characteristics of GRS. Through the experimental data analysis, the characteristics of the hysteresis curve and the variation of parameters are revealed. It supports the comprehensive understanding of the dynamic properties of GRS. In addition, the dynamic shear modulus and dynamic damping ratio models established based on the dynamic triaxial experimental data can provide the prediction of dynamic parameters of residual soil considering the influence of D-W cycles and fines content. It is helpful to carry out dynamic analysis of railway subgrade stability and provide theoretical support for the construction of the railroad subgrade in the residual soil area.

Test Materials
The GRS used in this study was extracted from the edge of a railway subgrade, representing typical soils therein. The GRS was dried for 48 h, pulverized and then sieved through a 5.0 mm sieve to generate the prepared soils.
The specific soil gravity (G s ) is 2.72. The particle size distribution of the soil obtained from sieve analysis exhibits a composition of 37.2% fines (<0.075 mm), 57.3% sand (0.075~4.75 mm), and 5.5% Gravel grain (>4.75 mm). The values of liquid limit (LL), plastic limit (PL) and PI are 46.6%, 23.8% and 22.8%, respectively. The maximum dry density and optimal moisture content (OMC) of the soil are 1.65 g/cm 3 and 20.3%. According to the Criteria for classification of soils (GBJ145-90) [45], the soil is classified as residual soil. In order to study the influence of fines content, the soil is divided into seven particle groups (5-2 mm, 2-1 mm, 1-0.5 mm, 0.5-0.25 mm, 0.25-0.1 mm, 0.1-0.075 mm and <0.075 mm) by the sieving method. In order to meet the specification requirements for high-speed railroad bed fillers [46], the fines content varies from 0%, 10%, 20%, to 30% by adding fine particles (particle size < 0.075 mm). The percentages refer to the ratio of the dry mass of fine particles to the dry mass of soil. The Standard Proctor tests were conducted to obtain the optimal moisture content and maximum dry density of the adjusted soils. The particle size distributions and physical properties of soils with different fines contents are illustrated in Table 1.

Preparation of GRS Samples
The GRS samples were prepared through the standard procedure defined by the Specification of Soil Tests (GB/T50123-1999) [47]. The distilled water was added to the prepared soils to reach the optimal moisture content of each group, and the soil mixtures were kept in the enclosed bags without evaporation for 48 h to ensure the water could be uniformly distributed. Subsequently, the soil triaxial samples were compacted in five equal layers on the bottom platen inside a split mold with a height of 100 mm and a diameter of 50 mm. The weight and height of each layer were measured to obtain a homogeneous dry density of 95% of ρ dmax . The acceptable error for quality variation within the same set of soil samples is set at less than 1%. If the variation exceeds 1%, the soil samples will be re-prepared. It ensures the consistency and reliability of the soil samples.

Drying-Wetting Cycle Process
To explore the effect of the D-W cycles on the dynamic properties of GRS in the railway subgrade, the soil samples subjected to D-W cycles were used to conduct the dynamic triaxial tests. According to the environmental parameters, the severe D-W cycles of the railway subgrade of Xiamen in Fujian province caused by the summer monsoon climate conditions were simulated, which included the highest temperature of about 55 • C (that is close to that at the ground surface in summer) and the moisture content of approximately 10% (which is equal to the moisture content of subgrade in-situ). The prepared soil samples in OMC were all wetted to saturation and were then subjected to the D-W cycles. Each D-W cycle includes three phases, namely, drying, wetting and saturation. The first D-W cycle process for soil samples involves three steps: initial moisture content to saturation, drying to a 10% moisture content state, and then saturation again. The subsequent D-W cycles of a soil sample are considered to have undergone a D-W cycle as long as it dries from the saturated state to 10% moisture content and then wets to saturation. To find an optimum maximum number of D-W cycles, a series of preliminary tests were conducted with respect to 5, 7, 9, and 11 D-W cycles under the same testing conditions. The results indicated that their stress-strain curves of 7, 9, and 11 D-W cycles are very close. Thus, seven D-W cycles were used as the maximum D-W cycles from a time-saving point of view. In the present study, GRS samples with different D-W cycles were prepared via the following procedures (see Figure 1): Appl. Sci. 2023, 13, 6660 5 of 26 to saturation, drying to a 10% moisture content state, and then saturation again. The subsequent D-W cycles of a soil sample are considered to have undergone a D-W cycle as long as it dries from the saturated state to 10% moisture content and then wets to saturation.
To find an optimum maximum number of D-W cycles, a series of preliminary tests were conducted with respect to 5, 7, 9, and 11 D-W cycles under the same testing conditions. The results indicated that their stress-strain curves of 7, 9, and 11 D-W cycles are very close. Thus, seven D-W cycles were used as the maximum D-W cycles from a timesaving point of view. In the present study, GRS samples with different D-W cycles were prepared via the following procedures (see Figure 1): (1) Drying process. The soil samples were put into the oven at 55 °C for about 6 h (to simulate real drying field conditions), then taken out and left for 24 h to cool down naturally to room temperature. Finally, the samples were weighed to determine the change in moisture content over time. The trial test results showed that the change of moisture content is less than 0.5%, and the moisture content remains around 10%, indicating that the samples at 10% moisture content are in completely dry conditions.
(2) Wetting process. The samples were first pasted by the wet filter papers, then sprayed uniformly with water every 10 min using an electric sprayer. The wetting process is continuously monitored for soil sample quality and moisture content, and it is terminated when neither changes.
(3) Saturation process. The saturation process includes vacuum saturation and backpressure saturation. In vacuum saturation, triaxial saturators filled with the soil samples were placed in a barrel for vacuum saturation. The triaxial saturator is shown in Figure 1. Each triaxial saturator had a permeable stone at each end, and the specimen was placed between two permeable stones. During the saturation process, the specimen was saturated by exhaust and water absorption through two permeable stones. The procedure was as follows: First, the soil samples were vacuumed in the saturator with a vacuum pump for 2 h, then the vacuum pump valve was closed, and the water filling valve was opened to add water. When the water level reached 3/4 of the soil sample height, the valve was closed, and the vacuum pump was opened for 2 h. Then the soil samples were kept under a vacuum of one standard atmosphere pressure for 24 h. The back-pressure saturation was used to install the soil samples in a dynamic triaxial instrument by applying back pressure (1) Drying process. The soil samples were put into the oven at 55 • C for about 6 h (to simulate real drying field conditions), then taken out and left for 24 h to cool down naturally to room temperature. Finally, the samples were weighed to determine the change in moisture content over time. The trial test results showed that the change of moisture content is less than 0.5%, and the moisture content remains around 10%, indicating that the samples at 10% moisture content are in completely dry conditions.
(2) Wetting process. The samples were first pasted by the wet filter papers, then sprayed uniformly with water every 10 min using an electric sprayer. The wetting process is continuously monitored for soil sample quality and moisture content, and it is terminated when neither changes.
(3) Saturation process. The saturation process includes vacuum saturation and backpressure saturation. In vacuum saturation, triaxial saturators filled with the soil samples were placed in a barrel for vacuum saturation. The triaxial saturator is shown in Figure 1. Each triaxial saturator had a permeable stone at each end, and the specimen was placed between two permeable stones. During the saturation process, the specimen was saturated by exhaust and water absorption through two permeable stones. The procedure was as follows: First, the soil samples were vacuumed in the saturator with a vacuum pump for 2 h, then the vacuum pump valve was closed, and the water filling valve was opened to add water. When the water level reached 3/4 of the soil sample height, the valve was closed, and the vacuum pump was opened for 2 h. Then the soil samples were kept under a vacuum of one standard atmosphere pressure for 24 h. The back-pressure saturation was used to install the soil samples in a dynamic triaxial instrument by applying back pressure to monitor the pore water pressure coefficient B, and when the value of B reached above 95%, the saturation process finished.

Dynamic Triaxial Test
The dynamic triaxial test was carried out using the British GDS (Global Digital Systems) standard dynamic triaxial test system. A static load of 15 kPa σ s was applied to simulate the action of the superstructure. Considering the effect of D-W cycles and the limited depth of the train load, the confining consolidation pressure σ 3 was set to 60 kPa.
Generally, the dynamic loads of high-speed trains are between 30 kPa and 70 kPa [48], and the loading frequency is between 1 Hz and 5 Hz, with actual data measured as 1.39-1.85 Hz [49]. Therefore, the dynamic load's amplitude σ d was set to 60 kPa, and the loading frequency was set to 2 Hz. The experiment was over as soon as the axial strain of the soil sample exceeded 10% or the number of vibration cycles reached 10,000 [50].

Hysteresis Curve Characteristics of GRS
The hysteresis curve refers to the dynamic stress-strain relationship throughout a load cycle. Its morphological characteristics can be used to describe the dynamic deformation, viscosity, stiffness change, energy loss, and other significant dynamic aspects of soil. The shape of the hysteresis curve can be affected by the confining pressure, load frequency, load amplitude, stress history, and naturally occurring environmental changes in the soil. Figure 2 illustrates the trend of the hysteresis curve variations for different fines contents without D-W cycles after 500 cyclic loadings. The corresponding observations are summarized as follows: (1) The hysteresis curve changes from the initial shuttle-like to a willow shape as the number of load cycles increases. (2) The largest hysteresis loop lies in the first cycle, and it gradually becomes smaller and denser in the following cycles. GRS exhibits an elastoplastic behavior, and the irreversible plastic strain gradually enlarges with increasing cyclic numbers. (3) As the fines content increases, the axial strain increases significantly. Comparing the axial strain of C0 with that of C30, the latter is nearly three times the former. The reasons behind these observations could be that the fine particles play a lubricating role on the surface of the coarse particles when the fines content increases from 0% to 10%. Moreover, when the fines content further increases to 30%, the fine particles may block the internal drainage channel of pores, resulting in the accumulation of excess pore water pressure. As the pore water pressure accumulates to the limit, the fine particles will be pushed out, and the axial strain will increase. Consequently, the axial strain becomes greater and greater with the increase of the fines content. Figure 3 plots the hysteresis curve changes of soil samples without adding fine particles for the 1st, 3rd, 5th, and 7th D-W cycles under 500 loading cycles. It can be observed that the hysteresis curve changes show a non-monotonic pattern. The axial strain decreases first and then increases significantly when the number of D-W cycles increases. With the results of Dou et al. [51] and Fu et al. [52], it is inferred that such a result in this paper is probably due to the new skeleton formed by fine and coarse particles being denser within three D-W cycles and the axial strain decreases. However, when the number of D-W cycles continuously increases to five, the hydrophilic capacity of clay minerals decreases, and the adhesion between particles weakens. At the same time, the soil structure becomes loose under the combined action of dynamic load and D-W cycles, which eventually leads to significant growth of axial strain. Appl

Hysteresis Curve Parameters
To a certain extent, the hysteresis curve's evolution law reveals the soil's internal structure change rule. By extracting the parameters of the shape change from the hysteresis curve, the influence of external conditions on soil damage can be obtained. This study uses the hysteresis curve area S, the slope of the long axis k, the center offset of hysteresis curve d and residual strain ε sp to analyze the hysteresis curve characteristics of GRS. The calculation methods of S, k, d and ε sp and their variations will be introduced separately below.
(1) Hysteresis curve area S The variation of S normally indicates the degree of soil compaction. A large S shows that the soil structure is loose and easily compacted and that a large amount of energy is lost during compaction. On the contrary, a small S indicates that the soil structure is dense and difficult to compact, and less energy is lost during compaction.
The area of the arc from the initial point of loading, to the final point of unloading is the area of the hysteresis curve. The hysteresis curve area S is expressed as follows: where σ di and ε i denote the deviator stress and the axial strain under the ith loading cycle.
(2) Major axis gradient k The stiffness and elasticity modulus of the soil can be reflected by the major axis gradient k of the hysteresis curve. When k grows, the soil stiffness and elasticity modulus increase and vice versa. The major axis gradient k is calculated by: where σ N,max , ε N,max , σ N,min and ε N,min are the maximum and the minimum stress and strain of the Nth hysteresis curve, respectively.
(3) Center offset d The center offset d is the distance between the centers of the previous and subsequent hysteresis curves, which reflects the soil deformation and internal structural damage. The larger the d, the greater the irrecoverable deformation of the soil sample and the greater the damage to the internal structure of the soil, and vice versa. Since the graph enclosed by the hysteresis curve is not a regular shape, it is not trivial to pinpoint the exact location of the center point. This work primarily studies the trend of its center distance with loading and D-W cycles. The difference between the strain maximum of the previous hysteresis curve and the latter hysteresis curve is proportional to the center distance. Therefore, the center offset d is computed as follows: where ε N,max and ε N−1,max represent the maximum strain of the Nth and the (N − 1)th hysteresis curve, respectively.

(4) Residual strain ε sp
The residual strain is the strain difference between the end point of loading and the starting point in one cycle, namely, the degree of non-closure or the size of the opening of the hysteresis curve. Thus, the residual strain ε sp is given as follows: where ε s and ε e are the strain values corresponding to the starting and ending points of the Nth loading, respectively.

Effect of Number of Loading Cycles on the Hysteresis Curve Parameters
The morphology of the hysteresis curve of GRS changes dramatically with the increase of loading cycles, resulting in significant changes in its morphological parameters. As shown in Figure 4, it is observed that the higher the number of loads, the smaller the area and the degree of opening, while the larger the degree of tilt.

Effect of Number of Loading Cycles on the Hysteresis Curve Parameters
The morphology of the hysteresis curve of GRS changes dramatically with the increase of loading cycles, resulting in significant changes in its morphological parameters. As shown in Figure 4, it is observed that the higher the number of loads, the smaller the area and the degree of opening, while the larger the degree of tilt.   [32,53]. The ratio of the decrease under the first 20 loads to the total decrease can be found in Table 2. The average reduction ratio for S, d and εsp are 77.8%, 92.5% and 94.8%, respectively, with a maximum of 96.6% and a minimum of 72.9%. (2) Unlike the other three parameters, as shown in Figure 5b, k increases sharply in the early load stage and then increases slowly, and moreover, k approximately changes linearly with LogN.    [32,53]. The ratio of the decrease under the first 20 loads to the total decrease can be found in Table 2. The average reduction ratio for S, d and ε sp are 77.8%, 92.5% and 94.8%, respectively, with a maximum of 96.6% and a minimum of 72.9%. (2) Unlike the other three parameters, as shown in Figure 5b, k increases sharply in the early load stage and then increases slowly, and moreover, k approximately changes linearly with LogN. The reason for these observations might be that in the early stage of loading (N ≤ 20), the soil is rapidly compacted and the deformation of the soil is predominantly plastic [54]. At this stage, the internal structural damage and residual strain of the soil are most obvious, while in the late stage of loading, the soil compactness is enhanced with filled pores and increasing contact between soil particles, and the soil deformation gradually transitions from plastic to elastic, so the residual strain decreases and the soil stiffness gradually increases. The reason for these observations might be that in the early stage of loading (N ≤ 20), the soil is rapidly compacted and the deformation of the soil is predominantly plastic [54]. At this stage, the internal structural damage and residual strain of the soil are most obvious, while in the late stage of loading, the soil compactness is enhanced with filled pores and increasing contact between soil particles, and the soil deformation gradually transitions from plastic to elastic, so the residual strain decreases and the soil stiffness gradually increases.

Effect of Number of Fines Content on the Hysteresis Curve Parameters
The fines content in the GRS affects the shape of the hysteresis curve. Figure 6 depicts the variation of four hysteresis curve parameters with the fines content. According to the amount of fines content, it is divided into three ranges, i.e., low fines content (C0 to C10, Range I), medium fines content (C10 to C20, Range II) and high fines content (C20 to C30, Range III). When the fines content increases, S, d, and εsp of the hysteresis curve increase. They increase quickly in Range I and Range III but slowly in Range II. Additionally, the influence of fines content is greater at low-load cycles than at high-load cycles. Taking S as an example, the values of S under cyclic load are listed in Table 3. Under low load cycles (N = 1-20), the average area of SC30 is 6.04, 2.5 times that of SC0. However, the average area of SC30 is only 1.4 times that of SC0 under high-loading cycles (N = 1000-10000). The soil group of C0 has the maximum k, and the increase of fines content makes k decrease. k has a noticeable decrease in Range I and Range III, but it has a small change at low load cycles and fluctuates at high loading cycles in Range II.

Effect of Number of Fines Content on the Hysteresis Curve Parameters
The fines content in the GRS affects the shape of the hysteresis curve. Figure 6 depicts the variation of four hysteresis curve parameters with the fines content. According to the amount of fines content, it is divided into three ranges, i.e., low fines content (C0 to C10, Range I), medium fines content (C10 to C20, Range II) and high fines content (C20 to C30, Range III). When the fines content increases, S, d, and ε sp of the hysteresis curve increase. They increase quickly in Range I and Range III but slowly in Range II. Additionally, the influence of fines content is greater at low-load cycles than at high-load cycles. Taking S as an example, the values of S under cyclic load are listed in Table 3. Under low load cycles (N = 1-20), the average area of S C30 is 6.04, 2.5 times that of S C0 . However, the average area of S C30 is only 1.4 times that of S C0 under high-loading cycles (N = 1000-10,000). The soil group of C0 has the maximum k, and the increase of fines content makes k decrease. k has a noticeable decrease in Range I and Range III, but it has a small change at low load cycles and fluctuates at high loading cycles in Range II.

Effect of D-W Cycles on the Hysteresis Curve Parameters
The morphological parameters for four groups of soils with different fines contents soil under seven D-W cycles were plotted in Figures 7-10 to illustrate in detail the effect of D-W cycles. Since the compaction effect of soil is obvious in the first 20 loading cycles, these parameters change do not fully reflect soil properties. In addition, excessive loading numbers also affect soil properties (see the green and grey areas in Figure 7a,b), and the center offset d and residual strain ε sp mainly occurs at N < 100, after which they are close to 0. Therefore, N = 20 to 90 (blue area) is the load-interval that best represents the influence of the D-W cycle on soil properties.

Effect of D-W Cycles on the Hysteresis Curve Parameters
The morphological parameters for four groups of soils with different fines contents soil under seven D-W cycles were plotted in Figures 7-10 to illustrate in detail the effect of D-W cycles. Since the compaction effect of soil is obvious in the first 20 loading cycles, these parameters change do not fully reflect soil properties. In addition, excessive loading numbers also affect soil properties (see the green and grey areas in Figure 7a,b), and the center offset d and residual strain εsp mainly occurs at N < 100, after which they are close to 0. Therefore, N = 20 to 90 (blue area) is the load-interval that best represents the influence of the D-W cycle on soil properties.   From Figures 7, 8, 9 and 10b, it can be observed that the k fluctuates at T0~3, while there is a significant decrease at T3~5, followed by a gradual stabilization at T5~7. Compared to S, d and ε sp , k shows an opposite variation as the D-W cycles increase. The major axis gradient k ave under different D-W cycles are calculated and listed in Table 4. The value of k ave first increases and then decreases. Usually, the range of k ave is larger at the later load stage than that at the early load stage. Appl      From Figures 7-10b, it can be observed that the k fluctuates at T0~3, while there is a significant decrease at T3~5, followed by a gradual stabilization at T5~7. Compared to S, d and εsp, k shows an opposite variation as the D-W cycles increase. The major axis gradient kave under different D-W cycles are calculated and listed in Table 4. The value of kave first increases and then decreases. Usually, the range of kave is larger at the later load stage than that at the early load stage. The above variation may be explained as follows. As shown by the schematic diagram of the initial state of GRS in Figure 11a, the grain sizes of GRS are mainly between 0.5 mm and 0.075 mm, and the total amounts of coarse and fine particles are close, while   N3  N1  N3  N1  N3  N1  N3  N1  N3  N1  N3  N1  N3  N1  N3  N1  N3   C0  397 652  399  614  429  620  391  580  356  530  267  472  275  436  267  418  162  233  C10  277 456  291  480  276  442  268  431  285  432  229  392  239  380  216  352  75  127  C20  280 535  301  463  230  409  333  491  290  468  237  358  230  390  265  392  102  176  C30  189 457  213  392  220  408  222  377  220  392  209  359  215  386  220  364  33  97 Note: k ave means the average of k under different numbers of loads. N1 means N = 1, 2, 3, . . . , 9; N3 means N = 1000, 2000, 3000, . . . , 9000.
The above variation may be explained as follows. As shown by the schematic diagram of the initial state of GRS in Figure 11a, the grain sizes of GRS are mainly between 0.5 mm and 0.075 mm, and the total amounts of coarse and fine particles are close, while the content of medium particles is quite low. Due to the small content of medium particles, the macrospores are larger and are filled by fine particles and pore water. At the beginning of D-W cycles (T2~3, see Figure 11b), the fine particles are dispersed, and the cementitious substances in the fine particles are dissolved in the water and partly lost with the migration of water [1,55]. However, during the drying process, the fine particles re-combine and aggregate to form a denser soil skeleton [56], resulting in a soil energy loss, a reduction of plastic deformation and an increase in stiffness. In the middle of D-W cycles (T3~5, see Figure 11c), the intrusion of water makes fine particles expand and disperse into finer particles, weakening the adhesion among the fine particles in the wetting process. On the other hand, finer particles continue to be lost with the water migration in the drying process [57]. Ultimately, the soil structure gradually loosens, so it shows more +energy loss and plastic strain characteristics. At the later D-W process (T5-7, see Figure 11d), some coarse particles are broken into fine particles after undergoing the previous D-W cycles, and with the further dispersion of fine particles, the fine contents filled in the pores increase.
More fine particles wrap the coarse particles, resulting in a more stable skeleton structure. As a result, the variation of parameters is stabilized.
plastic deformation and an increase in stiffness. In the middle of D-W cycles (T3~5, see Figure 11c), the intrusion of water makes fine particles expand and disperse into finer particles, weakening the adhesion among the fine particles in the wetting process. On the other hand, finer particles continue to be lost with the water migration in the drying process [57]. Ultimately, the soil structure gradually loosens, so it shows more +energy loss and plastic strain characteristics. At the later D-W process (T5-7, see Figure 11d), some coarse particles are broken into fine particles after undergoing the previous D-W cycles, and with the further dispersion of fine particles, the fine contents filled in the pores increase. More fine particles wrap the coarse particles, resulting in a more stable skeleton structure. As a result, the variation of parameters is stabilized. Figure 11. Effect mechanism of D-W cycles. Table 5 tabulates the influence of the number of load cycles, fines content and the number of D-W cycles on each parameter. It can be observed that the number of load cycles has the strongest influence on the four parameters, while the number of D-W cycles has the weakest influence on the parameters. The variation of ΔSave by the number of load cycles is greater than 1.39, while the variation of ΔSave by fines content and D-W cycles is less than 1.39 in most cases, especially in N3, where the variation is all less than 0.94 (see column ΔSave' in Table 5). The number of D-W cycles has the least influence on the parameters. For example, Δεspave changes by fines content reach 4% in more than half of the cases, while Δεspave changes by D-W cycle periods surpass 4% only once.  Table 5 tabulates the influence of the number of load cycles, fines content and the number of D-W cycles on each parameter. It can be observed that the number of load cycles has the strongest influence on the four parameters, while the number of D-W cycles has the weakest influence on the parameters. The variation of ∆S ave by the number of load cycles is greater than 1.39, while the variation of ∆S ave by fines content and D-W cycles is less than 1.39 in most cases, especially in N3, where the variation is all less than 0.94 (see column '∆S ave ' in Table 5). The number of D-W cycles has the least influence on the parameters. For example, ∆ε spave changes by fines content reach 4% in more than half of the cases, while ∆ε spave changes by D-W cycle periods surpass 4% only once.

Summary
In addition, there is a coupling between the factors that will change the degree of influence of the parameters. For example, the degree of influence of the loading cycles on four parameters is the greatest at T0-C30 (∆S ave = −5.3, ∆k ave = 267.9, ∆d ave = −8.05, ∆ε spave = −9.82), the degree of influence of fine particle contents and number of D-W cycles on S, d, ε sp is the least at N1-T5 (∆S ave = 0.89, ∆d ave = 1.65%, ∆ε spave = 1.91%) and N1-C10 (∆S ave = 1.02, ∆d ave = 1.88%, ∆ε spave = 2.25%), respectively. In addition, there is a significant decrease in the degree of influence of the number of loading cycles after one D-W cycle.

Dynamic Parameters of GRS
(1) Dynamic shear modulus G d Since the dynamic shear modulus (G d ) can not be directly determined in the dynamic triaxial tests, it is usually calculated indirectly, as follows: where E d is elasticity modulus, µ is Poisson's ratio. According to [58], the Poisson's ratio of GRS ranges from 0.3 to 0.33, thus µ is taken as 0.3 in this study. The dynamic elasticity modulus (E d ) is obtained from the dynamic triaxial test. It can be obtained as: where σ A , σ B and ε A , ε B are the stress and strain values corresponding to the intersection points of A and B, respectively.
(2) Damping ratio λ The damping ratio is the ratio of the dissipated energy and the elastic strain energy in one load cycle. According to the Standard for geotechnical testing method (GB/T 50123-2019) [59], the formula is given as follows: where A z is the area of the hysteresis curve and A s is the area of the triangle.

Characteristics of Dynamic Parameters of GRS
By taking soil that has not undergone any D-W cycle as an example, Figure 12 depicts the dynamic shear modulus G d and damping ratio λ with the number of load cycles for GRS. It is easy to see that the variation of dynamic parameters with the number of loading cycles has two stages, i.e., the initial loading stage (N < 20) and the late loading stage (N ≥ 20). In the initial loading stage (N < 20), with the increases in the number of loads, the dynamic shear modulus G d grows rapidly while the damping ratio λ decreases rapidly. It exhibits obvious strain-hardening characteristics [60]. This may be due to the significant compaction effect of soil, and the soil deformation is dominated by plastic deformation. Therefore, the change of dynamic parameters at this stage does not reflect the soil characteristics. In the late loading stage (N ≥ 20), both the increasing rate of dynamic shear modulus and the decreasing rate of damping ratio slows down and afterward reach stabilization (N > 1000). The dynamic shear modulus growth and the damping ratio reduction at this stage are mainly attributable to the particle dislocation and pore compression in the soil, resulting in pore water pressure and further compaction of the soil. The soil enters the elastoplastic deformation stage. When the number of load cycles continues to increase (N > 1000), the soil becomes more dense, and the soil deformation is dominated by elastic deformation. At this stage, dynamic parameters tend to be fixed. Therefore, the change of dynamic parameters at 20~1000 loading cycles is analyzed below.
where A  , B  and A  , B  are the stress and strain values corresponding to the intersection points of A and B, respectively.
(2) Damping ratio λ The damping ratio is the ratio of the dissipated energy and the elastic strain energy in one load cycle. According to the Standard for geotechnical testing method (GB/T 50123-2019) [59], the formula is given as follows: where Az is the area of the hysteresis curve and As is the area of the triangle.

Characteristics of Dynamic Parameters of GRS
By taking soil that has not undergone any D-W cycle as an example, Figure 12 depicts the dynamic shear modulus Gd and damping ratio λ with the number of load cycles for GRS. It is easy to see that the variation of dynamic parameters with the number of loading cycles has two stages, i.e., the initial loading stage (N < 20) and the late loading stage (N ≥ 20). In the initial loading stage (N < 20), with the increases in the number of loads, the dynamic shear modulus Gd grows rapidly while the damping ratio λ decreases rapidly. It exhibits obvious strain-hardening characteristics [60]. This may be due to the significant compaction effect of soil, and the soil deformation is dominated by plastic deformation. Therefore, the change of dynamic parameters at this stage does not reflect the soil characteristics. In the late loading stage (N ≥ 20), both the increasing rate of dynamic shear modulus and the decreasing rate of damping ratio slows down and afterward reach stabilization (N > 1000). The dynamic shear modulus growth and the damping ratio reduction at this stage are mainly attributable to the particle dislocation and pore compression in the soil, resulting in pore water pressure and further compaction of the soil. The soil enters the elastoplastic deformation stage. When the number of load cycles continues to increase (N > 1000), the soil becomes more dense, and the soil deformation is dominated by elastic deformation. At this stage, dynamic parameters tend to be fixed. Therefore, the change of dynamic parameters at 20~1000 loading cycles is analyzed below.  The effect of fines content on the dynamic parameters of GRS is closely related to its water content. When the soil is in the saturated state, the larger the fines contents, the smaller the dynamic shear modulus and the higher the damping ratio [61]. Figure 13 plots the evolution of dynamic parameters with the variation of fines contents of GRS. The dynamic shear modulus shows a trend of decreasing, then increasing, and finally decreasing again with the increase of fines content (see Figure 13a). The damping ratio generally shows an increasing trend with the increase of fines content, and the growth rate slows down as the number of load cycles increases (see Figure 13b). modulus without D-W cycles, the dynamic shear modulus after several D-W cycles exhibit a decreasing trend (see Table 7, GdT7 < GdT0).
The relationship between the damping ratio of GRS with D-W cycles is depicted in Figure 15. Compared with the influence of the number of cyclic loading on the damping ratio, the influence of D-W cycles on it is negligible. Table 8. shows the range of damping ratios among seven D-W cycles. It is observed that the differences are very small. It can be concluded that the damping ratio of GRS is less affected by D-W cycles.   Figure 14). It is observed that the influence of the number of D-W cycles on the dynamic shear modulus is greater when the fine particle content is low (see Figure 14a-d). Table 6 shows the range of dynamic shear modulus among seven D-W cycles. The range G d for C0 is nearly three or five times that for C30. Meanwhile, compared with the dynamic shear modulus without D-W cycles, the dynamic shear modulus after several D-W cycles exhibit a decreasing trend (see Table 7, G dT7 < G dT0 ).        T0  T7  T0  T7  T0  T7  T0  T7  T0  T7  T0  T7  T0  The relationship between the damping ratio of GRS with D-W cycles is depicted in Figure 15. Compared with the influence of the number of cyclic loading on the damping ratio, the influence of D-W cycles on it is negligible. Table 8. shows the range of damping ratios among seven D-W cycles. It is observed that the differences are very small. It can be concluded that the damping ratio of GRS is less affected by D-W cycles.

Dynamic Parameter Models of GRS
Some scholars have proposed formulas for fitting the dynamic elastic modulus E d and the number of loading cycles N [33,44]. The form of the proposed formulas is the same (see Equation (8)). Mu et al. [44] established the relationship between the dynamic elastic modulus E d and the confining pressure σ 3 . Xu et al. [33] explored the effects of frequency f, dynamic load amplitude σ d , initial partial stress σ s , and confining pressure σ 3 on E d .
where A and B are fitting parameters, and N is the number of load cycles. The model fitting parameters proposed above are related to the load and not to the soil properties. Therefore, the two models are presented to predict the dynamic shear modulus and damping ratio considering soil properties (D-W cycles and fine particle content).
(1) Dynamic shear modulus model From Figure 16, it can be observed that the relationship between the dynamic shear modulus and the number of load cycles (N > 20) can be approximated by a straight line through the transformation of . Therefore, all the data in Figure 16 are fitted according to the following formula: where G d20 and G d1000 are the dynamic shear modulus after 1 and 1000 loading cycles under C0 and T0. 'a G ' and 'b G ' are the fitting parameters and N is number of load cycles. The fitting results are shown in Figure 16 and Table 9.  It is important to highlight that prior to utilizing Equation (9), a preliminary test needs to be conducted to determine the dynamic shear modulus for N = 20 and N = 1000 without the DW cycle of GRS.
(2) Damping ratio model It is important to highlight that prior to utilizing Equation (9), a preliminary test needs to be conducted to determine the dynamic shear modulus for N = 20 and N = 1000 without the DW cycle of GRS.
(2) Damping ratio model As can be observed from Figure 17, the damping ratio and the number of loadings follow an exponential relationship, so the fitted equation is defined as: where a λ and b λ are the fitting parameters and N is the number of loading cycles. The values fitted from equation and the laboratory data are represented in Figure 17. Table 10 describes the values for the parameters as well as the determination coefficients R 2 . In general, the experimental data are close to the continuous lines of the models. In this sense, both the dynamic shear modulus and damping ratio models are very accurate, since in all the cases the determination coefficients (R 2 ) yield a result of over 94% (Tables 9 and 10). However, the coupling effect of fines content and D-W cycles are obvious, it is difficult to determine the relationship between the fitting parameters, D-W cycles, and the fines content, but the laws presented in the model can be applied to engineering. Table 10 describes the values for the parameters as well as the determination coefficients R 2 . In general, the experimental data are close to the continuous lines of the models. In this sense, both the dynamic shear modulus and damping ratio models are very accurate, since in all the cases the determination coefficients (R 2 ) yield a result of over 94% (Tables 9 and  10). However, the coupling effect of fines content and D-W cycles are obvious, it is difficult to determine the relationship between the fitting parameters, D-W cycles, and the fines content, but the laws presented in the model can be applied to engineering.

Conclusions
The dynamic behaviors of GRS in Southern China were experimentally studied. The effects of load cycles, fines content, and D-W cycles were investigated through the dynamic triaxial tests, where four fines contents and seven D-W cycles were considered. The main conclusions are as follows: (1) The morphology of the hysteresis curve changes with the numbers of load and D-W cycles, as well as the fines content. The hysteresis curve area S, center offset d, and residual strain ε sp of GRS vary similarly with the increase of load cycles. They all decrease sharply in the early stage(N < 20), and the reduction rate slows down in the later load stage (20 ≤ N ≤ 1000) and finally stabilizes in the end(N > 1000). The major axis gradient k increases with the number of load cycles and approximately changes linearly with log N. At the same time, as the fines content increases, S, d and ε sp increase and k decreases. All of these properties fluctuate with the increase of D-W cycles and tend to be stable after about seven D-W cycles.
(2) The growth of dynamic shear modulus and the reduction of damping ratio with the number of load cycles can be divided into the rapid change stage(N < 20) and the slow change stage (20 ≤ N ≤ 1000). Moreover, the dynamic shear modulus and damping ratio show a good linear relationship with the logarithmic number of load cycles and two prediction equations were obtained by the data fitting. With the increase of fines content, the dynamic shear modulus shows a trend of decreasing while the damping ratio shows an increasing trend. The dynamic shear modulus changes with D-W cycles, and it decreases after several D-W cycles. However, the damping ratio is less affected by D-W cycles.
(3) The dynamic characteristics of GRS are influenced by the number of load cycles, fines content, as well as the number of D-W cycles. Among them, the number of load cycles has the dominant influence and the number of D-W cycles has the least influence. The influence of fines content is more pronounced for the low number of load cycles than for the high number of load cycles. In addition, the effects of fines content and the number of D-W cycles are not independent.
(4) Fitting models are proposed to estimate the dynamic shear modulus and damping ratio based on the number of load cycles. These models can predict dynamic shear modulus and damping ratio considering the effects of fines content and D-W for dynamic stability analysis of railway subgrade.

Conflicts of Interest:
The authors declare that there is no conflict of interest regarding the publication of this paper.