Due to the high-order B-spline basis functions utilized in isogeometric analysis(IGA)and the repeatedly updating global stiffness matrix of topology optimization,Isogeometric topology optimization(ITO)intrinsically su...Due to the high-order B-spline basis functions utilized in isogeometric analysis(IGA)and the repeatedly updating global stiffness matrix of topology optimization,Isogeometric topology optimization(ITO)intrinsically suffers from the computationally demanding process.In this work,we address the efficiency problem existing in the assembling stiffness matrix and sensitivity analysis using B˙ezier element stiffness mapping.The Element-wise and Interaction-wise parallel computing frameworks for updating the global stiffness matrix are proposed for ITO with B˙ezier element stiffness mapping,which differs from these ones with the traditional Gaussian integrals utilized.Since the explicit stiffness computation formula derived from B˙ezier element stiffness mapping possesses a typical parallel structure,the presented GPU-enabled ITO method can greatly accelerate the computation speed while maintaining its high memory efficiency unaltered.Numerical examples demonstrate threefold speedup:1)the assembling stiffness matrix is accelerated by 10×maximumly with the proposed GPU strategy;2)the solution efficiency of a sparse linear system is enhanced by up to 30×with Eigen replaced by AMGCL;3)the efficiency of sensitivity analysis is promoted by 100×with GPU applied.Therefore,the proposed method is a promising way to enhance the numerical efficiency of ITO for both single-patch and multiple-patch design problems.展开更多
基金supported by the National Key R&D Program of China(2023YFB2504601)National Natural Science Foundation of China(52205267).
文摘Due to the high-order B-spline basis functions utilized in isogeometric analysis(IGA)and the repeatedly updating global stiffness matrix of topology optimization,Isogeometric topology optimization(ITO)intrinsically suffers from the computationally demanding process.In this work,we address the efficiency problem existing in the assembling stiffness matrix and sensitivity analysis using B˙ezier element stiffness mapping.The Element-wise and Interaction-wise parallel computing frameworks for updating the global stiffness matrix are proposed for ITO with B˙ezier element stiffness mapping,which differs from these ones with the traditional Gaussian integrals utilized.Since the explicit stiffness computation formula derived from B˙ezier element stiffness mapping possesses a typical parallel structure,the presented GPU-enabled ITO method can greatly accelerate the computation speed while maintaining its high memory efficiency unaltered.Numerical examples demonstrate threefold speedup:1)the assembling stiffness matrix is accelerated by 10×maximumly with the proposed GPU strategy;2)the solution efficiency of a sparse linear system is enhanced by up to 30×with Eigen replaced by AMGCL;3)the efficiency of sensitivity analysis is promoted by 100×with GPU applied.Therefore,the proposed method is a promising way to enhance the numerical efficiency of ITO for both single-patch and multiple-patch design problems.
