Limit Analysis, Numerical and Physical Modelling of Pile Stabilized Slopes

Vast researches have been performed in the eld of earth slope stability analysis including limit equilibrium, strength reduction, and limit analysis methods. All the available methods present slope safety factors in a range with a bit of difference and conrm each other. Validation of analytical results performs with instrumentation in actual slopes existing in the eld. Also, another approach that uses for validating results is experimental modeling. The physical modeling requires manufacturing of the intended model in the laboratory concerning the reducing effect of dimensions on the other parameters. This paper investigates slope safety factors against sliding by simulating the slopes in the laboratory and with an image processing system. The test container has dimensions of 1.5×1.5×2 m. The results have illustrated the crest displacement in reinforced slope, with increasing slope angle 13 degrees, is 30 mm in the experimental test and 13 mm in numerical modeling. In the unreinforced slopes, when the slope angle increased by 8 degrees, the experiment test failed, and the factor of safety in the numerical modeling is less than one. be considered. For modeling an actual concrete pile with a diameter of 800 mm, a section must be selected so that the moment of inertia resulted from applying the intended coefficient is equivalent to the actual model. In the laboratory model, a PVC tube with a diameter of 100 mm and wall thickness of 2 mm placed sandy soil is equivalent to a pile with a diameter of 800 mm in the actual model.


1-Introduction
Slope stability was investigated by different methods. These methods have been categorized into three groups consist of Limit Equilibrium (L.E.), Limit Analysis (L.A.), and Strength Reduction method (S.R.M.).
In the limit equilibrium method, the wedge is divided into several pieces, and the equilibrium of forces and moments, are written in each piece, and nally, the whole wedge balance and safety factor against sliding are determined [1][2][3]. The use of limit analysis methods in actual complex problems is of limitation [4,5]. An upper and lower limit are considered in the limit analysis, which in the upper limit, the assumptions of ow principle and application of energy balance in inter-point vertical and inclined surfaces are considered along the sliding surface [4][5][6][7]. The shear strength reduction method also has been utilized by many researchers for slope stability analysis [8][9][10][11]. Also, wide researches have been conducted in the eld of earth slope stability by using a row of piles and various scientists have studied the behavior of exerted forces to the piles and the slope factor of safety by using analytical relations and numerical modeling [12][13][14][15][16][17][18][19][20][21][22].
Also, many studies perform concerning the experimental investigation of earth slopes. Some cases are mention in the following. Vakil et al. [23] have studied the bearing capacity and settlement of a foundation constructed on an earth slope reinforced with stone columns by experimental modeling and numerical analyses. They used soft clay in their laboratory model and compacted the slope layers in thicknesses of 50 mm. Orenes et al. [24] have studied the factors affecting failure initiation. For this purpose, a series of experimental modeling performs on sandy slopes, and they investigated failure conditions in the slopes by using the arti cial load method and seepage. Lee et al. [31], Shara and Sojoudi [32], Hajiazizi et al. [33], and Mazaheri [33] have been studied the effect of slope stabilization on the concrete pile. Image processing systems have been used for several years in the researcher's studies, especially in studying uid dynamics [25][26]. However, the use of this system is less discusses in Civil and Geotechnical engineering works. White et al. [27] measured soil deformation by conducting the image processing system. Also, they have published a guideline for using this system in soil mechanics [28].
Ould Baba and Peth [29] have examined the suitability of Particle Image Velocimetry (PIV) in investigating creep on slopes by experimental tests. These methods use to validate the modeling results.
As it was mentioned, various researches have been performed in the eld of investigating earth slopes by the use of numerical analyses, analytical relations, and experimental modeling; however, it is not used from the image processing system to monitor slopes movement in the performed studies, up until the present time. Furthermore, previous laboratory tests have been done in boxes with small dimensions, so that the use of large boxes and image processing to examine the slope displacement is remarkable. In the present study, a slope with primary characteristics designed, and the rest of the investigations are proceeding based on the proposed slopes.

2-Statement Of The Problem
This paper is an attempt to make a comparison between numerical analyses, and experimental studies by the employment of laboratory equipment and the image processing system. For this purpose, an earth slope is rst designed based on the existing soil properties that are summarized in Table 1. For determining sand properties in Table 1 Standard test method of ASTM D3080, ASTM D1556, and ASTM C127 have been used. This slope has been depicted in Fig. 1.
The given slope having the dimensions presented in Fig. 1 and the properties summarized in Table 1 is a model in Plaxis and Geo-Slope software. Initial safety factors of 1.35 and 1.37 have been obtained from Plaxis and Slope/w software, respectively. Therefore, it is obvious that the introduced slope is stable in its primary state.

3-Experimental studies
To do experimental studies of the earth slope, a box container with dimensions of 1.5×1.5×2 m has been manufactured. It was used from a 10 mm thick glass in the two large sides of the box to monitor the earth's slope behavior. Regarding the box dimensions, earth slopes with heights up to 12 m can be model with a scale of 1/10. The preparation process of the model is described in the following, step by step.
At rst, according to the soil properties that were presented in Table 1, the earth slope sample was prepared by placing 100 mm layers with dry pluviation technique, and by repeating this procedure and changing the height of pluviating soil material, the required compaction was provided. To reduce friction between the soil material; and the glass wall, a transparent sheet layer was attached to the glass, and it was also used as a lubricant to reduce friction between the transparent sheet and the glass. Moreover, it should be noted that the transparent sheet has been applied strati ed because uni ed applying of the sheet would prevent proper behavior of the soil. In Table 2, the number of tests depicted.
During the preparation of the earth slope, the soil is placed in 100 mm layers with dry pluviation technique, and to control the slope movement, it was used from wooden cylinders with 20 mm of diameter and 50 mm of height in the container, which their base have been grinded and painted with red color. The height of 50 mm is used with the aim that the woods are involved with the soil, completely and their behavior is quite consistent with the soil behavior. These wooden pieces are located in pre-speci ed places marked on the side glass with 100 mm distance from each other. The testing box container and the used wooden cylinders have been illustrated in Fig. 2.
Regarding the large size of the box, the use of weight to create overstress to increase force and reduce safety factor is di cult. From the other side, with regards to being of glass the box sides, the use of weight might be risky work. For this purpose, to reduce the slope safety factor and destabilize, it was decided to increase the slope angle. In this regard, a base along with a crane was improvised on a side of the slope and the slope angle became changeable during the test with the help of a cable attached to the slope side and the crane.
The performed experiments can be divide into two major groups of reinforced and unreinforced earth slopes. In the reinforced slopes, there are using from piles with different distances to stabilize the slopes.
To determine the pile's diameter and stiffness in manufacturing laboratory models, it is required to perform dimension analysis and the moment of inertia and stiffness are required to be selected so that they have an adequate consistency with the actual slope in nature. In this study, the ratio equal to 1/10 was adopted in the preparation of the laboratory model.
According to the available references [30], in the determination of pile stiffness in the laboratory model, the material modulus of elasticity changed by the ratio of , where is introduced equal to 0.5 and 1 for sand and clay soils, respectively. To determine the pile moment of inertia in the laboratory model, it is required that the ratio of be considered. For modeling an actual concrete pile with a diameter of 800 mm, a section must be selected so that the moment of inertia resulted from applying the intended coefficient is equivalent to the actual model. In the laboratory model, a PVC tube with a diameter of 100 mm and wall thickness of 2 mm placed sandy soil is equivalent to a pile with a diameter of 800 mm in the actual model.

3-1-Unreinforced Slope
To investigate the behavior of earth slopes, the earth slope prepares by the above-described procedure, and the wooden cylinders are employed in speci c points with a distance of 100 mm from each other. There is a typical slope shown in Fig. 3. To determine the failure wedge in the slope, and the sliding path, a side of the slope was raised with the help of the crane. During increasing of the slope angle, the whole

4-Numerical Modeling
Two software including Plaxis and Slope/w have been used for numerical analysis of the reinforced and unreinforced slopes prepared in the laboratory. The slope characteristics in the numerical model were considered equal to the laboratory model parameters presented in Table 1. A typical model created in the software is illustrated in Fig. 4.

5-Limit Analysis
Limit Analysis (LA) theorem can be used to examine the safety of earth slope in Geotechnical Engineering. To solve slope stability problems, the use of limit analysis has almost exclusively concentrated on the kinematic theorem [35-37], because under certain assumptions, this is generally simpler to use than the static approach.
Application of the kinematic theorem requires equating the rate of work done by tractions and body forces to the internal energy dissipation rate, for any assumed strain rate field which is governed by the normality rule and is compatible with the velocities at the boundary of the failing soil mass (kinematically admissible failure mechanism). This can be expressed by the following work equation. The rotational log-spiral collapse mechanism, which was earlier examined by Chen [35] and many other researchers, is adopted herein. The geometry of the failure surface ( Fig. 1) is described by the log-spiral equation. movement is recorded by using a camera that is installed in front of the slope. An increase in the slope angle will be continued until the slope collapses.

3-2-Slope Reinforced by Pile
As it was described, it could be used from a 100 mm diameter tube for reinforcing the slope in the prepared laboratory model. In the modeling of slopes reinforced by piles, the distance between piles and the xity of pile end are the most important parameters. To x the pile end, a foam layer is embedded at 200 mm at bottom of the box. This foam layer has two proper functions; one is providing xity of the pile end and the other is preventing stress concentration and avoiding pressure to be exerted on the glass sides during increasing of the slope angle.
In the pile-reinforced slopes, pile distances ranging from 1 to 5 times the pile diameter have been examined, and the whole movement is video recorded during an increase in the slope angle. During all tests, in addition to monitoring the general path of movement, displacement in 5 speci c points including the slope crest, slope heel, slope middle on the ground surface and the slope middle-upper and lower than the sliding surface has been assessed and taken into comparison.
where r 0 = radius of the log-spiral with respect to angle .
The most effective position of the piles within the slope is where the stabilizing force needed to increase the safety factor to the desired value takes the minimum value. The effects of pile location on the required Fp, the force exerted on unit width of sliding mass by the piles, show that if the retaining force that a row of piles can provide is large enough, the piles should be installed near the toe of the slope where the stabilizing force can produce maximum stabilization results. It also indicates that the most economic location for piles in slope stabilization is near the toe of the slope [16].
In this research, a comparison between laboratory model, limit analysis, and numerical analysis with an increase in the slope angle has been depicted in Table 2.

6-Results And Examinations
As it was described in the previous sections, an earth slope was manufactured and analyzed both numerically and experimentally. In this section, the results obtained from investigations are presented in several sub-sections.

6-1-studying slope displacement by the help of image processing system
In the experimental studies, all steps of the work were being photographed. After performing the tests, the recorded images have been processed by using MATLAB software, and the overall displacements of the speci ed points have been evaluated for both reinforced and unreinforced slopes. In Fig. 6, an example of the images processed in the software is shown for different testing times. There are ten images in the gure, which were captured with increasing the earth slope angle. As is obvious in Fig. 6 (image 10), the sliding wedge of the slope could be determined.
Images were also taken during increments of the slope angle for the pile-reinforced earth slope. A comparative illustration presents in Fig. 7 between displacements of the points corresponding to both reinforced and unreinforced earth slopes.
As is clear from the comparison of Figs. 7(a) and 7(b), reinforcement of slopes has a signi cant in uence on reducing slope displacements. Figure 7(a) represents displacements of points in the unreinforced slopes with an increase of 8 degrees in the slope angle. Figure 7(b) demonstrates displacements of points in the pile-reinforced slope with pile distances of 3 times its diameter and an increase of 13 degrees in the slope angle. As is clear, by increasing the reinforced slope angle up to 13 degrees, there is not only no falling occurred, but also its points' displacements are much less than the corresponding displacement of the unreinforced slope with an 8-degree increase of the angle. In the unreinforced slope, a displacement of about 90 mm occurs in the points inside the wedge with a slope rotation equal to 8 degrees; however, in the slope reinforced by 2 piles, a maximum displacement of 30 mm takes place with a slope rotation equal to 13 degrees.

6-2-Comparison of slope safety factors of laboratory and numerical models
In the laboratory model, a slope with before-described characteristics is prepared and its safety factor has been calculated in the software. The slope angle is increased gradually and along with it, displacements are measured in special points at the middle of the slope. Also, in the numerical analysis, slope angle variation imitates the experiment and nally, the obtained results have been taken into comparison.
As it can be observed in Fig. 8(b), displacement at 5 different points in the numerical model (including the slope crest, slope heel, slope middle on the ground surface, the slope middle inside the wedge and near the sliding surface, and the slope middle out of the failure wedge) have been measured. These displacements are shown in Fig. 9 when the slope angle is increased by10 degrees.
As it can be seen in Fig. 9, the slope safety factor reaches a value less than 1 when the slope angle is increased by 10 degrees. This causes the slope to be failed and increases displacement in different points of the slopes. The maximum displacement is related to point A on the top of the slope crest is equal to 80 mm. The difference between displacements of points D and E implies that point D is located inside the wedge, while point E is out of it.

6-3-Comparison of results of laboratory model and numerical analysis in the unreinforced slope
In the following, to investigate and compare results from laboratory model and numerical analysis, the slope displacement and its safety factor have been determined by changing the slope angle. The results present in Table 3. 6-4-Comparison of results of laboratory model and numerical analysis in the reinforced slope As it was described, there have been used from Poly Vinyl Chloride, P.V.C., tubes with a diameter of 100 mm and wall thickness of 2 mm in the laboratory modeling to reinforce the slopes. According to the previous studies [13], the best place of the slope for pile installation is at the middle area and the most appropriate pile performance is when it is xed-end. Therefore, it was used from a 200mm thick foam layer, which the pile end goes down into it, to simulate a xed-end type pile. Thus, in the laboratory models reinforced by piles, all the piles are installed at the middle of the slope and are plunged into the foam layer. The models have been tested with pile distances of 1D, 2D, 3D, and 4D, where D is the pile diameter. The obtained results are compared with the results from numerical analysis of the slopes. Figure 11 shows a typical laboratory model.
As is speci ed in Fig. 11, slopes with different conditions of reinforcement have been prepared and evaluated. Similar to the previous case, the angle of the reinforced slope was also increased to destabilize it; meanwhile, all displacements of the slope have been photographed. The same procedure is also ful lled in the software and their results are compared with each other. In the slopes reinforced by pile, an increase of the slope angle has been continued up to 13 degrees. The results are summarized in Table 4.
It could be found from Table 4 that the minimum displacement of different points of the slope (including the slope crest) occurs when the piles are installed with a distance equal to their diameter.

6-5-Evaluating Slope Stabilization in a Complex Landslides by Concrete Pile; Case Study in the west of Iran
The Iran Gas Trunkline (IGAT) is a series of large diameter pipelines constructed from gas re neries in the south of Iran (Khuzestan and Bushehr provinces) n order to transfer natural gas to consumption centers across the country.
IGAT6, 56 inches (1,420 mm) in diameter, transfers natural gas produced in South Pars phases 6 to 10 from Asalouyeh to Khūzestān Province to be consumed there, in the west of the country, and Iraq. Figure 12 indicates a series of instabilities according to the installing of the underground pipeline. In this gure, there are three sliding slopes. As depicted in Fig. 12, Slope 1 is possible to stabilize by concrete pile.

7-Conclusions
As it was shown in this paper, the use of an image processing system could be of help in studying the behaviors of earth slopes. The major outcomes are presented in the following: When the slope angle increasing 8 degrees, the crest displacement in the unreinforced slope is 90 and 76 mm in experiment test and numerical modeling, respectively. When the slope angle increased 10 degrees in the unreinforced slope, the model was failed in the experimental test, and the factor of safety is 0.98 in numerical modeling.
In the reinforced slope, when the pile distance is ve times of pile diameter, increasing the slope angle 13 degrees the crest displacement is equal to 30 and 13 mm in experiment test and numerical modeling respectively.
Displacements that occurred in different points of the laboratory model have a suitable consistency with numerical modeling; however, in some cases, they got a little greater value. It could correspond to the shaking of the physical model during an increase in the slope angle.    Figure 1 Schematic illustration of the given earth slope studied in this paper Page 13/18 The testing box and pieces used for controlling movement path   Controlling slope displacement using image processing system by the help of MATLAB software  Displacement of the desired points in the numerical analysis   A Typical reinforced slope manufactured in laboratory Figure 12 A Typical proposed reinforced slope by concrete pile in a The Iran Gas Trunkline (IGAT6)