Image based simulation of one-dimensional compression tests on carbonate sand

High factors of safety and conservative methods are commonly used on foundation design on shelly carbonate soils. A better understanding of the behavior of this material is, thus, critical for more sustainable approaches for the design of a number of offshore structures and submarine pipelines. In particular, understanding the physical phenomena taking place at the microscale has the potential to spur the development of robust computational methods. In this study, a one-dimensional compression test was performed inside an X-ray scanner to obtain 3D images of the evolving internal structure of a shelly carbonate sand. A preliminary inspection of the images through five loading increments has shown that the grains rearrange under loading and in some cases cracks develop at the contacts. In order to replicate of the experiments in the numerical domain, the 3D image of the soil prior to loading was imported into a micro Finite Element (µFE) framework. This image-based modelling tool enables measurements of the contact force and stress map inside the grains while making use of the real microstructure of the soil. The potential of the µFE model to contribute insights into yield initiation within the grain is demonstrated here. This is of particular interest to better understand the breakage of shelly grains underpinning their highly compressive behavior

On the tensile strength of soil grains in Hertzian response

The breakage initiation of soil grains is controlled by its tensile capacity. Despite the importance of tensile strength, it is often disregarded due to difficulties in measurement. This paper presents an experimental and numerical investigation on the effect of tensile strength on Hertzian response of a single soil grain. Hertz theory is commonly used in numerical simulation to present the contact constitutive behaviour of a purely elastic grain under normal loading. This normal force:displacement comes from stress distribution and concentration inside the grain. When the stress reaches the tensile capacity, a crack initiates. A series of numerical tests have been conducted to determine the sensitivity of Hertzian response to the selected tensile strength used as an input data. An elastic-damage constitutive model has been employed for spherical grains in a combined finite-discrete element framework. The interpretation of results was enriched by considering previous theoretical work. In addition, systematic experimental tests have been carried out on both spherical glass beads and grains of two different sands, i.e. Leighton Buzzard silica sand and coarse carbonate sand from Persian Gulf. The preliminary results suggest that lower tensile strength leads to a softer response under normal loading. The wider range of responses obtained for the carbonate sand, are believed to be related to the large variety of grain shape associated with bioclastic origin of the constituent grains

On the torsional loading of elastoplastic spheres in contact

The mechanical interaction between two bodies involves normal loading in combination with tangential, torsional and rotational loading. This paper focuses on the torsional loading of two spherical bodies which leads to twisting moment. The theoretical approach for calculating twisting moment between two spherical bodies has been proposed by Lubkin [1]. Due to the complexity of the solution, this has been simplified by Deresiewicz for discrete element modelling [2]. Here, the application of a simplified model for elastoplastic spheres is verified using computational modelling. The single grain interaction is simulated in a combined finite discrete element domain. In this domain a grain can deform using a finite element formulation and can interact with other objects based on discrete element principles. For an elastoplastic model, the contact area is larger in comparison with the elastic model, under a given normal force. Therefore, the plastic twisting moment is stiffer. The results presented here are important for describing any granular system involving torsional loading of elastoplastic grains. In particular, recent research on the behaviour of soil has clearly shown the importance of plasticity on grain interaction and rearrangement

Comparison between a μFE model and DEM for an assembly of spheres under triaxial compression

This paper presents a simple case of a Face Centred Cubic (FCC) array of 2,000 spheres under triaxial compression to compare the results obtained using the Discrete Element Method (DEM) and a micro finite element model (μFE). This μFE approach was developed so that the internal structure of the soil can be obtained using x-ray computed tomography and converted into a numerical fabric. The individual grains are represented as continuum deformable bodies and the inter-granular interaction based on the defined contact laws. In order to demonstrate the simple contact constitutive behaviour used in this μFE model, the response for two contacting elastic spheres is compared with theoretical equations. The strength at failure of the packing of 2,000 spheres is seen to yield similar values for DEM, μFE and the analytical solution. When comparing the evolving void ratio, a good agreement between the two numerical models was observed for very small strains but as the strain increases, the values start to diverge, which is believed to be related with the rigidity of the grains used in DEM

CFD modelling of the effect of capillary pressure on retention behaviour of water menisci at inter-particle contacts

This paper presents a Computational Fluid Dynamics (CFD) model on the effect of capillary pressure on the retention behaviour of a granular material. The model proposes an unprecedented CFD insight into the onset of liquid menisci at the inter-particles contact under varying hydraulic conditions. The present work models the material grains as smooth spherical particles that define a porous network filled by two interstitial fluids: air and silicon oil. The numerical model has been subsequently validated against experimental measurements of the degree of saturation at different capillary pressures taken by Dullien et al. [F.A. Dullien, C. Zarcone, I.F. MacDonald, A. Collins, R.D. Bochard. J. Colloid Interface Sci. 127, 2 (1989)] in a system of smooth glass beads flooded with silicon oil. Results from the numerical simulations confirm the good capability of the model to reproduce the experimental retention behaviour of the granular material. Finally, the present paper laid the basis for future CFD studies on the effect of various factors (e.g. hydraulic hysteresis, surface roughness and/or grain shape) on the capillary pressure acting at the interparticle contact

Softening-based interface model and nonlinear load-settlement response analysis of piles in saturated and unsaturated multi-layered soils

This work presents a simplified method for the nonlinear analysis of the load–displacement response of piles in multi-layered soils. As a starting step, a new interface model based on the disturbed state concept (DSC) is put forth to simulate the interface shear stress-displacement relationship by considering the nonlinear hardening–softening behaviour. In the new model, input parameters can be conveniently calibrated using conventional interface shear tests or on-site tests. The good agreement between predictions and experimental data from interface direct shear tests validated the performance of the proposed DSC model. The DSC model performed better in terms of predictions when compared to the hyperbolic one. Next, the soil-structure interface model and bearing capacity theory are coupled to provide a theoretical framework for the analysis of pile load-transfer in saturated and unsaturated multi-layered soils, where the DSC model is employed to represent base resistance as well as skin friction. This work also discusses the profile of steady-state in-situ matric suction, soil–water characteristic curve, and pore-water pressure of unsaturated soils. The proposed method has the advantage of being used in practice as it is simple to obtain input parameters from laboratory tests, as well as Standard Penetration or Cone Penetration Tests. The proposed framework is finally applied to the analysis of five well-documented case studies. The proposed approach and the static load test results from the field measurements are found to be in satisfactory agreement, indicating that the proposed method performs well. The proposed method is suggested to be utilised for preliminary analysis, planning a suitable programme of loading tests, as well as optimizing the pile design by back analysis of the load test results.</p

Characterisation of physical and mechanical properties of seven particulate materials proposed as traction enhancers

Abstract Particulate materials are utilised in many applications to manipulate the friction between surfaces. This dataset provides the characteristics of particulates used as rail sand in the train’s wheel/rail interface (via an on-board system) to facilitate the train’s acceleration and deceleration. Seven materials are studied including Austrian rail sand, standard Great British rail sand, waste glass beads, recycled crushed glass, non-coated alumina, coated alumina, and dolomite. The main objective of this research is to provide a physical and mechanical characterisation of these granular materials in terms of their density, bulk behaviour, particle size, particle shape, hardness, reduced modulus, and mineralogical properties. In particular, three-dimensional raw and post-processed micro-computed tomography images of more than 1200 particles are shared. The results provide a detailed dataset which can be used in ongoing and future experimental and numerical investigations studying the role of particulates in the wheel/rail interface

Imaging the root–rhizosphere interface using micro computed tomography: quantifying void ratio and root volume ratio profiles

Root growth alters soil fabric and consequently its mechanical and physical properties. Recent studies show that roots induce compaction of soil in their immediate vicinity, a region that is central for plant health. However, high quality quantification of root influence on the soil fabric, able to inform computational models is lacking from the literature. This study quantifies the relationship between soil physical characteristics and root growth, giving special emphasis on how roots in early stage formation influence the physical architecture of the surrounding soil structure. High-resolution X-ray micro-Computed Tomography (µCT) is used to acquire three dimensional images of two homogeneously-packed samples. It is observed that the void ratio profile extending from the soil-root interface into the bulk soil is altered by root growth. The roots considerably modify the immediate soil physical characteristics by creating micro cracks at the soil-root interface and by increasing void ratio. This paper presents the mechanisms that led to the observed structure as well as some of the implications that it has in such a dynamic zone

Imaging the root–rhizosphere interface using micro computed tomography: quantifying void ratio and root volume ratio profiles

Root growth alters soil fabric and consequently its mechanical and physical properties. Recent studies show that roots induce compaction of soil in their immediate vicinity, a region that is central for plant health. However, high quality quantification of root influence on the soil fabric, able to inform computational models is lacking from the literature. This study quantifies the relationship between soil physical characteristics and root growth, giving special emphasis on how roots in early stage formation influence the physical architecture of the surrounding soil structure. High-resolution X-ray micro-Computed Tomography (µCT) is used to acquire three dimensional images of two homogeneously-packed samples. It is observed that the void ratio profile extending from the soil-root interface into the bulk soil is altered by root growth. The roots considerably modify the immediate soil physical characteristics by creating micro cracks at the soil-root interface and by increasing void ratio. This paper presents the mechanisms that led to the observed structure as well as some of the implications that it has in such a dynamic zone