基于固定网格和拓扑导数的结构拓扑优化自适应泡泡法

为继承传统拓扑优化泡泡法变量少、精度高等优点,并克服其网格重划频繁、孔洞合并操作繁琐等不足,提出了一种基于固定网格和拓扑导数的自适应泡泡方法.该方法的主要特点是:(1)采用有限胞元固定网格分析方法计算结构力学响应,在优化过程中无需网格更新和重划分,就能保证较高的分析精度;(2)根据拓扑导数信息指导结构区域中孔洞的引入,不仅消除了优化结果对孔洞初始布局的依赖性,还能有效控制设计变量的数量;(3)引入拓扑导数阈值和孔洞影响区域新概念,实现了孔洞引入频次和位置的自适应调节,保证了拓扑优化过程的数值计算稳定性;(4)采用光滑变形隐式曲线描述孔洞边界,不仅设计参数少、变形能力强,而且便于处理孔洞间的融合/分离操作以及与固定网格分析方法的有机结合.理论分析和数值算例表明,改进后的自适应泡泡法能够消除传统泡泡法因采用拉格朗日网格和参数化B样条曲线模型而存在的实施困难,采用很少的设计变量就可获得边界光滑清晰的优化结果.     

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

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

基于移动可变形孔洞方法的超弹性结构拓扑优化

现有隐式拓扑优化方法在进行超弹性结构拓扑优化设计时，具有设计变量多、中间设计有限元分析存在严重的收敛性和设计结果无法直接导入CAD/CAE系统等问题。为解决这些问题，提出了一种基于移动可变形孔洞的显式拓扑优化方法来进行承受大变形的超弹性结构设计，材料本构采用常用的Mooney-Rivlin模型。首先，介绍了移动可变形孔洞方法的基本思想和可变形孔洞的显式描述方法；其次，构造了基于移动可变形孔洞方法的超弹性结构拓扑优化的数学列式，给出了相应的灵敏度结果；最后，通过数值算例验证了本方法的有效性。数值结果表明，该方法可以通过较少的设计变量和非常稳健的优化过程，给出边界由B样条曲线描述且可与CAD/CAE软件无缝连接的超弹性结构设计。     

积分背景网格对自适应EFG法拓扑优化的影响研究

采用应力能量范数作为误差指标,探讨了EFG法中积分背景网格对计算精度的影响,得到了合理划分背景网格的建议;建立了以节点密度为设计变量、以最小化柔度为优化日标的拓扑优化模型。采用以节点密度值为加点判据的自适应规则加点方案,开展了连续体结构的拓扑优化研究,该加点方案能有效地减少设计变量的个数,探讨了背景网格对拓扑优化结果的影响。算例结果表明,采用合适的背景网格不仅能进一步减少设计变量的个数,而且能够改善拓扑优化结果的光滑性,使计算效率和精度得到提高。     

基于拓扑描述函数的连续体结构拓扑优化方法

提出了一种利用拓扑描述函数(TDF)作为拓扑设计变量求解连续体结构拓扑优化问题 的新方法. 优化问题的目标函数是结构的整体柔顺性，约束条件为对于可利用材料的体积限 制. 这种方法不仅可以消除拓扑优化中经常出现的棋盘格式等数值不稳定现象，而且能够有 效地抑制传统算法处理此类优化问题时所引发的边界扩散效应. 与其它的基于水平集描述函 数的拓扑优化方法相比，所提出的算法不仅无需求解控制水平集函数演化的双曲守恒方 程，而且合理地考虑了目标函数的拓扑导数信息，因而使得算法的计算效率有了显著的提高.     

基于移动可变形孔洞方法的超弹性结构拓扑优化

现有隐式拓扑优化方法在进行超弹性结构拓扑优化设计时,具有设计变量多、中间设计有限元分析存在严重的收敛性和设计结果无法直接导入CAD/CAE系统等问题。为解决这些问题,提出了一种基于移动可变形孔洞的显式拓扑优化方法来进行承受大变形的超弹性结构设计,材料本构采用常用的Mooney-Rivlin模型。首先,介绍了移动可变形孔洞方法的基本思想和可变形孔洞的显式描述方法;其次,构造了基于移动可变形孔洞方法的超弹性结构拓扑优化的数学列式,给出了相应的灵敏度结果;最后,通过数值算例验证了本方法的有效性。数值结果表明,该方法可以通过较少的设计变量和非常稳健的优化过程,给出边界由B样条曲线描述且可与CAD/CAE软件无缝连接的超弹性结构设计。     

基于拓扑描述函数的特定性能复合材料设计

提出了一种基于拓扑描述函数进行特定性能复合材料设计的新方法. 采用拓扑描 述函数作为设计变量，把复合材料微结构的设计问题转化为一个在周期性单胞上的拓扑优 化问题. 拓扑描述函数以及相应正则化机制的引入不仅可以有效消除棋盘格式等数值不稳定 现象，而且能够有效地抑制传统算法处理此类优化问题时所引发的边界扩散效应. 数值结果 表明所提出的方法可以稳定高效地实现具有特定性能的复合材料微结构设计.     

考虑最小尺寸精确控制的SIMP和MMC混合拓扑优化方法

拓扑优化作为一种先进设计方法, 已被成功用于多个学科领域优化问题求解, 但从拓扑优化结果到其工程应用之间仍存在诸多阻碍, 如在结构设计中存在难以制造的小孔或边界裂缝和单铰链连接等. 在拓扑优化设计阶段考虑结构最小尺寸控制是解决上述问题的一种有效手段. 在最小尺寸控制的结构拓扑优化方法中, 通用性较强的固体各向同性材料惩罚法SIMP优化结果边界模糊不光滑, 包含精确几何信息的移动变形组件法MMC对初始布局具有较强依赖性. 本文提出一种考虑最小尺寸精确控制的SIMP和MMC混合拓扑优化方法. 所提方法继承了二者优势, 避免了各自缺点. 在该方法中, 首先采用活跃轮廓算法ACWE获取SIMP输出的拓扑结构边界轮廓数据, 提出了SIMP优化结果到MMC组件初始布局的映射方法. 其次, 通过引入组件的3个长度变量, 建立了半圆形末端的多变形组件拓扑描述函数模型. 最后, 以组件厚度变量为约束, 构建了考虑结构最小尺寸控制的拓扑优化模型. 采用最小柔度问题和柔性机构问题验证了所提方法的有效性. 数值结果表明, 所提方法在无需额外约束的条件下, 仅通过组件厚度变量下限设置, 可实现整体结构的最小尺寸精确控制, 并获得了具有全局光滑的拓扑结构边界.      

固定网格下基于光滑变形隐式曲线的孔形优化设计方法

孔形优化设计是减缓开孔结构孔边应力集中的有效手段，其对应力分析精度和孔边曲线的几何表达能力都有着较高的要求．论文针对现有固定网格孔形优化方法分析精度较低和/或孔边曲线变形能力有限等不足，采用有限胞元方法在固定网格下对开孔结构进行高精度的力学性能分析，采用能够自由光滑变形的隐式曲线描述待优化的孔洞边界，进而建立优化模型并推导出相应的解析灵敏度计算公式，形成了固定网格下开孔结构孔形优化设计新框架．通过对不同载荷边界条件下的开孔平板结构进行孔形优化设计，表明本方法无需网格更新、灵敏度推导简便、力学分析精度高且优化设计空间大，能够有效降低孔边应力集中程度．     

基于节点独立变量的连续体结构动态拓扑优化

提出基于节点独立变量的连续体结构动态拓扑优化方法.以动态结构响应量最小或最大为目标,体积比为约束,建立了动态结构拓扑优化模型.基于数字图像处理过滤技术得到清晰、边界光滑和体现网格无关性的优化结果.通过二维结构数值算例对理论方法进行验证.结果表明,该方法在连续体结构动态拓扑优化设计中具有可行性和有效性.     

基于网格退化和重组技术的类桁架拓扑优化方法

发展了一种基于有限元网格退化和重组技术的类桁架拓扑优化方法,可在不改变设计域包络的情况下(如维持流型、艺术外观和附加功能等)对结构进行减重。首先,分别对二维问题和三维问题发展相应的有限元网格退化算法,并在刚度等效的意义上对网格进行重组,建立了具有杆元拓扑特征的有限元模型。其次,以全局种子网格的长度尺寸和杆元横截面积为优化变量,构造了域内双层驱动拓扑优化问题,得到具有最优体分比的杆元拓扑结构。数值算例表明,所提方法可获得新型式的结构拓扑优化方案,并可将结构拓扑优化理论推向工程化应用。     

Optimum shape and topology design using the boundary element method

A procedure is developed for simultaneous shape and topology design optimization of linear elastic two-dimensional continuum structures. An intuitive approach is presented to treat such topological problems whereby material is eliminated from within the structure (by introducing holes at regions of low stress) through a sequence of shape optimization processes. A mathematical programming technique coupled with the boundary element (BE) method of response and sensitivity analyses enables the optimal positioning of these holes plus optimization of the overall structural shape. The analytical derivative BE formulation is explained together with the use of appropriate design velocity fields, and example problems are solved to demonstrate the optimization procedure.     

Structural weight minimization under stress constraints and multiple loading

In this paper we present an approach for structural weight minimization under von Mises stress constraints and multiple load-cases. The minimization problem is solved by using the topological derivative concept, which allows the development of efficient and robust topology optimization algorithms. Since we are dealing with multiple loading, the resulting sensitivity is obtained as a sum of the topological derivatives associated with each load-case. The derived result is used together with a level-set domain representation method to devise a topology design algorithm. Several numerical examples are presented showing the effectiveness of the proposed approach.     

三维有限元六面体网格几何自适应再生成方法

基于几何自适应的六面体网再生成方法,提出了一种以栅格法为基本方法,对有限元网格再划分过程中的网格再划分标准、几何形状的获取以及新旧网格之间物理量的传递等关键问题进行了研究.重点介绍了基于几何自适应的六面体网格再生成算法,首先对旧网格实体模型进行识别,根据实体模型几何特征建立加密源点信息场;然后采用栅格法生成核心网格并对...     

Shape optimization and its extension to topological design based on isogeometric analysis

In most of structural optimization approaches, finite element method (FEM) has been employed for structural response analysis and sensitivity calculation. However, the approaches generally suffer certain drawbacks. In shape optimization, cumbersome parameterization of design domain is required and time consuming remeshing task is also necessary. In topology optimization, design results are generally restricted on the initial design space and additional post-processing is required for communication with CAD systems. These drawbacks are due to the use of different mathematical languages in design or geometric modeling and numerical analysis: spline basis functions are used in design and geometric modeling whereas Lagrangian and Hermitian polynomials in analysis. Isogeometric analysis is a very attractive and promising alternative to overcome the limitations resulting from the use of the conventional FEM in structural optimization. In isogeometric analysis, the same spline information such as control points and spline basis functions which represent geometries in CAD systems are also used in numerical analysis. Such unification of the mathematical languages in CAD, analysis and design optimization can resolve the issues mentioned above. In this work, structural shape optimization using isogeometric analysis is studied on 2D and shell problems. The proposed framework is extended to topology optimization using trimming techniques. New inner fronts are introduced by trimming spline curves in topology optimization. Trimmed surface analysis which was recently proposed to analyze arbitrary complex topology problems is employed for topology optimization. Some benchmarking problems in shape and topology optimization are treated using the proposed approach.     

基于双向演化的局部水平集组合算法用于拓扑优化

本文基于双向渐进结构优化(BESO)方法和局部水平集方法(LLSM),提出局部水平集组合算法,在拓扑优化中实现双向演化。为改进LLSM的孔洞成核能力,新算法以所提离散水平集函数为节点设计变量,拓扑导数为灵敏度,按照BESO优化准则进行双向演化得到稳定拓扑解。然后通过迭代求解距离正则化方程(DRE)来组合LLSM获得最终拓扑。在LLSM中,DRE用来代替重生成步骤,并构造条件稳定差分格式求解DRE。最终给出典型实例验证所提算法的收敛性和数值稳定性。     

Second order topological sensitivity analysis

The topological derivative provides the sensitivity of a given cost function with respect to the insertion of a hole at an arbitrary point of the domain. Classically, this derivative comes from the second term of the topological asymptotic expansion, dealing only with infinitesimal holes. However, for practical applications, we need to insert holes of finite size. Therefore, we consider one more term in the expansion which is defined as the second order topological derivative. In order to present these ideas, in this work we apply the topological-shape sensitivity method as a systematic approach to calculate first as well as second order topological derivative for the Poisson’s equations, taking the total potential energy as cost function and the state equation as constraint. Furthermore, we also study the effects of different boundary conditions on the hole: Neumann and Dirichlet (both homogeneous). Finally, we present some numerical experiments showing the influence of the second order topological derivative in the topological asymptotic expansion, which has two main features: it allows us to deal with hole of finite size and provides a better descent direction in optimization process.