Alternating active-phase algorithm for multimaterial topology optimization problems: a 115-line MATLAB implementation

A new algorithm for the solution of multimaterial topology optimization problems is introduced in the present study. The presented method is based on the splitting of a multiphase topology optimization problem into a series of binary phase topology optimization sub-problems which are solved partially, in a sequential manner, using a traditional binary phase topology optimization solver; internal solver. The coupling between these incomplete solutions is ensured using an outer iteration strategy based on the block coordinate descend method. The presented algorithm provides a general framework to extend the traditional binary phase topology optimization solvers for the solution of multiphase topology optimization problems. Interesting features of the presented algorithm are:generality, simplicity and the ease of implementation. The presented algorithm is used to solve multimaterial minimum structural and thermal compliance topology optimization problems based on the classical optimality criteria method. Details of MATLAB implementation are presented and the complete program listings are provided as the Supplementary Materials. The success and performance of the presented method are studied through several two dimensional numerical examples.

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

[2]  Yi Min Xie,et al.  Evolutionary Topology Optimization of Continuum Structures: Methods and Applications , 2010 .

[3]  Hongchao Zhang,et al.  A nonmonotone spectral projected gradient method for large-scale topology optimization problems , 2010 .

[4]  P. Tseng.  Convergence of a Block Coordinate Descent Method for Nondifferentiable Minimization , 2001 .

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

[6]  G. Allaire,et al.  SHAPE OPTIMIZATION OF A SODIUM FAST REACTOR CORE , 2012 .

[7]  Michael Yu Wang,et al.  3D Multi-Material Structural Topology Optimization with the Generalized Cahn-Hilliard Equations , 2006 .

[8]  A. Chambolle,et al.  Design-dependent loads in topology optimization , 2003 .

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

[10]  David G. Luenberger,et al.  Linear and nonlinear programming , 1984 .

[11]  Tony F. Chan,et al.  A Multiphase Level Set Framework for Image Segmentation Using the Mumford and Shah Model , 2002, International Journal of Computer Vision.

[12]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[13]  M. Wang,et al.  Piecewise constant level set method for structural topology optimization , 2009 .

[14]  Xiaoming Wang,et al.  Color level sets: a multi-phase method for structural topology optimization with multiple materials , 2004 .

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

[16]  G. Allaire,et al.  Optimization of nuclear fuel reloading by the homogenization method , 2002 .

[17]  Mei Yulin,et al.  A level set method for structural topology optimization and its applications , 2004 .

[18]  Paul Tseng,et al.  A coordinate gradient descent method for linearly constrained smooth optimization and support vector machines training , 2010, Comput. Optim. Appl..

[19]  O. Sigmund,et al.  Multiphase composites with extremal bulk modulus , 2000 .

[20]  S. Lucidi,et al.  Decomposition Algorithm Model for Singly Linearly-Constrained Problems Subject to Lower and Upper Bounds , 2009 .

[21]  Shiwei Zhou,et al.  Multimaterial structural topology optimization with a generalized Cahn–Hilliard model of multiphase transition , 2006 .

[22]  Pablo Pedregal,et al.  Optimal design of 2D conducting graded materials by minimizing quadratic functionals in the field , 2005 .

[23]  Michael Yu Wang,et al.  Design of piezoelectric actuators using a multiphase level set method of piecewise constants , 2009, J. Comput. Phys..

[24]  S. Torquato,et al.  Composites with extremal thermal expansion coefficients , 1996 .

[25]  E. M. L. Beale,et al.  Nonlinear Programming: A Unified Approach. , 1970 .

[26]  G. Allaire,et al.  Structural optimization using sensitivity analysis and a level-set method , 2004 .

[27]  Martin P. Bendsøe,et al.  Optimization of Structural Topology, Shape, And Material , 1995 .

[28]  Michael Yu Wang,et al.  Synthesis of shape and topology of multi-material structures with a phase-field method , 2004 .

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

[30]  R. Howard,et al.  Local convergence analysis of a grouped variable version of coordinate descent , 1987 .

[31]  Ole Sigmund,et al.  Design of materials with extreme thermal expansion using a three-phase topology optimization method , 1997, Smart Structures.

[32]  T. Chan,et al.  A Variational Level Set Approach to Multiphase Motion , 1996 .

[33]  Y. Xie,et al.  Bi-directional evolutionary topology optimization of continuum structures with one or multiple materials , 2009 .

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

[35]  S. Torquato,et al.  Design of smart composite materials using topology optimization , 1999 .

[36]  Giampaolo Liuzzi,et al.  On the convergence of a Jacobi-type algorithm for singly linearly-constrained problems subject to simple bounds , 2011, Optim. Lett..

[37]  S. Shtrikman,et al.  A variational approach to the theory of the elastic behaviour of multiphase materials , 1963 .

[38]  G. Allaire,et al.  Shape optimization by the homogenization method , 1997 .

[39]  Michael Yu Wang,et al.  A level-set based variational method for design and optimization of heterogeneous objects , 2005, Comput. Aided Des..

[40]  Ole Sigmund.  Recent developments in extremal material design , 2001 .