Data-driven topology optimization of spinodoid metamaterials with seamlessly tunable anisotropy

We present a two-scale topology optimization framework for the design of macroscopic bodies with an optimized elastic response, which is achieved by means of a spatially-variant cellular architecture on the microscale. The chosen spinodoid topology for the cellular network on the microscale (which is inspired by natural microstructures forming during spinodal decomposition) admits a seamless spatial grading as well as tunable elastic anisotropy, and it is parametrized by a small set of design parameters associated with the underlying Gaussian random field. The macroscale boundary value problem is discretized by finite elements, which in addition to the displacement field continuously interpolate the microscale design parameters. By assuming a separation of scales, the local constitutive behavior on the macroscale is identified as the homogenized elastic response of the microstructure based on the local design parameters. As a departure from classical FE2-type approaches, we replace the costly microscale homogenization by a data-driven surrogate model, using deep neural networks, which accurately and efficiently maps design parameters onto the effective elasticity tensor. The model is trained on homogenized stiffness data obtained from numerical homogenization by finite elements. As an added benefit, the machine learning setup admits automatic differentiation, so that sensitivities (required for the optimization problem) can be computed exactly and without the need for numerical derivatives – a strategy that holds promise far beyond the elastic stiffness. Therefore, this framework presents a new opportunity for multiscale topology optimization based on data-driven surrogate models.

[1]  T. E. Bruns,et al.  Topology optimization of non-linear elastic structures and compliant mechanisms , 2001 .

[2]  J. -F. M. Barthelemy,et al.  Automatic differentiation as a tool in engineering design , 1992 .

[3]  M. Bendsøe,et al.  Generating optimal topologies in structural design using a homogenization method , 1988 .

[4]  H. E. Cook,et al.  Brownian motion in spinodal decomposition , 1970 .

[5]  M. Bendsøe Optimal shape design as a material distribution problem , 1989 .

[6]  Seung Chul Han,et al.  A Microscopic Shell Structure with Schwarz’s D-Surface , 2017, Scientific Reports.

[7]  Jian Su,et al.  Automatic differentiation in robust optimization , 1996 .

[8]  Chang Liu,et al.  Machine Learning-Driven Real-Time Topology Optimization Under Moving Morphable Component-Based Framework , 2018, Journal of Applied Mechanics.

[9]  M. Bendsøe,et al.  Shape optimization of structures for multiple loading conditions using a homogenization method , 1992 .

[10]  M. Bendsøe,et al.  Topology Optimization: "Theory, Methods, And Applications" , 2011 .

[11]  L. Valdevit,et al.  The mechanical response of cellular materials with spinodal topologies , 2019, Journal of the Mechanics and Physics of Solids.

[12]  Ole Sigmund,et al.  A 99 line topology optimization code written in Matlab , 2001 .

[13]  M. Ashby,et al.  Cellular Materials in Nature and Medicine , 2010 .

[14]  D. Kochmann,et al.  Microstructural patterns with tunable mechanical anisotropy obtained by simulating anisotropic spinodal decomposition , 2018, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[15]  H. Rodrigues,et al.  Hierarchical optimization of material and structure , 2002 .

[16]  Michael Stingl,et al.  Efficient Two-Scale Optimization of Manufacturable Graded Structures , 2012, SIAM J. Sci. Comput..

[17]  Howon Lee,et al.  Ultralight, ultrastiff mechanical metamaterials , 2014, Science.

[18]  Andrej Cherkaev,et al.  Microstructures of Composites of Extremal Rigidity and Exact Bounds on the Associated Energy Density , 1997 .

[19]  A. Karma,et al.  Evolution of nanoporosity in dealloying , 2001, Nature.

[20]  F. Zok,et al.  Effects of nodal fillets and external boundaries on compressive response of an octet truss , 2018 .

[21]  O. Sigmund,et al.  Topology optimization approaches , 2013, Structural and Multidisciplinary Optimization.

[22]  Julia R. Greer,et al.  Reexamining the mechanical property space of three-dimensional lattice architectures , 2017 .

[23]  O. Sigmund Tailoring materials with prescribed elastic properties , 1995 .

[24]  Juliane Junker Topics In The Mathematical Modelling Of Composite Materials , 2016 .

[25]  Daniel A. White,et al.  Simple, accurate surrogate models of the elastic response of three-dimensional open truss micro-architectures with applications to multiscale topology design , 2019, Structural and Multidisciplinary Optimization.

[26]  Vikram Deshpande,et al.  Concepts for enhanced energy absorption using hollow micro-lattices , 2010 .

[27]  C. M. Portela,et al.  Impact of node geometry on the effective stiffness of non-slender three-dimensional truss lattice architectures , 2018, Extreme Mechanics Letters.

[28]  J. E. Hilliard,et al.  Free Energy of a Nonuniform System. I. Interfacial Free Energy , 1958 .

[29]  Andreas Griewank,et al.  Evaluating derivatives - principles and techniques of algorithmic differentiation, Second Edition , 2000, Frontiers in applied mathematics.

[30]  N. Kikuchi,et al.  General topology optimization method with continuous and discrete orientation design using isoparametric projection , 2015 .

[31]  Axel Voigt,et al.  Faceting of Equilibrium and Metastable Nanostructures: A Phase-Field Model of Surface Diffusion Tackling Realistic Shapes , 2015 .

[32]  Abhinav Sharma,et al.  Effect of anisotropic diffusion on spinodal decomposition. , 2018, Soft matter.

[33]  Levent Burak Kara,et al.  A data-driven investigation and estimation of optimal topologies under variable loading configurations , 2014, Comput. methods Biomech. Biomed. Eng. Imaging Vis..

[34]  Zhu Jiejiang,et al.  Topology Optimization Design of 3D Continuum Structure with Reserved Hole Based on Variable Density Method , 2016 .

[35]  Axel Voigt,et al.  A new phase-field model for strongly anisotropic systems , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[36]  D. Kochmann,et al.  An assessment of numerical techniques to find energy‐minimizing microstructures associated with nonconvex potentials , 2019, International Journal for Numerical Methods in Engineering.

[37]  Li Zheng,et al.  Inverse-designed spinodoid metamaterials , 2020, npj Computational Materials.

[38]  Ole Sigmund,et al.  Homogenization-based stiffness optimization and projection of 2D coated structures with orthotropic infill , 2019, Computer Methods in Applied Mechanics and Engineering.

[39]  E. Favvas,et al.  What is spinodal decomposition , 2008 .

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

[41]  M. Bendsøe,et al.  Material interpolation schemes in topology optimization , 1999 .

[42]  A. Hamza,et al.  Scaling equation for yield strength of nanoporous open-cell foams , 2006 .

[43]  Daniel A. White,et al.  Multiscale topology optimization using neural network surrogate models , 2019, Computer Methods in Applied Mechanics and Engineering.

[44]  R. Adler,et al.  Random Fields and Geometry , 2007 .

[45]  P. Pedersen On optimal orientation of orthotropic materials , 1989 .

[46]  In Gwun Jang,et al.  Deep learning for determining a near-optimal topological design without any iteration , 2018, Structural and Multidisciplinary Optimization.

[47]  Peter D. Dunning,et al.  Simultaneous material and structural optimization by multiscale topology optimization , 2016 .

[48]  J. Greer,et al.  Strong, lightweight, and recoverable three-dimensional ceramic nanolattices , 2014, Science.

[49]  Erik Lund,et al.  Discrete material optimization of general composite shell structures , 2005 .

[50]  Yang Gao,et al.  Computer Methods in Biomechanics and Biomedical Engineering: Imaging & Visualization , 2015 .

[51]  E. Thomas,et al.  Micro‐/Nanostructured Mechanical Metamaterials , 2012, Advanced materials.

[52]  Martin P. Bendsøe,et al.  Free material optimization via mathematical programming , 1997, Math. Program..

[53]  Stewart,et al.  Spinodal decomposition of a crystal surface. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[54]  Douglas E. Smith,et al.  Continuous Fiber Angle Topology Optimization for Polymer Composite Deposition Additive Manufacturing Applications , 2019, Fibers.

[55]  Pauli Pedersen,et al.  On thickness and orientational design with orthotropic materials , 1991 .

[56]  Bruder,et al.  Spinodal decomposition in thin films of a polymer blend. , 1992, Physical review letters.

[57]  C. M. Portela,et al.  Extreme mechanical resilience of self-assembled nanolabyrinthine materials , 2020, Proceedings of the National Academy of Sciences.

[58]  B. D. Nguyen,et al.  Design of the P-surfaced shellular, an ultra-low density material with micro-architecture , 2017 .

[59]  Zafer Gürdal,et al.  Combined topology and fiber path design of composite layers using cellular automata , 2005 .

[60]  John W. Cahn,et al.  Phase Separation by Spinodal Decomposition in Isotropic Systems , 1965 .

[61]  P. Breitkopf,et al.  Concurrent topology optimization design of material and structure within FE2 nonlinear multiscale analysis framework , 2014 .

[62]  Dirk Mohr,et al.  3D Plate‐Lattices: An Emerging Class of Low‐Density Metamaterial Exhibiting Optimal Isotropic Stiffness , 2018, Advanced materials.

[63]  Tong Gao,et al.  A bi‐value coding parameterization scheme for the discrete optimal orientation design of the composite laminate , 2012 .

[64]  Ivan V. Oseledets,et al.  Neural networks for topology optimization , 2017, Russian Journal of Numerical Analysis and Mathematical Modelling.

[65]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[66]  O. Sigmund Morphology-based black and white filters for topology optimization , 2007 .

[67]  Jun Wu,et al.  Design and Optimization of Conforming Lattice Structures , 2019, IEEE Transactions on Visualization and Computer Graphics.

[68]  Helder C. Rodrigues,et al.  A hierarchical model for concurrent material and topology optimisation of three-dimensional structures , 2008 .

[69]  Yan Zhang,et al.  Concurrent topology optimization for cellular structures with nonuniform microstructures based on the kriging metamodel , 2018, Structural and Multidisciplinary Optimization.

[70]  N. Kikuchi,et al.  A homogenization method for shape and topology optimization , 1991 .

[71]  Glaucio H. Paulino,et al.  Bridging topology optimization and additive manufacturing , 2015, Structural and Multidisciplinary Optimization.

[72]  Julien Yvonnet,et al.  Numerically explicit potentials for the homogenization of nonlinear elastic heterogeneous materials , 2009 .

[73]  T. Shi,et al.  Optimization of composite structures with continuous spatial variation of fiber angle through Shepard interpolation , 2017 .

[74]  P. Breitkopf,et al.  Multiscale structural topology optimization with an approximate constitutive model for local material microstructure , 2015 .

[75]  A. Mohraz,et al.  Bicontinuous Macroporous Materials from Bijel Templates , 2010, Advanced materials.

[76]  H. Wadley,et al.  Mechanical metamaterials at the theoretical limit of isotropic elastic stiffness , 2017, Nature.

[77]  Yong-Wei Zhang,et al.  Discrete‐Continuum Duality of Architected Materials: Failure, Flaws, and Fracture , 2018, Advanced Functional Materials.

[78]  Jimmy Ba,et al.  Adam: A Method for Stochastic Optimization , 2014, ICLR.

[79]  Isabelle Charpentier,et al.  On higher-order differentiation in nonlinear mechanics , 2012, Optim. Methods Softw..

[80]  Ole Sigmund,et al.  On the Design of Compliant Mechanisms Using Topology Optimization , 1997 .

[81]  J. Petersson,et al.  Numerical instabilities in topology optimization: A survey on procedures dealing with checkerboards, mesh-dependencies and local minima , 1998 .

[82]  P. Pedersen Bounds on elastic energy in solids of orthotropic materials , 1990 .

[83]  Tam H. Nguyen,et al.  A computational paradigm for multiresolution topology optimization (MTOP) , 2010 .

[84]  Rashid K. Abu Al-Rub,et al.  Multifunctional Mechanical Metamaterials Based on Triply Periodic Minimal Surface Lattices , 2019, Advanced Engineering Materials.

[85]  Steve Marschner,et al.  Microstructures to control elasticity in 3D printing , 2015, ACM Trans. Graph..

[86]  Liang Gao,et al.  Topology optimization for multiscale design of porous composites with multi-domain microstructures , 2019, Computer Methods in Applied Mechanics and Engineering.