This study presents a fully resolved numerical method coupling framework for capturing the interaction between debris flows and flexible barriers. The proposed coupling method combines a simplified flexible wire representation, smoothed particle hydrodynamics (SPH), and discrete element method (DEM). Specifically, the flexible wire representation is used to represent a flexible barrier that intercepts debris flows. The coupled SPH-DEM simulates a two-phase debris flow containing arbitrarily shaped boulders. Additionally, a modified dynamic boundary method (mDBC) is employed to couple SPH with DEM. These three methods are strongly coupled to simulate the behavior of boulders in debris flows and improve the efficiency and accuracy of the interaction between flexible barriers and debris flows. The proposed coupling framework is evaluated through a series of verification instances, which are consistent with experimental observations. Subsequently, the impact dynamics of the two-phase debris flow on the flexible barrier is modeled and analyzed. Overall, the framework holds great potential for the design and optimization of flexible barriers for debris flow mitigation in engineering practice.
{"title":"A fully resolved SPH-DEM for simulation of debris flows with arbitrary particle shapes impacting flexible barriers","authors":"Hao Xiong, Mengjie Hao, Debo Zhao, Xiaoxuan Gan, Zhen-Yu Yin, Xiangsheng Chen","doi":"10.1007/s11440-025-02544-8","DOIUrl":"10.1007/s11440-025-02544-8","url":null,"abstract":"<div><p>This study presents a fully resolved numerical method coupling framework for capturing the interaction between debris flows and flexible barriers. The proposed coupling method combines a simplified flexible wire representation, smoothed particle hydrodynamics (SPH), and discrete element method (DEM). Specifically, the flexible wire representation is used to represent a flexible barrier that intercepts debris flows. The coupled SPH-DEM simulates a two-phase debris flow containing arbitrarily shaped boulders. Additionally, a modified dynamic boundary method (mDBC) is employed to couple SPH with DEM. These three methods are strongly coupled to simulate the behavior of boulders in debris flows and improve the efficiency and accuracy of the interaction between flexible barriers and debris flows. The proposed coupling framework is evaluated through a series of verification instances, which are consistent with experimental observations. Subsequently, the impact dynamics of the two-phase debris flow on the flexible barrier is modeled and analyzed. Overall, the framework holds great potential for the design and optimization of flexible barriers for debris flow mitigation in engineering practice.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1403 - 1430"},"PeriodicalIF":5.6,"publicationDate":"2025-02-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143571021","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-27DOI: 10.1007/s11440-025-02540-y
Qiusheng Wang, Songnan Qiu, Hong Zheng, Ruitao Zhang
O’Kelly et al. critique the Wang et al. undrained shear strength prediction model paper in terms of its utility and accuracy. Both criticisms suffer from the miscomprehension of the model comparison method. Utility of the three-parameter liquidity index (IL)—undrained shear strength (Su) models has been enhanced through optimizing some empirical equation.
{"title":"Response to Comment on “Undrained shear strength prediction of clays using liquidity index”","authors":"Qiusheng Wang, Songnan Qiu, Hong Zheng, Ruitao Zhang","doi":"10.1007/s11440-025-02540-y","DOIUrl":"10.1007/s11440-025-02540-y","url":null,"abstract":"<div><p>O’Kelly et al. critique the Wang et al. undrained shear strength prediction model paper in terms of its utility and accuracy. Both criticisms suffer from the miscomprehension of the model comparison method. Utility of the three-parameter liquidity index (<i>I</i><sub><i>L</i></sub>)—undrained shear strength (<i>S</i><sub><i>u</i></sub>) models has been enhanced through optimizing some empirical equation.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1499 - 1501"},"PeriodicalIF":5.6,"publicationDate":"2025-01-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143571013","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-23DOI: 10.1007/s11440-024-02454-1
Takayuki Sakai, Masaki Nakano
In this study, we automatically estimated the parameters of the modified Cam-Clay model, a representative constitutive model for soil. The estimation was carried out by minimizing the objective function using the dynamic multiswarm particle swarm optimization (DMS-PSO) algorithm, which is an improvement over the original PSO. The objective function was newly defined by quantifying the discrepancy between the targeted results and the model calculations in q-p′-v space. DMS-PSO divides particles into several islands to search globally and prevent local solutions, and even particles that fall into a local solution can be relocated. To evaluate the automatic estimation performance of DMS-PSO, we examined whether model parameters could be correctly estimated from the calculation results (Consideration (1)) and whether the DMS-PSO algorithm could consistently obtain the same parameter values when reproducing the experimental results (Consideration (2)). Regarding Consideration (1), the objective function was consistently smaller than 1.0 × 10–6 when the number of particles was greater than 400 and the number of islands was greater than 40. At this time, the parameter values could be estimated to the fifth decimal place. When two experiments were conducted, the estimation was obtained approximately 1.5 times faster than when only one was conducted. Regarding Consideration (2), the coefficient of variation of the parameters obtained from 100 estimations was at most 1%, and the parameter values were estimated to be approximately the same each time. In addition, narrowing the solution search range based on soil physical properties could reduce the variation in parameters by approximately 10%. Additionally, the parameters could be accurately estimated by data from at least two mechanical experiments.
{"title":"Efficient automatic estimation of soil constitutive model parameters via particle swarm optimization","authors":"Takayuki Sakai, Masaki Nakano","doi":"10.1007/s11440-024-02454-1","DOIUrl":"10.1007/s11440-024-02454-1","url":null,"abstract":"<div><p>In this study, we automatically estimated the parameters of the modified Cam-Clay model, a representative constitutive model for soil. The estimation was carried out by minimizing the objective function using the dynamic multiswarm particle swarm optimization (DMS-PSO) algorithm, which is an improvement over the original PSO. The objective function was newly defined by quantifying the discrepancy between the targeted results and the model calculations in <i>q</i>-<i>p</i>′-<i>v</i> space. DMS-PSO divides particles into several islands to search globally and prevent local solutions, and even particles that fall into a local solution can be relocated. To evaluate the automatic estimation performance of DMS-PSO, we examined whether model parameters could be correctly estimated from the calculation results (Consideration (1)) and whether the DMS-PSO algorithm could consistently obtain the same parameter values when reproducing the experimental results (Consideration (2)). Regarding Consideration (1), the objective function was consistently smaller than 1.0 × 10<sup>–6</sup> when the number of particles was greater than 400 and the number of islands was greater than 40. At this time, the parameter values could be estimated to the fifth decimal place. When two experiments were conducted, the estimation was obtained approximately 1.5 times faster than when only one was conducted. Regarding Consideration (2), the coefficient of variation of the parameters obtained from 100 estimations was at most 1%, and the parameter values were estimated to be approximately the same each time. In addition, narrowing the solution search range based on soil physical properties could reduce the variation in parameters by approximately 10%. Additionally, the parameters could be accurately estimated by data from at least two mechanical experiments.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1001 - 1017"},"PeriodicalIF":5.6,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11440-024-02454-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143570969","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}
Pub Date : 2025-01-22DOI: 10.1007/s11440-024-02518-2
Linyuan Han, Xiaoping Zhou
A coupled hydro-mechanical field-enriched finite element method (HM-FE-FEM) is developed for simulating the hydraulic fracture process of rocks subjected to in situ stresses. The staggered iterative scheme based on the Newton–Raphson iterative algorithm is utilized to address the coupled hydro-mechanical problems. The opening process of the hydraulic fracture and natural fracture subjected to in situ stresses can be reproduced by the proposed method. The accuracy of the proposed method is verified through stress intensity factor, fracture opening and experiment results. For coalescence patterns between hydraulic and natural fractures, the influences of horizontal stress difference and intersection angle are analyzed and discussed. Additionally, for coalescence patterns between hydraulic fracture and natural karst caves, the influences of cave location, cave radius and horizontal stress difference are investigated and discussed. The numerical results illustrate the effectiveness of HM-FE-FEM in handling the complex coalescence behaviors in the hydraulic fracture process of underground rock reservoirs with different defects. This study has a certain guiding significance for optimizing the fracturing operation parameters.
{"title":"A coupled hydro-mechanical field-enriched finite element method for simulating the hydraulic fracture process of rocks subjected to in situ stresses","authors":"Linyuan Han, Xiaoping Zhou","doi":"10.1007/s11440-024-02518-2","DOIUrl":"10.1007/s11440-024-02518-2","url":null,"abstract":"<div><p>A coupled hydro-mechanical field-enriched finite element method (HM-FE-FEM) is developed for simulating the hydraulic fracture process of rocks subjected to in situ stresses. The staggered iterative scheme based on the Newton–Raphson iterative algorithm is utilized to address the coupled hydro-mechanical problems. The opening process of the hydraulic fracture and natural fracture subjected to in situ stresses can be reproduced by the proposed method. The accuracy of the proposed method is verified through stress intensity factor, fracture opening and experiment results. For coalescence patterns between hydraulic and natural fractures, the influences of horizontal stress difference and intersection angle are analyzed and discussed. Additionally, for coalescence patterns between hydraulic fracture and natural karst caves, the influences of cave location, cave radius and horizontal stress difference are investigated and discussed. The numerical results illustrate the effectiveness of HM-FE-FEM in handling the complex coalescence behaviors in the hydraulic fracture process of underground rock reservoirs with different defects. This study has a certain guiding significance for optimizing the fracturing operation parameters.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1315 - 1339"},"PeriodicalIF":5.6,"publicationDate":"2025-01-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143570922","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-17DOI: 10.1007/s11440-025-02533-x
Zhiyu Gong, Guoliang Dai, Weiwei Xu, Xinsheng Chen, Hongbo Liu
The stiffened deep cement mixing (SDCM) pile is a new composite pile type, which is a stiffened core inserted into the cemented soil, compared to the DCM pile structure. The SDCM pile not only combines the advantages of DCM piles, but also improves the corresponding bearing capacity. In this paper, the bearing characteristics (the ultimate bearing capacity, the pile end resistance and the pile side friction) of PHC pipe pile and the long-core SDCM piles with number of different piles are investigated through field tests. The ultimate bearing capacities of PHC pipe pile, the long-core SDCM pile and the double long-core SDCM pile are 5.82, 8.23 and 14.67 MN. Combined with finite element analysis, it was found that the long-core SDCM pile is a friction pile in sandy soil, the cemented soil can provide lateral friction resistance and its role in bearing the vertical load is cannot be ignored. Parametric analyses were carried out for parameters: the length ratio of the cemented soil and PHC pipe pile (LC/LP), the radius ratio of the cemented soil and PHC pipe pile (RC/RP) and the modulus of elasticity of PHC pipe pile (EP) on the long-core SDCM pile with the objective of ultimate bearing capacity. The ultimate bearing capacity will be increased with increasing in LC/LP and RC/RP, and changes in EP will not affect the ultimate bearing capacity. Finally, the bearing capacity equations of the long-core SDCM pile were adjusted with the results of the study.
{"title":"Field and 3D numerical investigation on bearing characteristics of the long-core SDCM piles under vertical load in sandy soil","authors":"Zhiyu Gong, Guoliang Dai, Weiwei Xu, Xinsheng Chen, Hongbo Liu","doi":"10.1007/s11440-025-02533-x","DOIUrl":"10.1007/s11440-025-02533-x","url":null,"abstract":"<div><p>The stiffened deep cement mixing (SDCM) pile is a new composite pile type, which is a stiffened core inserted into the cemented soil, compared to the DCM pile structure. The SDCM pile not only combines the advantages of DCM piles, but also improves the corresponding bearing capacity<b>.</b> In this paper, the bearing characteristics (the ultimate bearing capacity, the pile end resistance and the pile side friction) of PHC pipe pile and the long-core SDCM piles with number of different piles are investigated through field tests. The ultimate bearing capacities of PHC pipe pile, the long-core SDCM pile and the double long-core SDCM pile are 5.82, 8.23 and 14.67 MN. Combined with finite element analysis, it was found that the long-core SDCM pile is a friction pile in sandy soil, the cemented soil can provide lateral friction resistance and its role in bearing the vertical load is cannot be ignored. Parametric analyses were carried out for parameters: the length ratio of the cemented soil and PHC pipe pile (<i>L</i><sub>C</sub>/<i>L</i><sub>P</sub>), the radius ratio of the cemented soil and PHC pipe pile (<i>R</i><sub>C</sub>/<i>R</i><sub>P</sub>) and the modulus of elasticity of PHC pipe pile (<i>E</i><sub>P</sub>) on the long-core SDCM pile with the objective of ultimate bearing capacity. The ultimate bearing capacity will be increased with increasing in <i>L</i><sub>C</sub>/<i>L</i><sub>P</sub> and <i>R</i><sub>C</sub>/<i>R</i><sub>P</sub>, and changes in <i>E</i><sub>P</sub> will not affect the ultimate bearing capacity. Finally, the bearing capacity equations of the long-core SDCM pile were adjusted with the results of the study.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1341 - 1362"},"PeriodicalIF":5.6,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143571156","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-17DOI: 10.1007/s11440-024-02440-7
Hejian Zhu, Andrew J. Whittle, Roland J.-M. Pellenq
The mechanical properties of clay minerals are largely dependent upon the chemical compositions and the mesoscale fabrics of the constituent particles. This paper describes results of a series of mesoscale molecular dynamics simulations of the hydrostatic compression and shear strain behavior for initially randomly oriented assemblies of 103 illite primary particles. The particles are simulated as rigid-body ellipsoids that interact through the single-site, Gay–Berne potential function. This corresponds to a coarse-grained model based on prior atomistic scale computation of the potential of mean force for water-mediated interactions between pairs of particles through free energy perturbation method. We investigate the mesoscale fabrics of the NPT-equilibrated assemblies for confining pressures ranging from 1.0 to 125 atm, including path dependence associated with unloading and reloading. We analyze and quantify the geometric arrangement including particle orientation, specific surface area, properties of particle stacks/aggregates, and interstack pair correlation functions. The compression of each particle assembly is associated with large irrecoverable changes in void ratio, while unloading and reloading involves much smaller, largely recoverable volumetric strains. The results are qualitatively similar to macroscopic compression behavior reported in laboratory tests. We simulate the uniaxial and shear behavior at each of the equilibrated pressure states through a series of strain-controlled steps, allowing full relaxation of the virial stresses computed at each step. The simulations investigate directional and path dependence of the shear behavior for strain deviations up to 0.2%. The results show the onset on nonlinear stiffness properties at strain levels (sim)0.01% and hysteretic behavior upon unloading and reloading. Small-strain stiffness properties of the particle assemblies are qualitatively in good agreement with quasi-static, elastic stiffness properties reported for illitic clays.
{"title":"Mesoscale simulation of the compression and small-strain elastic shear behavior of illite nanoparticle assemblies","authors":"Hejian Zhu, Andrew J. Whittle, Roland J.-M. Pellenq","doi":"10.1007/s11440-024-02440-7","DOIUrl":"10.1007/s11440-024-02440-7","url":null,"abstract":"<div><p>The mechanical properties of clay minerals are largely dependent upon the chemical compositions and the mesoscale fabrics of the constituent particles. This paper describes results of a series of mesoscale molecular dynamics simulations of the hydrostatic compression and shear strain behavior for initially randomly oriented assemblies of 10<sup>3</sup> illite primary particles. The particles are simulated as rigid-body ellipsoids that interact through the single-site, Gay–Berne potential function. This corresponds to a coarse-grained model based on prior atomistic scale computation of the potential of mean force for water-mediated interactions between pairs of particles through free energy perturbation method. We investigate the mesoscale fabrics of the NPT-equilibrated assemblies for confining pressures ranging from 1.0 to 125 atm, including path dependence associated with unloading and reloading. We analyze and quantify the geometric arrangement including particle orientation, specific surface area, properties of particle stacks/aggregates, and interstack pair correlation functions. The compression of each particle assembly is associated with large irrecoverable changes in void ratio, while unloading and reloading involves much smaller, largely recoverable volumetric strains. The results are qualitatively similar to macroscopic compression behavior reported in laboratory tests. We simulate the uniaxial and shear behavior at each of the equilibrated pressure states through a series of strain-controlled steps, allowing full relaxation of the virial stresses computed at each step. The simulations investigate directional and path dependence of the shear behavior for strain deviations up to 0.2%. The results show the onset on nonlinear stiffness properties at strain levels <span>(sim)</span>0.01% and hysteretic behavior upon unloading and reloading. Small-strain stiffness properties of the particle assemblies are qualitatively in good agreement with quasi-static, elastic stiffness properties reported for illitic clays.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 2","pages":"781 - 802"},"PeriodicalIF":5.6,"publicationDate":"2025-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143362109","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.1007/s11440-024-02511-9
Gaurav Verma, Jyoti K. Chetri, Krishna R. Reddy, Stefan J. Green
Fugitive emissions of methane (CH4), carbon dioxide (CO2), and hydrogen sulfide (H2S) from municipal solid waste landfills are major environmental concerns. To address this, a biogeochemical cover (BGCC) system is developed to mitigate these emissions. This study evaluates the effectiveness of the BGCC in comparison with a conventional soil cover (SC) system using a new large-scale laboratory setup that simulates near-field scale conditions. Both cover systems were exposed to synthetic landfill gas (LFG) across five phases, featuring varying gas compositions and influx rates. Surface emission rates and gas concentrations were continuously monitored. Post-termination of the experiments, both cover systems were dismantled, and samples were collected from different depths and locations to analyze spatial variations in physico-chemical properties. Select samples from the biocover layer of both the cover systems and basic oxygen furnace (BOF) slag layer of BGCC were subjected to batch tests to measure potential CH4 oxidation rates and residual carbonation capacity, respectively. The results showed that both cover systems achieved their highest CH4 removal efficiency at moderate influx rates (23.9–25.5 g CH4/m2-day), with BGCC's CH4 removal ranging from 74.7 to 79.7% and SC's from 83.5 to 99.8%. Complete H2S removal occurred in the biocover layer of both systems. The highest average CH4 oxidation rates were 277.9 µg CH4/g-day at 50 cm below-ground surface (bgs) in BGCC and 260.2 µg CH4/g-day at 70 cm bgs in SC, with the lowest oxidation rates observed at deeper regions (at 85 cm bgs) of both covers. The breakthrough of CO2 occurred after 156 days of continuous exposure and could be attributed to the desiccation of the BOF slag layer. Overall, the BGCC system effectively mitigated CH4, CO2, and H2S emissions, whereas the SC system only mitigated CH4 and H2S at moderate flux rates, indicating that BGCC provides a comprehensive solution for LFG mitigation.
{"title":"Large-scale testing to assess the performance of biogeochemical and soil cover systems for landfill gas mitigation","authors":"Gaurav Verma, Jyoti K. Chetri, Krishna R. Reddy, Stefan J. Green","doi":"10.1007/s11440-024-02511-9","DOIUrl":"10.1007/s11440-024-02511-9","url":null,"abstract":"<div><p>Fugitive emissions of methane (CH<sub>4</sub>), carbon dioxide (CO<sub>2</sub>), and hydrogen sulfide (H<sub>2</sub>S) from municipal solid waste landfills are major environmental concerns. To address this, a biogeochemical cover (BGCC) system is developed to mitigate these emissions. This study evaluates the effectiveness of the BGCC in comparison with a conventional soil cover (SC) system using a new large-scale laboratory setup that simulates near-field scale conditions. Both cover systems were exposed to synthetic landfill gas (LFG) across five phases, featuring varying gas compositions and influx rates. Surface emission rates and gas concentrations were continuously monitored. Post-termination of the experiments, both cover systems were dismantled, and samples were collected from different depths and locations to analyze spatial variations in physico-chemical properties. Select samples from the biocover layer of both the cover systems and basic oxygen furnace (BOF) slag layer of BGCC were subjected to batch tests to measure potential CH<sub>4</sub> oxidation rates and residual carbonation capacity, respectively. The results showed that both cover systems achieved their highest CH<sub>4</sub> removal efficiency at moderate influx rates (23.9–25.5 g CH<sub>4</sub>/m<sup>2</sup>-day), with BGCC's CH<sub>4</sub> removal ranging from 74.7 to 79.7% and SC's from 83.5 to 99.8%. Complete H<sub>2</sub>S removal occurred in the biocover layer of both systems. The highest average CH<sub>4</sub> oxidation rates were 277.9 µg CH<sub>4</sub>/g-day at 50 cm below-ground surface (bgs) in BGCC and 260.2 µg CH<sub>4</sub>/g-day at 70 cm bgs in SC, with the lowest oxidation rates observed at deeper regions (at 85 cm bgs) of both covers. The breakthrough of CO<sub>2</sub> occurred after 156 days of continuous exposure and could be attributed to the desiccation of the BOF slag layer. Overall, the BGCC system effectively mitigated CH<sub>4</sub>, CO<sub>2</sub>, and H<sub>2</sub>S emissions, whereas the SC system only mitigated CH<sub>4</sub> and H<sub>2</sub>S at moderate flux rates, indicating that BGCC provides a comprehensive solution for LFG mitigation.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1471 - 1494"},"PeriodicalIF":5.6,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143571152","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.1007/s11440-025-02537-7
Ni An, Guanqi Wang, Di Wang, Gang Ma, Xiaolin Chang, Wei Zhou
The discrete element method (DEM) is proving to be a reliable tool for studying the behavior of granular materials and has been increasingly used in recent years. The accuracy of a DEM model depends heavily on the accuracy of the particle property parameters chosen which is of vital importance for studying the mechanical properties of granular materials. However, the existing DEM parameter calibration methods are limited in terms of applicability, and the trial-and-error method remains the most common way for DEM parameter calibration. This paper presents a novel calibration method for DEM parameters using the multi-objective tree-structured parzen estimator algorithm based on prior physical information (MOTPE-PPI). The MOTPE-PPI does not rely on the training datasets and may optimize with every single test, significantly reducing the computational efforts for DEM simulation. Moreover, MOTPE-PPI is suitable for a variety of contact models and damping parameters in DEM simulation, showing robust applicability and practical feasibility. Taking an example, the DEM parameters of sandy gravel material collected from Dashixia rockfill dam in China are calibrated using MOTPE-PPI in the paper. The prior physical information is obtained through a series of triaxial loading–unloading tests, single-particle crushing tests, and literature research. Seven parameters in the rolling resistance linear contact model and breakage model are considered, and the optimization process takes only 25 iterations. Through quantitative comparison with existing parameter calibration methods, the high efficiency and wide applicability of the DEM parameter calibration method proposed in this study. The calibrated DEM parameters are used to investigate the hysteretic behavior and deformation characteristics of the granular material, revealing that the accumulation of plastic strain and resilient modulus is related to confining pressure, stress level, and the number of cycles.
{"title":"DEM parameter calibration based on multi-objective Bayesian optimization and prior physical information","authors":"Ni An, Guanqi Wang, Di Wang, Gang Ma, Xiaolin Chang, Wei Zhou","doi":"10.1007/s11440-025-02537-7","DOIUrl":"10.1007/s11440-025-02537-7","url":null,"abstract":"<div><p>The discrete element method (DEM) is proving to be a reliable tool for studying the behavior of granular materials and has been increasingly used in recent years. The accuracy of a DEM model depends heavily on the accuracy of the particle property parameters chosen which is of vital importance for studying the mechanical properties of granular materials. However, the existing DEM parameter calibration methods are limited in terms of applicability, and the trial-and-error method remains the most common way for DEM parameter calibration. This paper presents a novel calibration method for DEM parameters using the multi-objective tree-structured parzen estimator algorithm based on prior physical information (MOTPE-PPI). The MOTPE-PPI does not rely on the training datasets and may optimize with every single test, significantly reducing the computational efforts for DEM simulation. Moreover, MOTPE-PPI is suitable for a variety of contact models and damping parameters in DEM simulation, showing robust applicability and practical feasibility. Taking an example, the DEM parameters of sandy gravel material collected from Dashixia rockfill dam in China are calibrated using MOTPE-PPI in the paper. The prior physical information is obtained through a series of triaxial loading–unloading tests, single-particle crushing tests, and literature research. Seven parameters in the rolling resistance linear contact model and breakage model are considered, and the optimization process takes only 25 iterations. Through quantitative comparison with existing parameter calibration methods, the high efficiency and wide applicability of the DEM parameter calibration method proposed in this study. The calibrated DEM parameters are used to investigate the hysteretic behavior and deformation characteristics of the granular material, revealing that the accumulation of plastic strain and resilient modulus is related to confining pressure, stress level, and the number of cycles.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1379 - 1401"},"PeriodicalIF":5.6,"publicationDate":"2025-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143571151","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.1007/s11440-025-02539-5
Meng Sun, Pengfei Liu, Yuxuan Chen, Bate Bate, Fan Xue
The inclusion of calcite precipitates (CaCO3) in soft soil can improve the mechanical properties. Understanding the variability in sand stiffness due to heterogeneous precipitates is crucial for stiffness evaluation and prediction. A novel discrete element-Monte Carlo (DE-MC) method was proposed to quantify the sand stiffness variability induced by stochastic distributions of calcite precipitates, specifically focusing on shear wave velocity (Vs) as an indicator of soil stiffness. A total of 1972 samples were constructed to simulate stochastic spatial distributions of calcite precipitates. Through joint stochastic analysis, the preferential paths formed by calcite clusters were identified as significant contributors to Vs variability. The normalized connectivity per unity distance contact weight (Cd,n) exhibited the most correlated relation with Vs. Two weight selection methods were applicable for using Cd,n to characterize and predict Vs. The results suggest that the DE-MC method has the potential to assess the variability in sand stiffness quantitatively.
{"title":"Sand stiffness variability induced by stochastic distributions of calcite precipitates: a Monte Carlo-DEM study","authors":"Meng Sun, Pengfei Liu, Yuxuan Chen, Bate Bate, Fan Xue","doi":"10.1007/s11440-025-02539-5","DOIUrl":"10.1007/s11440-025-02539-5","url":null,"abstract":"<div><p>The inclusion of calcite precipitates (CaCO<sub>3</sub>) in soft soil can improve the mechanical properties. Understanding the variability in sand stiffness due to heterogeneous precipitates is crucial for stiffness evaluation and prediction. A novel discrete element-Monte Carlo (DE-MC) method was proposed to quantify the sand stiffness variability induced by stochastic distributions of calcite precipitates, specifically focusing on shear wave velocity (<i>V</i><sub>s</sub>) as an indicator of soil stiffness. A total of 1972 samples were constructed to simulate stochastic spatial distributions of calcite precipitates. Through joint stochastic analysis, the preferential paths formed by calcite clusters were identified as significant contributors to <i>V</i><sub>s</sub> variability. The normalized connectivity per unity distance contact weight (<i>C</i><sub><i>d,n</i></sub>) exhibited the most correlated relation with <i>V</i><sub>s</sub>. Two weight selection methods were applicable for using <i>C</i><sub><i>d,n</i></sub> to characterize and predict <i>V</i><sub>s</sub>. The results suggest that the DE-MC method has the potential to assess the variability in sand stiffness quantitatively.</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 3","pages":"1363 - 1377"},"PeriodicalIF":5.6,"publicationDate":"2025-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143571096","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-10DOI: 10.1007/s11440-024-02447-0
Antonia Nousiou, Georgios Anagnostou
Anhydritic claystones are widely distributed in the Gypsum Keuper formation. Their swelling is associated with the chemical process of anhydrite to gypsum transformation and has caused extensive damages in tunnels. Even though this problem has attracted great scientific interest, an adequate mathematical description of the swelling of anhydritic rocks is still missing. The present paper contributes towards closing this gap by formulating a coupled chemo-mechanical constitutive model, which considers anhydritic rock as an elastoplastic porous medium according to the principle of effective stresses, with a Mohr–Coulomb yield criterion, a non-associated flow rule and an additional, chemically induced strain component. The volumetric chemical strain is equal to the sum of the changes of the volume of the solids and of the pore volume. The change of the volume of the solids depends on the stoichiometry of the chemical reaction and is proportional to the mass of the transformed anhydrite. The pore volume may increase or decrease during the anhydrite to gypsum transformation, depending on how gypsum grows. The pore volume increases if the gypsum crystals crack and expand the matrix, and decreases if the gypsum crystals precipitate within the available pore space. The proposed model considers experimental results according to which the higher the stresses and porosity, the lower the increase in pore volume. In addition, the model assumes that the chemical strains are coaxial with the principal stresses and that the volumetric chemical strain in each principal direction is inversely proportional to the corresponding principal stress. The model is calibrated with results of tests on artificial anhydrite-kaolin specimens and achieves a very high correlation degree (R2 = 0.92).
{"title":"A chemo-mechanical model of the swelling of anhydritic claystones","authors":"Antonia Nousiou, Georgios Anagnostou","doi":"10.1007/s11440-024-02447-0","DOIUrl":"10.1007/s11440-024-02447-0","url":null,"abstract":"<div><p>Anhydritic claystones are widely distributed in the Gypsum Keuper formation. Their swelling is associated with the chemical process of anhydrite to gypsum transformation and has caused extensive damages in tunnels. Even though this problem has attracted great scientific interest, an adequate mathematical description of the swelling of anhydritic rocks is still missing. The present paper contributes towards closing this gap by formulating a coupled chemo-mechanical constitutive model, which considers anhydritic rock as an elastoplastic porous medium according to the principle of effective stresses, with a Mohr–Coulomb yield criterion, a non-associated flow rule and an additional, chemically induced strain component. The volumetric chemical strain is equal to the sum of the changes of the volume of the solids and of the pore volume. The change of the volume of the solids depends on the stoichiometry of the chemical reaction and is proportional to the mass of the transformed anhydrite. The pore volume may increase or decrease during the anhydrite to gypsum transformation, depending on how gypsum grows. The pore volume increases if the gypsum crystals crack and expand the matrix, and decreases if the gypsum crystals precipitate within the available pore space. The proposed model considers experimental results according to which the higher the stresses and porosity, the lower the increase in pore volume. In addition, the model assumes that the chemical strains are coaxial with the principal stresses and that the volumetric chemical strain in each principal direction is inversely proportional to the corresponding principal stress. The model is calibrated with results of tests on artificial anhydrite-kaolin specimens and achieves a very high correlation degree (<i>R</i><sup><i>2</i></sup> = 0.92).</p></div>","PeriodicalId":49308,"journal":{"name":"Acta Geotechnica","volume":"20 2","pages":"823 - 841"},"PeriodicalIF":5.6,"publicationDate":"2025-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11440-024-02447-0.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143361990","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}