Next Article in Journal
Dynamic Monitoring of Steel Beam Stress Based on PMN-PT Sensor
Next Article in Special Issue
Numerical Study on Mechanical Characteristics of Tower Sections with Main Member Disconnection
Previous Article in Journal
Evaluation Method and the Influence of Visual Comfort of Ceramic Tiles in Indoor Environment—A Study Based on the Delphi and AHP
Previous Article in Special Issue
Examination of Damage Evolution in Slurry Masonry Schist Subjected to Biaxial Compressive Stresses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Settlement of a Pile Foundation Considering Linear and Rheological Properties of Soils

by
Armen Z. Ter-Martirosyan
*,
Lyubov Yu. Ermoshina
* and
George O. Anzhelo
Department of Soil Mechanics and Geotechnical Engineering, National Research Moscow State Civil Engineering University, 26, Yaroslavskoye Shosse, 129337 Moscow, Russia
*
Authors to whom correspondence should be addressed.
Buildings 2024, 14(9), 2830; https://doi.org/10.3390/buildings14092830
Submission received: 2 July 2024 / Revised: 2 September 2024 / Accepted: 5 September 2024 / Published: 8 September 2024
(This article belongs to the Special Issue Building Foundation Analysis: Soil–Structure Interaction)

Abstract

:
Despite numerous studies of single piles and practical experience with their application, methods for calculating settlements of pile foundations remain limited. The existing objective need for specialized methods of pile foundation settlement calculation that take into account the rheological properties of the base soils is becoming more and more important, especially in the construction of unique objects in complex ground conditions. When predicting the stress–strain state of the pile–raft-surrounding soil mass system, it is allowed to consider not the entire pile foundation as a whole, but only a part of it—the computational cell. In the present work, we have solved the problems of determining the strains of the computational cell consisting of the pile, the raft and the surrounding soil according to the column pile scheme and hanging pile scheme, on the basis of the Kelvin–Voigt rheological model, which is a model of a viscoelastic body consisting of parallel connected elements: Hooke’s elastic spring and Newtonian fluid. According to our results, we obtained graphs of the dependence of strains of the computational cell on time at different pile spacing and different values of coefficients of viscosity of the surrounding soil, and a formula for calculating the reduced modulus of deformation of the pile. The results of the present study can significantly improve the accuracy of calculations during construction on clayey soils with pronounced rheological properties and, as a result, increase the reliability of pile structures in general.

1. Introduction

For many years, the pile foundation has been one of the most popular types of foundations used in practice due to its reliability, speed of construction and the possibility of construction in almost any engineering–geological and climatic conditions. At the same time, the widespread use of pile foundations is due to the increasing number of stories and height of buildings, as well as the increasing weight of technological equipment.
When an external load is applied to a pile foundation, in addition to the settlements of the piles themselves, there is settlement of the soil in the inter-pile space, causing the friction forces developed on the lateral surfaces of the piles to not fully develop, unlike in single piles. When increasing the spacing between piles to 3–4 pile diameters, the development of friction forces on the lateral surface will only be partially realized.
It should be noted that the stress–strain state of the active zone in pile foundations significantly differs from the stress–strain state of the soil mass around a single pile, as in pile foundations there is a mutual influence of piles on each other. Therefore, it is not permissible to use the results of static tests of single piles and analytical methods to determine their settlements when calculating the settlements of a pile foundation.
Despite a large number of scientific works devoted to the study of single piles (Luga A.A. [1,2], Ogranovich A.B. [3], Golubkov V.N. [4,5], etc.), and the availability of many years of experience in the use of piles in the construction of various buildings and structures, for quite a long time there was no single engineering method for calculating the settlement of pile foundations and pile groups.
In the 1960s, the method by Egorov K.E. began to be used to calculate the settlement of pile groups [6], based on the calculation of a rectangular foundation placed on a natural base. The specific feature of this calculation method was that the pile foundation was considered as a massive foundation, and the place of application of the calculated pressure on the soil under the foundation bottom was dependent on the ratio of the length of the piles to the width of the foundation. If this ratio was less than or equal to one, then the point of application of the calculated pressure was taken as the surface located at the tip of the pile. If the ratio of pile length to foundation width exceeded unity, the design pressure was distributed to the surface at a depth equal to the foundation width [7].
At present, the prediction of the rate of settlement of foundations over time, formed mainly by weak clay soils, is quite relevant in the design of various types of foundations, including pile foundations. Since most settlements are uneven and can develop over years and even decades, it is necessary to know the rate of development of these settlements over time to ensure the normal operation of buildings and structures erected on weak clay soils.
Foundation settlement over time can be predicted by taking into account the long-term complex physical processes that occur during load transfer to soils that have rheological properties (creep properties). To describe the stress–strain state of the base foundation system, taking into account the viscoelastic and viscoelastic–plastic properties of the base soils, it is necessary to take into account many parameters, the values of which may change during the operation of the structure [8]. Based on the above, it is important to note that the time factor must be taken into account when calculating foundations on bases formed by soils possessing rheological properties.
Numerous experimental and theoretical studies aimed at studying the rheological properties of clay soils, presented in the works of Tertsagi K. [9], Gersevanov N.M. [10], Florin V.A. [11,12], Cytovich N.A. [13], Meschyan S.R. [14], Vyalov S.S. [15], Maslov N.N. [16], Ter-Martirosyan Z.G. [17] and many others, allowed us to develop methods for calculating foundations on a natural base over time.
Many researchers have studied the interaction of piles with clay soils, including Vyalov S.S. [15], Bartolomei, A.A. [18], Dalmatov B.I. and Lapshin F.K. [19], Ter-Martirosyan Z.G. [17], Omelchak I.M. [8] and many others.
Liu J. et al. [20] proposed a tri-linear softening model to describe the interaction between a single pile and the surrounding soil, and a tri-linear plastic model that describes the interaction between the pile and the underlying soil. Based on the mechanical properties of soils and the theory of elasticity, the authors developed equations that describe the dependence of the settlement of a single pile on the vertical load in layered soils. In performing the calculations, the pile was considered as a flexible shaft with a constant cross-section. The authors of the study have demonstrated the effectiveness of the proposed calculation method for a single pile and provided recommendations on using the proposed models in different soil layers for a correct analysis of the relationship between the pile settlement and the applied vertical load.
The study [21] presents an analytical solution for determining the stress–strain state of a soil mass interacting with underground structures, including with piles.
The interaction of a single long pile with a multilayer soil mass, taking into account the elastic and rheological properties of the soil mass, was investigated in [22,23]. In [22], the solution of the problem by analytical method taking into account linear and nonlinear properties of soils is presented. For describing nonlinear shear deformations, the elasto-plastic Timoshenko S.P. model was used. The authors performed a comparison of the results of the analytical solution to the elastic problem with the results of numerical modeling in the Plaxis 3D software package, and also obtained an expression for determining the reduced shear modulus for a multilayer soil mass. The study [23] considers the solution to the problem of determining the displacements of a long pile in a surrounding multilayer base. It was found that the rate of stress change under the pile heel depends on the viscosity of the soil, while the rheological hardening coefficient significantly influences the time of pressure stabilization under the pile heel and the time of pile settlement stabilization.
Based on the modified Burgers model, in which the spring is replaced by a hyperbolic model, a nonlinear approach to calculate the settlement over time of a vertically loaded single pile as well as a group of piles in a layered soft soil was developed in [24]. In this work, parametric studies were carried out to investigate the influence of the parameters included in the calculation on the time-dependent settlement of a single pile. Based on the conducted research, it was shown that the theoretical results agree well with the measurement results, indicating a sufficiently high accuracy of the approach proposed in the work.
For calculating the long-term settlement of a single pile embedded in viscoelastic base soils, a three-dimensional viscoelastic model with a fractional-order derivative was proposed by Li, Zhang and Liu to describe the rheological behavior of the soils around the pile [25].
A simple approach to analyzing the elasto-plastic behavior of a single pile in layered soils, which allows for a quick assessment of the settlement of a single pile as well as a group of piles, is presented in [26]. The proposed approach is based on the application of two models, one of which describes the nonlinear relationship between the shear stresses occurring on the lateral surface of the pile and the displacement of the soil around the pile. The second model is based on a bilinear hardening model to simulate the relationship between the settlement of the pile base and the resistance of the pile tip.
In [27], a model was proposed for determining the settlement of a vertically loaded single hanging pile embedded in layered base soils.
To evaluate the behavior of a single vertical pile over time, embedded in soil with elasto-plastic properties, the authors of [28] conducted an analysis using the finite element method in a two-dimensional setting. For modeling the soil around the pile, a linear elastic perfectly plastic Mohr–Coulomb model was used, and for modeling the pile material, a linear elastic model was employed. Based on the calculations performed, it was found that a pile in non-cohesive soil has greater resistance under rapid loading than under long-term loading, whereas the opposite results were obtained for a pile in cohesive soil. The geotechnical model presented in this work allows the determination of lateral deformations as well as lateral soil stress and its variation over time.
An analytical method for evaluating the behavior of a single tapered pile and a group of piles from applied static axial compressive loads was proposed in [29]. Based on the proposed method, an iterative computer program was developed to calculate the settlement and bearing capacity of a single tapered pile.
An analytical approach for predicting the settlement of vertically loaded pile and pile groups was proposed by Xia and Zou [30]. A piecewise function was applied to investigate the nonlinear relationship between surface friction (friction on the side surface of the pile) and relative displacement of the pile–soil system. In this paper, a comparative analysis of the results obtained using the proposed calculation methodology with field test data and results obtained by other researchers has been carried out.
In the study by Long Li and Yousheng Deng [31], an algorithm is proposed to calculate the settlement of a pile foundation in a soil base having linear–viscous–elastic properties. The calculation algorithm presented in this paper is based on the analysis of the change in linear viscoelastic settlement of the soil depending on the geometric parameters of the pile (pile spacing, pile length-to-diameter ratio) and the deformation characteristics of the soil (elastic modulus and Poisson’s ratio). The authors of the work performed a comparative analysis of the obtained results of the study with the results of field tests of piles.
A methodology based on the use of stress coefficients is presented by Gendy [32] to determine the settlement of hanging piles located in clay soil. The stress coefficients were obtained from the Mindlin equation by eliminating the Poisson’s ratio from this equation. Based on the results of the conducted research, equations were derived that allow predicting the settlement of a single pile, a pile group and a pile–raft foundation on clay soils.
It should be emphasized that taking into account the rheological properties of the base soils is important not only in static calculations, but also in dynamic calculations of foundations. This statement is supported by a number of studies and scientific papers that demonstrate the influence of rheological characteristics on the behavior of foundations subjected to dynamic loads [33,34]. The correct interpretation of these properties allows for a more accurate assessment of the stability and durability of structures, especially under variable and dynamic loads such as seismic impacts or vibrations induced by machinery.
Based on the analysis of the scientific and technical literature on the research topic, the following main conclusions can be drawn:
-
Most scientific papers focus on determining the settlement and bearing capacity of a single pile, rather than the pile foundation as a whole;
-
When determining the settlement of pile foundations, most scientific works do not take into account the rheological properties of the base soils;
-
The existing scientific papers present limited data on the influence of one of the most important rheological parameters—the coefficient of soil viscosity—on the settlement of pile foundations in general;
-
A number of papers present calculations that take into account the rheological properties of base soils in relation to dynamic calculations.
Based on the above, the purpose of this paper is to determine the settlement of a pile foundation, taking into account the linear and rheological properties of the base soils based on the Kelvin–Voigt viscoelastic rheological model of the viscoelastic body according to the column pile and hanging pile schemes.

2. Materials and Methods

It is known that when a uniformly distributed load is applied to a pile foundation, it will be distributed among the piles, the raft and the surrounding soil mass located in the inter-pile space. Moreover, this distribution significantly depends on the diameter, length and spacing of the piles, as well as on the physical–mechanical and rheological properties of the soils and the pile material. With such a load transfer mechanism (using piles and a raft), a complex, heterogeneous stress–strain state is formed in the soil mass, which can change in space and time, especially when considering the rheological properties of the soils.
The important and determining design parameters when calculating a pile foundation are the settlement and bearing capacity of both the pile foundation as a whole and the individual piles within it (in the case of uneven distribution of piles in the pile foundation).
When predicting the stress–strain state of the pile–raft-surrounding soil mass system, it is allowed to consider not the entire pile foundation as a whole, but only a part of it—the computational cell.
The computational cell is a thick-walled soil cylinder of limited dimensions L ,   2 b 1 , accommodating the pile, the raft and the surrounding soil mass. The geometric dimensions of the computational cell are chosen based on the load on the foundation and the engineering–geological conditions of the construction site.
Figure 1 shows the design scheme of the pile foundation with the delineation of the boundaries of the computational cell L ,   2 b 1 (dashed line).
When calculating piles as part of a pile foundation, the main thing is the quantitative assessment of the settlement of the pile foundation as a whole, which is largely determined by the magnitude of the load on the raft and the correct assessment of the stress–strain state of the pile–raft-surrounding soil mass system.
Depending on the ratio of the stiffness of the surrounding soil ( G 1 ) and the underlying soil ( G 2 ), the calculation of the stress–strain state of the pile foundation should be carried out either according to the column pile scheme when G 2 G 1 , or according to the hanging pile scheme when G 2 > G 1 .
According to SP 24.13330.2021 [35], a column pile is a pile resting on rocky or weakly deformed soils and transmitting the load to the base mainly through the heel.
A hanging pile is a pile that transmits the load to the base through the lateral surface and the heel [35].
Figure 2 shows the design scheme of the interaction of the pile as part of the pile foundation with the surrounding and underlying soils and the raft according to two schemes.
Since, in the following, all main calculations will be based on the Kelvin–Voigt rheological model, the main features of this model are presented below for a better understanding of the results obtained.
The Kelvin–Voigt model is a rheological model of a viscoelastic body consisting of parallel connected elements: an elastic Hooke spring (H) and a Newtonian fluid (N) (Figure 3). When two elements are connected in parallel, the strain of each element will be the same and the total stress will be defined as the sum of the stresses on each element.
The dependence of strains ( ε ) on stresses ( σ ) according to the Kelvin–Voigt rheological model will be as follows:
σ = ε · E + ε ˙ · η
where E is the modulus of deformation; ε ˙ is the strain rate; and η is the coefficient of soil viscosity.
Note that the Kelvin–Voigt model is designed to analyze the mechanical properties of continuous media. In the context of our study, the ground can be classified as a continuous medium, which allows us to use this model to describe its behavior. This is confirmed by the ability of the model to take into account both elastic and viscous properties of the material, which makes it convenient for analyzing the stress–strain state of the soil mass under various loads.
A more detailed description of the Kelvin–Voigt rheological model can be found in the following references: [17,36,37].

3. Results

3.1. Solution of the Problem in a Linear Formulation According to the Column Pile Scheme

Let a long pile with a given diameter of 2 a 1 and length l 1 in a pile foundation be embedded in the soil mass, with its lower end resting on an underlying relatively dense soil layer ( G 2 G 1 ). The computational domain of the massif containing the pile and the foundation is a thick-walled soil cylinder of limited dimensions L , 2 b 1 .
When a static uniformly distributed load ( σ N = c o n s t ) is applied to the raft, it will be distributed between the piles and the surrounding soil within each computational cell. Considering the compressive deformation of the pile shaft and the surrounding soil mass (in accordance with the problem conditions), we can write the following equilibrium equation:
σ N = σ p · ω + σ s · ( 1 ω ) ,
where σ N is the uniformly distributed load applied to the raft; σ p is the stress acting in the pile shaft at the level of the pile head; σ s is the stress acting in the soil mass at the contact with the raft; ω is the dimensionless coefficient equal to the ratio of the cross-sectional area of the pile shaft to the cross-sectional area of the entire computational cell, i.e., ω = a 1 2 / b 1 2 ; and a 1 and b 1 are the radius of the pile and the radius of the computational cell, respectively.
Let us also write down the condition of equality of settlements of the raft, the pile and the surrounding soil mass:
S N = S p = S s
Equation (3) can be represented in the following form:
m r · σ N = m p · σ p = m s · σ s
where m p and m s are the coefficients of relative compressibility of the pile and the soil, respectively.
The joint solution of Equations (2) and (4) will determine σ p and σ s , respectively:
σ p = σ N · E p E p · ω + E s · 1 ω ;   σ s = σ N · E s E p · ω + E s · 1 ω
The reduced modulus of deformation of the computational soil cell can be found from condition (2), i.e.:
E r = E p · ω + E s · ( 1 ω )
From Formula (6), it follows that as ω increases, E r also increases, and conversely. Consequently, other things being equal, by changing the ratio ω = a 1 2 / b 1 2 , it is possible to adjust both the value of the reduced modulus of deformation of the computational cell E r and the calculated value of the settlement of the cell as a whole.
Based on the above, the settlement of the pile foundation can be determined according to the following formula:
S p r = σ N E r · 0.8 · l 1
where l 1 is the length of the pile.
When calculating by Formula (7), the settlement of the underlying layer should be neglected.
Based on the known dependency S = ε · l , the solution to the settlement problem of a pile foundation cell can be obtained when the length ( l ) and strain ( ε ) are known. Since length is an input parameter, the following problems will be considered in order to determine the strains of the computational cell.

3.2. The Stress–Strain State of the Cell in an Elasto-Viscous Formulation According to the Column Pile Scheme Based on the Kelvin–Voigt Model

Let us write down the equilibrium condition characterizing the process of development of settlement of the raft over time, assuming in the first approximation that it coincides with the equilibrium Equation (2), then we obtain:
σ N = σ p ( t ) · ω + σ s ( t ) · ( 1 ω )
The stresses acting in the pile shaft at the level of the pile head ( σ p ) and the stresses acting in the soil mass at the contact with the raft ( σ s ) can be represented by the following relationships:
σ p t = ε p · E p ;   σ s t = ε s · E s + ε s ˙ · η s
where ε p is the strain of the pile; E p is the modulus of deformation of the pile; ε s is the strain of the surrounding soil; E s is the modulus of deformation of the surrounding soil; ε s ˙ is the rate of development of strains of the surrounding soil; and η s is the coefficient of viscosity of the surrounding soil.
In Figure 4, the calculated model of the interaction of a pile in a pile foundation and the surrounding soil mass under compressive deformation is presented based on the Kelvin–Voigt model.
Substituting into the equilibrium condition (8) the dependences obtained in (9), we obtain:
σ N = ε p · E p · ω + ( ε s · E s + ε s ˙ · η s ) · ( 1 ω )
Based on the condition of equality of strains of the pile and the surrounding soil ε p = ε s = ε , we obtain:
σ N = ε · ( E p · ω + E s · 1 ω ) + ε ˙ · η s · ( 1 ω )
where ε = ε p = ε s , ε ˙ = ε s ˙ .
Let us represent the differential Equation (11) in the following form:
ε ˙ + ε · P = Q
P = E p · ω + E s · 1 ω η s · ( 1 ω ) ; Q = σ N η s · 1 ω
The solution of the differential Equation (12) is known [38] and has the following form:
ε ( t ) = e P d t [ Q · e P d t d t + C ]
By performing the corresponding transformations in (14), we obtain the following equation:
ε ( t ) = Q P + C · e P · t
We determine the integration constant from the initial condition when ε 0 = 0 at t = 0 :
C = Q P
Substituting the integration constant C rom Equation (16) into Equation (15), we obtain:
ε ( t ) = Q P · ( 1 e P · t )
It follows from Equation (17) that:
-
at t 0   ε 0 = 0 ;
-
at t   ε = Q P = σ N E p · ω + E s · 1 ω .
Table 1 shows the characteristics of the pile and surrounding soil used in the calculation using Formula (17).
E p a 1 2 b 1 2 b 2 2 b 3 2 b 4 E s η s 1 η s 2 η s 3 Figure 5, Figure 6 and Figure 7 show the graphs of the dependence of strains of the computational cell on time ε t , obtained with different pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 and with different coefficients of viscosity of the surrounding soil η r 1 < η r 2 < η r 3 based on Formula (17).
Figure 8 shows the graphs of the dependence of strains of the computational cell on time ε t , obtained with different coefficients of viscosity of the surrounding soil η s 1 < η s 2 < η s 3 and with pile spacing 2 b 2 , based on Formula (17).

3.3. The Stress–Strain State of the Cell in an Elasto-Viscous Formulation According to the Hanging Pile Scheme Based on the Kelvin–Voigt Model

The analysis of the solution of the problem in the elastic–viscous formulation according to the column pile scheme showed that if we replace the modulus of deformation of the pile with the reduced modulus of deformation, taking into account the local settlement under the pile heel, we can find the solution of the problem according to the hanging pile scheme. For this purpose, it is sufficient to change the stiffness of the pile in such a way that the following condition is satisfied:
σ p E p · l 1 = σ p E p · l 1 + σ p · π · a 1 · 1 ν 2 2 2 · E 2
E p = E p 1 + E p · π · a 1 · 1 ν 2 2 2 · E 2 · l 1
where E p is the reduced modulus of deformation of the pile; E p is the modulus of deformation of the pile; E 2 is the modulus of deformation of the underlying soil mass; ν 2 is the Poisson’s ratio of the underlying soil mass; a 1 is the radius of the pile; and l 1 is the length of the pile.
Table 2 shows the characteristics of the pile and surrounding soil used in the calculation using Formula (17), taking into account Formula (19).
Based on the above, we can conclude that by replacing the modulus of deformation of the pile E p with the reduced modulus of deformation E p , in the rheological Equation (11), we obtain the stress–strain state of the cell according to the hanging pile scheme (Figure 9a).
As a result, we can predict the strains of the computational cell both using the column pile scheme ε 1 t and the hanging pile scheme ε 2 t (Figure 10) based on Formula (17), taking into account (19).

4. Discussion

When the load is transferred to the pile foundation, in addition to the settlement of the piles themselves, the settlement of the surrounding soil located in the inter-pile space occurs. When considering the stress–strain state of the active zone in pile foundations, it should be noted that it differs significantly from the stress–strain state of the surrounding soil mass around a single pile, because in pile foundations there is a mutual influence of the piles on each other, as well as a time distribution of the external load between the raft, piles and the surrounding soil. However, despite this, numerous studies [20,21,22,23,25,26,27,28,29] have focused on determining the settlement of a single pile rather than the pile foundation as a whole.
The formation of a complex and heterogeneous stress–strain state of the system pile–raft-surrounding soil mass system depends on:
-
geometric characteristics (cross-sectional and longitudinal profile shape, length) and pile material, as well as pile spacing;
-
engineering–geological and hydrogeological conditions of the construction site;
-
mechanical characteristics of the surrounding soil, which include deformation ( E , ν ) , strength ( c , φ ) and rheological ( η ) properties.
Indeed, the factors presented above have a significant influence on the formation of the stress–strain state of the active zone in the calculation of the pile foundation as a whole, which has been reviewed and noted by Long Li and Yousheng Deng in [31]. This study investigated the change in linear viscoelastic settlement of soil depending on the pile spacing, the ratio of pile length to diameter, soil elasticity modulus and Poisson’s ratio. However, unfortunately, this study did not pay sufficient attention to the influence of the rheological properties of clay soils (including soil viscosity) on the obtained values of settlement of the pile foundation over time.
Taking into account the rheological properties of clay soils is necessary in order to predict over time the final settlements of pile foundations constructed on clay soils, to prevent in advance the negative phenomena associated with creep of clay soils and, if necessary, to optimize the design solutions of foundations.
In the present study, the problem of the interaction between the pile, the raft and the surrounding soil mass was solved in a linear formulation using the column pile scheme. As a result, Formula (6) was obtained for calculating the reduced deformation modulus of the computational soil cell, as well as Formula (7) for calculating the settlement of the pile foundation.
To take into account the rheological properties of the base soils, a problem was solved using an elastic–viscous formulation in the column pile scheme and in the hanging pile scheme based on the well-known Kelvin–Voigt rheological model. This model allows us to describe creep of clay soils over time quite well because it includes a viscosity parameter.
Let us consider in more detail the main advantages of this model for a better understanding of the results obtained in this study.
  • The structure of the model: as mentioned earlier in this paper, the Kelvin–Voigt model consists of an elastic element (spring) and a viscous element (Newtonian fluid) connected in parallel. This allows this model to describe the processes of elastic deformation and viscous flow, which is characteristic of the behavior of clayey soils under the action of a constant load. We also note that the Kelvin–Voigt model takes into account both the rate and magnitude of deformation, describing the behavior of the soil under the action of constant or temporary loads. For example, when the load is suddenly increased, the elastic deformation will occur immediately, while the viscous part will determine how the soil will continue to deform over time.
  • The ability to take creep into account: in contrast to simple elastic models, the Kelvin–Voigt model allows us to take into account the creep phenomenon inherent in clayey soils. This means that when a constant load is applied, deformations will slowly continue to develop in the clayey soil over a long period of time. These deformations can be quantified using this model, which represents an important aspect in the design of various foundations and enclosing structures, since the final settlement can significantly affect the project and operation of the structure.
  • Taking into account the time factor: the Kelvin–Voigt model makes it possible to estimate deformations of soil with rheological properties both in the short-term (a few days) and long-term (several decades) time range, which is especially important in the design of foundations and various enclosing structures.
  • Flexibility and adaptability: the Kelvin–Voigt model can be used to calculate both single-layer and multilayer bases.
  • Ease of use: the mathematical description of the model is relatively simple, which will allow its future implementation in numerical calculation methods.
  • Empirical validity: the model agrees well with experimental data obtained from tests of clayey soils. This makes it a reliable tool for engineering analysis.
Based on the above, we conclude that the Kelvin–Voigt model is a powerful tool for calculating and predicting long-term deformations of soil bases with rheological properties. The model’s ability to take into account both elastic and viscous components allows engineers to more accurately estimate long-term deformations and make informed decisions in the design of various foundations and enclosing structures. Thus, the application of this model provides an opportunity for more reliable and efficient design and minimization of potential risks for the objects under construction.
Based on the solution to the problem using the column pile scheme based on the Kelvin–Voigt model, the following results were obtained:
-
It is shown that the vertical strains of the computational cell also increase with increasing pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 (Figure 5, Figure 6 and Figure 7), because with increasing pile spacing, the load applied to the raft starts to be supported more by the surrounding soil mass located in the inter-pile space, which has a modulus of deformation much lower than that of the pile itself;
-
At low values of the coefficient of viscosity of the surrounding soil ( η s 1 = 10 9 Poise) (Figure 5), comparable to the liquid medium, there is a peak value of vertical strains exceeding the residual vertical strains, which is not observed in more viscous media ( η s 2 = 10 11 Poise, η s 3 = 10 13 Poise) (Figure 6 and Figure 7). This phenomenon indicates significant differences in soil behavior depending on its rheological properties;
-
In accordance with the character of location of the graphs of dependence of strains of the computational cell on time ε t presented in Figure 8, it can be seen that the time of stabilization of vertical strains of the computational cell is directly proportional to the value of the coefficient of viscosity of the surrounding soil, i.e., the higher the value of the coefficient of viscosity, the more time will be needed in order to stabilize the vertical strains of the computational cell from the applied load. However, it should be noted that changing the coefficient of viscosity of the surrounding soil, both increasing and decreasing, will affect only the time of stabilization of vertical strains of the computational cell, and will not change their value. The results obtained confirm the importance of considering time factors in assessing the stability and durability of pile structures, especially under variable loads.
Based on the solution to the problem using the hanging pile scheme based on the Kelvin–Voigt model, the following results were obtained:
-
Formula (19) for calculating the reduced modulus of deformation of the pile ( E p ) was obtained, which makes it possible to take into account the strains of the underlying soil layer and considerably simplify the procedure for calculating the stress–strain state of the computational cell according to the hanging pile scheme. The resulting equation takes into account both the stiffness of the pile itself ( E p ) and the stiffness of the underlying soil mass under the pile heel ( E 2 );
-
Building graphs of dependence of vertical strains of the computational cell over time according to the column pile scheme ε 1 t and according to the hanging pile scheme ε 2 t , presented in Figure 10, showed that they are significantly different. This difference shows that taking into account the local settlement under the pile heel in the design diagram leads to an increase in settlement, which corresponds to the “hanging pile” design diagram (Figure 9a);
-
When the modulus of deformation of the surrounding soil ( E s ), the modulus of deformation of the underlying soil mass ( E 2 ) and the Poisson’s ratio of the underlying soil mass ( ν 2 ) decrease, the vertical strains of the computational cell in the hanging pile scheme will increase, and conversely, when these parameters increase, the vertical strains of the computational cell will decrease;
-
Change in the coefficient of viscosity of the surrounding soil ( η s ) does not influence the value of vertical strains of the computational cell both for the column pile and hanging pile schemes. However, this parameter significantly influences the time of stabilization of vertical strains, which can be important in the design of foundations, especially under conditions of dynamic loads on building structures.
The analysis of the obtained results allows us to conclude that when designing pile foundations, it is necessary to take an individual approach to each specific construction case, taking into account the type of soil, its rheological properties and its mechanical behavior. The use of the proposed Kelvin–Voigt viscoelastic body model and the obtained formula for calculating the reduced modulus of deformation can significantly increase the accuracy of calculations and, consequently, the reliability of pile structures in general. Importantly, further research in this area will help to optimize existing design standards, leading to safer and more efficient solutions in engineering practice.

5. Conclusions

The acceleration of scientific and technical progress in foundation engineering is inextricably linked to the development of new calculation methods that allow for a more reliable description of the interaction between different types of foundations and base soils during the transfer of loads from buildings, structures and technological equipment [8]. Numerous observations of the settlement of various buildings and structures over the years have shown that a promising direction of research into the interaction of various foundations and soil bases is to take into account the rheological properties of soils (consideration of viscous deformation under the action of constant and variable loads).
Taking into account the fact that many buildings and structures are constructed on bases composed of weak clay soils, the settlement of which, as a rule, is uneven and can develop over years and even decades, it is necessary to consider the rheological properties of clay soils when calculating and designing various types of foundations, including pile foundations. The widespread use of pile foundations is due to their reliability, speed of construction and possibility of construction in almost any engineering–geological and climatic conditions.
The following main results were obtained, meeting the purpose of this study as outlined in the introduction:
-
The dependence (17) was obtained, which allows us to determine the vertical strains of the computational cell (and, at known pile length, it is possible to obtain the settlement of the pile foundation as a whole) both for the column pile scheme and for the hanging pile scheme;
-
The Formula (19) for calculation of the reduced modulus of deformation of the pile ( E p ) was obtained, which allows us to take into account the local settlement of the underlying soil under the heel of the pile and thus to pass from the solution of the problem according to the column pile scheme to the solution of the problem according to the hanging pile scheme (Figure 10 clearly shows the difference in the values of vertical strains of the computational cell obtained when solving the problem using the column pile scheme ε 1 t and the hanging pile scheme, taking into account the reduced modulus of deformation of the pile ε 2 t );
-
It was found that the coefficient of viscosity of the surrounding soil has no influence on the values of vertical strains of the computational cell. This parameter only influences the stabilization time of the vertical strains, i.e., how quickly the vertical strains stabilize. Figure 8 effectively demonstrates how the stabilization time of vertical strains of the computational cell varies significantly with different values of the coefficients of viscosity of the surrounding soil η s 1 < η s 2 < η s 3 .
The practical significance of this study consists in solving the actual problem of determining the settlement of a pile foundation, taking into account the linear and rheological properties of soils under the column pile scheme and under the hanging pile scheme on the basis of the known viscoelastic Kelvin–Voigt model.
As future research, it is planned to further develop methods of quantitative assessment of the stress–strain state of the base-foundation system, taking into account rheological properties of the clay soil in order to exclude excessive development of settlements, tilts and displacements of various buildings and structures, as well as to optimize the adopted design solutions of foundations.

Author Contributions

Conceptualization, methodology, A.Z.T.-M.; software, validation, L.Y.E. and G.O.A.; formal analysis, writing—review and editing, all authors; investigation, A.Z.T.-M., L.Y.E. and G.O.A., supervision, project administration, A.Z.T.-M.; resources, A.Z.T.-M. and G.O.A.; visualization, L.Y.E. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the National Research Moscow State University of Civil Engineering (grant for fundamental and applied scientific research, project No. 02-392/130).

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Luga, A.A. Methodical Instructions on Calculation of Settlements of Single Piles; TsNIIS: Moscow, Russia, 1963. [Google Scholar]
  2. Luga, A.A. Settlement of Pile Foundations Working in Conditions of Absence of Cluster Effect VNITS. Investigations of Bearing Capacity of Foundation Bases; Transizdat: Moscow, Russia, 1965. [Google Scholar]
  3. Ogranovich, A.B. To the question of determining the settlement of a single pile. Bases Found. Soil Mech. 1963, 1. [Google Scholar]
  4. Golubkov, V.N.; Tugaenko, Y.F.; Khutoryansky, B.O. Investigation of the process of deformation zone formation in single pile foundations. Bases Found. 1971, 4, 9–13. [Google Scholar]
  5. Golubkov, V.N. On the nature of joint work of piles and soil. Bases Found. 1972, 5, 29–35. [Google Scholar]
  6. Egorov, K.E. To the question of permissible settlement of foundations of structures. Bases Found. 1952, 18. [Google Scholar]
  7. Bartolomei, A.A. Calculation of Settlements of Strip Pile Foundations; Stroyizdat: Moscow, Russia, 1972; 128p. [Google Scholar]
  8. Omelchak, I.M. Fundamentals of Calculation of Pile Foundations Taking into Account Rheological Properties of Foundation Soils. Ph.D. Thesis, Perm State Technical University, Perm, Russia, 2003. [Google Scholar]
  9. Tertsagi, K. Theory of Soil Mechanics; Gostroyizdat: Moscow, Russia, 1961; 507p. [Google Scholar]
  10. Gersevanov, N.M. Theoretical Foundations of Soil Mechanics and Their Practical Applications; Strojizdat: Moscow, Russia, 1948; 248p. [Google Scholar]
  11. Florin, V.A. Fundamentals of Soil Mechanics. Volume 1; Strojizdat: Leningrad, Russia, 1959; 356p. [Google Scholar]
  12. Florin, V.A. Fundamentals of Soil Mechanics. Volume 2; Gosstroyizdat: Leningrad, Russia, 1961; 543p. [Google Scholar]
  13. Cytovich, N.A. Soil Mechanics; State Publishing House of Literature on Construction, Architecture, and Building Materials: Moscow, Russia, 1963; 437p. [Google Scholar]
  14. Meschyan, S.R. Creep of Clay Soils; Publishing House of the Academy of Sciences of the Armenian SSR: Yerevan, Armenia, 1967; 320p. [Google Scholar]
  15. Vyalov, S.S. Rheological Foundations of Soil Mechanics: Textbook for Construction Universities; Vysshaya Shkola: Moscow, Russia, 1978; 447p. [Google Scholar]
  16. Maslov, N.N. Long-Term Stability and Deformation of Displacement Retaining Walls; Energiya: Moscow, Russia, 1968; 467p. [Google Scholar]
  17. Ter-Martirosyan, Z.G. Rheological Parameters of Soils and Foundation Calculations; Strojizdat: Moscow, Russia, 1990; 200p. [Google Scholar]
  18. Bartolomei, A.A.; Omelchak, I.M.; Yushkov, B.S. Pile Foundation Settlement Prediction; Stroyizdat: Moscow, Russia, 1994; 384p. [Google Scholar]
  19. Dalmatov, B.I.; Lapshin, F.K.; Rossikhin, Y.V. Design of Pile Foundations in Conditions of Weak Soils; Stroyizdat: Leningrad, Russia, 1975; 240p. [Google Scholar]
  20. Liu, J.; Xiao, H.B.; Tang, J.; Li, Q.S. Analysis of load-transfer of single pile in layered soil. Comput. Geotech. 2004, 31, 127–135. [Google Scholar] [CrossRef]
  21. Ter-Martirosyan, Z.G.; Ter-Martirosyan, A.Z. Foundation and Pile Interaction with Swelling Soil. Soil Mech. Found. Eng. 2018, 55, 291–297. [Google Scholar] [CrossRef]
  22. Ter-Martirosyan, Z.G.; Akuletskii, A.S. Interaction between a long pile and multi-layer underlying soils. Vestn. MGSU 2021, 16, 168–175. [Google Scholar] [CrossRef]
  23. Ter-Martirosyan, Z.G.; Akuletskii, A.S. Interaction between a long pile and the multilayered soil body with account for elastic and rheological properties as well as soil stabilization. Vestn. MGSU 2021, 16, 608–614. [Google Scholar] [CrossRef]
  24. Feng, S.; Li, X.; Jiang, F.; Lei, L.; Chen, Z. A nonlinear approach for time-dependet settlement analysis of a single pile and pile groups. Soil Mech. Found. Eng. 2017, 54, 7–16. [Google Scholar] [CrossRef]
  25. Li, X.-M.; Zhang, Q.-Q.; Liu, S.-W. Semianalytical solution for long-term settlement of a single pile embedded in fractional derivative viscoelastic soils. Int. J. Geomech. 2021, 21, 04020246. [Google Scholar] [CrossRef]
  26. Zhang, Q.-Q.; Zhang, Z.-M.; He, J.-Y. A simplified approach for settlement analysis of single pile and pile groups considering interaction between identical piles in multilayered soils. Comput. Geotech. 2010, 37, 969–976. [Google Scholar] [CrossRef]
  27. Liu, Q.-J.; Wang, J.-B.; Ma, J.-J.; Gao, W.-H. Vertically-loaded single floating pile in layered soils by thin annulus element method. Chin. J. Geotech. Eng. 2019, 41, 748–754. [Google Scholar] [CrossRef]
  28. Abbas, M.I.; Chik, H.Z.; Taha, R.M.; Shafiqu, S.M.Q. Time-dependent lateral response of pile embedded in elasto-plastic soil. J. Cent. South Univ. Technol. (Engl. Ed.) 2010, 17, 372–380. [Google Scholar] [CrossRef]
  29. Li, Y.; Li, W. Load-displacement behaviour of tapered piles: Theoretical modelling and analysis. Geomech. Eng. 2021, 26, 1–11. [Google Scholar] [CrossRef]
  30. Xia, Z.-Q.; Zou, J.-F. Simplified approach for settlement analysis of vertically loaded pile. J. Eng. Mech. 2017, 143, 04017124. [Google Scholar] [CrossRef]
  31. Li, L.; Deng, Y. Analysis of Settlement of Group Pile Foundation in Linear Viscoelastic Soil. Adv. Civ. Eng. 2023, 2023, 3207304. [Google Scholar] [CrossRef]
  32. El Gendy, M. Deriving equations for analyzing friction piles in clay soil. Ain Shams Eng. J. 2007, 42, 1–27. [Google Scholar]
  33. Chen, Q.; Li, Z.; Huang, X.; Zhang, C. Study of rheological-mechanical properties and vibration mechanics bandgap of row pile foundation. Eur. J. Comput. Mech. 2023, 32, 31–52. [Google Scholar] [CrossRef]
  34. Ren, X.; Wang, L. Study on the time domain semi analytical method for horizontal vibration of pile in saturated clay. Appl. Sci. 2024, 14, 778. [Google Scholar] [CrossRef]
  35. SP 24.13330.2021 «Pile Foundations»; Minstroy of Russia: Moscow, Russia, 2021; p. 113.
  36. Voigt, W. Ueber die Innere Reibung der Festen Körper, Insbesondere der Krystalle; Abhandlungen der Königlichen Gesellschaft von Wissenschaften zu Göttingen: German, 1890; pp. 3–47. [Google Scholar]
  37. Bulíček, M.; Málek, J.; Rajagopal, K.R. On Kelvin-Voigt model and its generalizations. Evol. Equ. Control Theory 2012, 1, 17–42. [Google Scholar] [CrossRef]
  38. Bronshtein, I.N.; Semendyaev, K.A. Handbook of Mathematics; State Publishing House of Physical and Mathematical Literature: Moscow, Russia, 1962; 608p. [Google Scholar]
Figure 1. Delineation of the computational cell boundaries in a pile foundation.
Figure 1. Delineation of the computational cell boundaries in a pile foundation.
Buildings 14 02830 g001
Figure 2. Design scheme of interaction of a pile as part of the pile foundation with the surrounding and underlying soils and the raft according to the following schemes: (a) hanging pile; (b) column pile: 1—pile; 2—raft; 3—surrounding soil mass; 4—underlying soil mass.
Figure 2. Design scheme of interaction of a pile as part of the pile foundation with the surrounding and underlying soils and the raft according to the following schemes: (a) hanging pile; (b) column pile: 1—pile; 2—raft; 3—surrounding soil mass; 4—underlying soil mass.
Buildings 14 02830 g002
Figure 3. The Kelvin–Voigt rheological model.
Figure 3. The Kelvin–Voigt rheological model.
Buildings 14 02830 g003
Figure 4. The calculated model of the interaction of a pile in a pile foundation and the surrounding soil mass under compressive deformation based on the Kelvin–Voigt model.
Figure 4. The calculated model of the interaction of a pile in a pile foundation and the surrounding soil mass under compressive deformation based on the Kelvin–Voigt model.
Buildings 14 02830 g004
Figure 5. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 and with a coefficient of viscosity of the surrounding soil η s 1 .
Figure 5. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 and with a coefficient of viscosity of the surrounding soil η s 1 .
Buildings 14 02830 g005
Figure 6. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 and with a coefficient of viscosity of the surrounding soil η s 2 .
Figure 6. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 and with a coefficient of viscosity of the surrounding soil η s 2 .
Buildings 14 02830 g006
Figure 7. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 and with a coefficient of viscosity of the surrounding soil η s 3 .
Figure 7. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different pile spacing 2 b 1 < 2 b 2 < 2 b 3 < 2 b 4 and with a coefficient of viscosity of the surrounding soil η s 3 .
Buildings 14 02830 g007
Figure 8. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different coefficients of viscosity of the surrounding soil η s 1 < η s 2 < η s 3 and with pile spacing 2 b 2 .
Figure 8. Graphs of the dependence of strains of the computational cell on time ε t , obtained with different coefficients of viscosity of the surrounding soil η s 1 < η s 2 < η s 3 and with pile spacing 2 b 2 .
Buildings 14 02830 g008
Figure 9. Design diagram for determining the settlement of: (a) hanging pile S 2 ; (b) column pile S 1 .
Figure 9. Design diagram for determining the settlement of: (a) hanging pile S 2 ; (b) column pile S 1 .
Buildings 14 02830 g009
Figure 10. Graphs of dependence of strains of the computational cell on time by the column pile scheme ε 1 t and by the hanging pile scheme ε 2 t .
Figure 10. Graphs of dependence of strains of the computational cell on time by the column pile scheme ε 1 t and by the hanging pile scheme ε 2 t .
Buildings 14 02830 g010
Table 1. Characteristics of the pile and surrounding soil (according to the column pile scheme).
Table 1. Characteristics of the pile and surrounding soil (according to the column pile scheme).
Name of the ParameterSymbolValueUnits
Characteristics of the pile
Modulus of deformation of the pile E p 3 × 107kPa
Radius of the pile a 1 0.30m
Pile spacing * 2 b 1 1.80m
2 b 2 2.40m
2 b 3 3.00m
2 b 4 3.60m
Characteristics of the surrounding soil
Modulus of deformation of the surrounding soil E s 2 × 104kPa
Coefficient of viscosity of the surrounding soil η s 1 1 × 109Poise
η s 2 1 × 1011Poise
η s 3 1 × 1013Poise
*—in the present study, the pile spacing is equal to the diameter of the computational cell (Figure 1).
Table 2. Characteristics of the pile and surrounding soil.
Table 2. Characteristics of the pile and surrounding soil.
Name of the ParameterSymbolValueUnits
Characteristics of the pile
Modulus of deformation of the pile E p 3 × 107kPa
Radius of the pile a 1 0.30m
Pile spacing * 2 b 2 2.40m
Length of the pile l 1 30m
Characteristics of the soil
Modulus of deformation of the surrounding soil E s 2 × 104kPa
Modulus of deformation of the underlying soil mass E 2 3 × 104kPa
Poisson’s ratio of the underlying soil mass ν 2 0.36-
Coefficient of viscosity of the surrounding soil η s 1 × 1012Poise
*—in the present study, the pile spacing is equal to the diameter of the computational cell (Figure 1).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ter-Martirosyan, A.Z.; Ermoshina, L.Y.; Anzhelo, G.O. Settlement of a Pile Foundation Considering Linear and Rheological Properties of Soils. Buildings 2024, 14, 2830. https://doi.org/10.3390/buildings14092830

AMA Style

Ter-Martirosyan AZ, Ermoshina LY, Anzhelo GO. Settlement of a Pile Foundation Considering Linear and Rheological Properties of Soils. Buildings. 2024; 14(9):2830. https://doi.org/10.3390/buildings14092830

Chicago/Turabian Style

Ter-Martirosyan, Armen Z., Lyubov Yu. Ermoshina, and George O. Anzhelo. 2024. "Settlement of a Pile Foundation Considering Linear and Rheological Properties of Soils" Buildings 14, no. 9: 2830. https://doi.org/10.3390/buildings14092830

APA Style

Ter-Martirosyan, A. Z., Ermoshina, L. Y., & Anzhelo, G. O. (2024). Settlement of a Pile Foundation Considering Linear and Rheological Properties of Soils. Buildings, 14(9), 2830. https://doi.org/10.3390/buildings14092830

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