Post-liquefaction undrained monotonic behaviour of sands: experiments and DEM simulations

This paper presents the results of a series of servo-controlled cyclic triaxial tests and numerical simulations using the three- dimensional discrete element method (DEM) on post-liquefaction undrained monotonic strength of granular materials. In a first test series,undrained monotonic tests were carried out after dissipating the excess pore water pressure developed during liquefaction. The influence of different parameters such as amplitude of axial strain,relative density and confining pressure prior to liquefaction on the post-liquefaction undrained response have been investigated.The results obtained highlight an insignificant influence of amplitude of axial strain, confining pressure and a significant influence of relative density on the post-liquefaction undrained monotonic stress-strain response.In the second series, undrained monotonic tests were carried out on similar triaxial samples without dissipating the excess pore water pressure developed during liquefaction. The results highlight that the amplitude of axial strain prior to liquefaction has a significant influence on the post-liquefaction undrained monotonic response.In addition,DEM simulations have been carried out on an assembly of spheres to simulate post-liquefaction behaviour.The simulations were very similar to the experiments with an objective to understand the behaviour of monotonic strength of liquefied samples from the grain scale. The numerical simulations using DEM have captured qualitatively all the features of the post-liquefaction undrained monotonic response in a manner similar to that of the experiments.In addition,a detailed study on the evolution of micromechanical parameters such as the average coordination number and induced anisotropic coefficients has been reported during the post-liquefaction undrained monotonic loading.

Impact resistance and evaluation of retained strength on geotextiles

Over the last few decades, geotextiles have progressively been incorporated into geotechnical applications, especially in the field of coastal engineering. Geotextile materials often act as separator and a filter layer between rocks laid above and subgrade beneath. This versatile material has gradually substituted traditional granular materials because of its ease of installation, consistent quality and labour costefficiency. However, geotextiles often suffer damage during installation due to high dynamic bulk loading of rock placement. This can degrade geotextiles' mechanical strength. The properties considered in this paper include the impact resistance and retained strength of geotextiles. In general, the greater the impact energy applied to geotextiles, the greater the potential for damage. Results highlight the inadequacy of using index derived values as an indicator to determine geotextile performance on site because test results shows that geotextiles (staple fibre (SF) and continuous filament (CF)) with better mechanical properties did not outperform lower mechanical strength materials. The toughest CF product with a CBR index value of 9696N shows inferior impact resistance compared to SF product with the least CBR strength (2719N) given the same impact energy of 9.02 kJ. Test results also indicated that the reduction of strength for CF materials were much greater (between 20 and 50%) compared to SF materials (between 0 and 5%) when subjected to the same impact energy of 4.52 kJ.

Effect of relative density and confining pressure on Poisson ratio from bender and extender elements tests

By using the bender and extender elements tests, together with measurements of the travel times of shear (S) and primary (P) waves, the variation of Poisson ratio (nu) was determined for dry sands with respect to changes in relative densities and effective confining pressures (sigma(3)). The tests were performed for three different ranges of particle sizes. The magnitude of the Poisson ratio decreases invariably with an increase in both the relative density and the effective confining pressure. The effect of the confining pressure on the Poisson ratio was found to become relatively more significant for fine-grained sand as compared with the coarse-grained sand. For a given material, at a particular value of sigma(3), the magnitude of the Poisson ratio decreases, almost in a linear fashion, with an increase in the value of maximum shear modulus (G(max)). The two widely used correlations in literature, providing the relationships among G(max), void ratio (e) and effective confining pressure (sigma(3)), applicable for angular granular materials, were found to compare reasonably well with the present experimental data for the fine- and medium-grained sands. However, for the coarse-grained sand, these correlations tend to overestimate the values of G(max).

Dynamics of sheared inelastic dumbbells

We study the dynamical properties of the homogeneous shear flow of inelastic dumbbells in two dimensions as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are modelled as smooth fused disks characterized by the ratio of the distance between centres (L) and the disk diameter (D), with an aspect ratio (L/D) varying between 0 and 1 in our simulations. Area fractions studied are in the range 0.1-0.7, while coefficients of normal restitution (e(n)) from 0.99 to 0.7 are considered. The simulations use a modified form of the event-driven methodology for circular disks. The average orientation is characterized by an order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axes of all dumbbells in the same direction). We investigate power-law fits of S as a function of (L D) and (1 - e(n)(2)) There is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased. The order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. The mean energy of the velocity fluctuations in the flow direction is higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocities are Gaussian to a very good approximation. The pressure is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation, even at the lowest area fractions studied here. The mean angular velocity of the particles is equal to half the vorticity at low area fractions, but the magnitude systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.

Effect of Dilatation on the Elasto-Plastic Response of Bulk Metallic Glasses Under Indentation

Elasto-plastic response of bulk metallic glasses (BMGs) follows closely the response of granular materials through pressure dependent (or normal stress) yield locus and shear stress induced material dilatation. On a micro-structural level, material dilatation is responsible for stress softening and formation of localized shear band, however its influence on the macro-scale flow and deformation is largely unknown. In this work, we systematically analyze the effect of material dilatation on the gross indentation response of Zr-based BMG via finite element simulation. The strengthening/softening effect on the load-depth response and corresponding stress-strain profiles are presented in light of differences in elastic-plastic regimes under common indenters. Through comparison with existing experimental results, we draw conclusions regarding selection of suitable dilatation parameters for accurately predicting the gross response of BMGs

Microstructure, shape, stability and melting inembeddednanoparticles

The present article reviews some of the current work on a new class of materials which are nanoscale granular materials. We shall discuss in this paper two phase granular materials where one of the phases having nanometric dimension is embedded in a matrix of larger dimension. Known as nanoembedded materials, nanocomposites or ultrafine granular materials, this class of materials has attracted attention because of the opportunity of basic studies on the effect of size and embedding matrix on transformation behaviors as well as some novel properties, which include structural, magnetic and transport properties. These are in addition to the tremendous interests in what is known as quantum structures(embedded particles size less than 5 nm) for the case of semiconductors, which will not be discussed here. We shall primarily review the work done on metallic systems where the dispersed phases have low melting points and borrow extensively from the work done in our group. The phase transformations of the embedded particles show distinctive behavior and yield new insights. We shall first highlight briefly the strategy of synthesis of these materials by non-equilibrium processing techniques, which will be followed by examples where the effect of length scales on phase transformation behaviors like melting and solidification are discussed.

Taming of large diameter triaxial setup

Triaxial tests are essential to estimate the shear strength properties of the soil or rock. Normally triaxial tests are carried out on samples of 38 mm diameter and 76 mm height. Granular materials, predominantly used in base/sub-base construction of pavements or in railways have size range of 60-75 mm. Determination of shear strength parameters of those materials can be made possible only through triaxial tests on large diameter samples. This paper describes a large diameter cyclic triaxial testing facility set up in the Geotechnical Engineering lab of Indian Institute of Science. This setup consists of 100 kN capacity dynamic loading frame, which facilitates testing of samples of up to 300 mm diameter and 600 mm height. The loading ram can be actuated up to a maximum frequency of 10 Hz, with maximum amplitude of 100 mm. The setup is capable of carrying out static as well as dynamic triaxial tests under isotropic, anisotropic conditions with a maximum confining pressure of 1 MPa. Working with this setup is a difficult task because of the size of the sample. In this paper, a detailed discussion on the various problems encountered during the initial testing using the equipment, the ideas and solutions adopted to solve them are presented. Pilot experiments on granular sub-base material of 53 mm down size are also presented.

A generalized model for the effective thermal conductivity of porous media based on self-similarity

A generalized model for the effective thermal conductivity of porous media is derived based on the fact that statistical self-similarity exists in porous media. The proposed model assumes that porous media consist of two portions: randomly distributed non-touching particles and self-similarly distributed particles contacting each other with resistance. The latter are simulated by Sierpinski carpets with side length L = 13 and cutout size C = 3, 5, 7 and 9, respectively, depending upon the porosity concerned. Recursive formulae are presented and expressed as a function of porosity, ratio of areas, ratio of component thermal conductivities and contact resistance, and there is no empirical constant and every parameter has a clear physical meaning. The model predictions are compared with the existing experimental data, and good agreement is found in a wide range of porosity of 0.14-0.80, and this verifies the validity of the proposed model.

Dynamic characterization of micro-particle systems

Ordered granular systems have been a subject of active research for decades. Due to their rich dynamic response and nonlinearity, ordered granular systems have been suggested for several applications, such as solitary wave focusing, acoustic signals manipulation, and vibration absorption. Most of the fundamental research performed on ordered granular systems has focused on macro-scale examples. However, most engineering applications require these systems to operate at much smaller scales. Very little is known about the response of micro-scale granular systems, primarily because of the difficulties in realizing reliable and quantitative experiments, which originate from the discrete nature of granular materials and their highly nonlinear inter-particle contact forces.

In this work, we investigate the physics of ordered micro-granular systems by designing an innovative experimental platform that allows us to assemble, excite, and characterize ordered micro-granular systems. This new experimental platform employs a laser system to deliver impulses with controlled momentum and incorporates non-contact measurement apparatuses to detect the particles’ displacement and velocity. We demonstrated the capability of the laser system to excite systems of dry (stainless steel particles of radius 150 micrometers) and wet (silica particles of radius 3.69 micrometers, immersed in fluid) micro-particles, after which we analyzed the stress propagation through these systems.

We derived the equations of motion governing the dynamic response of dry and wet particles on a substrate, which we then validated in experiments. We then measured the losses in these systems and characterized the collision and friction between two micro-particles. We studied wave propagation in one-dimensional dry chains of micro-particles as well as in two-dimensional colloidal systems immersed in fluid. We investigated the influence of defects to wave propagation in the one-dimensional systems. Finally, we characterized the wave-attenuation and its relation to the viscosity of the surrounding fluid and performed computer simulations to establish a model that captures the observed response.

The findings of the study offer the first systematic experimental and numerical analysis of wave propagation through ordered systems of micro-particles. The experimental system designed in this work provides the necessary tools for further fundamental studies of wave propagation in both granular and colloidal systems.

Orientation, Flow, and Clogging in a Two-Dimensional Hopper: Ellipses vs. Disks

Two-dimensional (2D) hopper flow of disks has been extensively studied. Here, we investigate hopper flow of ellipses with aspect ratio $\alpha = 2$, and we contrast that behavior to the flow of disks. We use a quasi-2D hopper containing photoelastic particles to obtain stress/force information. We simultaneously measure the particle motion and stress. We determine several properties, including discharge rates, jamming probabilities, and the number of particles in clogging arches. For both particle types, the size of the opening, $D$, relative to the size of particles, $\ell$ is an important dimensionless measure. The orientation of the ellipses plays an important role in flow rheology and clogging. The alignment of contacting ellipses enhances the probability of forming stable arches. This study offers insight for applications involving the flow of granular materials consisting of ellipsoidal shapes, and possibly other non-spherical shapes.

Universal Quake Statistics: From Compressed Nanocrystals to Earthquakes.

Slowly-compressed single crystals, bulk metallic glasses (BMGs), rocks, granular materials, and the earth all deform via intermittent slips or "quakes". We find that although these systems span 12 decades in length scale, they all show the same scaling behavior for their slip size distributions and other statistical properties. Remarkably, the size distributions follow the same power law multiplied with the same exponential cutoff. The cutoff grows with applied force for materials spanning length scales from nanometers to kilometers. The tuneability of the cutoff with stress reflects "tuned critical" behavior, rather than self-organized criticality (SOC), which would imply stress-independence. A simple mean field model for avalanches of slipping weak spots explains the agreement across scales. It predicts the observed slip-size distributions and the observed stress-dependent cutoff function. The results enable extrapolations from one scale to another, and from one force to another, across different materials and structures, from nanocrystals to earthquakes.

Effect of machine dependence on wall failure property measurements, Jenike versus an annular shear tester

Annular, ring or torsional shear testers are commonly used in bulk solids handling research for the purpose of powder characterisation or equipment design. This paper reports from a DEFRA sponsored project which aims to develop an industrial powder flow-ability tester, (based on the annular shear tester) that is economic to buy and quick and easy to use in trained but unskilled hands. This paper compares the wall failure loci measured with an annular shear cell with measurements obtained using the accepted standard wall friction tester, the Jenike shear cell. These wall failure loci have been measured for several bulk solids which range from fine cohesive powders to free-flowing granular materials, on a stainless steel 304 2B wall surface.

A DEM study of silo discharge of a cohesive solid

Bulk handling of powders and granular solids is common in many industries and often gives rise to handling difficulties especially when the material exhibits complex cohesive behaviour. For example, high storage stresses in a silo can lead to high cohesive strength of the stored solid, which may in turn cause blockages such as ratholing or arching near the outlet during discharge. This paper presents a Discrete Element Method study of discharge of a granular solid with varying levels of cohesion from a flat-bottomed silo. The DEM simulations were conducted using the commercial EDEM code with a recently developed DEM contact model for cohesive solids implemented through an API. The contact model is based on an elasto-plastic contact with adhesion and uses hysteretic non-linear loading and unloading paths to model the elastic-plastic contact deformation. The adhesion parameter is a function of the maximum contact overlap. The model has been shown to be able to predict the stress history dependent behaviour depicted by a flow function of the material. The effects of cohesion on the discharge rate and flow pattern in the silo are investigated. The predicted discharge rates are compared for the varying levels of cohesion and the effect of adhesion is evaluated. The ability of the contact model to qualitatively predict the phenomena that are present in the discharge of a silo has been shown with the salient feature of mixed flow from a flat bottomed hopper identified in the simulation.

Experimental and discrete element modelling of cohesive iron ore fines

In this study, the behaviour of iron ore fines with varying levels of adhesion was investigated using a confined compression test and a uniaxial test. The uniaxial test was conducted using the semi-automated uniaxial EPT tester in which the cohesive strength of a bulk solid is evaluated from an unconfined compression test following a period of consolidation to a pre-defined vertical stress. The iron ore fines were also tested by measuring both the vertical and circumferential strains on the cylindrical container walls under vertical loading in a separate confined compression tester - the K0 tester, to determine the lateral pressure ratio. Discrete Element Method simulations of both experiments were carried out and the predictions were compared with the experimental observations. A recently developed DEM contact model for cohesive solids, an Elasto-Plastic Adhesive model, was used. This particle contact model uses hysteretic non-linear loading and unloading paths and an adhesion parameter which is a function of the maximum contact overlap. The model parameters for the simulations are phenomenologically based to reproduce the key bulk characteristics exhibited by the solid. The simulation results show a good agreement in capturing the stress history dependent behaviour depicted by the flow function of the cohesive iron ore fines while also providing a reasonably good match for the lateral pressure ratio observed during the confined compression K0 tests. This demonstrates the potential for the DEM model to be used in the simulation of bulk handling applications.

Desenvolvimento de um sistema misto de pasta à base de geopolímero e cimento portland para correção de perda de circulação

The materials engineering includes processes and products involving several areas of engineering, allowing them to prepare materials that fulfill the needs of various new products. In this case, this work aims to study a system composed of cement paste and geopolymers, which can contribute to solving an engineering problem that directly involves the exploitation of oil wells subject to loss of circulation. To correct it, has been already proposed the use of granular materials, fibers, reducing the drilling fluid or cement paste density and even surface and downhole mixed systems. In this work, we proposed the development of a slurry mixed system, the first was a cement-based slurry and the second a geopolymer-based slurry. The cement-based slurry was formulated with low density and extenders, 12.0 ppg (1.438 g/cm ³), showing great thixotropic characteristics. It was added nano silica at concentrations of 0.5, 1.0 and 1.5 gps (66.88, 133.76 and 200.64 L/m3) and CaCl2 at concentrations of 0.5, 1, 0 and 1.5%. The second system is a geopolymer-based paste formulated from molar ratios of 3.5 (nSiO2/nAl2O3), 0.27 (nK2O/nSiO2), 1.07 (nK2O/nAl2O3) and 13.99 (nH2O/nK2O). Finally, we performed a mixture of these two systems, for their application for correction of circulation lost. To characterize the raw materials, XRD, XRF, FTIR analysis and titration were performed. The both systems were characterized in tests based on API RP10B. Compressive strength tests were conducted after curing for 24 hours, 7 and 28 days at 58 °C on the cement-based system and the geopolymer-based system. From the mixtures have been performed mixability tests and micro structural characterizations (XRD, SEM and TG). The results showed that the nano silica, when combined with CaCl2 modified the rheological properties of the cement slurry and from the concentration of 1.5 gpc (200.64 L / m³) it was possible to obtain stable systems. The system mixture caused a change in the microstructure of the material by favoring the rate of geopolymer formation to hinder the C3S phase hydration, thus, the production of CSH phases and Portlandite were harmed. Through the mixability tests it can be concluded that the system, due to reduced setting time of the mixture, can be applied to plug lost circulation zones when mixed downhole