Computational Particle Mechanics最新文献

{"title":"Study on coal wall spalling characteristics and stability control of steeply inclined coal seam mining face","authors":"Shengli Yang,&nbsp;Qiang Li,&nbsp;Hao Yue,&nbsp;Shuai Yang,&nbsp;Fengqi Liu","doi":"10.1007/s40571-024-00897-7","DOIUrl":"10.1007/s40571-024-00897-7","url":null,"abstract":"<div><p>Coal wall spalling characteristics and stability control mechanisms in steeply inclined coal seam (SICS) mining are more complex and dynamic due to the imbalanced mechanical constraint environment of the coal wall. This study investigates the characteristics of coal wall spalling in SICS mining using field measurements, theoretical analysis, and numerical simulations. The internal relationship between coal wall spalling and different influencing factors is also established, and the corresponding coal wall stability control measures are put forward. Results indicate that the primary form of coal wall spalling is wedge-shaped. Under the combined action of roof pressure and gravity, the wedge body produces two components along the dip and strike directions, forming two structural planes, a and b. The a-plane predominantly experiences shear slip failure, while the b-plane primarily undergoes tensile failure. The failure criteria for the a-plane and b-plane of the wedge body are determined, with mining height, coal strength, roof pressure, and coal seam dip angle being the main influencing factors in SICS mining. The stability relationship of the coal wall under various advancing distances is investigated using the 3DEC numerical simulation. The majority of coal wall spalling occurs as an asymmetric wedge body, and the relationship between different influencing factors and coal wall spalling is examined. A comprehensive preventive and control technique for coal wall spalling in SICS mining is proposed. The engineering application demonstrates positive results and provides theoretical and technical guidance for the prevention and control of coal wall spalling in SICS mining.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 3","pages":"1729 - 1750"},"PeriodicalIF":2.8,"publicationDate":"2025-01-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145162782","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Research on coal gangue transport and shield beam bearing pressure law based on FLAC3D-PFC3D coupled simulation and physical similar simulation","authors":"Zhenghan Qin,&nbsp;Yong Yuan,&nbsp;Zhenbin Mao,&nbsp;Xin Xu,&nbsp;Yong Li,&nbsp;Libao Li,&nbsp;Zhongshun Chen,&nbsp;Bo Li","doi":"10.1007/s40571-024-00899-5","DOIUrl":"10.1007/s40571-024-00899-5","url":null,"abstract":"<div><p>Top coal movement and real-time access to the caving state of top coal is the core scientific problem of intelligent longwall top coal caving, focusing on the pressure change law of the shield beam of hydraulic support during the process of top coal crushing and moving, the use of theoretical analysis, numerical simulation and physical simulation methods, research the rule of top coal movement and its impact on the change of the pressure of the shield beam. The stress state of shield beam under static state and coal caving state of hydraulic support is analyzed, and the results show that the force of the shield beam at rest is directly proportional to the support force of the column, and the moving state is positively correlated with the bearing capacity of the top beam. Based on the FDM-DEM coupling method, a FLAC^3D-PFC^3D coupling physical model was established, and the optimal initial coal discharge position was determined to be three steps away from the top coal caving, and the morphology surface of the top coal caving body gradually evolved from an ellipse with the long axis of 3.54 m and the short axis of 1.95 m to an ellipse with the long axis of 7.51 m and the short axis of 4.25 m, and the ellipse’s long axis was deflected by 18° in the direction of the pointer; and the position of the lowest point of line of demarcation between coal and gangue was linear and gradually deflected downward. Linear relationship and gradually downward offset, the offset point is located directly above the support shield beam. This phenomenon leads to the overall decrease of shield beam pressure, and the pressure of shield beam decreases rapidly and then increases slowly during the process of initial top coal caving and periodic coal caving. Finally, through the self-developed large-size physical similarity simulation platform for top coal caving, it was verified that in the cycle coal caving, the shield beam pressure showed the characteristic of rapid decrease and then slow increase, and the pressure before and after coal caving was generally reduced by 10%–30%. In the two research methods, the ratio of rapid change time and slow change time of shield beam pressure from the start of coal caving to the time when the gangue rate reaches 20% is approximately 5:2, so when the pressure of the shield beam is monitored to decrease during coal caving, it can be stopped when the time from the start of coal caving to the time of the lowest point of the pressure drop is 0.4 times, and at this time, the gangue rate is about 20%.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 3","pages":"1773 - 1793"},"PeriodicalIF":2.8,"publicationDate":"2025-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145162443","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"SPH-FE coupling for the simulation of confined flow through permeable deformable membranes","authors":"Matthias Brugger,&nbsp;Roland Traxl,&nbsp;Roman Lackner","doi":"10.1007/s40571-024-00892-y","DOIUrl":"10.1007/s40571-024-00892-y","url":null,"abstract":"<div><p>We present an extension of smoothed particle hydrodynamics (SPH) toward fluid flows involving the interaction with permeable deformable membranes. For this purpose, a coupled SPH-FE method based on a variational formulation of the immersed boundary (IB) method is developed. In the proposed method, weakly compressible SPH is used for the discretization of the fluid and a finite element (FE) method for thin structures for the discretization of the membrane. We consider confined flow in a two-dimensional fluid domain, with the membrane being represented as an elastic beam. Adopting the framework available for the IB method, the flux through the permeable membrane as described by Darcy’s law is considered. Finally, the proposed SPH-FE method is applied to two benchmark problems, i.e., the contraction of a circular membrane and the deformation of a membrane in a channel flow, comparing the numerical results with available analytical solutions.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 3","pages":"1665 - 1682"},"PeriodicalIF":2.8,"publicationDate":"2025-01-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s40571-024-00892-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161620","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"SPH simulation for 3D non-isothermal injection molding filling process using GPU acceleration","authors":"Yunpu Liu,&nbsp;Mengke Ren,&nbsp;Junfeng Gu,&nbsp;Zheng Li,&nbsp;Shilun Ruan,&nbsp;Changyu Shen","doi":"10.1007/s40571-024-00880-2","DOIUrl":"10.1007/s40571-024-00880-2","url":null,"abstract":"<div><p>The nature of high viscosity and pressure in the injection molding process poses a great challenge for numerical simulation in terms of numerical stability, especially when using particle-based meshless methods. In the present work, 3D filling stage of injection molding is simulated using smoothed particle hydrodynamics (SPH) method. To counter the instability caused by high viscosity and pressure, various methods including a new non-penetration boundary treatment, modified low-dissipation Riemann solver, kernel gradient correction and particle shift technique are applied. GPU parallel computing is achieved by using Taichi language to boost computing efficiency. 3D non-isothermal injection molding process is performed for rectangular cavity, tensile test specimen and a customized transparent injection mold which we intend to perform visual injection experiment to verify our simulation in future work. The properties of flow field such as pressure and velocity are shown and compared with Moldflow simulation. The results of our simulation show good agreement with Moldflow.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 2","pages":"1319 - 1333"},"PeriodicalIF":2.8,"publicationDate":"2025-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143919026","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Investigation of fracture energy rate for the combined finite-discrete element method","authors":"Peitao Li,&nbsp;Quansheng Liu,&nbsp;Lidan Fan,&nbsp;Yongqiang Yu,&nbsp;Feng Gao","doi":"10.1007/s40571-024-00809-9","DOIUrl":"10.1007/s40571-024-00809-9","url":null,"abstract":"<div><p>In the combined finite-discrete element method, the crack element was an important bridge to realize the transition from continuous to discontinuous deformation. The crack element was assumed as a non-thickness cohesive element, which could not be tested directly. The fracture energy rate for the crack element was an essential parameter determining the deformation and fracture process during the FDEM simulation. However, the common parameter calibration method required experience and time. Thus, a new estimation method of fracture energy rate was derived based on the energy mechanism. The proposed estimation method was consistent with the Griffith fracture criterion but also can reflect the influence of mechanical and inherent properties. Then, the influence of fracture energy rate on rock strength and fracture characteristic was analysed. The results showed that both the tensile strength and uniaxial compressive strength were increased with the fracture energy rate. A smaller tensile fracture energy rate would lead to a larger ratio of tensile fracture and a decrease in shear fracture. It was opposite for the shear fracture energy rate. Based on the rock strength and fracture characteristics, the initial coefficient of fracture energy rate should be 1.0. The fracture energy rate determined by the proposed estimation method and initial coefficient was close to the optimal value, which was well verified by the laboratory test result and tunnel engineering. Finally, the improved calibration steps of fracture energy rate for the complex material were put forward, which was expected for the fast and accurate calibration of model parameters.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 1","pages":"413 - 436"},"PeriodicalIF":2.8,"publicationDate":"2025-01-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143553843","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"DEM studies of cracking behavior of jointed rocks with different matching states under impact loading","authors":"Xiao Huaiguang,&nbsp;Yan Yatao,&nbsp;Wang Siwei","doi":"10.1007/s40571-024-00889-7","DOIUrl":"10.1007/s40571-024-00889-7","url":null,"abstract":"<div><p>Joints assume crucial roles in the propagation of stress waves for dynamic disasters. The SHPB system and jointed rocks with two matching states were numerically constructed. The mismatched joint comprises a triangular-shaped rough surface and a flat surface, whereas fully matched joint consists of rough surface and its biting surface. Subsequently, dynamic cracking behaviors and stress distribution of these jointed rocks subjected to different loading rates were examined. The simulated rock parameters are obtained based on cement mortar experiments. The results indicate that the dynamic strength of jointed rock is proportional to loading rate. Meanwhile, the strength and average modulus of mismatched jointed rocks (MJR) are significantly lower than those of fully matched jointed rocks (FJR) under the same loading. Stress concentration occurs at the tips of the mismatched joint and is readily influenced by the joint matching state. Owing to these rough joint surfaces, the MJR cracked first near the triangular teeth and expanded from the root of the teeth to both sides. However, the FJR was less affected by the joints and resembled that of intact rocks. Additionally, the effect of the joint undulation angle and the <i>P</i>-wave duration on the cracking behavior of jointed rock was also deliberated. The higher the joint undulation angle, the greater the concentrated stress.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 3","pages":"1617 - 1631"},"PeriodicalIF":2.8,"publicationDate":"2024-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145169842","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A high-resolution pseudo-polygon discrete element model for regional sea ice","authors":"Reshvar Kuppurangi,&nbsp;Min Wang","doi":"10.1007/s40571-024-00891-z","DOIUrl":"10.1007/s40571-024-00891-z","url":null,"abstract":"<div><p>This work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 3","pages":"1653 - 1664"},"PeriodicalIF":2.8,"publicationDate":"2024-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145168594","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Insights into the flow characteristics during hydraulic fracturing","authors":"Rezvan Abdi,&nbsp;Marek Krzaczek,&nbsp;Meisam Abdi","doi":"10.1007/s40571-024-00862-4","DOIUrl":"10.1007/s40571-024-00862-4","url":null,"abstract":"<div><p>This paper presents a numerical model to study fracture propagation during water-based hydraulic fracturing. To address the computational challenges associated with the numerical model, the proposed approach employs a set of overlapping spheres arranged in a monolayer to construct a porous specimen containing pre-existing cracks. The fluid-filled cracks represent various stages of initiation and propagation of fluid-driven fracture. The high-pressure fluid flow within the fractures is considered under isothermal conditions. Unlike the conventional focus on rock fracture analysis, the presented approach focuses on flow characteristics during fracture growth. The main objective of the presented study is to provide a detailed description of the computational fluid dynamics (CFD) aspects of fracture propagation during hydraulic fracturing to aid in calibration and validation of simplified discrete element method (DEM) models coupled with CFD representing this phenomenon. Experimental validations performed in previous studies support the model's reliability, making it useful in particular for calibration and validation of coupled 2D DEM-CFD models constructed from one layer of spheres. Obtaining experimental data for such cases is practically challenging, and the proposed model addresses the lack of reliable experimental data for hydraulic fracturing. To achieve this, representative specimens are designed, accurate simulations are conducted and precise assessments of the results are performed. Key variables such as density, pressure, velocity, porosity, and permeability were measured to facilitate the validation and calibration of future DEM-CFD studies.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 2","pages":"1139 - 1153"},"PeriodicalIF":2.8,"publicationDate":"2024-12-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s40571-024-00862-4.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143919002","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Experimental study and particle flow code numerical simulation of crack propagation and failure characteristics of prefabricated double-fissured rock samples","authors":"Zhaoyu Li,&nbsp;Bin Yang,&nbsp;Zhiguo Xia,&nbsp;Changxiang Wang,&nbsp;Tianqi Jiang,&nbsp;Zengxiang Lu,&nbsp;Jie Zhang","doi":"10.1007/s40571-024-00882-0","DOIUrl":"10.1007/s40571-024-00882-0","url":null,"abstract":"<div><p>Primary fissures within rock masses in underground engineering severely affect their mechanical properties. The propagation and coalescence of cracks are the key factors influencing the stability of engineering rock masses, thereby posing challenges to retaining the long-term stability of underground engineering rock masses. In this work, Brazilian splitting tests were carried out on sandstone samples with prefabricated fissures with different inclination angles. The strain field cloud map of the sample was obtained via digital speckle technology. The crack propagation evolution law, displacement field and stress field distribution characteristics of the fissured rock were studied from a microscopic perspective, and the failure mode of the sample was analysed in combination with a strain field cloud map. The results show that the mechanical parameter curve presents a “W” shape as the fissure inclination increases. The presence of fissures reduces the bearing capacity of rock samples, making the strength of fissured rock less than that of intact rock. The crack initiation position around the inclined fissure (F1) progressively transitions from the midpoint to the tip as the fissure inclination angle increases. The compressive stress and tensile stress concentration areas are primarily distributed near the fissure tip. As the fissure inclination angle increases, the maximum displacement value tends to decrease initially, followed by an increase. The rock bridge coalesces in the form of tensile cracks, and the fissure inclination angle alters the failure mode of the fissured rock.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 3","pages":"1565 - 1578"},"PeriodicalIF":2.8,"publicationDate":"2024-12-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145167194","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Towards the estimation of wall shear stress in smoothed particle hydrodynamics","authors":"Sumanta Laha,&nbsp;Georgios Fourtakas,&nbsp;Prasanta Kumar Das,&nbsp;Amir Keshmiri","doi":"10.1007/s40571-024-00879-9","DOIUrl":"10.1007/s40571-024-00879-9","url":null,"abstract":"<div><p>Over the past few decades, smoothed particle hydrodynamics (SPH) has emerged as an alternative computational fluid dynamics (CFD) technique, yet the estimation of wall shear stress lacks adequate standardisation. Wall shear stress is a critical metric in numerous applications, and hence, this is the focus of this paper. The present study proposes a novel SPH-based method for estimating wall shear stress using velocity data from the fluid particles adjacent to the wall. Wall shear stress is then calculated at the wall based on the wall shear stress data of the neighbouring fluid particles. For laminar flow, wall shear stress is estimated directly from velocity gradients, while for turbulent flow, the Smagorinsky large eddy simulation (LES) model with eddy viscosity is used. The results obtained from the model are rigorously validated against experimental, simulation and analytical data, confirming its effectiveness across different flow conditions. This validation highlights the reliability of the proposed model for fluid dynamics and bio-fluid mechanics research.</p></div>","PeriodicalId":524,"journal":{"name":"Computational Particle Mechanics","volume":"12 2","pages":"1309 - 1317"},"PeriodicalIF":2.8,"publicationDate":"2024-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s40571-024-00879-9.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143919068","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}