On barrier and modified barrier multigrid methods for 3d topology optimization

One of the challenges encountered in optimization of mechanical structures, in particular in what is known as topology optimization, is the size of the problems, which can easily involve millions of variables. A basic example is the minimum compliance formulation of the variable thickness sheet (VTS) problem, which is equivalent to a convex problem. We propose to solve the VTS problem by the Penalty-Barrier Multiplier (PBM) method, introduced by R.\ Polyak and later studied by Ben-Tal and Zibulevsky and others. The most computationally expensive part of the algorithm is the solution of linear systems arising from the Newton method used to minimize a generalized augmented Lagrangian. We use a special structure of the Hessian of this Lagrangian to reduce the size of the linear system and to convert it to a form suitable for a standard multigrid method. This converted system is solved approximately by a multigrid preconditioned MINRES method. The proposed PBM algorithm is compared with the optimality criteria (OC) method and an interior point (IP) method, both using a similar iterative solver setup. We apply all three methods to different loading scenarios. In our experiments, the PBM method clearly outperforms the other methods in terms of computation time required to achieve a certain degree of accuracy.

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

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

[3]  Michael Zibulevsky,et al.  Penalty/Barrier Multiplier Methods for Convex Programming Problems , 1997, SIAM J. Optim..

[4]  Marc Teboulle,et al.  Hidden convexity in some nonconvex quadratically constrained quadratic programming , 1996, Math. Program..

[5]  Roman A. Polyak,et al.  Modified barrier functions (theory and methods) , 1992, Math. Program..

[6]  Volker Schulz,et al.  Interior point multigrid methods for topology optimization , 2000 .

[7]  Michael Stingl,et al.  On the solution of large-scale SDP problems by the modified barrier method using iterative solvers , 2009, Math. Program..

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

[9]  Stephen J. Wright Primal-Dual Interior-Point Methods , 1997, Other Titles in Applied Mathematics.

[10]  Florian Jarre,et al.  Optimal Truss Design by Interior-Point Methods , 1998, SIAM J. Optim..

[11]  Michal Kočvara,et al.  Mechanical design problems with unilateral contact , 1998 .

[12]  Michal Kočvara Chapter 11: Truss Topology Design by Conic Linear Optimization , 2017 .

[13]  Michael Stingl,et al.  PENNON: A code for convex nonlinear and semidefinite programming , 2003, Optim. Methods Softw..

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

[15]  Michal Kocvara,et al.  Primal-Dual Interior Point Multigrid Method for Topology Optimization , 2016, SIAM J. Sci. Comput..

[16]  R. Polyak Smooth optimization methods for minimax problems , 1988 .

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

[18]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[19]  M. Saunders,et al.  Solution of Sparse Indefinite Systems of Linear Equations , 1975 .

[20]  Martin P. Bendsøe Topology Optimization , 2009, Encyclopedia of Optimization.

[21]  Roland W. Freund,et al.  A QMR-based interior-point algorithm for solving linear programs , 1997, Math. Program..

[22]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .

[23]  Martin P. Bendsøe,et al.  A New Method for Optimal Truss Topology Design , 1993, SIAM J. Optim..