THE PARTIAL PROJECTION METHOD IN THE FINITE ELEMENT DISCRETIZATION OF THE REISSNER-MINDLIN PLATE MODEL

In the paper a linear combination of both the standard mixed formulation and the displacement one of the Reissner-Mindlin plate theory is used to enhance stability of the former and to remove ''locking'' of the later. For this new stabilized formulation, a unified approach to convergence analysis is presented for a wide spectrum of finite element spaces. As long as the rotation space is appropriately enriched, the formulation is convergent for the finite element spaces of sufficiently high order. Optimal-order error estimates with constants independent of the plate thickness are proved for the various lower order methods of this kind.

A Dual-Horizon Nonlocal Diffusion Model and Its Finite Element Discretization

In this paper,we present a dual-horizon nonlocal diffusion model,in which the influence area at each point consists of a standard sphere horizon and an irregular dual horizon and its geometry is determined by the distribution of the varying horizon parameter.We prove the mass conservation and maximum principle of the proposed nonlocal model,and establish its well-posedness and convergence to the classical diffusion model.Noticing that the dual horizon-related term in fact vanishes in the corresponding variational formof themodel,we then propose a finite element discretization for its numerical solution,which avoids the difficulty of accurate calculations of integrals on irregular intersection regions between the mesh elements and the dual horizons.Various numerical experiments in two and three dimensions are also performed to illustrate the usage of the variable horizon and demonstrate the effectiveness of the numerical scheme.

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.

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.

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.

Dynamic Characteristics of Irregular Ice Floes Based on Polyhedral Discrete Element Method

In polar regions, floating ice exhibits distinct characteristics across a range of spatial scales. It is well recognized that the irregular geometry of these ice formations markedly influences their dynamic behavior. This study introduces a polyhedral Discrete Element Method (DEM) tailored for polar ice, incorporating the Gilbert-Johnson-Keerthi (GJK) and Expanding Polytope Algorithm (EPA) for contact detection. This approach facilitates the simulation of the drift and collision processes of floating ice, effectively capturing its freezing and fragmentation. Subsequently, the stability and reli ability of this model are validated by uniaxial compression on level ice fields, focusing specifically on the influence of compression strength on deformation resistance. Additionally, clusters of ice floes nav igating through narrow channels are simulated. These studies have qualitatively assessed the effects of Floe Size Distribution (FSD), initial concentration, and circularity on their flow dynamics. The higher power-law exponent values in the FSD, increased circularity, and decreased concentration are each as sociated with accelerated flow in ice floe fields. The simulation results distinctly demonstrate the con siderable impact of sea ice geometry on the movement of clusters, offering valuable insights into the complexities of polar ice dynamics.

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.

True triaxial behavior of sandy soils under both drained and undrained conditions:A discrete element perspective

An advanced discrete element servomechanism that can simultaneously and independently control the evolution equations of six stress and strain components without introducing severe stress concentration is implemented.Such a comprehensive series of discrete element method simulations of both drained and undrained behavior of transversely isotropic sandy soils are successfully conducted in the true triaxial setting.During the simulation process,the evolution patterns of the load-bearing structure of the granular specimen are tracked using a contact-normal-based fabric tensor.The simulation results show that sandy soils exhibit more significant non-coaxiality between the loading direction and the major principal direction of the fabric tensor under extension than under compression.Therefore,the fabric of the sandy soils under extension has a stronger tendency to evolve toward the loading direction than that under compression,causing a more significant disturbance to the load-bearing structure.Consequently,compared with the extension loading condition,the transversely isotropic specimen under compression exhibits a higher shear strength and stronger dilatancy under drained conditions and a stronger liquefaction resistance under undrained conditions.

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.

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.

Discrete element and finite element coupling simulation and experiment of hot granule medium pressure forming

The granule medium of discreteness is supposed to be continuous（Drucker-Prager model） in the existing finite element simulation analysis on the hot granule medium pressure forming（HGMF） process, so the granule medium may produce tensile stress in the process of pressure-transferring and flowing, which does not coincide with the reality. The analysis method, discrete element and finite element（DE-FE） coupling simulation, is proposed to solve the problem. The material parameters of simulation model are obtained by the pressure-transfer performance test of granule medium and the hot uniaxial tensile test of sheet metal. The DE-FE coupling simulation platform is established by adopting Visual Basic language. The features in the process that AA7075-T6 conical parts are formed by the HGMF process are analyzed and verified by the process test. The studies show that the results of DE-FE coupling simulation coincide well with the test results, which provides a new analysis method to solve the mechanics problem in the coupling of discrete and continuum.

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.

Discrete element modeling of migration and evolution rules of coarse aggregates in the static compaction process

To investigate migration and evolution rules of coarse aggregates in the static compaction process, an algorithm of generating digital coarse aggregates that can reflect real morphology（ such as shape, size and fracture surface） of aggregate particles, is represented by polyhedral particles based on the discrete element method（ DEM）. A digital specimen comprised of aggregates and air voids is developed. In addition,a static compaction model consisting of a digital specimen and three plates is constructed and a series of evaluation indices such as mean contact force σMCF, wall stress in direction of zcoordinate σWSZZ, porosity and coordination numbers are presented to investigate the motion rules of coarse aggregates at different compaction displacements of 7. 5, 15 and 30 mm. The three-dimensional static compaction model is also verified with laboratory measurements. The results indicate that the compaction displacements are positively related to σMCF and σWSZZ, which increase gradually with the increase in iterative steps. When the compaction proceeds, the digital specimen porosity decreases, but the coordination number increases. The variation ranges of these four indices are different at different compaction displacements. This study provides a method to analyze the compaction mechanism of particle materials such as asphalt mixture and graded broken stone.

Discrete element modeling of ice loads on ship hulls in broken ice fields

Ice loads on a ship hull affect the safety of the hull structure and the ship maneuvering performance in ice-covered regions. A discrete element method （DEM） is used to simulate the interaction between drifting ice floes and a moving ship. The pancake ice floes are modelled with three-dimensional （3-D） dilated disk elements considering the buoyancy, drag force and additional mass induced by the current. The ship hull is modelled with 3D disks with overlaps. Ice loads on the ship hull are determined through the contact detection between ice floe element and ship hull element and the contact force calculation. The influences of different ice conditions （current velocities and directions, ice thicknesses, concentrations and ice floe sizes） and ship speeds are also examined on the dynamic ice force. The simulated results are compared qualitatively well with the existing field data and other numerical results. This work can be helpful in the shil3 structure design and the navigation securitv in ice-covered fields.

Discrete Element Modeling of Debris Avalanche Impact on Retaining Walls

In China,gravity retaining walls are widely used as protection structures against rockfalls,debris flows and debris avalanches along the roads in mountainous areas.In this paper,the Discrete Element Method(DEM) has been used to investigate the impact of granular avalanches and debris flows on retaining walls.The debris is modeled as two dimensional circular disks that interact through frictional sliding contacts.The basic equations that control the deformation and motion of the particles are introduced.A series of numerical experiments were conducted on an idealized debris slide impacting a retaining wall.The parametric study has been performed to examine the influences of slope geometry,travel distance of the sliding mass,wall position,and surface friction on the impact force exerted on the wall.Results show that:1) the force achieves its maximum value when slope angle is equal to 60°,as it varies from 30° to 75°;2) an approximate linear relationship between the impact force and the storage area length is determined.

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.

Evaluating the mechanical properties of anisotropic shale containing bedding and natural fractures with discrete element modeling

Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.