Pub Date : 2025-09-08DOI: 10.1016/j.gete.2025.100738
Jianping Zuo , Bo Lei , Genshui Wu , Haiyan Liu , Massimo Coli
Research on the failure behavior of Longmaxi shale is vital for shale reservoir reconstruction. Shale inherently contains some initial micro-cracks, which significantly affect its strength and failure behavior. In this paper, a refined boundary multi-level parallel bonded grain-based model (multi-level PB-GBM) in Particle Flow Code (PFC2D) was developed, and the effect of inherent initial damage on shale strength and failure behavior was quantitatively investigated. The results showed that inherent initial damage significantly influences the failure pattern and mechanical properties of shale. The newly generated cracks of the initially damaged samples are significantly self-organized compared with those of the undamaged samples, indicating that the inherent initial damaged cracks induce the orientation and aggregation of micro-cracks. High initially damaged samples mainly demonstrate by splitting-shear coupled fracture as a result of the co-evolution of primary and secondary micro-cracks. Generally, rock strength gradually decreases as the initial damage increases. When the inherent initial damage within the sample is low, the rock strength is greatly influenced by confining pressure, whereas when the initial damage is high enough, the initial damage contributes more to the rock strength.
{"title":"Investigation on mechanical response and fracture behavior of initially damaged shale based on multi-level PB-GBM method","authors":"Jianping Zuo , Bo Lei , Genshui Wu , Haiyan Liu , Massimo Coli","doi":"10.1016/j.gete.2025.100738","DOIUrl":"10.1016/j.gete.2025.100738","url":null,"abstract":"<div><div>Research on the failure behavior of Longmaxi shale is vital for shale reservoir reconstruction. Shale inherently contains some initial micro-cracks, which significantly affect its strength and failure behavior. In this paper, a refined boundary multi-level parallel bonded grain-based model (multi-level PB-GBM) in Particle Flow Code (PFC2D) was developed, and the effect of inherent initial damage on shale strength and failure behavior was quantitatively investigated. The results showed that inherent initial damage significantly influences the failure pattern and mechanical properties of shale. The newly generated cracks of the initially damaged samples are significantly self-organized compared with those of the undamaged samples, indicating that the inherent initial damaged cracks induce the orientation and aggregation of micro-cracks. High initially damaged samples mainly demonstrate by splitting-shear coupled fracture as a result of the co-evolution of primary and secondary micro-cracks. Generally, rock strength gradually decreases as the initial damage increases. When the inherent initial damage within the sample is low, the rock strength is greatly influenced by confining pressure, whereas when the initial damage is high enough, the initial damage contributes more to the rock strength.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"44 ","pages":"Article 100738"},"PeriodicalIF":3.7,"publicationDate":"2025-09-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145106909","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-09-02DOI: 10.1016/j.gete.2025.100736
Jiacun Liu , Junjie Xiao , Ying Xu , Xing Li , Kaiwen Xia , Gang Han
Under the influence of high three-dimensional geostress, rocks transition into the ductile domain, undergoing continuous plastic hardening and volumetric contraction. Accurately describing the three-dimensional anisotropic deformation of rocks within ductile domain is of great significance for deep underground engineering. Therefore, a three-dimensional fractional elastoplastic constitutive within ductile domain is proposed in this study, including yield function and fractional flow rule. The ductile yield function is based on the modified Mohr-Coulomb criterion and generalized Matsuoka-Nakai deviatoric function. The deviatoric stress of yield surface is negatively correlated to hydrostatic pressure, but positively correlated to Lode angle. The yield surfaces in both meridian and deviatoric planes evolve with the plastic internal variable, accurately capturing the stress state during hardening. Two different fractional orders are used to control the plastic flow direction within meridian and deviatoric planes, represented by dilation angle and plastic deflection angle, respectively. These fractional orders are determined based on the relationship between plastic shear strain and volumetric strain, and they vary with the plastic internal variable, effectively capturing the plastic flow direction throughout hardening. The proposed model is validated using green sandstone data from hydrostatic compression and true-triaxial tests. The effect of fractional orders on the dilation angle and plastic deflection angle is discussed. Under the influence of fractional orders, both dilation angle and plastic deflection angle range from to . Besides, a comparison between the non-orthogonality and orthogonality flow rules is made. These results indicate that the fractional flow rule significantly improves the applicability and accuracy of constitutive model.
{"title":"A three-dimensional fractional elastoplastic constitutive model for rocks within ductile domain","authors":"Jiacun Liu , Junjie Xiao , Ying Xu , Xing Li , Kaiwen Xia , Gang Han","doi":"10.1016/j.gete.2025.100736","DOIUrl":"10.1016/j.gete.2025.100736","url":null,"abstract":"<div><div>Under the influence of high three-dimensional geostress, rocks transition into the ductile domain, undergoing continuous plastic hardening and volumetric contraction. Accurately describing the three-dimensional anisotropic deformation of rocks within ductile domain is of great significance for deep underground engineering. Therefore, a three-dimensional fractional elastoplastic constitutive within ductile domain is proposed in this study, including yield function and fractional flow rule. The ductile yield function is based on the modified Mohr-Coulomb criterion and generalized Matsuoka-Nakai deviatoric function. The deviatoric stress of yield surface is negatively correlated to hydrostatic pressure, but positively correlated to Lode angle. The yield surfaces in both meridian and deviatoric planes evolve with the plastic internal variable, accurately capturing the stress state during hardening. Two different fractional orders are used to control the plastic flow direction within meridian and deviatoric planes, represented by dilation angle and plastic deflection angle, respectively. These fractional orders are determined based on the relationship between plastic shear strain and volumetric strain, and they vary with the plastic internal variable, effectively capturing the plastic flow direction throughout hardening. The proposed model is validated using green sandstone data from hydrostatic compression and true-triaxial tests. The effect of fractional orders on the dilation angle and plastic deflection angle is discussed. Under the influence of fractional orders, both dilation angle and plastic deflection angle range from <span><math><msup><mrow><mn>0</mn></mrow><mo>∘</mo></msup></math></span> to <span><math><mrow><mo>−</mo><msup><mrow><mn>90</mn></mrow><mo>∘</mo></msup></mrow></math></span>. Besides, a comparison between the non-orthogonality and orthogonality flow rules is made. These results indicate that the fractional flow rule significantly improves the applicability and accuracy of constitutive model.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"44 ","pages":"Article 100736"},"PeriodicalIF":3.7,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145005361","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-09-01DOI: 10.1016/j.gete.2025.100735
Jian Xu , Liangkun Ding , Zihan Li , Jiayuan Li
This study employed the microbially induced calcium carbonate precipitation (MICP) technique to investigate the mechanism of desert sand stabilization through a multiscale approach, ranging from macro to micro levels. A multi-objective optimization model was created to enhance surface strength, CaCO3 content, and solidified layer thickness using a comprehensive analysis of multiple factors. The solidification effect was validated with wind tunnel and water retention tests. Microstructural mechanisms were examined through XRD, SEM, and PCAS. Results indicate that the optimum parameters for MICP technology are the 1:2.12 mix ratio, the 1.895 mol/L cementation solution concentration, and 4 treatment cycles. There was also a clear correlation between the performance indexes after solidification. The parameters optimized by the response surface method were essentially the same as those obtained from the experiments, with a difference of less than 5 % between the repeated test results and the optimized results. Under conditions of high CSC (single treatment cycle) or low CSC (multiple treatment cycles), MICP-treated desert sands can achieve highly efficient sand fixation and long-lasting water retention. Microanalysis revealed that increasing CSC and Tc altered the mode of particle contact from point to surface, and a significant negative correlation was observed between pore parameters and surface strength. This proves that it improves the water retention and mechanical strength of desert sand.
{"title":"MICP-enhanced wind erosion resistance of desert sand: process parameter optimization and microstructural mechanism","authors":"Jian Xu , Liangkun Ding , Zihan Li , Jiayuan Li","doi":"10.1016/j.gete.2025.100735","DOIUrl":"10.1016/j.gete.2025.100735","url":null,"abstract":"<div><div>This study employed the microbially induced calcium carbonate precipitation (MICP) technique to investigate the mechanism of desert sand stabilization through a multiscale approach, ranging from macro to micro levels. A multi-objective optimization model was created to enhance surface strength, CaCO<sub>3</sub> content, and solidified layer thickness using a comprehensive analysis of multiple factors. The solidification effect was validated with wind tunnel and water retention tests. Microstructural mechanisms were examined through XRD, SEM, and PCAS. Results indicate that the optimum parameters for MICP technology are the 1:2.12 mix ratio, the 1.895 mol/L cementation solution concentration, and 4 treatment cycles. There was also a clear correlation between the performance indexes after solidification. The parameters optimized by the response surface method were essentially the same as those obtained from the experiments, with a difference of less than 5 % between the repeated test results and the optimized results. Under conditions of high <em>CSC</em> (single treatment cycle) or low <em>CSC</em> (multiple treatment cycles), MICP-treated desert sands can achieve highly efficient sand fixation and long-lasting water retention. Microanalysis revealed that increasing <em>CSC</em> and <em>T</em><sub><em>c</em></sub> altered the mode of particle contact from point to surface, and a significant negative correlation was observed between pore parameters and surface strength. This proves that it improves the water retention and mechanical strength of desert sand.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"43 ","pages":"Article 100735"},"PeriodicalIF":3.7,"publicationDate":"2025-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144925645","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-09-01DOI: 10.1016/j.gete.2025.100706
Philip J. Vardon , Anne-Catherine Dieudonné , John. S. McCartney , Jean-Michel Pereira , David Smeulders , Guillermo Narsilio
{"title":"Accelerating the Energy Transition with Energy Geotechnics: editorial","authors":"Philip J. Vardon , Anne-Catherine Dieudonné , John. S. McCartney , Jean-Michel Pereira , David Smeulders , Guillermo Narsilio","doi":"10.1016/j.gete.2025.100706","DOIUrl":"10.1016/j.gete.2025.100706","url":null,"abstract":"","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"43 ","pages":"Article 100706"},"PeriodicalIF":3.7,"publicationDate":"2025-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145048917","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-08-29DOI: 10.1016/j.gete.2025.100731
A. Gajo
The energy and entropy balance equations of porous media saturated by one or more non-miscible fluids have been investigated by many Authors leading to expressions which are based on different thermodynamic potentials and include various simplifying assumptions. Thus the various approaches often appear to be unrelated with respect to each other. In this work, two thermodynamic potentials recently proposed in the literature for porous media saturated by one or two non-miscible and compressible pore fluids are exploited for reconsidering different and perfectly equivalent expressions of the energy balance equations given in terms of internal energies, entropies and enthalpies, without simplifying assumptions. In particular, the entropy fluxes and the dissipation functions are presented for a simple case of irreversible response of the solid skeleton, involving neither irreversibility of the solid grain response, nor elastoplastic coupling nor frozen inelastic energy. Some comparisons with the formulations proposed in the literature are discussed.
{"title":"Energy and entropy balance laws for porous media saturated by one or two non-miscible pore fluids at different temperatures","authors":"A. Gajo","doi":"10.1016/j.gete.2025.100731","DOIUrl":"10.1016/j.gete.2025.100731","url":null,"abstract":"<div><div>The energy and entropy balance equations of porous media saturated by one or more non-miscible fluids have been investigated by many Authors leading to expressions which are based on different thermodynamic potentials and include various simplifying assumptions. Thus the various approaches often appear to be unrelated with respect to each other. In this work, two thermodynamic potentials recently proposed in the literature for porous media saturated by one or two non-miscible and compressible pore fluids are exploited for reconsidering different and perfectly equivalent expressions of the energy balance equations given in terms of internal energies, entropies and enthalpies, without simplifying assumptions. In particular, the entropy fluxes and the dissipation functions are presented for a simple case of irreversible response of the solid skeleton, involving neither irreversibility of the solid grain response, nor elastoplastic coupling nor frozen inelastic energy. Some comparisons with the formulations proposed in the literature are discussed.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"44 ","pages":"Article 100731"},"PeriodicalIF":3.7,"publicationDate":"2025-08-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145333440","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-08-24DOI: 10.1016/j.gete.2025.100734
Yan-jie Feng , Cheng-zhi Qi , Fa Zhao , Tao Li , An-sen Gao , Xiao-yu Ma
The deformation and failure of rock are affected by the strain rate. However, the effect of the strain rate on rock failure has not been well studied at the microscale level. This study focuses on this complexity by constructing a three-scale wing-shaped crack propagation model that accounts for the interactions between cracks. Based on this model, we determined the initiation and coalescence times of 3 scale levels cracks to reveal the effect of the strain rate on the propagation and failure patterns of multiscale cracks. In addition, we analyzed the critical strain rate (strain rate required for simultaneous coalescence of adjacent scale-level cracks) that leads to failure. The results show that both crack initiation and coalescence times decrease significantly with increasing strain rate and that an increase in initial crack length leads to earlier crack initiation. For the multiscale crack model, as the strain rate increases, large-scale cracks (1-st-scale level cracks) coalesce first, followed sequentially by medium- and small-scale cracks (2-nd and 3-d-scale level cracks). Furthermore, we observed that the critical strain rate initially increases and then decreases with increasing initial crack concentration. Moreover, both an increase in the initial crack length and size decrease factor (ratio of the length of a specific scale-level crack to the length of adjacent larger scale-level cracks) lead to a decrease in the critical strain rate, further confirming the influence of crack size on the failure properties of rock. A comparison with the existing theoretical model shows that the proposed theoretical model is reasonable.
{"title":"Influence of the strain rate on the deformation and failure of rocks with multiscale cracks","authors":"Yan-jie Feng , Cheng-zhi Qi , Fa Zhao , Tao Li , An-sen Gao , Xiao-yu Ma","doi":"10.1016/j.gete.2025.100734","DOIUrl":"10.1016/j.gete.2025.100734","url":null,"abstract":"<div><div>The deformation and failure of rock are affected by the strain rate. However, the effect of the strain rate on rock failure has not been well studied at the microscale level. This study focuses on this complexity by constructing a three-scale wing-shaped crack propagation model that accounts for the interactions between cracks. Based on this model, we determined the initiation and coalescence times of 3 scale levels cracks to reveal the effect of the strain rate on the propagation and failure patterns of multiscale cracks. In addition, we analyzed the critical strain rate (strain rate required for simultaneous coalescence of adjacent scale-level cracks) that leads to failure. The results show that both crack initiation and coalescence times decrease significantly with increasing strain rate and that an increase in initial crack length leads to earlier crack initiation. For the multiscale crack model, as the strain rate increases, large-scale cracks (1-st-scale level cracks) coalesce first, followed sequentially by medium- and small-scale cracks (2-nd and 3-d-scale level cracks). Furthermore, we observed that the critical strain rate initially increases and then decreases with increasing initial crack concentration. Moreover, both an increase in the initial crack length and size decrease factor (ratio of the length of a specific scale-level crack to the length of adjacent larger scale-level cracks) lead to a decrease in the critical strain rate, further confirming the influence of crack size on the failure properties of rock. A comparison with the existing theoretical model shows that the proposed theoretical model is reasonable.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"43 ","pages":"Article 100734"},"PeriodicalIF":3.7,"publicationDate":"2025-08-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144894788","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This paper presents a numerical comparison of the vertical pull-out capacity of square and circular anchors in purely cohesive soils (i.e. clays in undrained conditions). For simplicity, ultrathin, infinitely rigid anchors are considered and to isolate the effect of anchor shape, comparisons are made between anchors of equal area and embedment depth. Finite Element Limit Analyses (FELA) are used to compute upper and lower bound values of the break-out factor over the full range of embedment ratios, and the associated failure mechanisms are identified. The results show for the first time (to the best of the authors’ knowledge) that square anchors exhibit slightly higher efficiency at shallow embedment ratios due to their larger perimeter, while at greater depths, circular anchors become more efficient as a result of the different failure mechanisms involved. The study also investigates the influence of anchor inclination and shows that inclined anchors have a higher pull-out capacity in vented conditions due to elongated failure mechanisms. Under attached conditions, the deep failure mechanism is obtained in most cases with the corresponding constant break-out factor. In addition, the paper analyses the influence of anchor spacing in anchor groups, identifying optimal spacing to avoid capacity reduction due to interaction effects. For shallow depths, a spacing of about two times the anchor width is sufficient, while deeper installations require larger spacings due to the extended failure zone. Once the deep failure mechanism is reached, spacing requirements decrease again, less than two times the anchor width. Overall, the presented numerical simulations offer insights for the design of plate anchors in cohesive soils, contributing to the advancement of offshore foundation technologies.
{"title":"Numerical comparison between square and circular plate anchors in clay","authors":"Mohammadreza Jahanshahinowkandeh, Marina Miranda, Jorge Castro","doi":"10.1016/j.gete.2025.100733","DOIUrl":"10.1016/j.gete.2025.100733","url":null,"abstract":"<div><div>This paper presents a numerical comparison of the vertical pull-out capacity of square and circular anchors in purely cohesive soils (i.e. clays in undrained conditions). For simplicity, ultrathin, infinitely rigid anchors are considered and to isolate the effect of anchor shape, comparisons are made between anchors of equal area and embedment depth. Finite Element Limit Analyses (FELA) are used to compute upper and lower bound values of the break-out factor over the full range of embedment ratios, and the associated failure mechanisms are identified. The results show for the first time (to the best of the authors’ knowledge) that square anchors exhibit slightly higher efficiency at shallow embedment ratios due to their larger perimeter, while at greater depths, circular anchors become more efficient as a result of the different failure mechanisms involved. The study also investigates the influence of anchor inclination and shows that inclined anchors have a higher pull-out capacity in vented conditions due to elongated failure mechanisms. Under attached conditions, the deep failure mechanism is obtained in most cases with the corresponding constant break-out factor. In addition, the paper analyses the influence of anchor spacing in anchor groups, identifying optimal spacing to avoid capacity reduction due to interaction effects. For shallow depths, a spacing of about two times the anchor width is sufficient, while deeper installations require larger spacings due to the extended failure zone. Once the deep failure mechanism is reached, spacing requirements decrease again, less than two times the anchor width. Overall, the presented numerical simulations offer insights for the design of plate anchors in cohesive soils, contributing to the advancement of offshore foundation technologies.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"43 ","pages":"Article 100733"},"PeriodicalIF":3.7,"publicationDate":"2025-08-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144911968","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-08-20DOI: 10.1016/j.gete.2025.100732
Satar Mahdevari , Pedram Bakhtiari Haftlang
Coal remains a cornerstone of global energy supply, driving the need for more efficient and technologically advanced extraction methods. This study introduces a numerical framework that couples the Smoothed Particle Hydrodynamics (SPH) with the Finite Element Method (FEM) to model the dynamic response of coal under waterjet-assisted cutting—an emerging technique recognized for its applicability, minimal stress disturbance, and safe working conditions in underground mining. Implemented in LS-DYNA, the model captures two-phase fluid–solid interactions, including jet-induced fracture initiation, propagation, and material removal. A detailed parametric investigation evaluates the effects of jet velocity, nozzle diameter, impingement angle, and cutting duration on coal fragmentation behavior. Model predictions were rigorously validated through controlled laboratory experiments, achieving reliable correlation with empirical results—showing mean absolute errors of 7.2 % in Cutting Depth (CD) and 5.8 % in Cutting Volume (CV). To address the performance constraints of Pure Water Jet (PWJ) systems, extended simulations were conducted for Abrasive Water Jet (AWJ) and Ice Abrasive Water Jet (IAWJ) techniques. The AWJ configuration enhanced CD and CV by 51 % and 66 %, respectively, while IAWJ achieved up to 20 % improvement over PWJ. Stress field analysis further revealed that increased jet velocity is significantly more effective than nozzle enlargement in maximizing cutting efficiency. These findings not only validate the SPH–FEM model as a predictive tool but also offer actionable insights for optimizing next-generation waterjet systems in deep coal mining applications.
{"title":"Coupled SPH–FEM modeling of waterjet-assisted coal cutting: Numerical simulation and experimental validation","authors":"Satar Mahdevari , Pedram Bakhtiari Haftlang","doi":"10.1016/j.gete.2025.100732","DOIUrl":"10.1016/j.gete.2025.100732","url":null,"abstract":"<div><div>Coal remains a cornerstone of global energy supply, driving the need for more efficient and technologically advanced extraction methods. This study introduces a numerical framework that couples the Smoothed Particle Hydrodynamics (SPH) with the Finite Element Method (FEM) to model the dynamic response of coal under waterjet-assisted cutting—an emerging technique recognized for its applicability, minimal stress disturbance, and safe working conditions in underground mining. Implemented in LS-DYNA, the model captures two-phase fluid–solid interactions, including jet-induced fracture initiation, propagation, and material removal. A detailed parametric investigation evaluates the effects of jet velocity, nozzle diameter, impingement angle, and cutting duration on coal fragmentation behavior. Model predictions were rigorously validated through controlled laboratory experiments, achieving reliable correlation with empirical results—showing mean absolute errors of 7.2 % in Cutting Depth (CD) and 5.8 % in Cutting Volume (CV). To address the performance constraints of Pure Water Jet (PWJ) systems, extended simulations were conducted for Abrasive Water Jet (AWJ) and Ice Abrasive Water Jet (IAWJ) techniques. The AWJ configuration enhanced CD and CV by 51 % and 66 %, respectively, while IAWJ achieved up to 20 % improvement over PWJ. Stress field analysis further revealed that increased jet velocity is significantly more effective than nozzle enlargement in maximizing cutting efficiency. These findings not only validate the SPH–FEM model as a predictive tool but also offer actionable insights for optimizing next-generation waterjet systems in deep coal mining applications.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"43 ","pages":"Article 100732"},"PeriodicalIF":3.7,"publicationDate":"2025-08-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144886021","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-08-05DOI: 10.1016/j.gete.2025.100730
Lei Qin , Jiawei Li , Haifei Lin , Shugang Li , Miao Mu , Shiyin Lv , Niandong Chen
Geological variability results in coal seams with diverse ranks exhibiting distinct physical properties, critically influencing liquid nitrogen (LN₂) fracturing outcomes. We systematically assess the mechanical and acoustic damage induced by single LN2 freeze-thaw (LNSFT) and repeated freeze-thaw cycles (LNCFT) in three representative coal ranks—lignite, bituminite, and anthracite—via ultrasonic measurement, uniaxial compression, and acoustic emission (AE) techniques. Results demonstrate that initially, pore water solidification enhances coal strength and acoustic integrity; subsequently, crack initiation and propagation induced by frost heave, thermal stress, and LN₂ expansion progressively weaken these properties. This balance between strengthening and weakening is primarily governed by coal pore structure, fissures, and moisture content. AE patterns under loading distinctly follow steady-state, activation, and attenuation phases, with both the freezing and thawing phases promoting shear-oriented fracture development. Damage indices (D), computed from ultrasonic P-wave velocity (v), peak strength (σ), and elastic modulus (E), reveal an inverse correlation between freeze-thaw damage severity and coal rank, indicating that higher-rank coals exhibit greater structural stability and freeze-thaw resistance. Furthermore, under equivalent cumulative freezing durations, LNCFT cause significantly greater damage than LNSFT, highlighting a cumulative damage effect. These insights provide critical guidance for optimizing LN₂ fracturing techniques aimed at enhancing coal seam permeability.
{"title":"Damage and acoustic characteristics of water-saturated coals with different ranks under liquid nitrogen freezing and thawing treatments","authors":"Lei Qin , Jiawei Li , Haifei Lin , Shugang Li , Miao Mu , Shiyin Lv , Niandong Chen","doi":"10.1016/j.gete.2025.100730","DOIUrl":"10.1016/j.gete.2025.100730","url":null,"abstract":"<div><div>Geological variability results in coal seams with diverse ranks exhibiting distinct physical properties, critically influencing liquid nitrogen (LN₂) fracturing outcomes. We systematically assess the mechanical and acoustic damage induced by single LN<sub>2</sub> freeze-thaw (LNSFT) and repeated freeze-thaw cycles (LNCFT) in three representative coal ranks—lignite, bituminite, and anthracite—via ultrasonic measurement, uniaxial compression, and acoustic emission (AE) techniques. Results demonstrate that initially, pore water solidification enhances coal strength and acoustic integrity; subsequently, crack initiation and propagation induced by frost heave, thermal stress, and LN₂ expansion progressively weaken these properties. This balance between strengthening and weakening is primarily governed by coal pore structure, fissures, and moisture content. AE patterns under loading distinctly follow steady-state, activation, and attenuation phases, with both the freezing and thawing phases promoting shear-oriented fracture development. Damage indices (<em>D</em>), computed from ultrasonic P-wave velocity (<em>v</em>), peak strength (<em>σ</em>), and elastic modulus (<em>E</em>), reveal an inverse correlation between freeze-thaw damage severity and coal rank, indicating that higher-rank coals exhibit greater structural stability and freeze-thaw resistance. Furthermore, under equivalent cumulative freezing durations, LNCFT cause significantly greater damage than LNSFT, highlighting a cumulative damage effect. These insights provide critical guidance for optimizing LN₂ fracturing techniques aimed at enhancing coal seam permeability.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"43 ","pages":"Article 100730"},"PeriodicalIF":3.7,"publicationDate":"2025-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144781379","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Energy piles present an innovative energy-saving technology that can fulfill two critical building needs of structural support and energy supply. In practice, developing simple but efficient methods to predict the thermomechanical response of energy pile foundations is essential for geotechnical engineers. In this study, a practical method was proposed for the thermomechanical analyses of energy pile foundations. The proposed method could effectively describe the interactions between the grouped energy piles, the surrounding soil, and the stiff soil strata underlying the pile tip. Based on this method, parametric analyses were performed to evaluate the effects of several aspects, including the foundation geometries and ground properties, on the pile displacement behavior and the pile-to-pile interaction. Further, the proposed method was used for the displacement analysis for a square pile group containing sixteen energy piles under thermomechanical loads. Comparisons with results obtained through the experimental investigations and finite-element methods prove that the proposed method is capable of capturing the displacement response of energy pile foundations with reasonable accuracy. The aim of this study is to offer a practical method and a reliable reference to geotechnical engineers during the design of energy pile foundations.
{"title":"Displacement analysis for energy pile foundations under thermomechanical loads","authors":"Jincheng Fang , Shijin Feng , Yong Zhao , Hongxin Chen","doi":"10.1016/j.gete.2025.100725","DOIUrl":"10.1016/j.gete.2025.100725","url":null,"abstract":"<div><div>Energy piles present an innovative energy-saving technology that can fulfill two critical building needs of structural support and energy supply. In practice, developing simple but efficient methods to predict the thermomechanical response of energy pile foundations is essential for geotechnical engineers. In this study, a practical method was proposed for the thermomechanical analyses of energy pile foundations. The proposed method could effectively describe the interactions between the grouped energy piles, the surrounding soil, and the stiff soil strata underlying the pile tip. Based on this method, parametric analyses were performed to evaluate the effects of several aspects, including the foundation geometries and ground properties, on the pile displacement behavior and the pile-to-pile interaction. Further, the proposed method was used for the displacement analysis for a square pile group containing sixteen energy piles under thermomechanical loads. Comparisons with results obtained through the experimental investigations and finite-element methods prove that the proposed method is capable of capturing the displacement response of energy pile foundations with reasonable accuracy. The aim of this study is to offer a practical method and a reliable reference to geotechnical engineers during the design of energy pile foundations.</div></div>","PeriodicalId":56008,"journal":{"name":"Geomechanics for Energy and the Environment","volume":"43 ","pages":"Article 100725"},"PeriodicalIF":3.7,"publicationDate":"2025-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144826998","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}