基于GPU加速的等几何拓扑优化高效多重网格求解方法

doi:10.3969/j.issn.1004-132X.2024.04.004

中国机械工程 ›› 2024, Vol. 35 ›› Issue (04): 602-613.DOI: 10.3969/j.issn.1004-132X.2024.04.004

基于GPU加速的等几何拓扑优化高效多重网格求解方法

杨峰1;罗世杰1;杨江鸿1;王英俊1,2

1.华南理工大学机械与汽车工程学院聚合物新型成型装备国家工程研究中心，广州，510641
2.华中科技大学数字制造装备与技术国家重点实验室，武汉，430074

出版日期:2024-04-25 发布日期:2024-05-24
通讯作者: 王英俊（通信作者），男，1984年生，教授、博士研究生导师。研究方向为拓扑优化、等几何分析、CAD/CAE一体化。E-mail：wangyj84@scut.edu.cn。
作者简介:杨峰，男，1998年生，硕士研究生。研究方向为拓扑优化、GPU并行计算。
基金资助:
国家自然科学基金（52075184）；数字制造装备与技术国家重点实验室开放基金（DMETKF2021020）

A GPU-accelerated High-efficient Multi-grid Algorithm for ITO

YANG Feng1;LUO Shijie1;YANG Jianghong1;WANG Yingjun1,2

1.National Engineering Research Center of Novel Equipment for Polymer Processing,School of
Mechanical and Automotive Engineering,South China University of Technology,Guangzhou,510641
2.State Key Laboratory of Digital Manufacturing Equipment and Technology,Huazhong University
of Science and Technology,Wuhan,430074

Online:2024-04-25 Published:2024-05-24

摘要/Abstract

摘要： 针对大规模等几何拓扑优化(ITO)计算量巨大、传统求解方法效率低的问题，提出了一种基于样条h细化的高效多重网格方程求解方法。该方法利用h细化插值得到粗细网格之间的权重信息，然后构造多重网格方法的插值矩阵，获得更准确的粗细网格映射信息，从而提高求解速度。此外，对多重网格求解过程进行分析，构建其高效GPU并行算法。数值算例表明，所提出的求解方法与线性插值的多重网格共轭梯度法、代数多重网格共轭梯度法和预处理共轭梯度法相比分别取得了最高1.47、11.12和17.02的加速比。GPU并行求解相对于CPU串行求解的加速比高达33.86，显著提高了大规模线性方程组的求解效率。

关键词: 等几何拓扑优化, 方程组求解, h细化, 多重网格法, GPU并行计算

Abstract: An efficient multi-grid equation solving method was proposed based on the h-refinement of splines to address the challenges posed by large-scale ITO computation and low efficiency of traditional solving methods. By the proposed method, the weight information obtained through h-refinement interpolation between coarse and fine grids was used to construct the interpolation matrix of the multi-grid method, thereby enhancing the accuracy of mapping information for both coarse and fine grids and improving computational efficiency. Additionally, a comprehensive analysis of the multi-grid solving process was conducted, culminating in the development of an efficient GPU parallel algorithm. Numerical examples illustrate that the proposed method outperforms existing methods, demonstrating speedup ratios of 1.47, 11.12, and 17.02 in comparison to the linear interpolation multi-grid conjugate gradient method algebraic multi-grid conjugate gradient method, and pre-processing conjugate gradient method respectively. Furthermore, the acceleration rate of GPU parallel solution surpasses that of CPU serial solution by 33.86 times, which significantly enhances the efficiency of solving large-scale linear equations.

Key words: isogeometric topology optimization(ITO), system of equations, h-refinement, multi-grid method, GPU parallel computing

中图分类号:

TH122

杨峰, 罗世杰, 杨江鸿, 王英俊, . 基于GPU加速的等几何拓扑优化高效多重网格求解方法[J]. 中国机械工程, 2024, 35(04): 602-613.

YANG Feng, LUO Shijie, YANG Jianghong, WANG Yingjun, . A GPU-accelerated High-efficient Multi-grid Algorithm for ITO[J]. China Mechanical Engineering, 2024, 35(04): 602-613.

导出引用管理器 EndNote|Ris|BibTeX

链接本文: https://www.cmemo.org.cn/CN/10.3969/j.issn.1004-132X.2024.04.004

https://www.cmemo.org.cn/CN/Y2024/V35/I04/602

参考文献

［1］杨雨豪,郑伟,王英俊.一种自由度缩减和收敛加速的高效等几何拓扑优化方法［J］.中国机械工程,2022,33(23):2811-2821.
YANG Yuhao, ZHENG Wei, WANG Yingjun. An Efficient Isogeometric Topology Optimization Methods Using DOF Reduction and Convergence Acceleration［J］. China Mechanical Engineering, 2022, 33(23):2811-2821.
［2］BENDSE M P. Optimal Shape Design as a Material Distribution Problem［J］. Structural Optimization, 1989, 1(4):193-202.
［3］BENDSE M P, KIKUCHI N. Generating Optimal Topologies in Structural Design Using a Homogenization Method［J］. Computer Methods in Applied Mechanics and Engineering, 1988, 71(2):197-224.
［4］ALLAIRE G, JOUVE F, TOADER A-M. ALevel-set Method for Shape Optimization［J］. Comptes Rendus Mathematique, 2002, 334(12):1125-1130.
［5］QUERIN O M, STEVEN G P, XIE Y M. Evolutionary Structural Optimisation(ESO) Using a Bidirectional Algorithm［J］. Engineering Computations, 1998, 15(8):1031-1048.
［6］GUO X, ZHANG W, ZHANG J, et al. ExplicitStructural Topology Optimization Based on Moving Morphable Components(MMC) with Curved Skeletons［J］. Computer Methods in Applied Mechanics and Engineering, 2016, 310:711-748.
［7］WANG Chuang, ZHU Jihong, WU Manqiao, et al. Multi-scale Design and Optimization for Solid-lattice Hybrid Structures and Their Application to Aerospace Vehicle Components［J］. Chinese Journal of Aeronautics, 2021, 34(5):386-398.
［8］GAO B, MENG D, SHI W, et al. Topology Optimization and the Evolution Trends of Two-speed Transmission of EVs［J］. Renewable and Sustainable Energy Reviews, 2022, 161:112390.
［9］BRANDO P, SILVA P T, PARENTE M, et al. μSmartScope—Towards a Low-cost Microscopic Medical Device for Cervical Cancer Screening Using Additive Manufacturing and Optimization［J］. Proceedings of the Institution of Mechanical Engineers, Part L:Journal of Materials:Design and Applications, 2022, 236(2):267-279.
［10］LI Y, LAI Y, LU G, et al. Innovative Design of Long-span Steel-concrete Composite Bridge Using Multi-material Topology Optimization［J］. Engineering Structures, 2022, 269:114838.
［11］丁延冬,罗年猛,杨奥迪，等.Bézier单元刚度映射下的高效多重网格等几何拓扑优化方法［J］.中国机械工程,2022,33(23):2801-2810.
DING Yandong, LUO Nianmeng, YANG Aodi, et al. Efficient Multigrid Isogeometric Topology Optimization under Bézier Element Stiffness Mapping［J］. China Mechanical Engineering, 2022, 33(23):2801-2810.
［12］HUGHES T J R, COTTRELL J A, BAZILEVS Y. Isogeometric Analysis:CAD, Finite Elements, NURBS, Exact Geometry and Mesh Refinement［J］. Computer Methods in Applied Mechanics and Engineering, 2005, 194(39/41):4135-4195.
［13］SEO Y D, KIM H J, YOUN S K. IsogeometricTopology Optimization Using Trimmed Spline Surfaces［J］. Computer Methods in Applied Mechanics and Engineering, 2010, 199(49/52):3270-3296.
［14］KUMAR A V, PARTHASARATHY A. Topology Optimization Using B-spline Finite Elements［J］. Structural and Multidisciplinary Optimization, 2011, 44(4):471-481.
［15］GAO J, WANG L, LUO Z, et al. IgaTop:an Implementation of Topology Optimization for Structures Using IGA in MATLAB［J］. Structural and Multidisciplinary Optimization, 2021, 64(3):1669-1700.
［16］WANG Y, WANG Z, XIA Z, et al. Structural Design Optimization Using Isogeometric Analysis:a Comprehensive Review［J］. Computer Modeling in Engineering & Sciences, 2018, 117(3):455-507.
［17］GAO J, XIAO M, ZHANG Y, et al. A Comprehensive Review of Isogeometric Topology Optimization:Methods, Applications and Prospects［J］. Chinese Journal of Mechanical Engineering, 2020, 33(6):34-47.
［18］NAGY A P, ABDALLA M M, GURDAL Z. Isogeometric Sizing and Shape Optimization of Beam Structures［J］. Computer Methods in Applied Mechanics and Engineering, 2010,199:1216-1230.
［19］NGUYEN D M，ANTON E，ALLAN R G，et al. Isogeometric Shape Optimization of Vibrating Membranes［J］. Computer Methods in Applied Mechanics and Engineering, 2011, 200:1343-1353.
［20］YOUN D H, KIM M G, KIM H S. Shape Design Optimization of SPH Fluid-structure Interactions Considering Geometrically Exact Interfaces［J］. Structural and Multidisciplinary Optimization, 2011, 44:319-336.
［21］刘宏亮,祝雪峰,杨迪雄.基于等几何分析的结构优化设计研究进展［J］.固体力学学报,2018,39(3):248-267.
LIU Hongliang, ZHU Xuefeng, YANG Dixiong. Research Advances in Isogeometric Analysis-based Optimum Design of Structure［J］. Chinese Journal of Solid Mechanics, 2018, 39(3):248-267.
［22］AAGE N, ANDREASSEN E, LAZAROV B S, et al.Giga-voxel Computational Morphogenesis for Structural Design［J］. Nature, 2017, 550(7674):84-86.
［23］THOMAS S J, ANANTHAN S, YELLAPANTULA S, et al. AComparison of Classical and Aggregation-based Algebraic Multigrid Preconditioners for High-fidelity Simulation of Wind Turbine Incompressible Flows［J］. SIAM Journal on Scientific Computing, 2019, 41(5):S196-S219.
［24］刘石,陈德祥,冯永新，等.等几何分析的多重网格共轭梯度法［J］.应用数学和力学,2014,35(6):630-639.
LIU Shi, CHEN Dexiang, FENG Yongxin, et al. A Multigrid Preconditioned Conjugate Gradient Method for Isogeometric Analysis［J］. Applied Mathematics and Mechanics, 2014, 35(6):630-639.
［25］WANG Y, LIAO Z, YE M, et al. An Efficient Isogeometric Topology Optimization Using Multilevel Mesh, MGCG and Local-update Strategy［J］. Advances in Engineering Software, 2020, 139:102733.
［26］邓亮,徐传福,刘巍，等.交替方向隐式CFD解法器的GPU并行计算及其优化［J］.计算机应用,2013,33(10):2783-2786.
DENG Liang, XU Chuanfu, LIU Wei, et al. Parallelization and Optimization of Alternating Direction Implicit CFD Solver on GPU［J］. Journal of Computer Applications, 2013, 33(10):2783-2786.
［27］XIA Z, WANG Y, WANG Q, et al. GPU Parallel Strategy for Parameterized LSM-based Topology Optimization Using Isogeometric Analysis［J］. Structural and Multidisciplinary Optimization, 2017, 56:413-434.
［28］韩琪,蔡勇.基于GPU的大规模拓扑优化问题并行计算方法［J］.计算机仿真,2015,32(4):221-226.
HAN Qi, CAI Yong. An Efficient Parallel Implementation of Large-scale Topology Optimization Problems Using GPU［J］. Computer Simulation, 2015, 32(4):221-226.
［29］AVNAT O, YAVNEH I. On the Recursive Structure of Multigrid Cycles［J］. SIAM Journal on Scientific Computing, 2022(V3):S103-S126.
［30］MADAMS A, SIFAKIS E, TERAN J. A Parallel Multigrid Poisson Solver for Fluids Simulation on Large Grids［C］∥Proceedings of the 2010 ACM SIGGRAPH/Eurographics Symposium on Computer Animation. Madrid,2010:65-74.
［31］LIU Y, SCHMIDT B. Lightspmv:Faster Cuda-compatible Sparse Matrix-vector Multiplication Using Compressed Sparse Rows［J］. Journal of Signal Processing Systems, 2018, 90:69-86.
［32］LIU Y, GONZLEZ-DOMNGUEZ J, SCHMIDT B. Faster Compressed Sparse Row(CSR)-based Sparse Matrix-vector Multiplication Using CUDA［C］∥GPU Technology Conference. Silicon Valley, 2015:17-20.
［33］BECK R. Algebraic Multigrid by Component Splitting for Edge Elements on Simplicial Triangulations［R］. ZIB-Report(SC-99-40),1999.
［34］NOTAY Y. Aggregation-based Algebraic Multigrid for Convection-diffusion Equations［J］. SIAM Journal on Scientific Computing, 2012, 34(4):A2288-A2316.
［35］BERNASCHI M, DAMBRA P, PASQUINI D. AMG Based on Compatible Weighted Matching for GPUs［J］. Parallel Computing, 2020, 92:102599.
［36］WANG Y, GAO L, QU J, et al. Isogeometric Analysis Based on Geometric Reconstruction Models［J］. Frontiers of Mechanical Engineering, 2021,16:782-797.

基于GPU加速的等几何拓扑优化高效多重网格求解方法

A GPU-accelerated High-efficient Multi-grid Algorithm for ITO

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 1

编辑推荐

Metrics