Damage Mechanism of Ultra-thin Asphalt Overlay(UTAO) based on Discrete Element Method

Aiming to analyze the damage mechanism of UTAO from the perspective of meso-mechanical mechanism using discrete element method(DEM),we conducted study of diseases problems of UTAO in several provinces in China,and found that aggregate spalling was one of the main disease types of UTAO.A discrete element model of UTAO pavement structure was constructed to explore the meso-mechanical mechanism of UTAO damage under the influence of layer thickness,gradation,and bonding modulus.The experimental results show that,as the thickness of UTAO decreasing,the maximum value and the mean value of the contact force between all aggregate particles gradually increase,which leads to aggregates more prone to spalling.Compared with OGFC-5 UTAO,AC-5 UTAO presents smaller maximum and average values of all contact forces,and the loading pressure in AC-5 UTAO is fully diffused in the lateral direction.In addition,the increment of pavement modulus strengthens the overall force of aggregate particles inside UTAO,resulting in aggregate particles peeling off more easily.The increase of bonding modulus changes the position where the maximum value of the tangential force appears,whereas has no effect on the normal force.

An Innovative Coupled Common-Node Discrete Element Method-Smoothed Particle Hydrodynamics Model Developed with LS-DYNA and Its Applications

In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.

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.

Investigation of the block toppling evolution of a layered model slope by centrifuge test and discrete element modeling

Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.

Discrete Element Modelling of Damage Evolution of Concrete Considering Meso-Structure of ITZ

The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.

Three dimensional discrete element modelling of the conventional compression behavior of gas hydrate bearing coal

To analyze the relationship between macro and meso parameters of the gas hydrate bearing coal(GHBC)and to calibrate the meso-parameters,the numerical tests were conducted to simulate the laboratory triaxial compression tests by PFC3D,with the parallel bond model employed as the particle contact constitutive model.First,twenty simulation tests were conducted to quantify the relationship between the macro–meso parameters.Then,nine orthogonal simulation tests were performed using four meso-mechanical parameters in a three-level to evaluate the sensitivity of the meso-mechanical parameters.Furthermore,the calibration method of the meso-parameters were then proposed.Finally,the contact force chain,the contact force and the contact number were examined to investigate the saturation effect on the meso-mechanical behavior of GHBC.The results show that:(1)The elastic modulus linearly increases with the bonding stiffness ratio and the friction coefficient while exponentially increasing with the normal bonding strength and the bonding radius coefficient.The failure strength increases exponentially with the increase of the friction coefficient,the normal bonding strength and the bonding radius coefficient,and remains constant with the increase of bond stiffness ratio;(2)The friction coefficient and the bond radius coefficient are most sensitive to the elastic modulus and the failure strength;(3)The number of the force chains,the contact force,and the bond strength between particles will increase with the increase of the hydrate saturation,which leads to the larger failure strength.

Correlations between mineral composition and mechanical properties of granite using digital image processing and discrete element method

This study investigated the correlations between mechanical properties and mineralogy of granite using the digital image processing(DIP) and discrete element method(DEM). The results showed that the X-ray diffraction(XRD)-based DIP method effectively analyzed the mineral composition contents and spatial distributions of granite. During the particle flow code(PFC2D) model calibration phase, the numerical simulation exhibited that the uniaxial compressive strength(UCS) value, elastic modulus(E), and failure pattern of the granite specimen in the UCS test were comparable to the experiment. By establishing 351 sets of numerical models and exploring the impacts of mineral composition on the mechanical properties of granite, it indicated that there was no negative correlation between quartz and feldspar for UCS, tensile strength(σ_(t)), and E. In contrast, mica had a significant negative correlation for UCS, σ_(t), and E. The presence of quartz increased the brittleness of granite, whereas the presence of mica and feldspar increased its ductility in UCS and direct tensile strength(DTS) tests. Varying contents of major mineral compositions in granite showed minor influence on the number of cracks in both UCS and DTS tests.

Computational fluid dynamics-discrete element method simulation of stirred tank reactor for graphene production

Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.

Capability of discrete element method to investigate the macro-micro mechanical behaviours of granular soils considering different stress conditions and morphological gene mutation

Discrete element method(DEM)has been widely utilised to model the mechanical behaviours of granular materials.However,with simplified particle morphology or rheology-based rolling resistance models,DEM failed to describe some responses,such as the particle kinematics at the grain-scale and the principal stress ratio against axial strain at the macro-scale.This paper adopts a computed tomography(CT)-based DEM technique,including particle morphology data acquisition from micro-CT(mCT),spherical harmonic-based principal component analysis(SH-PCA)-based particle morphology reconstruction and DEM simulations,to investigate the capability of DEM with realistic particle morphology for modelling granular soils’micro-macro mechanical responses with a consideration of the initial packing state,the morphological gene mutation degree,and the confining stress condition.It is found that DEM with realistic particle morphology can reasonably reproduce granular materials’micro-macro mechanical behaviours,including the deviatoric stressevolumetric straineaxial strain response,critical state behaviour,particle kinematics,and shear band evolution.Meanwhile,the role of multiscale particle morphology in granular soils depends on the initial packing state and the confining stress condition.For the same granular soils,rougher particle surfaces with a denser initial packing state and a higher confining stress condition result in a higher degree of shear strain localisation.

Effect of heterogeneity on mechanical and micro-seismic behaviors of sandstone subjected to multi-level cyclic loading: A discrete element method investigation

In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.This will lead to a poor reproduction of rock’s behaviors in terms of stress-strain relationship and micro-seismic characteristics in numerical simulation.This work aims to analyze and reveal the impact of parameter heterogeneity on the rock’s fatigue and micro-seismic properties based on PFC3D.Two distribution patterns(uniform and Weibull distributions),are implemented to assign four critical parameters(i.e.tensile strength,cohesion,parallel bond stiffness and linear stiffness)for 32 sets of numerical schemes.The results show that the models with high heterogeneity of tensile strength and cohesion can better reproduce the stress-strain relationship as well as the patterns of cumulative AE counts and energy magnitude.The evolution of the proportion of three-level AE events in the laboratory test is consistent with the numerical results when the highly heterogeneous tensile strength and cohesion are distributed.The numerical results can provide practical guidance to the PFC-based modeling of rock heterogeneity when exposed to multi-level cyclic loading and AE monitoring.

Virtual rutting test of asphalt mixture using discrete element method

In order to investigate the permanent deformation behavior of asphalt mixtures from discontinuity, the virtual rutting test of asphalt mixtures is developed by the discrete element method （DEM）. A digital specimen generation procedure considering aggregate gradation and irregular shape is developed based on the probability theory and the Monte Carlo method. The virtual rutting test is then conducted based on the generated digital specimen. In addition, on the basis of the time-temperature superposition （TTS） principle, a calculation method is used to reduce the computation time of the virtual rutting test. The simulation results are compared with the laboratory measurements. The results show that the calculation method based on the TFS principle in the discrete element （DE） viscoelastic model can significantly reduce the computation time. The deformation law of asphalt mixtures in the virtual rutting test is similar to the laboratory measurements, and the deformation and the dynamic stability of the virtual rutting test are slightly greater than the laboratory measurements. The two-dimensional virtual rutting test can predict the permanent deformation performance of asphalt mixtures.

Virtual fracture test of asphalt mixture based on discrete element method

In order to study fracture behaviors of asphalt mixtures, virtual tests of the two-dimensional（2D） microstructure based on the discrete element method（ DEM） are designed. The virtual structure of the 2D digital specimen of asphalt mixture is generated based on a particle generation program, in which the gradation and the irregular shapes of aggregates are considered. With the 2D digital specimens, a DEM-based mixture model is established and center-point beam fracture simulation tests are conducted by the DEM. Meanwhile, a series of calibration tests are carried out in laboratory to evaluate the DEM model and validate the methods of virtual fracture tests. The test results indicate that the fracture intensity of asphalt mixtures predicted by the DEM matches very well with the intensity obtained in laboratory. It is concluded that the microstructural virtual tests can be used as a supplemental tool to evaluate fracture properties of asphalt mixtures.

A Consistent Time Level Implementation Preserving Second-Order Time Accuracy via a Framework of Unified Time Integrators in the Discrete Element Approach

In this work,a consistent and physically accurate implementation of the general framework of unified second-order time accurate integrators via the well-known GSSSS framework in the Discrete Element Method is presented.The improved tangential displacement evaluation in the present implementation of the discrete element method has been derived and implemented to preserve the consistency of the correct time level evaluation during the time integration process in calculating the algorithmic tangential displacement.Several numerical examples have been used to validate the proposed tangential displacement evaluation;this is in contrast to past practices which only seem to attain the first-order time accuracy due to inconsistent time level implementation with different algorithms for normal and tangential directions.The comparisons with the existing implementation and the superiority of the proposed implementation are given in terms of the convergence rate with improved numerical accuracy in time.Moreover,several schemes via the unified second-order time integrators within the framework of the GSSSS family have been carried out based on the proposed correct implementation.All the numerical results demonstrate that using the existing state-of-the-art implementation reduces the time accuracy to be first-order accurate in time,while the proposed implementation preserves the correct time accuracy to yield second-order.

Assessment of strain bursting in deep tunnelling by using the finite-discrete element method

Rockbursting in deep tunnelling is a complex phenomenon posing significant challenges both at the design and construction stages of an underground excavation within hard rock masses and under high in situ stresses. While local experience, field monitoring, and informed data-rich analysis are some of the tools commonly used to manage the hazards and the associated risks, advanced numerical techniques based on discontinuum modelling have also shown potential in assisting in the assessment of rockbursting. In this study, the hybrid finite-discrete element method(FDEM) is employed to investigate the failure and fracturing processes, and the mechanisms of energy storage and rapid release resulting in bursting, as well as to assess its utility as part of the design process of underground excavations.Following the calibration of the numerical model to simulate a deep excavation in a hard, massive rock mass, discrete fracture network(DFN) geometries are integrated into the model in order to examine the impact of rock structure on rockbursting under high in situ stresses. The obtained analysis results not only highlight the importance of explicitly simulating pre-existing joints within the model, as they affect the mobilised failure mechanisms and the intensity of strain bursting phenomena, but also show how the employed joint network geometry, the field stress conditions, and their interaction influence the extent and depth of the excavation induced damage. Furthermore, a rigorous analysis of the mass and velocity of the ejected rock blocks and comparison of the obtained data with well-established semi-empirical approaches demonstrate the potential of the method to provide realistic estimates of the kinetic energy released during bursting for determining the energy support demand.

THE APPLICATION OF DISCRETE ELEMENT METHOD IN SOLVING THREE-DIMENTIONAL IMPACT DYNAMICS PROBLEMS

A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.

Geo-engineered buffer capacity of two-layered absorbing system under the impact of rock avalanches based on Discrete Element Method

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.

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.

Understanding roof deformation mechanics and parametric sensitivities of coal mine entries using the discrete element method

Although conventional coal mine designs are conservative regarding pillar strength,local failures such as roof-falls and pillar bursts still affect mine safety and operations.Previous studies have identified that discontinuous,layered roof materials have some self-supporting capacity.This research is a preliminary step towards understanding these mechanics in coal-measure rocks.Although others have considered broad conceptual models and simplified analogs for mine roof behavior,this study presents a unique numerical model that more completely represents in-situ roof conditions.The discrete element method(DEM)is utilized to conduct a parametric analysis considering a range of in-situ stress ratios,material properties,and joint networks to determine the parameters controlling the stability of single-entries modeled in two-dimensions.Model results are compared to empirical observations of roof-support effectiveness(ARBS)in the context of the coal mine roof rating(CMRR)system.Results such as immediate roof displacement,overall stability,and statistical relationships between model parameters and outcomes are presented herein.Potential practical applications of this line of research include:(1)roof-support optimization for a range of coal-measure rocks,(2)establishment of a relationship between roof stability and pillar stress,and(3)determination of which parameters are most critical to roof stability and therefore require concentrated evaluation.

Discrete element method of improved performance of railway ballast bed using elastic sleeper

With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.

Predicting the Dynamic Behavior of Asphalt Concrete Using Three-dimensional Discrete Element Method

A user-defined three-dimensional （3D） discrete element model was presented to predict the dynamic modulus and phase angle of asphalt concrete （AC）. The 3D discrete element method （DEM） model of AC was constructed employing a user-defined computer program developed using the ＂Fish＂ language in PFC3D. Important microstructural features of AC were modeled, including aggregate gradation, air voids and mastic. The irregular shape of aggregate particle was modeled using a clump of spheres. The developed model was validated through comparing with experimental measurements and then used to simulate the cyclic uniaxial compression test, based on which the dynamic modulus and phase angle were calculated from the output stress- strain relationship. The effects of air void content, aggregate stiffness and volumetric fraction on AC modulus were further investigated. The experimental results show that the 3D DEM model is able to accurately predict both dynamic modulus and phase angle of AC across a range of temperature and loading frequencies. The user- defined 3D model also demonstrated significant improvement over the general existing two-dimensional models.