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

维普资讯 第 29卷 第 2期 固 体 力 学 学 报 Vo1．29NO．2 2008年 6月 CHINESEJOURNALOF SOLID MECHANICS June2008 基于半隐式格式的水平集法连续体结构形状和拓扑优化方法 罗俊召 陈立平 罗 震 (华 中科技大学机械学 院CAD中心 ，武汉 ，430074) 摘 要 提 出一种改进的基于水平集方法 的结构拓扑优化方法．将半隐式的加性算子分裂方法 (AOS)引入 传统的水平集方程求解 ，使得原来 的Hamilton—Jacobi偏微分方程摆脱 了差分法 中CFL条件对时间步长的严格限 制 ，差分格式变得高效且稳定．在求解水平集方程过程 中不用再对高维的水平集函数进行耗时的周期性的初始化 ， 这样解决了传统水平集方法在优化过程 中不能生成新孔 的问题．通过典型算例验证 了该文算法 的有效性． 关键词 拓扑优化 ，水平集法 ，加性算子分裂 ，有限差分法 0 弓l言 为了获得最优拓扑，必须经过很多迭代步的演化．同 水平集法最初是 Osher和 Sethian_l在研究 曲 时为了保证差分格式的求解精度 ，需要在优化过程 线 (或曲面)以曲率相关的速度演化 时用于描述曲线 中不断地初始化水平集函数．全局的重新初始化以 (或曲面)的演化过程 的．其基本思想就是将 曲线或 及步长小导致在优化过程 中不能产生新 的孔洞 ，只 曲面隐式地表达成为一个高维水平集函数的零水平 能通过边界合并生成孔 ，使得优化结果严重依赖初 集．由于水平集是一种描述边界的隐式表达法 ，它能 始值．为解决不能生成新孔和设计结果严重依赖初 够灵活地同时描述结构边界形状和拓扑变化．Seth— 始值的问题，Burger等_7将拓扑导数引入到水平集 Jan和 Wiegmann_2最早将水平集法应用到结构优 法中，使其在优化过程中能够生成新孔 ，但是该种方 化设计 ，形状和拓扑的变化根据设计边界上 的等效 法很难控制拓扑形状导数和形状导数进行有效的转 应力来演化．之后 wang等 和 Allaire等 ]扩展 了 换．Wei等 ]采用 AOS方法求解水平集方程的一个 水平集拓扑优化 的方法 ，将水平集方法和形状导数 变形形式 ，也能够在优化过程中生成新孔 ，但是该方 结合起来 ，通过结构边界 的法向速度建立 了Hamil— 法稳定性较差．本文采用 A0S方法来求解水平集 ton-Jacobi方程 同形状 导数之 间的关系 ，采用逆风 拓扑优化 中的水平集方程，因为 AOS格式为无条 策略来求解水平集方程 ，实现结构边界的演化直至 件稳定，因此极大提高了优化过程的收敛速度，同时 获得最优解． 不用在优化过程 中对水平集 函数进行频繁初始化 ， 水平集拓扑优化方法是一种基于隐式边界描述 使得优化过程能够产生新孔 ，减轻了优化结果对初 的方法 ，同均匀化方法和 SIMP方法 相比具有多 始值的依赖性． 个优点 ，例如 ：形状和拓扑优化可以同时进行并且保 1 基于水平集法的优化模型 持结构边界 的光滑 ，没有棋盘格现象和中间密度单 如 图1(a)，对于一个给定的具有光滑边界 的设 元等．但同时也存在收敛速度慢 ，优化结果受初始值 计域 ，假设存在一个函数 声(z)满足 ： 影响大 ，优化过程中只能实现孔洞的合并而不能产 f声(z)0∞ Vz∈~／a2／ 生新孔等缺点． 声(z)