Next Article in Journal
Acyclic Nucleic Acids with Phosphodiester Linkages—Synthesis, Properties and Potential Applications
Next Article in Special Issue
Evaluating the Design Criteria for Light Embankment Piling: Timber Piles in Road and Railway Foundations
Previous Article in Journal
Clinical Interpretation of Working Volume and Weight Support in Upper Limb Robotic Neurorehabilitation after Stroke
Previous Article in Special Issue
Study on the Earth Pressure during Sinking Stage of Super Large Caisson Foundation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Settlement and Bearing Capacity of Rectangular Footing in Reliance on the Pre-Overburden Pressure of Soil Foundation

by
Zaven G. Ter-Martirosyan
,
Armen Z. Ter-Martirosyan
and
Huu H. Dam
*
Department of Soil Mechanics and Geotechnical Engineering, National Research Moscow State Civil Engineering University, 26, Yaroslavskoye Shosse, 129337 Moscow, Russia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(24), 12124; https://doi.org/10.3390/app112412124
Submission received: 16 November 2021 / Revised: 6 December 2021 / Accepted: 8 December 2021 / Published: 20 December 2021
(This article belongs to the Special Issue Advances in Geotechnical Engineering Ⅱ)

Abstract

:
This article presents a solution for the quantitative evaluation of the stress–strain state (SSS) and the bearing capacity of rectangular foundations, factoring in the unit weight of the soil mass and different values of pre-overburden pressure (POP). In order to assess the SSS of the soil subgrade below a rigid rectangular footing under a uniformly distributed load, the authors applied the Boussinesq basic solution for an elastic half-space subjected to a vertical point load on its surface. As a result, the formulas for vertical stress, mean stress, shear strain, and volumetric strain for any point in Cartesian coordinates (x, y, z) and foundation settlement were determined. Additionally, the application of Hencky’s system of physical equations, with non-linear dependencies between mean stress and volumetric strain as well as deviator stress and shear strain, along with the experimental curves, depicts the relationships between bulk modulus and volume stress, and shear modulus and shear stress. The authors point out the non-linear behavior of the subgrade soil and propose a method for estimating the bearing capacity of a rigid rectangular foundation.

1. Introduction

The solutions for the problem of settlement of rigid rectangular footings have been proposed by a variety of authors [1,2,3,4,5,6] with various approaches. However, the limitation of which is a lack of explanation for the failure mechanism of soil medium underneath a rigid rectangular footing. By considering the vertical displacement as a result of components due to shear strain and volumetric strain of a loaded soil element, the authors denote not only the non-linear dependence between applied load and strain components but also the approach to the bearing capacity of rigid rectangular footings. To calculate shear strain and volumetric strain, the elastic–perfectly plastic model of Timoshenko [7] and the non-linear model proposed by Grigoryan are applied to indicate the non-linear behavior of the soil subgrade below a rectangular footing subjected to a uniformly distributed load, thereby achieving a dependency of deformation on load with a double curvature. So far, the non-linear analysis of soil subgrade under a rectangular foundation has solely been conducted using the finite element method. Therefore, it is essential to consider the non-linear soil deformation properties in analytical solutions based on the results from compression and triaxial tests. As a result, the varying values of shear modulus and bulk modulus depending on the deviator stress are elucidated clearly at any point in Cartesian coordinates. In addition, linear behavior can also be observed in the non-linear solution, when considering the shear modulus and bulk modulus as constant values. Consequently, the non-linear analysis turns out to be a classic linear analysis, based on Hooke’s law (G = const and K = const).
In this article, for the purpose of estimating the settlement and bearing capacity, we applied Hencky’s system of equations, which enables engineers to consider both linear and non-linear behaviors of soil media under a rectangular footing as well as to thoroughly separate the total vertical strain ε z into shear strain ε γ and volumetric strain ε v as follows:
ε z = ε γ + ε v = σ z σ m G ( σ m , τ i / τ i ) + σ m K ( σ m )
in which G ( σ m , τ i / τ i ) and K ( σ m ) are, respectively, shear modulus and bulk modulus, the values of which are relevant to the mean stress σ m as well as the ratio of mobilized shear stress τ i to maximum shear stress τ i , which can be found using the Mohr–Coulomb criterion at the considered points.
In a particular case, when G = const and K = const, Equation (1) turns into the Hooke’s law equation [8] that is presented as:
ε z = σ z σ m 2 G + σ m K
G ( σ m , τ i / τ i ) and K ( σ m ) can be retrieved from standard triaxial tests, the general view of which is observed in Figure 1 below.
A significant advantage of using Hencky’s system of equations is that it allows for the determination of the typical characteristics of total strain, volume, and shear strain of loaded soil elements. Figure 1 accurately depicts real soil behavior, with a decaying dependency between volumetric strain and mean stress, and a persistent dependency between shear strain and deviator stress. The final state following the latter dependency leads soil elements to failure. In other words, the progression of shear strain due to the increment of deviator stress induces an unlimited increase in the vertical deformation of soil samples in triaxial tests and the appearance of failure surfaces in the soil subgrade under rectangular and strip footings, while the volumetric strain is solely increased to stable values. This means that the volume or the porosity of soil elements will reach constant values when the critical state is achieved.

2. Methods and Analysis

2.1. Theoretical Basics for Predicting the Settlement of a Rectangular Foundation in a Linear Formulation

As a computational geo-mechanical model for describing the stress state, a linearly deformable half-space under the influence of a constantly concentrated load q is considered, calculation schema of which is depicted in Figure 2. According to Boussinesq [9], the equations for the Cartesian normal stress and the mean stress are presented as follows:
σ p . z ( x , y , z ) = 3 q 2 π z 3 R 5
σ p . m ( x , y , z ) = q 3 π ( 1 + ν ) z ( x 2 + y 2 + z 2 ) 3
σ p . x = q 2 [ 3 x 2 z R 5 ( 1 2 ν ) ( x 2 y 2 ( x 2 + y 2 ) R ( z + R ) + y 2 z ( x 2 + y 2 ) R 3 ) ]
σ p . y = q 2 [ 3 y 2 z R 5 ( 1 2 ν ) ( y 2 x 2 ( x 2 + y 2 ) R ( z + R ) + x 2 z ( x 2 + y 2 ) R 3 ) ]
in which R = ( x 2 + y 2 + z 2 ) .
Based on Figure 3, in order to determine the distribution of all the stresses above at a depth beneath a rectangular footing, integrations Figure of expressions (3), (4), (5), and (6) with respect to η (from −a to a) and ξ (from −b to b) are carried out. However, in order to optimize the analytical calculations, the authors propose solely integrating Equations (3) and (4). The results are denoted as the following:
σ z ( x , y , z ) = 3 q z 3 2 π a a b b d η d ξ ( ( x η ) 2 + ( y ξ ) 2 + z 2 ) 5 = = q 2 π { arctg [ ( x a ) ( y b ) z ( x a ) 2 + ( y b ) 2 + z 2 ] + arctg [ ( x + a ) ( y + b ) z ( x + a ) 2 + ( y + b ) 2 + z 2 ] arctg [ ( x a ) ( y + b ) z ( x a ) 2 + ( y + b ) 2 + z 2 ] arctg [ ( x + a ) ( y b ) z ( x + a ) 2 + ( y b ) 2 + z 2 ] } + q 2 π { z ( x + a ) ( y + b ) [ ( x + a ) 2 + ( y + b ) 2 + 2 z 2 ] [ ( x + a ) 2 + z 2 ] [ ( y + b ) 2 + z 2 ] ( x + a ) 2 + ( y + b ) 2 + z 2 z ( x + a ) ( y b ) [ ( x + a ) 2 + ( y b ) 2 + 2 z 2 ] [ ( x + a ) 2 + z 2 ] [ ( y b ) 2 + z 2 ] ( x + a ) 2 + ( y b ) 2 + z 2 + + z ( x a ) ( y b ) [ ( x a ) 2 + ( y b ) 2 + 2 z 2 ] [ ( x a ) 2 + z 2 ] [ ( y b ) 2 + z 2 ] ( x a ) 2 + ( y b ) 2 + z 2 z ( x a ) ( y + b ) [ ( x a ) 2 + ( y + b ) 2 + 2 z 2 ] [ ( x a ) 2 + z 2 ] [ ( y + b ) 2 + z 2 ] ( x a ) 2 + ( y + b ) 2 + z 2 }
σ m ( x , y , z ) = q z 3 π ( 1 + ν ) a a b b d η d ξ ( ( x η ) 2 + ( y ξ ) 2 + z 2 ) 3 = q ( 1 + ν ) 3 π { arctg [ ( x a ) ( y b ) z ( x a ) 2 + ( y b ) 2 + z 2 ] + arctg [ ( x + a ) ( y + b ) z ( x + a ) 2 + ( y + b ) 2 + z 2 ] arctg [ ( x a ) ( y + b ) z ( x a ) 2 + ( y + b ) 2 + z 2 ] arctg [ ( x + a ) ( y b ) z ( x + a ) 2 + ( y b ) 2 + z 2 ] }
While considering the obtained expressions (7) and (8) on the vertical axis z (x = 0; y = 0), their brief forms can be expressed as follows:
σ z = 2 q π [ arctg a b z a 2 + b 2 + z 2 + a b z ( a 2 + b 2 + 2 z 2 ) ( a 2 + z 2 ) ( b 2 + z 2 ) a 2 + b 2 + z 2 ]
σ m = 4 q ( 1 + ν ) 3 π [ arctg ( a b z a 2 + b 2 + z 2 ) ]
Under the generalized Hooke’s Law, the total strain ε z of a soil element can be considered through two components, volumetric strain ε v and shear strain ε γ , which are dependent on vertical stress, mean stress, and the stiffness of the soil.
ε z = ε γ + ε v = σ z σ m 2 G + σ m K
The obtained expressions above reveal that, for the linear elastic soil medium, the shear strain and volumetric strain of soil elements at any point in Cartesian coordinates (x, y, z) can be determined using expressions of σ z and σ m along with the stiffness parameters. In the particular case of a rigid rectangular footing, in order to evaluate its settlement, the expression (11) requires an integration with respect to z (from 0 to H, where H is the depth of the compressed thickness of the soil massif under the footing). Therefore, the total settlement is presented as follows:
S = S γ + S v = 0 H ε γ d z + 0 H ε v d z = 0 H σ z σ m 2 G d z + 0 H σ m K d z
in which S γ and S v are the values of settlement due to shear strain and volumetric strain, respectively.
For the purpose of convenience, the authors propose an analytical method of total settlement assessment using expressions of σ z and σ m on the axis that passes through the center (0, 0, 0) of a rigid rectangular footing. The expressions of S γ and S v are illustrated as:
S γ = 0 H σ z σ m 2 G d z = 0 H [ q ( 1 2 ν ) 3 π G arctg a b z a 2 + b 2 + z 2 + q π G a b z ( a 2 + b 2 + 2 z 2 ) ( a 2 + z 2 ) ( b 2 + z 2 ) a 2 + b 2 + z 2 ] d z
S v = 0 H σ m K d z = 4 q ( 1 + ν ) 3 π K 0 H arctg ( a b z a 2 + b 2 + z 2 ) d z
The analyzed expression of settlement S γ due to shear strain consists of three components A B C as follows:
S γ = A + B + C
in which
A = q ( 1 2 ν ) 3 π G h   arctg a b h a 2 + b 2 + h 2
B = q ( 1 2 ν ) 3 π G { H 1 ln | 8 H 1 3 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 + h 2 8 H 1 3 + 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 + h 2 | + H 2 ln | 8 H 2 3 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 + h 2 8 H 2 3 + 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 + h 2 | H 1 ln | 8 H 1 3 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 8 H 1 3 + 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 | H 2 ln | 8 H 2 3 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 8 H 2 3 + 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 | }
C = q 2 π G [ a ln | a 2 + b 2 + h 2 b a 2 + b 2 + h 2 + b | a ln | a 2 + b 2 b a 2 + b 2 + b | + + b ln | a 2 + b 2 + h 2 a a 2 + b 2 + h 2 + a | | b ln | a 2 + b 2 a a 2 + b 2 + a | | ]
H 1 = ( a 2 + b 2 ) ( a 2 + b 2 ) 2 4 a 2 b 2 8 ;   H 2 = ( a 2 + b 2 ) + ( a 2 + b 2 ) 2 4 a 2 b 2 8
The settlement S v due to the volumetric strain caused by mean stress can be procured by analyzing expression (12) with the form of:
S v = 4 q ( 1 + ν ) 3 π K h arctg a b h a 2 + b 2 + h 2 + + 4 q ( 1 + ν ) 3 π K H 1 ln | 8 H 1 3 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 + h 2 8 H 1 3 + 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 + h 2 | 4 q ( 1 + ν ) 3 π K H 1 ln | 8 H 1 3 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 8 H 1 3 + 2 ( a 2 + b 2 ) H 1 + a b a 2 + b 2 | + + 4 q ( 1 + ν ) 3 π K H 2 ln | 8 H 2 3 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 + h 2 8 H 2 3 + 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 + h 2 | 4 q ( 1 + ν ) 3 π K H 2 ln | 8 H 2 3 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 8 H 2 3 + 2 ( a 2 + b 2 ) H 2 + a b a 2 + b 2 |
Considering the example with the following inputs: q1 = 1000 kPa; a1 = 0.5 m; b1 = 1.5 m; ν 1 = 0.3; H1 = 20 m; G1 = 10,000 kPa; K1 = 65,000 kPa; q2 = 1500 kPa; a2 = 1.5 m; b2 = 3 m; ν 2 = 0.3; H2 = 25 m; G2 = 15,000 kPa; K2 = 97,500 kPa; q3 = 1500 kPa; a3 = 2 m; b3 = 4 m; ν 3 = 0.3; H3 = 30 m; G3 = 20,000 kPa; and K3 = 130,000 kPa. By using math software for engineering calculations (Mathcad), we obtain: S γ 1 = 0.06105 m; S v 1 = 0.01436 m; S γ 2 = 0.1180 m; S v 2 = 0.02911 m; S γ 3 = 0.1556 m; and S v 3 = 0.03859 m. It is substantially important to note that the values of S γ 1 , S γ 2 , and S γ 3 are 4.25, 4.06, and 4.032 times higher than the values of S v 1 , S γ 2 , and S v 3 ( S γ 1 = 4.25 S v 1 ; S γ 2 = 4.06 S v 2 ; S γ 3 = 4.032 S v 3 ), respectively. Thus, it can be concluded that in linear elastic half-space, the settlement of a loaded soil subgrade mainly occurs due to the shear strain, which is over 80% of the total value, while nearly 20% of the total value comes from the volume change. The distributions of shear strain ε γ ( z ) and volumetric strain ε v ( z ) , calculated using expression (11), with a depth corresponding to the first group of parameters (q1, a1, b1, ν 1 , H1, G1, K1) are illustrated in Figure 4 and Figure 5.
As presented in expression (11), the shear ε γ ( z ) and volumetric strains ε v ( z ) are calculated through the deviator stress σ z σ m and volumetric stress σ m , which are derived from expressions (7) and (8). As a result, the variation of ε γ ( z ) , and ε v ( z ) with depth is accompanied by a change of σ z σ m and σ m , respectively. Therefore, the decrease in values of σ z σ m and σ m as the depth increases leads to the dip in ε γ ( z ) and ε v ( z ) .

2.2. Hencky’s System of Physical Equations

Apart from the solution for the linear relationships above, the usage of Hencky’s system of physical equations [10] allows modeling of non-linear relationships between stress and strain (see also [11]), which are mathematically described as follows:
ε x = χ ( σ x σ m ) + χ σ m ; ε γ . xy = 2 χ . τ xy ε y = χ ( σ y σ m ) + χ σ m ; ε γ . yz = 2 χ . τ yz ε z = χ ( σ z σ m ) + χ σ m ; ε γ . zy = 2 χ . τ zy
in which
χ = γ i 2 τ i = f ( τ i , σ m , μ σ ) 2 τ i   and   χ = ε m σ m = f ( τ i , σ m , μ σ ) σ m
τ i is the value of current shear stress of a loaded soil element and μ σ is the Lode–Nadai coefficient, which characterizes the forms of the stress state.
It is noted that, when χ and χ take the forms of 1 2 G and 1 K , respectively, Hencky’s system of physical equations turns into the Hooke’s law formulas for soil elements, in which G = E 2 ( 1 + ν ) and K = E ( 1 2 ν ) are applied as the primary stiffness parameters.
In order to model the non-linear relationship between the volumetric strain and mean stress of loaded soil elements, we used an exponential function dependency proposed by Grigoryan [12]. Its form is as follows:
ε m ( σ m ) = ε ( 1 e α σ m )
where ε is the limiting value of volumetric deformation, ε m ( σ m ) is the average volumetric deformation with varying values depending on mean stress σ m , and α is a non-linearity parameter, which can be determined experimentally.
By considering the elementary definition of volume modulus observed by the ratio of volumetric strain to volume stress, the varying values of secant bulk modulus K, which depend on the stress state of the loaded soil elements, can be represented by the following form:
ε m σ m = 1 K = ε ( 1 e α σ m ) σ m K = σ m ε m
Being close to infinite, ε m will eventually reach ε when α = 0 and ε = ε m , and expression (24) will have the form of Hooke’s law for linear elastic materials.
To express the elastic–plastic behaviors of soil elements subjected to shear stress, the formula proposed by Tymoshenko [13,14] can be applied.
ε γ . i = τ i G c e τ i τ i τ i
where τ i is the maximum value of shear stress, which is calculated by the following equation:
τ i = ( σ m + σ g ) t g ϕ i + c i
in which σ m is the mean stress of loaded soil elements; ϕ i and c i are the angle of internal friction and the cohesive strength, respectively, determined by the Mohr diagram and failure envelopes; and σ g is the vertical earth pressure, which in the condition of normally consolidated soil is represented by the form:
σ g = γ h
where γ is the unit weight of soil and h is the depth at which τ i is estimated.
In the case of over-consolidated soil with the maximum effective vertical overburden stress that the soil element has sustained in the past, the formula of σ g takes the following form:
σ g = γ h + σ P
where σ P is pre-overburden pressure, which can be determined through pre-consolidation pressure using Casagrande’s graphical method [15] based on the compression test results. The instructions for the method can be observed in the Russian Federation’s national standard on soils, under laboratory method for determining the over-consolidation characteristics. Equation (25) implies that:
G = G e ( 1 τ i τ i )
in which G e is the value of elastic shear modulus of the soil.
It follows from expression (29) that for τ i τ i , we procure G 0 and ε γ . i . Figure 6a–c demonstrate respectively the colour shadings of mean stress σ m , vertical stress σ z , and deviator stress σ z σ m that are derived from Equations (7) and (8) in the XOZ-plane.

3. Results and Discussion

3.1. Calculation of Total Strain in Soil Layers Taking into Account the Non-Linear Deformability and Their Various Mechanical Parameters

Throughout this article, it has been emphasized that when estimating the settlement of a rectangular footing, Hencky’s system of physical equations prominently manifests its special feature of analyzing the total deformations as shear and volumetric strains, taking into consideration the varying values of soil stiffness parameters, which depends on the stress-state of the soil elements. The equation for total strain is briefly displayed as follows:
ε z = ε γ + ε v = σ z σ m G ( σ m , τ i / τ i ) + σ m K ( σ m )
Putting Equations (24) and (29) into the procured expression (30) above, we get:
ε z = σ z σ m 2 G e ( 1 τ i τ i ) + ε ( 1 e α σ m )
in which the definitions of G e , τ i , ε , σ z , and σ m are detailed above.
According to [16], Mohr’s circle reveals that the maximum shear stress in soil elements takes a value that is equal to half of the difference between the maximum principal stress σ 1 and σ 3 . On the axis that passes through the center (0, 0, 0) of a rigid rectangular footing, the vertical stress σ z and the two components σ x and σ y act as the principal stress σ 1 , σ 2 , and σ 3 , respectively. Therefore, the value of shear stress that occurs in a cubic soil element τ i can be accepted as:
τ i = σ 1 σ 3 2 = σ z σ x 2
where σ x is calculated by integrating σ p . x , which corresponds to the case of point load with respect to η (from −a to a) and to ξ (from −b to b).
σ x = a a b b q 2 [ 3 ( x η ) 2 z R 1 5 ( 1 2 ν ) ( ( x η ) 2 ( y ξ ) 2 ( ( x η ) 2 + ( y ξ ) 2 ) R 1 ( z + R 1 ) + ( y ξ ) 2 z ( ( x η ) 2 + ( y ξ ) 2 ) R 1 3 ) ] d η d ξ
R 1 = ( x η ) 2 + ( y ξ ) 2 + z 2
With a rigid square footing, thanks to its symmetrical shape, the stress–strain state (SSS) of the soil subgrade under the footing can be considered as the SSS of a soil sample in a symmetric triaxial test, in which σ 3 and σ 2 are the same. Therefore, the equation τ i can be written in the following form (see also [17]):
τ i = σ z σ x 2 = σ z 3 σ m σ z 2 2 = 3 2 ( σ z σ m )
where σ z and σ m are calculated by Equations (9) and (10).
From analysis of Equation (31), it can be concluded that the component deformations grow according to the loading increment. The curve of the shear strain ε γ takes a parabolic form and has a tendency of infinity, while the progress of the volumetric strain ε v tends to be dampened. As a sum of the two components ε γ and ε v , the total strain of a loaded soil element possesses a double curvature and turns out to be infinite at τ i τ i . Dependency curves between shear strain, volumetric strain, total strain, and loading values with various values of stiffness and strength parameters can be observed using Equation (30). Figure 7 and Figure 8 indicate the PTC Mathcad results for Equation (31).

3.2. Calculation of Total Settlement of Soil Layers Taking into Account the Non-Linear Deformability and Their Various Mechanical Parameters

The dependency curve of total settlement S(q) on surface load q is derived by integrating Equation (31) with respect to z (from 0 to H, where H is the depth of the compressed thickness of the soil massif under the footing). This is represented as follows:
S ( q ) = S γ ( q ) + S v ( q ) = 0 H σ z ( q ) σ m ( q ) 2 G e [ 1 1 2 σ z ( q ) σ x ( q ) ( σ m ( q ) + γ h + σ P ) t g ϕ + c ] d z + 0 H ε ( 1 e α σ m ( q ) ) d z
in which S γ ( q ) and S v ( q ) are the shear and volume components, respectively, that are determined as follows:
S γ ( q ) = 0 H σ z ( q ) σ m ( q ) 2 G e [ 1 1 2 σ z ( q ) σ x ( q ) ( σ m ( q ) + γ h + σ P )     t g ϕ + c ] d z
S v ( q ) = 0 H ε ( 1 e α σ m ( q ) ) d z
Along with the increment in values of surface load q, the rise in settlement S(q) will arrive at an infinite value based on the strength parameters of the soil. The value of surface load q corresponding to the infinity state can be observed as the bearing capacity of the rigid rectangular footing. Figure 9 shows dependency plots for S γ ( q ) 1 , S v ( q ) 1 , and S ( q ) 3 , determined by Equation (35), while Figure 10 shows dependency plots for S ( q ) at different values of pre-overburden pressure σ P , and Figure 10 demonstrates dependency plots for S ( q ) at various values of foundation depth z.

4. Conclusions

  • The settlement and bearing capacity of rectangular footings are regarded as the main calculation parameters when designing foundations used in construction.
  • These calculations show that geo-mechanical models for foundation, including geometric parameters of the footing (length, width, depth, initial and boundary conditions) and computational soil models, such as linear and non-linear models with different systems of physical equations (Hooke and Hencky), have a significant influence on the evaluation of settlement and bearing capacity.
  • The computational soil model used in this article, with the possibility of lateral expansion ( ε x 0 , and ε y 0 ), along with the elastic–plastic model for shear strain–stress dependency and the non-linear model for volumetric strain–stress dependency as parts of Hencky’s system of physical equations, allows us to present the total deformation as the sum of the volumetric and shear components ( ε z = ε γ + ε v ). Only in this case do the strain–stress curve and the settlement–load curve both have decaying and persistent (double curvature) characteristics.
  • The computational soil model proposed here indicates that, thanks to the choice of different combinations of stiffness (Ge, ν , ε , α ) and strength parameters ( ϕ and c), the curve of dependency between settlement and load can be observed as decaying or persistent with various bending degrees of the curvature.
  • This article presents a substantial dependency of total deformation on stress as well as total settlement on load, depending on the various values of pre-overburden pressure σ P . The growth of pre-overburden pressure σ P leads to an increase in the bearing capacity of the foundation soil, and a decrease of total settlement compared to the case with lower values of σ P .

Author Contributions

Conceptualization, methodology, investigation, software, Z.G.T.-M. and A.Z.T.-M.; formal analysis, H.H.D. and A.Z.T.-M.; writing—original draft preparation, Z.G.T.-M. and A.Z.T.-M.; writing—review and editing, all authors; visualization, H.H.D.; supervision, Z.G.T.-M. and A.Z.T.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Ministry of Science and Higher Education (grant # 075-15-2021-686). All tests were carried out using research equipment of The Head Regional Shared Research Facilities of the Moscow State University of Civil Engineering.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are included within the article. The original details of the data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dempsey, J.P.; Li, H. A rigid rectangular footing on an elastic layer. Geotechnique 1989, 39, 147–152. [Google Scholar] [CrossRef]
  2. Butterfield, R.; Banerjee, P.K. A rigid disc embedded in an elastic half space. Geotech 1971, 2, 35–52. [Google Scholar]
  3. Mullan, S.J.; Sinclair, G.B.; Brothers, P. Stresses for an elastic half-space uniformly indented by a rigid rectangular footing. Int. J. Numer. Anal. Methods Geomech. 1980, 4, 277–284. [Google Scholar] [CrossRef]
  4. Nwoji, C.; Onah, H.N.; Mama, B.; Ike, C.C. Solution of Elastic Half Space Problem using Boussinesq Displacement Potential Functions. Asian J. Appl. Sci. 2017, 5. [Google Scholar]
  5. Pantelidis, L. Elastic Settlement Analysis for Various Footing Cases Based on Strain Influence Areas. Geotech. Geol. Eng. 2020, 38, 4201–4225. [Google Scholar]
  6. Pantelidis, L.; Gravanis, E. Elastic Settlement Analysis of Rigid Rectangular Footings on Sands and Clays. Geosciences 2020, 13, 491. [Google Scholar] [CrossRef]
  7. Timoshenko, S.P. Theory of Elasticity (Teoriya Elastichnosti); Nauka: Moscow, Russia, 1975. [Google Scholar]
  8. Whitlow, R. Basic Soil Mechanics; Pearson Prentice Hall: New York, NY, USA, 2001. [Google Scholar]
  9. Boussinesq, J. Application des Potentiels à L’étude de L’équilibre et du Mouvement des Solides Élastiques: Principalement au Calcul des Déformations et des Pressions que Produisent, dans ces Solides, des Efforts Quelconques Exercés sur une Petite Partie de Leur Interieur; Gauthier-Villars: Paris, France, 1885. [Google Scholar]
  10. Hencky, H. Zur Theorie plastischer Deformationen und der hierdurch im Material hervorgerufenen Nachspannungen. Zammzeitschrift Fur Angew. Math. Und Mech. 1924, 4, 323–334. [Google Scholar] [CrossRef]
  11. Ter-Martirosyan, Z.G. Predicting the settlement of foundation based taking into account induced anisotropy. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1015, 012048. [Google Scholar] [CrossRef]
  12. Florin, V.A. Fundamentals of Soil Mechanics, Volume I; Gosstroyizdat: Moscow, Russia, 1959. (In Russian) [Google Scholar]
  13. Vyalov, S.S. Rheological Basis of Soil Mechanics (Reologicheskie Osnovy Mekhaniki Gruntov); Vyshaya shkola: Moscow, Russia, 1978. [Google Scholar]
  14. Trofimov, V.T. Soil Science (Gruntovedenie), 6th ed.; Nauka: Moscow, Russia, 2005. [Google Scholar]
  15. Casagrande, A. The determination of pre-consolidation load and itʼs practical significance. In Proceeding of the 1st International Conference on Soil Mechanics, Cambridge, MA, USA, 22–26 June 1936. [Google Scholar]
  16. Labuz, J.F.; Zang, A. Mohr–Coulomb Failure Criterion. Rock Mech. Rock Eng. 2012, 45, 975–979. [Google Scholar] [CrossRef] [Green Version]
  17. Ter-Martirosyan, Z.G. Soil Mechanics (Mekhanika Gruntov); ASV Publisher: Moscow, Russia, 2009. [Google Scholar]
Figure 1. Schematic representation of the results of standard triaxial soil tests under a kinematic loading mode or along the path of crush.
Figure 1. Schematic representation of the results of standard triaxial soil tests under a kinematic loading mode or along the path of crush.
Applsci 11 12124 g001
Figure 2. Stress in an elastic half-space caused by a vertical point load on the surface.
Figure 2. Stress in an elastic half-space caused by a vertical point load on the surface.
Applsci 11 12124 g002
Figure 3. Stress in an elastic half-space caused by a uniformly distributed load on the surface.
Figure 3. Stress in an elastic half-space caused by a uniformly distributed load on the surface.
Applsci 11 12124 g003
Figure 4. Distribution of ε γ ( z ) with depth z.
Figure 4. Distribution of ε γ ( z ) with depth z.
Applsci 11 12124 g004
Figure 5. Distribution of ε v ( z ) with depth z.
Figure 5. Distribution of ε v ( z ) with depth z.
Applsci 11 12124 g005
Figure 6. Color shadings of σ m (a), σ z (b), and σ z σ m (c) (at a = 0.75 m; b = 1.5 m; q = 1000 kPa) and the calculation scheme for calculating shear and volumetric strain based on Hencky’s system of physical equations (d).
Figure 6. Color shadings of σ m (a), σ z (b), and σ z σ m (c) (at a = 0.75 m; b = 1.5 m; q = 1000 kPa) and the calculation scheme for calculating shear and volumetric strain based on Hencky’s system of physical equations (d).
Applsci 11 12124 g006aApplsci 11 12124 g006b
Figure 7. The dependency curve between ε γ 1 , ε v 2 , ε z 3 , calculated using Equation (31) and load q of a soil element at a depth of 5 m.
Figure 7. The dependency curve between ε γ 1 , ε v 2 , ε z 3 , calculated using Equation (31) and load q of a soil element at a depth of 5 m.
Applsci 11 12124 g007
Figure 8. The dependency curve between ε z 3 , calculated using Equation (31) and load q of a soil element at a depth of 5 m with various values of stiffness ( G e , ν , ε , α ) and strength ( ϕ , and c) parameters.
Figure 8. The dependency curve between ε z 3 , calculated using Equation (31) and load q of a soil element at a depth of 5 m with various values of stiffness ( G e , ν , ε , α ) and strength ( ϕ , and c) parameters.
Applsci 11 12124 g008
Figure 9. The dependency curve between S γ ( q ) 1 , S v ( q ) 1 , and S ( q ) 3 , calculated using Equation (35), and load q.
Figure 9. The dependency curve between S γ ( q ) 1 , S v ( q ) 1 , and S ( q ) 3 , calculated using Equation (35), and load q.
Applsci 11 12124 g009
Figure 10. The dependency curve between total settlement S ( q ) , calculated by Equation (35), and load q with various values of pre-overburden pressure σ P and foundation depth z.
Figure 10. The dependency curve between total settlement S ( q ) , calculated by Equation (35), and load q with various values of pre-overburden pressure σ P and foundation depth z.
Applsci 11 12124 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ter-Martirosyan, Z.G.; Ter-Martirosyan, A.Z.; Dam, H.H. Settlement and Bearing Capacity of Rectangular Footing in Reliance on the Pre-Overburden Pressure of Soil Foundation. Appl. Sci. 2021, 11, 12124. https://doi.org/10.3390/app112412124

AMA Style

Ter-Martirosyan ZG, Ter-Martirosyan AZ, Dam HH. Settlement and Bearing Capacity of Rectangular Footing in Reliance on the Pre-Overburden Pressure of Soil Foundation. Applied Sciences. 2021; 11(24):12124. https://doi.org/10.3390/app112412124

Chicago/Turabian Style

Ter-Martirosyan, Zaven G., Armen Z. Ter-Martirosyan, and Huu H. Dam. 2021. "Settlement and Bearing Capacity of Rectangular Footing in Reliance on the Pre-Overburden Pressure of Soil Foundation" Applied Sciences 11, no. 24: 12124. https://doi.org/10.3390/app112412124

APA Style

Ter-Martirosyan, Z. G., Ter-Martirosyan, A. Z., & Dam, H. H. (2021). Settlement and Bearing Capacity of Rectangular Footing in Reliance on the Pre-Overburden Pressure of Soil Foundation. Applied Sciences, 11(24), 12124. https://doi.org/10.3390/app112412124

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop