Analysis of Soil around Micropile Composite Structures under Horizontal Load

The current approach, which is based on conformal transformation, is to map micropile holes in comparison with unit circle domain. The stress field of soil around a pile plane, as well as the plane strain solution to displacement field distribution, can be obtained by adopting complex variable functions of elastic mechanics. This paper proposes an approach based on Winkler Foundation Beam Model, with the assumption that the soil around the micropiles stemmed from a series of independent springs. The rigidity coefficient of the springs is to be obtained from the planar solution. Based on the deflection curve differential equation of Euler-Bernoulli beams, one can derive the pile deformation and internal force calculationmethod ofmicropile composite structures under horizontal load. In the end, we propose reinforcing highway landslides with micropile composite structure and conducting on-site pile pushing tests. The obtained results from the experiment were then compared with the theoretical approach. It has been indicated through validation analysis that the results obtained from the established theoretical approach display a reasonable degree of accuracy and reliability.


Introduction
Generally, the diameter of a micropile is about 70-300 mm in a small diameter filling pile [1].The slenderness ratio is relatively big.Its preliminary application and exploitation were explored by Fondedile in Italy [2].Micropile composite structures refer to antislide structures that are composed of several, miniature, and single piles with a cap lid at the pile tip, which jointly bears the horizontal load [3].The structure adeptly adapts to shifting terrain during construction with small vibration and noise caused by the construction.It is characterized by a small pile diameter, rapid construction, and flexible piles.Thus, it has been widely used in building reinforcements, shake-proof, foundation underpinning, foundation excavation support, landslide control, and other types of engineering found in buildings [4][5][6].
Previous research on micropile structure mainly discussed ground stabilization, building and rectification, and so forth at vertical load bearing, which was specific to internal force deformation calculations, an analysis of micropiles, and the internal force calculation of a combination of micropile groups [7][8][9].The initial research yielded some achievements.Cantoni et al. [10] proposed a design and calculation method based on reticular micropiles, working under the assumption that the retaining structure would need to be complex; when Macklin designed the anchorage retaining wall, he simplified it as gravity retaining walls in order to analyze the internal force of the micropile [11]; Feng et al. [12] proposed an interaction analysis model for the pile-rock soil and masspiles found in flat micropile systems; they also established a mechanical model to calculate the internal force and deformation of a micropile system using a finite element method; Juran et al. [1] proposed a design approach using a mesh micropile reinforced slope, by assuming that the dense composite strengthening body formed by reticular micropiles, the internal soil mass, and the internal system were not subjected to tensile stress; furthermore, Brown and Shie [13] calculated a pile's internal force by applying a nonlinear, elastic-plastic subgrade reaction method (called the - curve method).However, there was less research on micropile structures under horizontal load.At present, most of the previously discussed engineering designs of micropiles adopt a calculation approach of specific to normal, antislide piles.However, the force model and design calculation theory of micropile composite structure have not been perfected yet.Establishing a calculation mode suitable for micropile composite structure and proposing a reasonable calculation method are urgently needed.Thus, this paper further explores research that capitalizes on this missed opportunity.
This paper discusses an analytical solution to stress and displacement distribution under horizontal load, based on the mechanics theory of two-dimensional elastic complex functions.Using the Winkler Foundation Beam Model, this paper assumes that the soil around a micropile stems from a series of independent springs.The rigidity coefficient of a spring can be obtained using the planar solution.After that, based on the deflection curve of the differential equation of an Euler-Bernoulli beam, the pile deformation and internal force calculation methods of a micropile composite structure under horizontal load can be derived using two modes, namely, by fixing one end, with the other end sliding, as well as fixing both ends.In the end, the paper suggests reinforcing highway landslides using micropile composite structures and conducting on-site pile pushing tests.The results obtained from the experiment have been compared to the theoretical approach to verify the accuracy and reliability of the theoretical approach.

Establishment of the Plane Strain Solution
Model of Micropiles 2.1.Description of the Problem.Then, researching the effects of horizontal load on the micropiles, it is important to consider the internal force and deformation analysis of rigid disc hole structure around the semi-infinite space surface at the  plane (rectangular coordinate system).As shown in Figure 1, the  region is a region of semi-infinite space, except for the disc-structure, and  is the distance between the soil around the pile and the center of the pile;  is the radius of the disc-structure; and  is the horizontal load applied in the center of the pile (  is  component of direction;   is  component of direction).It is assumed that  0 is the displacement under horizontal load applied in the center of the pile.

Fundamental Assumption of the Mechanical Model.
While obtaining the stress field of the soil around the pile and the plane strain solution of the displacement field distribution by adopting the complex variable functions of plane elastic mechanics, the following assumptions are made: (1) While establishing a strain model of two-dimensional complex elastic mechanics, since the rigidity of the microsteel pile is bigger than that around the pile, assume that mini micropiles and cross section around the piles are rigid disc hole structures.(2) Assume that the deformation of the soil around the pile presents a tendency of elastic variation under horizontal load.
(3) Compare the deformation of the micropile and the soil around the pile under horizontal load.

Basic Control Equation.
According to the complex variable function of plane elastic mechanics, the complex analytic function of stress and displacement components at  region are to be expressed in  1 () and  1 (); then [14], where  is any point on the micropile hole,   ,   , and   are the stress and strain component of any point, Re separates the real part and imaginary part in  = 3 − 4, Plane strain; (3 − )/(1 + ), Plane stress, and  = /2(1 + ), where  is modulus of elasticity and  is the Poisson ratio.
The analytic function and expression  1 () and  1 () under stress at infinity bounded conditions are as follows: where ,   , and   , respectively, represent the depth of the tunnel, the bulk density of the surrounding rock, and the relationship of the lateral pressure coefficient and the surrounding rock stress.In this equation,  is positive integer.
According to the analytic function under complex variables functions at bounded conditions, the displacement and stress boundary conditions are as follows: The condition of stress boundary: 2.4.Conformal Mapping.One must map the micropile orifice to unit annulus.The semi-infinite space orifice within  plane can be converted to within the plane  after conformal mapping as shown in Figure 2.Then, the mapping function of zone  after conformal mapping is shown in [15].Therefore, where any point  located in the  plane is expressed as polar coordinates; that is,  =   and  is a parameter determined by  and : (5) After transforming the function at plane  to the function of , the stress and displacement component will change to [16] where  = (), Φ() =   1 (), and Ψ() =   1 ().From formula (6), it is detected that, in order to get the stress and displacement values through using the fundamental equation, tone has to solve the complex function and obtain the solution to  0 () and  0 ().

Solution to Displacement 𝑢 0 under Horizontal Load. (1)
According to the plane displacement boundary conditions of micropile-soil structure, the following can be seen from the assumed conditions.
The displacement boundary conditions can be divided as follows: infinity to pile core position, namely,  = ; so the soil mass will not cause direct impact on pile deformation, whereas displacement exists at the contact surface of the pile hole boundary and soil mass, namely, | + | = .The strain of the micropiles soil structure to plane  = 3−4 can be seen in formula (3).In other words, the displacement boundary is as follows: where  1 is distance between the soil around the pile and the center of the pile and  0 is the boundary curve of the pile hole.
When the function of the  plane is transformed into a function of the  plane,  = (), Φ() =   1 (), and Ψ() =   1 () are substituted into formula (2) to get the following: where  0 () = ∑  =0    − and  0 () = ∑  =0    − .From the given situation, one can obtain ∑   = , ∑   = 0, and  =  =   =   = 0.When the initial conditions are substituted into formula (8), one can get the following: By substituting formula (9) into the displacement boundary condition, formula (7), it can be derived that where Through the simultaneous application of the boundary conditions of the two equations found in formula (10), the expression of  0 () and  0 () can be derived (including the unknown displacement of  0 ).
Similarly, according to formula (12), the value of  0 () can be obtained by taking conjugation at both sides and then applying the Cauchy integral operator.Thus, the information can be simultaneously obtained by combining formula (10), to get the value of  0 .
2.6.Solution to  0 () and  0 ().To obtain the boundary conditions of a miniature pile, make  =  and unfold () and () in a Laurent series form; namely [17], Depending on formula (4), one can get the following: Using formulas ( 14) and ( 15) substituted into boundary conditions  1 of formula (7), one can get the following: Make () = () = 2( + V) and  * () = (1 − )() = ∑ ∞      .Formulas ( 14), (15), and ( 17)-( 19) can be substituted into the boundary conditions  1 of formula (7); by eliminating   and   , one finds   and   : We simultaneously solved the four formulas in (17) to get all of the coefficients except for  0 , since  0 represents rigid body displacement; thus, no stress will be generated and it can be deemed as 0. By  times of iteration, one can obtain   ; thus, to get  0 () and  0 (), one must substitute  0 () and  0 () in the complex function of the fundamental equation, to ascertain the stress and displacement field of the soil around a pile of micropiles under horizontal load.According to the assumed conditions specified in Sections 2.1 and 2.2, the relationship between horizontal load and the horizontal displacement can be obtained from  = / 0 .

Establishment of Pile-Soil Mechanics Model Based on Winkler Foundation Beam
In most cases, while ministeel tub piles under horizontal load are applied in landslide control and slope reinforcement, a pile groups' layout will be adopted.Especially for multistage slopes, platforms are available for each grade.Miniature piles can be set at the platform for reinforcement.The multiple, miniature piles that are exposed on the platform can be fixed with a top beam to enhance their sliding resistance force.The micropile layout is shown in Figure 3.

Model Assumption.
While miniature piles are applied for landslide reinforcement, the mechanism can be used to address bigger shearing resistance in miniature pile landslides.Pile-soil structures will resist landslide thrust formed by micropiles as well as the soil mass.For internal calculation, the Winkler low econometric model is adopted.This research assumes that the soil around the micropile has a series of discrete springs and a rigidity coefficient of  = / 0 .As shown in Figure 3, since the micropile is fixed at the bed rock below the sliding surface, it is assumed that the sliding surface is as fixed constraint.Due to the fact that the micropile is fixed and connected through a top beam, as compared to miniature pile, the top beam can be deemed as a rigid member.Under the effects of horizontal load, the top beam only experiences horizontal displacement, so the displacement at each pile cap will be the same [18].Thus, the internal force calculation of the model is shown in Figure 4.

Internal Force Calculation.
According to the fundamental theory of elastic mechanics, the stress and calculation model of a single micropile can be obtained as shown in Figure 5.It is composed of two calculation models, namely, with one fixed and one sliding end under a concentrated load (Figure 5

Solution to Calculation Model with One Fixed and
One Sliding End under a Concentrated Load.According to the internal force calculation model specified in Figure 5(a), when considering micropile structures with one fixed end and one sliding end under a concentrated load, the boundary  condition on the top of the micropile is  = 0; the bending moment is  = 0; and  is horizontal shear.The boundary condition on the bottom of the micropile is  = , the angle is  = 0, and the horizontal displacement is  = 0.The flexural differential equation of an Euler-Bernoulli beam with one fixed and one sliding end is [19]   4   4 +  () = 0.
In the formula above,  is the elasticity modulus of the miniature pile,  is the inertia moment of a micropile cross section,  is the rigidity coefficient of an assumed spring, and () is the horizontal displacement of the soil around the pile along the pile body.By substituting the boundary conditions in the flexural differential equation ( 18), the analytical solution can be obtained: where  = (/4) 1/4 and  1 ,  2 ,  3 , and  4 are the integral constant.where () is the uniform load of the soil around the pile.To solve this equation,  = (/4) 1/4 can be substituted into formula (19) to get the following:

Solution to Calculation Model with Both Ends
By substituting the boundary conditions into the flexural differential equation ( 18), the analytical solution can

Steel bar welded to
The type of bracket  ( Figure 9: The horizontal static load test of micropiles.

The Determination of the Elastic Modulus (𝐸).
In the previous calculation,  is a crucial parameter, usually determined by laboratory experiment and experimental relationship of previous practice.However, as laboratory experiment, it is very difficult to reproduce the soil pile in live load effect on the stress path, and experiential relationship cannot determine the value of  either as it can be influenced obviously by human interference.Due to the above reason, in this paper, the field test of single pile at the scene of the elastic deformation range was applied.The elastic model value is the inversion analysis through the stress-strain curve, with the basic process described as follows.
Under the load, the micropile is loaded step by step, using the same load increment.Record each time load displacement deformation, therefore, according to the stress-strain curve obtained from test to calculate the value of elastic modulus.

Engineering Practice Analysis and Verification of Pile-Soil Mechanics Model
4.1.Engineering Examples.In this example, the landslide of a highway is medium type, about 80-100 m wild and 160 m long.The front is about 5 m thick, whereas the central part is 10-15 m thick and the back part is 3.5-8.0m thick.The average thickness of the landslide is about 10 m and the total volume of landslide is about 160,000 m 3 .The landslide demonstrates typical factors with the slope facing the empty and the steep slope.The features of landslide appearance are obvious: the back has tensile cracks and the two sides have pinnate cracks.The front is flanked by a shear seam with a ballooning extrusion crack forming and radioactive cracks.Specific to the characteristics of this landslide, it is proposed to adopt a light, retaining structure with a ministeel tub piles composite structure for reinforcement.

Establishment of Experiment Process and Numerical Simulation.
According to the design requirements, we performed on-site horizontal static load tests within the landslide reinforcement range.We then simulated an action mechanism and stress distribution in the micropile composite structure by adopting a jack to provide load.Through earth pressure cells that were installed before and after the micropile, we were able to monitor the earth pressure variation and the landslide thrust while checking the status of each row of piles.The stress of the miniature piles is measured through a reinforcement meter that was welded on a miniature steel pipe.The pile body deformation was indirectly measured through an inclinometer that was installed on the pile sides.
Figure 6 shows the model demonstration diagram.

Experiment Process.
The in situ test uses the micropile to determine the rate of reinforcement.The grouting coagulation of the soil strength grade is C25.The micropile length indicates the landslide segment, which is 8 m long.The pile diameter is 150 mm and the tube diameter is 50 mm.The main reinforcement pile contains 3 roots, made up of 28 reinforced steel pipes.The microcap sets the C30 concrete capping beam with a beam that is 0.5 m high and 1.5 m wide.In accordance with the requirements for the load test, the test is conducted with a grade 11 effective load, using two gauges to record the data.The average is taken as the final result.On the 12th level (96 t) load, the counterforce device becomes damaged, indicating the end of the test.The load of the destruction is the horizontal limit load.At the moment when the maximum amount of bending occurs, the steel of the yielding tensile zone is the corresponding load.The micropile section and test process are shown in Figures 7, 8, and 9.

Result Analysis.
The rate of pile displacement of each grade of load was obtained through on-site experiments and analysis, as shown in Figures 10 and 11.In Figures 10  and 11, it can be detected that, under horizontal load, the displacement of the micropile composite structure above the sliding surface (3 m) is more obvious than that below the sliding surface, which indicates that the micropile composite structure presents a tendency to lean forward.Since the pile only leaned 2 m forward at the base below the sliding surface, the horizontal displacement is basically 0, which indicates that the anchorage effects at the anchorage section are comparatively better.When consistently increasing the horizontal load, the variable quantities of the displacement of pile top of each row of piles are the same.This consistency in the variable quantities is caused by the lid cap contracting to make the micropile composite structure act as a whole for antisliding by adopting a jack to provide step load.The bending moment and the displacement distribution laws of the three-row piles are similar.A bending moment above 0.25 m is 0. No bending deflection of the micropile is generated, due to the constraint     of the cap lid.The point of contraflexure occurred at 0.9 m above the sliding surface and 2.5 m below the sliding surface for the three rows of piles.The maximum sagging moment is at 1 m above the sliding surface and the maximum hogging moment is at 0.5 m below the sliding surface.The value of bending moment increases along with the horizontal load.
As the horizontal load increased to grades 9-120 KPa (close to the upper limit of the horizontal load of an antislide structure of the micropile combination), the increment of the maximum bending moment (the absolute value) at the third row will be the maximum, followed by that of the second row and then the first row.If each row of piles adopts the same bending strength design without considering the impact of plastic failure on the soil mass between the piles, under the effects of landslide thrust beyond the upper limit of the horizontal load of the micropile combined mechanism, then the sequence for each row of piles is the third row, followed by the second row and the first row.11 and 13, it can be seen that, according to the pile-soil response theory calculation method under horizontal load, the pile displacement and bending moment are similar to the results found in simulated field experiments under each grade of load, which shows that the theoretical approach demonstrated in this paper is feasible.From the figures, it can be detected that the results of the experiment are smaller than the bending moment and the rate of pile body displacement.This result occurs because, in most cases, while searching for solution to the Winkler Foundation Beam Model, the sheer force of the soil between piles is generally ignored.But for engineering design, the solution from theoretical approach adopted by the paper is simply safe; thus, it can satisfy design accuracy.

Conclusion
(1) This paper discusses an analytical solution to stress around a micropile.Displacement distribution under horizontal load is obtained according to the planar complex function theory.The defects found in the internal force analysis after applying uniform sections were addressed.The stress and displacement distribution were analyzed by utilizing the proposed analytical solution in the soil around the micropiles.
(2) Based on the Winkler Foundation Beam Model, we assumed that the soil around the micropiles stemmed from a series of independent springs.The rigidity coefficient of a spring is obtained using a planar solution.After that, based on the deflection curve differential equation of an Euler-Bernoulli beam, the pile deformation and internal force calculation methods of micropile composite structures under horizontal loads can be derived using two modes, which have provided theoretical guidance for engineering designs.On the one hand, one end is fixed with the other end sliding; on the other hand, both ends are fixed.
(3) By comparing the results obtained from on-site pile pushing tests, it is determined that the results obtained of the established theoretical approach are reasonably accurate and reliable.

Figure 3 :
Figure 3: Micropile composite structure set on step of slope.
(a)) and with both ends fixed under a uniformly distributed load (Figure5(b)).Specific to the calculation model found in Figure4, for three-row micropile structure, the internal force of the micropile AD is obtained through superposition of Figures5(a) and 5(b).The internal forces of micropile BE and CF are derived from Figure 5(a).

Figure 6 :
Figure 6: The schematic drawing of a test model of micropiles composite structure.
The third row of piles

Figure 10 : 3
Figure 10: Comparative curves of distribution of pile deflection.

Figure 11 :
Figure 11: Comparative curves of distribution of pile model.

Figure 12 :
Figure 12: Comparative curves of the distribution of pile deflection.
Theoretical Calculation.The proposed approach incorporated the theoretical calculation model of Sections 2 and 3 and the mechanical parameters and conditions of the experiment in order to calculate the displacement of each row of micropiles under 26.6 KPa, 79.8 KPa, and 119.7 KPa, as well as the bending moment under 26.6 KPa and 119.7 KPa, as shown in Figures12 and 13.By comparing Figures 10 and 12 to Figures