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 %.展开更多
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.展开更多
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 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.展开更多
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.展开更多
A numerical method for simulating the stability of particle-packing structures is presented. The packing structures were modeled on the basis of face-centered cubic (fcc) and body-centered cubic (bcc) structures, and ...A numerical method for simulating the stability of particle-packing structures is presented. The packing structures were modeled on the basis of face-centered cubic (fcc) and body-centered cubic (bcc) structures, and the stability of these structures was investigated using the distinct element method. The interaction between the particles was simplified by considering repulsive, adhesive, and damping forces, and the stability against the gravitational force was simulated. The results under a certain set of parameters showed characteristic deformation when the particles were arranged in an fcc array. Focusing on the local structure, the resulting model was divided into several domains: The bottom base, four top corners, and intermediate domains. The bottom base notably became a body-centered tetragonal (bct) structure, which corresponds to a uniaxially compressed bcc structure. Conversely, the models based on the bcc arrangement were structurally stable, as no specific deformation was observed, and a monotonously compressed bct structure was obtained. Consequently, the bcc arrangement is concluded to be more stable against uniaxial compression, such as the gravitational force, in a particle-packing system.展开更多
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.展开更多
Reliable estimates of slope stability are essential for safe design and planning of road cut hill slopes which accommo-date a number of tourist destinations around the world. The failure of cut slopes along these hill...Reliable estimates of slope stability are essential for safe design and planning of road cut hill slopes which accommo-date a number of tourist destinations around the world. The failure of cut slopes along these hills puts human life in grave danger and it is also disastrous for the economy. In the present study, a section of 100 m high jointed basalt hill slope has been analyzed numerically in a distinct element code, which is apt for simulating the behavior for jointed rock. The analysis was carried out for both the dry and saturated conditions. The distinct element analysis of the hill slope demonstrates it to be marginally stable under dry condition, while for the saturated condition, the hill slope fails along well defined joint planes.展开更多
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.展开更多
降压开采会造成影响深海能源土宏观力学特性发生显著变化,而胶结型是天然气水合物赋存于深海能源土储层的主要类型,因此降压开采条件下胶结型深海能源土力学特性研究对于提高天然气水合物开采效率,预防海洋灾害具有重要意义。研究了胶...降压开采会造成影响深海能源土宏观力学特性发生显著变化,而胶结型是天然气水合物赋存于深海能源土储层的主要类型,因此降压开采条件下胶结型深海能源土力学特性研究对于提高天然气水合物开采效率,预防海洋灾害具有重要意义。研究了胶结型深海能源土降压开采时的胶结和损伤影响,构建了三维离散元(Discrete Element Method,DEM)数值模型,并分析了围压、饱和度和损伤对于降压开采条件下胶结型深海能源土力学性能影响。研究结论对于深海能源土储层中天然气水合物开采研究具有较好的借鉴意义。展开更多
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.展开更多
文摘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 %.
文摘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.
基金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 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.
文摘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.
文摘A numerical method for simulating the stability of particle-packing structures is presented. The packing structures were modeled on the basis of face-centered cubic (fcc) and body-centered cubic (bcc) structures, and the stability of these structures was investigated using the distinct element method. The interaction between the particles was simplified by considering repulsive, adhesive, and damping forces, and the stability against the gravitational force was simulated. The results under a certain set of parameters showed characteristic deformation when the particles were arranged in an fcc array. Focusing on the local structure, the resulting model was divided into several domains: The bottom base, four top corners, and intermediate domains. The bottom base notably became a body-centered tetragonal (bct) structure, which corresponds to a uniaxially compressed bcc structure. Conversely, the models based on the bcc arrangement were structurally stable, as no specific deformation was observed, and a monotonously compressed bct structure was obtained. Consequently, the bcc arrangement is concluded to be more stable against uniaxial compression, such as the gravitational force, in a particle-packing system.
基金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.
文摘Reliable estimates of slope stability are essential for safe design and planning of road cut hill slopes which accommo-date a number of tourist destinations around the world. The failure of cut slopes along these hills puts human life in grave danger and it is also disastrous for the economy. In the present study, a section of 100 m high jointed basalt hill slope has been analyzed numerically in a distinct element code, which is apt for simulating the behavior for jointed rock. The analysis was carried out for both the dry and saturated conditions. The distinct element analysis of the hill slope demonstrates it to be marginally stable under dry condition, while for the saturated condition, the hill slope fails along well defined joint planes.
基金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.
文摘降压开采会造成影响深海能源土宏观力学特性发生显著变化,而胶结型是天然气水合物赋存于深海能源土储层的主要类型,因此降压开采条件下胶结型深海能源土力学特性研究对于提高天然气水合物开采效率,预防海洋灾害具有重要意义。研究了胶结型深海能源土降压开采时的胶结和损伤影响,构建了三维离散元(Discrete Element Method,DEM)数值模型,并分析了围压、饱和度和损伤对于降压开采条件下胶结型深海能源土力学性能影响。研究结论对于深海能源土储层中天然气水合物开采研究具有较好的借鉴意义。
基金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.
