Computing elastic moduli on 3-D X-ray computed tomography image stacks

Abstract A numerical task of current interest is to compute the effective elastic properties of a random composite material by operating on a 3D digital image of its microstructure obtained via X-ray computed tomography (CT). The 3-D image is usually sub-sampled since an X-ray CT image is typically of order 10003 voxels or larger, which is considered to be a very large finite element problem. Two main questions for the validity of any such study are then: can the sub-sample size be made sufficiently large to capture enough of the important details of the random microstructure so that the computed moduli can be thought of as accurate, and what boundary conditions should be chosen for these sub-samples? This paper contributes to the answer of both questions by studying a simulated X-ray CT cylindrical microstructure with three phases, cut from a random model system with known elastic properties. A new hybrid numerical method is introduced, which makes use of finite element solutions coupled with exact solutions for elastic moduli of square arrays of parallel cylindrical fibers. The new method allows, in principle, all of the microstructural data to be used when the X-ray CT image is in the form of a cylinder, which is often the case. Appendix A describes a similar algorithm for spherical sub-samples, which may be of use when examining the mechanical properties of particles. Cubic sub-samples are also taken from this simulated X-ray CT structure to investigate the effect of two different kinds of boundary conditions: forced periodic and fixed displacements. It is found that using forced periodic displacements on the non-geometrically periodic cubic sub-samples always gave more accurate results than using fixed displacements, although with about the same precision. The larger the cubic sub-sample, the more accurate and precise was the elastic computation, and using the complete cylindrical sample with the new method gave still more accurate and precise results. Fortran 90 programs for the analytical solutions are made available on-line, along with the parallel finite element codes used.

[1]  J. van Dommelen,et al.  X-Ray Computed Tomography Based Modelling of Polymeric Foams , 2010 .

[2]  J. Watt,et al.  Hashin-Shtrikman bounds on the effective elastic moduli of polycrystals with orthorhombic symmetry , 1979 .

[3]  Julián Bravo-Castillero,et al.  Closed-form expressions for the effective coefficients of a fiber-reinforced composite with transversely isotropic constituents – I. Elastic and square symmetry , 2001 .

[4]  Michael Ortiz,et al.  Nanomechanics of Defects in Solids , 1998 .

[5]  Stephen A. Langer,et al.  OOF: an image-based finite-element analysis of material microstructures , 2001, Comput. Sci. Eng..

[6]  S. Torquato,et al.  Photographic granularity : mathematical formulation and effect of impenetrability of grains , 1990 .

[7]  Guohua Cao,et al.  Compressive sampling based interior reconstruction for dynamic carbon nanotube micro-CT. , 2009, Journal of X-ray science and technology.

[8]  M. S. Bartlett,et al.  The expected number of clumps when convex laminae are placed at random and with random orientation on a plane area , 1954, Mathematical Proceedings of the Cambridge Philosophical Society.

[9]  Thorpe,et al.  Percolation properties of random ellipses. , 1988, Physical review. A, General physics.

[10]  G. Milton.  The Theory of Composites , 2002 .

[11]  A. Ortona,et al.  Evaluation of a simple finite element method for the calculation of effective electrical conductivity of compression moulded polymer–graphite composites , 2013 .

[12]  L. Rayleigh,et al.  LVI. On the influence of obstacles arranged in rectangular order upon the properties of a medium , 1892 .

[13]  Matthieu Faessel,et al.  3D Modelling of random cellulosic fibrous networks based on X-ray tomography and image analysis , 2005 .

[14]  E. Garboczi,et al.  User manual for finite element and finite difference programs:: a parallel version of NIST IR 6269 , 2003 .

[15]  E. Garboczi,et al.  Intrinsic Viscosity and the Polarizability of Particles Having a Wide Range of Shapes , 2007 .

[16]  V. Kushch.  Multipole Expansion Method in Micromechanics of Composites , 2013 .

[17]  Day,et al.  Universal conductivity curve for a plane containing random holes. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[18]  David R. McKenzie,et al.  Transport properties of regular arrays of cylinders , 1979, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[19]  Leslie Greengard,et al.  A renormalization method for the evaluation of lattice sums , 1994 .

[20]  Edward J. Garboczi,et al.  Institute of Physics Publishing Modelling and Simulation in Materials Science and Engineering Linear Elastic Properties of 2d and 3d Models of Porous Materials Made from Elongated Objects , 2022 .

[21]  Edward J. Garboczi,et al.  Elastic Properties of Model Porous Ceramics , 2000, cond-mat/0006334.

[22]  E. Maire,et al.  Microstructural analysis of alumina chromium composites by X-ray tomography and 3-D finite element simulation of thermal stresses , 2003 .

[23]  J. Watt,et al.  Clarification of the Hashin‐Shtrikman bounds on the effective elastic moduli of polycrystals with hexagonal, trigonal, and tetragonal symmetries , 1980 .

[24]  A. Sangani,et al.  Elastic coefficients of composites containing spherical inclusions in a periodic array , 1987 .

[25]  J. Quintanilla,et al.  Local volume fraction fluctuations in random media , 1997 .

[26]  L. Mishnaevsky,et al.  Meso cell model of fiber reinforced composite: Interface stress statistics and debonding paths , 2008 .

[27]  Edward J. Garboczi,et al.  An algorithm for computing the effective linear elastic properties of heterogeneous materials: Three-dimensional results for composites with equal phase poisson ratios , 1995 .

[28]  N. Nicorovici,et al.  Green's tensors and lattice sums for electrostatics and elastodynamics , 1997, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[29]  R. Hill.  Elastic properties of reinforced solids: some theoretical principles , 1963 .

[30]  A. P. Roberts,et al.  Computation of the linear elastic properties of random porous materials with a wide variety of microstructure , 2002, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[31]  S. Torquato,et al.  Random Heterogeneous Materials: Microstructure and Macroscopic Properties , 2005 .

[32]  M Stampanoni,et al.  Finite element 3D reconstruction of the pulmonary acinus imaged by synchrotron X-ray tomography. , 2008, Journal of applied physiology.

[33]  V. Kushch.  Computation of the effective elastic moduli of a granular composite material of regular structure , 1987 .

[34]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[35]  S. Shtrikman,et al.  A variational approach to the theory of the elastic behaviour of multiphase materials , 1963 .

[36]  E. Garboczi,et al.  Finite Element and Finite Difference Programs for Computing the Linear Electric and Elastic Properties of Digital Images of Random Materials | NIST , 1998 .

[37]  John C. Russ,et al.  The Image Processing Handbook , 2016, Microscopy and Microanalysis.

[38]  E. Garboczi,et al.  Elastic moduli of composites containing a low concentration of complex-shaped particles having a general property contrast with the matrix , 2012 .

[39]  Edward J. Garboczi,et al.  A hybrid finite element-analytical method for determining the intrinsic elastic moduli of particles having moderately extended shapes and a wide range of elastic properties , 2006 .

[40]  S. L. Crouch,et al.  Elastic fields and effective moduli of particulate nanocomposites with the Gurtin–Murdoch model of interfaces , 2013 .