Research Article | Open Access
Fang Wei, Wang Jiao, Zhang Wan-dongxing, Shi Li-jun, "Lower-Bound Solution of Foundation-Bearing Capacity under Circular Uniformly Distributed Load", Advances in Civil Engineering, vol. 2020, Article ID 8890450, 6 pages, 2020. https://doi.org/10.1155/2020/8890450
Lower-Bound Solution of Foundation-Bearing Capacity under Circular Uniformly Distributed Load
The semispatial static allowable stress field was constructed by three-dimensional stress columns, and the analytical lower-bound solution of bearing capacity of Mohr–Coulomb foundation beneath circular uniformly distributed load was put forward for the first time. The influence of the amount of stress columns and soil shear strength parameters on the lower-bound solution of the foundation-bearing capacity is analyzed. The research showed that (1) when the number of stress columns is more than 9 (n > 4) and the friction angle is less than 30°, the relative error between the lower-bound solution of the bearing capacity and exact solution is less than 6.6%, (2) the suggested analytical solution is applicable for common clayey materials; thus, the application scope of the analytical lower-bound solution is further expanded, and (3) the influence of cohesion on bearing capacity is linear. However, with the increase of the friction angle, the bearing capacity increases faster and faster under all cohesion levels. The correctness and effectiveness of the proposed method were verified by literature comparison.
In the 1950s, Drucker and Prager combined the static field with the kinematical field and put forward the limit analysis method for geotechnical problems, which provided a new tool for solving the limit load of foundations. The limit analysis method includes upper-bound and lower-bound theorems. However, comparing with the kinematically allowable velocity field in the upper-bound method, it is more difficult to establish a static allowable stress field . Fortunately, finite-element technique can provide a powerful tool in solving some lower-bound problems. For instance, Li and Liu [2, 3] established the static allowable field and calculated the ultimate load by means of the finite-element method and linear programming. Singh and Basudhar  dealt with the prediction of lower-bound-bearing capacity of smooth (surface) and rough (embedded) strip footings based on finite elements and nonlinear programming. Huang  derived the lower-limit solution of the ultimate bearing capacity of a rigid strip footing based on finite element by using second-order cone programming. Furthermore, Sutcliffe, Yu, and Sloan  applied numerical limit analysis to evaluate the bearing capacity of strip footings for a jointed rock mass. Jiang and Sun  obtained the numerical solution of the ultimate bearing capacity of the composite pile foundation by using the lower-limit finite-element method. Zhang and Luan  deduced the lower-limit solution of the bearing capacity of a strip foundation under the composite loading mode. The solutions of the lower-bound-bearing capacity of reinforced soil-retaining wall, as well as shallow strip foundations in the vicinity of slopes, were successfully investigated [9, 10]. It could be concluded that research studies on the ultimate bearing capacity in lower-bound analysis concentrated on two-dimensional foundation, and three-dimensional cases were rarely reported [11–14]. However, the spatial effect of circular foundation cannot be ignored. Yang  obtained the analytical solution of the lower-limit-bearing capacity of circular foundation for Tresca materials, which is still one step distant from the ordinary foundation. By adopting Hoek–Brown failure criterion, Chakraborty and Kumar  evaluated the ultimate bearing capacity of a circular footing, placed over rock mass, by using the lower-bound theorem of limit analysis in conjunction with finite elements and nonlinear optimization.
In order to acquire the lower-bound solution of bearing capacity of Mohr–Coulomb foundation beneath circular footing, it is necessary to extend the two-dimensional stress plane to the three-dimensional stress space, which was proved to be feasible for pure cohesive soil . Furthermore, Chen  pointed out that since there was no new and physical complexity, it would be essentially possible to build a spatial stress field for soil foundation with φ > 0°, and the bearing capacity of three-dimensional footing on weightless foundation was sufficiently evaluated.
Based on the existing literature, stress columns were adopted to construct the static allowable stress field which conformed to the principle of the lower-bound theorem. Subsequently, the analytical lower-bound solution of the bearing capacity of Mohr–Coulomb foundation under the circular uniform load was presented. Finally, the results of examples were compared with acknowledged solutions and the reliability of the proposed analytical solution was verified.
2. Research Background and Lower-Bound Theorem
The bearing capacity of foundations is one of the three classical problems in soil mechanics. Figure 1(a) displays a typical load-deformation curve of loading tests on surficial footings. The curve contains four sections, which successively represent elastic, elastoplastic, plastic deformation, and strain hardening stage. In order to determine the ultimate bearing capacity (Pc) of the foundation as an elastoplastic material, two methods are usually adopted: one is to investigate the whole evolution process of the foundation, from the elastic deformation state to the plastic limit state, which is based on the theory of plastoelasticity mechanics (Figure 1(b)). Due to the complexity of practical engineering problems, it is usually only suitable for simple working conditions. The second is to ignore the process of elastic deformation, regard the material as a rigid-plastic body, and focus on the behavior of the structure in the plastic limit state, namely, the limit analysis method (Figure 1(c)). It can effectively simplify the analysis process and reflect the most essential content in the plastic deformation, so it is widely used in the engineering field.
If there exists a static allowable stress field throughout the whole object without any yield, which can balance with the load acting on the stress boundary, the object will never fail. When the deformation of the object reaches the limit state, the power of the real surface force in the given velocity field is always greater than (or equal to) the power of the corresponding surface force in the same velocity field of any other static allowable field. The lower-bound theorem indicates that (1) the ideal object can adjust itself to take on the potential external load and (2) among all the loads corresponding to the static allowable stress field, the ultimate load is the largest.
In general, the lower-bound theorem can be cooperated with the upper-bound theorem to define the actual interval of the ultimate load. If the upper limit is equivalent to the lower limit, then the ultimate load is a complete solution.
3. Spatial Stress Field and Lower-Bound Solution
According to the theoretical background of the lower-bound method [1, 15], following assumptions are adopted: (1) the foundation bears a circular uniformly distributed load; (2) the foundation is weightless; (3) the foundation is a rigid-plastic body, subjecting to the Mohr–Coulomb yield criterion.
3.1. A Single Stress Column
Similar to stress bands adopted in the two-dimensional lower-bound theorem , the concept of the three-dimensional stress column is firstly put forward and applied by Yang . It could be utilized to build the spatial static allowable stress field in lower-bound problems. The stress column is defined as an infinitely long “pile” suffering from uniaxial compression. For the simplest case, only one stress column beneath loading area is taken into consideration (Figure 2). Therefore, the lower-bound solution of bearing capacity of the foundation (σz) could be given byin which c and φ are the cohesion and the friction angle, respectively.
3.2. Five Stress Columns
Subsequently, a static allowable spatial stress field is shown in Figure 3, which is constructed by 5 stress columns. The stress column TO stands beneath the circular uniform load, and the horizontal projection of the cross section of the column totally coincides with the load distribution. The others, TA, TB, TC, and TD, are evenly distributed on a conical surface with an apex angle of 2α, subjecting to sin α = cos φ/2 [1, 15].
The axial compressive stresses of all columns are assumed to be P. In order to ensure that the stress state of everywhere in the foundation does not violate the material yield criterion after the superposition of stress fields, a horizontal stress column is added in the foundation with an axial compression of R = 2c cos φ/(1 − sin φ). The resultant stress field has omnidirectional symmetry in the horizontal direction. According to stress superposition principle, the resultant stresses along axes x, y, and z would be obtained by
Introducing the Mohr–Coulomb yield criterion,
Substituting (2), then
As a result, the lower-bound-bearing capacity could be stated as follows:
3.3. Nine Stress Columns
A more sophisticated static allowable stress field composed of 9 stress columns is shown in Figure 4. Similarly, the stress columns except for TO are evenly distributed on a conical surface with an apex angle of 2α. The axial compressive stresses of all columns are P, and a horizontal stress column with an axial compression of R is also supplemented.
The superposed stresses could be given by
According to the Mohr–Coulomb criterion,
Consequently, the lower-bound solution corresponding to 9 stress columns would be
3.4. Myriad Stress Columns
Above processes could be expanded to stress field constructed by 2n+1 (n is an even number) stress columns, including a vertical column and 2n equiangularly distributed columns on a conical surface with an apex angle of 2α. Also, a horizontal stress column with a compression of R is also supplemented. The superposed stresses could be calculated by the following equation:
Considering Mohr–Coulomb failure criterion,
So, we get
When the number of stress columns increases unlimitedly, as shown in Figure 5, the exact lower-bound solution of bearing capacity could be obtained by
As we know, the slip-line method and the limit equilibrium method can also be employed to solve geotechnical stability problems. The slip-line method tries to derive the plastic equilibrium differential equations and then obtains the corresponding solution by constructing the slip-line network. It is always complicated, and it relies on engineer’s experience and judgment to a great extent. On the other hand, a typical failure mode is specified in the limit equilibrium method, and the critical slip surface could be obtained by searching for the minimum safety factor. The deformation or velocity field is not considered in the limit equilibrium method, and the stress distribution beside the slip surface also remains lacking of enough validation. The most significant advantage of limit analysis is that no matter how complex the geometry and load of the structure are, a practical ultimate load can always be obtained conveniently. As presented above, the limit analysis method is competent to obtain the ultimate bearing capacity with such simplicity.
4. Verification and Evaluation
4.1. Literature Comparison
(1)By extending the 2D static allowable stress field to 3D environment, 5c was considered to be a lower-bound-bearing capacity of surficial rectangular or square footing with pure cohesive foundation (φ = 0°). By further improving the stress field, it would also be an eligible solution for all convex footings (such as circular, pentagon, or hexagon base) , not far from the solution of this research (4.8c).(2)Shield  investigated the bearing capacity of circular footing on soil layers with φ = 0°, as shown in Figure 6. B and H are the diameter of circular footing and the layer thickness below the base. As a result, the recommended upper-bound solution could be given by For semi-infinite homogeneous foundation space, B/H = 0; then, Consequently, the upper-bound-bearing capacity of the foundation with circular footing is slightly higher than the suggested lower-bound solution.(3)Yang  performed research on the analytical lower-bound solution of bearing capacity of circular footing upon Tresca foundation, and the comparative results with different stress column numbers are shown in Table 1. It shows no differences between these research studies; thus, the original solution is only the special case of this method at φ = 0°. However, Yang’s solution is not applicable in ordinary Mohr–Coulomb foundation (φ ≠ 0°), so the proposed solution in this research would have a larger scope of application.
4.2. The Effect of Stress Column Amount
From (4) and (5), the influence of the number of stress columns on lower-bound-bearing capacity could be easily evaluated. As shown in Figure 7, with the increase of stress column amount, the lower-bound solution gradually rises and finally converges to the exact solution determined by (12). For cases of φ = 0°, 10°, 20°, and 30°, the exact solution of σz/c is 4.800, 7.466, 12.209, and 21.192, respectively. Accordingly, Table 2 lists the numerical values and the relative errors of lower-bound solutions regarding variant amount of stress columns. It could be seen that for cases of n > 4 and φ < 30°, the relative errors of lower-bound solution of foundation-bearing capacity would be less than 6.6%. And the relative error will be reduced to less than 2% while n > 16 and φ < 30°. It means that a better lower-bound solution can be obtained by constructing a more sophisticated static allowable stress field. However, the improvement of lower-bound solution might be inefficient if the stress column amount is large enough.
4.3. The Effect of Shear Strength Parameters
According to (12), it is evident that the cohesion has a linear effect on lower-bound-bearing capacity. The effect of the friction angle is shown in Figure 8, in which the cohesion varies from 10 kPa to 40 kPa and the friction angle locates in the range of 0°–30°. The lower-bound-bearing capacity of the foundation with shear strength parameters of (c = 10 kPa, φ = 0°), (c = 40 kPa, φ = 0°), (c = 10 kPa, φ = 30°), and (c = 40 kPa, φ = 30°) would be 48 kPa, 192 kPa, 211.92 kPa, and 847.69 kPa, respectively. It is also indicated that for all cohesion levels, the growth rate of the bearing capacity gradually accelerates as the friction angle increases. The bearing capacity of φ = 30° is about 4.4 times that of φ = 0°. As Mohr–Coulomb materials, the bearing capacities of ordinary clayey foundations are much higher than those of Tresca material or pure cohesive soil. It can also partly explain why the base built in clay is more reliable than that in silt or sand.
Three-dimensional stress columns were adopted to build the half-spacial stress field for foundation beneath circular uniform load. The analytical lower-bound solution of bearing capacity of Mohr–Coulomb foundation, which greatly expanded the application scope of the original solution, was presented for the first time. Both effects of the stress column amount and shear strength parameters were evaluated, and the validity of the proposed solution was verified by the previous literature.
This study further enriches the relevant research results of the limit analysis method in the field of foundation-bearing capacity. Furthermore, the circular uniformly distributed load in this paper can be extended to any centrosymmetric distribution shape so long as the resultant force passes through the centroid vertically.
A part of the raw data cannot be shared at this time as the data also form part of an ongoing study.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
This work was funded by the Natural Science Foundation, China, grant number: 51408059.
- W. F. Chen, Limit Analysis and Soil Plasticity, Elsevier Science, Amsterdam, 1975.
- L. Li and B. C. Liu, “Lower bound limit analysis on bearing capacity of slope and its reliability,” Chinese Journal of Rock Mechanics and Engineering, vol. 20, no. 4, pp. 508–513, 2001.
- L. Li and B. C. Liu, “Ultimate bearing capacity on slope by lower bound analysis method,” China Railway Science, vol. 22, no. 1, pp. 79–83, 2001.
- D. N. Singh and P. K. Basudhar, “Optimal lower bound bearing capacity of strip footings,” Soils and Foundations, vol. 33, no. 4, pp. 18–25, 1993.
- Q. W. Huang, M. S. Huang, and G. H. Wang, “Calculation of bearing capacity of strip footings using lower bound limit method,” Chinese Journal of Geotechnical Engineering, vol. 29, no. 4, pp. 572–579, 2007.
- D. J. Sutcliffe, H. S. Yu, and S. W. Sloan, “Lower bound solutions for bearing capacity of jointed rock,” Computers and Geotechnics, vol. 31, no. 1, pp. 23–36, 2004.
- G. Sun, J. Sun, and X. D. Wang, “Finite element lower bound limit solutions for ultimate bearing capacity of subsoil under raft of composite piled foundations,” Chinese Journal of Geotechnical Engineering, vol. 33, no. 1, pp. 128–134, 2011.
- Q. Y. Zhang, M. T. Luan, and S. F. Zhao, “Lower bound limit analysis of bearing capacity of strip footing due to combined loading,” Journal of China University of Mining & Technology, vol. 37, no. 3, pp. 364–368, 2008.
- D. N. Singh and P. K. Basudhar, “Determination of the optimal lower-bound-bearing capacity of reinforced soil-retaining walls by using finite elements and non-linear programming,” Geotextiles and Geomembranes, vol. 12, no. 7, pp. 665–686, 1993.
- F. K. Ahmad, I. Ardavan, and J. C. Reza, “Seismic bearing capacity of shallow strip foundations in the vicinity of slopes using the lower bound finite element method,” Soils and Foundations, vol. 59, no. 6, pp. 1891–1905, 2019.
- J. X. Wang, Y. F. Wu, and Z. Li, “Lower bound theory of 3-D plastic limit analysis and its application,” Rock and Soil Mechanics, vol. 25, no. 10, pp. 1627–1631, 2004.
- R. T. Shield, “The plastic indentation of a layer by a flat punch,” Quarterly of Applied Mathematics, vol. 13, no. 1, pp. 27–46, 1955.
- R. T. Shield, “On Coulomb’s law of failure in soils,” Journal of the Mechanics & Physics of Solids, vol. 4, no. 1, pp. 10–16, 1955.
- H. J. Yang, Z. J. Shen, and J. H. Wang, “3D lower bound bearing capacity of smooth rectangular surface footings,” Mechanics Research Communications, vol. 30, no. 5, pp. 481–492, 2003.
- X. L. Yang and L. Li, “Discussion about lower-bound analytical solution of circular foundation bearing capacity using limit analysis method,” Journal of the China Railway Society, vol. 23, no. 6, pp. 77–80, 2001.
- M. Chakraborty and J. Kumar, “Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis,” Computers and Geotechnics, vol. 70, pp. 138–149, 2015.
Copyright © 2020 Fang Wei 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.