A Hybrid Parallel Strategy for Isogeometric Topology Optimization via CPU/GPU Heterogeneous Computing

This paper aims to solve large-scale and complex isogeometric topology optimization problems that consumesignificant computational resources. A novel isogeometric topology optimization method with a hybrid parallelstrategy of CPU/GPU is proposed, while the hybrid parallel strategies for stiffness matrix assembly, equationsolving, sensitivity analysis, and design variable update are discussed in detail. To ensure the high efficiency ofCPU/GPU computing, a workload balancing strategy is presented for optimally distributing the workload betweenCPU and GPU. To illustrate the advantages of the proposedmethod, three benchmark examples are tested to verifythe hybrid parallel strategy in this paper. The results show that the efficiency of the hybrid method is faster thanserial CPU and parallel GPU, while the speedups can be up to two orders of magnitude.

Full-Scale Isogeometric Topology Optimization of Cellular Structures Based on Kirchhoff-Love Shells

Cellular thin-shell structures are widely applied in ultralightweight designs due to their high bearing capacity and strength-to-weight ratio.In this paper,a full-scale isogeometric topology optimization(ITO)method based on Kirchhoff-Love shells for designing cellular tshin-shell structures with excellent damage tolerance ability is proposed.This method utilizes high-order continuous nonuniform rational B-splines(NURBS)as basis functions for Kirchhoff-Love shell elements.The geometric and analysis models of thin shells are unified by isogeometric analysis(IGA)to avoid geometric approximation error and improve computational accuracy.The topological configurations of thin-shell structures are described by constructing the effective density field on the controlmesh.Local volume constraints are imposed in the proximity of each control point to obtain bone-like cellular structures.To facilitate numerical implementation,the p-norm function is used to aggregate local volume constraints into an equivalent global constraint.Several numerical examples are provided to demonstrate the effectiveness of the proposed method.After simulation and comparative analysis,the results indicate that the cellular thin-shell structures optimized by the proposed method exhibit great load-carrying behavior and high damage robustness.

Explicit Isogeometric Topology Optimization Method with Suitably Graded Truncated Hierarchical B-Spline

This work puts forward an explicit isogeometric topology optimization(ITO)method using moving morphable components(MMC),which takes the suitably graded truncated hierarchical B-Spline based isogeometric analysis as the solver of physical unknown(SGTHB-ITO-MMC).By applying properly basis graded constraints to the hierarchical mesh of truncated hierarchical B-splines(THB),the convergence and robustness of the SGTHB-ITOMMC are simultaneously improved and the tiny holes occurred in optimized structure are eliminated,due to the improved accuracy around the explicit structural boundaries.Moreover,an efficient computational method is developed for the topological description functions(TDF)ofMMC under the admissible hierarchicalmesh,which consists of reducing the dimensionality strategy for design space and the locally computing strategy for hierarchical mesh.We apply the above SGTHB-ITO-MMC with improved efficiency to a series of 2D and 3Dcompliance design problems.The numerical results show that the proposed SGTHB-ITO-MMC method outperforms the traditional THB-ITO-MMCmethod in terms of convergence rate and efficiency.Therefore,the proposed SGTHB-ITO-MMC is an effective way of solving topology optimization(TO)problems.

Multi-Patch Black-White Topology Optimization in Isogeometric Analysis

Topological optimization plays a guiding role in the conceptual design process.This paper conducts research on structural topology optimization algorithm within the framework of isogeometric analysis.For multi-component structures,the Nitsche’smethod is used to glue differentmeshes to performisogeometricmulti-patch analysis.The discrete variable topology optimization algorithm based on integer programming is adopted in order to obtain clear boundaries for topology optimization.The sensitivity filtering method based on the Helmholtz equation is employed for averaging of curved elements’sensitivities.In addition,a simple averaging method along coupling interfaces is proposed in order to ensure the material distribution across coupling areas is reasonably smooth.Finally,the performance of the algorithm is demonstrated by numerical examples,and the effectiveness of the algorithm is verified by comparing it with the results obtained by single-patch and ABAQUS cases.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

A Comprehensive Review of Isogeometric Topology Optimization:Methods,Applications and Prospects

Topology Optimization(TO)is a powerful numerical technique to determine the optimal material layout in a design domain,which has accepted considerable developments in recent years.The classic Finite Element Method(FEM)is applied to compute the unknown structural responses in TO.However,several numerical deficiencies of the FEM significantly influence the effectiveness and efficiency of TO.In order to eliminate the negative influence of the FEM on TO,IsoGeometric Analysis(IGA)has become a promising alternative due to its unique feature that the Computer-Aided Design(CAD)model and Computer-Aided Engineering(CAE)model can be unified into a same mathematical model.In the paper,the main intention is to provide a comprehensive overview for the developments of Isogeometric Topology Optimization(ITO)in methods and applications.Finally,some prospects for the developments of ITO in the future are also presented.

T-Splines Based Isogeometric Topology Optimization with Arbitrarily Shaped Design Domains

In this paper,a new isogeometric topology optimization(ITO)method is proposed by using T-splines based isogeometric analysis(IGA).The arbitrarily shaped design domains,directly obtained from CAD,are represented by a single T-spline surface which overcomes the topological limitations of Non-Uniform Rational B-Spline(NURBS).The coefficients correlated with control points are directly used as design variables.Therefore,the T-spline basis functions applied for geometry description and calculation of structural response are simultaneously introduced to represent the density distribution.Several numerical examples show that the proposed approach leads to a coherent workflow to handle design problems of complicated structures.The optimized results are free of checkerboard patterns without additional stabilization and filtering techniques due to the properties of T-splines,which also simplified the post-processing.In addition,through performing local refinement,we can easily achieve multiresolution optimization and infill optimization within the T-splines based framework.In general,the proposed method provides a possibility to design,analyze,and optimize engineering structures in a uniform model,which has the potential to improve design efficiency and reduce the cost of product development.

Multiscale Isogeometric Topology Optimization with Unified Structural Skeleton

This paper proposes a multiscale isogeometric topology optimization(ITO)method where the configuration and layout of microstructures are optimized simultaneously.At micro scale,a shape deformation method is presented to transform a prototype microstructure(PM)for obtaining a series of graded microstructures(GMs),where microstructural skeleton based on the level set framework is applied to retain more topology features and improve the connectability.For the macro scale calculation,the effective mechanical properties can be estimated by means of the numerical homogenization method.By adopting identical non-uniform rational basis splines(NURBS)as basis functions for both parameterized level set model and isogeometric calculation model,the isogeometric analysis(IGA)is integrated into the level set method,which contributes to improving the accuracy and efficiency.Numerical examples demonstrate that,the proposed method is effective in improving the performance and manufacturability.

A Combined Shape and Topology Optimization Based on Isogeometric Boundary Element Method for 3D Acoustics

A combined shape and topology optimization algorithm based on isogeometric boundary element method for 3D acoustics is developed in this study.The key treatment involves using adjoint variable method in shape sensitivity analysis with respect to non-uniform rational basis splines control points,and in topology sensitivity analysis with respect to the artificial densities of sound absorption material.OpenMP tool in Fortran code is adopted to improve the efficiency of analysis.To consider the features and efficiencies of the two types of optimization methods,this study adopts a combined iteration scheme for the optimization process to investigate the simultaneous change of geometry shape and distribution of material to achieve better noise control.Numerical examples,such as sound barrier,simple tank,and BeTSSi submarine,are performed to validate the advantage of combined optimization in noise reduction,and to demonstrate the potential of the proposed method for engineering problems.

ALTERNATING OPTIMIZATION METHOD FOR ISOGEOMETRIC TOPOLOGY OPTIMIZATION WITH STRESS CONSTRAINTS

Topology optimization (TO) has developed rapidly recently. However, topology optimization with stress constraints still faces many challenges due to its highly non-linear properties which will cause inefficient computation, iterative oscillation, and convergence guarantee problems. At the same time, isogeometric analysis (IGA) is accepted by more and more researchers, and it has become one important tool in the field of topology optimization because of its high fidelity. In this paper, we focus on topology optimization with stress constraints based on isogeometric analysis to improve computation efficiency and stability. A new hybrid solver combining the alternating direction method of multipliers and the method of moving asymptotes (ADMM-MMA) is proposed to solve this problem. We first generate an initial feasible point by alternating direction method of multipliers (ADMM) in virtue of the rapid initial descent property. After that, we adopt the method of moving asymptotes (MMA) to get the final results. Several benchmark examples are used to verify the proposed method, and the results show its feasibility and effectiveness.

Structural Design Optimization Using Isogeometric Analysis: A Comprehensive Review

Isogeometric analysis(IGA),an approach that integrates CAE into conventional CAD design tools,has been used in structural optimization for 10 years,with plenty of excellent research results.This paper provides a comprehensive review on isogeometric shape and topology optimization,with a brief coverage of size optimization.For isogeometric shape optimization,attention is focused on the parametrization methods,mesh updating schemes and shape sensitivity analyses.Some interesting observations,e.g.the popularity of using direct(differential)method for shape sensitivity analysis and the possibility of developing a large scale,seamlessly integrated analysis-design platform,are discussed in the framework of isogeometric shape optimization.For isogeometric topology optimization(ITO),we discuss different types of ITOs,e.g.ITO using SIMP(Solid Isotropic Material with Penalization)method,ITO using level set method,ITO using moving morphable com-ponents(MMC),ITO with phase field model,etc.,their technical details and applications such as the spline filter,multi-resolution approach,multi-material problems and stress con-strained problems.In addition to the review in the last 10 years,the current developmental trend of isogeometric structural optimization is discussed.

An XBi-CFAO Method for the Optimization of Multi-Layered Variable Stiffness Composites Using Isogeometric Analysis

This paper presents an effective fiber angle optimization method for two and multi-layered variable stiffness composites.A gradient-based fiber angle optimization method is developed based on isogeometric analysis(IGA).Firstly,the element densities and fiber angles for two and multi-layered composites are synchronously optimized using an extended Bi-layered continuous fiber angle optimization method(XBi-CFAO).The densities and fiber angles in the base layer are attached to the control points.The structure response and sensitivity analysis are accomplished using the non-uniform rational B-spline(NURBS)based IGA.By the benefit of the B-spline space,this method is free from checkerboards,and no additional filtering is needed to smooth the sensitivity numbers.Then the curved fiber paths are generated using the streamline method and the discontinuous fiber paths are smoothed using a partitioned selection process.The proposed method in the paper can alleviate the phenomenon of fiber discontinuity,enhance information retention for the optimized fiber angles of the singular points and save calculating resources effectively.

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

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

Noise Pollution Reduction through a Novel Optimization Procedure in Passive Control Methods

This paper proposes a novel optimization framework in passive control techniques to reduce noise pollution.The geometries of the structures are represented by Catmull-Clark subdivision surfaces,which are able to build gap-free Computer-Aided Design models and meanwhile tackle the extraordinary points that are commonly encountered in geometricmodelling.The acoustic fields are simulated using the isogeometric boundary elementmethod,and a density-based topology optimization is conducted to optimize distribution of sound-absorbing materials adhered to structural surfaces.The approach enables one to perform acoustic optimization from Computer-Aided Design models directly without needingmeshing and volume parameterization,thereby avoiding the geometric errors and time-consuming preprocessing steps in conventional simulation and optimization methods.The effectiveness of the present method is demonstrated by three dimensional numerical examples.

基于等几何分析的板壳结构形状拓扑协同优化

板壳结构轻量化设计是工程中的常见问题,采用有限元法对板壳结构进行拓扑优化时,难以获得高质量的网格以精准描述几何模型,且单独使用拓扑优化受制于初始设计空间。提出了一种基于等几何分析的形状拓扑协同优化方法,该方法执行先形状优化、后拓扑优化的步骤,在最佳结构形状的基础上实现最佳材料布局的优化目标。相比于其他组合形式的优化机制,极大地提高了计算精度与计算效率。采用3个数值算例验证了本方法的有效性与高效性,结果表明,与经典等几何拓扑优化方法相比,可得到更高性能的优化结构,有助于进一步拓宽结构优化的应用范围。

基于自适应泡泡法的薄壳结构拓扑优化设计

为有效解决薄壳结构拓扑优化设计难题,并满足其对分析模型精度和优化结果质量的高要求,结合等几何壳体分析方法提出一种基于自适应泡泡法的新型拓扑优化设计框架.等几何分析技术在薄壳分析方面具有天然的优势:一方面可为薄壳结构建立起精确的NURBS分析模型,避免了模型转换操作及误差;另一方面还可保证待分析物理场的高阶连续性,无需设置转角自由度等.为了在给定壳面上实现结构的拓扑演化,借助NURBS曲面(即等几何分析中的薄壳中面)的映射关系,仅需在规则的二维参数区域内改变结构拓扑即可.鉴于此,采用自适应泡泡法在壳面参数区域内开展拓扑优化,该方法包含孔洞建模、孔洞引入和固定网格分析3个模块,其在当前工作中分别基于闭合B样条、拓扑导数理论和有限胞元法实现.其中,闭合B样条兼具参数和隐式两种表达形式,参数形式便于在CAD系统中直接生成精确的结构模型;隐式形式不仅便于开展孔洞的融合/分离操作,还能与有限胞元法有机结合以替代繁琐的修剪曲面分析方法.理论分析和数值算例表明,所提优化设计框架将复杂的薄壳结构拓扑优化问题转化为简单的二维结构拓扑优化问题,在保证足够分析精度的基础上使用相对很少的设计变量就可得到具有清晰光滑边界且便于导入到CAD系统的优化结果.

基于等几何方法的功能梯度Mindlin板结构拓扑优化设计

提出一种基于等几何分析的功能梯度Mindlin板结构拓扑优化方法,实现了沿着板面内方向梯度材料渐变与结构拓扑协同设计。基于等几何分析中的非均匀有理B样条(non-uniform rational B-spline,NURBS)基函数,保证了板结构建模、仿真分析和优化设计的参数化统一。由于等几何方法基函数的高阶连续性、非负性、规范性等特点,建立了合理的设计变量与梯度材料相对密度的分布关系。为了克服Mindlin板仿真分析的数值闭锁问题,采用NURBS等几何分析单元的缩减积分策略。数值算例展示了该方法的有效性,且获得了具有清晰细节特征的功能梯度板新型设计方案。结果表明,相比于预先定义材料配比,该方法可以实现功能梯度材料属性与结构拓扑的协同优化设计。

从计算机辅助设计(CAD)到人辅助设计(HAD)--一种等几何拓扑优化方法

本文提出了一种新的设计模式——人辅助设计,以取代传统的计算机辅助设计。在人辅助设计中,计算机可以通过一种新的等几何拓扑优化自动完成整个产品设计,而人类仅需协助轻微修改设计以满足要求。文中提出了一种嵌入域等几何拓扑优化用于设计具有不规则设计域的复杂模型,并且可以基于分层等几何拓扑优化结果自动生成优化结果的可编辑几何模型。测试了三个算例以验证所提出的等几何拓扑优化模式,包括一个具有规则设计域的3D悬臂梁,一个具有不规则设计域的汽车零件和一个具有多尺度结构的MBB梁。结果表明,所提出的等几何拓扑优化模式可以自动生成高质量的优化模型。因此,该技术具有成为革命性技术的巨大潜力,将当前设计模式由计算机辅助设计转变为人辅助设计。

Isogeometric topology optimization based on energy penalization for symmetric structure

We present an energy penalization method for isogeometric topology optimization using moving morphable components(ITO–MMC),propose an ITO–MMC with an additional bilateral or periodic symmetric constraint for symmetric structures,and then extend the proposed energy penalization method to an ITO–MMC with a symmetric constraint.The energy penalization method can solve the problems of numerical instability and convergence for the ITO–MMC and the ITO–MMC subjected to the structural symmetric constraint with asymmetric loads.Topology optimization problems of asymmetric,bilateral symmetric,and periodic symmetric structures are discussed to validate the effectiveness of the proposed energy penalization approach.Compared with the conventional ITO–MMC,the energy penalization method for the ITO–MMC can improve the convergence rate from 18.6%to 44.5%for the optimization of the asymmetric structure.For the ITO–MMC under a bilateral symmetric constraint,the proposed method can reduce the objective value by 5.6%and obtain a final optimized topology that has a clear boundary with decreased iterations.For the ITO–MMC under a periodic symmetric constraint,the proposed energy penalization method can dramatically reduce the number of iterations and obtain a speedup of more than 2.

Level set-based isogeometric topology optimization for maximizing fundamental eigenfrequency

Maximizing the fundamental eigenfrequency is an efficient means for vibrating structures to avoid resonance and noises.In this study,we develop an isogeometric analysis(IGA)-based level set model for the fonnulation and solution of topology optimization in cases with maximum eigenfrequency.The proposed method is based on a combination of level set method and IGA technique,which uses the non-uniform rational B-spline(NURBS),description of geometry,to perfonn analysis.The same NURBS is used for geometry representation,but also for IGA-based dynamic analysis and parameterization of the level set surface,that is,the level set function.The method is applied to topology optimization problems of maximizing the fundamental eigenfrequency for a given amount of material.A modal track method,that monitors a single target eigenmode is employed to prevent the exchange of eigenmode order number in eigenfrequency optimization.The validity and efficiency of the proposed method are illustrated by benchmark examples.