Geotechnics, Volume 3, Issue 3 (September 2023) – 22 articles

Rutting is one of the most common types of distress in flexible pavement structures. There are two fundamental methods of designing pavement structures: conventional empirical methods and analytical approaches. Many analytical and empirical design procedures assume that rutting is mostly of asphalt origin and can be reduced by limiting the vertical deformation or stress at the top of the subgrade, but they do not quantify the rutting depth itself. Mechanistic–empirical models to predict the permanent deformations of unbound pavement layers have been well investigated and are rather common in North America; however, they are not widely utilized in the rest of the world. To date, there is no generally accepted, widely recognized, and documented procedure for calculating permanent deformations and thus for determining the rutting depth in flexible pavement courses originating from the unbound granular layers. This paper presents a layered calculation method with which the deformation of soil layers (base, subbase, and subgrade courses) under flexible pavements due to repeated traffic load can be determined. In the first step, the cyclic strain amplitude is calculated using a nonlinear material model that is based on particle size distribution parameters (d₅₀ and C_U) and dependent on the mean normal stress, relative density, and actual strain level. In the second step, the HCA (High Cycle Accumulation) model is used to calculate the residual settlement of each sublayer as a function of the number of cycles. It is shown that the developed model is suitable for describing different types of subgrades and pavement cross-sections. It is also demonstrated with finite element calculations that the developed model describes both the elastic and plastic strains sufficiently accurately. The developed model can predict the settlement and rutting of pavement structures with sufficient accuracy based on easily available particle size distribution parameters without the need for complex laboratory and finite element tests. Full article

The indirect tension test is an important laboratory test for rock characterization. The presence of rock fabric, such as schistosity, complicates the assessment of test results. One hundred and forty-five indirect tension tests were conducted on mica schist specimens to investigate the effect of schistosity orientation on failure mode and tensile strength. Tensile strength results did not provide a clear relationship between schistosity orientation and tensile strength, so the failure patterns were investigated. A new naming scheme for failure modes was developed, incorporating fracture patterns observed in the specimen faces and edges. The Single Mode failure group specimens had only one failure pattern that appeared on both specimen faces, either axial failure (seventy-three specimens), schistosity failure (six specimens), or out-of-plane failure (seven specimens). The Mixed Mode failure group had thirty-two specimens that exhibited one failure pattern on one face and another on the other. The Hybrid Mode failure group had twenty-seven specimens with multiple failure patterns on both specimen faces. It was noted that Mixed Mode and Hybrid Mode specimens with components of axial failure had higher indirect tensile strengths than specimens without elements of axial failures. Statistical analyses of the tensile strength data using Levene’s Test for equal variances and two-sample t-tests showed no statistical difference between the Mixed Mode and Hybrid Mode failure groups. However, there was a statistical difference between the tensile strengths of the Single Mode axial failure specimens and the combined Mixed Mode and Hybrid Mode failure groups. These results clearly emphasize that indirect tensile strength should be assessed using schistosity orientation and failure mode. Full article

This paper assesses the performance of an embankment constructed in 2010 with a stabilised expansive soil. Two types of treatment were employed at construction time: 4% lime and a mix of 2% lime and 3% cement. A sampling campaign was carried out in 2021 to evaluate the long-term performance of the stabilised soil properties. To assess the compressibility of the soil, oedometer tests were carried out on samples from different parts of the embankment. The results were compared to the compression curve of the untreated soil, also sampled in the same embankment. Complementary shrinkage tests were performed to investigate the effect of the treatment on swelling and shrinkage. The obtained results show that the yield stress of the material from the outer part was inferior to 100 kPa, similarly to the yield stress of the untreated soil, demonstrating a strong alteration in the effect of both treatments over time. This alteration was noticeable to a distance of approximately 2 m from the external surface. Beyond this distance, the performance of the soil was comparable to the behaviour of recently treated soil, with yield stresses close to 1000 kPa. These observations, similar for each treatment dosage, raise questions as to the durability of the treatment on the outer part of the backfill. Full article

Geotechnical engineering relies heavily on predicting soil strength to ensure safe and efficient construction projects. This paper presents a study on the accurate prediction of soil strength properties, focusing on hydrated-lime activated rice husk ash (HARHA) treated soil. To achieve precise predictions, the researchers employed two grey-box machine learning models—classification and regression trees (CART) and genetic programming (GP). These models introduce innovative equations and trees that readers can readily apply to new databases. The models were trained and tested using a comprehensive laboratory database consisting of seven input parameters and three output variables. The results indicate that both the proposed CART trees and GP equations exhibited excellent predictive capabilities across all three output variables—California bearing ratio (CBR), unconfined compressive strength (UCS), and resistance value (R_value) (according to the in-situ cone penetrometer test). The GP proposed equations, in particular, demonstrated a superior performance in predicting the UCS and R_value parameters, while remaining comparable to CART in predicting the CBR. This research highlights the potential of integrating grey-box machine learning models with geotechnical engineering, providing valuable insights to enhance decision-making processes and safety measures in future infrastructural development projects. Full article

The design of fiber-reinforced soil structures, such as embankments and pavements, can be carried out using the results of unconsolidated, undrained triaxial compression tests conducted on specimens at their “as-compacted” water content and analyzed in terms of total stresses. The effects of soil and fiber type on the mechanical behavior of fiber-reinforced soils have not been methodically or adequately examined in the past under these conditions, and the effects of fiber length and content on the shear strength parameters of fiber-reinforced soils need further experimental documentation. Accordingly, five soils ranging from “excellent” to “poor” materials for use in earthwork structures were tested in the present study, in combination with five types of polypropylene fibers having lengths ranging from 9 to 50 mm. Unconsolidated undrained triaxial compression tests were conducted on specimens at their “as-compacted” water content, with fiber contents ranging from 0.5 to 2% by weight of dry soil. Fiber reinforcement reduces the stiffness and increases the deformability of the soil. The fiber-reinforced soils exhibit a more ductile behavior in comparison with the unreinforced soils. A Mohr–Coulomb type linear failure criterion satisfactorily describes the shear strength behavior of fiber-reinforced soils in total stress terms. The cohesion values of the fiber-reinforced soils range between 61 kPa and 301 kPa and increase up to seven times in comparison with the cohesion values of the unreinforced soils. The variations of the angle of internal friction of soils due to fiber reinforcement are generally limited to ±25%. The cohesion improvement due to fiber reinforcement is increased with increasing fiber content and fiber length up to 30 mm and is inversely proportional to the fine-grained fraction and the cohesion of the unreinforced soil. Full article

The analysis of sands, and the foundation systems with which they interact, are largely dependent on macroscale behavioral parameters that represent the aggregated response of several microscale characteristics. This research paper examines the influence of surface texture, or smoothness, on the behavior of sands. The challenge of estimating or measuring smoothness, due to its microscale feature domain, is addressed through an examination of six artificially graded sand specimens. These specimens are evaluated both visually and numerically to characterize their surface smoothness. The first approach described is a simple visual method that uses a smoothness scale consistent with those of roundness and sphericity. This method, which can be performed with a tool as simple as a hand lens, evaluates a group of representative particles collectively. The second approach is also a visual evaluation, but it utilizes images obtained via scanning electronic microscopy, traditional optical microscopy, and newer low-cost digital microscopes that can be rapidly connected to a smartphone or laptop. To validate these visual estimates, a novel third approach is introduced. This approach is a more objective numerical analysis measurement technique that enables rapid and economic quantification of smoothness. This technique may assist both practitioners and academics in their understanding of the macroscale response of coarse-grained soils. In addition to the visual methods, this research also conducted several laboratory index tests to observe the mechanical behavior of the specimens, considering their particle shape and surface smoothness properties. The results indicate that angular sands have greater minimum and maximum void ratios, a larger difference between the minimum and maximum void ratios, greater critical state friction angles, and greater flow rates through an orifice of fixed size. When adjusted for surface smoothness using the proposed approach, the behavior of the sands—particularly the limit void ratio results—appears to be more predictable in some cases. These results provide additional evidence of particle smoothness contributing to the strength behavior of sand, which may be particularly useful in the domains of slope stability, land reclamation, soil–structure interaction, and soil dynamics. Full article

This article introduces a novel numerical scheme within the finite element method (FEM) to study soil heterogeneity, specifically focusing on the root–soil matrix in fracture treatments. Material properties, such as Young’s modulus of elasticity, cohesion, and the friction angle, are considered as randomly distributed variables. To address the inherent uncertainty associated with these distributions, a Monte Carlo simulation is employed. By incorporating the uncertainties related to material properties, particularly the root component that contributes to soil heterogeneity, this article provides a reliable estimation of the factor of safety, failure surface, and slope deformation, all of which demonstrate a progressive behavior. The probability distribution curve for the factor of safety (FOS) reveals that an increase in the root area ratio (RAR) results in a narrower range and greater certainty in the population mean, indicating reduced material variation. Moreover, as the slope angle increases, the sample mean falls within a wider range of the probability density curve, indicating an enhanced level of material heterogeneity. This heterogeneity amplifies the level of uncertainty when predicting the factor of safety, highlighting the crucial importance of accurate information regarding heterogeneity to enhancing prediction accuracy. Full article

This paper reviews numerical methods used to simulate desiccation cracks in clayey soils. It examines five numerical approaches: Finite Element (FEM), Lattice Boltzmann (LBM), Discrete Element (DEM), Cellular Automaton (CAM), and Phase Field (PFM) Methods. The paper presents a simplified description of the methods, including their basic numerical formulations. Several factors such as the multiphase nature of soils, heterogeneity, nonlinearities, coupling, scales of analysis, and computational aspects are discussed. The review highlights the characteristics, strengths, and limitations of each method. FEM shows a good capacity to deal with the thermo-hydromechanical behavior of clays when drying that complement well with the ability of DEM to deal with particle interactions as well as LBM, PFM, and CAM to deal with complex crack patterns. The article concludes by reviewing the integration of multiple numerical methods to enhance the simulation of desiccation cracks in clayey soils and proposing what is the best option to continue improving the study of this problem. Full article

Mudrocks are fine-grained clay-rich rocks that comprise different lithotypes forming more than 60% of all sedimentary rocks, and thus, they occur frequently in engineering projects either as natural ground or as made ground. These rocks may display a range of engineering behaviours controlled mostly by their composition and structural features. Due to rapid breakdown and susceptibility to volume changes, they may cause problems both during and after construction. Research into the susceptibility of mudrocks to breakdown aims to predict problematic behaviour and provide guidance for avoiding or mitigating these effects. Low-durability materials that disintegrate during sampling and testing can be especially difficult to assess. The paper reviews laboratory techniques for mudrock characterization as well as describes geological and engineering geological classification schemes generally used to describe and classify these materials. The value of some of the tests and determinations in the evaluation of a series of mudrock data taken from the literature is presented. Full article

To effectively apply various soil types for embankments, understanding their compaction characteristics is crucial. One crucial factor affecting compaction is suction, which plays a significant role as it is typically performed under unsaturated conditions. Suction varies with soil density, water content, and fines content. This study directly measures suction after soil compaction using the triaxial apparatus, unlike the Soil water characteristic curve (SWCC), assessing its impact on compaction characteristics. Immediate suction measurement after compaction provides apparent suction, resembling on-site conditions with open pore air pressure. Comparing SWCC with apparent suction at each compacted state reveals that suction and air entry value increase with initial density, positively impacting compaction. Notably, apparent suction aligns better with wetting process suction from the SWCC due to added water during specimen preparation. Empirical equations are derived to obtain suction contours across various density and saturation ranges, aiding in understanding suction variations on the compaction curve. Even slight variations in saturation causes noticeable changes in apparent suction during higher compaction efforts, affecting soil compaction characteristics. Therefore, the precise control of saturation control is needed to achieve desired properties of compacted soil, especially at higher compaction efforts and with various soil types. This understanding significantly impacts the mechanical behavior of unsaturated soils. Full article

Severe climatic and environmental conditions warrant the use of stabilization agents in aid of compaction for sustainable improvement in engineering properties of clays. Physicochemical agents are a viable option because they are cost effective, environmentally friendly, and offer improved long-term performance of treated soils. This research developed a fundamental understanding of the clay–water–electrolyte admixtures relations. Based on a comprehensive literature review, the effect of nanomaterials, biopolymers, and geopolymers on the behavior of compacted clays was investigated. It was found that all of these admixtures facilitate the development of an aggregated soil microstructure through unique mechanisms. Biopolymers have the highest water adsorption capacity followed by geopolymers and then by nanomaterials. The effect of admixtures on optimum compaction properties follows a decreasing trend similar to untreated clays (S = 80% ± 20%). The variation of hydraulic conductivity, compression index, and compressive strength are largely within the family of curves identified by typical relationships for compacted clays. These preliminary findings indicate that not all engineering properties are improved to the same level by the different types of physicochemical admixtures. The specific nature of geotechnical engineering (soil type and site conditions) as well as the wide range of admixture types and potential biodegradation of some of the reagents are the major shortcoming of using this class of materials. Full article

Conventional geophysical methods are suitable for estimating the thicknesses of subsoil layers. By combining several geophysical methods, the uncertainties can be assessed. Hence, the reliability of the results increases with a more accurate engineering solution. To estimate the base of an abandoned landfill, we collected data using classical approaches: high-resolution seismic reflection and refraction, with more modern methods including passive surface wave analysis and horizontal-to-vertical spectral ratio (HVSR) measurements. To evaluate the thickness of the landfill, three different datasets were acquired along each of the two seismic lines, and five different processing methods were applied for each of the two arrays. The results of all the classical methods indicate very consistent correlations and mostly converge to clear outcomes. However, since the shear wave velocity of the landfill is relatively low (<150 (m/s)), the uncertainty of the HVSR results is significant. All these methods are engineering-oriented, environmentally friendly, and relatively low-cost. They may be jointly interpreted to better assess uncertainties and therefore enable an efficient solution for environmental or engineering purposes. Full article

The interaction between groundwater and civil engineering works is a key aspect in geotechnical design. In the case of excavations confined in sheet pile walls, steel sheeting, diaphragm walls, cut-off walls, or cofferdams, this design requires the estimation, among other soil mechanics properties, of the groundwater flowing into the excavation (seepage) caused by piezometry depletion. Numerical methods, graphical solutions, and analytical procedures are the methodologies traditionally used to solve this issue, solutions of which require an understanding of basic soil mechanical properties, hydraulic conditions and structure geometry. In this work, the discriminated non-dimensionalization technique is applied to obtain, for the first time, the dimensionless groups that govern the seepage, in anisotropic conditions, in large-scale scenarios where groundwater flow is not conditioned by impervious bedrock or the length of the back of the wall: ${\mathsf{\pi }}_1=\frac{\mathrm{a}}{\mathrm{b}},{\mathsf{\pi }}_2=\frac{{\mathrm{k}}_{\mathrm{x}}{\mathrm{b}}^2}{{\mathrm{k}}_{\mathrm{y}}{\mathrm{c}}^2}$and, ${\mathsf{\pi }}_3=\mathrm{T}/\mathrm{b}$. Numerical simulations are carried out to check the validity of dimensionless groups and to develop three sets of type curves that relate to these groups. Once the physical and geometrical data are known, the seepage (Q), the characteristic depth (T*) and the characteristic horizontal extension (L*) can be directly and easily calculated from these abacuses. The influence of anisotropy on the characteristic lengths is also addressed. Full article

This study focuses on understanding the historical tsunami events in Eastern Indonesia, specifically the Banda Sea region, by extracting information from the limited and challenging-to-interpret historical records. The oldest detailed account of a tsunami in Indonesia dates back to 1674, documented in the book Waerachtigh Verhael Van de Schlickelijcke Aerdbebinge by Rumphius. The study aims to comprehend the primary source of the tsunami and analyze its characteristics to facilitate future tsunami risk reduction. The methodology includes collecting topography and bathymetry data, conducting landslide scenario analysis, employing a two-layer wave propagation model, and performing spectral analysis. The study utilizes comprehensive datasets, investigates potential landslide scenarios, simulates tsunami propagation, and analyzes frequency characteristics using the fast Fourier transform. The 1674 event yielded a runup height of approximately 50–100 m, whereas this study underestimated the actual runup. To illustrate the tsunami wave along the bay’s coastline, a Hovmöller diagram was employed. By analyzing the Hovmöller diagram, the power spectral density was computed, revealing five prominent period bands: 6.96, 5.16, 4.1, 3.75, and 3.36 min. The integration of these components provides a rigorous approach to understanding tsunami dynamics and enhancing risk assessment and mitigation in the study area. Full article

The rainfall-induced landslide disasters in July 2018 in Southwestern Japan yet again exemplified the severity of slope failure-related damage and the need for improvement of early warning systems. The Japanese Meteorological Agency (JMA) uses a method based on a threshold value of soil water index (SWI), a conceptual measurement that represents saturation of slope soil. The current SWI early warning system uses 60-min rainfall data on a 5-km² mesh and does not take into consideration other landslide conditioning factors such as slope angle and geology. This study calculates SWI values during the July 2018 disasters in Kure City (Hiroshima Prefecture) using 1-min XRAIN rainfall data in a 250-m mesh to investigate the relationship between SWI and landslide occurrence. It was found that the SWI threshold of 124 mm used in the JMA early warning system for the area was surpassed in all cells. A new SWI threshold calculation method taking slope angle and geology into consideration and produced with machine learning is proposed, comprising power lines for different geological units at a two-dimensional graph where points located above the threshold line represent landslide risk. It is judged that this method would provide a more accurate early warning system for landslide disasters. Full article

Microbially induced calcite precipitation (MICP) is a green bio-inspired soil solidification technique that depends on the ability of urease-producing bacteria to form calcium carbonate that bonds soil grains and, consequently, improves soil mechanical properties. Meanwhile, different treatment methods have been adopted to tackle the key challenges in achieving effective MICP treatment. This paper proposes the combined method as a new MICP treatment approach, aiming to develop the efficiency of MICP treatment methods and simulate naturally cemented soil. This method combines the premixing, percolation, and submerging MICP methods. The strength outcomes of Portland-cemented and MICP-cemented sand using the percolation and combined methods were compared. For Portland-cemented sand, the UCS values varied from 0.6 MPa to 17.2 MPa, corresponding to cementation levels ranging from 5% to 30%. For MICP-cemented sand, the percolation method yielded UCS values ranging from 0.5 to 0.9 MPa, while the combined method achieved 3.7 MPa. The strength obtained by the combined method is around 3.7 times higher than that of the percolation method. The stiffness of bio-cemented samples varied between 20 and 470 MPa, while for Portland-cemented sand, it ranged from 130 to 1200 MPa. In terms of calcium carbonate distribution, the percolation method exhibited higher concentration at the top of the sample, while the combined method exhibited more precipitation at the top and perimeter, with less concentration in the central bottom region, equivalent to 10% of a half section’s area. Full article

The thermal conductivity of materials is a crucial property with diverse applications, particularly in engineering. Understanding soil thermal conductivity is crucial for designing efficient geothermal systems, predicting soil temperatures, and assessing soil contamination. This paper aimed to predict quartz sand thermal conductivity by using four mathematical models: multiple linear regression (MLR), artificial neural network (ANN), classification and regression random forest (CRRF), and genetic programming (GP). A grey-box AI method, GP, was used for the first time in this topic. Seven inputs affecting thermal conductivity were evaluated in the study, including sand porosity, degree of saturation, coefficient of uniformity, coefficient of curvature, mean particle size, and minimum and maximum void ratios. In predicting thermal conductivity, the MLR model performed poorly, with a coefficient of determination R² = 0.737 and a mean absolute error MAE = 0.300. Both ANN models using the Levenberg–Marquardt algorithm and the Bayesian Regularization (BR) algorithm outperformed the MLR model with an accuracy of R² = 0.916 and an error of MAE = 0.151. In addition, the CRRF model had the best accuracy of R² = 0.993 and MAE = 0.045. In addition, GP showed acceptable performance in predicting sand thermal conductivity. The R² and MAE values of GP were 0.986 and 0.063, respectively. This paper presents the best GP equation for evaluating other databases. Additionally, the porosity and saturation of the sand were found to have the greatest impact on the model results, while coefficients of curvature and uniformity had the least influence. Overall, the results of this study demonstrate that grey-box artificial intelligence models can be used to accurately predict quartz sand thermal conductivity. Full article

This study employs a 3D discrete element method (DEM) to simulate cone penetration tests (CPTs) in granular soils, taking into account the effect of temperature. A coupled thermal mechanical model is developed to allow for heat transfer and storage in the granular materials. The CPT simulations are conducted on granular samples prepared at various temperatures, with the specific heat and velocity of thermal conductivity being identified as two critical factors that influence sample heating time. Additionally, the thermal expansion coefficient is a crucial parameter that is closely related to the porosity of the sample. As the sample temperature increases, the particles expand, resulting in an increase in cone resistance. Full article

This research investigates the convenience of structural identification tools to detect the rocking motion tendency, using as input the structural response to ambient vibrations. The rocking ratio and rocking spectrum are proposed as original tools to highlight the rocking motion and its frequency content. The proposed procedure allows the detection and quantification of rocking using only building vertical motion records in both cases of ambient vibration and earthquake. First, three-dimensional finite element models of reinforced concrete buildings are adopted to simulate the structural response to white noise vibration. Different low- and high-rise buildings are studied, having framed structure and frame–wall system, regular and irregular structure, shallow foundation and underground floors. The structural response obtained numerically is analyzed using different signal processing tools to obtain the dynamic features of buildings, and the rocking motion tendency is identified by comparison with a reference fixed base condition. Then, the reliability of the proposed methodology to detect rocking motion attitude, using only the structural motion, is verified and quantified using the proposed tools. Finally, the same approach is applied to real structural motion records of a high-rise reinforced concrete building. Full article

The current study investigates the capacity of the proposed meta-material layout for the blast protection of above-ground steel pipes against explosions. The philosophy of the meta-material layout’s design is described adequately, and the 1D periodic structures’ theory is adopted for the analytical prediction of the layout’s band-gaps. The special characteristics of the blast loading are explained, and specific time-related parameters are calculated. The layout is tested numerically for nine explosion scenarios of various magnitude via the finite element program ABAQUS, and the CONWEP model is selected for the simulation of the explosions. The results demonstrate a significant reduction in the maximum displacements developed on the pipe’s spring line and crown within a blast loading. This study composes an extension of the author’s previous research on buried steel pipes and surface explosion, advancing now the applicability of the meta-material layouts for the cases of above-ground steel pipes towards explosions and blast hazards. The outer goal is the investigation and the further spreading of the beneficial exploitation of meta-materials concepts for the scope of the pipelines’ effective blast protection, readdressing that this way is a major hazard for this type of structure and a gap in the current literature. Full article

Electromagnetic (EM) waves, traditionally used for purposes such as geophysical characterization, impact properties to be measured. This paper describes the effects of radio frequency (RF) waves on the hydraulic conductivity of glass beads and natural sand. A series of tests was conducted using a customized, rigid-wall, cylindrical permeameter inside a resonant cavity made of Plexiglas covered with electrically conductive transparent films. Constant-head ASTM-D2434 tests were performed to measure the samples’ hydraulic conductivity. RF stimulation was performed using a magnetically coupled loop antenna at various frequencies and input RF-power levels. The hydraulic conductivity of both natural sand and glass-bead samples increased with RF stimulation. Furthermore, the measurement of the electric field component of RF waves was also performed to illustrate the pattern of the electric field, as well as evaluate RF’s impact on the hydraulic conductivity tests. The electric field was numerically simulated and validated against experimentally measured electric fields. A finite-difference numerical model was developed in MATLAB to analyze the seepage flow, which was then validated against the experimental results. An optimization scheme was then used to develop a governing equation for RF’s impact on hydraulic conductivity. Full article

This paper describes the relationship between bedrock depth (D) and site fundamental frequency (f₀) in the Nakdonggang delta region in the southeastern part of the Korean peninsula. We collected borehole logs to confirm the thickness of the sediments and estimated the f₀ at over 200 locations across the delta using the horizontal-to-vertical spectral ratio (HVSR) method. We developed an f₀ map of the study area by spatially interpolating the f₀ values using the Ordinary Kriging method. The bedrock depth in the main delta showed a power-law dependence on the f₀. The derived f₀–D model predicted much shallower bedrock depths compared with similar studies from other parts of the world. This was attributed to the fact that the Nakdonggang delta region is composed of relatively low V_s Holocene sediments. With an f₀ map, the derived model could enable a quick estimation of the bedrock depth, which could help to determine the site class in the Nakdonggang delta region according to the Korean Seismic Design Standard (KDS 17 10 00). Full article