Research Article | Open Access

# Deformation Analysis of Granular Soils under Dynamic Compaction Based on Stochastic Medium Theory

**Academic Editor:**Giovanni Garcea

#### Abstract

Dynamic compaction (DC) is widely used to improve the mechanical properties of soils and other granular fill material upon which foundations or other structures are to be built. To calculate the inner deformation induced by DC, a computational model based on stochastic medium theory was developed to deduce the amount of deformation in the fill from the geometry of the DC crater. For this model, the tamper-soil system was simplified to an axisymmetric geometry and the probability of any unit of material below the crater being deformed can be calculated. Subsequently, the deformation at any point can be obtained by integration. Input parameters for the model were established by back analysis using the results from a DC field test and published data. Comparing the model results to results from actual and simulated DC programs shows that the computational model is useful for calculating deformation induced by DC without the need to consider the constitutive model of the soil.

#### 1. Introduction

Compacting the soil that will be under a foundation is an essential procedure to ensure the stability of the structure. Improving the soil beneath the foundation can be accomplished by a number of different methods including static or dynamic rollers and fiber reinforcement [1–3]. However, static or dynamic rollers are limited and can only compact the soil to a relatively shallow depth (0.2–1.0 m) [4] and fiber reinforcement is expensive and not environment-friendly for large construction projects [5]. Compared with the above methods, dynamic compaction (DC) is a suitable method with the following advantages: (1) high density compaction to a considerable depth; (2) wide applicability in soils with different soil properties and different water contents; (3) simplicity of operation using widely available and simple equipment; (4) low cost. Menard [6] first introduced DC for construction and conducted systematically research on the method. Dynamic compaction simply involves dropping a tamper (10–40 t) from a height of 10–20 m onto the ground. Several blows are repeated until the ground under the tamper reaches the desired degree of soil improvement.

The extent of soil improvement is one of the most important concerns for a DC project; it is the essential indicator of the project’s success. The soil improvement can be divided into two parts: the depth of improvement (*DI*) and the lateral extent of improvement (*LI*).

The* DI* is a major issue facing the DC program designers and contractors [7]. The first equation for predicting* DI* was proposed by Menard [6] and it can be expressed as , where* M* is the tamper weight in tones and* H* is the falling height in meters. Mayne [8] summarized 120 sites with diverse soil properties and modified this equation to the form:where* n* (commonly equal to 0.3–0.8) is an empirical constant related to soil properties. Subsequent researches by Leonards [9], Lukas [10], Rollins [11], Slocombe [12], and Luoguo [13] calibrated and determined different values for* n* and proposed other forms for the equation based on the energy applied. However, these equations can only roughly evaluate the* DI *and do not consider factors like soil properties and tamper radius. Although researchers like Chow [14] and Smits [15] proposed computational models that did consider these factors, the parameters in their models are difficult to measure, and the values these factors have affect the results.

The* LI* is the maximum width the improvement extends horizontally away from the tamper. It is an essential parameter for determining the optimum distance between tamping points. There has been little research on how to calculate* LI*. Poran [16] discussed the factors that influence the degree of improvement and divided these factors into site-dependent factors and equipment-dependent factors. Two equations, one for the input energy and one for the dimensions of the plastic zone, were derived to evaluate the depth and width of improvement. The plastic zone was assumed to be semispherical shape.

In recent years, the dynamic mechanisms of DC have undergone further analysis and many simulations have been developed using finite or discrete element methods. For instance, Lee and Gu [17] proposed a method for estimating the degree and depth of improvement based on a two-dimensional finite element model. Ma [4] used the PFC2D/PFC3D software (Itasca Consulting Group, Inc., Minneapolis, MN, USA) to determine improvement and* DI* by analyzing porosity variations. Mostafa [18] conducted a comprehensive parametric study on the behavior of cohesive soils. Jia [19] investigated the behavior of granular soils during DC using PFC/FLAC software (Itasca Consulting Group). However, preparing these computational models and simulations is very complicated and time-consuming and it is not really feasible to use them in the field at construction sites.

During DC, the deformation can reflect the degree of soil improvement; this is a good indicator of how much the soil has been compacted. In most cases, when the soil in the crater is depressed less than a specified amounts per drop for the final two drops, it can be assumed that the soil’s strength has been improved enough to meet the project’s design requirements. The deformation in the soil can also be a reflection of the change in soil properties and used to determine how much the soil has been reinforced by the DC.

Stochastic medium theory was initially proposed by Litwiniszyn [20] and developed further by Liu et al. [21, 22]. The stochastic medium in question is granular material like sand, gravel, and rock fragments. The excavation of underground space often causes deformation of the overlying strata [23, 24] and the theory is mainly used to calculate subsidence after an opening has been excavated. The subsidence caused by the excavation will be transferred to the overlying strata (Figure 1(a)). For DC, the settling is caused by the impact of the tamper. Figure 1(b) shows soil particles compressed by tamping. Comparing Figure 1(a) with Figure 1(b) shows that the two kinds of settling are similar for noncohesive soils.

**(a) Mining subsidence**

**(b) Dynamic compaction**

In this paper, the deformation caused by the DC is modeled using stochastic medium theory. A computational model is proposed based on this theory and two parameters, a compression factor (determined by the type of DC equipment used) and an influence angle (a function of the soil’s properties). Finally, the applicability of this model is examined by comparing results calculated using the model with measurement made after a field test and with published data. This study can provide a method for determining how well a DC program is reinforcing the material upon which a foundation will be installed and will be valuable for the design and construction of DC projects.

#### 2. Model for Dynamic Compaction

##### 2.1. Introduction to Stochastic Medium Theory

In an Eulerian coordinate system, assume that the probability density for deformation of* dx* is for the* z*-axis symmetry of the deformation curve. The deformation probability* P* for* dx* can be expressed as

For an infinitesimal plane* ds *=* dxdy* centered at point* A*(*x*,*y*,*z*), as shown in Figure 2(a), the deformation probability is

**(a) Deformation of point**

**(b) The transformation of coordinate system**

If the coordinate system (*x*,*y*,*z*) is changed to (*x*_{1},*y*_{1},*z*) as shown in Figure 2(b), then (3) is transformed into

If the* x*_{1} axis passes through the point* A*, then

in which* c* is a constant.

The derivation of (7) will yield:

therefore

Equation (9) shows that the two sides of the equation are independent of* x* and* y*.

Let

Then, after integration, the following equation can be obtained:where* k*_{1} and* k*_{2} are parameters corresponding to the depth,* z*.

Therefore, the deformation probability of plane* ds* is

Assuming that the soil is not compacted, the deformation should be equal to the volume of the unit (*V*_{e}) in the crater, so

After solving, it can be seen that

Based on stochastic medium theory, the parameter can be expressed as* r*(*z*) (also called the influence radius) and* r*(*z*) is defined as where* z* is the depth of* ds* and *β* is the influence angle.

Therefore, the deformation at the point (*x*,* y*,* z*) induced by a unit in the crater can be expressed as

##### 2.2. Computational Model

Assuming that the foundation surface is smooth and the face of the tamper is circular, the crater formed by the tamper can be simplified as a frustum of a cone (Figure 3) and the DC can be simplified into an axisymmetric problem. To simplify the calculation, we only study the deformation in plane* zox* of the coordinate system. The deformation of point* A* induced by the unit in the crater can be divided into two directions: the vertical and the horizontal, which are* W*_{e} and in Figure 2(a), respectively. The vertical deformation is much easier to be detected than the horizontal. Therefore, in this paper, only vertical displacement is considered; horizontal displacement is ignored. Taking the central axis of the crater as the* z*-axis, for a unit* ds = dx*_{e}*dy*_{e}*dz*_{e} centered at (*x*_{e},* y*_{e},* z*_{e}) in the crater, the induced deformation at point* A*(*x*,* y*,* z*) (Figure 2) can be expressed by

The deformation induced in the crater can be obtained by integration, so

where Ω represents the whole region of the crater.

The crater’s coordinate system can be transformed into cylindrical coordinates and the unit in the crater changed from to . The geometric parameters of the crater (Figure 3) are the depth,* h*, and the depth within the (0,* h*) range represented by ; a vertical section through the crater is trapezoidal and the radius at different crater depths can be expressed as which can be calculated fromwhere and are the crater’s top and the bottom radii.

The deformation at any depth* z* can be calculated from

Dynamic compaction makes the soil denser. Therefore, a compression coefficient, , must be introduced to calculate the vertical movement of a unit of compacted material.

The influence angle can reflect the properties of the soil at depth* z*. However, if the soil at the construction site is considered to be homogeneous, the influence angle will be constant. For DC, the impact energy influence angle is the main factor affecting the amount of compaction. Therefore, the impact energy influence angle is used in the calculations in this paper. According to Yao [25], the DC impact influence angle for gravelly soil can be assumed to be 30° from the horizontal, so 30° is used as the influence angle in all the calculations presented below. The subsidence at a specific depth can be calculated using (20); the only unknown is , the compression coefficient.

#### 3. Model Verification

##### 3.1. Field Test Results

To acquire values for the variables in the model, a field test was conducted at the Chengde airport construction site, Hebei Province, China. The construction material was mainly crushed rock that complied with airport fill technical code MH/T5035-2017 [26]. Particle gradation curves for the eight different 1-m layers in the fill are shown in Figure 4.

The machine used for the DC was a 70-t crawler. The tamper weighed 16.8 t and was 2.5 m in diameter. The tamper imparted a DC energy of 1000 kNm with a fall of 6 m. To show the amount of settling in the different fill layers caused by the DC, thin plastic sheeting was laid down at three different horizons as the fill was being placed. The sheeting was installed at depths of 1.85 m, 2.8 m, and 4 m below the final fill surface as shown in Figures 5 and 6. For the calculations below, the geometry of the crater was simplified to make the crater a cylinder with a radius of 1.2 m and the parameters used in the calculations were = = 1.2 m,* h *= 0.8 m.

The most appropriate value of is the lowest least squares value for calculated using

where and are the calculated and observed results for deformation, respectively. The lower the , the better the match between the measured and the calculated result.

The observed and the calculated data are listed in Table 1. The calculated values are obviously larger than the observed values because the calculated values do not take the compression coefficient into account. The relationship between and was calculated using (21) and the results are shown in Figure 7. As shown in the figure, is minimized when = 0.65. The modified data listed in Table 1 were calculated with = 0.65. The final deformation is shown in Figure 8; the modified values agreed well with the observed values.

| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Note: Modified value = Calculated value × Compression Coefficient |

##### 3.2. Published Data (Ghassemi et al., [27])

The below-surface deformation induced by DC is rarely reported in field test results. Therefore, the results from this study are compared with the results from a DC simulation published by Ghassemi et al. [27]. Ghassemi et al. used finite element software to simulate sand compaction and verified the results of Oshima and Takada’s [28] centrifuge tests. The simulated DC energy was 4000 kNm, energy equivalent to a 20-t tamper dropping 20 m. Figure 1 in [27] shows the deformation after DC. The dimensions of the simulated crater in Ghassemi et al. were = 1.4 m, = 1.2 m,* h *= 1.2 m. Those crater dimensions were entered into (20) to calculate the soil deformation for that tamper configuration; the results are listed in Table 2. Figure 9 shows that for this DC simulated configuration is 0.7, so that value was used to calculate the “Modified” results listed in Table 2. As shown in Figure 10, the modified results show a good correlation with the data from Ghassemi et al.

| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Note: Modified value = Calculated value × Compression Coefficient |

#### 4. Discussion

##### 4.1. Critical Deformation

The computational model has shown its applicability for DC analysis. However, the primary shortcoming of this model is that the calculated deformation cannot be zero even for locations far removed from the tamped crater. This means that the zero deformation boundary cannot be determined using this model. To compensate for this shortcoming, the concept of a critical level of deformation has been formulated. The critical level of deformation is defined as deformation sufficient to significantly reinforce the soil. Deformation below this level will just modify the soil but not increase its strength. By comparing the Standard penetration test* N*_{120} (SPT) values from the fill at the Chengde airport construction site before and after DC (Figure 11), the rate of improvement can be calculated. Equation (22) can be used to quantify the amount of improvement.where* ω* is the rate of improvement.

Figure 12 is a graph showing the rate of improvement for a number of points in the fill plotted against the corresponding inner deformation. The fitted equation for the curve on that figure is

According to SPT’s performed at the construction site, the* DI* where the rate of improvement is 50% is about 5 m. The vertical deformation of points at the 5-m depth, from (22), is 0.04 m. Therefore, 0.04 m was determined to be the amount of deformation that defines the critical level of deformation boundary. Points within that boundary should have had their strength improved by at least 50%. The boundary’s location is defined by points with a vertical deformation of 0.04 m; a line drawn thorough those points is shown in Figure 13. The line shows that the reinforced area is bubble-shaped and the reinforced area extends nearly 3 m from the center of the crater. The area of reinforcement is consisted with the field test results.

##### 4.2. Crater Depth and the Reinforcement Extent

The* DI *and* LI *away from the tamping sites are two essential parameters that must be considered when a DC program is being designed. These parameters are commonly determined by using the Menard equation or practical field experience; both are largely based on the energy produced by the DC. However, for soils with different properties, the same impact energy produces different degrees of reinforcement. Because it is difficult to observe the change that have taken place in the soil after DC, the tamper-soil system is somewhat of a black box in that we only know the inputs to the system but cannot easily monitor the results. However, an “external interface” for this “black box” exists on the geometry of the crater. For the same DC energy, the shapes and the depths of craters are different for different soils. The crater is a reflection of the soil-tamper system and can be used to analyze the posttamping changes in the soil. In our stochastic-medium-based computational model, the crater is used as the main indicator of soil reinforcement. Figure 14 is the deformation of different depth with different crater depth calculated by the model. The deformation is decreased with the increase of the depth and radial distance. Meanwhile, the increase of the crater depth is corresponding to a larger reinforcement extent.

**(a) Vertical direction**

**(b) Horizontal direction**

Figures 15(a) and 15(b) present the relationships between* DI* and crater depth,* LI* and crater depth, respectively. The* DI* increases with the growth of the crater depth and the correlation is very nearly linear. The* LI* shows the same character with the* DI*, but the trend of growth slows down with the increase of crater depth. When the crater depth increases from 1 m to 2 m, the growth of* LI* is about 1.2 m. However, the growth of* LI* is just 0.6 m when the crater depth increases from 2 m to 3 m. Therefore, it can be concluded that, increasing the crater depth, the* DI* can be significantly expanded while* LI* increases consistently with a slowly decreasing rate (Figure 15(b)).

**(a) DI**

**(b) LI**

#### 5. Conclusions

A computational model based on stochastic medium theory has been developed to determine how much dynamic compaction (DC) reinforces the soil. This model is verified by comparing model results with the results from two compaction programs, one on crushed rock fill and the other simulating a sand fill. The compression coefficients for these two programs are 0.65 for the crushed rock and 0.7 for the sand. This method for establishing the degree of reinforcement has several advantages when compared with traditional methods. First, it takes the geometry of the crater into account. The shape and depth of the crater can reflect the DC energies and the soil properties and the crater’s dimension are easy to obtain at the construction site with no special equipment or personnel training required. Second, the model assesses the extent of the reinforcement both laterally and vertically instead of merely estimating the* DI*. Finally, the computation is simple and can be easily solved using commercially available programs. However, this model is a prototype and more details should be studied. Additional validation of the computational model is also needed before sufficient confidence on the model’s results can be established.

#### Data Availability

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

#### Conflicts of Interest

The authors declare no conflicts of interest.

#### Acknowledgments

This research was supported by the National Basic Research Program of China, grant number 2014CB047004. We thank David Frishman, Ph.D0, from Liwen Bianji, Edanz Group, China, for editing the English text of a draft of this manuscript.

#### References

- Y. Wang, P. Guo, F. Dai, X. Li, Y. Zhao, and Y. Liu, “Behavior and modeling of fiber-reinforced clay under triaxial compression by combining the superposition method with the energy-based homogenization technique,”
*International Journal of Geomechanics*, vol. 18, no. 12, Article ID 04018172, 2018. View at: Publisher Site | Google Scholar - Y. Wang, P. Guo, X. Li, H. Lin, Y. Liu, and H. Yuan, “Behavior of fiber-reinforced and lime-stabilized clayey soil in triaxial tests,”
*Applied Sciences*, vol. 9, no. 5, Article ID 900, 2019. View at: Publisher Site | Google Scholar - Y. X. Wang, P. P. Guo, H. Lin et al., “Numerical analysis of fiber-reinforced soils based on the equivalent additional stress concept,”
*International Journal of Geomechanics*, 2019. View at: Google Scholar - Z. Ma, F. Dang, and H. Liao, “Numerical study of the dynamic compaction of gravel soil ground using the discrete element method,”
*Granular Matter*, vol. 16, no. 6, pp. 881–889, 2014. View at: Publisher Site | Google Scholar - S.-J. Feng, K. Tan, W.-H. Shui, and Y. Zhang, “Densification of desert sands by high energy dynamic compaction,”
*Engineering Geology*, vol. 157, pp. 48–54, 2013. View at: Publisher Site | Google Scholar - L. Menard and Y. Broise, “Theoretical and practical aspects of dynamic consolidation,”
*Géotechnique*, vol. 25, no. 1, pp. 3–18, 1975. View at: Publisher Site | Google Scholar - W.-L. Zou, Z. Wang, and Z.-F. Yao, “Effect of dynamic compaction on placement of high-road embankment,”
*Journal of Performance of Constructed Facilities*, vol. 19, no. 4, pp. 316–323, 2005. View at: Publisher Site | Google Scholar - P. W. Mayne, J. S. Jones Jr., and J. C. Dumas, “Ground response to dynamic compaction,”
*Journal of Geotechnical Engineering*, vol. 110, no. 6, pp. 757–774, 1984. View at: Publisher Site | Google Scholar - G. A. Leonards, W. A. Cutter, and R. D. Holtz, “Dynamic compaction of granular soils,”
*Journal of the Geotechnical Engineering Division*, vol. 106, no. 1, pp. 35–44, 1980. View at: Google Scholar - R. G. Lukas, “Dynamic compaction for highway construction, Design and construction guidelines, Volume I,”
*Federal Highway Administration Report FHWA-RD-86-133*, 1986. View at: Google Scholar - K. M. Rollins and J. Kim, “Dynamic compaction of collapsible soils based on U.S. case histories,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 136, no. 9, pp. 1178–1186, 2010. View at: Publisher Site | Google Scholar - B. C. Slocombe,
*Dynamic Compaction. Ground Improvement*, vol. 2, Glasgow, Chapman & Hall. Cap., 1993. - V. Luongo, “Dynamic compaction: predicting depth of improvement,”
*Grouting, Soil Improvement and Geosynthetics, ASCE*, pp. 27–939, 1992. View at: Google Scholar - Y. K. Chow, D. M. Yong, K. Y. Yong, and S. L. Lee, “Dynamic compaction of loose sand deposits,”
*Soils and Foundations*, vol. 32, no. 4, pp. 93–106, 1992. View at: Publisher Site | Google Scholar - M. T. Smits and L. De Quelerij, “The effect of dynamic compaction on dry granular soils,” in
*Proceedings of the International Conference on Soil Mechanics and Foundation Engineering, 12th*, vol. 2, Rio de Janiero, Brazil, 1989. View at: Google Scholar - C. J. Poran and J. A. Rodriguez, “Finite element analysis of impact behavior of sand,”
*Soils and Foundations*, vol. 32, no. 4, pp. 68–80, 1992. View at: Publisher Site | Google Scholar - F. H. Lee and Q. Gu, “Method for estimating dynamic compaction effect on sand,”
*Journal of Geotechnical and Geoenvironmental Engineering*, vol. 130, no. 2, pp. 139–152, 2004. View at: Publisher Site | Google Scholar - K. F. Mostafa and R. Y. Liang, “Numerical modeling of dynamic compaction in cohesive soils,” in
*Proceedings of the Geo-Frontiers Congress 2011*, pp. 738–747, Dallas, Texas, United States, 2011. View at: Publisher Site | Google Scholar - M. Jia, Y. Yang, B. Liu, and S. Wu, “PFC/FLAC coupled simulation of dynamic compaction in granular soils,”
*Granular Matter*, vol. 32, no. 4, Article ID 76, 2018. View at: Publisher Site | Google Scholar - J. Litwiniszyn, “The theories and model research of movements of ground masses,” in
*Proceedings of the European Congress Ground Movement*, pp. 203–209, Leeds, UK, 1957. View at: Google Scholar - B. C. Liu, “Ground surface movements due to underground excavation in the People's Republic of China,”
*Excavation, Support and Monitoring*, pp. 781–817, 1993. View at: Google Scholar - J. S. Yang, B. C. Liu, and M. C. Wang, “Modeling of tunneling-induced ground surface movements using stochastic medium theory,”
*Tunnelling and Underground Space Technology*, vol. 19, no. 2, pp. 113–123, 2004. View at: Publisher Site | Google Scholar - P. Guo, X. Gong, and Y. Wang, “Displacement and force analyses of braced structure of deep excavation considering unsymmetrical surcharge effect,”
*Computers & Geosciences*, vol. 113, Article ID 103102, 2019. View at: Publisher Site | Google Scholar - Y. Wang, S. Shan, C. Zhang, and P. Guo, “Seismic response of tunnel lining structure in a thick expansive soil stratum,”
*Tunnelling and Underground Space Technology*, vol. 88, pp. 250–259, 2019. View at: Publisher Site | Google Scholar - Y. P. Yao and B. Z. Zhang, “Reinforcement range of dynamic compaction based on volumetric strain,”
*Rock and Soil Mechanics*, vol. 9, Article ID 032, 2016 (Chinese). View at: Google Scholar - MH/T5035-2017, “Technical code for high filling engineering of airport,”
*Civil Aviation Administration of China*, 2017 (Chinese). View at: Google Scholar - A. Ghassemi, A. Pak, and H. Shahir, “A numerical tool for design of dynamic compaction treatment in dry and moist sands,”
*Iranian Journal of Science & Technology, Transaction B: Engineering*, vol. 33, no. 4, pp. 313–326, 2009. View at: Google Scholar - A. Oshima and N. Takada, “Relation between compacted area and ram momentum by heavy tamping,” in
*Proceedings of the International Conference on Soil Mechanics and Foundation Engineering-International Society for Soil Mechanics and Foundation Engineering, AA BALKEMA*, vol. 3, pp. 1641–1644, 1997. View at: Google Scholar

#### Copyright

Copyright © 2019 Jifang Du et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.