Effects of Soil Fabric on Undrained Behavior of Sands

The undrained behavior of sands in monotonic triaxial compression and extension tests was simulated using the Distinct Element Method (DEM). Soil specimens were prepared at different initial soil fabrics but at similar void ratios and the effects of soil fabric on the undrained behavior were investigated. The DEM results show that soil fabric and its change have profound effects on the undrained response of sands. They also provide some insights for the interpretation of the published experimental data that show the effects of specimen reconstitution methods and of preshearing on the undrained behavior of sands

Three-dimensional shakedown solutions for anisotropic cohesive-frictional materials under moving surface loads

Previous work on three-dimensional shakedown analysis of cohesive-frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three-dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower-bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross-anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only

Influence of packing density and stress on the dynamic response of granular materials

Laboratory geophysics tests including bender elements and acoustic emission measure the speed of propagation of stress or sound waves in granular materials to derive elastic stiffness parameters. This contribution builds on earlier studies to assess whether the received signal characteristics can provide additional information about either the material’s behaviour or the nature of the material itself. Specifically it considers the maximum frequency that the material can transmit; it also assesses whether there is a simple link between the spectrum of the received signal and the natural frequencies of the sample. Discrete element method (DEM) simulations of planar compression wave propagation were performed to generate the data for the study. Restricting consideration to uniform (monodisperse) spheres, the material fabric was varied by considering face-centred cubic lattice packings as well as random configurations with different packing densities. Supplemental analyses, in addition to the DEM simulations, were used to develop a more comprehensive understanding of the system dynamics. The assembly stiffness and mass matrices were extracted from the DEM model and these data were used in an eigenmode analysis that provided significant insight into the observed overall dynamic response. The close agreement of the wave velocities estimated using eigenmode analysis with the DEM results confirms that DEM wave propagation simulations can reliably be used to extract material stiffness data. The data show that increasing either stress or density allows higher frequencies to propagate through the media, but the low-pass wavelength is a function of packing density rather than stress level. Prior research which had hypothesised that there is a simple link between the spectrum of the received signal and the natural sample frequencies was not substantiated

Calibration of linear contact stiffnesses in discrete element models using a hybrid analytical-computational framework

Efficient selections of particle-scale contact parameters in discrete element modelling remain an open question. The aim of this study is to provide a hybrid calibration framework to estimate linear contact stiffnesses (normal and tangential) for both two-dimensional and three-dimensional simulations. Analytical formulas linking macroscopic parameters (Young's modulus, Poisson's ratio) to mesoscopic particle parameters for granular systems are derived based on statistically isotropic packings under small-strain isotropic stress conditions. By taking the derived analytical solutions as initial approximations, the gradient descent algorithm automatically obtains a reliable numerical estimation. The proposed framework is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. The results show that this hybrid method practically reduces the time for artificial trials and errors to obtain reasonable stiffness parameters. The proposed framework can be extended to other parameter calibration problems in DEM