Zeitschriftendaten
Format
Zeitschrift
Erstveröffentlichung
09 Nov 2012
Erscheinungsweise
4 Hefte pro Jahr
Sprachen
Englisch
Open Access

Load capacity of the mixed bench and slab foundation. Numerical simulations and analytical calculation model

Akzeptiert: 03 Mar 2021
Zeitschriftendaten
Format
Zeitschrift
Erstveröffentlichung
09 Nov 2012
Erscheinungsweise
4 Hefte pro Jahr
Sprachen
Englisch
Introduction
Definition of the problem

Old, antique buildings are often founded on stone, masonry or concrete benches (Fig. 1). Modernization of the structure (to meet today's standards) often leads to situation where load capacity of the existing foundation benches is inadequate in comparison with loading. This is mainly due to load increase (when for example additional story is designed or way of the usage of the building is changed) or load capacity decreased (when basement is deepened). Sometimes these two situations occur simultaneously (Fig. 2).

In order to deal with such a problem, to increase load capacity of the existing foundation benches, several possible remedies are used (Fig. 3).

Dimensions of the existing foundations’ benches could be increased (both in horizontal and vertical direction). But such a technology requires digging out of the foundations, which is dangerous not only to the building structure itself, but also to the buildings in the vicinity. Limited available space (especially in urban environment) is also a serious problem in this technology. Another remedy is usage of the pressure injection or micropiles under the existing foundation benches. This method could be performed only by specialized geotechnical companies, which have a lot of experience in such a work. And still some danger for building structure itself and for buildings in the vicinity exists. Fourth possible method is to design additional reinforced concrete foundation slab between the existing foundation benches. Such a slab must be connected with existing benches, in order to form a united foundation. This is the easiest and cheapest method to increase the load capacity of the foundation. It does not require any additional excavations. But some design problems appear:

How to calculate the ultimate load of such a foundation when benches and slabs are on different levels, especially when soil conditions on the foundation level of the benches and slab are different?

How to calculate the internal forces (bending moments and shear forces) in the slab?

The second problem could be solved using the well-known slab on Winkler subsoil model. The first one is a main subject of this article.

Existing analytical calculations models for direct foundations’ ultimate load calculations

Problem of calculation of the shallow foundation ultimate load is very old. First works by Prandl (1920, 1921), later adopted by Reissner (1924), Keverling Buisman (1940) and Terzaghi (1943) lead to a simple analytical equation for the calculation of ultimate load. This approach, a bit modified by Meyerhof (1951, 1953, 1963) and Brinch Hansen (1970) is adopted by many codes, for example, by Eurocode-7 (EC-7). Outline of this approach (deterministic one) is presented below: $q=cNCbCsCiC+qNqbqsqiq+0.5γBNγbγsγiγ$ q = c{N_C}{b_C}{s_C}{i_C} + q{N_q}{b_q}{s_q}{i_q} + 0.5\gamma B{N_\gamma }{b_\gamma }{s_\gamma }{i_\gamma } $Nq=eπtgφtg2(45+φ/2)$ {N_q} = {e^{\pi tg\varphi }}t{g^2}\left( {45 + \varphi /2} \right) $Nc=(Nq−1)ctgφ$ {N_c} = \left( {{N_q} - 1} \right)ctg\varphi $Nγ=2(Nq−1)tgφ$ {N_\gamma } = 2\left( {{N_q} - 1} \right)tg\varphi where: q – ultimate load of the foundation [kPa], c – cohesion of the subsoil [kPa], Nq, Nc, Ng – bearing capacity coefficients [-], bc, bq, bγ – coefficients of load inclination, for vertical load =1 [-], sc, sq, sγ – coefficients of the foundation shape, for bench = 1 [-] $sq=1+B/Lsin(φ)$ {s_q} = 1 + B/L\sin \left( \varphi \right) $sγ=1−0.3B/L$ {s_\gamma } = 1 - 0.3B/L $sc=(sqNq−1)/(Nq−1)$ {s_c} = \left( {{s_q}{N_q} - 1} \right)/\left( {{N_q} - 1} \right) ic, iq, iγ – coefficients of the foundation inclination, for horizontal foundation =1 [-], q` – vertical stress on the foundation level outside of the foundation [kPa], ϕ – internal friction angle of the subsoil [∘], B – foundation width [m], L – foundation length [m], γ – soil volumetric weight [kN/m3].

It's worth noting that the equations for calculating the bearing capacity coefficients (as a functions of the internal friction angle) are still a subject of the investigations (Edgar Giovanny (2013), van Bars (2014)). Uncertainties of this model were discussed by Motra et al. (2016).

Wide discussion on the problem of partial material and load factors used in engineering calculations of direct foundations’ ultimate load is given by Bogusz and Godlewski (2019).

In the last few years, special attention was paid to the spatial variability of soil properties (mainly strength parameters—internal friction angle and cohesion), which leads to the usage of probabilistic approach in foundation safety investigation. In the work of Puła and Chwała (2015), soil strength parameters are considered as random fields, later averaged alongside the Prandtl's slip line. Random character of the slip line (as an effect of the random friction angle) is taken into account. Similar approach (but used to deal with seismic bearing capacity of strip footing) is presented by Johari et al. (2017). In the work of Johari et al. (2019), Young's modulus is treated as a random field in the analysis of differential settlement of strip footing. Stochastic finite element method (SFEM) is used. But linear elastic model is used to represent the soil behavior. This is far from the real soil behavior (nonlinearly-elastic in the small strain range and plastic close to the ultimate state).

But in the aforementioned methods for foundation's ultimate load estimation, a constant level of foundation (strictly speaking—constant vertical compressive stress on the foundation level outside of the foundation) is assumed, which is very often not true in the case of mixed bench and slab foundations. As far as the Author knows, there is no analytical solution for this problem in the literature.

Proposed calculation method

Numerical verification

In order to verify the above-proposed calculation method, numerical simulations were performed. FEM system ZSoil v18 was used. Elastic-plastic Mohr-Coulomb with cut-off condition (no tension) model was used to model the soil and elastic model for foundation. It is of course possible to use a more sophisticated model for soil (for example, Hardening Soil model). But it's not necessary when the main topic of interest is the ultimate load of the structure. Interface elements with no friction between foundation and soil were used to model the discontinued displacement field. For slab, beam elements were used in 2D plane strain analysis and shell elements in 3D analysis. Problem was run as a force-driven one. External load of the foundation was increased as long as a divergence of the iterative procedure occurs. Settlement – load curves were checked in order to judge if the divergence of the iterative process is caused by reaching of the ultimate load or only by some numerical instability (see Fig. 7).

First of all, the simple single bench ultimate load problem was considered, in order to compare the numerically obtained value of ultimate load with the results of Brinch Hansen (EC-7) approach. Foundation bench with 1 m width founded at 2 m below the terrain level on homogenous subsoil (medium capacity clay, see Table 1) was considered. 2D plane strain half-model was used to speed up the calculations (Fig. 5). Model consisted of 1957 Q4 2D-continuum elements, 2065 nodes and 19 interface elements for soil above the foundation level modelled as a continuum and 1271 Q4 2D-continuum elements, 1369 nodes and 4 interface elements for soil above the foundation level modelled as an equivalent load.

Parameters of soil used in the bearing capacity analysis.

Material E [MPa] γ [kN/m3] ϕ [°] c [kPa]
Medium capacity clay 15 20 10 10
Low capacity clay 5 18 5 5

Obtained value of ultimate load 280 kPa is significantly higher than the results of EC-7 calculations (187.5 kPa), difference is about 49%. Possible source of such a big difference is treating a soil above the foundation level as a source of the vertical compressive stress in the Brinch Hansen (EC-7) approach. So, the numerical simulation with soil above foundation level replaced by equivalent vertical load was performed. Obtained value of ultimate load 202 kPa is a bit higher than the results of EC-7 calculations (187.5 kPa), difference is about 8%. So, for simple problems, the numerical calculation gives values of an ultimate load quite close to those obtained from the Brinch Hansen (EC-7) approach, if soil above the foundation level is treated as a source of vertical stress only.

In the simulations presented in the later part of this paper, both approaches to the problem of soil above the foundation level (modelling as a Coulomb-Mohr continuum and as an equivalent load) are used, and the obtained results are compared.

Obtained failure mode is similar to those considered in the analytical solutions; triangle under the bench is visible in Fig. 6.

Then a more complicated case was analyzed. 2D plane strain half model of the benches and slab foundation (which corresponds to long foundation bench, B / L → 0) was built. 1 m width benches were founded on medium capacity clay 2 m below the terrain level, where 5.5 m width 30 cm thick slab—on low capacity clay 1 m below the terrain level (typical situation for antique buildings, see Table 1 and Fig. 8).

Model consisted of 2532 Q4 2D-continuum elements, 2694 nodes, 28 beam elements and 50 interface elements for soil above the foundation level modelled as a continuum and 1782 Q4 2D-continuum elements, 1915 nodes 28 beam elements and 35 interface elements for soil above the foundation level modelled as an equivalent load.

Obtained value of ultimate load 263.5 kPa (252 kPa with equivalent load approach) is a bit higher than the results of the proposed calculations method (221.5 kPa), difference is about 19% (14% with equivalent load approach). Obtained failure mode is illustrated in Fig. 9. Ultimate load of the slab itself (founded on a weak soil) calculated according to Brinch Hansen (EC-7) is 71.3 kPa; so such a conservative approach leads to great underestimation of the ultimate load, which is unacceptable from an economical point of view.

Also, 3D calculations of the ¼ of the real structure were performed (see numerical model shown in Fig. 10), in order to verify if the proposed calculations method works not only for long (B / L → 0) benches but also for situations where coefficients of the foundation shape have a significant influence on the obtained values of ultimate load.

Biggest of the used 3D model consisted of 27280 Q8 3D-continuum elements, 32296 nodes, 440 shell elements and 933 interface elements (soil above the foundation level modelled as a continuum, “long” foundation with B = 7.5 m and L = 28 m), where the smallest one consisted of 16040 Q8 3D-continuum elements, 19661 nodes, 104 shell elements and 693 interface elements (soil above the foundation level modelled as an equivalent load, “short” foundation with B = 7.5 m and L = 7.5 m).

Obtained result are illustrated in Fig. 11.

Obtained results show that treating soil above the foundation level as a Coulomb-Mohr continuum leads to obtain values of ultimate load about 17–19% higher than those obtained from the proposed approach. It's quite acceptable from an engineering and economical point of view. When using an equivalent load approach, differences drops down to 14% for “long” foundations and to about 5% for “short” (B/L = 1) ones. Obtained failure mode is similar to those observed in 2D simulations (Fig. 12).

Obtained load – settlement curves are almost identical for any points located on the bench. This is due to high stiffness of the bench. Maximal settlement difference is observed between the bench and the central point of the slab. The difference tends to vanish in the ultimate state (see Fig. 13). However, the values of the settlements presented here are only rough approximation of the reality. This is due to the simple constitutive Coulomb-Mohr model used and not taking into account the stiffness of the whole building. More sophisticated model, for example, Hardening Soil model, could be used to model the settlement process closer to the reality, but it will not change the obtained values of ultimate load.

Conclusions

Presented approach allows to adopt the Brinch-Hansen approach to the mixed bench and slab foundations. Calculated values of the ultimate load are a bit lower than those obtained from numerical simulations. Difference is between 17% and 19% when soil above the foundation level is treated (in numerical simulations) as a Coulomb-Mohr continuum and drops down to 5% for “short” (B/L = 1) foundations if equivalent load approach is used. Proposed approach works both for “short” (B/L = 1) and long (B/L = 0) foundations. So it can be used in engineering practice. Influence of the soil above foundation level on ultimate load of the structure is an additional margin of safety.

Parameters of soil used in the bearing capacity analysis.

Material E [MPa] γ [kN/m3] ϕ [°] c [kPa]
Medium capacity clay 15 20 10 10
Low capacity clay 5 18 5 5

Bogusz, W., Godlewski (2019), T. Philosophy of geotechnical design in civil engineering – possibilities and risks. Bulletin of the Polish Academy of Sciences Technical Sciences. 67(2), 289–306 BoguszW. GodlewskiT. 2019 Philosophy of geotechnical design in civil engineering – possibilities and risks Bulletin of the Polish Academy of Sciences Technical Sciences 67 2 289 306 Search in Google Scholar

Brinch Hansen, J. (1970) A revised extended formula for bearing capacity. Danish Geotechnical Institute Bulletin. 28, 5–11 Brinch HansenJ. 1970 A revised extended formula for bearing capacity Danish Geotechnical Institute Bulletin 28 5 11 Search in Google Scholar

Edgar Giovanny, D.S. (2013) Assessment of the range of variation of Nγ from 60 estimation methods for footings on sand. Canadian Geotechnical Journal. 50, 793–800 Edgar GiovannyD.S. 2013 Assessment of the range of variation of Nγ from 60 estimation methods for footings on sand Canadian Geotechnical Journal 50 793 800 Search in Google Scholar

Johari, A., Hosseini, S.M., Keshavarz A. (2017) Reliability analysis of seismic bearing capacity of strip footing by stochastic slip lines method. Computers and Geotechnics. 91, 203–217 JohariA. HosseiniS.M. KeshavarzA. 2017 Reliability analysis of seismic bearing capacity of strip footing by stochastic slip lines method Computers and Geotechnics 91 203 217 Search in Google Scholar

Johari, A., Sabzi A., Gholaminejad A. (2019) Reliability analysis of differential settlement of strip footings by stochastic response surface method (2019) Iran J Sci Technol Trans Civ Eng. 43, 37–48 JohariA. SabziA. GholaminejadA. 2019 Reliability analysis of differential settlement of strip footings by stochastic response surface method (2019) Iran J Sci Technol Trans Civ Eng 43 37 48 Search in Google Scholar

Keverling Buisman, A.S. (1940) Grondmechanica. Waltman, Delft, the Netherlands Keverling BuismanA.S. 1940 Grondmechanica Waltman Delft, the Netherlands Search in Google Scholar

Meyerhof, G.G. (1951) The ultimate bearing capacity of foundations. Géotechnique. 2, 301–332 MeyerhofG.G. 1951 The ultimate bearing capacity of foundations Géotechnique 2 301 332 Search in Google Scholar

Meyerhof, G.G. (1953) The bearing capacity of foundations under eccentric and inclined loads. In: Proceedings of the III International Conference on Soil Mechanics and Foundation Engineering, Zürich, Switzerland, 440–445 MeyerhofG.G. 1953 The bearing capacity of foundations under eccentric and inclined loads In: Proceedings of the III International Conference on Soil Mechanics and Foundation Engineering Zürich, Switzerland 440 445 Search in Google Scholar

Meyerhof, G.G. (1963) Some recent research on the bearing capacity of foundations. Canadian Geotechnical Journal. 1(1),16–26. MeyerhofG.G. 1963 Some recent research on the bearing capacity of foundations Canadian Geotechnical Journal 1 1 16 26 Search in Google Scholar

Motra, H.B., Stutz, H., Wuttke, F. (2016) Quality assessment of soil bearing capacity factor models of shallow foundations. Soils and Foundations. 56(2), 265–276 MotraH.B. StutzH. WuttkeF. 2016 Quality assessment of soil bearing capacity factor models of shallow foundations Soils and Foundations 56 2 265 276 Search in Google Scholar

Prandtl, L. (1920) Ueber die Haerte plastischer Koerper. Nachrichtender Gesellschaftder Wissenschaften. Berichte Mathem. -Phys. Kl, 74–85 PrandtlL. 1920 Ueber die Haerte plastischer Koerper. Nachrichtender Gesellschaftder Wissenschaften Berichte Mathem. -Phys Kl 74 85 Search in Google Scholar

Prandtl, L. (1921) Ueberdie Eindringfestigkeit (Haerte) plastischer Baustoffeund die Fes-tigkeit von Schneiden. Z. Angew. Math. Mech.1. Band1, 15–20 PrandtlL. 1921 Ueberdie Eindringfestigkeit (Haerte) plastischer Baustoffeund die Fes-tigkeit von Schneiden Z. Angew. Math. Mech.1. Band1 15 20 Search in Google Scholar

Puła, W., Chwała M. (2015) On spatial averaging along random slip lines in the reliability computations of shallow strip foundations. Computers and Geotechnics. 68, 128–136. PułaW. ChwałaM. 2015 On spatial averaging along random slip lines in the reliability computations of shallow strip foundations Computers and Geotechnics 68 128 136 Search in Google Scholar

Reissner, H. (1924) Zum Erddruckproblem. In: Biezeno, C.B., Burgers, J.M. (Eds.), Proceedings of the 1st International Congress for Applied Mechanics, Delft, The Netherlands: 295–311 ReissnerH. 1924 Zum Erddruckproblem In: BiezenoC.B. BurgersJ.M. (Eds.), Proceedings of the 1st International Congress for Applied Mechanics Delft, The Netherlands 295 311 Search in Google Scholar

Terzaghi, K. (1943) Theoretical Soil Mechanics. J. Wiley, New York. TerzaghiK. 1943 Theoretical Soil Mechanics J. Wiley New York Search in Google Scholar

van Baars, S. (2014) The inclination and shape factors for the bearing capacity of footings. Soils and Foundations. 54(5), 985–992 van BaarsS. 2014 The inclination and shape factors for the bearing capacity of footings Soils and Foundations 54 5 985 992 Search in Google Scholar

• Usefulness of the CPTU method in evaluating shear modulus G0 changes in the subsoil

Recommended articles from Trend MD