Numerical modeling of pit excavation with account for regional peculiarities of mechanical behavior of foundation soils
- Authors: Mangushev R.A.1, Bashmakov I.B.1, Paskacheva D.A.1
- Affiliations:
- Saint Petersburg University of Architecture and Civil Engineering
- Issue: Vol 15, No 3 (2024)
- Pages: 56–67
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/CG/article/view/4365
- DOI: https://doi.org/10.15593/2224-9826/2024.3.05
- Cite item
Abstract
The problems of undrained behavior of weak clayey soils of St. Petersburg are considered. The importance of taking into account the mechanism of shear hardening in calculations of excavations in weak clayey soils has been demonstrated. The problem of using undrained shear resistance as a strength parameter in full stress evaluation is shown. The modification of Y.I. Soloviev's theory of instantaneous strength has been presented with account for the formation of excess pore pressures under deviatoric load-ing in the plane formulation. With its help the calculated value of resistance to undrained shear through effective strength parameters for use in numerical modeling complexes is obtained. Numerical realization of the nonlinear model of soil with shear hardening in the software package implementing the finite element method has been performed. This model utilizes the proposed methodology for estimating undrained strength. Comparison of the calculation results of the proposed model and the Mohr-Coulomb model has been performed.
Full Text
IntroductionExpansion of the field of study of undrained behavior of soil foundation is associated with the practice of excavation pits calculation in the conditions of weak clayer soils. The study of properties of weak clayey soils is connected with their low mechanical characteristics with wide distribution in the active zone of the foundation during construction works.This is especially relevant for St. Petersburg, where the thickness of weak soils reaches an average of 15 m. With respect to the active development of underground construction in the city it is necessary to provide special approaches to assess the bearing capacity of weak clayey soils in order to make safe and cost-efficient engineering solutions during works.One of these approaches is the consideration of the weak clayey soils properties over time since they are characterized by low filtration coefficients. In this case, the possibility of formation and dispersion of excessive pore pressures is taken into ac-count, what is important for impacts on the foundation in a limited period of time (for example, the construction of a pit).Research in this area begins with the works of K. Terzaghi and is associated with the de-velopment of the theory of filtration consolidation. The influence of excessive pore pressures on soil strength was studied by A. Skempton [1–3], who introduced pore pressure coefficients and also conducted research for the examination of the clayer soils behavior. Significant research in the field of undrained behavior of soils was also carried out by A. Bishop, D. Henkel [4, 5], C. Ladd [6, 7] and others. Important studies in the field of undrained behavior were performed by the native scientists: A.G. Shashkin [8], V.M. Ulitsky [9], Y.I. Soloviev [10] and others.In the world practice, with reference to the possibility of forming excessive pore pressures of the foundation soils two approaches have been formed for the analysis:– analysis in effective stresses using effective strength parameters;– analysis in total stresses using full strength parameters.Evaluation of the bearing capacity of the foundation in total stresses is considered to be a more conservative approach. In the calculation it is used a well-known parameter – undrained shear resistance. Describing the shear resistance by only one parameter, it is impossible to take into account a number of dependencies strongly influencing its value:– resistance to undrained shear depends on the degree of disruption of the sample structure which is confirmed by A. Skempton [3], P. Mayne [11], M. Jamiolkowski [7];– type of stress-strain state of the sample during testing determines the value of undrained shear resistance which is confirmed by P. Mayne [11].Also, the analysis in total stresses does not provide the ability to perform time calculations taking into account the filtration processes in the foundation [12].During calculation in effective stresses, the angle of internal friction and cohesion are used as strength parameters. The use of effective strength parameters allows depending less on the quality of soil test specimens, what was proved by A. Skempton [3]. It also becomes possible to perform calculations in time. However, here an important role begins to play the laws of defor-mation embedded in the soil model.The purpose of this study is to develop a method for assessing the undrained strength of weak clayey soils during calculation at effective stresses and its implementation in the numeri-cal simulation software package.To achieve this goal, the following tasks have been completed:– modification of the theory of instantaneous strength by Y.I. Solovyov considering the pa-rameter which takes into account the development of excess pore pressures during shear defor-mation;– obtaining the calculated value of resistance to undrained shear Cu based on effective strength parameters and the undrained strength parameter Af;– development and numerical implementation of a nonlinear model realizing shear harden-ing with account for the obtained value of resistance to undrained shear.Experimental Methods and Theoretical ApproachesMathematical modeling of undrained soil behaviorThe special properties of clays play an important role in the study of undrained behavior of clayey soils. In the context of the pit excavation problem it is important to consider how these properties manifest themselves under shear deformation.This can be shown graphically by means of effective soil stress trajectories, which describe the change in the stress state of the specimen under loading. The theoretical description of the dependence of the stress trajectory on the soil model in numerical modeling and the role of the A. Skempton parameter are presented in the previous article by the authors [13].Thus, the main difficulty in using the analysis in effective stresses is the dependence of the ultimate strength realized in the calculation on the theoretical laws embedded in the soil model (see Fig. 1). Fig. 1. Trajectories of effective stresses [13]Рис. 1. Траектории эффективных напряжений [13]Method for Estimating Undrained Strength by Effective Strength Parameters under Conditions of Plane DeformationThe parameter Af can be used to obtain the value of undrained shear resistance Cu, which will take into account the peculiarities of deformation of weak foundation soils. However, de-termined under triaxial axisymmetric compression Af it will not correspond to the plane stress state, to which the problem of the excavation can be referred.For direct determination of excess pore pressure under shear deformation in plane defor-mation conditions it is possible to use the simple shear device (DSS), but in domestic practice it is found only in a limited number of scientific laboratories and testing on it is not encountered in the course of commercial engineering-geological surveys. This is facilitated by the lack of regulatory-technical base for testing in this device.In conditions of simple shear, the Skempton parameter can be determined by the following formula: (1)where – Skempton parameter for flat SSS; – pore pressure; – tangential stress change; – lowest main effective pressure; – the highest principal effective stress.To obtain the resistance to undrained shear in the case of aflat setting we modify the law of A. Skempton [1]: (2)where – Skempton Parameter for Flat SSS.To calculate the value of resistance to undrained shear, we will use the theory of instanta-neous strength of Yu.K. Solovyov [10]: (3)where – average stress in soil mass; – shear stresses in a soil mass; – effective specific cohesion; – effective angle of internal friction.Numerical Implementation of a Nonlinear Model with Shear HardeningTo conduct numerical experiments according to the proposed methodology we realize a nonlinear model with shear hardening. To describe the stress-strain relationship under devia-toric loading, we adopt the expression proposed by R.L. Kondner [14]. A similar expression was previously proposed by A.I. Botkin [15].Instead of the largest principal axial deformation as a variable, let us consider the parame-ter: (4)Then the expression for the relationship between stresses and strains under deviatoric load-ing will be as follows: (5)where – the inverse of the tangent modulus; – the inverse of the asymptotic stress deviator. Plastic component of deformation be found as: (6)where – elastic component of deformation.From the generalized law of R. Hooke we define : (7)where – Poisson's ratio; – modulus of elasticity.Then substituting expressions (5) and (7) into (6) we obtain: (8)Let us define the value as a value inverse to the modulus of elasticity and independent of stresses: (9)The value b is calculated as: (10)where – asymptotic value of the stress deviator; – stress deviator limit value; – frac-ture ratio.Let us substitute the ratio of the stress deviator to the stress deviator at fracture: (11)where – mobilized resistance to undrained shear.Let's substitute expressions (9), (10) and (11) into (8) and determine the value of : (12)As yield strength, we adopt Tresk's criterion to describe strength through undrained shear resistance: (13)Then is the hardening parameter. The law of plastic flow is assumed to be non-associated. As the plastic potential we take the Coulomb-Mohr law with the dilatancy angle : (14)The above relations are sufficient to realize the nonlinear solution. As an algorithm for nonlinear solution we adopt the method of initial stresses with step-by-step loading. The stiff-ness matrix is assumed to be variable depending on the achieved stress level. A detailed de-scription of the algorithm realization is described in A.B. Fadeev [16] and V.N. Paramonova [17].Results and discussionThe analytical solution obtained in the previous articles of the authors for calculating the parameter [13] and analytical solution for finding the pressure on the enclosing structure [18] show the necessity of correct consideration of the shear component of pore pressure in un-drained calculations of weak clayey soils. Also, it is noted in the works of C.P. Wroth [19], C. Surarak [20] and others.To develop a method for estimating undrained shear resistance based on effective strength parameters and deformation features of foundation soils, we modify Soloviev's theory of instan-taneous strength (3) using (2) to allow for shear deformation features. Here is a variant of the formula for the case of instantaneous loading of the soil mass, where : (15)where – natural average effective stresses in the soil; – coefficient of lateral soil pres-sure.By substituting the value, the solution of Y.I. Soloviev (3) can be obtained. This solution allows us to estimate the resistance to undrained shear with account for shear pore pressure, which develops significantly under deviatoric loading of weak clayey soils. To evaluate the in-fluence of shear hardening and the development of shear pore pressure on the deformations of the soil mass, we perform a numerical calculation of the excavation in three formulations:– undrained calculation in full stresses in the developed soil model with determination of strength parameters according to the proposed methodology (variant 1) – developed by PC;– undrained calculation in effective stresses in the Mohr-Coulomb model (variant 2) – Plaxis 2D;– drained calculation in effective stresses in Mohr-Coulomb model (variant 3) – developed by PC.The realization is carried out in a program complex based on the finite element method, developed by the scientific team which includes the authors [21–23].The calculation is performed in two stages: formation of the natural stress-strain state and pit digging. The general view of the calculation scheme is presented in Fig. 2.As soils we will take two engineering-geological elements characteristic for the thickness of weak soils of St. Petersburg districts. Characteristics are presented in the Table.As the resultsofthe calculation will present the isofields of displacements of the soil mass and the oints of plastic deformation in the developed software pachage in Fig. 3–5. Fig. 2. General view of the calculation scheme in the last phase of the calculationРис. 2. Общий вид расчетной схемы на последней фазе расчетаPhysical and mechanical characteristics of soils Физико-механические характеристики грунтов Characteris-tics The proposed model and methodology Mohr-Coulomb: Drained и Undrained ASoil № 1 (top layer)Physical , Mechanical , , Soil № 2 (bottom layer)Physical , Mechanical , , , 1 2 3Fig. 3. Isopoles of vertical deformations: variant 1 Uy = 51/15 mm, variant 2 Uy = 45/12 mm, variant 3 Uy = 83/2 mmРис. 3. Изополя вертикальных деформаций: вариант 1 Uy = 51/15 мм, вариант 2 Uy = 45/12 мм, вариант 3 Uy = 83/2 мм 1 2 3Fig. 4. Isopoles of horizontal deformations: variant 1 Ux = 19/0 mm, variant 2 Ux = 16/0 mm, variant 3 Ux = 11/4 mmРис. 4. Изополя горизонтальных деформаций: вариант 1 Ux = 19/0 мм, вариант 2 Ux = 16/0 мм, вариант 3 Ux = 11/4 мм 1 2 3Fig. 5. Plastic points: variant 1, variant 2, variant 3Рис. 5. Пластические точки: вариант 1, вариант 2, вариант 3The developed model shows larger displacements than the undrained calculation and smaller displacements than the drained one in the Mohr-Coulomb model. The drained calcula-tion shows larger zones of plastic deformation development. Comparison of undrained calcula-tions in the proposed model and in the Mohr-Coulomb model clearly illustrates the influence of plastic shear deformations on the calculation results.ConclusionThe influence of the type of SSS and deformation features of weak clayey soils on shear pore pressure requires either an adequate estimation of the magnitude of the undrained shear resistance or model selection with respect to the regional characteristics of the soils.It is convenient to estimate the magnitude of shear pore pressure using the Skempton pa-rameter Af. Determination of the parameter should be carried out in tests with pore pressure fix-ation with a calculation scheme realizing a SSS similar to reality.Comparison of soil mass deformations in the Mohr-Coulomb model and the proposed model shows the importance of the shear component of deformation for this type of problems.About the authors
R. A. Mangushev
Saint Petersburg University of Architecture and Civil Engineering
I. B. Bashmakov
Saint Petersburg University of Architecture and Civil Engineering
D. A. Paskacheva
Saint Petersburg University of Architecture and Civil Engineering
References
- Skempton A. The pore-pressure coefficients A and B. Geotechnique, 1954, vol. 4, pp. 143-147. doi: 10.1680/geot.1954.4.4.143
- Skempton A.W. The ø = 0 analysis of stability and its theoretical basis. 2nd International Conference on Soil Mechanics and Foundation Engineering, Proceedings, 1948, vol. 1.
- Skempton A.W., Sowa V.A. The behaviour of saturated clays during sampling and testing. Geotechnique, 1963, vol. 13, no. 4, pp. 269-290. doi: 10.1680/geot.1963.13.4.269
- Bishop A.W., Henkel D.J. The measurement of soil properties in the triaxial test. London, 1962, 228 p.
- Bishop A.W., Skinner A.E. The influence of high pore-water pressure on the strength of cohesionless soils. Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences, 1977, vol. 284, no. 1318, pp. 91-130. DOI: doi.org/10.1098/rsta.1977.0001
- Sheahan T.C., Ladd C.C., Germaine J.T. Rate-dependent undrained shear behavior of saturated clay. Journal of Geotechnical Engineering, 1996, vol. 122, no. 2, pp. 99-108. doi: 10.1061/(ASCE)0733-9410(1996)122:2(99)
- Jamiolkowski M. [et al.]. New developments in field and laboratory testing of soils. Proceedings of the eleventh international conference on soil mechanics and foundation engineering, San Francisco, 12-16 August 1985. Balkema (AA), 1985.
- Shashkin A.G. Visco-elastic-plastic model of the behavior of clay soil. Urban development and geotechnical construction, 2011, vol. 2, pp. 1-32.
- Ulitsky V.M., Shashkin A.G., Shashkin K.G., Shashkin V.A. Basics of joint calculations of buildings and foundations. Publishing house of the Institute "Georeconstruction", 2014.
- Solovyov Yu.I., Karaulov A.M., Vaganov P.S. The theory of instantaneous strength and its application in calculating the stability of consolidating soil masses. Design and study of foundations of hydraulic structures: Proceedings of the conf. and meeting in hydraulic engineering, 1980, pp. 104-105.
- Mayne P.W. [et al.]. Geomaterial behavior and testing. Proceedings of the 17th International Conference on Soil Mechanics and Geotechnical Engineering, 2009, vol. 4, pp. 2777-2872. doi: 10.3233/978-1-60750-031-5-2777
- Çinicioğlu S.F., Toğrol E. Embankment design on soft clays. Journal of geotechnical engineering, 1991, vol. 117, no. 11, pp. 1691-1705. doi: 10.1061/(ASCE)0733-9410(1991)117:11(1691)
- Mangushev R.A., Bashmakov I.B., Paskacheva D.A., Kvashuk A.V. Mathematical modeling of undrained behavior of soils. International Journal for Computational Civil and Structural Engineering, 2023, vol. 19, no. 1, pp. 97-111. doi: 10.22337/2587-9618-2023-19-1-97-111
- Kondner R.L. Hyperbolic stress-strain response: cohesive soils. Journal of the Soil Mechanics and Foundations Division, 1963, vol. 89, no. 1, pp. 115-143. doi: 10.1061/JSFEAQ.0000479
- Botkin A.I. On the strength of bulk and brittle materials. News of the Scientific Research Institute of Hydraulic Engineering, 1940, vol. 26, pp. 205-236.
- Fadeev A.B. Finite element method in geomechanics. Moscow, Nedra, 1987, 221 p.
- Paramonov V.N. Finite element method for solving nonlinear geotechnical problems. Saint Petersburg, Georeconstruction, 2012.
- Mangushev R.A., Diakonov I.P., Bashmakov I.B., Paskacheva D.A. Calculation method of determining the earth pressure on the diaphragm wall considering the undrained soil behavior. In Smart Geotechnics for Smart Societies, CRC Press, 2023, pp. 1015-1021. doi: 10.1201/9781003299127-142
- Wroth C.P. The interpretation of in situ soil tests. Geotechnique, 1984, vol. 34, no. 4, pp. 449-489. doi: 10.1680/geot.1984.34.4.449
- Surarak C., Likitlersuang S., Wanatowski D., Balasubramaniam A., Oh E., Guan H. Stiffness and strength parameters for hardening soil model of soft and stiff Bangkok clays. Soils and foundations, 2012, vol. 52, no. 4, pp. 682-697.
- Kolyukayev I.S., Solving the problem of slope stability using a numerical method. Proceedings of the National (All-Russian) Scientific and Technical Conference, St. Petersburg, April 10–13, 2023 Prospects for modern construction, 2023, Saint Petersburg, Saint Petersburg State University of Architecture and Civil Engineering, pp. 65-79.
- Polunin V.M., Kolyukayev I.S., Gorkina M.R. Analytical and numerical methods for determining the stress-strain state of a soil massif for solving a planar problem. Smart Geotechnics for Smart Societies, CRC Press, 2023, pp. 1991-1998. doi: 10.1201/9781003299127-304
- Polunin V.M. Analytical and numerical methods for determining the stressed state of a soil mass when solving a plane problem. Housing construction, 2023, no. 9, pp. 27-40. doi: 10.31659/0044-4472-2023-9-27-40
Statistics
Views
Abstract - 7
PDF (English) - 4
Refbacks
- There are currently no refbacks.