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

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

摘要/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

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

http://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.

[1]	张萌, 张松林, 刘玉为, 刘时成, 范鹏举. 可调谐半导体激光器压电驱动系统的优化设计[J]. 中国机械工程, 2024, 35(04): 656-665.
[2]	李荣启, 闫涛, 何智成, 米栋, 姜潮, 郑静. 流-热-力耦合的高性能结构拓扑优化设计方法[J]. 中国机械工程, 2024, 35(03): 487-497.
[3]	高进, 崔海冰, 樊涛, 李昂, 杜尊峰. 一种基于自适应Kriging集成模型的结构可靠性分析方法[J]. 中国机械工程, 2024, 35(01): 83-92.
[4]	程贤福, 章志宏, 王承辉, 潘逸飞. 产品架构演化及开放式设计策略[J]. 中国机械工程, 2024, 35(01): 109-124.
[5]	朱宗铭, 季苏强, 王浩, 唐蒲华, 梁亮. 基于流固耦合的介入机器人诊疗时血流动力学分析[J]. 中国机械工程, 2023, 34(21): 2592-2599.
[6]	屈小章, 张加贝, 翟方志. 高速列车散热离心风机性能灵敏性分析及优化[J]. 中国机械工程, 2023, 34(20): 2504-2512.
[7]	肖罡, 张斌, 李时春, 严惠军, 杨钦文. 基于流体压降控制的感应加热甲醇重整装置设计及其服役性能优化[J]. 中国机械工程, 2023, 34(17): 2048-2057,2076.
[8]	屈小章, 张加贝, 翟方志. 非随机载荷不确定的机车侧墙过滤系统两相流可靠性分析[J]. 中国机械工程, 2023, 34(15): 1881-1889.
[9]	耿雪晴, 吴孟武, 华林, . 轮腿式可变形车轮设计及整车控制研究[J]. 中国机械工程, 2023, 34(12): 1446-1452.
[10]	闫萌, 李涛, 杨晨, 王明宇, 杨东东. 面向机电产品绿色设计与评价协同的生命周期评价参数计算方法[J]. 中国机械工程, 2023, 34(12): 1453-1464.
[11]	陶亮, 张大山, 张小龙, 潘登, 占庆良. 智能轮胎开发平台专用轮辋总成设计与试验[J]. 中国机械工程, 2023, 34(09): 1111-1119.
[12]	刘基盛, 计良, 李威, 贾志新, 成金鑫, 方鹏程. 基于多代理模型的离心叶轮高效优化设计方法[J]. 中国机械工程, 2023, 34(08): 899-907.
[13]	王新愿, 周金宇, 谢里阳, 程锦翔. 结构可靠度求解的自适应细分重要抽样法[J]. 中国机械工程, 2023, 34(03): 300-306，313.
[14]	何宇凡, 孙江宏, 高锋, 李乃峥, 何雪萍, 王军见. 一种针对圆锥体外表面贴装的机械手设计分析与优化[J]. 中国机械工程, 2023, 34(01): 55-64.
[15]	丁延冬, 罗年猛, 杨奥迪, 王书亭, 朱浩然, 谢贤达. Bézier单元刚度映射下的高效多重网格等几何拓扑优化方法[J]. 中国机械工程, 2022, 33(23): 2801-2810.