Pub Date : 2025-01-16DOI: 10.1016/j.compgeo.2025.107067
Xuening Ma, Fan Bai
In the design of reinforced structure, the shear strength of the reinforcement-soil interface is a key indicator of reinforcement efficacy, which is significantly influenced by the geometric characteristics of geogrids. At present, the research on how the geogrid mesh shape and aperture ratio of geogrid affect the interface shear response is not sufficient. In this study, the traditional biaxial geogrid and three new types of geogrids were selected, through the pull-out test and coupled discrete element method (DEM) and finite difference method (FDM), the influences of the arrangement of transverse ribs and the change of aperture ratio on the effect and mechanism of the reinforced-soil interface were explored. The optimal mesh shape suitable for reinforced aeolian sand and the optimal aperture size ratio under various shapes of geogrids were analyzed. Based on the analysis of the displacement, shear band distribution and porosity evolution of aeolian sand reinforced by different geogrids during the drawing process, the reinforcement effect and mechanism were revealed from the mesoscopic level. The research results provide important reference for the design of new geogrids and the reinforcement design of desert roadbeds.
{"title":"The role of geogrid aperture shape and size in strengthening aeolian sands: Insights from a coupled DEM-FDM approach","authors":"Xuening Ma, Fan Bai","doi":"10.1016/j.compgeo.2025.107067","DOIUrl":"10.1016/j.compgeo.2025.107067","url":null,"abstract":"<div><div>In the design of reinforced structure, the shear strength of the reinforcement-soil interface is a key indicator of reinforcement efficacy, which is significantly influenced by the geometric characteristics of geogrids. At present, the research on how the geogrid mesh shape and aperture ratio of geogrid affect the interface shear response is not sufficient. In this study, the traditional biaxial geogrid and three new types of geogrids were selected, through the pull-out test and coupled discrete element method (DEM) and finite difference method (FDM), the influences of the arrangement of transverse ribs and the change of aperture ratio on the effect and mechanism of the reinforced-soil interface were explored. The optimal mesh shape suitable for reinforced aeolian sand and the optimal aperture size ratio under various shapes of geogrids were analyzed. Based on the analysis of the displacement, shear band distribution and porosity evolution of aeolian sand reinforced by different geogrids during the drawing process, the reinforcement effect and mechanism were revealed from the mesoscopic level. The research results provide important reference for the design of new geogrids and the reinforcement design of desert roadbeds.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107067"},"PeriodicalIF":5.3,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104938","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}
Pub Date : 2025-01-16DOI: 10.1016/j.compgeo.2024.107040
M.Q. Xu , Z.X. Yang
Understanding the response of sand to complex loading conditions is vital for practical geotechnical engineering. Circular rotational shear is a special stress path where the magnitudes of three principal stresses vary following a circular stress trajectory in the π-plane with their directions fixed. Although experimental studies under such stress paths are limited, the discrete element method appears to be an appealing approach to examine the response of granular materials to varying complex loading paths in numerical “virtual” tests. This study presents comprehensive numerical simulations of granular samples subjected to a circular stress path under varying conditions, including samples prepared with different bedding-plane angles and densities and subjected to different stress ratios. Both macroscopic and microscopic behaviors are presented and interpreted. A contact-normal-based fabric tensor is adopted in a detailed analysis to measure the internal structure of the granular assembly. The fabric, strain, and strain increment tensors are decomposed with respect to the stress tensor, and the evolutions of these components are presented along with the key influential factors. The results obtained in this study provide useful physical insight for the development of constitutive models for granular soils under general loading conditions.
{"title":"Three-dimensional discrete element simulations of inherently anisotropic granular materials subjected to circular rotational shear stress path","authors":"M.Q. Xu , Z.X. Yang","doi":"10.1016/j.compgeo.2024.107040","DOIUrl":"10.1016/j.compgeo.2024.107040","url":null,"abstract":"<div><div>Understanding the response of sand to complex loading conditions is vital for practical geotechnical engineering. Circular rotational shear is a special stress path where the magnitudes of three principal stresses vary following a circular stress trajectory in the π-plane with their directions fixed. Although experimental studies under such stress paths are limited, the discrete element method appears to be an appealing approach to examine the response of granular materials to varying complex loading paths in numerical “virtual” tests. This study presents comprehensive numerical simulations of granular samples subjected to a circular stress path under varying conditions, including samples prepared with different bedding-plane angles and densities and subjected to different stress ratios. Both macroscopic and microscopic behaviors are presented and interpreted. A contact-normal-based fabric tensor is adopted in a detailed analysis to measure the internal structure of the granular assembly. The fabric, strain, and strain increment tensors are decomposed with respect to the stress tensor, and the evolutions of these components are presented along with the key influential factors. The results obtained in this study provide useful physical insight for the development of constitutive models for granular soils under general loading conditions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107040"},"PeriodicalIF":5.3,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172637","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}
Pub Date : 2025-01-16DOI: 10.1016/j.compgeo.2025.107061
Meng Fan , Dong Su , Ning Zhang , Guojun Cai , Xiangsheng Chen
The bubble packing algorithm has been widely used in previous DEM studies to generate clump models that approximate concave particle shapes. However, the selection of optimal parameters—specifically, the ratio and distance —is often subjective or based solely on volume ratios. This study offers a more comprehensive quantitative analysis by evaluating both the geometrical and morphological accuracy of clump assemblies generated using various combinations of and . The analysis focuses on typical geotechnical particles, including Fujian sand, ballast, and cobble particles. The study establishes correlations for the approximation degrees of each geometrical and morphological parameter and the average number of sub-spheres per clump in relation to and . Using a 5 % threshold for approximation accuracy, the optimal parameter combination of and is identified and verified for each type of particles. The findings of this paper can provide a useful reference for future studies that utilize the bubble packing algorithm to generate clump models for approximating geotechnical particle shapes.
{"title":"Optimizing parameter combinations for clump models enabled by the bubble packing algorithm: Insights from geometrical and morphological approximation of typical geotechnical particles","authors":"Meng Fan , Dong Su , Ning Zhang , Guojun Cai , Xiangsheng Chen","doi":"10.1016/j.compgeo.2025.107061","DOIUrl":"10.1016/j.compgeo.2025.107061","url":null,"abstract":"<div><div>The bubble packing algorithm has been widely used in previous DEM studies to generate clump models that approximate concave particle shapes. However, the selection of optimal parameters—specifically, the ratio <span><math><msub><mi>ρ</mi><mi>r</mi></msub></math></span> and distance <span><math><mi>Φ</mi></math></span>—is often subjective or based solely on volume ratios. This study offers a more comprehensive quantitative analysis by evaluating both the geometrical and morphological accuracy of clump assemblies generated using various combinations of <span><math><msub><mi>ρ</mi><mi>r</mi></msub></math></span> and <span><math><mi>Φ</mi></math></span>. The analysis focuses on typical geotechnical particles, including Fujian sand, ballast, and cobble particles. The study establishes correlations for the approximation degrees of each geometrical and morphological parameter and the average number of sub-spheres per clump in relation to <span><math><msub><mi>ρ</mi><mi>r</mi></msub></math></span> and <span><math><mi>Φ</mi></math></span>. Using a 5 % threshold for approximation accuracy, the optimal parameter combination of <span><math><msub><mi>ρ</mi><mi>r</mi></msub></math></span> and <span><math><mi>Φ</mi></math></span> is identified and verified for each type of particles. The findings of this paper can provide a useful reference for future studies that utilize the bubble packing algorithm to generate clump models for approximating geotechnical particle shapes.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107061"},"PeriodicalIF":5.3,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172578","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}
Pub Date : 2025-01-15DOI: 10.1016/j.compgeo.2025.107065
Wei Qin , Chen Ye , Jiayu Gao , Guoliang Dai , Dong Wang , Youkou Dong
This study investigates the pore water pressure (PWP) behavior of soil around large-diameter open-ended thin-walled piles (LOTPs) during impact driving using a large deformation finite-element method. A comparative analysis of the PWP accumulation curves of the soil inside, outside, and below the LOTP tips with different diameters and wall thicknesses during impact driving is conducted under the same hammering solution. The PWP development is dependent on the absolute distance from the pile surface to the location of the soil and the dimensions of the LOTP. The excess pore water pressure (EPWP) accumulates and gradually dissipates, and its level decreases with increasing pile diameter. However, a negative excess pore water pressure (Ne-EPWP) is identified during hammering. Based on the above findings and analyses, a PWP prediction equation for LOTP during driving is proposed, and the predicted curves are compared with the numerical results. The influence of PWP accumulation after penetration of 2d (d is the LOTP internal diameter) does not increase significantly. This equation can provide the initial distribution field of PWP in saturated clay for LOTPs, thereby facilitating pile drivability analyses.
{"title":"Pore water pressure of clay soil around large-diameter open-ended thin-walled pile (LOTP) during impact penetration","authors":"Wei Qin , Chen Ye , Jiayu Gao , Guoliang Dai , Dong Wang , Youkou Dong","doi":"10.1016/j.compgeo.2025.107065","DOIUrl":"10.1016/j.compgeo.2025.107065","url":null,"abstract":"<div><div>This study investigates the pore water pressure (PWP) behavior of soil around large-diameter open-ended thin-walled piles (LOTPs) during impact driving using a large deformation finite-element method. A comparative analysis of the PWP accumulation curves of the soil inside, outside, and below the LOTP tips with different diameters and wall thicknesses during impact driving is conducted under the same hammering solution. The PWP development is dependent on the absolute distance from the pile surface to the location of the soil and the dimensions of the LOTP. The excess pore water pressure (EPWP) accumulates and gradually dissipates, and its level decreases with increasing pile diameter. However, a negative excess pore water pressure (Ne-EPWP) is identified during hammering. Based on the above findings and analyses, a PWP prediction equation for LOTP during driving is proposed, and the predicted curves are compared with the numerical results. The influence of PWP accumulation after penetration of 2<em>d</em> (<em>d</em> is the LOTP internal diameter) does not increase significantly. This equation can provide the initial distribution field of PWP in saturated clay for LOTPs, thereby facilitating pile drivability analyses.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107065"},"PeriodicalIF":5.3,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104495","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}
Pub Date : 2025-01-15DOI: 10.1016/j.compgeo.2025.107063
Chaozhe Zhang , Songyu Liu , Dingwen Zhang , Zhen Wang , Hongjiang Li , Guokai Zhang
The development of innovative and cost-effective pile-supported embankment technologies is an intrinsic demand driving the advancement of soft ground improvement techniques. A significant innovation within this approach is the penetrated precast concrete pile reinforced with cemented soil (PCCS) combined with partially penetrated deep mixing (DM) columns for layered ground treatment, in which the PCCS is created by embedding a precast concrete (cored) pile within a DM column. This study aims to develop a numerical unit model for this type of pile-supported embankment, denoted as a PCCS multi-pile composite foundation-supported embankment, providing insights into the soil arching effect within the embankment and elucidating the variation patterns of vertical stress and the horizontal earth pressure coefficient along the embankment height. Based on Hewlett’s hemispherical arch model, a refined composite soil-arching analytical model is proposed specifically for the multi-pile composite foundation-supported embankment. The findings reveal that the critical height of the embankment is 1.05 times the net pile spacing, with embankment height showing minimal influence on the height of internal and external arches. In contrast, an increase in the embankment filling modulus significantly affects the fan-shaped arches between adjacent piles, while its impact on spherical arches remains limited. The stress ratios between the PCCS-soil and DM column-soil derived from the three-dimensional composite arch model differ from the measured results by 7.3 % and 8.5 %, respectively, confirming the model’s accuracy.
{"title":"Analytical and numerical investigation of soil arching effect of the composite foundation reinforced by penetrated PCCSs and partially penetrated DM columns under embankment loading","authors":"Chaozhe Zhang , Songyu Liu , Dingwen Zhang , Zhen Wang , Hongjiang Li , Guokai Zhang","doi":"10.1016/j.compgeo.2025.107063","DOIUrl":"10.1016/j.compgeo.2025.107063","url":null,"abstract":"<div><div>The development of innovative and cost-effective pile-supported embankment technologies is an intrinsic demand driving the advancement of soft ground improvement techniques. A significant innovation within this approach is the penetrated precast concrete pile reinforced with cemented soil (PCCS) combined with partially penetrated deep mixing (DM) columns for layered ground treatment, in which the PCCS is created by embedding a precast concrete (cored) pile within a DM column. This study aims to develop a numerical unit model for this type of pile-supported embankment, denoted as a PCCS multi-pile composite foundation-supported embankment, providing insights into the soil arching effect within the embankment and elucidating the variation patterns of vertical stress and the horizontal earth pressure coefficient along the embankment height. Based on Hewlett’s hemispherical arch model, a refined composite soil-arching analytical model is proposed specifically for the multi-pile composite foundation-supported embankment. The findings reveal that the critical height of the embankment is 1.05 times the net pile spacing, with embankment height showing minimal influence on the height of internal and external arches. In contrast, an increase in the embankment filling modulus significantly affects the fan-shaped arches between adjacent piles, while its impact on spherical arches remains limited. The stress ratios between the PCCS-soil and DM column-soil derived from the three-dimensional composite arch model differ from the measured results by 7.3 % and 8.5 %, respectively, confirming the model’s accuracy.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107063"},"PeriodicalIF":5.3,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104496","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}
Pub Date : 2025-01-15DOI: 10.1016/j.compgeo.2024.107052
Jian Gong , Fei Wang , Liangbin Deng , Jiayan Nie , Hai Xu , Yi Zheng , Jie Jiang , Xiaoduo Ou , Xianwei Pang
The traditional stress-force-fabric relationship (SFF relationship) assumes that contact forces and branch vectors are independent, linking between the macro- and micro-scale behaviors of granular materials. However, both experimental results and numerical simulations indicate that this assumption does not hold for granular materials with high particle size non-uniformity. This study investigates the effect of non-uniformity on the SFF relationship for granular materials via 3D discrete element method (DEM). A series of drained triaxial tests was conducted on dense granular materials with continuous and gap gradations, varying the non-uniformity coefficient (Cu). The results show that as Cu increase, the stress ratios predicted by the traditional SFF relationship gradually deviate from that DEM data. This deviation is attributed to the correlation between contact forces and branch vectors, which can be characterized using a lower-order Fourier series expansion. Based on the traditional SFF relationship and incorporating this correlation, a modified SFF relationship is proposed. Verifications show that stress ratios from the modified SFF relationship align well with the DEM data. Finally, an anisotropic analysis provides insights into the microscopic mechanisms underlying the dependence of shear strength on Cu.
{"title":"Investigation of the effect of particle size non-uniformity on the stress-force-fabric relationship for granular materials","authors":"Jian Gong , Fei Wang , Liangbin Deng , Jiayan Nie , Hai Xu , Yi Zheng , Jie Jiang , Xiaoduo Ou , Xianwei Pang","doi":"10.1016/j.compgeo.2024.107052","DOIUrl":"10.1016/j.compgeo.2024.107052","url":null,"abstract":"<div><div>The traditional stress-force-fabric relationship (SFF relationship) assumes that contact forces and branch vectors are independent, linking between the macro- and micro-scale behaviors of granular materials. However, both experimental results and numerical simulations indicate that this assumption does not hold for granular materials with high particle size non-uniformity. This study investigates the effect of non-uniformity on the SFF relationship for granular materials via 3D discrete element method (DEM). A series of drained triaxial tests was conducted on dense granular materials with continuous and gap gradations, varying the non-uniformity coefficient (C<sub>u</sub>). The results show that as C<sub>u</sub> increase, the stress ratios predicted by the traditional SFF relationship gradually deviate from that DEM data. This deviation is attributed to the correlation between contact forces and branch vectors, which can be characterized using a lower-order Fourier series expansion. Based on the traditional SFF relationship and incorporating this correlation, a modified SFF relationship is proposed. Verifications show that stress ratios from the modified SFF relationship align well with the DEM data. Finally, an anisotropic analysis provides insights into the microscopic mechanisms underlying the dependence of shear strength on C<sub>u</sub>.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107052"},"PeriodicalIF":5.3,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104937","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}
Pub Date : 2025-01-14DOI: 10.1016/j.compgeo.2024.107025
Maozhu Peng , Huangcheng Fang , Yinghui Tian
An efficient soil-structure interaction interface is developed in this paper as a gap that can be filled with water to model the whole process of structure separation from soil (i.e. losing effective mechanical contact with suction sustained). This hydro-mechanical interface is formulated as node-to-face contact, using lubrication theory to describe the gap water motion, while both the structure and the soil are generalised as porous media obeying mass and momentum balances. Continuities of water pressure, flux, and stresses are enforced at the boundaries between the structure, gap, and soil. Incremental forms of governing equations are derived to accommodate material nonlinearity, and finite element formulations are detailed. The interface is validated by three examples. The first two consider a two-layer soil subjected to compression and tension, respectively, with the interface embedded in the middle. The numerical results agree well with the derived analytical solutions using Laplace transform. The third example models lifting a square footing off Modified Cam-Clay seabed at varying rates. Numerical results agree well with available centrifuge data.
{"title":"An H-M contact model to simulate soil-structure interactions under tension in offshore environments","authors":"Maozhu Peng , Huangcheng Fang , Yinghui Tian","doi":"10.1016/j.compgeo.2024.107025","DOIUrl":"10.1016/j.compgeo.2024.107025","url":null,"abstract":"<div><div>An efficient soil-structure interaction interface is developed in this paper as a gap that can be filled with water to model the whole process of structure separation from soil (i.e. losing effective mechanical contact with suction sustained). This hydro-mechanical interface is formulated as node-to-face contact, using lubrication theory to describe the gap water motion, while both the structure and the soil are generalised as porous media obeying mass and momentum balances. Continuities of water pressure, flux, and stresses are enforced at the boundaries between the structure, gap, and soil. Incremental forms of governing equations are derived to accommodate material nonlinearity, and finite element formulations are detailed. The interface is validated by three examples. The first two consider a two-layer soil subjected to compression and tension, respectively, with the interface embedded in the middle. The numerical results agree well with the derived analytical solutions using Laplace transform. The third example models lifting a square footing off Modified Cam-Clay seabed at varying rates. Numerical results agree well with available centrifuge data.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107025"},"PeriodicalIF":5.3,"publicationDate":"2025-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172636","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}
Pub Date : 2025-01-14DOI: 10.1016/j.compgeo.2024.107047
Zhaonan Wang , Gang Wang , Fenglei Han
In this study, we develop a novel framework by coupling the Updated Lagrangian Material Point Method (ULMPM) and the Total Lagrangian Material Point Method (TLMPM) to simulate the dynamic interaction between large-deformation soils and solid structures. The key contribution lies in enhancing the stabilized performance of TLMPM under nearly incompressible conditions through the implementation of a mixed-formulation scheme. Within the mixed-formulation scheme, linear momentum and the deformation gradient are treated as the primary unknowns, and the Petrov-Galerkin method is employed to derive the weak form of the governing equations. To explicitly characterize the contact behavior between ULMPM and TLMPM, a particle-based contact algorithm is introduced, where the contact forces are treated as external forces acting on different particle pairs. The accuracy and robustness of the proposed scheme are validated through the modeling of various quasi-static and dynamic benchmarks. By comparing the numerical results with those obtained from other methods and experimental data, the effectiveness of the framework in addressing nearly incompressible materials and soil-structure interaction problems is demonstrated.
{"title":"A mixed-formulated total Lagrangian material point method for dynamic interaction between structure and soil","authors":"Zhaonan Wang , Gang Wang , Fenglei Han","doi":"10.1016/j.compgeo.2024.107047","DOIUrl":"10.1016/j.compgeo.2024.107047","url":null,"abstract":"<div><div>In this study, we develop a novel framework by coupling the Updated Lagrangian Material Point Method (ULMPM) and the Total Lagrangian Material Point Method (TLMPM) to simulate the dynamic interaction between large-deformation soils and solid structures. The key contribution lies in enhancing the stabilized performance of TLMPM under nearly incompressible conditions through the implementation of a mixed-formulation scheme. Within the mixed-formulation scheme, linear momentum and the deformation gradient are treated as the primary unknowns, and the Petrov-Galerkin method is employed to derive the weak form of the governing equations. To explicitly characterize the contact behavior between ULMPM and TLMPM, a particle-based contact algorithm is introduced, where the contact forces are treated as external forces acting on different particle pairs. The accuracy and robustness of the proposed scheme are validated through the modeling of various quasi-static and dynamic benchmarks. By comparing the numerical results with those obtained from other methods and experimental data, the effectiveness of the framework in addressing nearly incompressible materials and soil-structure interaction problems is demonstrated.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107047"},"PeriodicalIF":5.3,"publicationDate":"2025-01-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172635","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}
Pub Date : 2025-01-13DOI: 10.1016/j.compgeo.2025.107055
Ze-Wei Ke , Sheng-Jie Wei , Shi-Yuan Yao , Si Chen , Yun-Min Chen , Yu-Chao Li
Solving the advection–diffusion equation (ADE) for contaminant transport in soil (forward problem) is of crucial importance in many environmental engineering topics, such as assessment of site contamination risks and design of engineered barriers. Although numerical techniques are widely used to solve the ADEs, they are not skilled at addressing inverse problems, such as identifying unknown parameters in the equations based on measurement data, especially when data are sparse or corrupted with noise. In this paper, forward and inverse problems of the contaminant transport in soils are solved using the newly developed physics-informed neural networks (PINN) incorporated with pre-training strategy, uncertainty quantification and domain decomposition method. Four cases are analyzed in detail to demonstrate the capability of the proposed approach. The results show that: (1) for forward problems, the proposed approach can provide spatio-temporal concentration distributions in a high agreement with analytical or numerical solutions, even for the two-dimensional case with layered soils; (2) for inverse problems, unknown parameters in the ADE can be accurately identified by the proposed approach based on a small amount of measured data, even for the case with two-parameter nonlinear adsorption model; (3) pre-training strategy can significantly enhance the training efficiency and prediction accuracy of PINN; (4) the uncertainty of the results can be effectively quantified by the proposed approach through incorporating the latent variable; and (5) the robustness against measured data noise can be ensured by the proposed approach. The proposed approach has the penitential to address contaminant transport problems under coupled multi-physics with multi-fidelity data.
{"title":"Pre-trained Physics-Informed Neural Networks for Analysis of Contaminant Transport in Soils","authors":"Ze-Wei Ke , Sheng-Jie Wei , Shi-Yuan Yao , Si Chen , Yun-Min Chen , Yu-Chao Li","doi":"10.1016/j.compgeo.2025.107055","DOIUrl":"10.1016/j.compgeo.2025.107055","url":null,"abstract":"<div><div>Solving the advection–diffusion equation (ADE) for contaminant transport in soil (forward problem) is of crucial importance in many environmental engineering topics, such as assessment of site contamination risks and design of engineered barriers. Although numerical techniques are widely used to solve the ADEs, they are not skilled at addressing inverse problems, such as identifying unknown parameters in the equations based on measurement data, especially when data are sparse or corrupted with noise. In this paper, forward and inverse problems of the contaminant transport in soils are solved using the newly developed physics-informed neural networks (PINN) incorporated with pre-training strategy, uncertainty quantification and domain decomposition method. Four cases are analyzed in detail to demonstrate the capability of the proposed approach. The results show that: (1) for forward problems, the proposed approach can provide spatio-temporal concentration distributions in a high agreement with analytical or numerical solutions, even for the two-dimensional case with layered soils; (2) for inverse problems, unknown parameters in the ADE can be accurately identified by the proposed approach based on a small amount of measured data, even for the case with two-parameter nonlinear adsorption model; (3) pre-training strategy can significantly enhance the training efficiency and prediction accuracy of PINN; (4) the uncertainty of the results can be effectively quantified by the proposed approach through incorporating the latent variable; and (5) the robustness against measured data noise can be ensured by the proposed approach. The proposed approach has the penitential to address contaminant transport problems under coupled multi-physics with multi-fidelity data.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107055"},"PeriodicalIF":5.3,"publicationDate":"2025-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143172634","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}
Pub Date : 2025-01-13DOI: 10.1016/j.compgeo.2025.107071
Xiaopeng Zhang , Kaiheng Hu , Hao Li , Shuang Liu , Junyu Liu
Dynamic numerical simulation of large-scale landslides with material point method (MPM) is becoming increasingly popular in the geotechnical engineering and natural hazard community. However, efficient and accurate simulation of the interaction between landslides and complex terrain which the sliding mass moving remains a significant challenge. This study uses a fixed Triangulated Irregular Network (TIN) model from a digital elevation model to represent complex terrain, and models the interaction between sliding mass and terrain with a contact algorithm between TIN meshes and material points. A novel local neighborhood search (LNS) scheme with bucket-sorting algorithm is proposed to improve the particle-to-face contact detection efficiency. The new search algorithm significantly enhances both the efficiency and accuracy of matching detection between the triangular mesh and material points. The correctness and capability in handling complex dynamic problems are demonstrated through two numerical cases. The catastrophic Xinmo landslide occurred in China in 2017 is also simulated to validate the method’s capability in large-scale landslide simulation. The results demonstrate that the proposed method is both accurate and computationally efficient. Furthermore, its adaptability and simplicity allow for seamless integration into other high-performance MPM frameworks, facilitating the simulation of large-scale gravity-driven flows on complex terrain.
{"title":"Efficient local neighborhood search scheme for three-dimensional MPM modeling of large-scale landslides on complex terrain","authors":"Xiaopeng Zhang , Kaiheng Hu , Hao Li , Shuang Liu , Junyu Liu","doi":"10.1016/j.compgeo.2025.107071","DOIUrl":"10.1016/j.compgeo.2025.107071","url":null,"abstract":"<div><div>Dynamic numerical simulation of large-scale landslides with material point method (MPM) is becoming increasingly popular in the geotechnical engineering and natural hazard community. However, efficient and accurate simulation of the interaction between landslides and complex terrain which the sliding mass moving remains a significant challenge. This study uses a fixed Triangulated Irregular Network (TIN) model from a digital elevation model to represent complex terrain, and models the interaction between sliding mass and terrain with a contact algorithm between TIN meshes and material points. A novel local neighborhood search (LNS) scheme with bucket-sorting algorithm is proposed to improve the particle-to-face contact detection efficiency. The new search algorithm significantly enhances both the efficiency and accuracy of matching detection between the triangular mesh and material points. The correctness and capability in handling complex dynamic problems are demonstrated through two numerical cases. The catastrophic Xinmo landslide occurred in China in 2017 is also simulated to validate the method’s capability in large-scale landslide simulation. The results demonstrate that the proposed method is both accurate and computationally efficient. Furthermore, its adaptability and simplicity allow for seamless integration into other high-performance MPM frameworks, facilitating the simulation of large-scale gravity-driven flows on complex terrain.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"180 ","pages":"Article 107071"},"PeriodicalIF":5.3,"publicationDate":"2025-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143104943","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}