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.

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.

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.

Dynamic failure process of expanded polystyrene particle lightweight soil under cyclic loading using discrete element method

Expanded polystyrene(EPS)particle-based lightweight soil,which is a type of lightweight filler,is mainly used in road engineering.The stability of subgrades under dynamic loading is attracting increased research attention.The traditional method for studying the dynamic strength characteristics of soils is dynamic triaxial testing,and the discrete element simulation of lightweight soils under cyclic load has rarely been considered.To study the meso-mechanisms of the dynamic failure processes of EPS particle lightweight soils,a discrete element numerical model is established using the particle flow code(PFC)software.The contact force,displacement field,and velocity field of lightweight soil under different cumulative compressive strains are studied.The results show that the hysteresis curves of lightweight soil present characteristics of strain accumulation,which reflect the cyclic effects of the dynamic load.When the confining pressure increases,the contact force of the particles also increases.The confining pressure can restrain the motion of the particle system and increase the dynamic strength of the sample.When the confining pressure is held constant,an increase in compressive strain causes minimal change in the contact force between soil particles.However,the contact force between the EPS particles decreases,and their displacement direction points vertically toward the center of the sample.Under an increase in compressive strain,the velocity direction of the particle system changes from a random distribution and points vertically toward the center of the sample.When the compressive strain is 5%,the number of particles deflected in the particle velocity direction increases significantly,and the cumulative rate of deformation in the lightweight soil accelerates.Therefore,it is feasible to use 5%compressive strain as the dynamic strength standard for lightweight soil.Discrete element methods provide a new approach toward the dynamic performance evaluation of lightweight soil subgrades.

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.

Prediction of constrained modulus for granular soil using 3D discrete element method and convolutional neural networks

To efficiently predict the mechanical parameters of granular soil based on its random micro-structure,this study proposed a novel approach combining numerical simulation and machine learning algorithms.Initially,3500 simulations of one-dimensional compression tests on coarse-grained sand using the three-dimensional(3D)discrete element method(DEM)were conducted to construct a database.In this process,the positions of the particles were randomly altered,and the particle assemblages changed.Interestingly,besides confirming the influence of particle size distribution parameters,the stress-strain curves differed despite an identical gradation size statistic when the particle position varied.Subsequently,the obtained data were partitioned into training,validation,and testing datasets at a 7:2:1 ratio.To convert the DEM model into a multi-dimensional matrix that computers can recognize,the 3D DEM models were first sliced to extract multi-layer two-dimensional(2D)cross-sectional data.Redundant information was then eliminated via gray processing,and the data were stacked to form a new 3D matrix representing the granular soil’s fabric.Subsequently,utilizing the Python language and Pytorch framework,a 3D convolutional neural networks(CNNs)model was developed to establish the relationship between the constrained modulus obtained from DEM simulations and the soil’s fabric.The mean squared error(MSE)function was utilized to assess the loss value during the training process.When the learning rate(LR)fell within the range of 10-5e10-1,and the batch sizes(BSs)were 4,8,16,32,and 64,the loss value stabilized after 100 training epochs in the training and validation dataset.For BS?32 and LR?10-3,the loss reached a minimum.In the testing set,a comparative evaluation of the predicted constrained modulus from the 3D CNNs versus the simulated modulus obtained via DEM reveals a minimum mean absolute percentage error(MAPE)of 4.43%under the optimized condition,demonstrating the accuracy of this approach.Thus,by combining DEM and CNNs,the variation of soil’s mechanical characteristics related to its random fabric would be efficiently evaluated by directly tracking the particle assemblages.

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.

Revisiting rolling and sliding in two-dimensional discrete element models

It has long been recognized that the rotation of single particles plays a very important role in simulations of granular flow using the discrete element method （DEM）. Many researchers have also pointed out that the effect of rolling resistance at the contact points should be taken into account in DEM simulations. However, even for the simplest case involving two-dimensional circular particles, there is no agreement on the best way to define roiling and sliding, and different definitions and calculations of rolling and sliding have been proposed. It has even been suggested that a unique rolling and sliding definition is not possible. In this paper we assess results from previous studies on rolling and sliding in discrete element models and find that some researchers have overlooked the effect of particles of different sizes. After considering the particle radius in the derivation of rolling velocity, all results reach the same outcome： a unique solution. We also present a clear and simple derivation and validate our result using cases of rolling. Such a decomposition of relative motion is objective, or independent of the reference frame in which the relative motion is measured.

Coupled two-dimensional discrete element and multibody dynamic modeling for interaction of the soil and rough shaped rigid bodies

The study of soil and rigid body system interactions is very important for the exploration of the Moon and Mars worldwide.The discrete element method(DEM)is a relatively accurate simulation method to study dry sand soil mechanical properties.However,it is not suitable for bodies that are in mutual contact,connected due to constraints or have complex inertia properties due to their geometry.An efficient combination of the two-dimensional discrete element and multibody dynamic modeling method is proposed to solve the problem,in which the contacts and frictions among the granular spheres and the multibody system,including the smooth and rough rigid bodies,are taken into account.In this work,the soil field is modeled by a two-dimensional DEM,and the dynamics of the constrained rigid body system are modeled by the Cartesian method.A detection algorithm is developed to address the interactions between spherical discrete elements and roughly shaped rigid bodies.The advantage of this coupled method is that it enables the simultaneous capture of both responses.Finally,the program is verified by simulation experiments of the three-ball collision and the collision among the rectangular bars and the three balls.Based on this,the movement of the toothed wheel in the granular matter is analyzed,and the results show that the wheel with six teeth and 30°inclination has the fastest forward speed.In extraterrestrial objects,the wheel grip worsens,but the forward speed first increases and then decreases with decreasing gravity acceleration and loads on wheels,which proves that the coupled two-dimensional DEM and multibody dynamic program is effective in solving engineering problems.

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.

Earth Pressure of the Trapdoor Problem Using Three-Dimensional Discrete Element Method

Load transformation from the yielding part of the soil to the adjacent part is known as the soil arching effect,which plays an important role in the design of various geotechnical infrastructures.Terzaghi’s trapdoor test was an importantmilestone in the development of theories on soil arching.The research on earth pressure of the trapdoor problem is presented in this paper using the three-dimensional(3D)discrete element method(DEM).Five 3D trapdoor models with different heights are established by 3DDEMsoftware PFC 3D.The variation of earth pressure on the trapdoor with the downward movement of the trapdoor,the distribution of vertical earth pressure along the horizontal direction,the distribution of vertical earth pressure along the vertical direction,the distribution of lateral earth pressure coefficient along the depth direction,the magnitude and direction of contact force chain are studied,respectively.Related research results show that the earth pressure on the trapdoor decreases rapidly after the downward movement of the trapdoor,and then reaches the minimum earth pressure.After that,the earth’s pressure will rise slightly,and whether this phenomenon occurs depends on the depth ratio.For the bottom soil,due to the stress transfer caused by the soil arching effect,the ratio of earth pressure in the loose area decreases,while the ratio of earth pressure in the stable area increases.With the trapdoor moving down,the vertical earth pressure along the depth in the stable zone is basically consistent with the initial state,which shows an approximate linear distribution.After the trapdoor moves down,the distribution of earth pressure along with the depth in the loose area changes,which is far less than the theoretical value of vertical earth pressure of its self-weight.Because of the compression of the soil on both sides,the lateral earth pressure coefficient of most areas on the central axis of the loose zone is close to the passive earth pressure coefficient Kp.The existence of a‘soil arch’can be observed intuitively from the distribution diagram of the contact force chain in the loose zone.

A Study on the Physical Properties of Banana Straw Based on the Discrete Element Method

To improve the application of discrete element models(DEM)to the design of agricultural crushers,in this study a new highly accurate model is elaborated.The model takes into account the fiber structure,porous nature of the material and the leaf sheath coating structure.Dedicated experimental tests are conducted to determine the required“intrinsic”and basic contact parameters of the considered banana straw materials.A large number of bonding parameters are examined in relation to the particle aggregation model in order to characterize different actual banana straws.Using the particle surface energy contact model,the viscosity characteristics of the crushed material are determined together with the related stacking angle(considered as the main response factor).Through single factor experiment analysis,it is found that when the surface energy is 0.9 J·m-2,the relative error between simulations and physical experiments is 5.288%.

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.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

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.

Discrete gap breathers in a two-dimensional diatomic face-centered square lattice

In this paper we study the existence and stability of two-dimensional discrete gap breathers in a two-dimensional diatomic face-centered square lattice consisting of alternating light and heavy atoms, with on-site potential and coupling potential. This study is focused on two-dimensional breathers with their frequency in the gap that separates the acoustic and optical bands of the phonon spectrum. We demonstrate the possibility of the existence of two-dimensional gap breathers by using a numerical method. Six types of two-dimensional gap breathers are obtained, i.e., symmetric, mirror-symmetric and asymmetric, whether the center of the breather is on a light or a heavy atom. The difference between one-dimensional discrete gap breathers and two-dimensional discrete gap breathers is also discussed. We use Aubry＇s theory to analyze the stability of discrete gap breathers in the two-dimensional diatomic face-centered square lattice.

Image encryption technique based on new two-dimensional fractional-order discrete chaotic map and Menezes–Vanstone elliptic curve cryptosystem

We propose a new fractional two-dimensional triangle function combination discrete chaotic map(2D-TFCDM)with the discrete fractional difference.Moreover,the chaos behaviors of the proposed map are observed and the bifurcation diagrams,the largest Lyapunov exponent plot,and the phase portraits are derived,respectively.Finally,with the secret keys generated by Menezes-Vanstone elliptic curve cryptosystem,we apply the discrete fractional map into color image encryption.After that,the image encryption algorithm is analyzed in four aspects and the result indicates that the proposed algorithm is more superior than the other algorithms.