Mesomechanical simulation of direct shear test on outwash deposits with granular discrete element method

The mechanical properties of outwash deposits which are taken as unconsolidated geo-materials with the characteristics of non-uniformity, heterogeneity and multiphase have attracted much attention in engineering. According to the results of laboratory direct shear test on the remolded samples, the soil particle parameters of numerical model based on in-situ particle size cumulative curves and 3D granular discrete element method were determined. Then, numerical experiments on different lithology, stone content and gradation composition were conducted. The results show that it is not a flat surface but a shear band that yields in the sample. The curve of particle velocity vs distance from the designed shear surface of test model that is taken as a datum plane in the vertical section of sample shows in "S" shape. The shear disturbance area is about twice the maximum diameter of stone blocks. The greater the stiffness of stone is, the rougher the shear surface is. The shear strength of outwash deposits is largely controlled by lithology and stone content, and the bite force between stone blocks is the root reason of larger friction angle. It is also shown that strain hardening and low shear dilatancy occur under high confining pressure as well as possibility of shear shrinkage. But it is easy to behave shear dilatation and strain softening under low confining pressure. The relationship between particle frictional coefficient and stone content presents an approximately quadratic parabola increase. The strain energy first increases and then drops with the increase of frictional energy. The cohesion increases with soil stiffness increasing but decreases with stone stiffness increasing. Numerical results are consistent with the laboratory test results of remolded samples, which indicate that this method can be a beneficial supplement to determine the parameters of engineering deposit bodies.

GPU-based discrete element simulation on flow stability of flat-bottomed hopper

In this study, the flow stability of the flat-bottomed hopper was investigated via GPU-based discrete element method(DEM) simulation. With the material height inside the hopper reducing, the fluctuation of the flow rate indicates an unstable discharge. The flow regions of the unstable discharge were compared with that of the stable discharge, a key transformation zone, where the voidage showed the largest difference between unstable and stable discharge, was revealed. To identify the relevance of the key transformation zone and the hopper flow stability, the voidage variation of the key transformation zone with material height reducing was studied.A sharp increase in the voidage in the key transformation zone was considered to be the standard for judging the unstable hopper flow, and the ‘Top–Bottom effect' of the hopper was defined, which indicated the hopper flow was unstable when the hopper only had the top area and the bottom area, because the voidage of particles in the top area and the bottom area were both variables.

Discrete element method study of shear-driven granular segregation in a slowly rotating horizontal drum

Segregation and mixing of granular materials are complex processes and are not fully understood. Motivated by industrial need, we performed a simulation using the discrete element method to study size segregation of a binary mixture of granular particles in a horizontal rotating drum. Particles of two dif- ferent sizes were poured into the drum until it was 50% full. Shear-driven segregation was induced by rotating the side-plates of the drum in the opposite direction to that of the cylindrical wall. We found that radial segregation diminished in these systems but did not completely vanish. In an ordinary rotating drum, a radial core of smaller particles is formed in the center of the drum, surrounded by larger revolving particles. In our system, however, the smaller particles were found to migrate toward the side-plates. The shear from anti-spinning side-plates reduces the voidage and increases the bulk density. As such, smaller particles in the mixer tend to move to denser regions. We varied the shear by changing the coefficient of friction on the side-plates to study the influence of shear rate on this migration. We also compared the extent of radial segregation with stationary side-plates and with side-plates moving in different angular directions.

A discrete element model for simulating saturated granular soil

A numerical model is developed to simulate saturated granular soil, based on the discrete element method. Soil particles are represented by Lagrangian discrete elements, and pore fluid, by appropriate discrete elements which represent alternately Lagrangian mass of water and Eulerian volume of space. Macroscale behavior of the model is verified by simulating undrained biaxial compression tests. Micro-scale behavior is compared to previous literature through pore pressure pattern visualization during shear tests. It is demonstrated that dynamic pore pressure patterns are generated by superposed stress waves. These pore-pressure patterns travel much faster than average drainage rate of the pore fluid and may initiate soil fabric change, ultimately leading to liquefaction in loose sands. Thus, this work demonstrates a tool to roughly link dvnamic stress wave patterns to initiation of liQuefaction nhenomena.

Numerical simulation of direct shear tests on mechanical properties of talus deposits based on self-adaptive PCNN digital image processing

The macro mechanical properties of materials with characteristics of large scale and complicated structural composition can be analyzed through its reconstructed meso-structures.In this work,the meso-structures of talus deposits that widely exist in the hydro-power engineering in the southwest of China were first reconstructed by small particles according to the in-situ photographs based on the self-adaptive PCNN digital image processing,and then numerical direct shear tests were carried out for studying the mechanical properties of talus deposits.Results indicate that the reconstructed meso-structures of talus deposits are more consistent with the actual situation because the self-adaptive PCNN digital image processing has a higher discrimination in the details of soil-rock segmentation.The existence and random distribution of rock blocks make the initial shear stiffness,the peak strength and the residual strength higher than those of the "pure soil" with particle size less than 1.25 cm apparently,but reduce the displacements required for the talus deposits reaching its peak shear strength.The increase of rock proportion causes a significant improvement in the internal friction angle of talus deposit,which to a certain degree leads to the characteristics of shear stress-displacement curves having a changing trend from the plastic strain softening deformation to the nonlinear strain hardening deformation,while an unconspicuous increase in cohesion.The uncertainty and heterogeneity of rock distributions cause the differences of rock proportion within shear zone,leading to a relatively strong fluctuation in peak strengths during the shear process,while movement features of rock blocks,such as translation,rotation and crossing,expand the scope of shear zone,increase the required shear force,and also directly lead to the misjudgment that the lower shear strength is obtained from the samples with high rock proportion.That,however,just explains the reason why the shear strength gained from a small amount of indoor test data is not consistent with engineering practice.

Simulation of particle flow in a bell-less type charging system of a blast furnace using the discrete element method

A three-dimensional model was established by the discrete element method （DEM） to analyze the flow and segregation of particles in a charging process in detail. The simulation results of the burden falling trajectory obtained by the model were compared with the industrial charging measurements to validate the applicability of the model. The flow behavior of particles from the weighing hopper to the top layer of a blast furnace and the heaping behavior were analyzed using this model. A radial segregation index （RSI） was used to evaluate the extent of the size segregation in the charging process. In addition, the influence of the chute inclination angle on the size segregation and burden profile during the charging process was investigated.

Effect of rolling friction on wall pressure,discharge velocity and outflow of granular material from a flat-bottomed bin

The present paper provides both experimental and DEM analyses of the filling and discharge of pea grains from a 3D flat-bottomed bin. In the DEM model, the fixed mean values of the experimentally determined single particle data, such as the particle density, Young＇s modulus, Poisson＇s ratio as well as the sliding and rolling friction coefficients were incorporated to analyse their effects on the macroscale indicators, such as the wall pressure, discharge velocities and material outflow parameters. The effect of rolling friction was studied based on the experimentally measured single particle rolling friction coefficient. This analysis is aimed at the quantitative prediction of flow parameters as related to the identification of material parameters.

Mechanics of granular column collapse in fluid at varying slope angles

This paper investigates the effect of initial volume fraction on the runout characteristics of collapse of granular columns on slopes in fluid. 2-D sub-grain scale numerical simulations are performed to understand the flow dynamics of granular collapse in fluid. The discrete element method（DEM） technique is coupled with the lattice Boltzmann method（LBM）, for fluid-grain interactions, to understand the evolution of submerged granular flows. The fluid phase is simulated using multiple-relaxation-time LBM（LBM-MRT） for numerical stability. In order to simulate interconnected pore space in 2-D, a reduction in the radius of the grains（hydrodynamic radius） is assumed during LBM computations. The collapse of granular column in fluid is compared with the dry cases to understand the effect of fluid on the runout behaviour. A parametric analysis is performed to assess the influence of the granular characteristics（initial packing） on the evolution of flow and run-out distances for slope angles of 0 °, 2.5°, 5 ° and 7.5 °. The granular flow dynamics is investigated by analysing the effect of hydroplaning, water entrainment and viscous drag on the granular mass. The mechanism of energy dissipation, shape of the flow front, water entrainment and evolution of packing density is used to explain the difference in the flow characteristics of loose and dense granular column collapse in fluid.

Particle shape consideration in numerical simulation of assemblies of irregularly shaped particles

The mechanical behavior of granular materials depends much on the shape of the constituent particles. Therefore appropriate modeling of particle, or grain, shape is quite important. This study employed the method of direct modeling of grain shape （Matsushima ＆ Saomto, 2002）, in which, the real shape of a grain is modeled by combining arbitrary number of overlapping circular elements which are connected to each other in a rigid way. Then, accordingly, a discrete-element program is used to simulate the assembly of grains. In order to measure the effects of grain shape on mechanical properties of assembly of grains, three types of grains-high angular grains, medium angular grains and round grains are considered where several biaxial tests are conducted on assemblies with different grain types, The results show that the angularity of grains greatly affects the behavior of granular soil.

Numerical analysis of enhanced mixing in a Gallay tote blender

The mixing performance of a multi-bladed baffle inserted into a traditional Gallay tote blender is explored by graphic processing unit-based discrete element method software. The mixing patterns and rates are investigated for a binary mixture, represented by two different colors, under several loading profiles. The baffle effectively enhances the convective mixing both in the axial and radial directions, because of the disturbance it causes to the initial flowing layer and solid-body zone, compared with a blender without a baffle. The axial mixing rate is affected by the gap between the baffle and the wall on the left and right sides, and an optimal blade length corresponds to the maximum mixing rate. However, the radial mixing rate increases with the blade length almost monotonically.