Study on the Calculation Method of the Diffusion Range of Grouting Slurry for a High-speed Railway Ballastless Track Subgrade

Grouting engineering is the most convenient and quick way to deal with postconstruction settlement. Based on the investigation of the Beijing Shenyang high-speed railway in the early stage, the microdestructive force of grout on the soil mass in the process of grouting with the method of fluid solid coupling theory was analyzed. -e numerical simulation results show that, in the grouting process, the greater the grouting pressure, the wider the diffusion range of the slurry and the greater the damage to the microscopic contact force of the slurry to the soil; in addition, when the grouting pressure is below 0.15MPa, the diffusion form of the slurry in the soil is osmotic grouting; when the pressure is about 0.3MPa, the diffusion form of the slurry in the soil is compacted grouting; on the basis of the simulation results, a deeper research is carried out, and the calculation formula of the slurry diffusion curve is obtained. It fills in the blank of the theoretical study of the slurry diffusion curve. In this paper, the diffusion formula of grouting slurry can be used in the theoretical study of clay grouting, and it can be used as a reference for practical engineering.


Introduction
In the recent years, the development of high-speed railway in China has been very rapid. With the development and improvement of the railway system, the problems exposed in the railway operation stage have also been valued by researchers. In the operational stage, high-speed railways will have different levels of subgrade settlement problems due to the different soil properties in different places. As a kind of roadbed settlement treatment method with convenient operation and excellent antisinking effect, grouting engineering has been widely used in the world. However, the research of grouting engineering by domestic and foreign researchers is far behind engineering practice.
is is because the concealment of grouting engineering and the complexity of soil properties have caused researchers to study grouting engineering difficulties. Many researchers have carried out different research studies on the microscopic mechanism of grouting engineering and the macroscopic antisinking effect and have obtained valuable research theories and experience in grouting [1][2][3]. Ho et al. used the finite element software Abaqus to study the mechanical properties of grouting materials under different confining pressures [4]; Rahimi et al. used neural network algorithms to study the waterproofing and penetrating function of the Bakhtiari dam [5]; Moayed et al. measured the shear strength of grouted stone bodies based on the laboratory test [6]; Zhu et al. used experiments to analyze the grouting reinforcement effect of different grouting materials [7]; Li et al. used the multiphysics coupling analysis finite element calculation software Comsol Multiphysics to study the influence of grouting on surface subsidence and groundwater seepage field during foundation pit excavation [8]. Li et al. focused on an improved prediction model to determine the limiting grouting pressure of compaction grouting considering the ground surface upheaval, which is caused by the three-dimensional conical shearing failure.
In summary, due to the limitations of current computing resources and the lack of experimental conditions, most of researchers only conduct research on a certain convenience of grouting engineering. In response to this research phenomenon, the research team carried out a series of research methods on the site inspection, laboratory test, numerical simulation, and field test of the typical section from the railway roadbed grouting project. Based on the work that the research team has carried out, this paper combines the study of the microscopic grouting mechanism with the calculation method of slurry diffusion and deeply studies the interaction between slurry and soil during grouting.

Discrete Element Contact Model.
e simulation can simulate the interaction between particles and fluid accurately EDEM-FLUENT coupling, not only considering the force of the fluid on the particles but also the resistance of the particles to the fluid. e two-way coupling of two software in grouting research has its own unique advantages.
Firstly, the soil contact model to study the damage of the mesocontact force of grouting is selected. EDEM provides a variety of particle contact models for selection, such as the Hertz-Mindlin contact model, Hertz-Mindlin bond model, and Hertz-Mindlin heat transfer. In this paper, the Hertz-Mindlin bonding model (Hertz-Mindlin with bonding model) was used to simulate the soil. e contact force and moment in the model can be decomposed into normal vectors (F n , M n ) and tangential vector (F s , M S ). In the parallel bond model, normal stiffness (N·m 3 ) and tangential stiffness (N·m 3 ) are (k n , k s ) and normal tensile strength (N·m − 2 ) and tangential shear strength (N·m − 2 ) are (σ c , τ c ). e radius of the bonding area is R.
Within one time step, the contact normal and tangential increments resulting from the contact displacement increment are ΔF n � − k n AΔu n , where (A � πR 2 ) is the equivalent disc area, Δu n is the normal direction of displacement, and Δu s is the tangential increment. e normal and tangential increments of the contact torque generated by the incremental rotation are ΔM n � − k n JΔθ n , where J � (1/2)πR 4 is the moment of inertia passing through the contact point for the equivalent disc cross section with the Δθ s direction as the axis, and Δθ n and Δθ s are the direction of rotation increment and tangential increment, respectively. According to the beam theory of material mechanics, the formula for calculating the maximum tensile stress and shear stress acting on the equivalent disk is (3) σ max ≥ σ or τ max ≥ τ, bonding damage. R is the contact radius of the particle, which is set to 3.5 mm in this manuscript.
In EDEM, the particle cohesion between particles by a critical normal stress, the tangential shear bond strength, and bonded disk radius between the particles are determined. If the bond strength between the particles reaches or exceeds this critical value during the simulation, the contact bond will be destroyed.

Fluid-Particle Interaction Force.
In CFD-DEM coupling algorithm, there are some mechanical models used to calculate the particle-fluid interaction force, such as the Drag Model, Lift Model, and Interfacial Heat Transfer Models [9,10]. At the same time, this manuscript focuses on the soil particle-slurry interaction forces which are mainly drag force. erefore, the Di Felice drag model [11] is selected to calculate the particle-fluid interaction force, and it can consider the interaction among particles by introducing volume fractions of single particle, which is used in other similar studies [12,13]. e Di Felice drag model is shown as follows: In the abovementioned formulas, d p is the particle diameter, C D is the drag coefficient of single particle, ρ f is the fluid phase density, ε s is the volume fraction of the particle phase, v f is the fluid phase velocity at the particle center, v s is the particle phase velocity, Re s is the particle Reynolds number.

Discrete Element Contact Model
In this microlevel, CFD-DEM coupling algorithm is used to simulate the whole grouting process. FLUENT is a popular computational fluid dynamics simulation software in the world, and EDEM can be widely used in many scientific research fields, such as mining, geotechnical, agriculture, and so on. e seamless coupling between fluent-EDEM and FLUENT-EDEM is also the world's leading CFD-DEM coupling calculation method. e coupling flow between the two software is shown in Figure 1

Grouting Microscopic Mechanism
Research Model

Calculation Model.
In EDEM, there are many soil particle models, in this paper, we used the Hertz-Mindlin bonding model to simulate soil. Due to the limitation of subgrade size and computer resources, it is impossible to simulate the soil with real graded soil particles in EDEM. erefore, in this simulation, the spherical particles were simulated with spherical particles with a particle size of 6 mm, and a simplified model of 60 cm × 60 cm × 30 cm was used to simulate the subgrade structure. As shown in Figure 2, 599044 spherical particles were filled.

Parameter Settings.
In the calculation process, due to the excessive size of the subgrade structure in this simulation, it is necessary to consider the surrounding rock pressure in the subgrade soil, so the micromechanical parameters are corrected in the simplified model of the subgrade [14]. e particle simulation parameters are shown in Table 1, and the slurry parameters in fluent are shown in Table 2.

Microscopic Grouting Mechanism Analysis.
e grouting is a complex physical process, and its principle is to form stone body by the diffusion of grouting slurry in the roadbed, so as to enhance the bearing capacity of soil and achieve the purpose of settlement control. It can be seen from Figure 3 that when the pressure is small, the movement of the slurry in the soil is in the form of infiltration; when the pressure is high, the slurry appears as compacted grouting in the soil; when the grouting pressure continues to increase, the mode of action of the slurry on the soil is cleft. With the increase of the calculated pressure, the damage degree of the slurry to the soil is also increasing. Figure 4 is an overview of the scope of the extrusion area. Each circle represents the scope of the      Figure 4 clearly depicts the regular increase of extrusion area with the increase of pressure. At 0.5 MPa, the extruded area reaches 32.8 cm, and the range of each extruded area is shown in Figure 5. Comprehensive analysis can be seen in Figures 3-5; the range of splitting increases with the increase of grouting pressure, and the slope of its curve is also increasing which means that the radius changes faster and faster, probably because that the slurry destroyed the soil structure when the grouting pressure was too high. Also, it is observed from Figure 5 that, during the grouting process, the influence of the slurry on the soil appears as an impact shape and is symmetrically distributed in the radial direction of the six grouting holes. Moreover, when the grouting pressure is low, the shape is not obvious, the reason of which is that the effect of the slurry on the soil is small while the pressure is low. At the same time, when the grouting pressure is large, the effect of the slurry on the soil becomes larger, and the response to the slurry on the microscopic level of the soil is also gradually increased. Figure 5 is fitted to the curve to obtain the mathematical curve of the splitting radius and grouting pressure as R � 64.53p 1.331 , and the fitting degree reaches 0.9755, which is a power exponential relationship between the splitting radius and the grouting pressure.

Slurry Diffusion Range Calculation Method.
From the abovementioned range of influence of the slurry on the soil, a very regular change can be obtained. Due to the limitation of space, this paper selects the contact failure map of the soil at the microscopic level of 0.4 MPa for analysis. e contact key damage envelope diagrams of the change graph are described, as shown in Figure 6. en, we take the center of the grouting pipe as the coordinate origin and take the upper part of the curve for analysis and observed that the envelope diagram can be divided into two parts: the first half is a logarithmic curve, and the second half is part of a circle. We select 20 coordinate points in the first half of the curve to export the coordinate points to the EXCEL for processing and obtain the coordinate curve as shown in Figure 7. In Figure 7, the curve can be divided into two parts with the point (X2, Y2) as the boundary.
Fitting the curve of the first half of the curve in Figure 7 is as shown in Figure 8. We got the first half of the curve equation as y � 2.2348 ln x − 1.3719, whose fitting degree reached 0.9803. e second half is like a circle and is fitted with a circular equation. We import the graphic into AUTOCAD and measure the radius of the circle to be 4.36 cm and the coordinates of the center of the circle to be 16.50 cm. So, the standard equation for the envelope curve of the second half is ((x − 16.5) 2 + y 2 � 4.36 2 ), and through the statistical analysis of the envelope curve, the ratio of the diffusion radius R 1 of the first half to the diffusion radius R 2 of the second half is about 3 : 1.
It can be observed from Figure 4 that the envelope curve of the damage range only increases as the pressure increases, but its shape does not change. Unknown constants in the expression are set as a, c, d, and r, so it can be assumed that     Advances in Civil Engineering the mathematical curve expression of the contact failure of the soil contact bond under different grouting pressure is y � a ln x + c (y − d) 2 + x 2 � r 2 , x > 0, y > 0, and the unknown coefficients are a, c, d, and r. erefore, the assumption is made that the unknown coefficients a, c, d, and r are related to the pressure of grouting, and each grouting pressure corresponds to a set of coefficients, corresponding to a set of slurry diffusion curves. So, the next problem that we must solve is finding the relationship between a, c, d, r, and the grouting pressure p.
From Section 5.1, the relationship between grouting pressure and rupture radius is R � 64.53p 1.331 , considering the radius of the grouting pipe and the origin of the coordinate at the center of the grouting pipe. So, the total radius R of the slurry diffusion range mathematical expression is expressed as R � 64.53p 1.331 + b, where b is the radius of the grouting pipe.
In summary, the following equations can be obtained: Among them, R is the total length of the slurry diffusion radius, R 1 is the first half of the slurry diffusion radius, R 2 is the second half of the slurry diffusion radius, b is the radius of the grouting pipe, p is the grouting pressure, and a, c, d, and r are the constants to be solved related to the grouting pressure.
As shown in Figure 9, we select two special points on the curve, select the connection point between the end point and the two-part connection curve, and set the coordinates to be (x1, y1), (x2, y2).
We bring (x1, y1) and (x2, y2) into mathematical equation (5). Since the curve in the graph is a continuous curve, the two parts of the curve are equal in the derivative of (x2, y2), so the following equations are obtained: In equation (6), we consider x 2 and y 2 as known constants and combine the first and second equations to solve a and c first and, then, combine the third equation with the fourth equation to solve for c and d.
In actual engineering, x 1 and y 1 can be directly derived, where y 1 is the radius of the shot hole on the grouting pipe, then x 1 can be calculated by the equation of the grouting pipe, and the calculation equation is x 2 + y 2 � b 2 . In fact, x 2 and y 2 are unknowns. So, we can solve x 2 and y 2 by the following equations combined into a system of equations (8). Also, x 2 , y 2 solved by it are shown in equation (9).
In the equations, erefore, a comprehensive analysis of all the abovementioned equations can be used to obtain the grouting diffusion curve equation as shown in the following equation: where y 1 is the radius of the shot hole, x 1 can be obtained by the equation x 2 + y 2 � b 2 , and x 2 and y 2 are the coordinates of the joint of the two curves, which can be calculated by equation (9). b is the radius of the grouting pipe. R is the radius of the contact failure range of the microscopic layer of the soil, which can be calculated by the equation R � 64.53p 1.331 + b. In this way, this paper obtained the mathematical expression of the grouting slurry diffusion curve under various grouting pressures.

Parameters Research
e relationship between grouting pressure and slurry diffusion has been obtained above. In order to confirm the applicability of the relationship, several pressures in the simulation are brought into the mathematical relationship to verify. is section takes 0.3 MPa as an example to show the whole process of its calculation. First, p � 0. 3  x > 0, y > 0.  Advances in Civil Engineering plotted in a graph, as shown in Figure 10. It can be seen that the diffusion curve is basically consistent with the results obtained in the simulation.

Conclusions
is paper takes a high-speed rail grouting section in actual engineering as an example. After the previous site research work and after discussion by the research group, FLUENT and EDEM coupled simulation were used to analyze the microscopic force of subgrade grouting under microscopic analysis.
e study analyzed the interaction between the slurry and the soil at the microscopic level and in-depth study of the calculation method of the diffusion range of slurry in soil. From the research work, this paper has the following conclusions: (1) During the grouting process, the damage of the slurry to the soil is mainly related to the pressure of grouting. e extent of the damage of the slurry to the microscopic force of the soil increases with the increase of grouting pressure. (2) In this simulation analysis test, when the grouting pressure is below 0.15 MPa, the diffusion form of the slurry in the soil appears as infiltration, and the influence on the soil surrounds about 3 cm around the grouting pipe; when the grouting pressure is about 0.3 MPa, the diffusion form of the slurry in the soil is compacted grouting, and the influence on the soil is about 15 cm around the grouting pipe; when the pressure of grouting exceeds 0.3 MPa, the damage of the slurry to the soil is in the form of splitting, which causes some damage to the structure of the soil. (3) In the study of grouting engineering, the study of its calculation formula is an indispensable part. In this paper, through the calculation and analysis of the grouting microscopic process, the calculation formula of the diffusion envelope curve of the slurry is deduced, and the formula of the slurry diffusion curve is obtained, which lays a foundation for further theoretical research of the subsequent grouting project.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors of this manuscript do not have any conflicts of interest regarding the publication of this article.