基于半隐式格式水平集法连续体结构形状和拓扑优化方法.pdfVIP

第29卷第2期 固体力学学报 V01．29No．2 CHINESE OFSOLIDMECHANICS 2008年6月 JOURNAL June2008 基于半隐式格式的水平集法连续体结构形状和拓扑优化方法。 罗俊召” 陈立平 罗 震 (华中科技大学机械学院cAD中心。武汉。430074) 摘 要 提出一种改进的基于水平集方法的结构拓扑优化方法．将半隐式的加性算子分裂方法(AOS)引入 制，差分格式变得高效且稳定．在求解水平集方程过程中不用再对高维的水平集函数进行耗时的周期性的初始化， 这样解决了传统水平集方法在优化过程中不能生成新孔的问题．通过典型算例验证了该文算法的有效性． 关键词拓扑优化，水平集法，加性算子分裂，有限差分法 O 引言 水平集法最初是0sher和Sethian[11在研究曲 线(或曲面)以曲率相关的速度演化时用于描述曲线 (或曲面)的演化过程的．其基本思想就是将曲线或 曲面隐式地表达成为一个高维水平集函数的零水平 集．由于水平集是一种描述边界的隐式表达法，它能 够灵活地同时描述结构边界形状和拓扑变化．Seth— ian和wiegmann[21最早将水平集法应用到结构优 化设计，形状和拓扑的变化根据设计边界上的等效 应力来演化．之后Wang等[33和Allaire等[41扩展了 水平集拓扑优化的方法，将水平集方法和形状导数 结合起来，通过结构边界的法向速度建立了Hami卜 ton—Jacobi方程同形状导数之间的关系，采用逆风 策略来求解水平集方程，实现结构边界的演化直至 获得最优解． 水平集拓扑优化方法是一种基于隐式边界描述 的方法，同均匀化方法和SIMP方法[5’6j相比具有多 个优点，例如：形状和拓扑优化可以同时进行并且保 持结构边界的光滑，没有棋盘格现象和中间密度单 元等．但同时也存在收敛速度慢，优化结果受初始值 影响大，优化过程中只能实现孔洞的合并而不能产 f拳(z)O甘Vz∈Q／an 生新孔等缺点． 庐(z)一。甘Vz∈an (1) 在传统的水平集拓扑优化方法中[3“]，水平集 【庐(z)O∞Vz∈D／n 方程通过显式的迎风差分格式来求解，为了保证数 值求解的稳定性，每次推进的时间步长必须满足 CFL条件，在空间中表现为每次演化的步长由最小 的网格尺寸决定．因而，对于一般密度的离散网格， *2007一05—10收到第1稿，2008一04—22收到修改稿． **通讯作者．Ejmail：junzha01029@gmail．com． 万方数据 · · 176 固体力学学报 其中n表示设计域，“表示位移域空间U中的一个 位移域，C为弹性张量，e为应变张量，p是体积力，厂 是应用在边界an的部分边界an。上的牵引力． 对上面的优化问题，有多种方法可以求解。本文 中采用增广拉格朗日的方法，将体积约束作为一个 惩罚项施加在目标函数上，得到如下的公式： (a)设计域 (b)水平集模型 H(拳)dn一 Min．：_厂(“，声)一．，。(“，声)+A(I (a)Designdomin (b)Levelsetm