The stability of slopes and tunnels is controlled by rock discontinuities,and the rock discontinuities roughness and the sliding direction play a signifcant role in shear failure.However,three-dimensional roughness ev...The stability of slopes and tunnels is controlled by rock discontinuities,and the rock discontinuities roughness and the sliding direction play a signifcant role in shear failure.However,three-dimensional roughness evaluation considering shear directions is scare,and the internal shear fracturing processes,micromechanical mechanisms and failure precursor of rock discontinuities are not well understood.Therefore,this study proposes a novel roughness evaluation index to quantitatively analyze the anisotropic characteristics of rock discontinuities.In conjunction with shear tests,a novel 3D-GBM modelling method considering the micromineral constituent and particle size distribution characteristics of granite as well as the geometric shape of discontinuities was realized.The strength,macro and micro-fracture characteristics,visual anisotropic shear evolution process and microfailure mechanism of granite discontinuities at diferent roughness and shear direction were investigated.Finally,the spatial and temporal evolutions of AE parameter b-value and magnitude M were further analyzed to reveal the shear fracture precursor of granite discontinuities.展开更多
Piled embankments are widely used in highway and railway engineering due to their economy and efficiency inovercoming several issues encountered in constructing embankments over weak soils. Soil arching, caused by the...Piled embankments are widely used in highway and railway engineering due to their economy and efficiency inovercoming several issues encountered in constructing embankments over weak soils. Soil arching, caused by the pile-subsoilrelative displacement （△s）, plays an important role in reducing the embankment load falling on weak soil, however, the funda-mental characteristics （e.g., formation and features） of soil arching remain poorly understood. In this study, a series of discreteelement method （DEM） modellings are performed to study the formation and features of soil arching with the variation of As inpiled embankments with or without geosynthetic reinforcement. Firstly, calibration for the modelling parameters is carried out bycomparing the DEM results with the experimental data obtained from the existing literature. Secondly, the analysis of the macro-and micro-behaviours is performed in detail. Finally, a parametric study is conducted in an effort to identify the influences of threekey factors on soil arching： the friction coefficient of the embankment fill （f）, the embankment height （h）, and the pile clear spacing（s-a）. Numerical results indicate that △s is a key factor governing the formation and features of soil arching in embankments. Tobe specific, soil arching gradually evolves from two inclined shear planes at a small △s to a hemispherical arch at a relatively largeAs. Then, with a continuous increase in △s, the soil arching height gradually increases and finally approaches a constant value of0.8（s-a） （i.e., the maximum soil arching height）. For a given case, the higher the soil arching height, the greater the degree of soilarching effect. The parametric study shows that the friction coefficient of the embankment fill has a negligible influence on theformation and features of soil arching. However, embankment height is a key factor governing the formation and features of soilarching. In addition, pile clear spacing has a significant effect on the formation of soil arching, but not on its features.展开更多
The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sand...The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sands were reproduced and the slump test and the direct shear test of the conditioned sands were implemented. A DEM equivalent model that can simulate the macro mechanical characteristic of the conditioned sands was proposed,and the corresponding numerical models of the slump test and the shear test were established. By selecting proper DEM model parameters,the errors of the slump values between the simulation results and the test results are in the range of 10.3%-14.3%,and the error of the curves between the shear displacement and the shear stress calculated with the DEM simulation is 4.68%-16.5% compared with that of the laboratory direct shear test. This illustrates that the proposed DEM equivalent model can approximately simulate the mechanical characteristics of the conditioned sands,which provides the basis for further simulation of the interaction between the conditioned soil and the chamber pressure system of the EPB machine.展开更多
A discrete element method was used to study the evolution of particle crushing in a rockfill sample subjected to triaxial shear. A simple procedure was developed to generate clusters with arbitrary shapes, which resem...A discrete element method was used to study the evolution of particle crushing in a rockfill sample subjected to triaxial shear. A simple procedure was developed to generate clusters with arbitrary shapes, which resembled real rockfill particles. A theoretical method was developed to define the failure criterion for an individual particle subjected to an arbitrary set of contact forces. Then, a series of numerical tests of large-scale drained triaxial tests were conducted to simulate the behaviors of the rockfill sample. Finally, we examined the development of micro-characteristics such as particle crushing, contact characteristics, porosity, deformation, movement, and energy dissipation. The simulation results were partially compared with the laboratory experiments, and good agreement was achieved, demonstrating that the particle crushing model proposed can be used to simulate the drained triaxial test ofrockfill materials. Based on a comparison of macro behaviors of the roekfill sample and micro structures of the particles, the microscopic mechanism of the rockfill materials subjected to triaxial shear was determined qualitatively. It is shown that the crushing rate, rather than the number of crushed particles, can be used to reflect the relationship between macro- and micro-mechanical characteristics of rockfill materials. These research results further develop our understanding of the deformation mechanism of rockfill materials.展开更多
The application of a combined finite-discrete element modeling approach to simu late the three-dimensional microscopic compaction behavior of single-layer met al powder system was described. The process was treated as...The application of a combined finite-discrete element modeling approach to simu late the three-dimensional microscopic compaction behavior of single-layer met al powder system was described. The process was treated as a static problem,wit h kinematical component being neglected. Due to ill condition,Cholesky’s metho d failed to solve the system equations,while conjugate gradient method was trie d and yielded good results. Deformation of the particles was examined and compar ed with the results of physical modeling experiments. In both cases,the inner p articles were deformed from sphere to polygonal column,with the edges turning f rom arc to straight line. The edge number of a particle was equal to the number of particles surrounding it. And the experiments show that the ductile metal par ticles can be densified only by their plastic deformation without the occurrence of rearrangement phenomenon.展开更多
Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were ex...Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.展开更多
Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid ...Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Discrete Element Method (DEM) is a powerful tool for simulating different types of mills. It also used for computing different types of particles such as rocks, grains, and molecules.</span></span><span...Discrete Element Method (DEM) is a powerful tool for simulating different types of mills. It also used for computing different types of particles such as rocks, grains, and molecules.</span></span><span style="white-space:normal;"><span style="font-family:""> </span></span><span style="white-space:normal;"><span style="font-family:"">DEM has been widely used in the field of rock mechanics. In the present work,</span></span><span style="white-space:normal;"><span style="font-family:""> </span></span><span style="white-space:normal;"><span style="font-family:"">DEM approach is applied to model the milling media (powder particles and balls) inside a planetary ball mill and to estimate the distribution of particles of a dry powder during milling. In fact, the efficiency of the DEM strongly depends on the input parameters. The DEM simulation results indicated that</span></span><span style="white-space:normal;"><span style="font-family:""> </span></span><span style="white-space:normal;"><span style="font-family:"">DEM is a promising tool for the simulation of the dynamic particles motion and interactions within planetary ball mill. These results could be utilized to further develop the synthesis performance, anticipate the reaction, and reduce the wear in the dry milling reactions.展开更多
Based on analysis of regularity of stacking coal,discrete element simultaneous simulation is adopted to predict the process of unloading coal,which is proved to be effcient in the prediction of ash content.The results...Based on analysis of regularity of stacking coal,discrete element simultaneous simulation is adopted to predict the process of unloading coal,which is proved to be effcient in the prediction of ash content.The results show that the altitude of new irregular coal is equal to the income coal volume divided by area of cabin.The distribution of infnitesimal flow velocity helps to induce the motion equation of infnitesimal element,which provides the mathematical model for computer simulation.Swarm,a computer programming language,is utilized in this study.Adaptive infnitesimal stacking algorithm helps settle the diffculties in attainment of infnitesimal elements.The result of simulation is similar to the actual situation,which can accurately predict the ash contents of current time and cumulative time.Coal movement in the cabin is a new project,the result of which can also be applied to other solid particles and the widespread of the result will be highly valued.展开更多
Many rock avalanches were triggered by the Wenchuan earthquake on May 12, 2008 in southwest China. Protection galleries covered with a single soil layer are usually used to protect against rockfall. Since one-layer pr...Many rock avalanches were triggered by the Wenchuan earthquake on May 12, 2008 in southwest China. Protection galleries covered with a single soil layer are usually used to protect against rockfall. Since one-layer protection galleries do not have sufficient buffer capacity, a two-layered absorbing system has been designed. This study aims to find whether an expanded poly-styrol （EPS） cushion, which is used in the soil-covered protection galleries for shock absorption, could be positioned under dynamic loadings. The dynamic impacts of the two-layered absorbing system under the conditions of rock avalanches are numerically simulated through a 2D discrete dement method. By selecting reasonable parameters, a series of numerical experiments were conducted to find the best combination for the two- layered absorbing system. The values of the EPS layer area as a percentage of the total area were set as 0% （Sl）, 22~ （S2）, and 70% （$3）. 22~ of the area of the EPS layer was found to be a reasonable value, and experiments were conducted to find the best position of the EPS layer in the two-layered absorbing system. The numerical results yield useful conclusions regarding the interaction between the impacting avalanches and the two-layered absorbing system. The soil layer can absorb the shock energy effectively and S2 （0.4-m thick EPS cushion covered with soil layer） is the most efficient combination, which can reduce the impact force, compared with the other combinations.展开更多
Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft...Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.展开更多
General purpose computing on GPU for scientific computing has been rapidly growing in recent years. We investigate the applicability of GPU to discrete element method (DEM) often used in particle motion simulation. NV...General purpose computing on GPU for scientific computing has been rapidly growing in recent years. We investigate the applicability of GPU to discrete element method (DEM) often used in particle motion simulation. NVIDIA provides a sample code for this type of simulation, which obtained superior performance than CPU in computational time. A computational model of the contact force in NVIDIA’s sample code is, however, too simple to use in practice. This paper modifies the NVIDIA’s simple model by replacing it with the practical model. The computing speed of the practical model on GPU is compared with the simple one on GPU and with the practical one on CPU in numerical experiments. The result shows that the practical model on GPU obtains the computing speed 6 times faster than the practical one on CPU while 7 times slower than that of the simple one on GPU. The effects of the GPU architectures on the computing speed are analyzed.展开更多
The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensio...The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.展开更多
1.Introduction and motivation Dry mixing is an important process step in the production of Li-ion battery electrodes.It is not only the first process step,but also the one with significant influence on the performance...1.Introduction and motivation Dry mixing is an important process step in the production of Li-ion battery electrodes.It is not only the first process step,but also the one with significant influence on the performance of the manufactured electrodes.Numerous authors show that the quality of manufactured electrodes correlate with the degree of comminution of the conductive additive carbon black(Bockholt et al.,2013;Mayer et al.,2022;Wenzel et al.,2014).展开更多
The efficiency of particle screening was studied over a range of vibrational parameters including amplitude, frequency and vibrational direction. The Discrete Element Method (DEM) was used to simulate the screening pr...The efficiency of particle screening was studied over a range of vibrational parameters including amplitude, frequency and vibrational direction. The Discrete Element Method (DEM) was used to simulate the screening process. A functional relationship between efficiency and the parameters, both singly and combined, is established. The function is a complicated exponential. Optimal amplitude and frequency values are smaller for particles near the mesh and larger for other particles. The optimum vibration angle is 45° for nearly all kinds of particles. A transverse velocity, V⊥, was defined and V⊥=0.2 m/s was identified to be the most efficient operating point by both simulation and experimental observation. Comparison of these results with those reported by others is included.展开更多
This paper presents numerical investigations into a ridged surface whose design is inspired by the geometry of a Farrer’sscallop.The objective of the performed research is to assess if the proposed Bioinspired Ridged...This paper presents numerical investigations into a ridged surface whose design is inspired by the geometry of a Farrer’sscallop.The objective of the performed research is to assess if the proposed Bioinspired Ridged Surface (BRS) can potentiallyimprove wear resistance of soil-engaging components used in agricultural machinery and to validate numerical simulationsperformed using software based on the Discrete Element Method (DEM).The wear performance of the BRS is experimentallydetermined and also compared with a conventional flat surface.Different size of soil particles and relative velocities between theabrasive sand and the testing surfaces are used.Comparative results show that the numerical simulations are in agreement withthe experimental results and support the hypothesis that abrasive wear is greatly reduced by substituting a conventional flatsurface with the BRS.展开更多
The accurate measurement of the fill level in the ball mill has not been resolved because of the interplay of many variable factors, which led the mill to be operated under the uneconomical condition and lost a lot of...The accurate measurement of the fill level in the ball mill has not been resolved because of the interplay of many variable factors, which led the mill to be operated under the uneconomical condition and lost a lot of energy. At present, some methods, such as vibration method and acoustic method, have been applied for measuring the fill level by the researchers. Aiming at the problem of the traditional methods for measuring the fill level, that is, the feature variables of the fill level suffer the influences of the ball load and the water content of the coal, a novel method to measure the fill level is proposed and a possible relation between the fill level and the angular position of the maximum vibration point on the mill shell is investigated. The angular positions of the maximum vibration point on the mill shell for different fill level cases are calculated theoretically under two assumptions, respectively. Meanwhile the charge motions of the mill for different fill level cases are simulated with the discrete element method （DEM）. And the simulation results are verified by comparing the motion trajectories of steel balls and power draft of the mill. The simulated movement trajectories of the outmost layer steel balls in the mill are monitored and analyzed to obtain the angular positions of the maximum vibration point on the mill shell. Both the results of the theoretical calculation and the 3D DEM simulation show that the position of the maximum vibration point on the mill shell moves to a lower angular positions as the fill level decreasing, which provides a new idea for measuring the filllevel accurately.展开更多
The discrete element method(DEM) has been widely used to simulate microscopic interactions between particles.Screening is a deeply complicated process when considering the law of motion for the particles,themselves.In...The discrete element method(DEM) has been widely used to simulate microscopic interactions between particles.Screening is a deeply complicated process when considering the law of motion for the particles,themselves.In this paper,a numerical model for the study of a particle screening process using the DEM is presented.Special attention was paid to the modeling of a vibrating screen that allows particles to pass through,or to rebound,when approaching the screen surface.Inferences concerning screen length and vibrating frequency as they relate to screening efficiency were studied.The conclusions were:three-dimensional simulation of screening efficiency along the screen length follows an exponential distribution;when the sieve vibrates over a certain frequency range the screening efficiency is stable;and,higher vibration frequencies can improve the handling capacity of the screening machine.展开更多
An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is bu...An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.展开更多
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.展开更多
基金Financial support to complete this study was provided in part by National Natural Science Foundation of China(52109119)the Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin(China Institute of Water Resources and Hydropower Research)(IWHR-SKL-202202)+1 种基金the China Postdoctoral Science Foundation Project(2022M723408)the Systematic Project of Guangxi Key Laboratory of Disaster Prevention and Engineering Safety(2020ZDK007)。
文摘The stability of slopes and tunnels is controlled by rock discontinuities,and the rock discontinuities roughness and the sliding direction play a signifcant role in shear failure.However,three-dimensional roughness evaluation considering shear directions is scare,and the internal shear fracturing processes,micromechanical mechanisms and failure precursor of rock discontinuities are not well understood.Therefore,this study proposes a novel roughness evaluation index to quantitatively analyze the anisotropic characteristics of rock discontinuities.In conjunction with shear tests,a novel 3D-GBM modelling method considering the micromineral constituent and particle size distribution characteristics of granite as well as the geometric shape of discontinuities was realized.The strength,macro and micro-fracture characteristics,visual anisotropic shear evolution process and microfailure mechanism of granite discontinuities at diferent roughness and shear direction were investigated.Finally,the spatial and temporal evolutions of AE parameter b-value and magnitude M were further analyzed to reveal the shear fracture precursor of granite discontinuities.
基金supported by the National Key Research and Development Program of China(2016YFC0800208)the National Natural Science Foundation of China(Nos.51278216,51478201,51308241,and 51608316)
文摘Piled embankments are widely used in highway and railway engineering due to their economy and efficiency inovercoming several issues encountered in constructing embankments over weak soils. Soil arching, caused by the pile-subsoilrelative displacement （△s）, plays an important role in reducing the embankment load falling on weak soil, however, the funda-mental characteristics （e.g., formation and features） of soil arching remain poorly understood. In this study, a series of discreteelement method （DEM） modellings are performed to study the formation and features of soil arching with the variation of As inpiled embankments with or without geosynthetic reinforcement. Firstly, calibration for the modelling parameters is carried out bycomparing the DEM results with the experimental data obtained from the existing literature. Secondly, the analysis of the macro-and micro-behaviours is performed in detail. Finally, a parametric study is conducted in an effort to identify the influences of threekey factors on soil arching： the friction coefficient of the embankment fill （f）, the embankment height （h）, and the pile clear spacing（s-a）. Numerical results indicate that △s is a key factor governing the formation and features of soil arching in embankments. Tobe specific, soil arching gradually evolves from two inclined shear planes at a small △s to a hemispherical arch at a relatively largeAs. Then, with a continuous increase in △s, the soil arching height gradually increases and finally approaches a constant value of0.8（s-a） （i.e., the maximum soil arching height）. For a given case, the higher the soil arching height, the greater the degree of soilarching effect. The parametric study shows that the friction coefficient of the embankment fill has a negligible influence on theformation and features of soil arching. However, embankment height is a key factor governing the formation and features of soilarching. In addition, pile clear spacing has a significant effect on the formation of soil arching, but not on its features.
基金Project (2007CB714006) supported by the National Basic Research Program of China
文摘The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sands were reproduced and the slump test and the direct shear test of the conditioned sands were implemented. A DEM equivalent model that can simulate the macro mechanical characteristic of the conditioned sands was proposed,and the corresponding numerical models of the slump test and the shear test were established. By selecting proper DEM model parameters,the errors of the slump values between the simulation results and the test results are in the range of 10.3%-14.3%,and the error of the curves between the shear displacement and the shear stress calculated with the DEM simulation is 4.68%-16.5% compared with that of the laboratory direct shear test. This illustrates that the proposed DEM equivalent model can approximately simulate the mechanical characteristics of the conditioned sands,which provides the basis for further simulation of the interaction between the conditioned soil and the chamber pressure system of the EPB machine.
基金supported by the National Key Basic Research Program of China (Grants No. 50879007 and 50979014)the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20090041110016)
文摘A discrete element method was used to study the evolution of particle crushing in a rockfill sample subjected to triaxial shear. A simple procedure was developed to generate clusters with arbitrary shapes, which resembled real rockfill particles. A theoretical method was developed to define the failure criterion for an individual particle subjected to an arbitrary set of contact forces. Then, a series of numerical tests of large-scale drained triaxial tests were conducted to simulate the behaviors of the rockfill sample. Finally, we examined the development of micro-characteristics such as particle crushing, contact characteristics, porosity, deformation, movement, and energy dissipation. The simulation results were partially compared with the laboratory experiments, and good agreement was achieved, demonstrating that the particle crushing model proposed can be used to simulate the drained triaxial test ofrockfill materials. Based on a comparison of macro behaviors of the roekfill sample and micro structures of the particles, the microscopic mechanism of the rockfill materials subjected to triaxial shear was determined qualitatively. It is shown that the crushing rate, rather than the number of crushed particles, can be used to reflect the relationship between macro- and micro-mechanical characteristics of rockfill materials. These research results further develop our understanding of the deformation mechanism of rockfill materials.
文摘The application of a combined finite-discrete element modeling approach to simu late the three-dimensional microscopic compaction behavior of single-layer met al powder system was described. The process was treated as a static problem,wit h kinematical component being neglected. Due to ill condition,Cholesky’s metho d failed to solve the system equations,while conjugate gradient method was trie d and yielded good results. Deformation of the particles was examined and compar ed with the results of physical modeling experiments. In both cases,the inner p articles were deformed from sphere to polygonal column,with the edges turning f rom arc to straight line. The edge number of a particle was equal to the number of particles surrounding it. And the experiments show that the ductile metal par ticles can be densified only by their plastic deformation without the occurrence of rearrangement phenomenon.
基金Supported by the China National Science and Technology Major Project(2016ZX05033002,2016ZX05033001).
文摘Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.
基金supported by the Natural Sciences and Engineering Research Council of Canada through Discovery Grant 341275 (G. Grasselli) and Engage EGP 461019-13
文摘Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
文摘Discrete Element Method (DEM) is a powerful tool for simulating different types of mills. It also used for computing different types of particles such as rocks, grains, and molecules.</span></span><span style="white-space:normal;"><span style="font-family:""> </span></span><span style="white-space:normal;"><span style="font-family:"">DEM has been widely used in the field of rock mechanics. In the present work,</span></span><span style="white-space:normal;"><span style="font-family:""> </span></span><span style="white-space:normal;"><span style="font-family:"">DEM approach is applied to model the milling media (powder particles and balls) inside a planetary ball mill and to estimate the distribution of particles of a dry powder during milling. In fact, the efficiency of the DEM strongly depends on the input parameters. The DEM simulation results indicated that</span></span><span style="white-space:normal;"><span style="font-family:""> </span></span><span style="white-space:normal;"><span style="font-family:"">DEM is a promising tool for the simulation of the dynamic particles motion and interactions within planetary ball mill. These results could be utilized to further develop the synthesis performance, anticipate the reaction, and reduce the wear in the dry milling reactions.
基金the financial support provided by the National Natural Science Foundation of China(No.51174202)Jiangsu Natural Science Foundation of China(No.20100095110013)
文摘Based on analysis of regularity of stacking coal,discrete element simultaneous simulation is adopted to predict the process of unloading coal,which is proved to be effcient in the prediction of ash content.The results show that the altitude of new irregular coal is equal to the income coal volume divided by area of cabin.The distribution of infnitesimal flow velocity helps to induce the motion equation of infnitesimal element,which provides the mathematical model for computer simulation.Swarm,a computer programming language,is utilized in this study.Adaptive infnitesimal stacking algorithm helps settle the diffculties in attainment of infnitesimal elements.The result of simulation is similar to the actual situation,which can accurately predict the ash contents of current time and cumulative time.Coal movement in the cabin is a new project,the result of which can also be applied to other solid particles and the widespread of the result will be highly valued.
基金financial support from the Project of National Science Foundation of China(Grant No.41272346)the National Outstanding Youth Funds(Grant No.41225011)+2 种基金financial support from the Science & Technology Research Plan of China Railway Eryuan Engineering Group CO.LTD (Grant No.13164196(13-15))the Project of National Science Foundation of China(Grant Nos. 41472293,91430105)"hundred talents" program of CAS
文摘Many rock avalanches were triggered by the Wenchuan earthquake on May 12, 2008 in southwest China. Protection galleries covered with a single soil layer are usually used to protect against rockfall. Since one-layer protection galleries do not have sufficient buffer capacity, a two-layered absorbing system has been designed. This study aims to find whether an expanded poly-styrol （EPS） cushion, which is used in the soil-covered protection galleries for shock absorption, could be positioned under dynamic loadings. The dynamic impacts of the two-layered absorbing system under the conditions of rock avalanches are numerically simulated through a 2D discrete dement method. By selecting reasonable parameters, a series of numerical experiments were conducted to find the best combination for the two- layered absorbing system. The values of the EPS layer area as a percentage of the total area were set as 0% （Sl）, 22~ （S2）, and 70% （$3）. 22~ of the area of the EPS layer was found to be a reasonable value, and experiments were conducted to find the best position of the EPS layer in the two-layered absorbing system. The numerical results yield useful conclusions regarding the interaction between the impacting avalanches and the two-layered absorbing system. The soil layer can absorb the shock energy effectively and S2 （0.4-m thick EPS cushion covered with soil layer） is the most efficient combination, which can reduce the impact force, compared with the other combinations.
基金Special Fund of Marine Commonweal Industry under contact Nos 201105016 and 201205007supported by National Marine Environment Forecasting Centrethe National Natural Science Foundation of China under contact No.41176012
文摘Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.
文摘General purpose computing on GPU for scientific computing has been rapidly growing in recent years. We investigate the applicability of GPU to discrete element method (DEM) often used in particle motion simulation. NVIDIA provides a sample code for this type of simulation, which obtained superior performance than CPU in computational time. A computational model of the contact force in NVIDIA’s sample code is, however, too simple to use in practice. This paper modifies the NVIDIA’s simple model by replacing it with the practical model. The computing speed of the practical model on GPU is compared with the simple one on GPU and with the practical one on CPU in numerical experiments. The result shows that the practical model on GPU obtains the computing speed 6 times faster than the practical one on CPU while 7 times slower than that of the simple one on GPU. The effects of the GPU architectures on the computing speed are analyzed.
文摘The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.
基金the Federal Ministry of Education and Research for funding this work in the project Sim4Pro(grant number 03XP0242B)within the ProZell Cluster.
文摘1.Introduction and motivation Dry mixing is an important process step in the production of Li-ion battery electrodes.It is not only the first process step,but also the one with significant influence on the performance of the manufactured electrodes.Numerous authors show that the quality of manufactured electrodes correlate with the degree of comminution of the conductive additive carbon black(Bockholt et al.,2013;Mayer et al.,2022;Wenzel et al.,2014).
基金the Special Topic of Key Science and Technology of Fujian Province Fund (No.2006HZ0002-2)
文摘The efficiency of particle screening was studied over a range of vibrational parameters including amplitude, frequency and vibrational direction. The Discrete Element Method (DEM) was used to simulate the screening process. A functional relationship between efficiency and the parameters, both singly and combined, is established. The function is a complicated exponential. Optimal amplitude and frequency values are smaller for particles near the mesh and larger for other particles. The optimum vibration angle is 45° for nearly all kinds of particles. A transverse velocity, V⊥, was defined and V⊥=0.2 m/s was identified to be the most efficient operating point by both simulation and experimental observation. Comparison of these results with those reported by others is included.
基金supported by the National High Technology Research and Development Program of China(863 Program)(Grant No.2009AA043603-4,2009AA043604-2)by National Foundations of Agricultural Technological Transformation of China(Grant No.2009GB23600507)+1 种基金by National Natural Science Foundation of China(Grant No.50675087,50635030)by "985 Project" of Jilin University
文摘This paper presents numerical investigations into a ridged surface whose design is inspired by the geometry of a Farrer’sscallop.The objective of the performed research is to assess if the proposed Bioinspired Ridged Surface (BRS) can potentiallyimprove wear resistance of soil-engaging components used in agricultural machinery and to validate numerical simulationsperformed using software based on the Discrete Element Method (DEM).The wear performance of the BRS is experimentallydetermined and also compared with a conventional flat surface.Different size of soil particles and relative velocities between theabrasive sand and the testing surfaces are used.Comparative results show that the numerical simulations are in agreement withthe experimental results and support the hypothesis that abrasive wear is greatly reduced by substituting a conventional flatsurface with the BRS.
基金supported by National Natural Science Foundation of China (Grant No. 50775035)New Doctor Teacher Foundation of Southeast University of China (Grant No. 9202000024)
文摘The accurate measurement of the fill level in the ball mill has not been resolved because of the interplay of many variable factors, which led the mill to be operated under the uneconomical condition and lost a lot of energy. At present, some methods, such as vibration method and acoustic method, have been applied for measuring the fill level by the researchers. Aiming at the problem of the traditional methods for measuring the fill level, that is, the feature variables of the fill level suffer the influences of the ball load and the water content of the coal, a novel method to measure the fill level is proposed and a possible relation between the fill level and the angular position of the maximum vibration point on the mill shell is investigated. The angular positions of the maximum vibration point on the mill shell for different fill level cases are calculated theoretically under two assumptions, respectively. Meanwhile the charge motions of the mill for different fill level cases are simulated with the discrete element method （DEM）. And the simulation results are verified by comparing the motion trajectories of steel balls and power draft of the mill. The simulated movement trajectories of the outmost layer steel balls in the mill are monitored and analyzed to obtain the angular positions of the maximum vibration point on the mill shell. Both the results of the theoretical calculation and the 3D DEM simulation show that the position of the maximum vibration point on the mill shell moves to a lower angular positions as the fill level decreasing, which provides a new idea for measuring the filllevel accurately.
基金Project 2006HZ0002-2 supported by the Special Topic Fund of Key Science and Technology of Fujian Province
文摘The discrete element method(DEM) has been widely used to simulate microscopic interactions between particles.Screening is a deeply complicated process when considering the law of motion for the particles,themselves.In this paper,a numerical model for the study of a particle screening process using the DEM is presented.Special attention was paid to the modeling of a vibrating screen that allows particles to pass through,or to rebound,when approaching the screen surface.Inferences concerning screen length and vibrating frequency as they relate to screening efficiency were studied.The conclusions were:three-dimensional simulation of screening efficiency along the screen length follows an exponential distribution;when the sieve vibrates over a certain frequency range the screening efficiency is stable;and,higher vibration frequencies can improve the handling capacity of the screening machine.
基金supported by the National Natural Science Foundation of China(Grant No.51105386)the National Deep-Sea Technology Project of Development and Research(Grant No.DYXM-115-04-02-01)the Fundamental Research Funds for the Central Universities(Grant No.2011QNZT058)
文摘An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.
基金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.
