Patterns of shear zones in granular bodies within a polar hypoplastic continuum

Summary The paper deals with numerical investigations on the patterning of shear zones in granular bodies. The behavior of dry sand during plane strain compression tests was numerically modelled with a finite element method using a hypoplastic constitutive relation within a polar (Cosserat) continuum. The constitutive relation was obtained through an extension of a non-polar one by polar quantities, viz. rotations, curvatures, couple stresses using the mean grain diameter as a characteristic length. This relation can reproduce the essential features of granular bodies during shear localisation. During FE-calculations, the attention was laid on the influence of boundary conditions and the distribution of imperfections in the granular specimen on the formation of patterns of shear zones.     

Behaviour of granular bodies in induced shear zones

This paper deals with the behaviour of granular bodies in induced shear zones. Shearing of an infinite narrow layer of sand between two very rough boundaries under conditions of free dilatancy is numerically modelled with a finite element method and a hypoplastic constitutive relation within a polar (Cosserat) continuum. The relation can reproduce the essential features of granular bodies during shear localisation. The material constants are easily determined from element test results and can be estimated from granulometric properties. The attention is laid on the influence of the initial void ratio, mean grain diameter, layer height, pressure level and grain roughness on the thickness of induced shear zones. In addition, the effect of dilatancy constraint on shear localisation is investigated. The results are also compared to solutions within a non-polar continuum. The FE-calculations demonstrate that polar effects manifested by the appearance of grain rotations and couple stresses are significant in shear zones and their thickness is mainly affected by the initial void ratio, the mean grain diameter and the layer height. Received: 16 November 1999     

DEM simulation of shear vibrational fluidization of granular material

Fluidization of dry granular material is the transition from a solid state to a liquid state when sufficient energy is applied during vibration. This behavior is important because it is closely related to deformations of geotechnical structures during an earthquake. The scientific challenge lies in the understanding on how strain localization is related to the fluidization zone during the entire shearing process. Despite the importance of the mechanical behavior of granular material during fluidization, it cannot be easily characterized using traditional direct shear test. In this paper, 2D DEM model is firstly conduct, shear vibrational fluidization is defined for dry granular material, and the discrete element method has been used to simulate the direct shear test on granular material under vibrational loading during shearing. The peak, residual and vibro-residual shear strength envelopes have been obtained from the numerical simulations. Three distinct zones have been identified in the upper shear box based on the observed changes in volumetric strain before vibration. During vibration, fluidization occurs in the three zones with the characteristics that the shear stress, porosity, volumetric strain, and the coordination number drop to relatively lower values. During vibration, material becomes denser than the critical state, and strain localization has been relieved. Densification of the material at the shear zone leads to a strengthening of the material which increases the shearing resistance after vibration. Furthermore, a comparison of the 2D and 3D simulations is performed. Results reveal that the motion of particles in the out-of-plane direction in the 3D simulations lead to smoother shear stress and more consistent with the experimental result.     

被动围压条件下岩石材料冲击压缩试验研究

为研究煤矿岩石材料被动围压条件下动态力学性能和变形破坏规律,利用Ø50mm变截面分离式Hopkinson压杆(SHPB)试验装置,对45#钢质套筒环向约束状态下煤矿岩石试件进行了不同加载速率冲击压缩试验。试验结果表明：被动围压条件下SHPB试验中,岩石试件的材料延性和抗破坏能力均得到增强,试件轴向应力是采用同种加载条件无围压SHPB试验时的1.2倍,破坏应变比无围压SHPB试验提高2～3倍,且径向应力随轴向应变增大总体呈上升趋势,试件破坏为压剪破坏模式,与无围压SHPB试验有所不同。     

DEM modeling of shear bands in crushable and irregularly shaped granular materials

A method of modeling convex or concave polygonal particles is proposed. DEM simulations of shear banding in crushable and irregularly shaped granular materials are presented in this work. Numerical biaxial tests are conducted on an identical particle assembly with varied particle crushability. The particle crushing is synchronized with the development of macroscopic stress, and the evolution of particle size distribution can be characterized by fractal dimension. The shear banding pattern is sensitive to particle crushability, where one shear band is clearly visible in the uncrushable assembly and X-shaped shear bands are evident in the crushable assembly. There are fewer branches of strong force chains and weak confinement inside the shear bands, which cause the particles inside the shear bands to become vulnerable to breakage. The small fragments with larger rotation magnitudes inside the shear bands form ball-bearing to promote the formation of shear bands. While there are extensive particle breakages occurring, the ball-bearing mechanism will lubricate whole assembly. With the increase of particle crushability the shear band formation is suppressed and the shear resistance of the assembly is reduced. The porosity inside the shear bands are related to the particle crushability.     

DEM analysis of swelling behaviour in granular media

Many granular materials change their volume as they absorb fluids. This phenomenon is called swelling and can be observed in a variety of solids, such as soils, wood, absorbent hygiene products (AHPs) and pharmaceutical excipients. Therefore, an in-depth understanding of grain swelling is of great importance. Since experimental investigations can often provide only limited information, while great insight could be gained from numerical modelling, rigorous numerical models for predicting particle swelling are required. Hence, the objective of this research is to develop and validate a Discrete Element Method (DEM) for swelling of particles. A first order kinetic model was employed to predict the volume expansion of a single grain and subsequently implemented in DEM. The validation of the model was accomplished by comparing the expansion with time of a packed bed made of super absorbent polymer (SAP) particles obtained numerically and experimentally. It was shown that the DEM model can accurately predict the bed expansion. The model was then employed to simulate the swelling of three different materials: superabsorbent polymer (SAP), rice and microcrystalline cellulose (MCC Avicel PH102). As expected, it is demonstrated that the material properties play a significant role on the swelling; the fastest to reach its maximum expansion is the granular bed made of MCC PH102, followed by SAP and rice. However, the highest swelling capacity is achieved with SAP. Moreover, a preliminary DEM analysis of the segregation in a swelling binary mixture is presented in this work. Results suggest that systems which contain a small number of particles, and thus are looser, are more prone to segregation. Future study could advance the developed model to analyse consequences of swelling phenomena in granular materials, such as segregation and heat generation.     

DEM simulation of the shear behaviour of breakable granular materials with various angularities

Particle shape is an important factor that affects particle breakage and the mechanical behaviour of granular materials. This report explored the effect of angularity on the mechanical behaviour of breakable granular materials under triaxial tests. Various angular particles are generated using the quasi-spherical polyhedron method. The angularity α is defined as the mean exterior angle of touching faces in a particle model. A breakable particle is constructed as an aggregate composed of coplanar and glued Voronoi polyhedra. After being prepared under the densest conditions, all assemblies were subjected to triaxial compression until a critical state was reached. The macroscopic characteristics, including the shear strength and dilatancy response, were investigated. Then, particle breakage characteristics, including the extent of particle breakage, breakage pattern and correlation between the particle breakage and energy input, were evaluated. Furthermore, the microscopic characteristics, including the contact force and fabric anisotropy, were examined to probe the microscopic origins of the shear strength. As α increases, the peak shear strength increases first and then remains constant, while the critical shear strength generally increases. Assemblies with larger angularity tend to cause more serious particle breakage. The relative breakage is linearly correlated with α under shear loading. Compared with unbreakable particles, the peak shear strength and the critical volumetric strain decline, and the degree of decline linearly increases with increasing α.     

FE-investigation of shear localization in granular bodies under high shear rate

Shear localization in granular materials under high shear rate is analysed with the finite element method and a micro-polar hypoplastic constitutive model enhanced by viscous terms. We consider plane strain shearing of an infinitely long and narrow granular strip of initially dense sand between two very rough walls under conditions of free dilatancy. The constitutive model can reproduce the essential features of granular materials during shear localization. The calculations are performed under quasi-static and dynamic conditions with different shear rates. In dynamic regime, the viscosity terms are formulated based on a modified Newtonian fluid and according to the formula by Stadler and Buggisch (Proceedings of the conference on Reliable flow of particulate solids, EFCE Pub. Series, vol 49. Chr. Michelsen Institute, Bergen, 1985). Emphasis is given to the influence of inertial and viscous forces on the shear zone thickness and mobilized wall friction angle.     

Tangential velocity profiles of granular material within horizontal rotating cylinders modelled using the DEM

Flow regimes of granular materials in horizontal rotating cylinders are industrially important since they have a strong influence on the rates of heat and mass transfer within these systems. The tangential velocity profile, which describes how the average particle velocity in the direction parallel to the surface of the bed varies along a radius perpendicular to the surface of the bed, has been examined for many experimental and simulated systems. This paper is concerned with tangential velocity profiles within rotating cylinders simulated using the discrete element method. For high fill levels good agreement is found between the simulated velocity profiles and the equation proposed by Nakagawa et al. (Exp Fluids 16:54–60, 1993) based on magnetic resonance measurements. At lower fill levels slip is observed between the cylinder wall and the particles in contact with it and also between the outer layer of particles and the bulk of the bed. It is demonstrated that this slip occurs when the particles in contact with the wall are able to rotate and that it may be prevented either by using non-spherical particles or by attaching “lifters” to the cylinder wall.     

Numerical analysis of critical state behaviors of granular soils under different loading conditions

Discrete element method is used to simulate granular assembly behaviors with different initial conditions under three different loading conditions—plain strain, conventional triaxial compression, and direct shear. Different deformation modes of specimens with different conditions are presented. Some important parameters of the critical state theory are investigated. Uniqueness of the critical state line is checked which shows that there is no a unique critical state line for specimens with different initial void ratios under different loading conditions. Frictional angles and dilation angles of specimens with different conditions at critical state are compared. Void ratios and coordination numbers of specimens at critical state are studied. Anisotropies of the particle orientation and normal contact force at initial state, critical state, as well as the evolutions during shearing are analyzed. The anisotropy is shown to have significant effects on the soil behaviors and is related to the non-uniqueness of the critical state line. The developed numerical models can be used to study the micromechanics and microstructure of the specimen subjected to different loading conditions in the future.     

Study of the shear behavior of binary granular materials by DEM simulations and experimental triaxial tests

This paper aims at studying the shear behavior of mixtures of fine and coarse particles by classical triaxial tests. The work is performed both on experimental tests and computer simulations by discrete element method. The comparisons between experimental and simulation results on monosized and binary samples show that the DEM model can reproduce deviatoric curves satisfactorily in experimental conditions. The shear behavior of monosized and binary systems with the same initial void ratio differs significantly, suggesting that the state of compaction of the system is more influential than the initial void ratio. Comparison between compacted and uncompacted samples confirms that compaction increases the shear strength of granular matter. At the particle scale, the coordination number decreases with the augmentation of the volume fraction of coarse particles. The average rotation velocity of fine particles is higher than coarse particles, but their particle stress tensor is smaller than coarse ones.     

DEM investigation on the evolution of microstructure in granular soils under shearing

The mechanical behaviors of granular soils at different initial densities and confining pressures in the drained and undrained triaxial tests are investigated micromechanically by three-dimensional discrete element method (DEM). The evolutions of the microstructure in the numerical specimen, including coordination number, contact force and anisotropies of contact normal and contact force, are monitored during the shearing. The typical shear behaviors of granular soils (e.g. strain softening, phase transformation, static liquefaction and critical state behavior) are successfully captured in the DEM simulation. It is found that the anisotropies of contact normal, normal and tangential contact forces comprise the shear resistance and show different evolution features during shearing. After large strain shearing, the microstructure of the soil will finally reach a critical state, although the evolution path depends on the soil density and loading mode. Similar to the macroscopic void ratio $e$ and deviatoric stress $q$ , the coordination number and anisotropies of contact normal and contact force at the critical state also depend on the mean normal effective stress $P^{\prime }$ at the critical state.     

DEM analysis of the size effects on the behavior of crushable granular materials

Four sets of individual-particle crushing tests were carried out on sandstone grains of different size with geometric similarity. The tensile strength was analyzed using Weibull statistics, and the size-hardening law was obtained. The experimental data also validated that the Weibull modulus is independent of the grain size. Considering both the shear and tensile fracture modes of the particle, the Mohr–Coulomb model with a tension cut-off was employed as the fracture criterion of a single particle. When the particle stresses satisfied the fracture criterion, three new fragments modeled by the ‘clump’ were generated to replace the broken particle. Nine spheres with four different sizes were released from the clump and allowed to continue crushing if the fragment stresses fulfilled the criterion again. Two polydisperse assemblies with different particle sizes but same initial fabrics were prepared. DEM simulations of triaxial shear tests with different grain sizes were carried out on the crushable granular material with varied confining pressures. The simulated stress–strain–dilation responses were in agreement with the experimental observations. The macro–micro responses of the two samples, including the stress–strain–dilation behavior, the particle crushing, and the normal contact force distribution, were discussed in detail. The cause of the size effect on the shear strength and deformation was thoroughly investigated through a variety of mechanism demonstrations and micromechanical analysis.     

Effect of periodic boundary conditions on granular motion in horizontal rotating cylinders modelled using the DEM

Applying periodic boundary conditions to DEM simulations of granular motion in horizontal, rotating cylinders can result in unexpected bulk axial flow. This axial flow is shown to contribute significantly to the mean square deviation in particle position and consequently must be considered in the analysis of processes governed by slow axial motion such as axial dispersion.     

Effect of vibration on solid-to-liquid transition in small granular systems under shear

The effect of vibration on the solid-to-liquid-like transition of a dense granular assembly under planar shear is studied numerically using soft particle molecular dynamics simulations in two dimensions. We focus on small systems in a thin planar Couette cell, examining the bistable region while increasing shear, with varying amounts of vertical vibration, and determine statistics of the shear required for fluidization. In the absence of vibration, the threshold value of the shear stress depends on the preparation of the system and has a broad distribution. However, adding periodic vibration both lowers the mean fluidization threshold value of the shear stress and decreased its variability. A previous study performed similar simulations using random noise; the results from these two studies exhibit excellent agreement with proper normalization over appropriate ranges of parameters.     

Micromechanical behavior of granular materials with inherent anisotropy under cyclic loading using 2D DEM

This paper presents the micromechanical behavior of granular materials due to different initial inherent anisotropic conditions during cyclic loading using the discrete element method (DEM). Oval particles were used to model the samples. Three samples, with three different inherent anisotropic conditions based on the particle’s bedding direction, were prepared and subjected to biaxial cyclic loading. The differences in the inherent anisotropic conditions of the samples affect the stress–strain-dilative behavior of granular materials. The width of the stress–strain cyclic loops decreases as the preferred bedding angle changes from vertical to horizontal. Contact fabric evolution is found to be dependent on the inherent anisotropic fabric of the sample during loading and unloading. The fabric anisotropy is dominant for horizontal particle bedding at the end of loading and for vertical particle bedding at the end of unloading. A change in fabric anisotropy is observed only for the first few loading–unloading cycles for the given conditions depicted in the present study.     

A numerical analysis of passive scalar transport in turbulent shear flows

Abstract

The development of the formation and vortex pairing process in a two‐dimensional shear flow and the associated passive scalar (mass concentration or energy) transport process was numerically simulated by using the Vortex‐in‐Cell (VIC) Method combined with the Upwind Finite Difference Method. The visualized temporal distributions of passive scalars resemble the vortex structures and the turbulent passive scalar fluxes showed a definite connection with the occurrence of entrainment during the formation and pairing interaction of large‐scale vortex structures. The profiles of spatial‐averaged passive scalar ø, turbulent passive scalar fluxes, u'ø’ and v'ø’, turbulent diffusivity of mean‐squared scalar fluctuation, v'ø‘ ², mean‐squared turbulent passive scalar fluctuation, √ø‘ ², skewness, and flatness factor of the probability density function of scalar fluctuation ø’ at three different times are calculated. With the lateral dimension scaled by the momentum thickness and the velocity scaled by the velocity difference across the shear layer, these profiles were shown to be self‐preserved. The probability density function of turbulent scalar fluctuation was found to be asymmetric and double‐peaked.     

In situ analysis of delayed fibre failure within water-aged GFRP under static fatigue conditions

A stress corrosion model has been applied to the microscopic analysis of the delayed fibre failure processes occurring within a water-aged unidirectional glass/epoxy composite under static fatigue loading (i.e. relaxation). By means of in situ microscopic observations, the individual fibre failures within an elementary volume located on the tensile side of the flexural specimens have been quantified as a function of time under various applied strain levels. It was found that the time dependence of the in situ fibre failure processes obeyed a stress corrosion model. From the microscopic observations, it was possible to assess consistent values of the parameters characterising the in situ fibre strength distribution and the subcritical crack propagation law. A comparison with separate static fatigue experiments using unimpregnated fibre bundles demonstrated that the specific physico-chemical environment encountered by the glass fibres within the aged epoxy matrix can induce significant changes in the subcritical crack propagation rates, as compared to stress corrosion cracking data collected in humid air.     

Analysis of the influence of crushing on the behavior of granular materials under shear

The behavior of granular materials mainly depends on the mechanical and engineering properties of particles in its structural matrix. Crushing or breakage of granular materials under compression or shear occurs when the energy available is sufficient to overcome the resistance of the material. Relatively little systematic research has been conducted regarding how to evaluate or quantify particle crushing and how it effects the engineering properties of the granular materials. The aim of this study is to investigate the effect of crushing on the bulk behavior of granular materials by using manufactured granular materials (MGM) rather than using a naturally occurring cohesionless granular material. MGM allow changing only one particle parameter, namely the “crushing strength”. Four different categories of MGM (with different crushing strength) are used to study the effect on the bulk shear strength, stiffness modulus, friction and dilatancy angle “engineering properties”. A substantial influence on the stress–strain behavior and engineering properties of granular materials is observed. Higher confining stress causes some non-uniformity (strong variations/jumps) in volumetric strain and a constant volumetric strain is not always observed under large shear deformations due to crushing, i.e. there is no critical state with flow regime (with constant volumetric strain).