Effects of Cyclic Freezing–Thawing on Dynamic Properties of Loess Reinforced with Polypropylene Fiber and Fly Ash

: Periodic freezing–thawing is recognized as a real threat to the mechanical properties of reinforced loess, which has been used in the recent construction of high-speed railways in northwest China; however, the performance of these materials under periodic freezing–thawing and dynamic loading has rarely been investigated. In this work, dynamic triaxial tests were conducted on ﬂy ash- and polypropylene ﬁber-reinforced loess with different blend ratios and freeze–thaw cycles. The dynamic shear modulus and damping ratio were investigated. The results revealed that cyclic freezing–thawing had a remarkable effect on the dynamic shear modulus and damping ratio, which demonstrated considerable reductions and increases, respectively, after cyclic freezing–thawing. Additionally, the dynamic shear modulus increased notably with the ﬂy ash content and conﬁning pressure and decreased with the water content. Meanwhile, the damping ratio increased with the ﬁber content and water content and decreased with the ﬂy ash content and conﬁning pressure. Com-paratively, the effects of polypropylene ﬁber on dynamic behavior were found to be not signiﬁcant. Furthermore, novel models were established to predict the dynamic shear modulus and damping ratio for reinforced loess. The results provide more information towards infrastructure design in seasonal frozen regions. and damping ratio were established by taking account of freeze–thaw cycle, FA content, PF content, initial water content, and confining pressure. The comparison of the experimental and predicted results illustrated that the newly established models are suitable to estimate the dynamic shear modulus and damping ratio of fly ash–polypropylene fiber-reinforced loess.


Introduction
Loess is widespread soil that stretches from seasonal frozen areas to permafrost areas in Northern China. Due to a lack of coarse-grained materials, it is often selected as a filling material for highways, railways, and high-speed railways (HSR) despite its unfavorable engineering characteristics such as its large pores and high sensitivity to water and temperature [1][2][3]. An effective method to mitigate these deficiencies and to improve the fill is to reinforce the soil with cementitious additives and reinforcing materials [4]. In practice, however, the freezing-thawing weathering process and dynamic traffic loads could strongly influence the long-term dynamic behaviors of the modified loess substructure and could threaten the operational safety of infrastructures [5][6][7]. Therefore, how to estimate the dynamic properties of reinforced loess in these complex environments has become an urgent problem that needs to be solved in cold regions.
Loess-reinforcing techniques to create better pavements based on conventional stabilizers, such as cement and lime, have become more mature due to long-term progress and practices [8][9][10]. However, the productions of these materials are energy-intensive, result in high CO 2 emissions, and are becoming increasingly expensive. Fly ash (FA) is a waste material from coal-firing power plants, and hundreds of millions of tons of it is generated annually [11]. In the context of energy-saving and sustainable development policies, many experimental research studies have been focused on FA-based reinforcing techniques for different types of soils. Most of the existing work on FA-reinforcing techniques focus on

Materials
The studied loess was collected from the southern region of Xi'an, Shaanxi Province, NW China. The key physical and index properties are detailed in Table 1. Figure 1 exhibits the gradation curves of the loess and fly ash (FA) in the present work. The FA and polypropylene fiber (PF) were supplied by a firm in Henan Province, China. FA is a fine, powdery material that is primarily composed of spherical glassy particles, and its chemical composition is listed in Table 2. With a calcium oxide content of 7.7%, it can be classified as Class F [27]. Table 3 provides the main PF parameters, which were provided by the manufacturer.
Before mixing the samples, all of the loess and additives were pulverized, oven-dried, and then measured by mass according to the predesignated mix ratio. Afterward, the samples were thoroughly hand-mixed with corresponding initial water content. In order to reduce the fiber's electrostatic effect for a more even mixture, all of the materials were divided into 10 parts prior to mixing. Every part of the mixture was manually stirred until visual uniformity was achieved. The prepared mixtures were cast in 10 sequential layers into a standard metallic mold (50 mm diameter, 100 mm high). Every layer surface was scraped prior to subsequent compaction to avoid the delamination effect. Finally, the mixtures were compacted by targeting the maximum dry density of T01 (1.76 g/cm 3 ). After that, the new molded specimens were sealed and placed in a curing chamber for 7 days at a constant temperature of 20.0 ± 2 • C. After curing, the samples were subjected to a cyclic F-T process in a programmable freeze-thaw chamber (Figure 2a). According to the test results derived by Orakoglu et al. [13], a uniform temperature inside the sample can be achieved within 12 h of the freezing/thawing process. Considering this, the specimens were conditioned at −10 • C for 12 h and were subsequently allowed to thaw at +20 • C for 12 h, which means that one F-T cycle took 24 h. This procedure was repeated 0, 1, 2, 3, 5, 7, and 10 times since a new equilibrium could be achieved after 5 to 10 F-T cycles for most soils [28]. Before mixing the samples, all of the loess and additives were pulverized, oven-dried, and then measured by mass according to the predesignated mix ratio. Afterward, the samples were thoroughly hand-mixed with corresponding initial water content. In order to reduce the fiber's electrostatic effect for a more even mixture, all of the materials were divided into 10 parts prior to mixing. Every part of the mixture was manually stirred until visual uniformity was achieved. The prepared mixtures were cast in 10 sequential layers into a standard metallic mold (50 mm diameter, 100 mm high). Every layer surface was scraped prior to subsequent compaction to avoid the delamination effect. Finally, the mixtures were compacted by targeting the maximum dry density of T01 (1.76 g/cm 3 ). After that, the new molded specimens were sealed and placed in a curing chamber for 7 days at a constant temperature of 20.0 ± 2 °C.
After curing, the samples were subjected to a cyclic F-T process in a programmable freeze-thaw chamber (Figure 2a). According to the test results derived by Orakoglu et al. [13], a uniform temperature inside the sample can be achieved within 12 h of the freezing/thawing process. Considering this, the specimens were conditioned at −10 °C for 12 h and were subsequently allowed to thaw at +20 °C for 12 h, which means that one F-T cycle took 24 h. This procedure was repeated 0, 1, 2, 3, 5, 7, and 10 times since a new equilibrium could be achieved after 5 to 10 F-T cycles for most soils [28].

Apparatus and Testing Procedure
GDS DYNTTS manufactured by Global Digital Systems Ltd. Instruments (Figure 2b) was utilized to conduct dynamic triaxial tests in this work. The strain-dependent dynamic behaviors of reinforced loess can be obtained by applying multi-stage dynamic stress control ( Figure 3). All of the specimens were subjected to about 15 different cyclic stress amplitudes, and each uniform amplitude consisted of 50 sinusoidal load cycles.

Apparatus and Testing Procedure
GDS DYNTTS manufactured by Global Digital Systems Ltd. Instruments (Figure 2b) was utilized to conduct dynamic triaxial tests in this work. The strain-dependent dynamic behaviors of reinforced loess can be obtained by applying multi-stage dynamic stress control ( Figure 3). All of the specimens were subjected to about 15 different cyclic stress amplitudes, and each uniform amplitude consisted of 50 sinusoidal load cycles.

Determination of Dynamic Parameters
The dynamic behavior of reinforced samples can be determined by their response to multi-stage dynamic loads. From the test results, the dynamic shear stress d τ and dy-

Determination of Dynamic Parameters
The dynamic behavior of reinforced samples can be determined by their response to multi-stage dynamic loads. From the test results, the dynamic shear stress τ d and dynamic strain γ d can be determined by where σ d is the dynamic axial stress; ε d is the dynamic axial strain; and µ is Poisson's ratio. Figure 4 presents the archetype of the τ d − γ d hysteresis loops of soils under dynamic axial loading. The backbone curves are obtained by connecting the loop peaks of different stress amplitudes. In the present work, the dynamic shear modulus G d was calculated as follows: where τ d and γ d are the amplitudes of the dynamic shear stress and strain, respectively.

Determination of Dynamic Parameters
The dynamic behavior of reinforced samples can be determined by their response to multi-stage dynamic loads. From the test results, the dynamic shear stress d τ and dynamic strain d γ can be determined by where d σ is the dynamic axial stress; d ε is the dynamic axial strain; and μ is Poisson's ratio. Figure 4 presents the archetype of the d d -τ γ hysteresis loops of soils under dynamic axial loading. The backbone curves are obtained by connecting the loop peaks of different stress amplitudes. In the present work, the dynamic shear modulus d G was calculated as follows: where d τ and d γ are the amplitudes of the dynamic shear stress and strain, respectively.  The damping ratio is a key parameter that reflects the energy dissipation characteristics of soils. According to Hardin and Drnevich [29], it can be calculated from the hysteresis loop using the following equation: where W D , which represents the energy loss, equals the area enclosed by the hysteresis loop, and W A , which represents the elastic strain energy, equals the area of triangle OAB. The hyperbolic model is applied in the backbone curve fitting as follows: where a and b are two positive parameters obtained from data fitting. Accordingly, the dynamic shear modulus can be obtained from: Water 2022, 14, 317 where γ dr is the reference strain. According to Equation (6), G d reaches its maximum value at γ d = 0, which is denoted as the maximum dynamic shear modulus, G dmax = 1/a; τ d approaching infinity is denoted as ultimate dynamic shear stress, τ dult = 1/b. Accordingly, γ dr is equal to a/b or τ dult /G dmax .
Hardin and Drnevich suggested that λ is a function of the normalized shear modulus G d /G dmax . In this study, the following relationship is adopted to describe strain-dependent behavior of the damping ratio: where λ max is the maximum damping ratio, and n is a fitting constant.

Effect of Freeze-Thaw Cycle
Figure 5a exhibits the backbone curves (τ d − γ d relationships) for various F-T cycles. Generally, all of the specimens present strain-hardening characteristics, which not only illustrates how τ d increases as the strain amplitude increases, but also the applicability of the Hardin model. The backbone curve gradually shifts downward with cyclic F-T, indicating the deterioration in the dynamic performance of the reinforced loess under F-T conditions. Moreover, this downward shift recedes with the number of F-T cycles. Specifically, the backbone curve shifts downward sharply when N ≤ 5, while this phenomenon becomes less obvious as the number of F-T cycles increases. Similar observations were also reported by Orakoglu et al. [13] and Jing et al. [30] for reinforced clayey soil. The reason for these phenomena lies in the deteriorating effects of F-T on the developed cementation among particles and fibers [14,31,32]. During the freezing stage, the water between particles freezes and its volume increase by 11%. This expansion results in the development of stress, micro-pore breakage, and finally, cracks being generated within the soil structure. During the thawing period, the void space and cracks in the soil structure propagate further, leading to soil particle rearrangement to create weaker interparticle bonds. Therefore, the reinforced loess endures a considerable loss in stiffness (lower d G ) while presenting high energy dissipation characteristics (higher λ ).

Effect of Fly Ash
The influence of FA on the dynamic properties, d τ , d G , and λ , of reinforced loess is demonstrated in Figure 6. It can be observed in Figure 6a that the backbone curve shifts upward as the FA content increases. dult τ increased from 35.79 kPa to 82.55 kPa as the FA content increased from 0% to 30%. The specimen without FA exhibited the lowest dmax G value of 19.81 MPa, whereas the specimen with a FA content of 30% reached the highest dmax G -value of 32.21 MPa, an increase of about 62.59% (Figure 6b). This is presumably related to the cementation effect caused by the pozzolanic reactions of FA [15]. This cementation effect provides a confinement between soil particles and the fiber, which results in higher stiffness and a lower damping ratio.  Figure 5b,c show the variations in G d and λ with shear strain γ d , respectively. As is shown, G d decreases with the number of F-T cycles, while λ increases. Notably, the specimen that does not experience F-T presented the best G dmax of approximately 48.52 MPa, while the specimen subjected to 5 F-T cycles achieved the smallest G dmax -value of 28.15 MPa (a reduction of 41.98%). Meanwhile, λ max increased by 52.58%, from 27.54% to 42.02%.
The reason for these phenomena lies in the deteriorating effects of F-T on the developed cementation among particles and fibers [14,31,32]. During the freezing stage, the water between particles freezes and its volume increase by 11%. This expansion results in the development of stress, micro-pore breakage, and finally, cracks being generated within the soil structure. During the thawing period, the void space and cracks in the soil structure propagate further, leading to soil particle rearrangement to create weaker interparticle bonds. Therefore, the reinforced loess endures a considerable loss in stiffness (lower G d ) while presenting high energy dissipation characteristics (higher λ).

Effect of Fly Ash
The influence of FA on the dynamic properties, τ d , G d , and λ, of reinforced loess is demonstrated in Figure 6. It can be observed in Figure 6a that the backbone curve shifts upward as the FA content increases. τ dult increased from 35.79 kPa to 82.55 kPa as the FA content increased from 0% to 30%. The specimen without FA exhibited the lowest G dmax -value of 19.81 MPa, whereas the specimen with a FA content of 30% reached the highest G dmax -value of 32.21 MPa, an increase of about 62.59% (Figure 6b). This is presumably related to the cementation effect caused by the pozzolanic reactions of FA [15]. This cementation effect provides a confinement between soil particles and the fiber, which results in higher stiffness and a lower damping ratio.
water between particles freezes and its volume increase by 11%. This expansion results in the development of stress, micro-pore breakage, and finally, cracks being generated within the soil structure. During the thawing period, the void space and cracks in the soil structure propagate further, leading to soil particle rearrangement to create weaker interparticle bonds. Therefore, the reinforced loess endures a considerable loss in stiffness (lower d G ) while presenting high energy dissipation characteristics (higher λ ).

Effect of Fly Ash
The influence of FA on the dynamic properties, d τ , d G , and λ , of reinforced loess is demonstrated in Figure 6. It can be observed in Figure 6a that the backbone curve shifts upward as the FA content increases. dult τ increased from 35.79 kPa to 82.55 kPa as the FA content increased from 0% to 30%. The specimen without FA exhibited the lowest dmax G value of 19.81 MPa, whereas the specimen with a FA content of 30% reached the highest dmax G -value of 32.21 MPa, an increase of about 62.59% (Figure 6b). This is presumably related to the cementation effect caused by the pozzolanic reactions of FA [15]. This cementation effect provides a confinement between soil particles and the fiber, which results in higher stiffness and a lower damping ratio. Generally, a decrease in λ was observed when the FA content increased (Figure 6c). This is because the addition of FA improves the rigidity of loess, reducing particle slippage and rearrangement, and thus, the reinforced loess exhibits lower energy dissipation characteristics [15]. Generally, a decrease in λ was observed when the FA content increased (Figure 6c). This is because the addition of FA improves the rigidity of loess, reducing particle slippage and rearrangement, and thus, the reinforced loess exhibits lower energy dissipation characteristics [15]. Figure 7 depicts the effects of PF on the τ d , G d , and λ of reinforced loess. It can be seen that both the backbone curve and G d − γ d curve present an upward-shifting trend first and then downward-shifting as the PF content increases, which indicates an optimal PF content. Specifically, compared with 0% PF content, τ dult of specimen with 0.25%, 0.5 % PF content increase by 20.63%and 47.67%, respectively. Nevertheless, it decreased by 15.07% as the PF content further increased from 0.5% to 1.0%. On the other hand, with the increase in the PF content from 0% to 0.5%, G dmax gradually increased by 20.89%, from 23.29 MPa to 28.15 MPa, and then declined by 8.22% as the PF content continued to increase (Figure 7b). The most noteworthy effect of PF was found on the specimen with a PF content of 0.5%. As for λ, it increased with the addition of PF, especially when γ d was larger than 0.1%.

Effect of Polypropylene Fiber
Generally, the addition of PF does not achieve a significant improvement in the stiffness and energy dissipation behavior of the reinforced loess compared to the effects of FA. The effects of PF can be explained by the cementing force formed between the PF and the mixture, resulting in significant increases in the shear stress as well as in the shear modulus. The disappointing reinforcement effect of a higher PF content may be due to the electrostatic interaction between PF that leads to the fibers forming groups. As a result, unevenly distributed fibers would decrease the reinforcement effect and destroy the soil structure, degrading the dynamic performance of reinforced soil.
The effects of PF can be explained by the cementing force formed between the PF and the mixture, resulting in significant increases in the shear stress as well as in the shear modulus. The disappointing reinforcement effect of a higher PF content may be due to the electrostatic interaction between PF that leads to the fibers forming groups. As a result, unevenly distributed fibers would decrease the reinforcement effect and destroy the soil structure, degrading the dynamic performance of reinforced soil.  Figure 8a,b present the d τ and d G of reinforced loess with initial water contents of 13%, 16%, 18.9%, 22%, and 25%, respectively. As expected, both of them exhibit a decreasing trend as the initial water content increases. To be specific, the dult τ of reinforced loess slides from 68.25 kPa to 54.71 kPa, and dmax G decreases from 34.69 MPa to 23.98 MPa with water contents ranging from 13% to 25%. However, there is an increment in λ as the initial water content increases (Figure 8c).

Effect of Initial Water Content
The effect of the initial water content can be attributed to two reasons. Firstly, the friction interaction inside the specimen is weakened due to the thickened bound water and increased free water [33]. Consequently, the soil particles in the reinforced loess can move easily under dynamic loads, which contributes to the deterioration of dynamic behavior. More than that, the deterioration effect of F-T is enhanced due to the increased amount of pore water, which turns into ice and separates the soil particles during F-T process.  Figure 8a,b present the τ d and G d of reinforced loess with initial water contents of 13%, 16%, 18.9%, 22%, and 25%, respectively. As expected, both of them exhibit a decreasing trend as the initial water content increases. To be specific, the τ dult of reinforced loess slides from 68.25 kPa to 54.71 kPa, and G dmax decreases from 34.69 MPa to 23.98 MPa with water contents ranging from 13% to 25%. However, there is an increment in λ as the initial water content increases (Figure 8c).

Effect of Confining Pressure
The effects of confining pressure on the d τ , d G , and λ of reinforced loess are presented in Figure 9. As is shown, d τ and d G increase as the confining pressure increase. An increase in dult τ of up to 70.96% and an increase in dmax G of 67.77% are observed with confining pressures ranging from 12.5 kPa to 400 kPa. Meanwhile, a marked decrease in λ is observed as the confining pressure increases. For example, at a strain amplitude of 0.5%, the value of λ is reduced from 32.18% to 24.00% when the confining pressure increases from 12.5 kPa to 400 kPa.  The effect of the initial water content can be attributed to two reasons. Firstly, the friction interaction inside the specimen is weakened due to the thickened bound water and increased free water [33]. Consequently, the soil particles in the reinforced loess can move easily under dynamic loads, which contributes to the deterioration of dynamic behavior. More than that, the deterioration effect of F-T is enhanced due to the increased amount of pore water, which turns into ice and separates the soil particles during F-T process.

Effect of Confining Pressure
The effects of confining pressure on the τ d , G d , and λ of reinforced loess are presented in Figure 9. As is shown, τ d and G d increase as the confining pressure increase. An increase in τ dult of up to 70.96% and an increase in G dmax of 67.77% are observed with confining pressures ranging from 12.5 kPa to 400 kPa. Meanwhile, a marked decrease in λ is observed as the confining pressure increases. For example, at a strain amplitude of 0.5%, the value of λ is reduced from 32.18% to 24.00% when the confining pressure increases from 12.5 kPa to 400 kPa. sented in Figure 9. As is shown, d τ and d G increase as the confining pressure increase. An increase in dult τ of up to 70.96% and an increase in dmax G of 67.77% are observed with confining pressures ranging from 12.5 kPa to 400 kPa. Meanwhile, a marked decrease in λ is observed as the confining pressure increases. For example, at a strain amplitude of 0.5%, the value of λ is reduced from 32.18% to 24.00% when the confining pressure increases from 12.5 kPa to 400 kPa. It can be observed that the dynamic behavior of the reinforced loess is significantly improved by the confining pressure. This can be explained by the enhanced friction force between the particles and by the particle rearrangement in the soil, which leads to less open void space, less partial micro-cracks, and stronger inter-particle contact [34]. In addition to this, this effect also leads to less energy dissipation during wave propagation, and leads to further decreases in λ .

Theoretical Analytical Formulations
According to Equation (7), dmax G and max λ are important parameters in estimating the dynamic properties of reinforced loess. Therefore, their relationships with the five key factors, N, FA C , PF C , w, and 3 σ , was evaluated first. Subsequently, novel models for the dynamic shear modulus and the damping ratio of reinforced loess in cold regions are constructed.
Before analysis, the standard maximum dynamic shear modulus dmax0 G and standard maximum damping ratio max0 λ are introduced to normalize the data and to evaluate It can be observed that the dynamic behavior of the reinforced loess is significantly improved by the confining pressure. This can be explained by the enhanced friction force between the particles and by the particle rearrangement in the soil, which leads to less open void space, less partial micro-cracks, and stronger inter-particle contact [34]. In addition to this, this effect also leads to less energy dissipation during wave propagation, and leads to further decreases in λ.

Theoretical Analytical Formulations
According to Equation (7), G dmax and λ max are important parameters in estimating the dynamic properties of reinforced loess. Therefore, their relationships with the five key factors, N, C FA , C PF , w, and σ 3 , was evaluated first. Subsequently, novel models for the dynamic shear modulus and the damping ratio of reinforced loess in cold regions are constructed.
Before analysis, the standard maximum dynamic shear modulus G dmax0 and standard maximum damping ratio λ max0 are introduced to normalize the data and to evaluate the influences of various factors. The G dmax and λ max values of sample T01, which are equal to 28.15 MPa and 34.71%, respectively, were selected. Figure 10 exhibits the relationships of G dmax and λ max with the number of F-T cycles. It can be observed G dmax sharply decreases from 48.52 MPa to 28.15 MPa and reaches its stable state after approximately five F-T cycles. Meanwhile, λ max increases almost linearly with the number of F-T cycles. Their relationship with the number of F-T cycles can be described by the following equations:

Normalized Maximum Shear Modulus and Damping Ratio
where g N (N) and d N (N) represent the normalized parameters, G dmax /G dmax0 and λ max / λ max0 , which are affected by the number of F-T cycles.    Figure 11 depicts the variations in dmax G and max λ caused by the FA content. As shown, dmax G keeps increasing near-linearly from 20.2 MPa to 32.5 MPa, achieving an increase of 60% when the FA content is increased from 0% to 30%. Similarly, max λ exhibits a linearly increasing tendency from 44.85% to 32.11% with the FA content. The following equations are used to describe these relationships:
where ( )   Figure 11 depicts the variations in G dmax and λ max caused by the FA content. As shown, G dmax keeps increasing near-linearly from 20.2 MPa to 32.5 MPa, achieving an increase of 60% when the FA content is increased from 0% to 30%. Similarly, λ max exhibits a linearly increasing tendency from 44.85% to 32.11% with the FA content. The following equations are used to describe these relationships: where g f (C FA ) and d f (C FA ) represent the normalized parameters, G dmax /G dmax0 and λ max /λ max0 , which are affected by the FA content.
Water 2022, 14, x FOR PEER REVIEW 11 of 17 Figure 11. Effects of fly ash content on dmax G and max λ . Figure 12 shows the effects of PF on the dmax G and max λ of reinforced loess. Generally, the effects of the PF content were found to not be significant in comparison to other factors. When the PF content increased from 0.0% to 1.0%, dmax G first experienced a moderate increase, and thereafter, it decreased slightly, while max λ increased almost linearly with the PF content. A quadratic polynomial and a linear model were proposed, respectively, for dmax G and max λ as follows: ( 10 00 where ( )   Figure 12 shows the effects of PF on the G dmax and λ max of reinforced loess. Generally, the effects of the PF content were found to not be significant in comparison to other factors. When the PF content increased from 0.0% to 1.0%, G dmax first experienced a moderate increase, and thereafter, it decreased slightly, while λ max increased almost linearly with the PF content. A quadratic polynomial and a linear model were proposed, respectively, for G dmax and λ max as follows: where g p (C PF ) and d p (C PF ) represent the normalized parameters, G dmax /G dmax0 and λ max /λ max0 , which are affected by the PF content. 10 00 where ( )   Figure 13 presents the effects of the initial water content on the dmax G and max λ of reinforced loess. As shown in the figure, a noticeable decrease in dmax G and an increase in max λ are observed when the initial water content increases from 13% to 25%. The effects of the initial water content can be described by the following equations: Figure 12. Effects of polypropylene fiber content on G dmax and λ max . Figure 13 presents the effects of the initial water content on the G dmax and λ max of reinforced loess. As shown in the figure, a noticeable decrease in G dmax and an increase in λ max are observed when the initial water content increases from 13% to 25%. The effects of the initial water content can be described by the following equations: where g w (w) and d w (w) represent the normalized parameters, G dmax /G dmax0 and λ max /λ max0 , which are affected by the initial water content.   28 87    Figure 14 generates G dmax and λ max versus the confining pressure on the semilogarithmic scale. As shown in the figure, an approximately linear increase in G dmax and an approximately linear decrease in λ max are observed as the confining pressure increases. Through regression analysis, G dmax and λ max are described by the following Equations: where g σ (σ 3 ) and d σ (σ 3 ) represent the normalized parameters, G dmax /G dmax0 and λ max / λ max0 , which are affected by the confining pressure; p a is the atmospheric pressure, 101.3 kPa.

Empirical Expression for Dynamic Shear Modulus
Considering the effects of F-T cycles, FA content, PF content, initial water content, and confining pressure, Equations (19) and (20) are proposed to predict the G dmax and λ max of reinforced loess in Equations (7) and (8): where the parameters and normalized functions of G dmax and λ max have been discussed in the former parts. Figure 15 presents the comparisons of the experimental G dmax and λ max with the predicted values from Equations (19) and (20). It can be observed that both the predicted G dmax and λ max exhibit an approximate linear relationship of 1:1 with the test data, which indicates that the G dmax and λ max of reinforced loess can be estimated from these equations well. Considering the effects of F-T cycles, FA content, PF content, initial water content, and confining pressure, Equations (19) and (20) are proposed to predict the dmax G and max λ of reinforced loess in Equations (7) and (8): where the parameters and normalized functions of dmax G and max λ have been discussed in the former parts. Figure 15 presents the comparisons of the experimental dmax G and max λ with the predicted values from Equations (19) and (20). It can be observed that both the predicted dmax G and max λ exhibit an approximate linear relationship of 1:1 with the test data, which indicates that the dmax G and max λ of reinforced loess can be estimated from these equations well. Generally, d G and λ are normalized by dmax G and max λ , respectively. As shown in Figure 16, the normalized dynamic shear modulus at different experimental conditions presents a similar variation tendency, and so does the normalized damping ratio. Equations (7) and (8) are utilized to fit these normalized data and result in Equations (21) and (22): where dr γ and n are normalized the parameters and are equal to 0.2132 and 0.6757, Generally, G d and λ are normalized by G dmax and λ max , respectively. As shown in Figure 16, the normalized dynamic shear modulus at different experimental conditions presents a similar variation tendency, and so does the normalized damping ratio. Equations (7) and (8) are utilized to fit these normalized data and result in Equations (21) and (22): where γ dr and n are normalized the parameters and are equal to 0.2132 and 0.6757, respectively.
Similarly, an empirical expression for the damping ratio can be deduced as: where d G and dmax0 G are expressed in MPa; FA C , PF C , and w are expressed as percentages (%); 3 σ and a p are expressed in kPa; and dmax0 G and max0 λ are equal to 28.15 MPa and 34.71%, respectively. Figure 17 presents the comparisons of the experimental values of d G and λ with the corresponding predicted values by Equations (23) and (24). It can be seen that the predicted results are in good agreement with the experimental values. The correlation coefficients are up to 0.952 and 0.973, respectively. It can be concluded that the prediction models can be confidently applied to estimate d G and λ in designing FA-PF-reinforced loess. Substituting Equation (19) into Equation (21), a comprehensive function for the G d of modified loess by considering multiple factors is established as follows:  (24) where G d and G dmax0 are expressed in MPa; C FA , C PF , and w are expressed as percentages (%); σ 3 and p a are expressed in kPa; and G dmax0 and λ max0 are equal to 28.15 MPa and 34.71%, respectively. Figure 17 presents the comparisons of the experimental values of G d and λ with the corresponding predicted values by Equations (23) and (24). It can be seen that the predicted results are in good agreement with the experimental values. The correlation coefficients are up to 0.952 and 0.973, respectively. It can be concluded that the prediction models can be confidently applied to estimate G d and λ in designing FA-PF-reinforced loess.

Conclusions
In this study, a series of dynamic triaxial tests were performed on loess that had been reinforced with FA and PF. The dynamic shear stress, dynamic shear modulus, and damping ratio were analyzed by considering five key factors for reinforced soil in cold regions. Finally, the impacts of these factors were incorporated into newly established models to determine the dynamic shear modulus and damping ratio. The following major conclusions can be drawn: (1) The reinforced samples exhibited strain-hardening behavior. The Hardin model was adopted to describe the dynamic shear stress-strain relationships and exhibited good agreement with the experimental data. (2) A considerable reduction in the dynamic shear modulus was observed with increasing freeze-thaw cycles and initial water content. On the contrary, the dynamic shear modulus increased as the FA content and confining pressure increases. The sample with 0.5% PP fiber showed the best dynamic performance, and generally, the impact of PP fiber is not significant by comparison. (3) The damping ratio decreased as the FA content and confining pressure increased.
However, it increased with cyclic freezing-thawing, PF content, and initial water content, which illustrates the better energy dissipation performance of the reinforced loess. (4) Novel empirical models for the dynamic shear modulus and damping ratio were established by taking account of freeze-thaw cycle, FA content, PF content, initial water content, and confining pressure. The comparison of the experimental and predicted results illustrated that the newly established models are suitable to estimate the dynamic shear modulus and damping ratio of fly ash-polypropylene fiber-reinforced loess.

Conclusions
In this study, a series of dynamic triaxial tests were performed on loess that had been reinforced with FA and PF. The dynamic shear stress, dynamic shear modulus, and damping ratio were analyzed by considering five key factors for reinforced soil in cold regions. Finally, the impacts of these factors were incorporated into newly established models to determine the dynamic shear modulus and damping ratio. The following major conclusions can be drawn: (1) The reinforced samples exhibited strain-hardening behavior. The Hardin model was adopted to describe the dynamic shear stress-strain relationships and exhibited good agreement with the experimental data. (2) A considerable reduction in the dynamic shear modulus was observed with increasing freeze-thaw cycles and initial water content. On the contrary, the dynamic shear modulus increased as the FA content and confining pressure increases. The sample with 0.5% PP fiber showed the best dynamic performance, and generally, the impact of PP fiber is not significant by comparison. (3) The damping ratio decreased as the FA content and confining pressure increased. However, it increased with cyclic freezing-thawing, PF content, and initial water content, which illustrates the better energy dissipation performance of the reinforced loess. (4) Novel empirical models for the dynamic shear modulus and damping ratio were established by taking account of freeze-thaw cycle, FA content, PF content, initial water content, and confining pressure. The comparison of the experimental and predicted results illustrated that the newly established models are suitable to estimate the dynamic shear modulus and damping ratio of fly ash-polypropylene fiber-reinforced loess.  Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.