Optimizing Bucket Elevator Performance through a Blend of Discrete Element Method, Response Surface Methodology, and Firefly Algorithm Approaches

This research introduces a novel approach to enhancing bucket elevator design and operation through the integration of discrete element method(DEM)simulation,design of experiments(DOE),and metaheuristic optimization algorithms.Specifically,the study employs the firefly algorithm(FA),a metaheuristic optimization technique,to optimize bucket elevator parameters for maximizing transport mass and mass flow rate discharge of granular materials under specified working conditions.The experimental methodology involves several key steps:screening experiments to identify significant factors affecting bucket elevator operation,central composite design(CCD)experiments to further explore these factors,and response surface methodology(RSM)to create predictive models for transport mass and mass flow rate discharge.The FA algorithm is then applied to optimize these models,and the results are validated through simulation and empirical experiments.The study validates the optimized parameters through simulation and empirical experiments,comparing results with DEM simulation.The outcomes demonstrate the effectiveness of the FA algorithm in identifying optimal bucket parameters,showcasing less than 10%and 15%deviation for transport mass and mass flow rate discharge,respectively,between predicted and actual values.Overall,this research provides insights into the critical factors influencing bucket elevator operation and offers a systematic methodology for optimizing bucket parameters,contributing to more efficient material handling in various industrial applications.

Distinct element modelling of fracture plan control in continuum and jointed rock mass in presplitting method of surface mining

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.

AN ESTIMATION OF FUZZY RELIABILITY OF DISTINCT ELEMENT METHOD IN PREDICTION OF SURFACE SUBSIDENCE DUE TO COAL MINING

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 %.

Distinct element method investigation on mechanical behavior within shear bands in granulates under the Earth and the Moon conditions

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.

Modelling of URM Walls Retrofitted with Cable: A Comparison Between a Basic Mechanical Model and Distinct Element Method

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.

Investigating Earth Reaction to Pull-Out Process of Frictional Rock Bolts Using Distinct Element Method

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.

Modeling and Simulation of Particle-Packing Structures and Their Stability Using the Distinct Element Method

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.

基于CFD-DEM耦合的梯级溜槽的设计与分析

针对传统物料转运过程中溜槽和输送带磨损严重、出口处粉尘浓度过高的问题,建立含臂架的梯级溜槽几何模型,采用基于计算流体力学与离散单元法(CFD-DEM)耦合的数值模拟方法,分析了臂架对转运溜槽的磨损以及对其出口处粉尘排放浓度的影响。仿真结果表明:含臂架的梯级溜槽可以有效控制物料流的速度和方向,降低对溜槽内表面的冲击磨损,降低出口处的粉尘量。

Influence of heterogeneity on rock strength and stiffness using discrete element method and parallel bond model

The particulate discrete element method(DEM) can be employed to capture the response of rock,provided that appropriate bonding models are used to cement the particles to each other.Simulations of laboratory tests are important to establish the extent to which those models can capture realistic rock behaviors.Hitherto the focus in such comparison studies has either been on homogeneous specimens or use of two-dimensional(2D) models.In situ rock formations are often heterogeneous,thus exploring the ability of this type of models to capture heterogeneous material behavior is important to facilitate their use in design analysis.In situ stress states are basically three-dimensional(3D),and therefore it is important to develop 3D models for this purpose.This paper revisits an earlier experimental study on heterogeneous specimens,of which the relative proportions of weaker material(siltstone) and stronger,harder material(sandstone) were varied in a controlled manner.Using a 3D DEM model with the parallel bond model,virtual heterogeneous specimens were created.The overall responses in terms of variations in strength and stiffness with different percentages of weaker material(siltstone) were shown to agree with the experimental observations.There was also a good qualitative agreement in the failure patterns observed in the experiments and the simulations,suggesting that the DEM data enabled analysis of the initiation of localizations and micro fractures in the specimens.

New Method to Measure the Fill Level of the Ball Mill I——Theoretical Analysis and DEM Simulation

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.

Three-dimensional distinct element modeling of fault reactivation and induced seismicity due to hydraulic fracturing injection and backflow

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.

A review of methods,applications and limitations for incorporating fluid flow in the discrete element method

The past decade has witnessed the substantial growth in research interests and progress on the subject of coupled hydro-mechanical processes in rocks and soils,driven mainly by the surge of research in unconventional hydrocarbon reservoirs and associated hazards.Many coupling techniques have been developed to include the effects of fluid flow in the discrete element method(DEM),and the techniques have been applied to a variety of geomechanical problems.Although these coupling methods have been successfully applied in various engineering fields,no single fluid/DEM coupling method is universal due to the complexity of engineering problems and the limitations of the numerical methods.For researchers and engineers,the key to solve a specific problem is to select the most appropriate fluid/DEM coupling method among these modeling technologies.The purpose of this paper is to give a comprehensive review of fluid flow/DEM coupling methods and relevant research.Given their importance,the availability or unavailability of best practice guidelines is outlined.The theoretical background and current status of DEM are introduced first,and the principles,applications,and advantages and disadvantages of different fluid flow/DEM coupling methods are discussed.Finally,a summary with speculation on future development trends is given.

Study on the particle breakage of ballast based on a GPU accelerated discrete element method

Breakage of particles will have greatly influence on mechanical behavior of granular material(GM)under external loads,such as ballast,rockfill and sand.The discrete element method(DEM)is one of the most popular methods for simulating GM as each particle is represented on its own.To study breakage mechanism of particle breakage,a cohesive contact mode is developed based on the GPU accelerated DEM code-Blaze-DEM.A database of the 3D geometry model of rock blocks is established based on the 3D scanning method.And an agglomerate describing the rock block with a series of non-overlapping spherical particles is used to build the DEM numerical model of a railway ballast sample,which is used to the DEM oedometric test to study the particles’breakage characteristics of the sample under external load.Furthermore,to obtain the meso-mechanical parameters used in DEM,a black-analysis method is used based on the laboratory tests of the rock sample.Based on the DEM numerical tests,the particle breakage process and mechanisms of the railway ballast are studied.All results show that the developed code can better used for large scale simulation of the particle breakage analysis of granular material.

NUMERICAL ANALYSIS OF SOIL PLUG INSIDE SUCTION FOUNDATIONS DURING SUCTION PENETRATION BY DISCRETE ELEMENT METHOD

The phenomenon of the soil plug usually rising inside the suction foundations during suction penetration was quantitatively described and predicted. The formation process of the soil plug was simulated and calculated by DEM （discrete element method） model. The seepage flow, the self-weight of soil, the friction on the chamber wall as well as the suction inside the chamber are considered as the main external forces in the process. The results are compared with a set of laboratory model tests performed by using three soil types （sand, silty clay and clay） in the Bohai Sea area. The heights of soil plug from numerical estimations are lower than those from model test results, mainly because the suction pressure and friction resistance are applied in an ideal way under the numerical simulation.

A new method to simulate irregular particles by discrete element method

Granular materials are ubiquitous in nature and important in various applications such as road and railway engineering. Granular materials exhibit complicated mechanical behaviors, which are affected significantly by the irregular shape of particles. Currently, the discrete element method (DEM) has been accepted as an effective approach to investigate the mechanical behaviors of granular materials. However, there are scarce simulations based on DEM in literatures considering the irregularity of particle shape. A new method is proposed to simulate individual real particle with irregular shape using clump constituted by overlapping spheres. First, the geometric model of real particle with surface nodes and inner nodes is established through digitally processing the computerized tomography (CT) scanning data. Second, a clump consisting of spheres is generated to simulate the real particle using a minimum distance criterion. The criterion is implemented by tree optimization algorithm. Influential factors are also introduced to balance the model accuracy and computing cost. Effects of the influential factors, including the density of geometric grid and the minimum distance, on simulations are discussed. Results show that this new method is simpler and more efficient than the previous methods in terms of the model accuracy and computing cost.

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method(DEM)has been intensively used to study the constitutive behaviour of granular materials.However,to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear.This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials.Three groups of models with spherical and clumped particles are investigated from four perspectives:(i)deviatoric stress and volumetric behaviour;(ii)critical state behaviour;(iii)stress-dilatancy relationship;and(iv)the evolution of principal stress ratio against axial strain.The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses,volumetric behaviour,and critical state behaviour observed in experiments for granular materials.On the other hand,some qualitative deviations between experiments and the investigated DEM simulations are also observed,in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain,which are proven to be critical for constitutive modelling.The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective.It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

Distinct Element Modelling of Unreinforced Masonry Wall Under Seismic Loads with and without Cable Retrofitting

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.

Comparison of undrained behaviors of granular media using fluidcoupled discrete element method and constant volume method

The fluid-coupled discrete element method(DEM)and the constant volume method as two types of discrete modeling methods for fundamental study of undrained responses of granular materials,have been discussed by many researchers.The fluid-coupled DEM,which couples the motions of discrete particles with pore fluid movements,is theoretically robust although it requires a large amount of computation time.As a substitution for the complex fluid-coupled DEM,the constant volume method simulates an undrained condition for a saturated granular material by simply preserving the total volume of a granular assembly without considering interactions between fluids and particles;hence,the validity of its results is questionable.In this paper,the undrained behaviors of granular assemblies simulated using the aforementioned two methods are compared.Based on a comparison of both macroscopic and microscopic responses given by the two methods,it is demonstrated that the constant volume method may reasonably simulate the responses of a loose saturated granular material with very coarse grains,which has a high permeability,and thus a rapid pore pressure equalization.However,it is ineffective in simulating the responses of a loose material with fine components due to its failure to capture the process of a slow dissipation of the excess pore pressure among the individual pores.With regard to the dense material adopted,similar behaviors at the early and intermediate shearing stages given by the two methods are displayed.

Distinct Element Modelling of Mahabaleshwar Road Cut Hill Slope

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.

An MPI parallel DEM-IMB-LBM framework for simulating fluid-solid interaction problems

The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.