Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
The thick layer and graded feeding technology of a belt roasting machine is an effective method for improving the production efficiency and quality index of pellet production,and a reasonable design of the mechanical ...The thick layer and graded feeding technology of a belt roasting machine is an effective method for improving the production efficiency and quality index of pellet production,and a reasonable design of the mechanical structure of the layer is the basis for optimizing the heat and mass transfer performance of the layer.Janssen effect and von Mises yield criterion were used to establish a simplified mathematical model describing the elastic and plastic deformation of the green pellet under the action of an external force.The mechanical characteristics of extrusion,contact,and elastic-plastic deformation between green pellet particles in the material layer of the belt roasting machine were modeled using EDEM software.For a green pellet size of 12 mm,as the layer height increases from 300 to 1000 mm,the maximum vertical pressure on the pellets increases from 11.64 to 24.01 N,and the porosity decreases from 27.04%to 22.01%.As the layer height increases,the contact between the green pellets becomes more intense,and the force chain structure of the layer becomes more stable;the Janssen effect is observed when the layer reaches 700 mm.The compressive strength of the green pellets is linearly related to the particle size,and the compressive strength increases with an increase in particle size.At a layer height of 600 mm,as the particle size of the green pellets increases from 8 to 20 mm,the maximum vertical pressure increases from 7.54 to 44.16 N,and the porosity increases from 23.20%to 31.47%,while the yield per unit of the layer decreased by 12.1%.Small particles have a more stable force chain structure,larger comparative area,and higher production efficiency;however,their compressive strength is lower.Large particles have higher compressive strength and good permeability in the layer,but the production efficiency is relatively low.In actual production,a variety of factors should be integrated to optimize the feeding,and a multi-granularity graded feeding is the most ideal feeding.展开更多
Despite the prevalence and validity of the universal distinct element code(UDEC)in simulations in geotechnics domain,water-weakening process of rock models remains elusive.Prior research has made positive contribution...Despite the prevalence and validity of the universal distinct element code(UDEC)in simulations in geotechnics domain,water-weakening process of rock models remains elusive.Prior research has made positive contributions to a presupposed link between modelling parameters and saturation degree,Sr.Nevertheless,this effort presents inaccurate results and limited implications owing to the misleading interpretation,that is,devoid of the basic logic in UDEC that modelling parameters should be calibrated by tested macroscopic properties in contrast to a presupposed relation with Sr.To fill this gap,a new methodology is proposed by coupling a computationally efficient parametric study with the simulation of water-weakening mechanisms.More specifically,tested macroscopic properties with different Sr values are input into parametric relations to acquire initial modelling parameters that are sequentially calibrated and modulated until simulations are in line with geomechanical tests.Illustrative example reveals that numerical water-weakening effects on macroscopic properties,mechanical behaviours,and failure configurations are highly consistent with tested ones with noticeable computational expediency,implying the feasibility and simplicity of this methodology.Furthermore,with compatibility across various numerical models,the proposed methodology substantially extends the applicability of UDEC in simulating water-weakening geotechnical problems.展开更多
This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on...This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on the geological data in Horn River Basin,Northeast British Columbia,Canada.The modeling results indicate that the maximum magnitude of seismic events appears at the fracturing stage.The increment of fluid volume in the fault determines the cumulative moment and maximum fault slippage,both of which are essentially proportional to the fluid volume.After backflow starts,the fluid near the joint intersection keeps flowing into the critically stressed fault,rather than backflows to the wellbore.Although fault slippage is affected by the changes of both pore pressure and ambient rock stress,their contributions are different at fracturing and backflow stages.At fracturing stage,pore pressure change shows a dominant effect on induced fault slippage.While at backflow stage,because the fault plane is under a critical stress state,any minor disturbance would trigger a fault slippage.The energy analysis indicates that aseismic deformation takes up a majority of the total deformation energy during hydraulic fracturing.A common regularity is found in both fracturing-and backflow-induced seismicity that the cumulative moment and maximum fault slippage are nearly proportional to the injected fluid volume.This study shows some novel insights into interpreting fracturing-and backflowinduced seismicity,and provides useful information for controlling and mitigating seismic hazards due to hydraulic fracturing.展开更多
To retrofit and strengthen existing unreinforced masonry (URM) structures to resist the potential earthquake damages has become an important issue in Australia. In order to secure the performance of URM under seismic ...To retrofit and strengthen existing unreinforced masonry (URM) structures to resist the potential earthquake damages has become an important issue in Australia. In order to secure the performance of URM under seismic loading in the future, a research project was carried out aimed at developing a simple and high strength seismic retrofitting technique for masonry structures. A series of experimental testing on URM walls retrofitted with an innovative technique by cable system have been conducted. The results indicated that both the strength and ductility of the tested speci-mens were significantly enhanced with the technique. An analytical model which is based on Dis-tinct Element Method (DEM) has also been developed to simulate the behaviour of URM walls be-fore and after retrofitting. The model is then further developed by applying a seismic wave to the wall to simulate the wall behavior under earthquake loads before and after retrofitting.展开更多
Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surfac...Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.展开更多
The quantitative evaluation of errors involved in a particular numerical modelling is of prime importance for the effectiveness and reliability of the method. Errors in Distinct Element Modelling are generated mainly ...The quantitative evaluation of errors involved in a particular numerical modelling is of prime importance for the effectiveness and reliability of the method. Errors in Distinct Element Modelling are generated mainly through three resources as simplification of physical model, determination of parameters and boundary conditions. A measure of errors which represent the degree of numerical solution 'close to true value' is proposed through fuzzy probability in this paper. The main objective of this paper is to estimate the reliability of Distinct Element Method in rock engineering practice by varying the parameters and boundary conditions. The accumulation laws of standard errors induced by improper determination of parameters and boundary conditions are discussed in delails. Furthermore, numerical experiments are given to illustrate the estimation of fuzzy reliability. Example shows that fuzzy reliability falls between 75%-98% when the relative standard errors of input data is under 10 %.展开更多
The Australian love of 'heritage' buildings (most of them are unreinforced masonry (URM)) means that greater attention is required to secure their performance under seismic or impact loading in the future. A r...The Australian love of 'heritage' buildings (most of them are unreinforced masonry (URM)) means that greater attention is required to secure their performance under seismic or impact loading in the future. A research project has been carried out to develop a new, economic and high strength retrofitting technique for masonry structures. A series of experimental testing on URM walls retrofitted with an innovative technique by cable system have been conducted. In this paper, an analytical model which is based on distinct element method (DEM) is developed to simulate the behaviour of retrofitted walls. In DEM, a solid is represented as an assembly of discrete blocks. Joints are modelled as interface between distinct bodies. It is a dynamic process and specially designed to model the behaviour of discontinuities. In order to assist the practising engineers to design this new retrofitted wall system, a simple mechanical model was also developed to predict the strength of the retrofitted walls. The results obtained from this simple mechanical model are compared with those from both experiments and distinct element model.展开更多
This letter mainly aims to investigate the mechanical behavior within shear bands in regolith both under the Earth and the Moon conditions via the distinct element method, in which a novel contact model considering in...This letter mainly aims to investigate the mechanical behavior within shear bands in regolith both under the Earth and the Moon conditions via the distinct element method, in which a novel contact model considering interparticle van der Waals forces and rolling resistance is employed. The results show that for regolith under both conditions the stress paths are almost identical inside and outside the shear bands but void ratio, average pure rotation rate, and strain paths are rather distinct with dilation, particle rotation and the second invariant of strain tensor mainly occurring within the bands. However, the regolith under the Moon condition has higher peak strength and more significant strain localization than those under the Earth condition.展开更多
The reaction of earth to pull-out process of frictional rock bolts was here modeled by the distinct element method (DEM). Ten frictional bolts were prepared;the expanding shells of five bolts included convex edges and...The reaction of earth to pull-out process of frictional rock bolts was here modeled by the distinct element method (DEM). Ten frictional bolts were prepared;the expanding shells of five bolts included convex edges and the others had the shells with concave bits. The strength of bolts was measured by applying a standard pull-out test;the results confirmed that the strength of shells with convex edges was remarkably more than the strength of other shells. Furthermore, a two-dimensional DEM model of the test was developed by a particle flow code;the obtained results showed that the reaction of rock particles to the contacts occurring between the convex edges and earth was considerably more than those of the concave bits. In the other words, the convex edges transferred the pull-out force into a large area of the surrounded rock, causing these bolts to have the highest resistance against earth movements.展开更多
Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P...Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.展开更多
The deformation and fracture evolution mechanisms of the strata overlying mines mined using sublevel caving were studied via numerical simulations.Moreover,an expression for the normal force acting on the side face of...The deformation and fracture evolution mechanisms of the strata overlying mines mined using sublevel caving were studied via numerical simulations.Moreover,an expression for the normal force acting on the side face of a steeply dipping superimposed cantilever beam in the surrounding rock was deduced based on limit equilibrium theory.The results show the following:(1)surface displacement above metal mines with steeply dipping discontinuities shows significant step characteristics,and(2)the behavior of the strata as they fail exhibits superimposition characteristics.Generally,failure first occurs in certain superimposed strata slightly far from the goaf.Subsequently,with the constant downward excavation of the orebody,the superimposed strata become damaged both upwards away from and downwards toward the goaf.This process continues until the deep part of the steeply dipping superimposed strata forms a large-scale deep fracture plane that connects with the goaf.The deep fracture plane generally makes an angle of 12°-20°with the normal to the steeply dipping discontinuities.The effect of the constant outward transfer of strata movement due to the constant outward failure of the superimposed strata in the metal mines with steeply dipping discontinuities causes the scope of the strata movement in these mines to be larger than expected.The strata in the metal mines with steeply dipping discontinuities mainly show flexural toppling failure.However,the steeply dipping structural strata near the goaf mainly exhibit shear slipping failure,in which case the mechanical model used to describe them can be simplified by treating them as steeply dipping superimposed cantilever beams.By taking the steeply dipping superimposed cantilever beam that first experiences failure as the key stratum,the failure scope of the strata(and criteria for the stability of metal mines with steeply dipping discontinuities mined using sublevel caving)can be obtained via iterative computations from the key stratum,moving downward toward and upwards away from the goaf.展开更多
This paper investigates the explicit use of rock reinforcement in a discontinuous stress analysis model.A series of numerical experiments was undertaken to evaluate the performance of local and global reinforcement mo...This paper investigates the explicit use of rock reinforcement in a discontinuous stress analysis model.A series of numerical experiments was undertaken to evaluate the performance of local and global reinforcement models implemented in universal distinct element code(UDEC).This was made possible by calibrating the reinforcement models to the laboratory behavior of a fully-grouted rebar bolt tested under pure pull and pure shear loading conditions.The model calibration focuses on matching different loading stages of the force-displacement curve including the initial elastic response,the hardening behavior and the bolt rupture.The paper concludes with a discussion on the suitability of the different reinforcement models in UDEC including their advantages and limitations.Finally,it addresses the choice of input parameters required for a realistic simulation of fully-grouted rebar bolts.展开更多
Based on the stress redistribution analysis of rock mass during the deep underground excavation, the unloading process of pre-flawed rock material was simulated by distinct element method (DEM). The effects of unloadi...Based on the stress redistribution analysis of rock mass during the deep underground excavation, the unloading process of pre-flawed rock material was simulated by distinct element method (DEM). The effects of unloading rate and flaw inclination angle on unloading strengths and cracking properties of pre-flawed rock specimens are numerically revealed. The results indicate that the unloading failure strength of pre-flawed specimen exhibits a power-function increase trend with the increase of unloading period. Moreover, combined with the stress state analysis on the flaws, it is found that the unloading failure strength increases with the increase of flaw inclination angle. The cracking distribution of pre-flawed specimens under the unloading condition closely depends on the flaw inclination angle, and three typical types of flaw coalescence are observed. Furthermore, at a faster unloading rate, the pre-flawed specimen experiences a sharper and quicker unloading failure process, resulting in more splitting cracks in the specimens.展开更多
From the geological structure of the columnar jointed rock mass, a visual model was established in software AUTOCAD by programming based on the algorithm of the Voronoi diagram. Furthermore, a program to convert the A...From the geological structure of the columnar jointed rock mass, a visual model was established in software AUTOCAD by programming based on the algorithm of the Voronoi diagram. Furthermore, a program to convert the AUTOCAD model into 3DEC (3-dimensional distinct element code) model was developed, and a numerical model was established in 3DEC. Moreover, the results of triaxial compression tests of columnar jointed rock masses were simulated numerically. The REV (representative element volume) scale was studied, and the result shows that the REV size is 3 m × 3 m. The proposed approach, the established model and the numerical simulation were applied to study the macro-mechanical properties and the equivalent strength parameters of the columnar jointed rock mass. The numerical simulation results are in good accordance with the in-situ test results.展开更多
Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrest...Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrestrial laser scanning(TLS)with distinct element method for rock mass characterization and stability analysis in tunnels.TLS records detailed geometric information of the surrounding rock mass by scanning and collecting the positions of millions of rock surface points without contact.By conducting a fuzzy K-means method,a discontinuity automatic identification algorithm was developed,and a method for obtaining the geometric parameters of discontinuities was proposed.This method permits the user to visually identify each discontinuity and acquire its spatial distribution features(e.g.occurrences,spac-ings,trace lengths)in great detail.Compared with hand mapping in conventional geotechnical surveys,the geometric information of discontinuities obtained by this approach is more accurate and the iden-tification is more efficient.Then,a discrete fracture network with the same statistical characteristics as the actual discontinuities was generated with the distinct element method,and a representative nu-merical model of the jointed surrounding rock mass was established.By means of numerical simulation,potential unstable rock blocks were assessed,and failure mechanisms were analyzed.This method was applied to detection and assessment of unstable rock blocks in the spillway and sand flushing tunnel of the Hongshiyan hydropower project after a collapse.The results show that the noncontact detection of blocks was more labor-saving with lower safety risks compared with manual surveys,and the stability assessment was more reliable since the numerical model built by this method was more consistent with the distribution characteristics of actual joints.This study can provide a reference for geological survey and unstable rock block hazard mitigation in tunnels subjected to complex geology and active rockfalls.展开更多
Based on a centrifuge model test and distinct element method(DEM), this study provides new insights into the uplift response of a shallow-buried structure and the liquefaction mechanism for saturated sand around the s...Based on a centrifuge model test and distinct element method(DEM), this study provides new insights into the uplift response of a shallow-buried structure and the liquefaction mechanism for saturated sand around the structure under seismic action. In the centrifuge test, a high-speed microscopic camera was installed in the structure model, by which the movements of particles around the structure were monitored. Then, a two-dimensional digital image processing technology was used to analyze the microstructure of saturated sand during the shaking event. Herein, a numerical simulation of the centrifuge experiment was conducted using a two-phase(solid and fl uid) fully coupled distinct element code. This code incorporates a particle-fl uid coupling model by means of a "fi xed coarse-grid" fl uid scheme in PFC3D(Particle Flow Code in Three Dimensions), with the modeling parameters partially calibrated based on earlier studies. The physical and numerical models both indicate the uplifts of the shallow-buried structure and the sharp rise in excess pore pressure. The corresponding micro-scale responses and explanations are provided. Overall, the uplift response of an underground structure and the occurrence of liquefaction in saturated sand are predicted successfully by DEM modeling. However, the dynamic responses during the shaking cannot be modeled accurately due to the restricted computer power.展开更多
The effect of intermediate stress(in situ tunnel axial)on a strainburst is studied with a threedimensional(3D)bonded block distinct element method(DEM).A series of simulations of strainbursts under true triaxial in si...The effect of intermediate stress(in situ tunnel axial)on a strainburst is studied with a threedimensional(3D)bonded block distinct element method(DEM).A series of simulations of strainbursts under true triaxial in situ stress conditions(i.e.high tangential stress,moderate intermediate stress and low radial stress)of near-boundary rock masses are performed.Compared with the experimental results,the DEM model is able to capture the stress-strain response,failure pattern and energy balance of strainbursts.The fracturing processes of strainbursts are also numerically reproduced.Numerical results show that,as the intermediate stress increases:(1)The peak strain of strainbursts increases,the yield stress increases,the rock strength increases linearly,and the ratio of yield stress to rock strength decreases,indicating that the precursory information on strainbursts is enhanced;(2)Tensile and shear cracks increase significantly,and slabbing and bending of rock plates are more pronounced;and(3)The stored elastic strain energy and dissipated energy increase linearly,whereas the kinetic energy of the ejected rock fragments increases approximately exponentially,implying an increase in strainburst intensity.By comparing the experimental and numerical results,the effect of intermediate stress on the rock strength of strainbursts is discussed in order to address three key issues.Then,the Mogi criterion is applied to construct new strength criteria for strainbursts by converting the one-face free true triaxial stress state of a strainburst to its equivalent true triaxial stress state.In summary,the effect of intermediate stress on strainbursts is a double-edged sword that can enhance the rock strength and the precursory information of a strainburst,but also increase its intensity.展开更多
To solve the problem of water loss during mining of shallow,buried coal seams,we have first analyzed the mechanism and suitability of solid-liquid coupling,i.e.,we used the FLUID-MECHANICS system of 3-Dimensional Dist...To solve the problem of water loss during mining of shallow,buried coal seams,we have first analyzed the mechanism and suitability of solid-liquid coupling,i.e.,we used the FLUID-MECHANICS system of 3-Dimensional Distinct Element Code(3DEC)in simulating dynamic water crannies in overlying strata,under mining conditions of a large longwall coalface.Next the dynamic initiation of a water cranny,its propagation and close phases were studied with 3DEC,along with the overlying strata breakage and recombination as the mining space of the shallow,buried coal seam increased.Combined with the change in the stress and displacement fields,the distribution features of the mining cranny were systematically studied.The effect of regularities and their effective measures of local filling and mine slicing technology in controlling mine crannies were investigated and the potential danger areas of water loss identified.Our results can be applied to decrease water loss during the exploitation of shallow,buried coal seams with a thin bedrock.The results also prove that 3DEC is a credible numerical analytical method to predict initiations of dynamic water crannies,their propagation,their closure phases and other concomitant hazards.展开更多
The section of shield tunnel of the Chengdu Metro line passes primarily through sandy cobble strata. There are many buildings with spread foundations along the lines. Shield tunnel construction will disturb the ground...The section of shield tunnel of the Chengdu Metro line passes primarily through sandy cobble strata. There are many buildings with spread foundations along the lines. Shield tunnel construction will disturb the ground, causing displacement or stress to adjacent spread foundations. Based on the similarity theory, a laboratory model test of shield tunnel driving was carried out to study the influence of shield tunnel excavation on the displace ment of adjacent spread foundation. The results show that foundation closer to the tunnel has greater displacement or settlement than that further away. The horizontal dis placement is small and is influenced greatly by the cutting face. The displacement along the machine driving direction is bigger and is significantly affected by the thrust force. Settlement occurs primarily when shield machine passes close to the foundation and is the greatest at that time. Uneven settlement at the bottom of the spread foundation reaches a maximum after the excavation ends. In a numerical simulation, a particle flow model was con structed to study the impact of shield tunnel excavation on the stresses in the ground. The model showed stress con centration at the bottom of the spread foundation. With the increasing ground loss ratio, a loose area appears in the tunnel dome where the contact force dropped. Above the loose area, the contact force increases, forming an arch shaped soil area which prevents the loose area from expanding to the ground surface. The excavation also changed the pressure distribution around spread foundation.展开更多
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
基金support from the Outstanding Youth Fund of Anhui Province(Grant No.2208085J19).
文摘The thick layer and graded feeding technology of a belt roasting machine is an effective method for improving the production efficiency and quality index of pellet production,and a reasonable design of the mechanical structure of the layer is the basis for optimizing the heat and mass transfer performance of the layer.Janssen effect and von Mises yield criterion were used to establish a simplified mathematical model describing the elastic and plastic deformation of the green pellet under the action of an external force.The mechanical characteristics of extrusion,contact,and elastic-plastic deformation between green pellet particles in the material layer of the belt roasting machine were modeled using EDEM software.For a green pellet size of 12 mm,as the layer height increases from 300 to 1000 mm,the maximum vertical pressure on the pellets increases from 11.64 to 24.01 N,and the porosity decreases from 27.04%to 22.01%.As the layer height increases,the contact between the green pellets becomes more intense,and the force chain structure of the layer becomes more stable;the Janssen effect is observed when the layer reaches 700 mm.The compressive strength of the green pellets is linearly related to the particle size,and the compressive strength increases with an increase in particle size.At a layer height of 600 mm,as the particle size of the green pellets increases from 8 to 20 mm,the maximum vertical pressure increases from 7.54 to 44.16 N,and the porosity increases from 23.20%to 31.47%,while the yield per unit of the layer decreased by 12.1%.Small particles have a more stable force chain structure,larger comparative area,and higher production efficiency;however,their compressive strength is lower.Large particles have higher compressive strength and good permeability in the layer,but the production efficiency is relatively low.In actual production,a variety of factors should be integrated to optimize the feeding,and a multi-granularity graded feeding is the most ideal feeding.
基金supported by the National Natural Science Foundation of China under Grant Nos.41977249 and 42090052the China Scholarship Council under file No.202204910040.
文摘Despite the prevalence and validity of the universal distinct element code(UDEC)in simulations in geotechnics domain,water-weakening process of rock models remains elusive.Prior research has made positive contributions to a presupposed link between modelling parameters and saturation degree,Sr.Nevertheless,this effort presents inaccurate results and limited implications owing to the misleading interpretation,that is,devoid of the basic logic in UDEC that modelling parameters should be calibrated by tested macroscopic properties in contrast to a presupposed relation with Sr.To fill this gap,a new methodology is proposed by coupling a computationally efficient parametric study with the simulation of water-weakening mechanisms.More specifically,tested macroscopic properties with different Sr values are input into parametric relations to acquire initial modelling parameters that are sequentially calibrated and modulated until simulations are in line with geomechanical tests.Illustrative example reveals that numerical water-weakening effects on macroscopic properties,mechanical behaviours,and failure configurations are highly consistent with tested ones with noticeable computational expediency,implying the feasibility and simplicity of this methodology.Furthermore,with compatibility across various numerical models,the proposed methodology substantially extends the applicability of UDEC in simulating water-weakening geotechnical problems.
基金supported by the Key Innovation Team Program of Innovation Talents Promotion Plan by Ministry of Science and Technology of China(Grant No.2016RA4059)National Natural Science Foundation of China(Grant Nos.41672268 and 41772286)。
文摘This paper presents a three-dimensional fully hydro-mechanical coupled distinct element study on fault reactivation and induced seismicity due to hydraulic fracturing injection and subsequent backflow process,based on the geological data in Horn River Basin,Northeast British Columbia,Canada.The modeling results indicate that the maximum magnitude of seismic events appears at the fracturing stage.The increment of fluid volume in the fault determines the cumulative moment and maximum fault slippage,both of which are essentially proportional to the fluid volume.After backflow starts,the fluid near the joint intersection keeps flowing into the critically stressed fault,rather than backflows to the wellbore.Although fault slippage is affected by the changes of both pore pressure and ambient rock stress,their contributions are different at fracturing and backflow stages.At fracturing stage,pore pressure change shows a dominant effect on induced fault slippage.While at backflow stage,because the fault plane is under a critical stress state,any minor disturbance would trigger a fault slippage.The energy analysis indicates that aseismic deformation takes up a majority of the total deformation energy during hydraulic fracturing.A common regularity is found in both fracturing-and backflow-induced seismicity that the cumulative moment and maximum fault slippage are nearly proportional to the injected fluid volume.This study shows some novel insights into interpreting fracturing-and backflowinduced seismicity,and provides useful information for controlling and mitigating seismic hazards due to hydraulic fracturing.
文摘To retrofit and strengthen existing unreinforced masonry (URM) structures to resist the potential earthquake damages has become an important issue in Australia. In order to secure the performance of URM under seismic loading in the future, a research project was carried out aimed at developing a simple and high strength seismic retrofitting technique for masonry structures. A series of experimental testing on URM walls retrofitted with an innovative technique by cable system have been conducted. The results indicated that both the strength and ductility of the tested speci-mens were significantly enhanced with the technique. An analytical model which is based on Dis-tinct Element Method (DEM) has also been developed to simulate the behaviour of URM walls be-fore and after retrofitting. The model is then further developed by applying a seismic wave to the wall to simulate the wall behavior under earthquake loads before and after retrofitting.
文摘Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.
文摘The quantitative evaluation of errors involved in a particular numerical modelling is of prime importance for the effectiveness and reliability of the method. Errors in Distinct Element Modelling are generated mainly through three resources as simplification of physical model, determination of parameters and boundary conditions. A measure of errors which represent the degree of numerical solution 'close to true value' is proposed through fuzzy probability in this paper. The main objective of this paper is to estimate the reliability of Distinct Element Method in rock engineering practice by varying the parameters and boundary conditions. The accumulation laws of standard errors induced by improper determination of parameters and boundary conditions are discussed in delails. Furthermore, numerical experiments are given to illustrate the estimation of fuzzy reliability. Example shows that fuzzy reliability falls between 75%-98% when the relative standard errors of input data is under 10 %.
文摘The Australian love of 'heritage' buildings (most of them are unreinforced masonry (URM)) means that greater attention is required to secure their performance under seismic or impact loading in the future. A research project has been carried out to develop a new, economic and high strength retrofitting technique for masonry structures. A series of experimental testing on URM walls retrofitted with an innovative technique by cable system have been conducted. In this paper, an analytical model which is based on distinct element method (DEM) is developed to simulate the behaviour of retrofitted walls. In DEM, a solid is represented as an assembly of discrete blocks. Joints are modelled as interface between distinct bodies. It is a dynamic process and specially designed to model the behaviour of discontinuities. In order to assist the practising engineers to design this new retrofitted wall system, a simple mechanical model was also developed to predict the strength of the retrofitted walls. The results obtained from this simple mechanical model are compared with those from both experiments and distinct element model.
基金supported by China National Funds for Distinguished Young Scientists (51025932)the National Natural Science Foundation of China (51179128)Program of Shanghai Academic Chief Scientist (11XD1405200)
文摘This letter mainly aims to investigate the mechanical behavior within shear bands in regolith both under the Earth and the Moon conditions via the distinct element method, in which a novel contact model considering interparticle van der Waals forces and rolling resistance is employed. The results show that for regolith under both conditions the stress paths are almost identical inside and outside the shear bands but void ratio, average pure rotation rate, and strain paths are rather distinct with dilation, particle rotation and the second invariant of strain tensor mainly occurring within the bands. However, the regolith under the Moon condition has higher peak strength and more significant strain localization than those under the Earth condition.
文摘The reaction of earth to pull-out process of frictional rock bolts was here modeled by the distinct element method (DEM). Ten frictional bolts were prepared;the expanding shells of five bolts included convex edges and the others had the shells with concave bits. The strength of bolts was measured by applying a standard pull-out test;the results confirmed that the strength of shells with convex edges was remarkably more than the strength of other shells. Furthermore, a two-dimensional DEM model of the test was developed by a particle flow code;the obtained results showed that the reaction of rock particles to the contacts occurring between the convex edges and earth was considerably more than those of the concave bits. In the other words, the convex edges transferred the pull-out force into a large area of the surrounded rock, causing these bolts to have the highest resistance against earth movements.
基金financially supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.
基金Financial support for this work was provided by the Youth Fund Program of the National Natural Science Foundation of China (No. 42002292)the General Program of the National Natural Science Foundation of China (No. 42377175)the General Program of the Hubei Provincial Natural Science Foundation, China (No. 2023AFB631)
文摘The deformation and fracture evolution mechanisms of the strata overlying mines mined using sublevel caving were studied via numerical simulations.Moreover,an expression for the normal force acting on the side face of a steeply dipping superimposed cantilever beam in the surrounding rock was deduced based on limit equilibrium theory.The results show the following:(1)surface displacement above metal mines with steeply dipping discontinuities shows significant step characteristics,and(2)the behavior of the strata as they fail exhibits superimposition characteristics.Generally,failure first occurs in certain superimposed strata slightly far from the goaf.Subsequently,with the constant downward excavation of the orebody,the superimposed strata become damaged both upwards away from and downwards toward the goaf.This process continues until the deep part of the steeply dipping superimposed strata forms a large-scale deep fracture plane that connects with the goaf.The deep fracture plane generally makes an angle of 12°-20°with the normal to the steeply dipping discontinuities.The effect of the constant outward transfer of strata movement due to the constant outward failure of the superimposed strata in the metal mines with steeply dipping discontinuities causes the scope of the strata movement in these mines to be larger than expected.The strata in the metal mines with steeply dipping discontinuities mainly show flexural toppling failure.However,the steeply dipping structural strata near the goaf mainly exhibit shear slipping failure,in which case the mechanical model used to describe them can be simplified by treating them as steeply dipping superimposed cantilever beams.By taking the steeply dipping superimposed cantilever beam that first experiences failure as the key stratum,the failure scope of the strata(and criteria for the stability of metal mines with steeply dipping discontinuities mined using sublevel caving)can be obtained via iterative computations from the key stratum,moving downward toward and upwards away from the goaf.
基金supported by the Natural Science and Engineering Council of Canada
文摘This paper investigates the explicit use of rock reinforcement in a discontinuous stress analysis model.A series of numerical experiments was undertaken to evaluate the performance of local and global reinforcement models implemented in universal distinct element code(UDEC).This was made possible by calibrating the reinforcement models to the laboratory behavior of a fully-grouted rebar bolt tested under pure pull and pure shear loading conditions.The model calibration focuses on matching different loading stages of the force-displacement curve including the initial elastic response,the hardening behavior and the bolt rupture.The paper concludes with a discussion on the suitability of the different reinforcement models in UDEC including their advantages and limitations.Finally,it addresses the choice of input parameters required for a realistic simulation of fully-grouted rebar bolts.
基金Projects(41630642,11472311)supported by the National Natural Science Foundation of ChinaProject(2017zzts181)supported by the Cultivating Excellent Ph Ds of Central South University,ChinaProject(201806370062)supported by the China Scholarship Council
文摘Based on the stress redistribution analysis of rock mass during the deep underground excavation, the unloading process of pre-flawed rock material was simulated by distinct element method (DEM). The effects of unloading rate and flaw inclination angle on unloading strengths and cracking properties of pre-flawed rock specimens are numerically revealed. The results indicate that the unloading failure strength of pre-flawed specimen exhibits a power-function increase trend with the increase of unloading period. Moreover, combined with the stress state analysis on the flaws, it is found that the unloading failure strength increases with the increase of flaw inclination angle. The cracking distribution of pre-flawed specimens under the unloading condition closely depends on the flaw inclination angle, and three typical types of flaw coalescence are observed. Furthermore, at a faster unloading rate, the pre-flawed specimen experiences a sharper and quicker unloading failure process, resulting in more splitting cracks in the specimens.
基金Projects(50911130366, 50979030) supported by the National Natural Science Foundation of China
文摘From the geological structure of the columnar jointed rock mass, a visual model was established in software AUTOCAD by programming based on the algorithm of the Voronoi diagram. Furthermore, a program to convert the AUTOCAD model into 3DEC (3-dimensional distinct element code) model was developed, and a numerical model was established in 3DEC. Moreover, the results of triaxial compression tests of columnar jointed rock masses were simulated numerically. The REV (representative element volume) scale was studied, and the result shows that the REV size is 3 m × 3 m. The proposed approach, the established model and the numerical simulation were applied to study the macro-mechanical properties and the equivalent strength parameters of the columnar jointed rock mass. The numerical simulation results are in good accordance with the in-situ test results.
基金support of the National Natural Science Foundation of China(Grant No.42102316)the Open Project of the Technology Innovation Center for Geological Environment Monitoring of Ministry of Natural Resources of China(Grant No.2022KFK1212005).
文摘Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrestrial laser scanning(TLS)with distinct element method for rock mass characterization and stability analysis in tunnels.TLS records detailed geometric information of the surrounding rock mass by scanning and collecting the positions of millions of rock surface points without contact.By conducting a fuzzy K-means method,a discontinuity automatic identification algorithm was developed,and a method for obtaining the geometric parameters of discontinuities was proposed.This method permits the user to visually identify each discontinuity and acquire its spatial distribution features(e.g.occurrences,spac-ings,trace lengths)in great detail.Compared with hand mapping in conventional geotechnical surveys,the geometric information of discontinuities obtained by this approach is more accurate and the iden-tification is more efficient.Then,a discrete fracture network with the same statistical characteristics as the actual discontinuities was generated with the distinct element method,and a representative nu-merical model of the jointed surrounding rock mass was established.By means of numerical simulation,potential unstable rock blocks were assessed,and failure mechanisms were analyzed.This method was applied to detection and assessment of unstable rock blocks in the spillway and sand flushing tunnel of the Hongshiyan hydropower project after a collapse.The results show that the noncontact detection of blocks was more labor-saving with lower safety risks compared with manual surveys,and the stability assessment was more reliable since the numerical model built by this method was more consistent with the distribution characteristics of actual joints.This study can provide a reference for geological survey and unstable rock block hazard mitigation in tunnels subjected to complex geology and active rockfalls.
基金National Natural Science Foundation of China under Grant Nos.41272296 and 51208294
文摘Based on a centrifuge model test and distinct element method(DEM), this study provides new insights into the uplift response of a shallow-buried structure and the liquefaction mechanism for saturated sand around the structure under seismic action. In the centrifuge test, a high-speed microscopic camera was installed in the structure model, by which the movements of particles around the structure were monitored. Then, a two-dimensional digital image processing technology was used to analyze the microstructure of saturated sand during the shaking event. Herein, a numerical simulation of the centrifuge experiment was conducted using a two-phase(solid and fl uid) fully coupled distinct element code. This code incorporates a particle-fl uid coupling model by means of a "fi xed coarse-grid" fl uid scheme in PFC3D(Particle Flow Code in Three Dimensions), with the modeling parameters partially calibrated based on earlier studies. The physical and numerical models both indicate the uplifts of the shallow-buried structure and the sharp rise in excess pore pressure. The corresponding micro-scale responses and explanations are provided. Overall, the uplift response of an underground structure and the occurrence of liquefaction in saturated sand are predicted successfully by DEM modeling. However, the dynamic responses during the shaking cannot be modeled accurately due to the restricted computer power.
基金We acknowledge the funding support from the National Natural Science Foundation of China(Grant Nos.52009016 and 52179118)the Fundamental Research Funds for the Central Universities(Grant No.2022QN1032).
文摘The effect of intermediate stress(in situ tunnel axial)on a strainburst is studied with a threedimensional(3D)bonded block distinct element method(DEM).A series of simulations of strainbursts under true triaxial in situ stress conditions(i.e.high tangential stress,moderate intermediate stress and low radial stress)of near-boundary rock masses are performed.Compared with the experimental results,the DEM model is able to capture the stress-strain response,failure pattern and energy balance of strainbursts.The fracturing processes of strainbursts are also numerically reproduced.Numerical results show that,as the intermediate stress increases:(1)The peak strain of strainbursts increases,the yield stress increases,the rock strength increases linearly,and the ratio of yield stress to rock strength decreases,indicating that the precursory information on strainbursts is enhanced;(2)Tensile and shear cracks increase significantly,and slabbing and bending of rock plates are more pronounced;and(3)The stored elastic strain energy and dissipated energy increase linearly,whereas the kinetic energy of the ejected rock fragments increases approximately exponentially,implying an increase in strainburst intensity.By comparing the experimental and numerical results,the effect of intermediate stress on the rock strength of strainbursts is discussed in order to address three key issues.Then,the Mogi criterion is applied to construct new strength criteria for strainbursts by converting the one-face free true triaxial stress state of a strainburst to its equivalent true triaxial stress state.In summary,the effect of intermediate stress on strainbursts is a double-edged sword that can enhance the rock strength and the precursory information of a strainburst,but also increase its intensity.
基金Projects 50374065,50490273,50474068 supported by the National Natural Science Foundation of ChinaNCET-05-0480 by the New Century Excellent Talentsin University+1 种基金07KF09 by the 2007 Research Fund of the State Key Laboratory of Coal Resources and Mine SafetyCUMT and 2005B002 by the Scientific Re-search Foundation of China University of Mining&Technology
文摘To solve the problem of water loss during mining of shallow,buried coal seams,we have first analyzed the mechanism and suitability of solid-liquid coupling,i.e.,we used the FLUID-MECHANICS system of 3-Dimensional Distinct Element Code(3DEC)in simulating dynamic water crannies in overlying strata,under mining conditions of a large longwall coalface.Next the dynamic initiation of a water cranny,its propagation and close phases were studied with 3DEC,along with the overlying strata breakage and recombination as the mining space of the shallow,buried coal seam increased.Combined with the change in the stress and displacement fields,the distribution features of the mining cranny were systematically studied.The effect of regularities and their effective measures of local filling and mine slicing technology in controlling mine crannies were investigated and the potential danger areas of water loss identified.Our results can be applied to decrease water loss during the exploitation of shallow,buried coal seams with a thin bedrock.The results also prove that 3DEC is a credible numerical analytical method to predict initiations of dynamic water crannies,their propagation,their closure phases and other concomitant hazards.
基金funded by the National Natural Science Foundation of China (Nos. 51278422 and 50925830)the National 973 Plan Topics of China (No. 2010CB732105)+1 种基金the National Science and Technology Pillar Program of China (No. 2012BAG05B03)the Sichuan Youth Science and Technology Foundation, China (No. 2012JQ0021)
文摘The section of shield tunnel of the Chengdu Metro line passes primarily through sandy cobble strata. There are many buildings with spread foundations along the lines. Shield tunnel construction will disturb the ground, causing displacement or stress to adjacent spread foundations. Based on the similarity theory, a laboratory model test of shield tunnel driving was carried out to study the influence of shield tunnel excavation on the displace ment of adjacent spread foundation. The results show that foundation closer to the tunnel has greater displacement or settlement than that further away. The horizontal dis placement is small and is influenced greatly by the cutting face. The displacement along the machine driving direction is bigger and is significantly affected by the thrust force. Settlement occurs primarily when shield machine passes close to the foundation and is the greatest at that time. Uneven settlement at the bottom of the spread foundation reaches a maximum after the excavation ends. In a numerical simulation, a particle flow model was con structed to study the impact of shield tunnel excavation on the stresses in the ground. The model showed stress con centration at the bottom of the spread foundation. With the increasing ground loss ratio, a loose area appears in the tunnel dome where the contact force dropped. Above the loose area, the contact force increases, forming an arch shaped soil area which prevents the loose area from expanding to the ground surface. The excavation also changed the pressure distribution around spread foundation.
