Abstract
We propose an efficient implementation of a new hybrid topology optimization algorithm based on multigrid approach that combines the parallelization strategy of CPU using OpenMP and heavily multithreading capabilities of modern Graphics Processing Units (GPU). In addition to that, significant computational efficiency in memory requirement has been achieved using homogenization strategy. The algorithm has been integrated with versatile computing platform of MATLAB for ease of use and customization. The bottlenecking repetitive solution of the state equation has been solved using an optimized geometric multigrid approach along with CUDA parallelization enabling an order of magnitude faster in computational time than current state of the art implementations. The main novelty lies in the efficient implementation wherein on the fly computation of auxiliary matrices in the multigrid scheme and modification in interpolation schemes using homogenization strategy removes memory limitation of GPUs. Memory hierarchy of GPU has also been exploited for further optimized implementations. All these enable solution of structures involving hundred millions of three dimensional brick elements to be accomplished in a standard desktop computer or a workstation. Performance of the proposed algorithm is illustrated using several examples including design dependent loads. Results obtained indicate the excellent performance and scalability of the proposed approach.
Similar content being viewed by others
Notes
A is computed on the fly and hence only Q is actually stored.
For numerical stability, we consider a very small stiffness of \(1/10^6\).
References
Bendsøe MP, Kikuchi N (1988) Generating optimal topologies in structural design using a homogenization method. Comput Methods Appl Mech Eng 71(2):197–224
Yoon GH, Choi H, Hur S (2018) Multiphysics topology optimization for piezoelectric acoustic focuser. Comput Methods Appl Mech Eng 332:600–623
Zhu X, Zhao C, Wang X, Zhou Y, Hu P, Ma Z-D (2019) Temperature-constrained topology optimization of thermo-mechanical coupled problems. Eng Optim
Kazemi H, Vaziri A, Norato JA (2020) Multi-material topology optimization of lattice structures using geometry projection. Comput Methods Appl Mech Eng 363:112895
Chakraborty S, Goswami S, Rabczuk T (2019) A surrogate assisted adaptive framework for robust topology optimization. Comput Methods Appl Mech Eng 346:63–84
Giraldo-Londoño O, Mirabella L, Dalloro L, Paulino GH (2020) Multi-material thermomechanical topology optimization with applications to additive manufacturing: design of main composite part and its support structure. Comput Methods Appl Mech Eng 363:112812
Liang X, To AC, Du J, Zhang YJ (2021) Topology optimization of phononic-like structures using experimental material interpolation model for additive manufactured lattice infills. Comput Methods Appl Mech Eng 377:113717
Sigmund O, Maute K (2013) Topology optimization approaches. Struct Multidiscip Optim 48(6):1031–1055
Kreissl S, Maute K (2012) Levelset based fluid topology optimization using the extended finite element method. Struct Multidiscip Optim 46(3):311–326
Wang Y, Wang Z, Xia Z, Poh LH (2018) Structural design optimization using isogeometric analysis: a comprehensive review. Comput Model Eng Sci 117(3):455–507
Gai Y, Zhu X, Zhang YJ, Hou W, Hu P (2020) Explicit isogeometric topology optimization based on moving morphable voids with closed B-spline boundary curves. Struct Multidiscip Optim 61(3):963–982
Sigmund O (2020) EML webinar overview: topology optimization-status and perspectives. Extreme Mech Lett 39:100855
Wadbro E, Berggren M (2009) Megapixel topology optimization on a graphics processing unit. SIAM Rev 51(4):707–721
Haase G, Liebmann M, Douglas CC, Plank G (2010) A parallel algebraic multigrid solver on graphics processing units. High performance computing and applications. Springer, Berlin, pp 38–47
Gandham R, Esler K, Zhang Y (2014) A GPU accelerated aggregation algebraic multigrid method. Comput Math Appl 68(10):1151–1160
Schmidt S, Schulz V (2011) A 2589 line topology optimization code written for the graphics card. Comput Vis Sci 14(6):249–256
Martínez-Frutos J, Martínez-Castejón PJ, Herrero-Pérez D (2015) Fine-grained GPU implementation of assembly-free iterative solver for finite element problems. Comput Struct 157:9–18
Xia Z, Wang Y, Wang Q, Mei C (2017) GPU parallel strategy for parameterized LSM-based topology optimization using isogeometric analysis. Struct Multidiscip Optim 56(2):413–434
Dick C, Georgii J, Westermann R (2011) A real-time multigrid finite hexahedra method for elasticity simulation using CUDA. Simul Model Pract Theory 19(2):801–816
Martínez-Frutos J, Martínez-Castejón PJ, Herrero-Pérez D (2017) Efficient topology optimization using GPU computing with multilevel granularity. Adv Eng Softw 106:47–62
Martínez-Frutos J, Herrero-Pérez D (2017) GPU acceleration for evolutionary topology optimization of continuum structures using isosurfaces. Comput Struct 182:119–136
Baiges J, Martínez-Frutos J, Herrero-Pérez D, Otero F, Ferrer A (2019) Large-scale stochastic topology optimization using adaptive mesh refinement and coarsening through a two-level parallelization scheme. Comput Methods Appl Mech Eng 343:186–206
Li H, Li H, Gao L, Li J, Li P, Yang Y (2021) Topology optimization of arbitrary-shape multi-phase structure with structured meshes based on a virtual phase method. Comput Methods Appl Mech Eng 387:114138
Xie X, Wang S, Xu M, Jiang N, Wang Y (2020) A hierarchical spline based isogeometric topology optimization using moving morphable components. Comput Methods Appl Mech Eng 360:112696
Ferrari F, Sigmund O (2020) A new generation 99 line Matlab code for compliance topology optimization and its extension to 3D. Struct Multidiscip Optim 62(4):2211–2228
Bendsøe MP (1989) Optimal shape design as a material distribution problem. Struct Optim 1(4):193–202
Bendsoe MP, Sigmund O (2003) Topology optimization: theory, methods, and applications. Springer, Berlin
Mlejnek H (1992) Some aspects of the genesis of structures. Struct Optim 5(1):64–69
Zhou M, Rozvany G (1991) The COC algorithm, Part II: topological, geometrical and generalized shape optimization. Comput Methods Appl Mech Eng 89(1–3):309–336
Sigmund O, Petersson J (1998) Numerical instabilities in topology optimization: a survey on procedures dealing with checkerboards, mesh-dependencies and local minima. Struct Optim 16(1):68–75
Bourdin B (2001) Filters in topology optimization. Int J Numer Meth Eng 50(9):2143–2158
Zeng S, Lee PS (2019) Topology optimization of liquid-cooled microchannel heat sinks: an experimental and numerical study. Int J Heat Mass Transf 142:118401
Sigmund O (1998) Topology optimization in multiphysics problems. In: 7th AIAA/USAF/NASA/ISSMO symposium on multidisciplinary analysis and optimization, p 4905
Wilson RB A simplicial algorithm for concave programming. Ph. D. Dissertation, Graduate School of Business Administration
Svanberg K (1987) The method of moving asymptotes-a new method for structural optimization. Int J Numer Meth Eng 24(2):359–373
Bendsøe MP, Sigmund O (1995) Optimization of structural topology, shape, and material, vol 414. Springer, Berlin
Papadrakakis M, Stavroulakis G, Karatarakis A (2011) A new era in scientific computing: domain decomposition methods in hybrid CPU-GPU architectures. Comput Methods Appl Mech Eng 200(13–16):1490–1508
Davis TA (2006) Direct methods for sparse linear systems. SIAM, Philadelphia
Hestenes MR, Stiefel E et al (1952) Methods of conjugate gradients for solving linear systems, vol 49. NBS, Washington, DC
Vassilevski PS (2008) Multilevel block factorization preconditioners: matrix-based analysis and algorithms for solving finite element equations. Springer, Berlin
CUDA C++ Programming Guide. NVIDIA Corporation, (2021)
Amir O, Aage N, Lazarov BS (2014) On multigrid-CG for efficient topology optimization. Struct Multidiscip Optim 49(5):815–829
Wu J, Dick C, Westermann R (2015) A system for high-resolution topology optimization. IEEE Trans Vis Comput Gr 22(3):1195–1208
Xie YM, Zuo ZH, Huang X, Black T, Felicetti P (2014) Application of topological optimisation technology to bridge design. Struct Eng Int 24(2):185–191
Cascone F, Faiella D, Tomei V, Mele E (2021) Stress lines inspired structural patterns for tall buildings. Eng Struct 229:111546
Sarkisian M, Lee P, Long E, Shook D (2010) Organic and natural forms in building design. Struct Congress 2010:2840–2851
Standard B Eurocode 1: actions on structures
Standard B (2004) Eurocode 2: design of concrete structures-, Part 1-1: general rules and rules for buildings, p 230
Otterness N, Anderson JH (2020) Amd gpus as an alternative to nvidia for supporting real-time workloads. In: 32nd Euromicro conference on real-time systems (ECRTS 2020), Schloss Dagstuhl-Leibniz-Zentrum für Informatik
Acknowledgements
APP acknowledge the financial support received from IIT Rookee in form of MHRD, Govt. Of India, fellowship. SC acknowledges the financial support received from Science and Engineering Research Board via Grant no. SRG/2021/000467.
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Replication of results
Matlab codes as well as the dependent CUDA C routines for the examples shown in this work are available at https://github.com/csccm-iitd/GPU-TO.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix 1: Three different system configurations used for testing
Appendix 1: Three different system configurations used for testing
See Table 11.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Padhi, A.P., Chakraborty, S., Chakrabarti, A. et al. Efficient hybrid topology optimization using GPU and homogenization-based multigrid approach. Engineering with Computers 39, 3593–3615 (2023). https://doi.org/10.1007/s00366-022-01771-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00366-022-01771-x