{"title":"A Riemann-based two-phase two-layer SPH method for simulating submarine landslide tsunamis","authors":"Yue Fang , Qiang Xu , Jianyun Chen","doi":"10.1016/j.compgeo.2025.107247","DOIUrl":"10.1016/j.compgeo.2025.107247","url":null,"abstract":"<div><div>Capturing dynamic water–soil variations and addressing water/soil boundary in numerical simulations of submarine landslide tsunamis remain major challenges. This study proposes a new two-phase two-layer Riemann-smoothed particle hydrodynamics (SPH) numerical model. In this model, the soil satisfies the Drucker–Prager yield criterion and water is treated as a weakly-compressible fluid. The SPH method based on a low-dissipation Riemann solver, is innovatively introduced into the mathematical framework of water–soil coupling to simulate the dynamic behaviors of two overlapping water and soil particle layers. The no-slip boundary treatment based on a one-sided Riemann format is proposed to achieve the water/soil–solid coupling in landslide simulations. The volume fraction is fully participated in the current SPH discretization process, and the Riemann solver is also applied to improve its calculation. Validations against two static and dynamic water–soil coupling cases demonstrate that present approach is effective. Then, the proposed model with different particle resolutions is further applied to explore the submarine landslide tsunamis under two different configurations. The landslide motions exhibit good agreement and convergence with experimental data, and dynamic water–soil evolutions during the landslide process are properly captured with reasonable dissipation control, which indicates the accuracy and usability of the current Riemann-based two-layer model.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107247"},"PeriodicalIF":5.3,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143791617","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}
Hao Liu , Yuan Gao , Jiaxuan Li , Xiaoyan Yang , Xin Liu , M.Hehsam El Naggar , Guoxiong Mei , Wenbing Wu
{"title":"A meta-damping model for soil-pile interaction in large-diameter pipe piles under horizontal dynamic loading","authors":"Hao Liu , Yuan Gao , Jiaxuan Li , Xiaoyan Yang , Xin Liu , M.Hehsam El Naggar , Guoxiong Mei , Wenbing Wu","doi":"10.1016/j.compgeo.2025.107233","DOIUrl":"10.1016/j.compgeo.2025.107233","url":null,"abstract":"<div><div>A comprehensive understanding of soil-pile interaction (SPI) is crucial for applications in earthquake engineering and offshore structures subjected to dynamic loads. The soil plug within large diameter pipe piles plays a significant role in dynamic horizontal loading. Previous models primarily focus on the soil plug’s resistance, neglecting its inertial effects, which leads to discrepancies in predicted natural frequencies and stiffness. This study introduces a <em>meta</em>-damping model to simultaneously account for both the resistance and inertial effects of the soil plug. A theoretical model for pipe piles embedded in saturated soil under horizontal dynamic loading is developed using this model. The model’s validity and accuracy are verified through comparisons with experimental data and existing theoretical models. The impact of the soil plug on the dynamic characteristics of the pile is explored, revealing that traditional plane strain models, which neglect the inertial effects, tend to overestimate the natural frequency and underestimate dynamic damping. The discrepancies are more pronounced for piles with larger diameters and poorer soil properties. This highlights the importance of considering the soil plug’s inertial effects for more accurate dynamic analysis.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107233"},"PeriodicalIF":5.3,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143791616","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":"Dissipative effects of baffles on the dynamics of debris flow and its impact on downstream structure","authors":"Trung-Kien Nguyen , Nhu H.T. Nguyen , Thanh-Trung Vo , Liuxin Chen","doi":"10.1016/j.compgeo.2025.107246","DOIUrl":"10.1016/j.compgeo.2025.107246","url":null,"abstract":"<div><div>Debris flows triggered by heavy rainfall constitute a mixture of slurry and earthen aggregates. These flows travel downhill at a high speed, causing significant damage to downstream facilities. Baffles are among the common protective structures installed along debris flow paths to decelerate the flows, impeding their impact on downstream structures. Previous research mainly focused on the interactions between baffle systems and debris flows, while the flow impact responses, especially on downstream structures, were rarely examined. Based on a coupled VOF-DEM method, this study investigates the impeding effects of baffles on the dynamics of debris mixture flows and their impact on a downstream barrier. We found that baffles are more effective in decelerating the flows with low Froude numbers generally occurring on gentle slopes. However, debris flow with high Froude numbers run up and overflow through baffle arrays, reducing interactions between baffles and the flows. Increasing the number of baffle rows in a baffle system can further decelerate the flows, thereby transitioning flow run-up to stable pipe-up at the downstream barrier, reducing dynamic impact force on the barrier. Remarkably, regardless of flow inertial conditions, baffle systems with more baffle rows and shorter baffles alleviate the flow impact on the downstream barrier better than those with fewer rows and higher baffles. This finding provides a reference for baffle designers to prioritise baffle systems with more rows rather than higher baffles.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107246"},"PeriodicalIF":5.3,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143785717","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}
Bing Bai , Bixia Zhang , Yanjie Ji , Yongchen Zong
{"title":"A thermodynamic multi-field model for unsaturated sulfate-saline soils considering crystallization process","authors":"Bing Bai , Bixia Zhang , Yanjie Ji , Yongchen Zong","doi":"10.1016/j.compgeo.2025.107251","DOIUrl":"10.1016/j.compgeo.2025.107251","url":null,"abstract":"<div><div>A novel thermo-hydro-mechanical-chemical (THMC) coupling model grounded in thermodynamic dissipation theory was established to unravel the intricate behavior of unsaturated sulfate-saline soils during cooling crystallization. The model quantifies energy transfer and dissipation during crystallization and introduces a method to calculate the amount of sulfate crystallization. It intricately captures the interdependencies between crystallization, pore water pressure, crystallization pressure and volumetric expansion, while also accounting for the dynamic feedback of latent heat from phase transitions on heat conduction. The reliability of the model was validated through experimental data. Numerical simulations explored the effects of cooling paths, thermal conductivity, initial salt content and initial porosity on the crystallization behavior and mechanical properties. The model provides theoretical support for optimizing the engineering design and facility maintenance of sulfate-saline soils.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107251"},"PeriodicalIF":5.3,"publicationDate":"2025-04-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143785716","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}
Xiaolin Huang , Weiqi Kang , Shengwen Qi , Xiaohui Zhang , Jiahu Du
{"title":"Comparative Analysis of microcracking behaviors and associated acoustic emission characteristics in sandstone subjected to compression-induced and direct tensile stresses","authors":"Xiaolin Huang , Weiqi Kang , Shengwen Qi , Xiaohui Zhang , Jiahu Du","doi":"10.1016/j.compgeo.2025.107236","DOIUrl":"10.1016/j.compgeo.2025.107236","url":null,"abstract":"<div><div>Rocks often undergo both direct tensile actions and indirect tensile actions induced by compression due to tectonic movements and human engineering activities, leading to the gradual accumulation of microcracks accompanied by acoustic emission (AE) events. However, to date, the differences in the microcracking behavior and associated AE of rocks under these two tensile conditions, as well as the underlying micromechanical mechanisms, have not been thoroughly understood. Here, we conducted a direct tensile test on a dog-bone-shaped sandstone sample and a Brazilian splitting test (representing compression-induced tension) on a disc-shaped sample. Corresponding nonlinear particle-based discrete element models were developed to simulate the tensile stress-displacement responses and macroscopic failure modes of the two sandstone samples. Results showed that sandstone under direct tension exhibits significantly lower strength than under compression-induced tension. In the disc-shaped model, microcracks appeared progressively at dispersed times, accompanied by high-frequency AE signals, mainly concentrated near the longitudinal symmetric axis, adjacent to the loading point. Conversely, microcracks in the dog-bone-shaped model appeared abruptly and were concentrated within a narrow range near the horizontal symmetric axis, with microcracks and energy release uniformly distributed, and fractures tending to be perpendicular to the loading direction. Compared to compression-induced tension, the dog-bone-shaped model exhibited smaller AE counts, energy release, and AE signal frequency, but a larger <em>b</em>-value, more than three times as high. Micromechanical analysis revealed that these differences are primarily due to significant variations in tensile/compressive contact force chains and strain energy distributions in particle assemblies under direct and compression-induced tension.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"184 ","pages":"Article 107236"},"PeriodicalIF":5.3,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143767856","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 non-orthogonal elastoplastic constitutive model incorporating cohesion degradation for frozen soil","authors":"Jingyu Liang , Haowen Li , Jilin Qi , Dechun Lu","doi":"10.1016/j.compgeo.2025.107231","DOIUrl":"10.1016/j.compgeo.2025.107231","url":null,"abstract":"<div><div>Cohesion provided by pore ice is a critical component influencing the mechanical behavior of frozen soil, as it not only cements soil particles together but also shares the external loads with them. In view the crucial role of cohesion in developing an elastoplastic model for frozen soil, this paper employs triaxial tensile strength (TTS) to characterize cohesion and proposes a TTS degradation expression driven by plastic shear strain. By directly incorporating TTS into the yield function, a framework for a Non-Orthogonal Elastoplastic (NOEP) constitutive model that accounts for cohesion degradation in frozen soil is developed. Furthermore, a hardening parameter incorporating TTS is introduced and used in conjunction with the modified yield function to determine the magnitude of the plastic strain increment. The non-orthogonal plastic flow rule is used to determine the direction of the plastic strain increment based on the modified yield function. Ultimately, by combining the elastic strain increment determined by Hooke’s rule, a NOEP constitutive model incorporating cohesion degradation for frozen soil is established. The validity and rationality of the proposed NOEP model in representing the stress–strain relationship of frozen soil are confirmed through comparisons with test results of frozen soil under the triaxial compression conditions. The proposed constitutive model provides a more comprehensive and precise representation of frozen soil’s response to external loading, enhancing the understanding of its shear deformation behavior and providing a robust theoretical foundation for engineering design and construction in cold regions.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107231"},"PeriodicalIF":5.3,"publicationDate":"2025-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143769105","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":"Nanoscale contact behavior of α-quartz asperities — A molecular dynamics approach","authors":"Sheng Li , Eiichi Fukuyama","doi":"10.1016/j.compgeo.2025.107239","DOIUrl":"10.1016/j.compgeo.2025.107239","url":null,"abstract":"<div><div>Rock surfaces are composed of asperities at various scales. Accurate modeling of the multiscale contact deformation of asperities is of great importance to study the earthquake source mechanics and hydraulic properties of fractures. In this study, to shed light on the nanoscale contact behavior of rock asperities that might not have been paid much attention to in previous studies, we simulated a series of nanoscale contact processes of <em>α</em>-quartz asperities under single-asperity contact employing molecular dynamics method. In addition to the investigation of the influence of contact configuration and asperity size on the asperity failure, we compared the simulation results with macroscopically elastic and elastoplastic contact models to evaluate the multiscale applicability. We observed that fracture occurs in the nanoscale <em>α</em>-quartz asperities during contact process, inconsistent with the traditional assumptions of elastic and elastoplastic deformation of rock asperities. Moreover, the asperity size and contact configuration significantly affect the failure mechanism of <em>α</em>-quartz asperities, specifically the transition from plasticity-dominant failure to fracture damage. Finally, a multiscale disparity exists, where macroscopic contact models are limited to predict the nanoscale contact behavior of <em>α</em>-quartz asperities. The results obtained here underscore the need to re-evaluate the assumptions of rock asperity deformation and emphasize the importance of considering contact configurations and multiscale effects of rock asperities.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107239"},"PeriodicalIF":5.3,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143760255","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}
Chih-Wei Lu , Minh-Tam Doan , Shi-Shuenn Chen , Jun-Yang Shi
{"title":"Probabilistic fragility assessment of shallow-founded buildings in liquefiable soils treated by enhanced drainage and densification considering short and long-term responses","authors":"Chih-Wei Lu , Minh-Tam Doan , Shi-Shuenn Chen , Jun-Yang Shi","doi":"10.1016/j.compgeo.2025.107213","DOIUrl":"10.1016/j.compgeo.2025.107213","url":null,"abstract":"<div><div>This study quantifies the seismic fragility assessment of shallow-founded buildings in liquefiable and treated soils, enhanced by drainage and densification, considering both short- and long-term behaviors. A conceptual framework is proposed for developing seismic fragility curves based on engineering demand parameters (EDPs) of buildings subjected to various earthquake magnitudes. The framework for establishing seismic fragility curves involves three essential steps. First, nonlinear dynamic analyses of soil-building systems are performed to assess both the short-term response, which occurs immediately following an earthquake, and the long-term response, when excess pore water pressure completely dissipates, and generate a dataset of building settlements. The seismic responses are compared in terms of excess pore water pressure buildup, immediate and residual ground deformation, and building settlement to explore the dynamic mechanisms of soil-building systems and evaluate the performance of enhanced drainage and densification over short- and long-term periods. Second, 38 commonly used and newly proposed intensity measures (IMs) of ground motions (GMs) are comprehensively evaluated using five statistical measures, such as correlation, efficiency, practicality, proficiency, and sufficiency, to identify optimal IMs of GMs. Third, fragility curves are developed to quantify probability of exceeding various capacity limit states, based on structural damage observed in Taiwan, for both liquefaction-induced immediate and residual settlements of buildings under different levels of IMs. Overall, this study proposes a rapid and straightforward probabilistic assessment approach for buildings in liquefiable soils, along with remedial countermeasures to enhance seismic resilience.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107213"},"PeriodicalIF":5.3,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143760254","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":"DEM multi-scale insights on the stress, energy and crack propagation in proppant-fractured rock collisions","authors":"Sarath C.R. Nallala , Huan He , Kostas Senetakis","doi":"10.1016/j.compgeo.2025.107234","DOIUrl":"10.1016/j.compgeo.2025.107234","url":null,"abstract":"<div><div>The collision behavior of an impactor (proppant simulant) colliding on weak clastic rock with pre-existing microcracks is examined in this work at multi-scales. Using discrete element-based analysis, the linear parallel bond contact model is used to simulate the rock base, while the Hertz-Mindlin contact model is used to capture the non-linear force–displacement relationship between grain-block contacts. A systematic parametric study is conducted to assess the impact of microcrack size, number, angle, and position on the damage potential of the impactor colliding the base rock. The results indicate that some of these parameters significantly influence the damage potential while others have minimal effects. Multiscale assessment is conducted by comparing stress and energy distribution results at the macroscale, penetration and contact formation information at the mesoscale and studying the discrete fracture network at the microscale. The findings of this study provide important insights into the collision behavior of proppants with clastic rocks, particularly during hydraulic fracturing operations and can potentially aid in the modeling of these systems at larger scales.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107234"},"PeriodicalIF":5.3,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143760256","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":"Interpreting random field parameters of spatially variable soil using ball penetrometer","authors":"Zhixuan Li, Ping Yi, Jun Liu","doi":"10.1016/j.compgeo.2025.107219","DOIUrl":"10.1016/j.compgeo.2025.107219","url":null,"abstract":"<div><div>Natural soil exhibits inherent spatial variability, which is crucial for site characterization and the reliability design of foundations. However, the interpretation accuracy of random field parameters for spatially variable soil, based on the undrained shear strength (<em>s</em><sub>u</sub>) profile obtained by the penetrometer, is hampered by the spatial averaging effect and weakest path effect. These effects become more pronounced as the observation scale (i.e., the ball diameter in this study) and the coefficient of variation of soil strength increase, respectively. Interpreting the spatial variability of actual soil strength based on <em>s</em><sub>u</sub> profile obtained by penetrometer remains challenging. In the present study, the random field parameters of the <em>s</em><sub>u</sub> profile obtained by the ball penetrometer are estimated using large deformation finite element analysis (LDFEA) and Markov chain Monte Carlo method. A filtering procedure, based on fast Fourier transform, is proposed to filter out the numerical oscillation of LDFEA on the <em>s</em><sub>u</sub> profile. Then the effect of the observation scale on the estimated random field parameters is investigated and it is proved that the estimated random field parameters are significantly influenced by the observation scale. The smaller the penetrometer dimension is, the closer the estimated random field parameters of the <em>s</em><sub>u</sub> profile approach those of actual soil strength. Therefore, a correction procedure is proposed and reliability analysis of a caisson and a circular foundation show that directly employing the estimated random field parameters in foundation design may lead to unsafe outcomes. Therefore, the estimated random field parameters should be corrected by the proposed interpretation procedure.</div></div>","PeriodicalId":55217,"journal":{"name":"Computers and Geotechnics","volume":"183 ","pages":"Article 107219"},"PeriodicalIF":5.3,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143760251","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}