Multi-grid reduced-order topology optimization

Additive printing allows the “single step” production of virtually any complex mechanical component. However, the manufacturing process involves a layer-by-layer deposition of material, which leads to an anisotropic mechanical behavior of the whole component. This would then entail a very fine 3D model to simulate the mechanical performance accurately. This simulation would also need to be integrated within an iterative design process in order to obtain the most efficient design. Both reasons explain the prohibitive number of calculations needed, which is currently beyond the capacities of existing software and computers. Recent research papers have opened promising pathways for integrating model reduction techniques within the overall topology optimization process. However, these approaches still present challenges such as choosing the minimum number and appropriate selection of the snapshots required to get accurate simulations. In this work, we present a methodology in the combined field of reduced-order modeling and topology optimization. The key idea consists of projecting the higher dimensional system of equations onto a lower dimensional space with the reduced basis vectors constructed using Proper Orthogonal Decomposition (POD). This reduced basis is updated in an incremental “on-the-fly” manner using alternatively costly high-fidelity and more rapid lower fidelity simulation snapshots. The variable-fidelity resolutions of successive approximations to the global system of equations are then integrated into the topology optimization process. The approaches are tested and computational savings and precision are compared, using both minimum compliance and compliant mechanism design benchmark problems.

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

[2]  Charbel Farhat,et al.  Nonlinear model order reduction based on local reduced‐order bases , 2012 .

[3]  Marian Brezina,et al.  Convergence of algebraic multigrid based on smoothed aggregation , 1998, Numerische Mathematik.

[4]  Krister Svanberg,et al.  A Class of Globally Convergent Optimization Methods Based on Conservative Convex Separable Approximations , 2002, SIAM J. Optim..

[5]  J. Marsden,et al.  Dimensional model reduction in non‐linear finite element dynamics of solids and structures , 2001 .

[6]  K. Svanberg The method of moving asymptotes—a new method for structural optimization , 1987 .

[7]  Ole Sigmund,et al.  Approximate reanalysis in topology optimization , 2009 .

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

[9]  Weihong Zhang,et al.  Extended Co-Kriging interpolation method based on multi-fidelity data , 2018, Appl. Math. Comput..

[10]  P Kerfriden,et al.  Bridging Proper Orthogonal Decomposition methods and augmented Newton-Krylov algorithms: an adaptive model order reduction for highly nonlinear mechanical problems. , 2011, Computer methods in applied mechanics and engineering.

[11]  David G. MacManus,et al.  A variable‐fidelity aerodynamic model using proper orthogonal decomposition , 2016 .

[12]  Wolfgang Hackbusch,et al.  Multi-grid methods and applications , 1985, Springer series in computational mathematics.

[13]  Piotr Breitkopf,et al.  Constrained Proper Orthogonal Decomposition based on QR-factorization for aerodynamical shape optimization , 2013, Appl. Math. Comput..

[14]  Ramana V. Grandhi,et al.  A survey of structural and multidisciplinary continuum topology optimization: post 2000 , 2014 .

[15]  Piotr Breitkopf,et al.  Model reduction for multidisciplinary optimization - application to a 2D wing , 2008 .

[16]  P. Villon,et al.  Nonlinear Shape-Manifold Learning Approach: Concepts, Tools and Applications , 2016, Archives of Computational Methods in Engineering.

[17]  G. Yoon Structural topology optimization for frequency response problem using model reduction schemes , 2010 .

[18]  Raúl Tempone,et al.  Galerkin Finite Element Approximations of Stochastic Elliptic Partial Differential Equations , 2004, SIAM J. Numer. Anal..

[19]  Habib N. Najm,et al.  A multigrid solver for two-dimensional stochastic diffusion equations , 2003 .

[20]  P. Breitkopf,et al.  Evolutionary topology optimization of elastoplastic structures , 2017 .

[21]  P. Breitkopf,et al.  Design of materials using topology optimization and energy-based homogenization approach in Matlab , 2015 .

[22]  Boyan Stefanov Lazarov,et al.  On multigrid-CG for efficient topology optimization , 2014 .

[23]  Piotr Breitkopf,et al.  Topology optimization of multiscale elastoviscoplastic structures , 2016 .

[24]  Weihong Zhang,et al.  Proper orthogonal decomposition with high number of linear constraints for aerodynamical shape optimization , 2014, Appl. Math. Comput..

[25]  Piotr Breitkopf,et al.  On-the-fly model reduction for large-scale structural topology optimization using principal components analysis , 2020, Structural and Multidisciplinary Optimization.

[26]  M. Zhou,et al.  The COC algorithm, Part II: Topological, geometrical and generalized shape optimization , 1991 .

[27]  Kevin Carlberg,et al.  Krylov-Subspace Recycling via the POD-Augmented Conjugate-Gradient Method , 2015, SIAM J. Matrix Anal. Appl..

[28]  C. Gogu,et al.  Efficient surrogate construction by combining response surface methodology and reduced order modeling , 2013, Structural and Multidisciplinary Optimization.

[29]  Piotr Breitkopf,et al.  Model reduction by CPOD and Kriging , 2010 .

[30]  Piotr Breitkopf,et al.  Towards surrogate modeling of material microstructures through the processing variables , 2017, Appl. Math. Comput..

[31]  R. Haftka,et al.  Review of multi-fidelity models , 2016, Advances in Computational Science and Engineering.

[32]  K. Meerbergen,et al.  Using Krylov‐Padé model order reduction for accelerating design optimization of structures and vibrations in the frequency domain , 2012 .

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

[34]  Thomas A. Manteuffel,et al.  Smoothed aggregation algebraic multigrid for stochastic PDE problems with layered materials , 2014, Numer. Linear Algebra Appl..

[35]  E. Sturler,et al.  Large‐scale topology optimization using preconditioned Krylov subspace methods with recycling , 2007 .

[36]  Ole Sigmund,et al.  Giga-voxel computational morphogenesis for structural design , 2017, Nature.

[37]  Xiangmin Jiao,et al.  A hybrid geometric + algebraic multigrid method with semi-iterative smoothers , 2014, Numer. Linear Algebra Appl..

[38]  Boyan Stefanov Lazarov,et al.  Topology optimization using PETSc: An easy-to-use, fully parallel, open source topology optimization framework , 2015 .

[39]  Christian Gogu,et al.  Improving the efficiency of large scale topology optimization through on‐the‐fly reduced order model construction , 2015 .

[40]  Piotr Breitkopf,et al.  Towards a space reduction approach for efficient structural shape optimization , 2013 .

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

[42]  Mary Frecker,et al.  Topology optimization of 2D continua for minimum compliance using parallel computing , 2006 .

[43]  James K. Guest,et al.  Reducing dimensionality in topology optimization using adaptive design variable fields , 2010 .

[44]  O. Sigmund,et al.  Efficient use of iterative solvers in nested topology optimization , 2010 .

[45]  B. Raghavan,et al.  Combined meso-scale modeling and experimental investigation of the effect of mechanical damage on the transport properties of cementitious composites , 2016 .