基于水平集方法及COMSOL的弹性结构拓扑优化.pdfVIP

基于水平集方法及 COMSOL的弹性结构拓扑优化 张绍华，向家伟 (桂林电子科技大学机电工程学院，广西 桂林 541004) 摘 要：针对现有弹性结构拓扑优化方法计算效率差、结构边界不光滑等缺点，提出了一种基于水平集方法及COMSOL 的弹性结构拓扑优化方法。优化方法利用水平集函数作为设计变量，结构的整体柔度最小为 目标函数 ，实体材料所 占的 体积比为约束条件，结合COMSOL偏微分方程 (PDE)模块中反应扩散方程有限元求解弹性结构拓扑优化问题。该方法 与密度惩罚法相比，得到了光滑的结构边界 ；与传统的水平集法相比，不用有限差分求解复杂的Hamilton-Jacobi方程 ， 摆脱了柯朗一弗里德里希斯 一列维 (CFL)条件的限制 ，提高了计算效率。数值算例结果表明该方法能够提高计算效率 并获得光滑的结构边界。 关键词：水平集方法；拓扑优化；COMSOL；有限元 中图分类号：TB12 文献标识码：B 文章编号：1672—545X(2015)01—0086—03 材料的有效利用一直是人类追求的目标，也是 出：要保证稳定性 ，数值解的依赖区域必须全部包含 许多研究领域不变的话题，并伴随着结构优化理论 解析解的依赖区域 ，这个条件称为 CFL条件。由于 和方法的产生而发展。拓扑优化是结构优化领域 内 CFL条件的限制，导致计算效率很低。 相对较新且发展很快的一种方法，因此近年来拓扑优 基于上述原因，提出了一种基于水平集方法及 化方法引起了广泛关注。对于弹性结构拓扑优化，目 COMSOL软件的弹性结构拓扑优化方法 ，该方法利 前采用较多的是均匀化方法lll、密度惩罚法 (SIMP)[2-3] 用水平集函数作为设计变量 ，结构的整体柔度最小 和水平集法 。均匀化方法是将微结构引入优化 问 为 目标函数，实体材料所 占的体积比为约束条件 ，结 题，采用微结构的尺寸参数作为拓扑设计变量，通过 合COMSOL软件中偏微分方程 (PDE)模块的反应扩 对其进行调整来改变材料的力学性能 ，从而得到结 散方程 ，以有限元求解反应扩散方程替代有限差分 构的最优拓扑，由于均匀化方法引入 了大量的设计 求解复杂的Hamilton—Jacobi方程，摆脱有限差分求 变量，造成计算量非常大。密度惩罚法把材料的密度 解水平集方程时 CFL条件的限制，从而提高计算效 作为拓扑设计变量 ，通过对中间密度进行惩罚，使优 率并得到光滑的结构边界。 化结果尽可能具有非 0即 1的密度分布 。密度惩罚 法虽然具有概念简单 ，实现容易等特点，但其优化结 1 水平集方法 果存在灰度单元 、棋盘格式、锯齿形边界等现象。水 平集法最初 Osher和Sethian[8]在研究曲线以曲率相 1．1基本原理 关的速度演化时用于描述 曲线的演化过程时提出来 水平集方法的基本思想是把材料边界隐性的表 的，2000年 Sethian和Wiegmann9[]首次把水平集法引 达为高一维水平集函数的零水平集，所有有关结构 入结构优化领域。传统的水平集方法是通过差分求 拓扑和结构边界的信息都体现在零水平集之中，由 解一个复杂的Hamilton—Jacobi方程来演化水平集函 满足优化条件的要求驱动零水平集的移动。对于一个 数的，即使采用本质无振荡的高分辨率格式，在数值 具有光滑边界的设计域，假设存在一个隐函数 () 求解过程中不可避免可能会遇到数值不稳定性 问 满足： 题 ，这是由于用显式差分方法数值求解双曲型方程 ． 1 ()o ∈Q。(材料)1 时网格剖分所造成的。Courant—Frideriche—I~wy[曾 ()=0 ∈n(边界)}I (1) 经分析了差分格式的稳定性和差分解 的收敛性 ，指