{"title":"Multiscale data-driven modeling of the thermomechanical behavior of granular media with thermal expansion effects","authors":"","doi":"10.1016/j.compgeo.2024.106789","DOIUrl":"10.1016/j.compgeo.2024.106789","url":null,"abstract":"<div><div>A multiscale data-driven (MSDD) methodology is proposed for simulating the thermomechanical behavior of granular materials subjected to thermal expansion. The macroscale is handled using a continuous model based on the Finite Volume Method (FVM), while the microscale response is captured at Representative Volume Elements (RVEs) with the Discrete Element Method (DEM). To significantly reduce the computational cost of the analyses, the microscale DEM computations are not performed online, <span><math><mrow><mi>i</mi><mo>.</mo><mi>e</mi><mo>.</mo></mrow></math></span>, simultaneously with the macroscale FVM ones, as generally done in standard multiscale approaches. Instead, they are performed in advance to create a comprehensive database of RVE solutions under different initial conditions and thermal strains. This dataset is then used to train an Artificial Neural Network (ANN), which serves as a surrogate model for the macroscale solver. The MSDD approach is validated against pure DEM solutions of problems with distinct thermal conditions. Remarkably, we demonstrate that with only three input parameters, namely porosity, fabric, and thermal strain, the surrogate model can predict the microstructure evolution, as well as the updated conductivity and Cauchy stress tensors of the granular assembly. This allows for a generally accurate simulation of transient thermomechanical analyses at a drastically lower computational cost than the pure DEM approach.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142322554","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Integrating machine learning techniques for predicting ground vibration in pile driving activities","authors":"","doi":"10.1016/j.compgeo.2024.106784","DOIUrl":"10.1016/j.compgeo.2024.106784","url":null,"abstract":"<div><div>This study focuses on the assessment of ground vibrations due to pile driving activities. Given the likelihood of excessive vibration due to the driving process, it is imperative to predict vibration levels during the design phase. The primary goal of this work is to integrate machine learning techniques, specifically Extreme Gradient Boosting (XGBoost) and Artificial Neural Networks (ANNs) for real-time vibration prediction. The training dataset was generated using a validated numerical model and the trained models were validated based on experimental results. This validation process highlights the efficiency and accuracy of Extreme Gradient Boosting in predicting the-free-field response of the ground.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142322555","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A cross-anisotropic constitutive framework for modelling gallery intersections excavated in the Callovo-Oxfordian claystone","authors":"","doi":"10.1016/j.compgeo.2024.106761","DOIUrl":"10.1016/j.compgeo.2024.106761","url":null,"abstract":"<div><div>This paper presents the development of a constitutive law tailored for cross-anisotropic rocks and its application on the intersection of two perpendicular unsupported galleries excavated in the the Callovo-Oxfordian (COx) claystone, considered as a potential geological formation to host an underground radioactive waste repository in France. An anisotropic Drucker-Prager elastoplastic law with shear strength hardening is developed to describe the rock matrix behaviour, and the fabric tensor method is used to account for the anisotropy of material strength related to its sedimentary layered structure. The effects of material anisotropy and in-situ stress state anisotropy on the short-term mechanical behaviour of the rock surrounding the intersection area are investigated, and the extent and shape of the plastic zone is presented for different gallery orientations and construction stages. In particular, the results indicate that the distribution of orthoradial stresses and plastic strains are mainly influenced by the anisotropic nature of the rock, leading to a preferential horizontal extension of the plastic zone regardless of the orientation of the galleries. Additionally, the zone affected by the excavation of the intersection is noticeable at a distance of maximum two diameters from the walls of the main gallery.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142319422","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Three-dimensional probabilistic site characterizations using multi-outputs sparse Bayesian learning","authors":"","doi":"10.1016/j.compgeo.2024.106757","DOIUrl":"10.1016/j.compgeo.2024.106757","url":null,"abstract":"<div><div>Three-dimensional probabilistic site characterization is the cornerstone of geotechnical digital transformation, because all engineering projects require an accurate understanding of subsurface geotechnical properties. Soil laboratory testing data or in-situ testing records are often used for data-driven site characterization. However, these site investigation data are often multivariate, uncertain, sparse, and spatially varying. In this paper, the existing sparse Bayesian learning method for three-dimensional (3D) probabilistic site characterization is extended to incorporate multiple soil properties, considering both the three-dimensional spatial variability and the cross-correlation among different soil properties. The proposed three-dimensional multiple-outputs sparse Bayesian learning (3D-MSBL) method is also capable of simulating multiple-correlated conditional random fields in 3D, with the benefit to quantify the statistical uncertainties of soil properties at unexplored locations. The proposed 3D-MSBL method is examined on three case studies. It is shown that the proposed method outperforms the existing single-output SBL method in geotechnical data-driven site characterization, giving more accurate predictions accuracy of soil properties at unexplored locations with smaller statistical uncertainties especially for sparse training data scenario.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142319452","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"In fluence of a firm stratum on probabilistic slope stability analyses with spatially variable undrained soil strength","authors":"","doi":"10.1016/j.compgeo.2024.106782","DOIUrl":"10.1016/j.compgeo.2024.106782","url":null,"abstract":"<div><div>The influence of a firm stratum on the stability of a slope under undrained conditions has long been of interest to geotechnical investigators, which has been studied in a number of previously important works in relation to slope stability analyses without considering soil spatial variability. This paper proposes another look at such a problem in the context of probabilistic slope stability analyses considering soil spatial variability. Here, the random field (RF) is used to simulate the spatially variable undrained soil strength. It is found that under stationary RF and non-stationary RF with the soil strength at the top ground surface (<em>s<sub>u0</sub></em>) larger than 0, the depth of the firm stratum (<em>H<sub>f</sub></em>) has a significant influence on the mean and standard deviation of factor of safety (i.e., <em>μ</em>[FS] and <em>σ</em>[FS], respectively). By contrast, under non-stationary RF with <em>s<sub>u0</sub></em> = 0, <em>H<sub>f</sub></em> has a slight influence on <em>μ</em>[FS], but its influence on <em>σ</em>[FS] is non-negligible. In addition, the autocorrelation distance is found to have an insignificant impact on the influential effect of <em>H<sub>f</sub></em> on <em>μ</em>[FS]. However, for <em>σ</em>[FS], this impact is not negligible. When the autocorrelation distance is smaller, the influence of <em>H<sub>f</sub></em> on <em>σ</em>[FS] would be more significant. Under non-stationary RF, the influence of <em>H<sub>f</sub></em> on <em>σ</em>[FS] would be slight if the autocorrelation distance is large enough. Furthermore, the impacts of slope ratio, <em>s<sub>u0</sub></em>, isotropic and anisotropic features on the influential effects of <em>H<sub>f</sub></em> are also investigated and discussed.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142319423","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Creep-fatigue interaction of rock salt using discrete element simulation","authors":"","doi":"10.1016/j.compgeo.2024.106771","DOIUrl":"10.1016/j.compgeo.2024.106771","url":null,"abstract":"<div><div>The rock surrounding underground salt caverns hosting compressed air energy storage (CAES) systems is subjected to simultaneous creep and fatigue loads. A proper design considering the mechanical response and damage evolution of rock salt under such loading conditions is crucial for ensuring long-term stability. In addition to the conventional laboratory creep-fatigue tests, numerical simulation can serve as a complementary and validating method for unveiling the underlying mechanisms of interaction occurrences. The integration of these methods offers significant value in advancing our understanding of the meso-mechanical properties of rock salt. This study delves deep into the creep-fatigue interactions of rock salt using discrete element method (DEM) simulations. By establishing a DEM model of rock salt, that employs a hybrid contact model of Burger’s and linear parallel bond (LPB) to represent creep and damage behavior, respectively, we unveil the S-shaped evolution of the strain–time curve, encompassing initial, steady, stable damage, and acceleration deformation stages. Our findings highlight that stress intervals significantly reduce fatigue and creep life compared to pure creep and fatigue conditions. Furthermore, we demonstrate the effectiveness of a bilinear cumulative damage rule in describing the creep-fatigue life of rock salt under uniaxial conditions, with longer stress intervals leading to increased cumulative crack numbers and earlier crack initiation times. Ultimately, our numerical salt model exhibits an inclined shear crack and numerous micro-cracks surrounding the macro-shear crack, aligning with fracture modes observed in laboratory tests and advancing our understanding of rock salt behavior in CAES system.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142319453","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A multifractal geometric model for estimating spontaneous imbibition in an unsaturated fractured core-scale network in a low-permeability reservoir","authors":"","doi":"10.1016/j.compgeo.2024.106746","DOIUrl":"10.1016/j.compgeo.2024.106746","url":null,"abstract":"<div><div>Several mathematical models have been established for describing the saturated seepage of fracture networks in low-permeability reservoirs, but there is still a lack of systematic research on unsaturated seepage models. Therefore, by combining these single connected fractures into the physical geometry of a parallel plate, a comprehensive spontaneous imbibition model at the fracture and core scales considering a rough wall geometry, streamline tortuosity, and heterogeneity of the fracture size distribution is derived by applying the classical cubic law describing flat flow and fractal geometry theory. Compared with existing models, our proposed model is verified by using experimental data on the spontaneous imbibition of a single fracture and fractured core available from the literature. Based on the established model, the mechanism through which various geometric elements control spontaneous imbibition is analysed at the fracture and core scales. The maximum fracture aperture plays an important role in controlling the wettability evolution behaviour at the core scale, and the fractal dimensions of fracture tortuosity and roughness are more sensitive to spontaneous imbibition in fractured cores than is the fractal dimension of the fracture size distribution. These research results are useful for evaluating the balance between reservoir damage and fracturing effects.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142315049","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A constitutive model considering the interaction between evolution of microstructure and hydro-mechanical behaviour of unsaturated soils","authors":"","doi":"10.1016/j.compgeo.2024.106775","DOIUrl":"10.1016/j.compgeo.2024.106775","url":null,"abstract":"<div><div>The evolution of microstructure induced by loading and unloading has a significant impact on the hydro-mechanical behaviour of soils, including volume change, shear strength, water retention and permeability. In this paper, a constitutive model based on the evolution of microstructure is established building on the approach of an existing mechanistic model. In this model, the evolution of microstructure is represented via changes in the pore size distribution (PSD) and assumed to be related solely to the change of void ratio induced by loading and unloading. A PSD-dependent Bishop’s effective stress coefficient χ*, which represents the coupled impact of PSD evolution on hydro-mechanical behaviour of soils, is used to replace the Bishop’s effective stress coefficient χ. The model can reproduce and predict the hydro-mechanical behaviour and evolution of microstructure and their interaction within a unified framework. It also has potential in studying the soil-water characteristic curve and multi-field-coupling of soils. Model response and sensitivity analysis are reported based on idealized parameters to give a primary evaluation on the model’s performance and feasibility of using PSDs from mercury intrusion porosimetry. It is found that whilst the model is sensitive to parameters representing inter-aggregate pore size distributions it can be satisfactorily applied to represent the hydro-mechanical behaviour and microstructural evolution of unsaturated soils.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0266352X24007146/pdfft?md5=2c48de659303c27e7e3ff87f0dc5928b&pid=1-s2.0-S0266352X24007146-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142315051","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Research on the static liquefaction failure of upstream tailings dams under continuous discharge conditions","authors":"","doi":"10.1016/j.compgeo.2024.106785","DOIUrl":"10.1016/j.compgeo.2024.106785","url":null,"abstract":"<div><div>Upstream tailings dams are extensively employed owing to their straightforward and economical construction methodologies worldwide. Nevertheless, they are susceptible to static liquefaction failures, wherein sustained tailings discharge (monotonic loading) stands as one of the primary triggering conditions for static liquefaction. This study centers on the Erlichong tailings dam and employs a multifaceted approach, combining field investigations, laboratory experiments, and fluid–solid coupling discrete element method. It mainly aims to ascertain the influence of tailings discharge rates, fine-grained interlayers, and drainage facilities on static liquefaction failure within the tailings dam. It is found that, higher discharge rates lead to accelerated accumulation of excess pore water pressure, rendering the tailings dam increasingly vulnerable to static liquefaction failure. Meanwhile, although in-situ mixed tailings exhibit greater susceptibility to static liquefaction than fine tailings with same initial void ratio under completely undrained conditions, the presence of fine-grained interlayers diminishes the vertical permeability of the tailings deposit in practical scenarios. Undrained conditions are established for the mixed tailings confined within the interlayers. Consequently, tailings dams containing fine-grained interlayers demonstrate heightened vulnerability to static liquefaction failure in comparison to homogeneous tailings dams. Furthermore, the incorporation of drainage facilities serves to mitigate the risk of static liquefaction, with the vertical-horizontal combined drainage scheme emerging as the optimal strategy for tailings dams with fine-grained interlayers. In summary, this study furnishes theoretical insights aimed at ensuring the secure operation of tailings dams.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142315050","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An anisotropic elastoplastic strong discontinuity model for shear failure in anisotropic rock masses","authors":"","doi":"10.1016/j.compgeo.2024.106762","DOIUrl":"10.1016/j.compgeo.2024.106762","url":null,"abstract":"<div><div>In this paper, we propose a noval anisotropic elastoplastic strong discontinuity-FEM (SD-FEM) for analyzing the complete progressive shear failure process in anisotropic rock masses. The deformation property preceding slip is described with an elastoplastic formulation incorporating the microstructure tensor approach. Anisotropic discontinuous bifurcation analysis is conducted to judge the initiation conditions and propagation direction of slip lines. Furthermore, a derived anisotropic stress-displacement relation on the discontinuity is utilized to describe the post-failure response asscociated with slip. Two numerical examples, namely the uniaxial compression test of anisotropic rock masses and the loading problem of the anisotropic rock slope, are used to demonstrate the remarkable capabilities of the anisotropic elastoplastic SD-FEM model. It is illustrated that this model can not only reflect the anisotropic mechanical characteristics of rock masses but also accurately simulate the complete progressive failure process, spanning from uniform deformation to sliding failure. Notably, it is observed that the magnitude of slip (i.e., discontinuous displacement) exhibits a linear increase with the vertical load, highlighting the elastic-brittle deformation characteristics of rock masses. Moreover, the eigenvalues of the global stiffness matrix remain positive even in the softening stage, which enables the numerical calculation to proceed and ensures the mesh-independent numerical solutions, indicating that the anisotropic SD-FEM can regularize the ill-posedness of the boundary value problem when strain softening occurs.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-09-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142315048","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}