A study of influence of gravity on bulk behavior of particulate solid

This paper examines the influence of gravity on the bulk responses of a granular solid. The loading scenarios in this study include confined compression, rod penetration into a granular medium and discharging through an orifice. Similar loading and flow conditions are likely to be encountered in the stress and deformation regimes that regoliths are subjected to in extraterrestrial exploration activities including in situ resource utilisation processes. Both spherical and non-spherical particles were studied using the discrete element method (DEM). Whilst DEM is increasingly used to model granular solids, careful validations of the simulation outcomes are rather rare. Thus in addition to exploring the effect of gravity, this paper also compares DEM simulations with experiments under terrestrial condition to verify whether DEM can produce satisfactory predictions.The terrestrial experiments were conducted with great care and simulated closely using DEM. The key mechanical and geometrical properties for the particles were measured in laboratory tests for use in the DEM simulations. A series of DEM computations were then performed under reduced gravity to simulate these experiments under extraterrestrial environment. It was found that gravity has no noticeable effect on the force transmission in the confined compression case; the loading gradient in the rod penetration is linearly proportional to the gravity; the mass flow rate in silo discharge is proportional to square root of the gravity and the angle of repose increases with reducing gravity. These findings are in agreement with expectation and existing scientific evidence.     

Three-dimensional discrete element models for the granular statics and dynamics of powders in cavity filling

Rapid granular flow from a moving container and angle of repose formation were investigated by numerical simulations using the discrete element method (DEM) and experiments. Grain models of various geometrical complexity were studied and their ability to reproduce the experiments in those regimes was explored. The predictive power of the most realistic model for gravity driven cavity filling was assessed. Good agreement between computed and measured density distributions within the filled cavities provides a basis for numerical process variations aiming at homogenized density distributions. The effect of numerical coarse graining was found to be negligible for all properties of interest provided that force laws are scaled properly and corrections for boundary effects are taken into account. The proposed scaling was tested for a certain set of force laws but could be applied to different DEM forces as well. An analytic mass flow law for powder discharge from a moving container was derived and verified by our DEM simulations.     

Discrete element parameter calibration and the modelling of dragline bucket filling

The Discrete Element Method (DEM) is useful for modelling granular flow. The accuracy of DEM modelling is dependent upon the model parameter values used. Determining these values remains one of the main challenges. In this study a method for determining the parameters of cohesionless granular material is presented. The particle size and density were directly measured and modelled. The particle shapes were modelled using two to four spheres clumped together. The remaining unknown parameter values were determined using confined compression tests and angle of repose tests. This was done by conducting laboratory experiments followed by equivalent numerical experiments and iteratively changing the parameters until the laboratory results were replicated. The modelling results of the confined compression tests were mainly influenced by the particle stiffness. The modelling results of the angle of repose tests were dependent on both the particle stiffness and the particle friction coefficient. From these observations, the confined compression test could be used to determine the particle stiffness and with the stiffness known, the angle of repose test could be used to determine the particle friction coefficient. Usually DEM codes do not solve the equations of motion for so-called walls (non-granular structural elements). However, in this study a dynamic model of a dragline bucket is developed and implemented in a commercial DEM code which allows the dynamics of the walls to be modelled. The DEM modelling of large systems of particles is still a challenge and procedures to simplify and speed up the modelling of dragline bucket filling are presented. Using the calibrated parameters, numerical results of bucket filling are compared to experimental results. The model accurately predicted the orientation of the bucket. The model also accurately predicted the drag force over the first third of the drag, but predicted drag forces too high for the subsequent part of the drag.     

Experimental and DEM analyses on wheel-soil interaction

In this paper, the wheel-soil interaction for a future lunar exploration mission is investigated by physical model tests and numerical simulations. Firstly, a series of physical model tests was conducted using the TJ-1 lunar soil simulant with various driving conditions, wheel configurations and ground void ratios. Then the corresponding numerical simulations were performed in a terrestrial environment using the Distinct Element Method (DEM) with a new contact model for lunar soil, where the rolling resistance and van der Waals force were implemented. In addition, DEM simulations in an extraterrestrial (lunar) environment were performed. The results indicate that tractive efficiency does not depend on wheel rotational velocity, but decreases with increasing extra vertical load on the wheel and ground void ratio. Rover performance improves when wheels are equipped with lugs. The DEM simulations in terrestrial environment can qualitatively reproduce the soil deformation pattern as observed in the physical model tests. The variations of traction efficiency against the driving condition, wheel configuration and ground void ratio attained in the DEM simulations match the experimental observations qualitatively. Moreover, the wheel track is found to be less evident and the tractive efficiency is higher in the extraterrestrial environment compared to the performance on Earth.     

不同重力场下颗粒冲击过程的离散元分析

开展不同重力场下颗粒材料冲击动力学研究有助于加深对颗粒运动机制的理解和深空探测任务的实施。本文采用离散元模拟对颗粒介质受球形冲击物的冲击过程进行了数值模拟，并与地球重力下冲击的试验结果进行对比验证。在此基础上，进一步研究了重力加速度对冲击物动力学的影响规律。计算结果表明，在所有重力加速度下，冲击物的穿透深度d与冲击速度v₀的关系可以用Poncelet模型表达；d与冲击物下落的总高度H表现为d～Hⁿ的幂律关系，当H<10 m时，d与H的幂率标度为0.322,而H>10 m时，d与H的幂率标度下降到0.211。此外，穿透深度小于冲击物半径时，重力加速度对冲击物减速过程无影响。在所有的重力加速度下，当冲击速度大于5 m/s时，冲击物的持续碰撞时间t_c是恒定的，且与重力的-1/2次方呈正比。     

Determining the angle of repose of sand under low-gravity conditions using discrete element method

This study is a comparative investigation of data, collected through experimental and numerical means, related to the flow of sand particles through a hopper under low-gravity conditions. During a parabolic airplane flight simulating low-gravity conditions, we determined effects of gravity on the angle of repose of sand pile particles by flowing dry sand from a hopper. The gravity effects on the angle of repose of the sand were negligible. Two-dimensional discrete element method (DEM) was used to simulate the angle of repose. Results were compared to observations made during the low-gravity experiments. Effects of varying parameters such as the friction coefficient and coefficient of rolling friction were determined by running various DEM simulations. Moreover, the effect of the elemental radius on the angle of repose was investigated using DEM. The angle of repose is influenced by certain changes in the friction coefficient and rolling friction values, but the elemental radius has only a negligible effect on the angle of repose within the range of variation. Results show that the DEM model used for this study might be applicable to determine terramechanical interactions under lunar surface gravity conditions, provided that parameters are adjusted and an extended period of simulation is achieved.     

Discrete element modeling of inherently anisotropic granular assemblies with polygonal particles

In the present article, we study the effect of inherent anisotropy, i.e., initial bedding angle of particles and associated voids on macroscopic mechanical behavior of granular materials, by numerical simulation of several biaxial compression tests using the discrete element method (DEM). Particle shape is considered to be irregular convex-polygonal. The effect of inherent anisotropy is investigated by following the evolution of mobilized shear strength and volume change during loading. As experimental tests have already shown, numerical simulations also indicate that initial anisotropic condition has a great influence on the strength and deformational behavior of granular assemblies. Comparison of simulations with tests using oval particles, shows that angularity influences both the mobilized shear strength and the volume change regime, which originates from the interlocking resistance between particles.     

Segregation of particulate solids: Experiments and DEM simulations

Segregation of granular materials is a complex phenomenon, difficult to measure quantitatively and to predict. Discrete element method (DEM) can be a useful tool to predict segregation effects and to support the industrial design. In this context, a very challenging idea is the characterization of the granular solids to provide the key parameters needed for a successful DEM simulation of segregation processes. Rolling friction, sliding friction and the coefficient of restitution are the critical parameters to be studied. These microscopic simulation parameters are calibrated by comparing the macroscopic behavior of granular matter in standard bulk experiments, which have the advantage of being highly repeatable and reliable.An experimental method is presented to characterize free surface segregation. The effects of different particle properties, particularly, shape and size, on segregation of cohesionless materials were investigated. From the experiments, particle size demonstrated a stronger effect on segregation than particle shape. Finally, the corresponding DEM simulations of the segregation experiments were presented. The parameters obtained by calibration were validated by the comparison of the modeled segregation behavior with the experimental results. Thus, calibrated DEM simulations are capable of predicting segregation effects.     

单轴压缩下页岩类脆性材料细观损伤演化的离散元分析

了解广泛存在的类似页岩的脆性材料各向异性对工程安全具有重要意义。本研究将页岩视为粘结颗粒材料,基于离散单元方法研究了横观各向同性脆性页岩的损伤演化。再现了不同层理角的页岩试样的破坏模式,并对比了实验和数值模拟的抗压强度和弹性模量。引入微裂纹的概念,通过定义裂纹密度函数,系统地研究了单轴压缩条件下,页岩层理角对细观结构的影响。此外基于平均配位数建立了配位数变化与细观损伤的联系,并根据配位数的变化与裂纹数量将加载过程分为三个阶段,分析了不同阶段配位数与裂纹数量的对应变化关系。研究表明,页岩的裂纹密度随着层理角的增加而增加,而试样的平均配位数在加载过程中先上升后剧烈下降,颗粒集合体在单轴压缩条件下的应力应变及裂纹数量曲线与平均配位数曲线有良好的一致性。该研究揭示了横观各向同性脆性岩石的破坏过程和内在机理,将为页岩类脆性材料的工程应用提供理论指导。     

Segregation of particulate solids＂ Experiments and DEM simulations

Segregation of granular materials is a complex phenomenon, difficult to measure quantitatively and to predict. Discrete element method （DEM） can be a useful tool to predict segregation effects and to support the industrial design. In this context, a very challenging idea is the characterization of the granular solids to provide the key parameters needed for a successful DEM simulation of segregation processes. Rolling friction, sliding friction and the coefficient of restitution are the critical parameters to be studied. These microscopic simulation parameters are calibrated by comparing the macroscopic behavior of granular matter in standard bulk experiments, which have the advantage of being highly repeatable and reliable. An experimental method is presented to characterize free surface segregation. The effects of different particle properties, particularly, shape and size, on segregation of cohesionless materials were investi- gated. From the experiments, particle size demonstrated a stronger effect on segregation than particle shape. Finally, the corresponding DEM simulations of the segregation experiments were presented. The parameters obtained by calibration were validated by the comparison of the modeled segregation behav- ior with the experimental results. Thus, calibrated DEM simulations are capable of predicting segregation effects.     

Numerical simulation study on particle breakage behavior of granular materials in confined compression tests

In order to study the fragmentation law, the confined compression experiment of granular assemblies has been conducted to explore the particle breakage characteristic by DEM approach in this work. It is shown that contact and contact force during the loading process gradually transform from anisotropy to isotropy. Meanwhile, two particle failure modes caused by different contact force states are analyzed, which are single-through-crack failure and multi-short-crack failure. Considering the vertical distribution of the number of cracks and the four characteristic stress distributions (the stress related to the maximum contact force, the major principal stress, the deviatoric stress and the mean stress), it is pointed out that the stress based on the maximum contact force and the major principal stress can reflect the distribution of cracks accurately. In addition, the size effect of particle crushing indicates that small size particles are prone to break. The lateral pressure coefficient of four size particles during the loading process is analyzed to explain the reason for the size effect of particle breakage.     

Effect of particle size distribution on micro- and macromechanical response of granular packings under compression

The role of particle size heterogeneity on micro- and macromechanical properties of assemblies of spherical particles was studied using DEM simulations. The response to an imposed load of a granular material composed of non-uniformly sized spheres subjected to uniaxial confined compression was investigated. A range of geometrical and micro-mechanical properties of granular packings (e.g., void fraction, contact force distribution, average coordination number and degree of mobilisation of friction at contacts between particles) were examined, and provided a more accurate interpretation of the macroscopic behaviour of mixtures than has previously been available. The macromechanical study included stress transmission, stiffness and angle of internal friction of the granular assemblies.The degree of polydispersity showed slight effect on both, the void fraction and the elastic properties of the system. The tendency for increase in the lateral-to-vertical pressure ratios was observed with an increasing degree of particle size heterogeneity; however, the different pressure ratios calculated for samples with various degrees of polydispersity lay within the range of data scatter.     

Numerical simulations of shock-induced load transfer processes in granular media using the discrete element method

By the discrete element method (DEM), we perform numerical simulations of shock-induced load transfer processes in granular layers composed of spherical particles packed in vertical channels. In order to isolate the load transfer through the grains’ contact points from the complicated load transfer processes, we simulate the shock wave interactions with granular layers having no permeability for gas. The shock loading is achieved by applying a downward step force on the top of the granular layers. Complex, three-dimensional load transfer processes in the granular media, which are extremely difficult to understand from experiments, are visualized based on the results from the present DEM simulation. The numerical results show that highly concentrated load transfer paths, through which shock loads are transferred mainly, exist in the granular media, and that the dimensions of the container of the granular media considerably affect the shock-induced load transfer processes. From a coarse-grained representation of intergranular stress, wave-like load transfer processes are clearly observed. For relatively deep granular layers, however, the wave fronts became unclear as they propagated.     

Hf基非晶合金夹芯结构长杆弹的侵彻行为

为研究Hf基非晶合金的变形行为及高速侵彻性能,分别开展了Hf基非晶合金材料静动态力学性能和Hf基非晶合金夹芯结构长杆弹高速侵彻45钢靶体试验研究,并与45钢夹芯长杆弹侵彻结果进行对比。研究发现:Hf基非晶合金具有较高的断裂强度,断裂时伴随有能量释放现象;Hf基非晶合金夹芯长杆弹侵彻钢靶过程可分为3个阶段:开坑、夹芯结构侵彻和剩余弹体侵彻。Hf非晶合金在侵彻过程中发生了明显的释能反应,显著地增强了弹体毁伤效应,扩大了侵彻弹孔直径,增加了弹体侵彻深度和弹孔体积。在高速冲击下,Hf基非晶合金夹芯长杆弹表现出优异的侵彻性能,可以为非晶合金材料在高效毁伤领域的应用提供新思路。     

ROTATIONAL RESISTANCE AND SHEAR-INDUCED ANISOTROPY IN GRANULAR MEDIA

This paper presents a micromechanical study on the behavior of granular materials under confined shear using a three-dimensional Discrete Element Method （DEM）. We consider rotational resistance among spherical particles in the DEM code as an approximate way to account for the effect of particle shape. Under undrained shear, it is found rotational resistance may help to increase the shear strength of a granular system and to enhance its resistance to liquefaction. The evolution of internal structure and anisotropy in granular systems with different initial conditions depict a clear bimodal character which distinguishes two contact subnetworks. In the presence of rotational resistance, a good correlation is found between an analytical stress-force-fabric relation and the DEM results, in which the normal force anisotropy plays a dominant role. The unique properties of critical state and liquefaction state in relation to granular anisotropy are also explored and discussed.     

Micro-scale behavior of granular materials during cyclic loading

This study presents the micro-scale behavior of granular materials under biaxial cyclic loading for differ- ent confining pressures using the two-dimensional （2D） discrete element method （DEM）. Initially, 8450 ovals were generated in a rectangular frame without any overlap. Four dense samples having confining pressures of 15, 25, 50, and 100 kPa were prepared from the initially generated sparse sample. Numeri- cal simulations were performed under biaxial cyclic loading using these isotropically compressed dense samples. The numerical results depict stress-strain-dilatancy behavior that was similar to that observed in experimental studies. The relationship between the stress ratio and dilatancy rate is almost indepen- dent of confining pressures during loading but significantly dependent on the confining pressures during unloading. The evolution of the coordination number, effective coordination number and slip coordina- tion number depends on both the confining pressures and cyclic loading. The cyclic loading significantly affects the microtopology of the granular assembly. The contact fabric and the fabric-related anisotropy are reported, as well. A strong correlation between the stress ratio and the fabric related to contact normals is observed during cyclic loading, irrespective of confining pressures.     

Calibration of granular material parameters for DEM modelling and numerical verification by blade-granular material interaction

The discrete element method (DEM) is a promising approach to model blade-granular material interactions. The accuracy of DEM models depends on the model parameters. In this study, a calibration process was developed to determine the parameter values. The particle size was the same as the real material and the particle shape was modelled using two spherical particles rigidly clumped together to form a single grain. Laboratory shear tests and compressions tests were used to determine the material internal friction angle and stiffness, respectively. These tests were replicated numerically using DEM models with different sets of particle friction coefficients and particle stiffness values. The shear test results are found to be dependent on both the particle friction coefficient and the particle stiffness. The compression test results show that it is only dependent on the particle stiffness. The combination of shear test and compression test results can be used to determine a unique set of particle friction and particle stiffness values. The calibration process was validated experimentally and numerically by modelling a blade moving through granular material. Results show that the forces acting on the blade can be accurately modelled with DEM and the maximum error is found to be 26%. The relative particle-blade displacements were used to predict the position and shape of the shear lines in front of the blade. A good qualitative correlation was achieved between the experiments and the DEM simulations.     

Micro–macro analysis of granular material behavior along proportional strain paths

When granular materials are subjected to proportional strain loading paths, they manifest a variety of behaviors depending on the initial void ratio of the specimen as well as the imposed dilatancy/contractancy rate. In some cases, the stress components may vanish over the duration of the test, and the specimen may progressively liquefy. To investigate this behavior, the authors have developed a kinematic approach to be deployed in two parts. First, numerical simulations are performed by means of a discrete element method. Secondly, two micromechanical models have corroborated the DEM results. The performance of these models may explain a number of microstructural mechanisms responsible for the macroscopic constitutive behavior.