Numerical Modeling in Geotechnical Engineering

A special issue of Geosciences (ISSN 2076-3263). This special issue belongs to the section "Geomechanics".

Deadline for manuscript submissions: closed (28 February 2021) | Viewed by 32541

Special Issue Editors


E-Mail Website
Guest Editor
Faculty of Engineering and Architecture, University of Enna “Kore”, 94100 Enna, Italy
Interests: geotechnical earthquake engineering; dynamic characterization of soil; liquefaction; site response analysis; foundation; landslide

E-Mail Website
Guest Editor
Department of Civil Engineering and Architecture, University of Catania, 95124 Catania, Italy
Interests: soil structure interaction; numerical modeling; earthquake engineering; tunnels; geotechnical innovative materials

E-Mail Website
Guest Editor
Department of Engineering, University of Palermo, 90133 Palermo, Italy
Interests: soil behaviour; foundation; tunnel

E-Mail Website
Guest Editor
ARISTOTLE University of Thessaloniki, 541 24 Thessaloniki, Greece
Interests: geotechnical earthquake engineering; soil-foundation-structure interaction; experimental soil-foundation-structure interaction; performance based design; structural dynamics; numerical analysis; foundation design and analysis; seismic behavior and rehabilitation of historical buildings and monuments; soil dynamics; soil mechanics

Special Issue Information

In recent decades, numerical simulation has increasingly been used for the analysis of stress–strain levels in many fields of geotechnical engineering toward better accounting for the possible initial and boundary conditions associated with micro- to megascale problems, including groundwater flow, fully coupled soil–structure systems, and environmentally friendly innovative materials.

In addition to the developments related to the methods itself (new constitutive models for soil and rock, new numerical procedures, and calculation methods), the role of the numerical approach has evolved from the research field into a daily engineering tool due to the increase of computer power.

This Special Issue will collect contributions on recent research advances and/or well-documented applications of numerical modeling in static and dynamic geotechnical engineering.

Contributions regarding tunnels and deep excavations in urban areas, soil–structure interaction issues, cultural heritage protection, retaining walls, foundations, dams, slope stability, and innovative materials are welcome. Applications of numerical modeling to dynamics, geotechnical earthquake engineering, and constitutive models including of unsaturated soil and soil–structure interface are also welcome.

Finally, we encourage the submission of contributions concerning operative applications with the experimental validation of the models.

Dr. Valentina Lentini
Dr. Maria Rossella Massimino
Dr. Maurizio Ziccarelli
Dr. Dimitris Pitilakis
Guest Editors

Manuscript Submission Information

Manuscripts should be submitted online at www.mdpi.com by registering and logging in to this website. Once you are registered, click here to go to the submission form. Manuscripts can be submitted until the deadline. All submissions that pass pre-check are peer-reviewed. Accepted papers will be published continuously in the journal (as soon as accepted) and will be listed together on the special issue website. Research articles, review articles as well as short communications are invited. For planned papers, a title and short abstract (about 100 words) can be sent to the Editorial Office for announcement on this website.

Submitted manuscripts should not have been published previously, nor be under consideration for publication elsewhere (except conference proceedings papers). All manuscripts are thoroughly refereed through a single-blind peer-review process. A guide for authors and other relevant information for submission of manuscripts is available on the Instructions for Authors page. Geosciences is an international peer-reviewed open access monthly journal published by MDPI.

Please visit the Instructions for Authors page before submitting a manuscript. The Article Processing Charge (APC) for publication in this open access journal is 1800 CHF (Swiss Francs). Submitted papers should be well formatted and use good English. Authors may use MDPI's English editing service prior to publication or during author revisions.

Keywords

  • constitutive models
  • numerical modeling
  • soil–structure interaction
  • earthquake engineering
  • liquefaction
  • slope stability
  • foundation
  • tunnels
  • environmentally friendly innovative materials
  • cultural heritage protection

Published Papers (11 papers)

Order results
Result details
Select all
Export citation of selected articles as:

Research

15 pages, 45034 KiB  
Article
Influence of a Thin Horizontal Weak Layer on the Mechanical Behaviour of Shallow Foundations Resting on Sand
by Maurizio Ziccarelli and Marco Rosone
Geosciences 2021, 11(9), 392; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11090392 - 16 Sep 2021
Cited by 1 | Viewed by 1603
Abstract
The presence of minor details of the ground, including soil or rock masses, occurs more frequently than what is normally believed. Thin weak layers, shear bands, and slickensided surfaces can substantially affect the behaviour of foundations, as well as that of other geostructures. [...] Read more.
The presence of minor details of the ground, including soil or rock masses, occurs more frequently than what is normally believed. Thin weak layers, shear bands, and slickensided surfaces can substantially affect the behaviour of foundations, as well as that of other geostructures. In fact, they can affect the failure mechanisms, the ultimate bearing capacity of footings, and the safety factor of the geotechnical system. In this research, numerically conducted through Finite Element Code Plaxis 2D, the influence of a horizontal thin weak layer on the mechanical behaviour of shallow footings was evaluated. The obtained results prove that the weak layer strongly influences both the failure mechanism and the ultimate bearing capacity if its depth is lower than two to four times the footing width. In fact, under these circumstances, the failure mechanisms are always mixtilinear in shape because the shear strains largely develop on the weak layer. However, the reduction in the ultimate bearing capacity is a function of the difference between the shear strength of the foundation soil and the layer. The presence of a thin weak layer decreases the ultimate bearing capacity up to 90%. In conclusion, this research suggests that particular attention must be paid during detailed ground investigations to find thin weak layers. Based on the obtained results, it is convenient to increase the soil volume investigation to a depth equal to four times the width of the foundation. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

18 pages, 17993 KiB  
Article
System Identification of Mosques Resting on Soft Soil. The Case of the Suleiman Mosque in the Medieval City of Rhodes, Greece
by Anna Karatzetzou, Dimitris Pitilakis and Stella Karafagka
Geosciences 2021, 11(7), 275; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11070275 - 30 Jun 2021
Cited by 4 | Viewed by 2225
Abstract
The present study focuses on the dynamic system identification of the Suleiman Mosque minaret in the medieval city of Rhodes, Greece. Suleiman Mosque was built in 1522 at the site of the destroyed Christian Church of the Apostles. First, we performed sets of [...] Read more.
The present study focuses on the dynamic system identification of the Suleiman Mosque minaret in the medieval city of Rhodes, Greece. Suleiman Mosque was built in 1522 at the site of the destroyed Christian Church of the Apostles. First, we performed sets of ambient vibration measurements at the minaret of the monument. Based on these data, we calculated the eigenproperties of the minaret. Next, we modeled the monument in three dimensions, using the finite element method. Six numerical models were considered. Model Ι is the simplest one (isolated, fixed base minaret). Model VI is the most complicated one (simulation of the whole mosque also considering soil–structure interaction and foundation flexibility). The calculated predominant periods and mode shapes of Models I–VI are validated against the microtremor field measurements, recorded on the minaret’s two floors and ground level. We elaborate on the reliability of finite element models for earthquake response evaluation, considering soil–structure interaction and foundation flexibility on the mode shape eigenfrequencies. Additionally, we discuss the seismic response of the minaret compared to the whole monument. We observed no significant difference in the first two modes of response, implying that the minaret’s dynamic behavior is slightly affected by the entire mosque’s presence. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

21 pages, 4957 KiB  
Article
The Use of Polyurethane Injection as a Geotechnical Seismic Isolation Method in Large-Scale Applications: A Numerical Study
by Michele Placido Antonio Gatto, Valentina Lentini, Francesco Castelli, Lorella Montrasio and Davide Grassi
Geosciences 2021, 11(5), 201; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11050201 - 04 May 2021
Cited by 22 | Viewed by 2911
Abstract
This paper analyses the effect of polyurethane injections on the seismic surficial response of cohesionless soils. For this purpose, dynamic finite element numerical analyses were performed through GiD + OpenSees. Both the soil and the composite material, resulted after the expansion of the [...] Read more.
This paper analyses the effect of polyurethane injections on the seismic surficial response of cohesionless soils. For this purpose, dynamic finite element numerical analyses were performed through GiD + OpenSees. Both the soil and the composite material, resulted after the expansion of the injected polyurethane, are modelled with a nonlinear hysteretic constitutive model. Based on the polyurethane percentage, a homogenisation of the characteristics was considered for the composite material: linear for density and damping, and exponential (experimentally calibrated) for the stiffness. An expansion coefficient quantifies how much the injected polyurethane expands: three expansion coefficients were considered, each of them related to a different polyurethane density. For the evaluation of the foam stiffness, a linear stiffness–density correlation was used, derived after impact tests. Results showed that polyurethane reduces the surficial accelerations proportionally to the ratio of its seismic impedance and volumetric percentage with respect to the soil seismic impedance and total volume. This is a preliminary indication for the design of polyurethane injections in cohesionless soils for seismic acceleration reduction. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

31 pages, 11664 KiB  
Article
Effects of Soil-Foundation-Interaction on the Seismic Response of a Cooling Tower by 3D-FEM Analysis
by Francesco Castelli, Salvatore Grasso, Valentina Lentini and Maria Stella Vanessa Sammito
Geosciences 2021, 11(5), 200; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11050200 - 03 May 2021
Cited by 32 | Viewed by 3524
Abstract
This paper reports on the results of soil-foundation numerical modelling and the seismic response of a cooling tower founded on piles of a petrochemical facility located in the city of Augusta (Sicily, Italy). The city was affected in the past by some destructive [...] Read more.
This paper reports on the results of soil-foundation numerical modelling and the seismic response of a cooling tower founded on piles of a petrochemical facility located in the city of Augusta (Sicily, Italy). The city was affected in the past by some destructive earthquakes (1693, 1848, and 1990) that damaged a large territory of Southeastern Sicily, which was characterized by a very high seismic hazard. With this aim, the paper reports the FEM modelling of the seismic behaviour of the coupled soil-structure system. To determine the soil profile and the geotechnical characteristics, laboratory and in situ investigations were carried out in the studied area. The seismic event occurred in January 1693 and has been chosen as a scenario earthquake. Moreover, a parametric study with different input motions has also been carried out. A Mohr-Coulomb model has been adopted for the soil, and structural elements have been simulated by means of an elastic constitutive model. Two different vertical alignments have been analysed, considering both the free-field condition and the soil-structure interaction. The dynamic response has been investigated in terms of accelerations, response spectra, and amplification functions. The results have also been compared with those provided by Italian technical regulations. Finally, the seismic response of the coupled soil-structure system has been further examined in terms of peak bending moments along the pile foundation, emphasizing the possibility of a kinematic interaction on piles induced by the seismic action. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

16 pages, 3234 KiB  
Article
Stability of Embankments Resting on Foundation Soils with a Weak Layer
by Maurizio Ziccarelli and Marco Rosone
Geosciences 2021, 11(2), 86; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11020086 - 13 Feb 2021
Cited by 3 | Viewed by 2341
Abstract
The presence of weak layers in geotechnical systems, including soil or rock masses, both natural and man-made, is more frequent than is normally believed. Weak layers can affect both failure mechanisms, in drained and in undrained conditions, as well as in static and [...] Read more.
The presence of weak layers in geotechnical systems, including soil or rock masses, both natural and man-made, is more frequent than is normally believed. Weak layers can affect both failure mechanisms, in drained and in undrained conditions, as well as in static and seismic conditions, and the safety factor. In the present study, conducted numerically using the finite-element method (FEM) Plaxis 2D code, the influence of a horizontal thin weak layer on stress and strain distribution, on failure mechanisms and on the overall stability of an embankment was evaluated. The results obtained prove that when the weak layer is located at a significant depth from the foundation plane, the failure mechanisms are normally mixtilinear in shape because the shear strains largely develop on the weak layer. As a result, the safety factor highly decreases compared to the same case without a weak layer. Then, in the presence of weak layers, even embankments that, if founded on homogeneous soils, would have very high global safety factors (higher than 2) can become unstable, i.e., the safety factor can become unitary. So particular attention must be paid during detail ground investigations to finding thin weak layers. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

13 pages, 4775 KiB  
Article
A Geospatial Approach for Mapping the Earthquake-Induced Liquefaction Risk at the European Scale
by Francesca Bozzoni, Roberta Bonì, Daniele Conca, Claudia Meisina, Carlo G. Lai and Elisa Zuccolo
Geosciences 2021, 11(1), 32; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11010032 - 08 Jan 2021
Cited by 10 | Viewed by 3392
Abstract
This paper presents a geospatial methodology for zoning the earthquake-induced soil liquefaction risk at a continental scale and set-up in a Geographic Information System (GIS) environment by coupling data-driven and knowledge-driven approaches. It is worth mentioning that liquefaction is a phenomenon of soil [...] Read more.
This paper presents a geospatial methodology for zoning the earthquake-induced soil liquefaction risk at a continental scale and set-up in a Geographic Information System (GIS) environment by coupling data-driven and knowledge-driven approaches. It is worth mentioning that liquefaction is a phenomenon of soil instability occurring at a very local spatial scale; thus, the mega-zonation of liquefaction risk at a continental scale is a hard facing challenge. Since the risk from natural disasters is the convolution of hazard, vulnerability, and exposure, the liquefaction risk mapping is based on the combination of geospatial explanatory variables, available at the continental scale, of the previously listed three assumed independent random variables. First, by applying a prediction model calibrated for Europe, the probability of liquefaction is mapped for the whole continent. Then, the Analytical Hierarchy Process (AHP) is adopted to identify areas that have a high risk of liquefaction, taking into account proxy data for exposure. The maps are computed for different levels of severity of ground shaking specified by three return periods (i.e., 475, 975, and 2475 years). A broad variety of stakeholders would benefit from the outcomes of this study, such as civil protection organizations, insurance and re-insurance companies, and infrastructure operators. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

25 pages, 14524 KiB  
Article
New Stress Reduction Factor for Evaluating Soil Liquefaction in the Coastal Area of Catania (Italy)
by Salvatore Grasso, Maria Rossella Massimino and Maria Stella Vanessa Sammito
Geosciences 2021, 11(1), 12; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences11010012 - 28 Dec 2020
Cited by 19 | Viewed by 3733
Abstract
In this paper, a study concerning the soil liquefaction potential in the city of Catania is presented. The stress-based liquefaction analysis framework for cohesionless soil includes a function that describes fundamental aspects of dynamic site response, i.e., the shear stress reduction coefficient, r [...] Read more.
In this paper, a study concerning the soil liquefaction potential in the city of Catania is presented. The stress-based liquefaction analysis framework for cohesionless soil includes a function that describes fundamental aspects of dynamic site response, i.e., the shear stress reduction coefficient, rd, which depends on several factors (depth; earthquake and ground motion characteristics; dynamic soil properties). Various relationships of rd are reported in literature because of the importance of assessment of CSR. Herein, new variations of rd with depth have been obtained using different deterministic earthquake scenarios as input motion. The relationships are based on large numbers of site response analyses for different site conditions. The new relationships obtained have been used for the evaluation of the liquefaction potential in the area of the Catania Harbour. The liquefaction resistance has been evaluated by the horizontal stress index (KD) from seismic dilatometer Marchetti tests (SDMTs). Various correlations were developed to estimate the CRR from KD, expressed in form of CRR-KD curves to differentiate between liquefiable and non-liquefiable zones. In this study three different CRR-KD curves have been used. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

11 pages, 3474 KiB  
Article
Numerical Evaluation of Natural Periods and Mode Shapes of Earth Dams for Probabilistic Seismic Hazard Analysis Applications
by Paolo Zimmaro and Ernesto Ausilio
Geosciences 2020, 10(12), 499; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences10120499 - 12 Dec 2020
Cited by 4 | Viewed by 2407
Abstract
The evaluation of natural periods and related mode shapes of earth dams represents a critical issue when performing structure-specific probabilistic seismic hazard analyses (PSHA). The identification of critical scenario events, using techniques such as disaggregation of the seismic hazard, and the calculation of [...] Read more.
The evaluation of natural periods and related mode shapes of earth dams represents a critical issue when performing structure-specific probabilistic seismic hazard analyses (PSHA). The identification of critical scenario events, using techniques such as disaggregation of the seismic hazard, and the calculation of a suitable target spectrum for ground motion selection and scaling procedures (e.g., the conditional mean spectrum), require at least the knowledge of the fundamental period of the system. This problem can be solved using analytical, numerical, and/or empirical techniques. We present several linear elastic modal analyses for an earth dam located in Southern Italy, using a numerical solution of the generalized eigenvalue problem obtained by the finite element method (FEM). Our numerical experiments are performed, testing various assumptions on boundary conditions, degree of saturation, and the distribution of geotechnical characteristics of the dam’s materials. We then compare our results against existing analytical solutions. We show that ignoring soil–structure interaction effects due to the flexibility of the dam foundation (i.e., under the assumption of fixed base) can lead to a substantial underestimation of the fundamental period of the dam. This effect should be carefully addressed when modal analysis results are used in PSHA-related applications. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

18 pages, 10411 KiB  
Article
Numerical Modelling of Structures Adjacent to Retaining Walls Subjected to Earthquake Loading
by Xiaoyu Guan and Gopal S. P. Madabhushi
Geosciences 2020, 10(12), 486; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences10120486 - 03 Dec 2020
Cited by 3 | Viewed by 2888
Abstract
In an urban environment, it is often necessary to locate structures close to existing retaining walls due to congestion in space. When such structures are in seismically active zones, the dynamic loading attracted by the retaining wall can increase. In a novel approach [...] Read more.
In an urban environment, it is often necessary to locate structures close to existing retaining walls due to congestion in space. When such structures are in seismically active zones, the dynamic loading attracted by the retaining wall can increase. In a novel approach taken in this paper, finite element-based numerical analyses are presented for the case of a flexible, cantilever sheet pile wall with and without a structure on the backfill side. This enables a direct comparison of the influence exerted by the structure on the dynamic behaviour of the retaining wall. In this paper, the initial static bending moments and horizontal stresses prior to application of any earthquake loading are compared to Coulomb’s theory. The dynamic behaviour of the retaining wall is compared in terms of wall-top accelerations and bending moments for different earthquake loadings. The dynamic structural rotation induced by the differential settlements of the foundations is presented. The accelerations generated in the soil body are considered in three zones, i.e., the free field, the active and the passive zones. The differences caused by the presence of the structure are highlighted. Finally, the distribution of horizontal soil pressures generated by the earthquake loading behind the wall, and in front of the wall is compared to the traditional Mononobe-Okabe type analytical solutions. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

22 pages, 18404 KiB  
Article
Innovative Seismic Microzonation Maps of Urban Areas for the Management of Building Heritage: A Catania Case Study
by Glenda Abate, Simone Bramante and Maria Rossella Massimino
Geosciences 2020, 10(12), 480; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences10120480 - 26 Nov 2020
Cited by 6 | Viewed by 2177
Abstract
Several urban areas in the Mediterranean have already been subjected to seismic microzonation studies aimed at determining the acceleration expected on the ground surface, therefore mitigating the associated seismic risks. These studies have been generally related to free-field conditions. The present paper shows [...] Read more.
Several urban areas in the Mediterranean have already been subjected to seismic microzonation studies aimed at determining the acceleration expected on the ground surface, therefore mitigating the associated seismic risks. These studies have been generally related to free-field conditions. The present paper shows innovative seismic microzonation maps based on a large-scale estimate of soil-structure interaction (SSI) effects on design accelerations for some areas characterized by a high seismic risk in Catania, Italy. The proposed procedure combined: (1) geotechnical characteristics; (2) building features; and (3) 1-D seismic response analyses in free-field conditions. The seismic hazard and site effects were evaluated using artificial inputs and inputs recorded recently in Catania. Structural fundamental periods and related spectral accelerations, considering both the fixed-base building configuration and flexible-base configuration, were mapped in the Google My Maps environment. These results showed that SSI often had a beneficial effect, but sometimes it had detrimental effects, especially for some masonry buildings. These maps provided important information for planning the seismic retrofitting of investigated buildings, which were based on more detailed analyses of SSI and the developed maps requiring them. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Figure 1

25 pages, 10790 KiB  
Article
Geotechnical Analysis and 3D Fem Modeling of Ville San Pietro (Italy)
by Rossella Bovolenta and Diana Bianchi
Geosciences 2020, 10(11), 473; https://0-doi-org.brum.beds.ac.uk/10.3390/geosciences10110473 - 22 Nov 2020
Cited by 5 | Viewed by 2593
Abstract
The paper describes the three-dimensional numerical model of Ville San Pietro, an Italian village subject to slope movements causing damage. The church (dating back to 1776), which is the most significant building of the area, is modelled too. The information from geotechnical and [...] Read more.
The paper describes the three-dimensional numerical model of Ville San Pietro, an Italian village subject to slope movements causing damage. The church (dating back to 1776), which is the most significant building of the area, is modelled too. The information from geotechnical and geophysical surveys on field are used to define the model geometry and the soil properties. A finite element code is adopted to simulate the slope behavior in occurrence of water table fluctuations, detected by piezometers, and to evaluate the slope displacements and stability. The validation of the model is carried out using the inclinometer and interferometry measures and by on-site inspections. The model demonstrated a good ability to simulate the slope behavior during the raising and lowering of the water table. The critical areas computed by the numerical code are in good accordance to the actual portions affected by soil displacements and damages. The modelling presented in this paper is crucial for future analyses that will take advantage of an innovative monitoring system, which will be installed on site. Full article
(This article belongs to the Special Issue Numerical Modeling in Geotechnical Engineering)
Show Figures

Graphical abstract

Back to TopTop