拓扑相关荷载作用下结构形状%2f拓扑优化的水平集方法.pdfVIP

郭旭等：拓扑相关荷载作用下结构形状／拓扑优化的水平集方法 627 拓扑相关荷载作用下结构形状／拓扑优化的 水平集方法+ 郭旭 赵康 工业装各结构分析国家重点实验室。大连理工大学力学系，大连116023 摘要 在本文中，我们发展了一种利用水平集技术求解拓扑相关荷载作用下结构拓扑优化问题的数 值方法。优化问题的目标函数是结构的整体柔顺性。约束条件为对于结构的体积限制。可以看到，将 水平集技术iJIA,优化问题，可以有效的解决荷载与拓扑相关给优化问题所带来的求解困难。文末的数 值算例表明了本文所提出的优化方法在处理此类问题时所具有的独到的优越性。 关键字拓扑优化：拓扑相关荷载；灵敏度分析；水平集 1引言 拓扑优化是结构优化领域内最具挑战性的研究课题。较之于尺寸优化和形状优化，拓扑优化能 够带来更大的优化收益。因此在近20年来得到了飞速的发展。对于连续体结构的拓扑优化，已经有 了大量的研究工作，其中基于微结构思想的均匀化化方法“。，直接借助材料密度进行刚度插值的SIMP 方法81以及利用启发式手段逐渐寻求改进设计的进化算法。1等都是比较有效的求解技术。有关这方面 童心 在以往连续体结构的拓扑优化研究中，往往都假设外荷载与设计无关，也就是不随结构的拓扑 变化而发生改变的。而对于拓扑相关荷载，意即其作用位置，大小和方向都要取决于当前结构拓扑 的荷载作用下的结构拓扑优化设计问题，就笔者所掌握的资料的来看，相关的研究工作还不是很多。 Hammer和Olho田言先研究了在静水压力作用下的二维结构拓扑优化问题。在他们的工作中，采用了 基结构的参数化方式和材料的密度作为拓扑设计变量。为了处理拓扑相关荷载，他们首先根据密度 等值线信息选取一些插值节点，然后利用这些节点用Bezier三次样条构造出一条光滑的荷载作用曲 线。最后再利用等效功原理把外荷载移置到有限单元的节点上。随后他们又将上述工作推广到了三 结构方法框架下研究了类似的拓扑优化问题。他们用模拟的变温荷载来模拟拓扑相关荷载，因而在 一定程度上减少了荷载作用曲线(曲面)提取所带来的额外工作量。 在本文中，我们基于水平集技术，提出了一种求解拓扑相关荷载作用下结构拓扑优化问题的新 方法。通过引入水平集函数，我们以隐含的方式对结构的拓扑和形状作了描述。这样就把拓扑优化 问题转化为了寻求最优水平集函数的问题。由于结构的边界可以用相应的零水平集来描述。因此利 用适当的数学变换技巧，我们可以方便地处理施加在结构上的拓扑相关荷载。文末的数值算例表明 了本文所提出的优化方法在处理此类问题时所具有的独到的优越性． 2结构拓扑优化的水平集方法 水平集方法是由seIhi∞和wie蛐首先引入结构优化领域的。采用水平集方法进行拓扑优化并 ·囡家自然科学基10032030留学归嗣人员启动基金壹助项目 628 郭旭等：拓扑相关荷载作用下结构形状／拓扑优化的水平集方法 使用响应泛函的灵敏度信息构造水平集发展所需要的速度场则是由Allaire以及Wang等人的工作开始 的。采用水平集方法求解结构拓扑优化问题的基本思路是引入一个水平集函数妒=妒(x，f)，然后采 用如下的方式对结构的拓扑形式加以描述。 onS 妒(x，f)=0 inQ ≯(x，f)0 inD／Q 妒(x，r)0 这里Q代表结构实体材料当前所占有的区域，S是其边界。D是一个计算区域，它包含了结构在优 化过程中可能占有的一切构型，参见附图l。 一旦我们用水平集函数庐隐式地描述了结构拓扑，那么结构拓扑优化的过程可以通过求解以下 Hamilton．Jacobi方程来实现。 詈+v．V舢 这里的V是一个适当选择的用来演化水平集函数的速度场，而时间变量则可对应梯