Research Papers

A 149 Line Homogenization Code for Three-Dimensional Cellular Materials Written in matlab

[+] Author and Article Information
Guoying Dong

Department of Mechanical Engineering,
McGill University,
817 Rue Sherbrooke Ouest, Room G53,
Montréal, QC H3A 0C3, Canada
e-mail: guoying.dong@mail.mcgill.ca

Yunlong Tang

Department of Mechanical Engineering,
McGill University,
817 Rue Sherbrooke Ouest, Room G53,
Montréal, QC H3A 0C3, Canada
e-mail: tang.yunlong@mail.mcgill.ca

Yaoyao Fiona Zhao

Department of Mechanical Engineering,
McGill University,
817 Rue Sherbrooke Ouest, Room 148,
Montréal, QC H3A 0C3, Canada
e-mail: yaoyao.zhao@mcgill.ca

1Corresponding author.

Contributed by the Materials Division of ASME for publication in the JOURNAL OF ENGINEERING MATERIALS AND TECHNOLOGY. Manuscript received October 25, 2017; final manuscript received May 24, 2018; published online July 10, 2018. Assoc. Editor: Toshio Nakamura.

J. Eng. Mater. Technol 141(1), 011005 (Jul 10, 2018) (11 pages) Paper No: MATS-17-1316; doi: 10.1115/1.4040555 History: Received October 25, 2017; Revised May 24, 2018

Cellular architectures are promising in a variety of engineering applications due to attractive material properties. Additive manufacturing has reduced the difficulty in the fabrication of three-dimensional (3D) cellular materials. In this paper, the numerical homogenization method for 3D cellular materials is provided based on a short, self-contained matlab code. It is an educational description that shows how the homogenized constitutive matrix is computed by a voxel model with one material to be void and another material to be solid. A voxel generation algorithm is proposed to generate the voxel model easily by the wireframe scripts of unit cell topologies. The format of the wireframe script is defined so that the topology can be customized. The homogenization code is then extended to multimaterial cellular structures and thermal conductivity problems. The result of the numerical homogenization shows that different topologies exhibit anisotropic elastic properties to a different extent. It is also found that the anisotropy of cellular materials can be controlled by adjusting the combination of materials.

Copyright © 2019 by ASME
Your Session has timed out. Please sign back in to continue.


Schaedler, T. A. , and Carter, W. B. , 2016, “ Architected Cellular Materials,” Annu. Rev. Mater. Res., 46(1), pp. 187–210. [CrossRef]
Queheillalt, D. T. , and Wadley, H. N. , 2005, “ Cellular Metal Lattices With Hollow Trusses,” Acta Mater., 53(2), pp. 303–313. [CrossRef]
Clough, E. C. , Ensberg, J. , Eckel, Z. C. , Ro, C. J. , and Schaedler, T. A. , 2016, “ Mechanical Performance of Hollow Tetrahedral Truss Cores,” Int. J. Solids Struct., 91, pp. 115–126. [CrossRef]
Schaedler, T. A. , Ro, C. J. , Sorensen, A. E. , Eckel, Z. , Yang, S. S. , Carter, W. B. , and Jacobsen, A. J. , 2014, “ Designing Metallic Microlattices for Energy Absorber Applications,” Adv. Eng. Mater., 16(3), pp. 276–283. [CrossRef]
Wadley, H. N. , and Queheillalt, D. T. , 2007, “ Thermal Applications of Cellular Lattice Structures,” Mater. Sci. Forum, 539–543, pp. 242–247. [CrossRef]
Hassani, B. , and Hinton, E. , 1998, “ A Review of Homogenization and Topology Optimization I—Homogenization Theory for Media With Periodic Structure,” Comput. Struct., 69(6), pp. 707–717. [CrossRef]
Bensoussan, A. , Lions, J. L. , and Papanicolaou, G. , 2011, Asymptotic Analysis for Periodic Structures, Vol. 374, Elsevier, Providence, RI.
Torquato, S. , 2013, Random Heterogeneous Materials: Microstructure and Macroscopic Properties, Vol. 16, Springer Science & Business Media, New York.
Bendsøe, M. P. , and Kikuchi, N. , 1988, “ Generating Optimal Topologies in Structural Design Using a Homogenization Method,” Comput. Methods Appl. Mech. Eng., 71(2), pp. 197–224. [CrossRef]
Arabnejad, S. , and Pasini, D. , 2013, “ Mechanical Properties of Lattice Materials Via Asymptotic Homogenization and Comparison With Alternative Homogenization Methods,” Int. J. Mech. Sci., 77, pp. 249–262. [CrossRef]
Hassani, B. , and Hinton, E. , 1998, “ A Review of Homogenization and Topology Optimization II: Analytical and Numerical Solution of Homogenization Equations,” Comput. Struct., 69(6), pp. 719–738. [CrossRef]
Andreassen, E. , and Andreasen, C. S. , 2014, “ How to Determine Composite Material Properties Using Numerical Homogenization,” Comput. Mater. Sci., 83, pp. 488–495. [CrossRef]
Chandrupatla, T. R. , Belegundu, A. D. , Ramesh, T. , and Ray, C. , 2002, Introduction to Finite Elements Engineering, Vol. 2, Prentice Hall, Upper Saddle River, NJ.
Lautrup, B. , 2005, Physics of Continuous Matter: Exotic and Everyday Phenomena in the Macroscopic World, CRC Press, Boca Raton, FL.
Xu, S. , Shen, J. , Zhou, S. , Huang, X. , and Xie, Y. M. , 2016, “ Design of Lattice Structures With Controlled Anisotropy,” Mater. Des., 93, pp. 443–447. [CrossRef]
Vaezi, M. , Chianrabutra, S. , Mellor, B. , and Yang, S. , 2013, “ Multiple Material Additive Manufacturing—Part 1: A Review,” Virtual Phys. Prototyping, 8(1), pp. 19–50. [CrossRef]
Xu, H. , and Pasini, D. , 2016, “ Structurally Efficient Three-Dimensional Metamaterials With Controllable Thermal Expansion,” Sci. Rep., 6(1), p. 34924. [CrossRef] [PubMed]
Yu, C.-H. , Haller, K. , Ingber, D. , and Nagpal, R. , 2008, “ Morpho: A Self-Deformable Modular Robot Inspired by Cellular Structure,” International Conference on Intelligent Robots and Systems (IROS 2008), Nice, France, Sept. 22–26, pp. 3571–3578.


Grahic Jump Location
Fig. 1

Wireframe of seven topologies

Grahic Jump Location
Fig. 2

An example of the voxel model of a lattice structure with cubic topology: (a) perspective view, (b) top view of each layers, and (c) 3D logical matrix

Grahic Jump Location
Fig. 3

Two cases when calculating the minimum distance from voxel center to the line segment of lattice wireframe

Grahic Jump Location
Fig. 4

Voxel models of seven topologies

Grahic Jump Location
Fig. 5

Local node numbers and the natural coordinate of hexahedron element

Grahic Jump Location
Fig. 6

An example of hexahedron mesh used to model the unit cell: (a) the ID of each element, (b) the nonperiodic node number, (c) global node number of element 1 compared to local node number, and (d) periodic node numbers

Grahic Jump Location
Fig. 7

An example of periodic cellular materials with randomized unit cell: (a) periodic Voronoi structure, (b) one unit cell, and (c) voxel model

Grahic Jump Location
Fig. 8

Effective Young's modulus of lattice structures with 10% relative density

Grahic Jump Location
Fig. 9

Effective Young's modulus of lattice structures with 30% relative density

Grahic Jump Location
Fig. 10

Effective Young's modulus of lattice structures with 50% relative density

Grahic Jump Location
Fig. 11

The relationship between the resolution of the voxel and the axial Young's modulus

Grahic Jump Location
Fig. 12

The convergence plot for PCG method with different preconditioners

Grahic Jump Location
Fig. 13

A three-material lattice cellular structure with cubic-center topologies

Grahic Jump Location
Fig. 14

The Young's modulus surface of cubic-center lattice with (a) multimaterials and (b) single material



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In