Research Article | Open Access

Ming Sun, Meng-fan Zong, Shao-jun Ma, Wen-bing Wu, Rong-zhu Liang, "Analytical Solution for One-Dimensional Consolidation of Soil with Exponentially Time-Growing Drainage Boundary under a Ramp Load", *Mathematical Problems in Engineering*, vol. 2018, Article ID 9385615, 10 pages, 2018. https://doi.org/10.1155/2018/9385615

# Analytical Solution for One-Dimensional Consolidation of Soil with Exponentially Time-Growing Drainage Boundary under a Ramp Load

**Academic Editor:**Enrico Conte

#### Abstract

By introducing the exponentially time-growing drainage boundary, this paper investigated the one-dimensional consolidation problem of soil under a ramp load. Firstly, the one-dimensional consolidation equations of soil are established when there is a ramp load acting on the soil surface. Then, the analytical solution of excess pore water pressure and consolidation degree is derived by means of the method of separation of variables and the integral transform technique. The rationality of this solution is also verified by comparing it with other existing analytical solutions. Finally, the consolidation behavior of soil is studied in detail for different interface parameters or loading scheme. The results show that the exponentially time-growing drainage boundary can reflect the phenomenon that the excess pore water pressure at the drainage boundaries dissipates smoothly rather than abruptly from its initial value to the value of zero. By adjusting the values of interface parameters and* c*, the presented solution can be degraded to Schiffman’s solution, which can compensate for the shortcoming that Terzaghi’s drainage boundary can only consider the two extreme cases of fully pervious and impervious boundaries. The significant advantage of the exponentially time-growing drainage boundary is that it can be applied to describe the asymmetric drainage characteristics of the top and bottom drainage surfaces of the actual soil layer by choosing the appropriate interface parameters and* c*.

#### 1. Introduction

The consolidation problem of soil has always been one of the key problems of soil mechanics, for it has close relation with the deformation, strength, stability, and seepage of soil mechanics. Since Terzaghi [1] published his consolidation theory in 1925, noteworthy process has been achieved in consolidation theory. In order to accurately describe the soil properties, a considerable amount of scholars has made a great effort to build the governing equation of consolidation problem for soil according to different soil constitutive models, such as Terzaghi-Rendulic diffusion equation [2], three-dimensional Biot’s consolidation equation [3], Maxwell model [4], Kelvin model [5, 6], and Merchant model [7]. Considering the significant influence of initial condition on the consolidation process of soil, a large number of investigators have also proposed various solutions for one-dimensional consolidation problem of soil under different loads, such as linear load [8], rectangular load [9], random load [10], cyclic loading [11–13], and ramp loading [13–15].

From the aforementioned literature review, it can be found that the drainage boundaries are normally treated as pervious or impervious, that is, Terzaghi’s drainage boundary. Nevertheless, the drainage at the boundaries of consolidating soil is impeded for most engineering applications [16]. With regard to this aspect, Gray [16] first proposed an impeded drainage boundary which was later pursued by Schiffman and Stein [17] to investigate the consolidation problem of soil with variable permeability and compressibility. Following Gray’s work, solutions of one-dimensional consolidation have been presented for two- [18] and multilayered soils [19, 20] with impeded drainage boundaries. However, Lee et al. [21] pointed out that the impeded drainage boundary conditions can be regarded as special cases of permeable, yet incompressible, top or bottom soil layer in multilayered system. Recently, Mei et al. [22] have put forward the exponentially time-growing drainage boundary condition (i.e., continuous drainage boundary as Mei’s definition), including permeable and impermeable boundary conditions as two extremities, to consistently describe the all drainage boundary conditions of saturated soil. Soon afterwards, studies on consolidation theory with exponentially time-growing drainage boundary condition have been conducted for saturated soils with a single layer [23] and multilayers [24] and unsaturated soil with a single layer [25]. Since the exponentially time-growing drainage boundary can allow the excess pore water pressure to dissipate smoothly rather than abruptly from its initial value given by the initial conditions to the value of zero, it may have board application prospect in engineering.

The objective of this article is to develop a general analytical solution for the one-dimensional consolidation of soil with exponentially time-growing drainage boundary (abbreviated as ETGD boundary for convenience) under a ramp load. Moreover, the detailed solution is obtained for the excess pore water pressure and overall average consolidation degree in terms of excess pore water pressure. Then, a comparison is made between the present solution and Schiffman’s solution [8] to verify the present solution. Compared with Conte’s solution [13] and Olson’s solution [14], it is also worth noting that the proposed solution draws closer to engineering reality for it can account for the variation of drainage boundaries with time during the consolidation process. In addition, the influence of the interface parameters of the ETGD boundary on the consolidation behavior of soil is discussed in detail based on the present solution.

#### 2. Problem Description

The thickness, coefficient of permeability, and coefficient of compressibility of soil are denoted as , , and , respectively. is half the thickness of the soil. As shown in Figure 1, a ramp load is applied on the top surface of soil layer, and the ultimate load and the loading time of the ramp load are represented as and , respectively. For a single stage load with constant loading rate, it can be expressed aswhere indicates the loading rate. For a constant ultimate load, the smaller means the larger loading rate.

According to Terzaghi’s consolidation theory, the governing equation of one-dimensional consolidation problem of soil can be expressed aswhere represents the excess pore water pressure of soil; is the downward vertical coordinate originated from the top boundary; is the loading time. indicates the consolidation coefficient of soil layer and denotes the unit weight of water, respectively.

Then, this work introduces the ETGD boundary as follows [20]:where and are interface parameters which reflect the drainage properties of the top and bottom drainage surfaces of soil layer, respectively. The two parameters can be obtained by experimental simulation or engineering measurement inversion.

The initial boundary is given as

#### 3. Analytical Solution for This Model

Letting and substituting it into (2) yieldwhere

Combining with the expression form of ramp load , can be rewritten aswhere

Then, the boundary conditions and initial condition can be rewritten as

The general solution of (5) can be set as

Substituting (12) into (5) giveswhere can be expressed according to Fourier series aswhere

In order to make (13) always true, the following equation can be derived:

For a single stage load with constant loading rate, it can be determined that . Moreover, it can also obtain that because .

Utilizing the Laplace transform technique to solve (17) yieldswhere

##### 3.1. Solution of Loading Stage

For loading stage (i.e., ), one can derive

Substituting (20) into (12) gives

Furthermore, the excess pore water pressure of soil during the loading stage can be derived as

##### 3.2. Solution of Constant Load Stage

For constant load stage (i.e., ), one can obtainwhere

Then, the excess pore water pressure of soil during the constant load stage can be obtained as

#### 4. Verification and Parametric Study

##### 4.1. Degeneration of Solution

In order to verify the effectiveness and accuracy of the presented solution, it is needed to degenerate the proposed solution and compare it with existing solutions. Therefore, this section conducts the verification for both loading stage and constant load stage.

###### 4.1.1. Degeneration of Solution at Loading Stage

When the interface parameters and approach infinity, (22) can be degenerated as

###### 4.1.2. Degeneration of Solution at Constant Load Stage

When the interface parameters and tend to infinity, (24) can be reduced as

Then, (25) can be degenerated as

Apparently, (26) and (28) are consistent with the corresponding equations proposed by Schiffman utilizing Terzaghi’s drainage boundary [8].

Then, a relevant computer program is also developed and applied to compare the presented solution with Schiffman’s solution [8]. For the degeneration of the ETGD boundary, the interface parameters, both and are set as sufficiently small values (i.e., ) and large enough values (i.e., ) which can be regarded as the impervious boundary and pervious boundary, respectively. The unit weight, elastic modulus, Poisson’s ratio, and permeability coefficient of soil are , , , and , respectively. Lastly, the time factor, which is expressed as , is also introduced for convenient comparison, and is the corresponding time factor with respect to the final loading time.

As shown in Figure 2, when the interface parameters and* c* and the time factor have smaller values, tends to 0.9 because the drainage velocity is small and the drainage time is very short. For the curve obtained by Schiffman’s solution, there are mutations at and which are different with the actual situations. When the interface parameters and* c* and the time factor have larger values, both the values of obtained by the presented solution and Schiffman’s solution are consistent and in symmetrical distribution along depth. The values of tend to zero because the drainage velocity is large enough and the drainage time is long enough. The reason for that the values of are not zero is that the consolidation process is at loading stage and the increment of applied loading should be first borne by the excess pore water pressure. The results show that the solution obtained by Terzaghi’s drainage boundary can not be reduced to the solution degenerated from the presented solution, and the curves obtained by the solution based on Terzaghi’s drainage boundary are not smooth and have mutations at the top and bottom drainage surfaces. In other words, whatever the value of time factor is, the excess pore water pressure based on Terzaghi’s drainage boundary at and is always zero. Actually, the pore water pressure can not be completely dissipated when the drainage time is zero or very small, and the exponentially time-growing drainage boundary can reflect this phenomenon that the excess pore water pressure at the drainage boundary dissipates smoothly rather than abruptly from its initial value given by the initial conditions to the value of zero.

Moreover, further comparison of the proposed solution with Schiffman’s solution is conducted by analyzing the influence of on excess pore water pressure for both loading stage and constant load stage. The smaller means the shorter loading time from loading stage to constant load stage and also implies the faster loading rate. As can be observed in Figure 3, at the loading stage (<), the excess pore water pressure increases with the increasing loading rate for the external load increases with the increase of loading rate at a given time. However, at the constant load stage ( > ) as shown in Figure 4, the excess pore water pressure at a given time decreases as the loading rate increases, for there is a longer time to dissipate the excess pore water pressure if the loading process is completed within a short period of time. In addition, it can also be seen from the diagrams that when the interface parameters and equal to 1000d^{−1}, the pore pressure curve obtained by the ETGD boundary coincides with the pore pressure curve obtained by Schiffman’s solution. This indicates that the solution obtained by the ETGD boundary can be reduced to Schiffman’s solution when and value are large enough, which further verifies the correctness of the present solution.

##### 4.2. Influence of Number of Term in Series on Calculated Results

Figure 5 indicates the influence of number of terms in series on calculated results. For the parameters of , it can be seen that the excess pore water pressure of loading stage is larger than that of constant load stage. The excess pore water pressure increases as the time factor increases during the loading stage. But for the constant load stage, the excess pore water pressure decreases with the increase of time factor because the pore water pressure dissipates gradually. When the time factor has smaller value, it needs more number of terms in series to meet the convergence of solution. However, if the time factor is large enough, the fluctuation of the curves is small and the convergence of solution can be satisfied by only taking two or three terms. Therefore, unless otherwise specified, the number of terms in series is set as five terms for efficient calculation in the following analysis.

##### 4.3. Influence of Interface Parameters and on the Excess Pore Water Pressure of Soil

Figure 6 illustrates the excess pore water pressure distribution curves at loading stage while* b*=*c*. When is equal to* c*, the top and bottom drainage surfaces of soil exhibit the same permeability which results in the fact that the excess pore water pressure curves are in symmetrical distribution about . With the gradual increase of time factor , the excess pore water pressures at the drainage surfaces and interior of soil gradually dissipate, instead of that the excess pore water pressures obtained by Schiffman’s solution are always zero at the drainage surfaces. It can be also observed from Figure 6 that the excess pore water pressure obtained by the presented solution will increase first and then decrease at the drainage surfaces. The reason for the abovementioned phenomena may be that at the beginning the excess pore water pressure at the drainage surfaces increases with the increase of the applied loading; after that the rate of dissipation becomes greater than the rate of loading increment. It is also worth noting that the solution based on Terzaghi’s drainage boundary can not reflect the above-mentioned phenomena.

As depicted in Figure 7, when is not equal to* c*, for example, , the upper part of the excess pore water pressure distribution curves at loading stage is obviously steeper than the lower part which means that the drainage velocity of the upper part of the soil is larger than that of the lower part of soil. Therefore, it can also be understood that the greater the value of the interface parameter, the larger the drainage velocity of the corresponding haft part of the soil. Terzaghi’s drainage boundary can not consider the difference of the drainage capacity between the upper and lower parts of soil layer, for the excess pore water pressure distribution curves at loading stage based on Terzaghi’s drainage boundary are always symmetrical about [8]. It can also be seen from Figure 7 that the significant advantage of the ETGD boundary is that it can be applied to approximately model the asymmetric drainage characteristics of the top and bottom drainage surfaces of actual soil layer by adjusting the interface parameters and* c*.

As demonstrated in Figure 8, when is equal to* c*, the top and bottom drainage surfaces of soil also exhibit the same permeability for the excess pore water pressure curves at constant load stage also distribute symmetrically about . It is also found that the excess pore water pressure at the drainage surfaces and interior of soil gradually dissipates as the time factor increases, but the excess pore water pressures at constant load stage obtained by Schiffman’s solution are also zero at the drainage surfaces.

As shown in Figure 9, when is not equal to* c*, for example , the lower part of the excess pore water pressure distribution curves at constant load stage is obviously steeper than the upper part which implies that the drainage velocity of the lower part of the soil is larger than that of the upper part of the soil. Meanwhile, the solution based on Terzaghi’s single-drainage condition can not be extended to the soil layer with the upper and bottom surfaces having different drainage capacity, for the excess pore water pressure distribution curves at constant load stage based on Terzaghi’s drainage boundary are always symmetrical about .

##### 4.4. Influence of Consolidation Coefficient on the Excess Pore Water Pressure of Soil

Figure 10 shows the influence of consolidation coefficient on the excess pore water pressure of soil. For the loading stage, for example, , the excess pore water pressure of soil decreases as the consolidation coefficient decreases. For the constant load stage, for example, , the excess pore water pressure of soil also decreases with the decrease of consolidation coefficient. The reason for these phenomena is that, for the same time factor , the smaller consolidation coefficient means the longer consolidation time, and the longer consolidation time will lead to more sufficient dissipation of the excess pore water pressure. It is also worth noting that the smaller consolidation coefficient means the weaker permeability of soil which results in smaller drainage velocity of soil. Therefore, the results show that the consolidation time has more significant influence on the excess pore water pressure than consolidation coefficient. As can be seen from Schiffman’s results [8], the change of consolidation coefficient for constant time factor has very weak impact on the excess pore water pressure which means that Terzaghi’s drainage boundary can not completely reflect the actual mechanism of consolidation process of soil influenced by time. It can also be observed from Figure 10 that the degree of the influence of consolidation coefficient on excess pore water pressure decreases with the decrease of consolidation coefficient.

##### 4.5. Influence of Loading Rate on the Excess Pore Water Pressure of Soil

Figure 11 indicates the influence of loading rate on excess pore water pressure when the other parameters of soil remain constant. As can be observed in Figure 11, the excess pore water pressure increases with the increase of loading rate because the excess pore water pressure can not be sufficiently dissipated within the shorter consolidation time when the loading rate is too fast. For the constant load stage, the excess pore water pressure under the same time factor decreases as the loading rate increases, for there is a longer time to dissipate the excess pore water pressure if the loading process is completed within a short period of time.

#### 5. Analysis of Consolidation Degree

According to Soil Mechanics, the overall average consolidation degree in terms of excess pore water pressure can be expressed as

Then, by substituting (22) into (29), the overall average consolidation degree in terms of excess pore water pressure at loading stage can be obtained as

Furthermore, by substituting (25) into (29), the overall average consolidation degree defined by excess pore water pressure at constant load stage can be derived as

Figure 12 depicts the influence of interface parameters and on the consolidation degree of soil. The consolidation degree obtained by the presented solution will approach that calculated by Schiffman’s solution if the interface parameters and are large enough, for the boundary surfaces with increasing interface parameters will tend to be pervious boundary. In contrast, if the interface parameters and are becoming smaller, the consolidation degrees obtained by the proposed solution will gradually deviate from that calculated by Schiffman’s solution. It can also be found that the difference between the consolidation degree obtained by the proposed solution and that calculated by Schiffman’s solution increases with the decrease of time factor. It is also worth noting from Figure 12 that the early consolidation process of Terzaghi’s drainage boundary is faster than that of the ETGD boundary, but the late consolidation process of the ETGD boundary is faster than that of Terzaghi’s drainage boundary. Furthermore, if the interface parameters and are large enough, the consolidation curves obtained by the ETGD boundary will essentially coincide with those obtained by Terzaghi’s drainage boundary at the late consolidation process. The results show that, for practical engineering, even though the early consolidation process of the soil with the boundary surfaces designed according to the ETGD boundary is slow, it can make more sufficient dissipation of the excess pore water pressure of the whole soil layer, and it will not lead to the slow drainage of the lower soil layer or even the discharge of water after the consolidation of the upper soil layer. At the same time, the ultimate consolidation time based on the ETGD boundary is almost the same as that based on Terzaghi’s drainage boundary. When utilizing the ETGD boundary to design the actual drainage boundaries, there is no need to deal with the drainage surfaces as fully pervious boundaries. It means that the overall economy of the design based on the ETGD boundary is clearly superior to that based on Terzaghi’s drainage boundary.

As shown in Figure 13, the smaller initial time factor implies the faster loading process of ramp load. It can be seen that the faster the loading process of ramp load, the faster the consolidation process. When the loading rate has larger value, the early consolidation process is faster but the late consolidation process is slower. In contrast, when the loading rate has slower value, the early consolidation process is slower but the late consolidation process is faster. For the ETGD boundary, when the consolidation degree is equal to 80%, the time factor is about 1.1, 1.3, and 8.8 when is equal to 0.1, 1, and 10, respectively. For Terzaghi’s drainage boundary, when consolidation degree is equal to 80%, the time factor is about 0.63, 1.2, and 8.4 when is equal to 0.1, 1, and 10, respectively. The results imply that the ultimate consolidation time based on the ETGD boundary is almost the same as that based on Terzaghi’s drainage boundary when the loading rate is slow enough.

#### 6. Conclusions

In this article, an analytical solution for the one-dimensional consolidation of soil layer under a ramp load is derived by adopting the ETGD boundary. The proposed analytical solution is also verified by comparing it with Schiffman’s solution [8]. Then, selective results are presented to illustrate the consolidation behavior of soil with ETGD boundary. Schiffman’s solution can be regarded as a special case of the present solution, and the excess pore water pressure and consolidation degree calculated by the presented solution will approach those calculated by Schiffman’s solution if the interface parameters and are large enough. The important advantage of the ETGD boundary is that it can be utilized to describe the asymmetric drainage characteristics of the top and bottom drainage surfaces of the actual soil layer by choosing the appropriate interface parameters and . Significantly, the present solution is just for homogeneous soil, and the consolidation problem of multilayered soil will be investigated by the authors’ team in near future.

#### 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 that they have no conflicts of interest.

#### Acknowledgments

This research is supported by the National Natural Science Foundation of China (Grants nos. 51578164, 51678547, 51878634, 51878185, and 41807262) and the China Postdoctoral Science Foundation Funded Project (Grants no. 2016M600711 and no. 2017T100664). The Research Funds provided by MOE Engineering Research Center of Rock-Soil Drilling & Excavation and Protection (Grant no. 201402) and the Fundamental Research Funds for the Central Universities-Cradle Plan for 2015 (Grant no. CUGL150411) are also acknowledged.

#### References

- K. Terzaghi,
*Erdbaumechanik and bodenphysikalischer grundlage*, Lpz. Deuticke, 1925. - L. Rendulic,
*Porenziffer und porenwasserdruck in Tonen*, vol. 7, Bauingenieur, 1936. - M. A. Biot, “General theory of three-dimensional consolidation,”
*Journal of Applied Physics*, vol. 12, no. 2, pp. 155–164, 1941. View at: Publisher Site | Google Scholar - T. K. Tan, “Secondary time effects and consolidation of clays,”
*Science in China Series A-Mathematics*, vol. 1, no. 11, pp. 1060–1075, 1958. View at: Google Scholar - R. E. Gibson and K. Y. Lo,
*A theory of consolidation for soils exhibiting secondary compression*, Norwegian Geotechnical Institute, Oslo, Norway, 1961. - K. Y. Lo, “Secondary compression of clays,”
*Journal of the Soil Mechanics and Foundation Division, ASCE*, vol. 87, no. 4, pp. 61–88, 1961. View at: Google Scholar - K. H. Xie and X. W. Liu, “A study on one dimensional consolidation of soils exhibiting rheological characteristics,” in
*International Symposium on Compression and Consolidation of Clayey Soils*, pp. 385–388, Rotterdam, The Netherlands, 1995. View at: Google Scholar - R. L. Schiffman, “Consolidation of soil under time-dependent loading and varying permeability,” in
*Proceedings of the Highway Research Board*, vol. 37, pp. 584–617, 1958. View at: Google Scholar - N. E. Wilson and M. M. Elgohary, “Consolidation of Soils Under Cyclic Loading,”
*Canadian Geotechnical Journal*, vol. 11, no. 3, pp. 420–423, 1974. View at: Publisher Site | Google Scholar - E. E. Alonso and R. J. Krizek, “Randomness of settlement rate under stochastic load,”
*Journal of the Geotechnical Engineering Divisio*, pp. 1211–1226, 1974. View at: Google Scholar - M. Favaretti and M. Soranzo, “A simplified consolidation theory in cyclic loading conditions,” in
*International Symposium on Compression and Consolidation of Clayey Soils*, A. A. Balkema, Ed., pp. 405–409, Rotterdam, The Netherlands, 1995. View at: Google Scholar - J. Z. Chen, X. R. Zhu, K. H. Xie, Q. Pan, and G. X. Zeng, “One dimensional consolidation of soft clay under trapezoidal cyclic loading,” in
*Proceedings of 2nd International Conference on Soft Soil Engineering*, pp. 211–216, Hohai University Press, Nanjing, 1996. View at: Google Scholar - E. Conte and A. Troncone, “One-dimensional consolidation under general time-dependent loading,”
*Canadian Geotechnical Journal*, vol. 43, no. 11, pp. 1107–1116, 2006. View at: Publisher Site | Google Scholar - R. E. Olson, “Consolidation under time-dependent loading,”
*Journal of the Geotechnical Engineering Division*, vol. 103, no. 1, pp. 55–60, 1977. View at: Google Scholar - E. Conte and A. Troncone, “Nonlinear consolidation of thin layers subjected to time-dependent loading,”
*Canadian Geotechnical Journal*, vol. 44, no. 6, pp. 717–725, 2007. View at: Publisher Site | Google Scholar - H. Gray, “Simultaneous consolidation of contiguous layers of unlike compressive soils,”
*Transactions of the American Society of Civil Engineers*, vol. 110, pp. 1327–1356, 1945. View at: Google Scholar - R. L. Schiffmann and J. R. Stein, “One-dimensional consolidation of layered systems,”
*Journal of the Soil Mechanics & Foundations Division*, vol. 96, no. 4, pp. 1499–1504, 1970. View at: Google Scholar - K.-H. Xie, X.-Y. Xie, and X. Gao, “Theory of one-dimensional consolidation of two-layered soil with partially drained boundaries,”
*Computers & Geosciences*, vol. 24, no. 4, pp. 265–278, 1999. View at: Publisher Site | Google Scholar - Y.-Q. Cai, X. Liang, and S.-M. Wu, “One-dimensional consolidation of layered soils with impeded boundaries under time-dependent loadings,”
*Applied Mathematics and Mechanics-English Edition*, vol. 25, no. 8, pp. 937–944, 2004. View at: Publisher Site | Google Scholar - R. P. Chen, W. H. Zhou, H. Z. Wang, and Y. M. Chen, “One-dimensional nonlinear consolidation of multi-layered soil by differential quadrature method,”
*Computers & Geosciences*, vol. 32, no. 5, pp. 358–369, 2005. View at: Publisher Site | Google Scholar - P. K. K. Lee, K. H. Xie, and Y. K. Cheung, “A study on one-dimensional consolidation of layered systems,”
*International Journal for Numerical & Analytical Methods in Geomechanics*, vol. 16, no. 11, pp. 815–831, 1992. View at: Google Scholar - G.-X. Mei, J. Xia, and L. Mei, “Terzaghi's one-dimensional consolidation equation and its solution based on asymmetric continuous drainage boundary,”
*Yantu Gongcheng Xuebao/Chinese Journal of Geotechnical Engineering*, vol. 33, no. 1, pp. 28–31, 2011. View at: Google Scholar - G.-X. Mei and Q.-M. Chen, “Solution of Terzaghi one-dimensional consolidation equation with general boundary conditions,”
*Journal of Central South University*, vol. 20, no. 8, pp. 2239–2244, 2013. View at: Publisher Site | Google Scholar - J.-C. Liu and G. H. Lei, “One-dimensional consolidation of layered soils with exponentially time-growing drainage boundaries,”
*Computers & Geosciences*, vol. 54, pp. 202–209, 2013. View at: Publisher Site | Google Scholar - L. Wang, D. Sun, and A. Qin, “Semi-analytical solution to one-dimensional consolidation for unsaturated soils with exponentially time-growing drainage boundary conditions,”
*International Journal of Geomechanics*, vol. 18, no. 2, 2018. View at: Google Scholar

#### Copyright

Copyright © 2018 Ming Sun 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.