Search

Particle morphology is of pivotal importance in granular materials at different length scales. This thesis begins with quantification and reconstruction of real particle shapes, then studies various kinds of granular material behaviour influenced by morphology features. The project background and related previous work are introduced in Chapters 1 and 2, respectively. Chapter 3 is firstly concentrated on reconstructing particle shapes using Spherical Harmonics (SH). DEM clumps are also used to approximate realistic particle shapes to simulate sand column collapse. The efficiency of SH analysis between different coordinate systems in FEM mesh qualities and reconstruction are also compared. Chapter 4 is on contact behaviour of rough spheres with different morphology features. Via FEM simulation, benchmarked with Hertzian solution, a semi-analytical model is proposed. In Chapter 5, to highlight the influences of contact curvature on single particle crushing behaviour, a rotational point loading in FDEM, benchmarked with an in-situ XCT experiment, is proposed. In Chapter 6, experiments of failure modes of cemented sands under different loading paths and the quantification of fracture fabric are conducted. FDEM simulations, benchmarked with a combined in-situ experiment with XCT and diffraction, of meso-scale concrete are also conducted. Effects of realistic aggregate morphology on overall concrete are investigated. In Chapter 7, permeability of uniformly graded 3D printed grains, is experimentally studied. Modified Kozney-Carman equation is also proposed. Chapter 8 concludes this work by summarising the findings and implications and provides an outlook on future works. This dissertation presents a new comprehensive understanding of granular behaviour influenced by its morphology features. Via the proposed framework combining both experimental and numerical information, it is helpful to design and optimise of various granular materials with specific morphology features.

Glacier surges are quasi-periodic episodes of rapid ice flow that arise from increases in slip rate at the ice–bed interface. The mechanisms that trigger and sustain surges are not well understood. Here, we develop a new model of incipient surge motion for glaciers underlain by sediments to explore how surges may arise from slip instabilities within a thin layer of saturated, deforming subglacial till. Our model represents the evolution of internal friction, porosity and pore water pressure within the till as functions of the rate and history of shear deformation, and couples the till mechanics to a simple ice-flow model. Changes in pore water pressure govern incipient surge motion, with less permeable till facilitating surging because dilation-driven reductions in pore water pressure slow the rate at which till tends towards a new steady state, thereby allowing time for the glacier to thin dynamically. The reduction of overburden (and thus effective) pressure at the bed caused by dynamic thinning of the glacier sustains surge acceleration in our model. The need for changes in both the hydromechanical properties of the till and the thickness of the glacier creates restrictive conditions for surge motion that are consistent with the rarity of surge-type glaciers and their geographical clustering.

The complexity and multiscale nature of material microstructures introduces significant intricacies to many mechanics problems for which we do not have a full theoretical understanding. Under loading, these microstructures can introduce significant nonlinearities that cannot be described sufficiently by current theories and models. This leads us to consider experiments we could perform to improve our understanding of such effects. This thesis describes the design of experiments exploring two aspects of material microstructure effects: (i) crack propagation and renucleation in soft brittle polymers and (ii) interparticle forces in granular materials. First, experimental and analysis methods are developed to study fracture mechanics in soft brittle polymers with the goal of developing a more detailed understanding of the effects of microstructural heterogeneities on crack propagation and renucleation in three-dimensions. To better understand these processes, experiments on crack propagation in thin soft polymers using confocal microscopy images are conducted. Traditional metrics associated with crack propagation including stress intensity factor (SIF, K) and energy release rate (ERR, G) are calculated by direct measurement of the crack tip opening displacement (CTOD, δt) on the sub-millimeter scale. Errors in these calculations are comparable to those reported in the literature for more traditional fracture experiment geometries. Fluorescent speckle images are captured using confocal microscopy imaging, a fast and low cost 3D optical imaging technique, to study crack geometry during propagation. Images of renucleation events are also captured allowing investigation of factors contributing to slow crack roughening observed by earlier researchers. The goal of this study is to provide an experimental method to enhance understanding of crack interactions with microstructural heterogeneities and of renucleation events, which can significantly improve our ability to design material toughness. To begin to understand the effects of engineered microstructural heterogeneities such as inclusions in materials, we must be able to produce such engineered systems and understand the interparticle interactions. To this end, a method to manufacture volumetrically speckled spheres in-house with controlled diameters was developed. Additionally, an experimental method combining confocal microscopy with digital volume correlation (DVC) was also used to study interparticle force transmission in 3D. Analysis of an in-plane 2D projection of volumetric surface data shows that three-dimensional effects play a significant role in the deformation of granular assemblies. Study of a single grain in 3D demonstrates progress in experimental capabilities and highlights the need for more studies to validate existing numerical models and theories for granular matter. Analysis of particle scale deformations and strains with the Granular Element Method (GEM) allows us to determine interparticle forces and understand the development and evolution of force chains in a granular assembly under a wide variety of loading conditions. These experiments can also lead to development of new understanding of the effects of inclusions on material properties, processes, and damage evolution.

The relative equilibria for the spherical, finite density three-body problem are identified. Specifically, there are 28 distinct relative equilibria in this problem which include the classical five relative equilibria for the point-mass three-body problem. None of the identified relative equilibria exist or are stable over all values of angular momentum. The stability and bifurcation pathways of these relative equilibria are mapped out as the angular momentum of the system is increased. This is done under the assumption that they have equal and constant densities and that the entire system rotates about its maximum moment of inertia. The transition to finite density greatly increases the number of relative equilibria in the three-body problem and ensures that minimum energy configurations exist for all values of angular momentum.

A method for modeling cubically anisotropic elasticity within the discrete element method is presented. The discrete element method (DEM) is an approach originally intended for modeling granular materials (sand, soil, and powders); however, recent developments have usefully extended it to model stochastic mechanical processes in monolithic solids which, to date, have been assumed to be elastically isotropic. The method presented here for efficiently capturing cubic elasticity in DEM is an important prerequisite for further extending DEM to capture the influence of elastic anisotropy on the mechanical response of polycrystals, composites, etc. The system demonstrated here uses a directionally assigned stiffness in the bonds between adjacent elements and includes separate schemes for achieving anisotropy with Zener ratios greater and smaller than one. The model framework is presented along with an analysis of the accessible space of elastic properties that can be modeled and an artificial neural network interpolation scheme for mapping input parameters to model elastic behavior.

The available studies in the literature which dealt with the scale effects of strip footings on different sand packing systematically still remain scarce. In this research, the variation of ultimate bearing capacity and deformation pattern of soil beneath strip footings of different widths under plane-strain condition on the surface of loose, medium-dense and dense sand have been systematically studied using experimental and noninvasive methods for measuring microscopic deformations. The presented analyses are based on model scale compression test analysed using Particle Image Velocimetry (PIV) technique. Upper bound analysis of the current study shows that the maximum vertical displacement of the sand under the ultimate load increases for an increase in the width of footing, but at a decreasing rate with relative density of sand, whereas the relative vertical displacement in the sand decreases for an increase in the width of the footing. A well agreement is observed between experimental results for different footing widths and relative densities. The experimental analyses have shown that there exists pronounced scale effect for strip surface footing. The bearing capacity factors Nγ rapidly decrease up to footing widths B=0.25 m, 0.35 m, and 0.65 m for loose, medium-dense and dense sand respectively, after that there is no significant decrease in Nγ. The deformation modes of the soil as well as the ultimate bearing capacity values have been affected by the footing widths. The obtained results could be used to improve settlement calculation of the foundation interacting with granular soil., {"references":["H. M. Jaeger, S. R Nagel, and R. P. Behringer, \"Granular solids, liquids, and gases,\" Rev. Mod. Phys., vol. 68, no.4, pp. 1259-1273, Oct.1996.","S. J. Antony, \"Link between single-particle properties and macroscopic properties in particulate assemblies: role of structures within structures,\" Phil. Trans. R. Soc. A, vol. 365, pp. 2879-2891, Sept. 2007.","J. E. Bowles, Foundation Analysis and Design, 5th ed., McGraw-Hill, Singapore, 1997.","B. M. Das, Shallow foundations: bearing capacity and settlement, 2nd ed., CRC Press, London, 2009.","C. Liu, and J. B. Evett, Soils and foundations, 6th ed., Pearson Prentice Hall, New Jersey, 2004.","K. Terzaghi, and R. B Peck, Soil mechanics in engineering practice. Wiley, London, 1967.","S. Hansbo, Foundation engineering, Elsevier, London, 1994.","J. H. Schmertmann, P. R. Brown, and J. P. Hartman, \"Improved strain influence factor diagrams,\" J. Geo. Eng. Div. vol. 104, no. GT8, pp. 1131-1135, Aug.1978.","R. J. Adrian, \"Particle-imaging techniques for experimental fluid mechanics\", Annu. Rev. Fluid Mech., vol. 23, pp. 261-304, 1991.\n[10]\tC. O'Loughlin, and B. Lehane, \"Nonlinear cone penetration test-based method for predicting footing settlements on sand, \"J. Geo. Geoenviron. Eng., vol. 136, no.3, pp. 409-416, Aug. 2010.\n[11]\tT.G. Murthy, E. Gnanamanickam, and S. Chandrasekar, \"Deformation field in indentation of a granular ensemble,\" Phys. Rev. E, vol. 85, no.061306, pp.1-11, June 2012.\n[12]\tZ. K. Jahanger, S. J., Antony, J., Richter, \"Displacement patterns beneath a rigid beam indenting on layered soil,\" in Pro. 8th Amer. Reg. Conf. Inter. Soc. Terrain-Vehicle Sys. Michigan, 2016, Paper No.67.\n[13]\tS. Albaraki, and S.J. Antony, \"How does internal angle of hoppers affect granular flow? Experimental studies using Digital Particle Image Velocimetry,\"Pow. Techn., vol. 268, pp. 253-260, Aug. 2014.\n[14]\tASTM, American Society for Testing and Materials, Soil and Rock, Building, Stores, Geotextiles, ASTM Standard, vol. 04.08, 1989.\n[15]\tK. Head, Manual of Soil Laboratory Test. Volume 1: soil Classification and Compaction Tests, 3rd ed., CRC Press, Boca Raton, FL, 2006.\n[16]\tJ. Liu, and M. Iskander, \"Adaptive cross correlation for imaging displacements in soils,\" J. Comput. Civil Eng., vol. 18, no.1, pp. 46-57, Jan. 2004.\n[17]\tJ. Dijkstra, D. J. White, and C. Gaudin, \"Comparison of failure modes below footings on carbonate and silica sands,\" Int. J. Phys. Model. Geotech., vol. 13, no. 1, pp. 1-12, Aug. 2013.\n[18]\tC.K. Lau, Scale effects in tests on footings, PhD thesis, University of Cambridge, UK, 1988.\n[19]\tA. Altaee, and B. H. Fellenius, \"Physical modeling in sand,\" Can. Geotech. J. vol. 31, no. 3, pp. 420-431, Feb. 1994.\n[20]\tG.P. Raymond, and F.E. Komos, \"Repeated load testing of a model plane strain footing,\" Can. Geotech. J., vol.15, no. 2, pp. 190-201, Nov. 1978.\n[21]\tD. White, and M. Bolton, \"Displacement and strain paths during plane-strain model pile installation in sand,\" Géotechnique, vol. 54, no.6, pp. 375-397, Apr. 2004.\n[22]\tJ. Kumar, and M.K Bhoi, \"Interference of two closely spaced strip footings on sand using model tests,\" J. Geo. Geoenviron. Eng., vol. 135, no. 4, pp. 595-604, Apr.2009.\n[23]\tE. Hamm, F. Tapia, and F. Melo, \"Dynamics of shear bands in a dense granular material forced by a slowly moving rigid body,\" Phys. Rev. E, vol. 84, no.041304, pp. 1-7, Oct. 2011.\n[24]\tS. O. Akbas, and F. H. Kulhawy, \"Axial compression of footings in cohesionless soils. I: Load-settlement behavior,\" J. Geotech. Geoenviron. Eng., vol.135, no.11. pp. 1562-1574, Nov. 2009.\n[25]\tA.S. Vesic, \"Analysis of ultimate loads of shallow foundations,\" Soil Mech. and Found. Div., ASCE, vol.99, no. SM1, pp. 45-73, Jan.1973.\n[26]\tE. E. De Beer, Bearing capacity and settlement of shallow foundations on sand, in Proc. of Symp. Bearing Capacity and Settlement of Foundation, Duke University, Durham, N.C. 1965, pp. 15-33.\n[27]\tJ. Lee, J. Eun, M. Prezzi, and R. Salgado, \"Strain influence diagrams for settlement estimation of both isolated and multiple footings in sand,\" J. Geo. Geoenviron. Eng., vol.134, no.4, pp. 417-427, Apr. 2008.\n[28]\tP. W. Mayne, and H. G. Poulos, \"Approximate displacement influence factors for elastic shallow foundations,\" J. Geo. Geoenviron. Eng., vol. 125, no. 6, pp. 453-460, June 1999.\n[29]\tW. F. Chen, Limit Analysis and Soil Plasticity, J. Ross Publishing, Fort Lauderdale, USA, 2008.\n[30]\tA. J. Lutenegger, and D. J. DeGroot, \"Settlement of shallow foundations on granular soils,\" University of Massachusetts Transportation Center, Amherst, MA 01003, report no. 6332, June 1995.\n[31]\tL. Bjerrum, A. and Eggestad, \"Interpretation of loading test on sand.\" in Proc. of European Conf. in Soil Mechanics, Weisbaden, West Germany 1, 1963, pp.199-203.\n[32]\tK. Terzaghi, R. B. Peck, and G., Mesri, Soil mechanics in engineering practice, third ed., John Wiley and sons, New York, 1996.\n[33]\tN. F. Ismael, and A. H. N. Ahmad, \"Bearing capacity of footings on calcareous sands.\" Soils Found., vol. 30, no. 3, pp.81-90, Sept. 1990. \n[34]\tB. Cerato, and A. J. Lutenegger, \"Scale effects of shallow foundation bearing capacity on granular material,\" J. Geotech. Geoenviron. Eng., vol. 133, no. 10, pp.1192-1202, Oct. 2007."]}

Non convex granular media are involved in many industrial processes as, e.g., particle calcination/drying in rotating drums or solid catalyst particles in chemical reactors. In the case of optimizing the shape of catalysts, the experimental discrimination of new shapes based on packing density and pressure drop proved to be difficult due to the limited control of size distribution and loading procedure. There is therefore a strong interest in developing numerical tools to predict the dynamics of granular media made of particles of arbitrary shape and to simulate the flow of a fluid (either liquid or gas) around these particles. Non-convex particles are even more challenging than convex particles due to the potential multiplicity of contact points between two solid bodies. In this work, we implement new numerical strategies in our home made high-fidelity parallel numerical tools: Grains3D for granular dynamics of solid particles and PeliGRIFF for reactive fluid/solid flows. The first part of this work consists in extending the modelling capabilities of Grains3D from convex to non-convex particles based on the decomposition of a non-convex shape into a set of convex particles. We validate our numerical model with existing analytical solutions and experimental data on a rotating drum filled with 2D cross particle shapes. We also use Grains3D to study the loading of semi-periodic small size reactors with trilobic and quadralobic particles. The second part of this work consists in extending the modelling capabilities of PeliGRIFF to handle poly-lobed (and hence non-convex) particles. Our Particle Resolved Simulation (PRS) method is based on a Distributed Lagrange Multiplier / Fictitious Domain (DLM/FD) formulation combined with a Finite Volume / Staggered Grid (FV/SG) discretization scheme. Due to the lack of analytical solutions and experimental data, we assess the accuracy of our PRS method by examining the space convergence of the computed solution in assorted flow configurations such as the flow through a periodic array of poly-lobed particles and the flow in a small size packed bed reactor. Our simulation results are overall consistent with previous experimental work.; Cette thèse porte sur l'étude numérique des écoulements fluide-particules rencontrés dans l'industrie. Ces travaux se situent dans le cadre de la compréhension des phénomènes qui se déroulent dans des tambours tournants et réacteurs à lit fixe en présence de particules de forme non convexe. En effet, la forme des particules influence de manière importante la dynamique de ces milieux. A cet effet, nous nous sommes servis de la plateforme numérique parallèle Grans3D pour la dynamique des milieux granulaires et PeliGRIFF pour les écoulements multiphasiques. Dans la première partie de cette thèse, nous avons développé une nouvelle stratégie numérique qui permet de prendre en compte des particules de forme arbitrairement non convexe dans le solveur Grains3D. Elle consiste à décomposer une forme non convexe en plusieurs formes convexes quelconques. Nous avons nommé cette méthode « glued-convex ». Le modèle a été validé avec succès sur des résultats théoriques et expérimentaux de tambours tournants en présence de particules en forme de croix. Nous avons aussi utilisé le modèle pour simuler le chargement de réacteurs à lits fixes puis des lois de corrélation sur les taux de vide ont été déduites de nos résultats numériques. Dans ces travaux, nous avons aussi testé les performances parallèles de nos outils sur des simulations numériques à grande échelle de divers systèmes de particules convexes. La deuxième partie de cette thèse a été consacrée à l'extension du solveur PeliGRIFF à pouvoir prendre en compte la présence de particules multilobées (non convexes) dans des écoulements monophasiques. Une approche du type Simulation Numérique Directe, basée sur les Multiplicateurs de Lagrange Distribués / Domaine Fictif (DLM/FD), a alors été adoptée pour résoudre l'écoulement autour des particules. Une série d'études de convergence spatiale a été faite basée sur diverses configurations et divers régimes. Enfin, ces outils ont été utilisés pour simuler des écoulements au travers de lits fixes de particules de forme multi-lobée dans le but d'étudier l'influence de la forme des particules sur l'hydrodynamique dans ces lits. Les résultats ont montré une consistance avec les résultats expérimentaux disponibles dans la littérature.

Using a granular-mechanics code, we study the impact of a highly porous granular body on a hard wall. The projectile consists of monodisperse adhesive micrometer-sized silica grains. For the impact velocities studied, v

Shear deformation of granular media leads to continual restructuring of particle contact network and mechanical interactions. These changes to the mechanical state include jamming of grains, collisions, and frictional slip of particles—all of which present abrupt perturbations of internal forces and release of strain energy. Such energy release events typically result in the generation of elastic waves in the kHz frequency range, termed acoustic emissions (AE). The close association between grain-scale mechanics and AE generation motivated the use of AE as surrogate observations to assess the mechanical state of complex materials and granular flows. The study characterizes AE generation mechanisms stemming from grain-scale mechanical interactions. Basic mechanisms are considered, including frictional slip between particles, and mechanical excitation of particle configurations during force network restructuring events. The intrinsic frequencies and energy content of generated AEs bear the signature of source mechanisms and of structural features of the grain network. Acoustic measurements in simple shear experiments of glass beads reveal distinct characteristics of AE associated with different source mechanisms. These findings offer new capabilities for non-invasive interrogation of micromechancial interactions and linkage to a stochastic model of shear zone mechanics. Certain statistical features of restructuring events and associated energy release during shearing were predicted with a conceptual fiber-bundle model (FBM). In the FBM the collective behavior of a large number of basic mechanical elements (representing e.g. grain contacts), termed fibers, reproduces the reaction of disordered materials to progressive loading. The failure of fibers at an individual threshold force corresponds to slipping of a particle contact or a single rearrangement event of the granular network. The energy release from model fiber breakage is the equivalent to elastic energy from abrupt grain rearrangement events and provides an estimate of the energy available for elastic wave generation. The coupled FBM–AE model was in reasonable agreement with direct shear experiments that were performed on large granular assemblies. The results underline the potential of using AE as a diagnostic tool to study micro-mechanical interactions, shear failure and mobilization in granular material., Granular Matter, 16 (5), ISSN:1434-5021, ISSN:1434-7636

International audience; The physical laws governing the morphogenesis of biological tissues remain largely misunderstood. In particular, the role of the mechanical interactions occurring in this process needs to be better understood and studied. Inner follicular cells surrounding the oocytes of Ciona intestinalis form an epithelial monolayer resulting from an accretion process (without mitosis or apoptosis). This epithelium is elementary and useful for morphogenesis studies: the cells exhibit polygon packing with a specific but non-systematically repeatable topology (i.e. the distribution of pentagons, hexagons and heptagons changes). To understand the role of mechanical forces in tissue formation, we propose an innovative ''2D spherical'' model based on the physics of divided media. This approach simulates the cellular mechanical behavior and epithelium structuration by allowing cells to adopt a large variety of shapes and to self-organize in response to mechanical interactions. The numerical parameters considered in the model are derived from experimental data in order to perform pertinent and realistic simulations. The results obtained are compared to biological observations using the same counting method to characterize epithelium topology. Numerical and experimental data appear close enough to validate the model. It is then used for exploratory studies dealing with ''Tissue Development Speed'' variation, which is not easily attainable by experimentation. We show that the formation speed of the tissue influences its topology and hence its packing organization.

We describe here how to design a triaxial shear setup suitable for both stress and structure measurements on granular packings. Our novel setup makes use of refractive index matched scanning techniques. It can thus be used to study any kind of material compatible with that method. We describe in particular how to use hydrogel spheres, which we have found to have properties which make them particularly well suited for stress and structure measurements of granular packings. We discuss here some mechanical properties of hydrogel spheres and methods to use them in refractive index matched scanning devices such as the triaxial shear setup discussed.

We have designed a novel experimental apparatus to probe the mechanics of sheared quasi two-dimensional frictional granular materials with tunable friction from the supporting base of the apparatus. The experiment consists of a floating layer of photoelastic disks, which is subject to deformation. Forces on the particles are measured at the particle scale, using their photoelastic properties. This novel setup makes the study of the role of basal friction on sheared granular media possible.

Complex grain-boundary structures, such as the 1-2 nm thick intergranular glassy films (IGF), play a prominent role in the failure behavior of nano-phased ceramics. The IGF plays the role of an imperfection and serves as the location of strain localization and failure. We have recently performed “theoretical” mechanical loading experiments on very large atomic models of IGF in silicon nitride using ab initio simulation to obtain their failure behavior. The ab initio simulations yield characteristic post-peak softening accompanied by strain localization zone. In this paper a micro-structural granular mechanics based higher-order continuum theory is applied to model the failure behavior of these types of material systems. The results obtained from the ab initio simulations are compared with those predicted by the higher order continuum theory.

Polymer processing entails several operations involving the handling of granular materials. We present here a review of some recent advances in the field which are of relevance to polymer processing. The stress distribution in static material is considered with application to particulates in bins and hoppers. A new model for description of the rheology of dense granular flows is reviewed. A model for surface granular flows with application to rotating cylinders is discussed. Finally the mixing and segregation is tumbling mixers is reviewed.