How to determine composite material properties using numerical homogenization

Abstract Numerical homogenization is an efficient way to determine effective macroscopic properties, such as the elasticity tensor, of a periodic composite material. In this paper an educational description of the method is provided based on a short, self-contained Matlab implementation. It is shown how the basic code, which computes the effective elasticity tensor of a two material composite, where one material could be void, is easily extended to include more materials. Furthermore, extensions to homogenization of conductivity, thermal expansion, and fluid permeability are described in detail. The unit cell of the periodic material can take the shape of a square, rectangle, or parallelogram, allowing for all kinds of 2D periodicities.

[1]  E. Hinton,et al.  A review of homogenization and topology opimization II—analytical and numerical solution of homogenization equations , 1998 .

[2]  Thomas J. R. Hughes,et al.  The Stokes problem with various well-posed boundary conditions - Symmetric formulations that converge for all velocity/pressure spaces , 1987 .

[3]  Anders Clausen,et al.  Efficient topology optimization in MATLAB using 88 lines of code , 2011 .

[4]  E. Hinton,et al.  A review of homogenization and topology optimization III—topology optimization using optimality criteria , 1998 .

[5]  J. Drummond,et al.  Laminar viscous flow through regular arrays of parallel solid cylinders , 1984 .

[6]  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.

[7]  J. Fish,et al.  Multiscale asymptotic homogenization for multiphysics problems with multiple spatial and temporal scales: a coupled thermo-viscoelastic example problem , 2002 .

[8]  E. Hinton,et al.  A review of homogenization and topology optimization I- homogenization theory for media with periodic structure , 1998 .

[9]  O. Sigmund.  Materials with prescribed constitutive parameters: An inverse homogenization problem , 1994 .

[10]  R. Cook,et al.  Concepts and Applications of Finite Element Analysis , 1974 .

[11]  S. Torquato,et al.  Design of materials with extreme thermal expansion using a three-phase topology optimization method , 1997 .

[12]  P R Fernandes,et al.  Permeability analysis of scaffolds for bone tissue engineering. , 2012, Journal of biomechanics.

[13]  James K. Guest,et al.  Optimizing multifunctional materials: Design of microstructures for maximized stiffness and fluid permeability , 2006 .

[14]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

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

[16]  K. Seffen,et al.  International Journal of Solids and Structures , 2015 .

[17]  N. Kikuchi,et al.  Preprocessing and postprocessing for materials based on the homogenization method with adaptive fini , 1990 .

[18]  Scott J. Hollister,et al.  Tailoring the mechanical properties of 3D-designed poly(glycerol sebacate) scaffolds for cartilage applications. , 2010, Journal of biomedical materials research. Part A.

[19]  Ole Sigmund,et al.  Saturated poroelastic actuators generated by topology optimization , 2011 .

[20]  S. Hollister.  Porous scaffold design for tissue engineering , 2005, Nature materials.

[21]  Alejandro R. Diaz,et al.  Designing materials with prescribed elastic properties using polygonal cells , 2003 .

[22]  Proceedings of the Royal Society (London) , 1906, Science.

[23]  Haimin Yao,et al.  Journal of the Mechanics and Physics of Solids , 2014 .

[24]  John K. Tomfohr,et al.  Lecture Notes on Physics , 1879, Nature.

[25]  A. Bensoussan,et al.  Asymptotic analysis for periodic structures , 1979 .