The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitati...The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitational potential of a prism. Discrete element method can be used to determine the gravitational potential of a prism, and can approximate the true gravitational potential values with sufficient accuracy (the smaller each element is, the more accurate the result is). Although Nagy's approach provided a closed expression, one does not know whether it is valid, due to the fact that this approach has not been confirmed in literatures. In this paper, a study on the comparison of Nagy's approach with discrete element method is presented. The results show that Nagy's formulas for determining the gravitational potential of a prism are valid in the domain both inside and outside the prism.展开更多
The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fractur...The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.展开更多
A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck ...A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.展开更多
A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in pol...A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geom...Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geometries of the fracture systems. However, no adequate efforts have been made to study this issue due to the current practical impossibility of laboratory tests with samples of large volumes con- taining many fractures, and the difficulty for controlling reliable initial and boundary conditions for large-scale in situ tests. Therefore, a reliable numerical predicting approach for evaluating anisotropy of fractured rock masses is needed. The objective of this study is to systematically investigate anisotropy of strength and deformability of fractured rocks, which has not been conducted in the past, using a nu- merical modeling method. A series of realistic two-dimensional （2D） discrete fracture network （DFN） models were established based on site investigation data, which were then loaded in different directions, using the code UDEC of discrete element method （DEM）, with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures. The results indicate that for design and safety assessments of rock engineering projects, the directional variations of strength and deformability of the fractured rock mass concerned must be treated properly with respect to the directions of in situ stresses. Traditional practice for simply positioning axial orientation of tunnels in association with principal stress directions only may not be adequate for safety requirements. Outstanding issues of the present study and su^zestions for future study are also oresented.展开更多
The electronic structures of BaMgF4 crystals containing an F colour centre are studied within the framework of the fully relativistic self-consistent Direc-Slater theory, using a numerically discrete variational （DV...The electronic structures of BaMgF4 crystals containing an F colour centre are studied within the framework of the fully relativistic self-consistent Direc-Slater theory, using a numerically discrete variational （DV-Xa） method. It is concluded from the calculated results that the energy levels of the F colour centre are located in the forbidden band. The optical transition energy from the ground state to the excited state for the F colour centre is about 5.12 eV, which corresponds to the 242-nm absorption band. These calculated results can explain the origin of the absorption bands.展开更多
The discrete element method is applied to investigate high-temperature spread in compacted metallic particle systems formed by high-velocity compaction. Assuming that heat transfer only occurs at contact zone between ...The discrete element method is applied to investigate high-temperature spread in compacted metallic particle systems formed by high-velocity compaction. Assuming that heat transfer only occurs at contact zone between particles, a discrete equation based on continuum mechanics is proposed to investigate the heat flux. Heat generated internally by friction between moving particles is determined by kinetic equations. For the proposed model, numerical results are obtained by a particle-flow-code-based program. Temperature profiles are determined at different locations and times. At a fixed location, the increase in temperature shows a logarithmic relationship with time. Investigation of three different systems indicates that the geometric distribution of the particulate material is one of the main influencing factors for the heat conduction process. Higher temperature is generated for denser packing, and vice versa. For smaller uniform particles, heat transfers more rapidly.展开更多
基金Supported by the National Natural Science Foundation of China (No.40637034, 40974015)the National 863 Program of China (No.2006AA12Z211)
文摘The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitational potential of a prism. Discrete element method can be used to determine the gravitational potential of a prism, and can approximate the true gravitational potential values with sufficient accuracy (the smaller each element is, the more accurate the result is). Although Nagy's approach provided a closed expression, one does not know whether it is valid, due to the fact that this approach has not been confirmed in literatures. In this paper, a study on the comparison of Nagy's approach with discrete element method is presented. The results show that Nagy's formulas for determining the gravitational potential of a prism are valid in the domain both inside and outside the prism.
文摘The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.
基金financial support from the National Natural Science Foundation of China (No. 51204181)the Research Fund for the Doctoral Program of Higher Education of China (No. 20110095120004)+2 种基金the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe Fundamental Research Funds for the Central Universities (Nos. 2011QNA10 and 2010QNB17)the China Postdoctoral Science Foundation (No. 20110491485)
文摘A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.
文摘A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
文摘Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geometries of the fracture systems. However, no adequate efforts have been made to study this issue due to the current practical impossibility of laboratory tests with samples of large volumes con- taining many fractures, and the difficulty for controlling reliable initial and boundary conditions for large-scale in situ tests. Therefore, a reliable numerical predicting approach for evaluating anisotropy of fractured rock masses is needed. The objective of this study is to systematically investigate anisotropy of strength and deformability of fractured rocks, which has not been conducted in the past, using a nu- merical modeling method. A series of realistic two-dimensional （2D） discrete fracture network （DFN） models were established based on site investigation data, which were then loaded in different directions, using the code UDEC of discrete element method （DEM）, with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures. The results indicate that for design and safety assessments of rock engineering projects, the directional variations of strength and deformability of the fractured rock mass concerned must be treated properly with respect to the directions of in situ stresses. Traditional practice for simply positioning axial orientation of tunnels in association with principal stress directions only may not be adequate for safety requirements. Outstanding issues of the present study and su^zestions for future study are also oresented.
基金supported by the Foundation of Shanghai Municipal Education Committee,China (Grant No. 09YZ210)the Shanghai Leading Academic Discipline Project (Grant No. S30502)
文摘The electronic structures of BaMgF4 crystals containing an F colour centre are studied within the framework of the fully relativistic self-consistent Direc-Slater theory, using a numerically discrete variational （DV-Xa） method. It is concluded from the calculated results that the energy levels of the F colour centre are located in the forbidden band. The optical transition energy from the ground state to the excited state for the F colour centre is about 5.12 eV, which corresponds to the 242-nm absorption band. These calculated results can explain the origin of the absorption bands.
基金S. Wang was supported by the Research Fund for Doctoral Program of Shandong Jianzhu University （Grant No. XNBS1338）and the National Natural Science Foundation of China （Grant No. 11471195）. Z. Zheng was supported by the National Natural Science Foundation of China （Grant Nos. 51174236 and 51134003）, the National Basic Research Program of China （Grant No. 2011 CB606306）, and the Opening Project of State Key Laboratory of Porous Metal Materials, China （Grant No. PMM-SKL-4-2012）.
文摘The discrete element method is applied to investigate high-temperature spread in compacted metallic particle systems formed by high-velocity compaction. Assuming that heat transfer only occurs at contact zone between particles, a discrete equation based on continuum mechanics is proposed to investigate the heat flux. Heat generated internally by friction between moving particles is determined by kinetic equations. For the proposed model, numerical results are obtained by a particle-flow-code-based program. Temperature profiles are determined at different locations and times. At a fixed location, the increase in temperature shows a logarithmic relationship with time. Investigation of three different systems indicates that the geometric distribution of the particulate material is one of the main influencing factors for the heat conduction process. Higher temperature is generated for denser packing, and vice versa. For smaller uniform particles, heat transfers more rapidly.
