Abstract
Non-water reacting polyurethane is a new type of grouting material for support reinforcement and emergency rescue, whose creep mechanical behavior is directly related to the stability and safety of the reinforcement project. To investigate the creep properties of non-water reacting polyurethane grouting material, a series of multistage loading creep tests were performed under uniaxial compression conditions. Then, the creep curves, creep strain rate, and creep failure modes of non-water reacting polyurethane grouting material were obtained. In addition, the effects of temperature and density on the creep properties of non-water reacting polyurethane grouting material were also analyzed. Finally, based on the finite difference program, the numerical simulation of the uniaxial compression creep test of non-water reacting polyurethane grouting material was conducted by employing the Burgers–Mohr creep model. The simulated creep curve was in good agreement with the test curve, which verified the feasibility of the model to describe the creep characteristics of non-water reacting polyurethane grouting material. The research results can provide a reference for studying the creep behavior of polyurethane grouting material and also contribute to the development and application of polyurethane material grouting technology.
1. Introduction
With the development of underground engineering, water inrush disasters occur frequently. The grouting technology is considered to be an effective means to prevent and control these engineering accidents [1–10]. Non-water reacting two-component polyurethane is a new type of grouting material, which is different from the traditional grouting material represented by cement mortar. It has the advantages of rapid hardening, large expansion, and a small permeability coefficient [11–15]. In recent years, it has gradually become a hot topic of international research. At present, it has been widely used in seepage prevention, infrastructure reinforcement, and as support materials in infrastructure engineering projects [16–24].
As non-water reacting polyurethane is a new type of material gradually emerging in recent years, the research on it is limited by test equipment and experimental data. Scholars have made slow progress in the study of its physical and mechanical properties, and the theoretical research is far behind the engineering practice [25–28]. Liu [29] conducted a uniaxial compression test and splitting test for low-density polyurethane materials (ρ < 0.3 g/cm3), and obtained the relationship between the density of polyurethane materials and the uniaxial compressive strength and splitting tensile strength. Pan [30] tested the expansion performance of non-water reacting polyurethane grouting material and obtained the law that the expansion force increases exponentially with the increase of density in fixed volume space. Shi et al. [31] studied the permeability, bending, and temperature characteristics of non-water reacting polyurethane grouting material and pointed out that the permeability of polyurethane materials decreases with the increase of density. Liu [32] tested the tensile properties of non-water reacting polyurethane grouting material by direct tensile test. The results showed that the tensile strength increased exponentially with the increase in density. Aiming at non-water reacting polyurethane grouting material used for repairing highways, settlement lifting for ballastless track of high-speed rail, and impermeable reinforcement of dams or dikes in a seismic region, Li et al. [33] investigated dynamic viscoelastic properties of polyurethane grouting material and proposed a viscoelastic constitutive model based on the generalized Maxwell model. Sham [34] studied the relationship between the viscoelasticity of polyurethane and frequency or temperature by using the Dynamic Mechanical Thermal Analysis method. Bagley and Torvik [35] studied creep and relaxation behaviors of polyurethane materials and proposed a fractional viscoelastic model.
However, despite the recent improvements in the research on the mechanical properties of non-water reacting polyurethane grouting material by many scholars, the current research mainly focuses on the conventional mechanical behavior of polyurethane materials, and there are only a few research reports on the creep mechanical behavior. In practical engineering, the creep phenomenon of materials widely exists in all kinds of engineering. A large number of field measurements and engineering practices show that the grouting materials in highway and railway tunnels, large slopes, mine roadways, and tunnel surrounding rocks exhibit creep characteristics [12–14]. Therefore, it is of great significance to study the creep characteristics of non-water reacting polyurethane materials. The research results can provide a scientific basis for the analysis of long-term deformation characteristics of subgrade and provide theoretical guidance for jacking settlement pipeline and tunnel support.
Due to the similarity of the mechanical properties of non-water reacting polyurethane grouting material with geotechnical materials and asphalt mixtures, this study employed the rheological mechanical properties of geotechnical materials and asphalt mixtures to study the creep mechanical properties of non-water reacting polyurethane grouting material. First, a series of multistage loading creep tests were carried out to obtain the creep curves of non-water reacting polyurethane grouting material. Then, the effects of density and temperature on the creep rate and creep failure mode of polyurethane grouting material were analyzed. Finally, the viscoelastic mechanical parameters of polyurethane grouting material were obtained by the fitting method, and the creep test curve of polyurethane material was simulated by a finite difference program. The reliability of the Burgers–Mohr creep model to describe the creep behavior of non-water reacting polyurethane grouting material was verified.
2. Creep Tests
2.1. Testing Equipment
In our study, the uniaxial compression creep tests were conducted using the CMT5105 series microcomputer-controlled electronic universal testing system, as shown in Figure 1. The test system has high accuracy and sensitivity for the measurement and control of load, deformation, and displacement. It can also carry out the automatic control test of constant velocity loading, constant velocity deformation, and constant velocity displacement, and has the functions of low-frequency load cycle, deformation cycle, and displacement cycle. It can be competent for the tensile, compression, bending, shear, peeling, tearing, and other tests of various solid material samples. The maximum test force of the test equipment is 100 kN and the resolution of the test force is 1/300000 FS (full scale). The deformation measurement range is 0.2 ∼ 100% FS and the relative error of deformation indication is ±0.5%. In order to study the temperature effect on creep characteristics of non-water reacting polyurethane, a high-low temperature test chamber was equipped with a temperature control range of −40°C to 150°C.

2.2. Sample Preparation
The non-water reacting polyurethane material is a two-component foam (polyols and isocyanate) that exhibits liquid properties before curing, with fluidity and expansion. The specimen used for the creep test was a cylinder with a diameter of 50 mm and a height of 50 or 100 mm. Therefore, the sample-making mold was specially developed, as shown in Figure 2.

Since the non-water reacting polyurethane grouting material can generate a large expansion force after mixing, the mold was made of steel to ensure sufficient rigidity. The grouting hole was set in the steel plate above the mold. In order to facilitate demolding, it is necessary to apply grease to the inner wall of the mold before grouting. It should be noted that the specimens with different densities were prepared by controlling the injection volume of polyols and isocyanate.
Due to the limitation of the testing equipment, some cylindrical samples with a height of 100 mm were cut into two cylindrical samples with a height of 50 mm by the numerical control punching machine. In order to ensure the smoothness of the end face of the specimen, the specimen with an uneven end face was polished with sandpaper. The surface deviation of the two ends was within a range of 0.05 mm. Figure 3 shows some of the prepared specimens.

The mass, height, and diameter of the specimen were measured by a digital display balance and vernier caliper with an accuracy of 0.1 g, and then the density of the polyurethane specimen was calculated. The uniaxial compressive strength of each specimen was calculated according to the empirical formula, as listed in Table 1. The empirical formula for uniaxial compressive strength is as follows [32]:where x denotes the density of the polyurethane specimen and y denotes the uniaxial compressive strength of the polyurethane specimen.
2.3. Creep Test Scheme
Previous studies have shown that density and temperature are important factors affecting the mechanical properties of polyurethane materials. Therefore, the effects of density and temperature were considered in this study. For creep tests to study the density effect, the multistep loading test method was used to reduce the influence of sample inhomogeneity. The axial stress was divided into several different levels. At present, no one has studied the creep mechanical properties of non-water reacting polyurethane grouting materials by laboratory tests. We also have no relevant literature for reference. Through a large number of tests, it is found that non-aqueous reactive polyurethane grouting materials have similar viscosity characteristics to asphalt mixture. Therefore, we searched the literature and relevant test procedures and referred to some test methods for the study of the rheological properties of asphalt mixture [25, 36–38]. The loading time for each step was tentatively determined to be approximately 1800 s (30 minutes). When the loading time of each step reached the predetermined value, the next step of the axial stress level was applied. The axial stress was increased step by step until the sample failed. In order to be consistent with the engineering environment as much as possible, the local average temperature (25°C) was taken as the reference temperature for creep tests. When considering the effect of temperature on the creep properties of the polyurethane grouting material, 25°C was taken as the initial temperature; the temperature increased or decreased by 10°C for each stage until the sample failed. It should be noted that the axial stress remained constant since loading throughout the creep experiment to study the temperature effect. After the temperature of each stage lasted for 1800 s, the temperature of the next stage was applied.
All of the creep test data were recorded during the test and the creep test scheme is outlined in Table 2. Before carrying out each creep test, a load of 50 N shall be pre-applied to the specimen at a loading speed of 1 mm/min, and the load shall be kept constant for 30 s.
3. Test Results Analysis
3.1. Creep Strain
Figure 4 illustrates the relationship curves of axial strain with time for the polyurethane grouting material. From Figure 4(a), it can be observed that the axial strain curves of polyurethane samples with different densities show typical creep characteristics, including attenuation creep, steady-state creep, and accelerated creep. When the axial stress was small, the polyurethane grouting material showed attenuation creep and steady-state creep behavior. With the gradual increase of axial stress, the creep of polyurethane grouting material gradually transited to the accelerated creep stage. In addition, polyurethane samples with higher density experienced a longer creep process. Taking the creep curves with densities of 0.141 and 0.357 g/cm3 as an example, the former creep loading time lasted for about 5800 s, while the latter lasted for nearly 11000 s, which was more than twice that of the former. Moreover, polyurethane samples with higher density experienced a higher creep failure strength. For lower density specimens (less than 0.233 g/cm3), the averaged creep failure strength is about 0.76 UCS. However, for the specimens with higher density (more than 0.233 g/cm3), their creep failure strength is basically above 0.9 UCS. This is because the polyurethane samples with higher densities showed larger strength, and thus the loading time was longer to reach higher creep failure strength. Therefore, in engineering practice, its service life can be improved by increasing the density of the polyurethane grouting material. Furthermore, it is also found that when the axial stress was less than 50% of the compressive strength, the polyurethane grouting material showed attenuated creep behavior. When it was 95% of the compressive strength, the polyurethane grouting material showed accelerated creep behavior. While it was 50% ∼ 95% of the compressive strength, the polyurethane grouting material showed steady-state creep behavior.

(a)

(b)

(c)
Polyurethane samples with densities of 0.408 and 0.567 g/cm3 were taken as examples. Table 3 lists the axial total strain εt, axial instantaneous strain εi, axial creep strain εc, and at instantaneous elastic modulus Ei different loading levels during the creep test.
From Table 3, the contribution of creep to axial total deformation increased with the increasing axial stress. In particular, the creep strain was less than 30% of the amount of axial total strain for the first loading level, when the axial creep strain exceeded 90% of the amount of axial total strain for the last loading level. In addition, it is also noticed that the elastic modulus of the polyurethane sample first increased and then decreased during the whole creep process. This is because the polyurethane grouting material had many internal voids, and with the increase of the axial stress, the elastic modulus of the polyurethane sample reached the maximum. Then, as the axial stress continued to increase, a certain degree of damage occurred inside the polyurethane sample, resulting in the deterioration of its stiffness. On the whole, the density had a great influence on the elastic modulus of polyurethane grouting material. The average elastic modulus of a polyurethane sample with a density of 0.567 g/cm3 was nearly twice larger than that of a polyurethane sample with a density of 0.408 g/cm3.
3.2. Creep Strain Rate
In order to study the law of creep strain rate, Figure 5 shows the variation curve of axial strain with time during each loading level for polyurethane samples with different densities. It is clearly seen that the polyurethane grouting material exhibits an obvious creep stress threshold, that is, the creep phenomenon occurs only after the axial stress reaches a certain value. Moreover, the lower the density, the lower the creep stress threshold of polyurethane samples. For specimens with lower density (less than 0.295 g/cm3), when the axial stress is lower than 0.4 UCS, no creep phenomenon will occur. However, for specimens with higher density (more than 0.295 g/cm3), the creep stress threshold is about 0.5 ∼ 0.6 UCS. In addition, it can also be seen that all specimens entered the accelerated creep stage before reaching the expected loading time (1800 s) during the final loading level. At the same time, the specimen did not fail when the axial load was applied. The final failure of the specimens occurred within 1800 s of the final loading level, i.e., creep failure. Moreover, the creep strain rate in the steady-state stage showed an increasing trend with the increase of density. However, when the density of the sample was too large, due to the compression of the voids in the sample, the deformation ability decreased, resulting in a decrease in the creep rate in the steady-state creep stage. Furthermore, with the increase of the axial stress, the creep strain rate in the steady-state creep stage also showed an increasing trend. The creep strain rate in the steady-state creep stage was usually 4% to 21% of the creep rate in the accelerated creep stage.

(a)

(b)

(c)

(d)

(e)

(f)
3.3. Creep Failure Characteristics
Figures 6∼10 illustrate the failure forms of non-water reacting polyurethane samples.

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)
It can be found that the creep failure modes of polyurethane specimens with different densities were significantly different. When the sample density was low, the sample was mainly a compressed deformation failure (see Figures 6 and 7). In addition, due to the difference in density and hardness, there were differences in the manifestations of compression failure of the specimens. The deformation of the two groups of samples with densities of 0.141 and 0.194 g/cm3 was only the wrinkle on the sample surface, without fracture failure on the sample surface (see Figure 7). For the sample with a density of 0.233 g/cm3, in addition to the same wrinkle, some wrinkles on the surface of the sample were broken (see Figure 7). When the density of the specimen was larger, there was a macro fracture surface. Vertical cracks appeared on the side surface of the sample with a density of 0.407 g/cm3. It indicates that the specimen occurred splitting tensile failure. For the samples with a density of 0.567 and 0.585 g/cm3, in addition to the vertical crack, there was a shear fracture surface with an inclination of 45°. When the density was large, the polyurethane sample showed tensile-shear composite failure.
3.4. Temperature Effect
Figure 11 shows the variation curves of axial strain with time for polyurethane samples with different densities under heating and cooling conditions. Clearly, the temperature had a significant effect on the creep behavior of polyurethane grouting material. When the temperature rose, the creep rate of the polyurethane sample was larger, and there was steady-state creep and accelerated creep stages. While the temperature decreased, the creep rate also decreased, and the polyurethane specimens only exhibited attenuated creep behavior. Therefore, it can be concluded that high temperature can enhance the creep of polyurethane grouting material, while low temperature can inhibit its creep behavior. The reason for this phenomenon is probably that the increase in temperature will increase the activity of the molecules in the polyurethane grouting material, causing the creep rate to increase significantly and thus inducing the accelerated creep behavior. In the low-temperature environment, the molecular activity inside the polyurethane grouting material was low, and thus the creep rate was greatly reduced, and only attenuated creep behavior occurs. In addition, by comparing Figures 11(a) and 11(b), it is found that the polyurethane specimen with a density of 0.231 g/cm3 had accelerated creep at 35°C, while the specimen with a density of 0.141 g/cm3 was only in the steady-state creep stage at 35°C, and the accelerated creep behavior occurred at 45°C. This further verified the above findings. The specimen with a smaller density contained less polyurethane grouting material, and the active energy of internal molecules was lower at the same temperature. Therefore, the temperature to reach the accelerated creep stage should be higher.

(a)

(b)

(c)
4. Numerical Simulation
4.1. Burgers–Mohr Creep Model
From the above analysis, the polyurethane grouting material had the three-stage characteristics of a typical creep curve, i.e., attenuation creep, steady-state creep, and accelerated creep. The Burgers creep model can well describe the viscoelastic characteristics of attenuation creep and steady-state creep. In order to better describe the accelerated creep characteristics, the FLAC3D finite difference program provided a Burgers-Mohr creep model which connects a Mohr–Coulomb plastic cell in series with the Burgers model [39], as shown in Figure 12. The Burgers part of the model was used to describe the first two stages of the creep curve. When the deviatoric stress exceeded yield strength and entered the accelerated creep stage, the Mohr–Coulomb plastic cell started to work. Next, the feasibility of using this model was explored to describe the creep mechanical behavior of polyurethane grouting material under uniaxial compression conditions. For the convenience of illustration, the polyurethane specimen with a density of 0.408 g/cm3 was taken as an example.

4.2. Model Parameter Solution
From Figure 12, it is found that the Burgers–Mohr creep model includes five mechanical parameters, i.e., Kelvin module E1, Kelvin visibility coefficient η1, Maxwell module E2, Maxwell visibility coefficient η2, and yield strength σs. Since this study was limited to uniaxial compression, the yield strength here was taken as uniaxial compressive strength (9.11 MPa). In addition, in order to describe the lateral deformation, the mechanical parameter Poisson’s ratio was also needed. According to the uniaxial compression test results, the Poisson’s ratio was 0.25.
Under the one-dimensional condition, the creep equation of the Burgers creep model can be written aswhere σ0 denotes the current stress level and t denotes the loading time.
Referring to the previous study [40], (2) was adopted to fit the axial strain curve under each level of stress level, and the fitted model parameters are listed in Table 4. The creep parameters of each stress level varied greatly, indicating that the creep parameters of polyurethane grouting material were non-stationary in the whole creep test process.
Figure 13 shows the comparison between the fitted curves and the test value. It is found that the fitted curves were in good agreement with the test data at each stress level, and the correlation coefficient was close to 1. It indicates that the Burgers creep model can well describe the viscoelastic characteristics of polyurethane grouting material under uniaxial compression creep.

(a)

(b)

(c)

(d)
4.3. Model Verification
According to the parameters in Table 4, the Burgers–Mohr creep model built-in FLAC3D was employed to simulate the uniaxial compression creep of the specimen with a density of 0.408 g/cm3. The loading level classification and creep time in the numerical simulation were completely consistent with the actual test. Figure 14 shows the comparison between the test curve and the simulated curve. It is observed that the simulated curve was in good agreement with the experimental curve, and the accelerated creep stage was successfully obtained. It has strongly proved that the Burgers–Mohr creep model can better describe the viscoelastic plastic creep characteristics of polyurethane grouting material under uniaxial compression conditions.

5. Conclusions
In this study, a series of uniaxial compression creep tests were performed to study the creep properties of polyurethane grouting material. The main conclusions are as follows:(1)Polyurethane grouting material showed typical creep characteristics under uniaxial compression. Density is the dominant factor affecting the creep behavior of polyurethane grouting material. Higher density can significantly increase the creep failure strength, stiffness, and creep failure time of polyurethane grouting material.(2)The contribution of creep to polyurethane grouting material deformation increased with the increasing axial stress level. The steady-state creep strain rate increased with the increase of the axial stress level.(3)The creep failure mode of polyurethane grouting material under uniaxial compression was significantly affected by the density. When the density was small (less than 0.2 g/cm3), the polyurethane sample showed the compression deformation failure mode with wrinkles on the surface. When the density was medium (0.2 ∼ 0.5 g/cm3), the polyurethane sample shows splitting failure. When the density was high (more than 0.5 g/cm3), the polyurethane sample showed tensile-shear composite failure.(4)The creep characteristics of polyurethane grouting material were significantly dependent on temperature. High temperature (more than 35°C) can enhance the creep of polyurethane grouting material, while low temperature (less than 35°C) can inhibit its creep behavior.(5)Burgers–Mohr creep model can better describe the attenuated creep, steady-state creep, and accelerate the creep of polyurethane grouting material under uniaxial compression conditions. It will provide a reference for the research on the creep behavior of polyurethane grouting material, and contribute to the development and application of polyurethane grouting material grouting technology.
Data Availability
The data presented in this study are available on request to the corresponding author. The data are not publicly available due to the data also form part of an ongoing study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was funded by the Scientific Research Project of Education Department of Henan Province of China (Grant no. 22A56008); Science and Technology Research Project of Department of Science and Technology of Henan Province of China (Grant no. 212102310965); and Young Backbone Teacher Training Project of Henan Urban Construction University (Grant no. YCJQNGGJS202202).