The Elastic Modulus and Damage Stress–Strain Model of Polypropylene Fiber and Nano Clay Modified Lime Treated Soil under Axial Load

Using polypropylene fiber (PPF) and nano clay modified lime treated soil (LS), the static and dynamic properties of fiber modified lime treated soil (FLS), nano clay modified lime treated soil (NLS), and fiber nano clay composite modified lime treated soil (NFLS) were studied. Through the unconfined compressive strength (UCS) test and dynamic triaxial test of FLS, NLS, and NFLS, the static and dynamic elastic modulus characteristics at 7 day curing age were explored, and the damage stress–strain model was established. The results show that: (1) Polypropylene fiber and nano clay can significantly enhance the mechanical properties of NFLS. Nano clay can promote the reaction between lime and soil to produce calcium silicate hydrate (C-S-H) and calcium aluminate hydrate (C-A-H), thus improving the strength of NFLS, and UCS can be increased by up to 103%. Polypropylene fiber can enhance the ductility of NFLS and increase the residual ductility strength, and the residual strength can be increased by 827%. (2) Nano clay can enhance the static and dynamic elastic modulus of modified lime treated soil. The static and dynamic elastic modulus of NLS, FLS, and NFLS are linear with the change of polypropylene fiber and nano clay content. The static and dynamic elastic modulus of NLS, FLS, and NFLS are linear, exponential, and logarithmic, respectively. (3) The mesoscopic random damage model can characterize the stress–strain relationship of NFLS. Polypropylene fiber and nano clay can improve the ductility and strength of modified LS, and the composite addition of polypropylene fiber and nano clay can improve the ability of modified LS to resist damage.


Introduction
Due to the demand of economic and urbanization development, the quality requirements of road infrastructure construction are increasing year by year. Using lime treated construction waste soil for road base is an effective measure to reduce construction waste. At the same time, it can replace sand, stone, and other materials, so as to reduce the irreversible damage to the ecological environment caused by mountain excavation and river excavation. The traffic load is mainly axial load. Due to the increase of road load, higher requirements are put forward for the mechanical properties of road base. Therefore, it is necessary to deeply explore the properties of lime treated soil and the main technical means to improve its mechanical properties.
Liu et al. [1] proved that after adding lime into the soil, the soil properties change, which can improve its unconfined compressive strength (UCS) to meet the actual engineering needs. Lime can not only significantly improve the UCS of the improved soil and reduce the expansion and contraction of the soil, but also react with the pozzolanic ash to form cementitious materials, which is an excellent material for stabilizing the soil [2,3]. Lemaire et al. [4] conducted UCS and X-ray diffraction (XRD) tests on lime treated soil. The results show that lime can make the soil form a "honeycomb" structure and significantly 88% and purity of 98%. The polypropylene fiber was a bundle monofilament, which was produced by Shaoxing Fiber High Tech Co., Ltd., Shaoxing City, Zhejiang Province, China, with a diameter of 18-46 µm and length of 6 mm, which has the characteristics of light weight and good stability, and can significantly improve the ductility and bending resistance of soil. The nano clay used was produced by Hubei Jinxi montmorillonite Co., Ltd., Wuhan City, Hubei Province, China, the density is 0.45 g/cm 3 , 99.9% of nano clay can pass the 200 mesh screen, and the content of montmorillonite is 96-98%, which is beige powder.

Test Scheme
LS was modified by nano clay, polypropylene fiber, composite nano clay and polypropylene fiber. The static and dynamic properties of the modified LS were studied by UCS test and dynamic triaxial test. The content of lime was 6%, and, on this basis, nano clay with the content of 0%, 2%, 4%, and 6% and polypropylene fiber with the content of 0%, 0.25%, 0.5%, and 0.75% was added respectively. The content of lime, nano clay, and polypropylene fiber was the mass ratio of dry soil. The mechanical properties of NFLS were investigated by adding 2%, 4%, and 6% nano clay in the optimal FLS and 0.25%, 0.5%, and 0.75% polypropylene fiber in the optimal NLS [20].
According to the Chinese test code for Inorganic Binder Stabilized Materials in Highway Engineering [22], the undisturbed soil of the required soil sample is crushed and dried in a 105 • C oven. After crushing, a 2 mm sieve is used to screen the soil. According to the test scheme, after weighing the corresponding quality of dry soil, quicklime, nano clay, and polypropylene fiber with an electronic scale, add water and stir evenly according to the optimal moisture content to make each admixture disperse and not agglomerate. Put the uniformly stirred aggregate into the sealing stuffs of plastic bags for 24 h, and make samples by static compaction, as shown in Figure 1. After demoulding, it should be cured for 7 days. In the first 6 days, it should be cured in a standard curing box with a temperature of 20 ± 2 • C and a relative humidity of more than 95%. On the last day, the sample should be taken out of the curing box and soaked in water at 20 ± 2 • C for curing. The water surface is approximately 2.5 cm above the sample. The test scheme is shown in Table 2. Note, L, N, and F respectively represent lime, nano clay, and polypropylene fiber. Lime content, nano clay content, and polypropylene content are the percentage of dry soil mass, and water content is the percentage of all dry mixtures. The UCS test sample is a cylinder with a diameter of 50 mm and a height of 50 mm. After curing, TKA-WCY-1F automatic unconfined compressive strength tester is selected for testing, as shown in Figure 2. Referring to the Chinese code for Highway Geotechnical Test [23], the loading rate of the instrument was 1 mm/min, and the stress-strain curve of the sample is obtained through the computer data acquisition system.
Because the traffic load has a long-term load impact on the road base, the dynamic elastic modulus of NFLS is studied by dynamic triaxial test to deeply explore the longterm traffic load bearing capacity of modified LS. After curing, GDS servo motor control system is selected for dynamic triaxial test. The size of the dynamic triaxial test sample is a cylinder with a diameter of 39.1 mm and a height of 80 mm. Because the surrounding earth pressure has little impact on the road base, the confining pressure is 0 kPa, and the displacement control method is used to vibrate for 50 times. The test is carried out with a sine wave with a frequency of 1 Hz and an amplitude of 0.1 mm. Data acquisition is carried out through microcomputer system, as shown in Figure 3.  The UCS test sample is a cylinder with a diameter of 50 mm and a height of 50 mm. After curing, TKA-WCY-1F automatic unconfined compressive strength tester is selected for testing, as shown in Figure 2. Referring to the Chinese code for Highway Geotechnical Test [23], the loading rate of the instrument was 1 mm/min, and the stress-strain curve of the sample is obtained through the computer data acquisition system. The UCS test sample is a cylinder with a diameter of 50 mm and a height of 50 mm. After curing, TKA-WCY-1F automatic unconfined compressive strength tester is selected for testing, as shown in Figure 2. Referring to the Chinese code for Highway Geotechnical Test [23], the loading rate of the instrument was 1 mm/min, and the stress-strain curve of the sample is obtained through the computer data acquisition system.
Because the traffic load has a long-term load impact on the road base, the dynamic elastic modulus of NFLS is studied by dynamic triaxial test to deeply explore the longterm traffic load bearing capacity of modified LS. After curing, GDS servo motor control system is selected for dynamic triaxial test. The size of the dynamic triaxial test sample is a cylinder with a diameter of 39.1 mm and a height of 80 mm. Because the surrounding earth pressure has little impact on the road base, the confining pressure is 0 kPa, and the displacement control method is used to vibrate for 50 times. The test is carried out with a sine wave with a frequency of 1 Hz and an amplitude of 0.1 mm. Data acquisition is carried out through microcomputer system, as shown in Figure 3.  Because the traffic load has a long-term load impact on the road base, the dynamic elastic modulus of NFLS is studied by dynamic triaxial test to deeply explore the long-term traffic load bearing capacity of modified LS. After curing, GDS servo motor control system is selected for dynamic triaxial test. The size of the dynamic triaxial test sample is a cylinder with a diameter of 39.1 mm and a height of 80 mm. Because the surrounding earth pressure has little impact on the road base, the confining pressure is 0 kPa, and the displacement control method is used to vibrate for 50 times. The test is carried out with a sine wave with a frequency of 1 Hz and an amplitude of 0.1 mm. Data acquisition is carried out through microcomputer system, as shown in Figure 3.  For SEM test, a Nissan high-low vacuum scanning electron microscope (JSM-6360LV, JEOL, TKY, Japan) was used. The electro-adhesive was evenly spread on the aluminum tray, and then the powder of the sample to be tested was glued to the electro-adhesive. After the required test sample was glued, a layer of platinum film was uniformly sprayed on the surface of the sample. After confirming uniform spraying, the aluminum tray was put into the electron microscope for SEM test.
FT-IR was tested by KBr compression method. The instrument used in the experiment was the Fourier Infrared spectrometer (NEXUS, Thermo Nicolet Corporation, American), with a wavelength range of 4000-400 cm −1 . The test sample was ground into powder and then added into potassium bromide powder. Grinding continued so that the sample powder and potassium bromide powder was fully mixed and ground into a fine powder. Then the powder was pressed into sheets by KBr method and put into the instrument for infrared spectrum analysis.

Stress-Strain Damage Model
According to the random damage mechanics of concrete, Li et al. [24,25] abstracted the concrete sample as a micro spring random damage system from the meso point of view, so as to establish the random damage constitutive model of concrete. According to the meso random damage theory, the damage stress-strain relationship of materials under compression is shown in Equation (1).
where E is the elastic modulus, D is the damage variable, and the mean and variance of D are shown in Equations (2) and (3), respectively.
The meso random damage model is mainly related to parameters λ, ζ, and ω [26,27]. E can be obtained from the dynamic triaxial test of NFLS, and then the parameters of λ, ζ, and ω can be identified according to the UCS test results, so as to establish the stressstrain damage model of NFLS. For SEM test, a Nissan high-low vacuum scanning electron microscope (JSM-6360LV, JEOL, Tokyo, Japan) was used. The electro-adhesive was evenly spread on the aluminum tray, and then the powder of the sample to be tested was glued to the electro-adhesive. After the required test sample was glued, a layer of platinum film was uniformly sprayed on the surface of the sample. After confirming uniform spraying, the aluminum tray was put into the electron microscope for SEM test.
FT-IR was tested by KBr compression method. The instrument used in the experiment was the Fourier Infrared spectrometer (NEXUS, Thermo Nicolet Corporation, American, Waltham, MA, USA), with a wavelength range of 4000-400 cm −1 . The test sample was ground into powder and then added into potassium bromide powder. Grinding continued so that the sample powder and potassium bromide powder was fully mixed and ground into a fine powder. Then the powder was pressed into sheets by KBr method and put into the instrument for infrared spectrum analysis.

Stress-Strain Damage Model
According to the random damage mechanics of concrete, Li et al. [24,25] abstracted the concrete sample as a micro spring random damage system from the meso point of view, so as to establish the random damage constitutive model of concrete. According to the meso random damage theory, the damage stress-strain relationship of materials under compression is shown in Equation (1).
where E is the elastic modulus, D is the damage variable, and the mean and variance of D are shown in Equations (2) and (3), respectively.
The meso random damage model is mainly related to parameters λ, ζ, and ω [26,27]. E can be obtained from the dynamic triaxial test of NFLS, and then the parameters of λ, ζ, and ω can be identified according to the UCS test results, so as to establish the stress-strain damage model of NFLS.

UCS Test Results
UCS and residual strength are important mechanical indexes for the study of NFLS. The peak value of stress-strain curve is taken as the UCS, and the stress corresponding to Polymers 2022, 14, 2606 6 of 18 peak strain plus 5% is selected as the residual strength. The UCS test results of NFLS are shown in Table 3.  Figure 4a,b show the stress-strain curve of NFLS with nano clay content N = 6%. The stress-strain curve of NFLS is a softening curve, that is, the stress first increases to the stress peak with strain, then decreases, and then tends to be gentle. Figure 4c,d show the UCS and residual strength of NFLS with nano clay content N = 6%. The increase of nano clay and fiber content can significantly improve the UCS of LS. The content of nano clay increases from 0% to 6%. The UCS of NFLS increases with the increase of nano clay content. When the content of nano clay is 6%, the UCS of NLS is 1947 kPa and the relative growth rate of UCS is 57%. When the optimal nano clay content is 6%, the UCS of NFLS increases with the increase of polypropylene fiber content, and when the fiber content is 0.75%, the UCS of NFLS is 2513 kPa, the growth rate relative to NLS is 29%, and the growth rate of UCS is the fastest. It shows that the strength of NFLS can be significantly improved by adding nano clay and fiber, and the optimal contents of nano clay and polypropylene fiber are 6% and 0.75% respectively. The residual strength of NLS increases first and then decreases with the increase of nano clay content, which is due to the increase of brittleness of NLS and the decrease of residual strength. When the content of nano clay is 6%, the residual strength of NFLS increases linearly with the increase of the content of polypropylene fiber. This is because the fiber can slow down the development of failure cracks of NFLS. With the increase of the content of fiber, the cracks develop more slowly and the bearing capacity is greater, indicating that polypropylene fiber can significantly enhance the ductility of NFLS and slow down its brittle failure. The interfacial force generated between the polypropylene fiber and soil can play a better role in the bonding of cracks. As shown in Figure 4e, under the same strain, with the increase of polypropylene fiber content, the fibers emerging from the crack of the sample increase significantly, the pull effect is significant, and the crack development slows down [28,29].  Figure 5a,b shows the stress-strain curve of FLS and NFLS with fiber content F = 0.75%, and Figure 5c,d shows the UCS and residual strength diagram of NFLS with fiber content F = 0.75%. With the increase of polypropylene fiber content, the UCS of FLS first increases and then decreases. When the polypropylene fiber content is 0.75%, the UCS of FLS shows a downward trend, which is due to the overhead phenomenon when the network structure is formed between the soil mass due to the excessive fiber content, and the relative compactness of the soil mass decreases, leading to a decrease in UCS growth. When the fiber content is 0.75%, the UCS of NFLS increases with the increase of nano clay content, and when the nano clay content is 6%, the UCS of NFLS increases the fastest. This is because nano clay can fill the pores between LS particles and make its internal structure more dense. Secondly, nano clay can promote the further reaction of LS, and on the other hand, it can have pozzolanic reaction with LS [30][31][32][33]. The volcanic ash reaction generates a sheet structure, as shown in Figure 6a. The sheet structure is cementitious substances such as hydrated calcium silicate C-S-H and hydrated calcium aluminate C-A-H, as  When the polypropylene fiber content is 0.75%, the UCS of FLS shows a downward trend, which is due to the overhead phenomenon when the network structure is formed between the soil mass due to the excessive fiber content, and the relative compactness of the soil mass decreases, leading to a decrease in UCS growth. When the fiber content is 0.75%, the UCS of NFLS increases with the increase of nano clay content, and when the nano clay content is 6%, the UCS of NFLS increases the fastest. This is because nano clay can fill the pores between LS particles and make its internal structure more dense. Secondly, nano clay can promote the further reaction of LS, and on the other hand, it can have pozzolanic reaction with LS [30][31][32][33]. The volcanic ash reaction generates a sheet structure, as shown in Figure 6a. The sheet structure is cementitious substances such as hydrated calcium silicate C-S-H and hydrated calcium aluminate C-A-H, as shown in Figure 6b,c. The generated cementitious material can fill the pores between the LS skeleton and make the LS particles agglomerate better. At the same time, the interface between fiber and soil particles is closer, the interfacial friction is greater, the damage degree of external force to the sample is reduced, and the strength of NFLS is further increased. However, the residual strength of NFLS increases slowly with the increase of nano clay content, indicating that polypropylene fiber contributes more to the residual strength of NFLS than nano clay. shown in Figure 6b,c. The generated cementitious material can fill the pores between the LS skeleton and make the LS particles agglomerate better. At the same time, the interface between fiber and soil particles is closer, the interfacial friction is greater, the damage degree of external force to the sample is reduced, and the strength of NFLS is further increased. However, the residual strength of NFLS increases slowly with the increase of nano clay content, indicating that polypropylene fiber contributes more to the residual strength of NFLS than nano clay. The above phenomena show that nano clay mainly provides strength for NFLS, but will enhance the brittle failure of LS. Polypropylene fiber can enhance the ductility of NFLS and improve its residual strength. shown in Figure 6b,c. The generated cementitious material can fill the pores between the LS skeleton and make the LS particles agglomerate better. At the same time, the interface between fiber and soil particles is closer, the interfacial friction is greater, the damage degree of external force to the sample is reduced, and the strength of NFLS is further increased. However, the residual strength of NFLS increases slowly with the increase of nano clay content, indicating that polypropylene fiber contributes more to the residual strength of NFLS than nano clay. The above phenomena show that nano clay mainly provides strength for NFLS, but will enhance the brittle failure of LS. Polypropylene fiber can enhance the ductility of NFLS and improve its residual strength. The above phenomena show that nano clay mainly provides strength for NFLS, but will enhance the brittle failure of LS. Polypropylene fiber can enhance the ductility of NFLS and improve its residual strength.

Variation Law of Static Elastic Modulus
The static elastic modulus is the ratio of stress to strain in the elastic deformation stage under compression. Through the unconfined compressive strength test data, the stress-strain curve of NFLS in the elastic stage under compression is fitted, and the straightline slope is the static elastic modulus. Figure 7a-d shows the variation curve of static elastic modulus under different content of nano clay and polypropylene fiber. Figure 7a shows the variation law of the static elastic modulus of NLS with the increase of nano clay content. When the content of nano clay is 6%, the static elastic modulus of NLS increases to 122 MPa, and the growth rate is 51% higher than that of LS, indicating that the addition of nano clay can improve the static elastic modulus of LS and enhance its deformation resistance. Figure 7b shows that when the content of nano clay is 6%, the static elastic modulus of NFLS decreases with the increase of the content of polypropylene fiber, because the stress of NFLS in the elastic stage decreases with the growth rate of strain. In Figure 7c, it can be found that the content of polypropylene fiber increases from 0.25% to 0.75%, and the static elastic modulus of FLS first increases and then decreases. When the content of polypropylene fiber is 0.5%, the elastic modulus of FLS is 69 MPa, while when the content of polypropylene increases by 0.75%, the static elastic modulus of FLS decreases to 55 MPa, indicating that on the one hand, this is due to the large amount of C-S-H and C-A-H produced by lime soil reaction attached to the surface of polypropylene fiber, so polypropylene fiber can enhance the adhesion of the interface between soil and fiber and make the samples bond together better, as shown in Figure 8. On the other hand, with the increase of the content of polypropylene fiber, the porosity between soil particles and fiber may increase, resulting in the decrease of its static elastic modulus. Figure 7d shows that when the content of polypropylene fiber is 0.75%, the static elastic modulus of NFLS first increases, then decreases, and finally significantly increases with the increase of nano clay content. When the content of nano clay is 6%, the static elastic modulus of NFLS is significantly increased to 90 MPa. This is because nano clay can fill the pores between fiber and soil, enhance the stiffness of NFLS, and reduce elastic deformation.
The relationship is fitted between the static elastic modulus and the content of nano clay and polypropylene fiber of NLS, FLS, and NFLS. The fitting formulas are shown in Equations (4)- (7). E (NLS) = 6.877y + 80.9, R 2 = 0.92 where E is the static elastic modulus (MPa); x is the content of polypropylene fiber (%); and y is the content of nano clay (%).

Hysteretic Curve
The dynamic elastic modulus of the material can be calculated from the slope of the hysteretic curve to characterize the stiffness of the soil [34]. Figure 9a-d shows the variation law of dynamic stress-strain hysteretic curves of NLS, FLS, and NFLS. Figure 9a,b shows that the hysteresis curves of NLS gradually approaches the stress axis with the increase of nano clay content. When the content of nano clay is 6%, different content of polypropylene fiber has a significant effect on the hysteretic curve. With the increase of polypropylene fiber content, the angle between hysteretic curve and stress axis increases. The dynamic elastic modulus of NFLS decreases with the increase of polypropylene fiber content. Figure 9c,d shows that the hysteresis curves of FLS fluctuates with the increase of polypropylene fiber content and its area and angle with the stress axis. The FLS area is the largest and the included angle with the stress axis is the smallest when the content of polypropylene fiber is 0.5%. When the content of polypropylene fiber is 0.75%, the inclination angles of hysteresis curves increase gradually with the content of nano clay.

Hysteretic Curve
The dynamic elastic modulus of the material can be calculated from the slope of the hysteretic curve to characterize the stiffness of the soil [34]. Figure 9a-d shows the variation law of dynamic stress-strain hysteretic curves of NLS, FLS, and NFLS. Figure 9a,b shows that the hysteresis curves of NLS gradually approaches the stress axis with the increase of nano clay content. When the content of nano clay is 6%, different content of polypropylene fiber has a significant effect on the hysteretic curve. With the

Hysteretic Curve
The dynamic elastic modulus of the material can be calculated from the slope of the hysteretic curve to characterize the stiffness of the soil [34]. Figure 9a-d shows the variation law of dynamic stress-strain hysteretic curves of NLS, FLS, and NFLS. Figure 9a,b shows that the hysteresis curves of NLS gradually approaches the stress axis with the increase of nano clay content. When the content of nano clay is 6%, different content of polypropylene fiber has a significant effect on the hysteretic curve. With the increase of polypropylene fiber content, the angle between hysteretic curve and stress axis increases. The dynamic elastic modulus of NFLS decreases with the increase of polypropylene fiber content. Figure 9c,d shows that the hysteresis curves of FLS fluctuates with the increase of polypropylene fiber content and its area and angle with the stress axis. The FLS area is the largest and the included angle with the stress axis is the smallest when the content of polypropylene fiber is 0.5%. When the content of polypropylene fiber is 0.75%, the inclination angles of hysteresis curves increase gradually with the content of nano clay.

Dynamic Elastic Modulus
According to the hysteretic curves, the dynamic elastic modulus of NFLS with different nano clay and polypropylene fiber content can be obtained according to Equation (8), as shown in Figure 10.  Figure 10a,b shows that the dynamic elastic modulus of NLS increases with the increase of nano clay content. The nano clay content increases from 0% to 6%, the dynamic elastic modulus of NLS increases to 189 MPa, and nano clay can significantly enhance the elastic modulus of LS. When the content of nano clay is 6%, the dynamic elastic modulus of NFLS decreases with the increase of the content of polypropylene fiber. Nano clay can enhance the dynamic elastic modulus of LS because nano clay can produce pozzolanic reaction by modifying LS to produce more cementitious materials and make its structure more compact. When the sample is pressed and does not break, the strain produced is very small, so the force on the fiber can be decomposed into a compressive force and a flexural force, which has almost no effect on the sample and mainly plays a role of reinforcement [35][36][37].

Dynamic Elastic Modulus
According to the hysteretic curves, the dynamic elastic modulus of NFLS with different nano clay and polypropylene fiber content can be obtained according to Equation (8), as shown in Figure 10. Figure 10a,b shows that the dynamic elastic modulus of NLS increases with the increase of nano clay content. The nano clay content increases from 0% to 6%, the dynamic elastic modulus of NLS increases to 189 MPa, and nano clay can significantly enhance the elastic modulus of LS. When the content of nano clay is 6%, the dynamic elastic modulus of NFLS decreases with the increase of the content of polypropylene fiber. Nano clay can enhance the dynamic elastic modulus of LS because nano clay can produce pozzolanic reaction by modifying LS to produce more cementitious materials and make its structure more compact. When the sample is pressed and does not break, the strain produced is very small, so the force on the fiber can be decomposed into a compressive force and a flexural force, which has almost no effect on the sample and mainly plays a role of reinforcement [35][36][37]. Figure 10c,d shows that the dynamic elastic modulus of FLS first decreases, then increases and then decreases with the increase of polypropylene fiber content. When the polypropylene fiber content is 0.5%, the dynamic elastic modulus of FLS is 120 MPa, while when the polypropylene fiber content increases to 0.75%, the dynamic elastic modulus decreases to 112 MPa. When the polypropylene fiber content is large, it is easy to agglomerate and interweave, reducing the compactness of soil. The capillary action of fiber will lead to soil erosion and reduce the dynamic elastic modulus of LS. When the content of polypropylene fiber is certain, the dynamic elastic modulus of NFLS decreases first and then increases significantly with the increase in the content of nano clay. When the content of nano clay is 6%, the dynamic elastic modulus of NFLS is 150 MPa.
where Ed is dynamic elastic modulus (MPa).  The relationship between dynamic elastic modulus of NLS, FLS, and NFLS and content of nano clay and polypropylene fiber was fitted. The fitting formulas are shown in Equations (9)-(12) respectively.
where E d is dynamic elastic modulus (MPa).

Relationship between Static and Dynamic Elastic Modulus
Comparing Figures 9 and 10, it can be found that the content of polypropylene fiber and nano clay has a significant impact on the static and dynamic elastic modulus of NFLS. Therefore, the prediction model of the impact of different content of polypropylene fiber and nano clay on the static and dynamic elastic modulus of modified LS can be established. Figure 11a,b shows the prediction model of the relationship between the content of polypropylene fiber, nano clay, and static and dynamic elastic modulus. It reflects the relationship between the static and dynamic elastic modulus of NFLS and the content of polypropylene fiber and nano clay. Under the same nano clay content, the predicted model surface shows a downward convex shape with the increase of polypropylene fiber content, indicating that the static and dynamic elastic modulus of NFLS first decreases and then increases with the change of polypropylene fiber content. Under the same polypropylene fiber content, the prediction model surface increases with the increase of nano clay content, indicating that the static and dynamic elastic modulus of NFLS increases linearly with the increase of nano clay content. The fitting formulas of static and dynamic elastic modulus of NFLS are shown in Equations (13) and (14): Therefore, the prediction model of the impact of different content of polypropylene fiber and nano clay on the static and dynamic elastic modulus of modified LS can be established. Figure 11a,b shows the prediction model of the relationship between the content of polypropylene fiber, nano clay, and static and dynamic elastic modulus. It reflects the relationship between the static and dynamic elastic modulus of NFLS and the content of polypropylene fiber and nano clay. Under the same nano clay content, the predicted model surface shows a downward convex shape with the increase of polypropylene fiber content, indicating that the static and dynamic elastic modulus of NFLS first decreases and then increases with the change of polypropylene fiber content. Under the same polypropylene fiber content, the prediction model surface increases with the increase of nano clay content, indicating that the static and dynamic elastic modulus of NFLS increases linearly with the increase of nano clay content. The fitting formulas of static and dynamic elastic modulus of NFLS are shown in Equations (13) and (14): From the comparison of static and dynamic elastic modulus prediction models according to Figure 11, it can be found that the change trend between the static and dynamic elastic modulus of NFLS and the content of nano clay and polypropylene fiber is basically the same. In the actual road engineering construction, the static elastic modulus is easy to obtain, but in the complex engineering environment, it is difficult to obtain the dynamic elastic modulus through the precise dynamic triaxial test. The static and dynamic elastic modulus test data of NFLS are fitted to establish the variation relationship of them, as shown in Figure 12. The static and dynamic elastic modulus of NLS, FLS, and NFLS conforms to the linear, exponential, and logarithmic relations respectively, and the fitting formula is shown in Equations (15)- (17). From the comparison of static and dynamic elastic modulus prediction models according to Figure 11, it can be found that the change trend between the static and dynamic elastic modulus of NFLS and the content of nano clay and polypropylene fiber is basically the same. In the actual road engineering construction, the static elastic modulus is easy to obtain, but in the complex engineering environment, it is difficult to obtain the dynamic elastic modulus through the precise dynamic triaxial test. The static and dynamic elastic modulus test data of NFLS are fitted to establish the variation relationship of them, as shown in Figure 12. The static and dynamic elastic modulus of NLS, FLS, and NFLS conforms to the linear, exponential, and logarithmic relations respectively, and the fitting formula is shown in Equations (15)- (17).

Stress-Strain Damage Model and Damage Evolution Law of NFLS
According to the UCS test data, the elastic modulus of NFLS is substituted into the objective function [25,26], and the parameters λ and ζ are identified by particle swarm optimization algorithm. The identification results are shown in Table 4. The mean value of damage variables µ D is calculated according to the identification results, and the stressstrain random damage model is established by using λ and ζ.  Figure 12. Relationship between static and dynamic elastic modulus.

Stress-Strain Damage Model and Damage Evolution Law of NFLS
According to the UCS test data, the elastic modulus of NFLS is substituted into the objective function [25,26], and the parameters λ and ζ are identified by particle swarm optimization algorithm. The identification results are shown in Table 4. The mean value of damage variables μD is calculated according to the identification results, and the stressstrain random damage model is established by using λ and ζ.
where σi and σ (εi) is the measured stress and model calculated stress with the strain is εi, respectively, and m is the number of measured data groups. Figure 13a-d shows the damage evolution curves of NFLS. When D = 0, the sample is not damaged. When D = 1, the sample is completely destroyed. When the strain is less than 1%, because the last day of curing is immersion curing, the early stage is mainly the compression process of sample pores, and there is no damage. With the increase of strain, the specimen is gradually compressed and destroyed until the damage variable D = 1 and the specimen is completely destroyed. Figure 13a shows that in the early stage of NLS test, before the strain reaches 1%, the sample is in the compaction stage, and the stress changes little with the strain. When the , where σ i and σ (ε i ) is the measured stress and model calculated stress with the strain is ε i , respectively, and m is the number of measured data groups. Figure 13a-d shows the damage evolution curves of NFLS. When D = 0, the sample is not damaged. When D = 1, the sample is completely destroyed. When the strain is less than 1%, because the last day of curing is immersion curing, the early stage is mainly the compression process of sample pores, and there is no damage. With the increase of strain, the specimen is gradually compressed and destroyed until the damage variable D = 1 and the specimen is completely destroyed. Figure 13a shows that in the early stage of NLS test, before the strain reaches 1%, the sample is in the compaction stage, and the stress changes little with the strain. When the strain increases from 1% to 5%, the specimen is mainly in the linear elastic stage and plastic stage, the stress changes obviously with the strain, and the damage accumulates gradually. When the strain is 5%, the damage variable D = 0.8, the sample begins to crack and plastic failure. When the strain is greater than 5%, the damage variable D tends to be stable with the increase of strain until the specimen is completely destroyed. Figure 13b shows the change of damage variable D of NFLS with the content of polypropylene fiber when the content of nano clay is 6%. When the strain is 8%, the damage variable of NFLS without fiber is D = 1 and the sample is completely destroyed. Then, with the increase of fiber content, the destruction speed of the sample slows down. Figure 13c reflects the change of damage variable D of FLS with polypropylene content. When the strain of the sample without fiber is 10%, the damage variable D = 1, and the sample is completely destroyed. With the increase of the content of polypropylene fiber, the failure rate of FLS slows down, and the sample is not completely destroyed until the strain is 12%. Figure 13d shows the curves of the damage variable D of NFLS with the content of nano clay when the content of fiber is fixed at 0.75%. With the increase of nano clay content, the damage curve is dense and the plastic damage is obvious. Due to the addition of polypropylene fiber, the damage of the sample slows down. When the strain is 10%, the sample is not completely damaged. On the one hand, nano clay can fill the pores between soil particles and delay the development of cracks. On the other hand, the fiber plays a reinforcing role, provides tension, and inhibits its plastic deformation.
ually. When the strain is 5%, the damage variable D = 0.8, the sample begins to crack and plastic failure. When the strain is greater than 5%, the damage variable D tends to be stable with the increase of strain until the specimen is completely destroyed. Figure 13b shows the change of damage variable D of NFLS with the content of polypropylene fiber when the content of nano clay is 6%. When the strain is 8%, the damage variable of NFLS without fiber is D = 1 and the sample is completely destroyed. Then, with the increase of fiber content, the destruction speed of the sample slows down. Figure 13c reflects the change of damage variable D of FLS with polypropylene content. When the strain of the sample without fiber is 10%, the damage variable D = 1, and the sample is completely destroyed. With the increase of the content of polypropylene fiber, the failure rate of FLS slows down, and the sample is not completely destroyed until the strain is 12%. Figure 13d shows the curves of the damage variable D of NFLS with the content of nano clay when the content of fiber is fixed at 0.75%. With the increase of nano clay content, the damage curve is dense and the plastic damage is obvious. Due to the addition of polypropylene fiber, the damage of the sample slows down. When the strain is 10%, the sample is not completely damaged. On the one hand, nano clay can fill the pores between soil particles and delay the development of cracks. On the other hand, the fiber plays a reinforcing role, provides tension, and inhibits its plastic deformation.

Conclusions
The elastic modulus and damage stress-strain model of NFLS are studied through UCS and dynamic triaxial test, and the following conclusions can be obtained.
(1) Polypropylene fiber and nano clay can significantly modify the strength of LS. The growth rate of UCS and residual strength of NFLS is the most significant when the content of polypropylene fiber is 0.75% and the content of nano clay is 6%. On the one hand, nano clay can promote the reaction between lime and soil, and on the other hand, it can react with LS to produce cementitious materials such as hydrated calcium silicate and hydrated calcium aluminate. Nano clay and the generated cementitious material can fill the pores between fiber and soil, enhance the interfacial friction between fiber and soil, and improve the strength of NFLS. UCS can be increased by up to 103%. Polypropylene fiber enhances the ductility of NFLS, slows down the development of cracks, and improves its residual strength. The residual strength can be increased by 827%.