A globally convergent method to accelerate topology optimization using on-the-fly model reduction

We present a globally convergent method to accelerate density-based topology optimization using projection-based reduced-order models (ROMs) and trust-region methods. To accelerate topology optimization, we replace the large-scale finite element simulation, which dominates the computational cost, with ROMs that reduce the cost of objective function and sensitivity evaluations by orders of magnitude. To guarantee convergence, we first introduce a trust-region method that employs generalized trust-region constraints and prove it is globally convergent. We then devise a class of globally convergent ROM-accelerated topology optimization methods informed by two theories: the aforementioned trust-region theory, which identifies the ROM accuracy conditions required to guarantee the method converges to a critical point of the original topology optimization problem; a posteriori error estimation theory for project-based ROMs, which informs ROM construction procedure to meet the accuracy conditions. This leads to trust-region methods that construct and update the ROM on-the-fly during optimization; the methods are guaranteed to converge to a critical point of the original, unreduced topology optimization problem, regardless of starting point. Numerical experiments on three different structural topology optimization problems demonstrate the proposed reduced topology optimization methods accelerate convergence to the optimal design by a factor of at least two.

[1]  J. Hesthaven,et al.  Certified Reduced Basis Methods for Parametrized Partial Differential Equations , 2015 .

[2]  Yongping Yu,et al.  The approximate reanalysis method for topology optimization under harmonic force excitations with multiple frequencies , 2017 .

[3]  David Amsallem,et al.  Fast local reduced basis updates for the efficient reduction of nonlinear systems with hyper-reduction , 2015, Advances in Computational Mathematics.

[4]  D. Rovas,et al.  A blackbox reduced-basis output bound method for noncoercive linear problems , 2002 .

[5]  Charbel Farhat,et al.  On the Use of Discrete Nonlinear Reduced-Order Models for the Prediction of Steady-State Flows Past Parametrically Deformed Complex Geometries , 2016 .

[6]  N. Nguyen,et al.  An ‘empirical interpolation’ method: application to efficient reduced-basis discretization of partial differential equations , 2004 .

[7]  Kurt Maute,et al.  Large-scale parallel topology optimization using a dual-primal substructuring solver , 2008 .

[8]  Ekkehard W. Sachs,et al.  Managing POD models by optimization methods , 2002, Proceedings of the 41st IEEE Conference on Decision and Control, 2002..

[9]  R. Carter.  On the global convergence of trust region algorithms using inexact gradient information , 1991 .

[10]  A. Patera,et al.  Reduced basis approximation and a posteriori error estimation for affinely parametrized elliptic coercive partial differential equations , 2007 .

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

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

[13]  M. Brand,et al.  Fast low-rank modifications of the thin singular value decomposition , 2006 .

[14]  H PaulinoGlaucio,et al.  Bridging topology optimization and additive manufacturing , 2016 .

[15]  L. Sirovich.  Turbulence and the dynamics of coherent structures. II. Symmetries and transformations , 1987 .

[16]  Sandra A. Santos,et al.  On the approximate reanalysis technique in topology optimization , 2018, Optimization and Engineering.

[17]  Ekkehard W. Sachs,et al.  Adaptive Trust-Region POD Methods in PIDE-Constrained Optimization , 2014 .

[18]  Karl Meerbergen,et al.  Accelerating Optimization of Parametric Linear Systems by Model Order Reduction , 2013, SIAM J. Optim..

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

[20]  Geoffrey Oxberry,et al.  Limited‐memory adaptive snapshot selection for proper orthogonal decomposition , 2017 .

[21]  B. Bourdin.  Filters in topology optimization , 2001 .

[22]  A. Patera,et al.  Reduced basis approximation and a posteriori error estimation for affinely parametrized elliptic coercive partial differential equations , 2007 .

[23]  Mathias Stolpe,et al.  An efficient second-order SQP method for structural topology optimization , 2016 .

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

[25]  Charbel Farhat,et al.  Progressive construction of a parametric reduced‐order model for PDE‐constrained optimization , 2014, ArXiv.

[26]  Lorenz T. Biegler,et al.  A trust-region framework for constrained optimization using reduced order modeling , 2013 .

[27]  Matthias Heinkenschloss,et al.  Inexact Objective Function Evaluations in a Trust-Region Algorithm for PDE-Constrained Optimization under Uncertainty , 2014, SIAM J. Sci. Comput..

[28]  P. A. Newman,et al.  Approximation and Model Management in Aerodynamic Optimization with Variable-Fidelity Models , 2001 .

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

[30]  Bart G. van Bloemen Waanders,et al.  A Trust-Region Algorithm with Adaptive Stochastic Collocation for PDE Optimization under Uncertainty , 2012, SIAM J. Sci. Comput..

[31]  Youhong Sun,et al.  An Efficient Reanalysis Method for Topological Optimization of Vibrating Continuum Structures for Simple and Multiple Eigenfrequencies , 2018, Mathematical Problems in Engineering.

[32]  E. Sachs,et al.  Trust-region proper orthogonal decomposition for flow control , 2000 .

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

[34]  L. Sirovich.  Turbulence and the dynamics of coherent structures. I. Coherent structures , 1987 .

[35]  Trenton Kirchdoerfer,et al.  Accelerating design optimization using reduced order models , 2019, ArXiv.

[36]  Luís N. Vicente,et al.  Analysis of Inexact Trust-Region SQP Algorithms , 2002, SIAM J. Optim..

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

[38]  Panos Y. Papalambros,et al.  Structural reanalysis for topological modifications – a unified approach , 2001 .

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

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

[41]  O. Sigmund,et al.  Filters in topology optimization based on Helmholtz‐type differential equations , 2011 .

[42]  Karen Willcox,et al.  A Certified Trust Region Reduced Basis Approach to PDE-Constrained Optimization , 2017, SIAM J. Sci. Comput..

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

[44]  Hung Nguyen-Xuan,et al.  Three-dimensional topology optimization of auxetic metamaterial using isogeometric analysis and model order reduction , 2019, Computer Methods in Applied Mechanics and Engineering.

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

[46]  Stephen J. Wright,et al.  Numerical Optimization (Springer Series in Operations Research and Financial Engineering) , 2000 .

[47]  Matthias Heinkenschloss,et al.  A Matrix-Free Trust-Region SQP Method for Equality Constrained Optimization , 2014, SIAM J. Optim..

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

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