Massively efficient filter for topology optimization based on the splitting of tensor product structure

Abstract

In this work, we put forward a massively efficient filter for topology optimization (TO) utilizing the splitting of tensor product structure. With the aid of splitting technique, the traditional weight matrices of B-splines and non-uniform rational B-spline implicit filters are decomposed equivalently into two or three submatrices, by which the sensitivity analysis is reformulated for the nodal design variables without altering the optimization process. Afterwards, an explicit sensitivity filter, which is decomposed by the splitting pipeline as that applied to implicit filter, is established in terms of the tensor product of the axial distances between adjacent element centroids, and the corresponding sensitivity analysis is derived for elemental design variables. According to the numerical results, the average updating time for the design variables is accelerated by two-order-of-magnitude for the decomposed filter compared with the traditional filter. In addition, the memory burden and computing time of the weight matrix are decreased by six- and three-order-of-magnitude for the decomposed filter. Therefore, the proposed filter is massively improved by the splitting of tensor product structure and delivers a much more efficient way of solving TO problems in the frameworks of isogeometric analysis and finite element analysis.

Graphical abstract

Keywords

topology optimization / isogeometric analysis / finite element analysis / tensor product structure / sensitivity analysis

References

[1] BendsøeM P, SigmundO. Topology Optimization: Theory, Methods, and Applications. 2nd ed. Berlin: Springer, 2003

[2] Bendsøe M P, Kikuchi N. Generating optimal topologies in structural design using a homogenization method. Computer Methods in Applied Mechanics and Engineering, 1988, 71(2): 197–224

[3] Bendsøe M P, Sigmund O. Material interpolation schemes in topology optimization. Archive of Applied Mechanics, 1999, 69(9): 635–654

[4] Sigmund O. A 99 line topology optimization code written in MATLAB. Structural and Multidisciplinary Optimization, 2001, 21(2): 120–127

[5] Xie Y M, Steven G P. A simple evolutionary procedure for structural optimization. Computers & Structures, 1993, 49(5): 885–896

[6] Allaire G, Jouve F, Toader A M. Structural optimization using sensitivity analysis and a level-set method. Journal of Computational Physics, 2004, 194(1): 363–393

[7] Mei Y L, Wang X M. A level set method for structural topology optimization and its applications. Advances in Engineering Software, 2004, 35(7): 415–441

[8] Guo X, Zhang W S, Zhong W L. Doing topology optimization explicitly and geometrically—a new moving morphable components based framework. Journal of Applied Mechanics, 2014, 81(8): 081009

[9] Zhang W H, Zhou Y, Zhu J H. A comprehensive study of feature definitions with solids and voids for topology optimization. Computer Methods in Applied Mechanics and Engineering, 2017, 325: 289–313

[10] Zhang W S, Jiang S, Liu C, Li D D, Kang P, Youn S K, Guo X. Stress-related topology optimization of shell structures using IGA/TSA-based moving morphable void (MMV) approach. Computer Methods in Applied Mechanics and Engineering, 2020, 366: 113036

[11] Zhou Y, Zhang W H, Zhu J H, Xu Z. Feature-driven topology optimization method with signed distance function. Computer Methods in Applied Mechanics and Engineering, 2016, 310: 1–32

[12] Xie X D, Wang S T, Xu M M, Jiang N, Wang Y J. A hierarchical spline based isogeometric topology optimization using moving morphable components. Computer Methods in Applied Mechanics and Engineering, 2020, 360: 112696

[13] Xie X D, Wang S T, Ye M, Xia Z H, Zhao W, Jiang N, Xu M M. Isogeometric topology optimization based on energy penalization for symmetric structure. Frontiers of Mechanical Engineering, 2020, 15(1): 100–122

[14] Xie X D, Yang A D, Wang Y J, Jiang N, Wang S T. Fully adaptive isogeometric topology optimization using MMC based on truncated hierarchical B-splines. Structural and Multidisciplinary Optimization, 2021, 63(6): 2869–2887

[15] Haber R B, Jog C S, Bendsøe M P. A new approach to variable-topology shape design using a constraint on perimeter. Structural Optimization, 1996, 11(1): 1–12

[16] Borrvall T. Topology optimization of elastic continua using restriction. Archives of Computational Methods in Engineering, 2001, 8(4): 351–385

[17] Bourdin B. Filters in topology optimization. International Journal for Numerical Methods in Engineering, 2001, 50(9): 2143–2158

[18] Sigmund O. Morphology-based black and white filters for topology optimization. Structural and Multidisciplinary Optimization, 2007, 33(4): 401–424

[19] Andreassen E, Clausen A, Schevenels M, Lazarov B S, Sigmund O. Efficient topology optimization in MATLAB using 88 lines of code. Structural and Multidisciplinary Optimization, 2011, 43(1): 1–16

[20] BrunsT E, Tortorelli D A. Topology optimization of non-linear elastic structures and compliant mechanisms. Computer Methods in Applied Mechanics and Engineering, 2001, 190(26‒27): 3443–3459

[21] Sigmund O, Maute K. Sensitivity filtering from a continuum mechanics perspective. Structural and Multidisciplinary Optimization, 2012, 46(4): 471–475

[22] Guest J K, Prévost J H, Belytschko T. Achieving minimum length scale in topology optimization using nodal design variables and projection functions. International Journal for Numerical Methods in Engineering, 2004, 61(2): 238–254

[23] Sigmund O. Manufacturing tolerant topology optimization. Acta Mechanica Sinica, 2009, 25(2): 227–239

[24] Xu S L, Cai Y W, Cheng G D. Volume preserving nonlinear density filter based on heaviside functions. Structural and Multidisciplinary Optimization, 2010, 41(4): 495–505

[25] Kawamoto A, Matsumori T, Yamasaki S, Nomura T, Kondoh T, Nishiwaki S. Heaviside projection based topology optimization by a PDE-filtered scalar function. Structural and Multidisciplinary Optimization, 2011, 44(1): 19–24

[26] Lazarov B S, Sigmund O. Filters in topology optimization based on Helmholtz-type differential equations. International Journal for Numerical Methods in Engineering, 2011, 86(6): 765–781

[27] Chen L L, Lu C, Lian H J, Liu Z W, Zhao W C, Li S Z, Chen H B, Bordas S P A. Acoustic topology optimization of sound absorbing materials directly from subdivision surfaces with isogeometric boundary element methods. Computer Methods in Applied Mechanics and Engineering, 2020, 362: 112806

[28] Wallin M, Ivarsson N, Amir O, Tortorelli D. Consistent boundary conditions for PDE filter regularization in topology optimization. Structural and Multidisciplinary Optimization, 2020, 62(3): 1299–1311

[29] Xie X D, Wang S T, Wang Y J, Jiang N, Zhao W, Xu M M. Truncated hierarchical B-spline-based topology optimization. Structural and Multidisciplinary Optimization, 2020, 62(1): 83–105

[30] Xie X D, Yang A D, Jiang N, Wang S T. Topology optimization using fully adaptive truncated hierarchical B-splines. Applied Mathematical Modelling, 2021, 96: 131–151

[31] Xie X D, Yang A D, Jiang N, Zhao W, Liang Z S, Wang S T. Adaptive topology optimization under suitably graded THB-spline refinement and coarsening. International Journal for Numerical Methods in Engineering, 2021, 122(20): 5971–5998

[32] Qian X P. Topology optimization in B-spline space. Computer Methods in Applied Mechanics and Engineering, 2013, 265: 15–35

[33] Wang M M, Qian X P. Efficient filtering in topology optimization via B-splines. In: Proceedings of ASME 2014 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Buffalo: ASME, 2014, V02BT03A011

[34] Costa G, Montemurro M, Pailhès J. A 2D topology optimisation algorithm in NURBS framework with geometric constraints. International Journal of Mechanics and Materials in Design, 2018, 14(4): 669–696

[35] Costa G, Montemurro M, Pailhès J. NURBS hyper-surfaces for 3D topology optimization problems. Mechanics of Advanced Materials and Structures, 2021, 28(7): 665–684

[36] Costa G, Montemurro M, Pailhès J. Minimum length scale control in a NURBS-based SIMP method. Computer Methods in Applied Mechanics and Engineering, 2019, 354: 963–989

[37] Costa G, Montemurro M, Pailhès J, Perry N. Maximum length scale requirement in a topology optimisation method based on NURBS hyper-surfaces. CIRP Annals, 2019, 68(1): 153–156

[38] Costa G, Montemurro M. Eigen-frequencies and harmonic responses in topology optimisation: a CAD-compatible algorithm. Engineering Structures, 2020, 214: 110602

[39] Rodriguez T, Montemurro M, Le Texier P, Pailhès J. Structural displacement requirement in a topology optimization algorithm based on isogeometric entities. Journal of Optimization Theory and Applications, 2020, 184(1): 250–276

[40] Montemurro M, Rodriguez T, Texier P L, Pailhès J. Multi-displacement requirement in a topology optimization algorithm based on non-uniform rational basis spline hyper-surfaces. In: Mariano P M, ed. Variational Views in Mechanics. Cham: Springer, 2021, 223–257

[41] RoinéT, MontemurroM, Pailhès J. Stress-based topology optimization through non-uniform rational basis spline hyper-surfaces. Mechanics of Advanced Materials and Structures, 2021, In press

[42] Montemurro M, Bertolino G, Roiné T. A general multi-scale topology optimisation method for lightweight lattice structures obtained through additive manufacturing technology. Composite Structures, 2021, 258: 113360

[43] Bertolino G, Montemurro M. Two-scale topology optimisation of cellular materials under mixed boundary conditions. International Journal of Mechanical Sciences, 2022, 216: 106961

[44] Montemurro M, Refai K. A topology optimization method based on non-uniform rational basis spline hyper-surfaces for heat conduction problems. Symmetry, 2021, 13(5): 888

[45] Montemurro M, Refai K, Catapano A. Thermal design of graded architected cellular materials through a CAD-compatible topology optimisation method. Composite Structures, 2022, 280: 114862

[46] Montemurro M. On the structural stiffness maximisation of anisotropic continua under inhomogeneous Neumann–Dirichlet boundary conditions. Composite Structures, 2022, 287: 115289

[47] Fernandez F, Puso M A, Solberg J, Tortorelli D A. Topology optimization of multiple deformable bodies in contact with large deformations. Computer Methods in Applied Mechanics and Engineering, 2020, 371: 113288

[48] Wang C, Zhang W H, Zhou L, Gao T, Zhu J H. Topology optimization of self-supporting structures for additive manufacturing with B-spline parameterization. Computer Methods in Applied Mechanics and Engineering, 2021, 374: 113599

[49] Kolda T G, Bader B W. Tensor decompositions and applications. SIAM Review, 2009, 51(3): 455–500

[50] Sidiropoulos N D, De Lathauwer L, Fu X, Huang K J, Papalexakis E E, Faloutsos C. Tensor decomposition for signal processing and machine learning. IEEE Transactions on Signal Processing, 2017, 65(13): 3551–3582

[51] Castellana D, Bacciu D. A tensor framework for learning in structured domains. Neurocomputing, 2022, 470: 405–426

[52] Boussé M, Vervliet N, Domanov I, Debals O, De Lathauwer L. Linear systems with a canonical polyadic decomposition constrained solution: algorithms and applications. Numerical Linear Algebra with Applications, 2018, 25(6): e2190

[53] Scholz F, Mantzaflaris A, Jüttler B. Partial tensor decomposition for decoupling isogeometric Galerkin discretizations. Computer Methods in Applied Mechanics and Engineering, 2018, 336: 485–506

[54] Mantzaflaris A, Jüttler B, Khoromskij B N, Langer U. Low rank tensor methods in Galerkin-based isogeometric analysis. Computer Methods in Applied Mechanics and Engineering, 2017, 316: 1062–1085

[55] de Boor C. On calculating with B-splines. Journal of Approximation Theory, 1972, 6(1): 50–62

[56] Svanberg K. The method of moving asymptotes—a new method for structural optimization. International Journal for Numerical Methods in Engineering, 1987, 24(2): 359–373

[57] Xia Q, Zhou T, Wang M Y, Shi T L. Shape and topology optimization for tailoring the ratio between two flexural eigenfrequencies of atomic force microscopy cantilever probe. Frontiers of Mechanical Engineering, 2014, 9(1): 50–57

[58] Xu M M, Wang S T, Xie X D. Level set-based isogeometric topology optimization for maximizing fundamental eigenfrequency. Frontiers of Mechanical Engineering, 2019, 14(2): 222–234

[59] Long K, Yang X Y, Saeed N, Tian R H, Wen P, Wang X. Topology optimization of transient problem with maximum dynamic response constraint using SOAR scheme. Frontiers of Mechanical Engineering, 2021, 16(3): 593–606

[60] Liu J K, Chen Q, Liang X, To A C. Manufacturing cost constrained topology optimization for additive manufacturing. Frontiers of Mechanical Engineering, 2019, 14(2): 213–221

[61] Wang Y J, Gao L, Qu J P, Xia Z H, Deng X W. Isogeometric analysis based on geometric reconstruction models. Frontiers of Mechanical Engineering, 2021, 16(4): 782–797

Acknowledgements

This work has been supported by the National Key R&D Program of China (Grant No. 2020YFB1708300) and China Postdoctoral Science Foundation (Grant No. 2021M701310). No conflicts of interest exist in this paper. The authors would like to thank the anonymous reviewers, whose suggestions helped to substantially improve this manuscript.

RIGHTS & PERMISSIONS

2022 Higher Education Press