结构拓扑优化的发展现状及未来

结构拓扑优化的发展现状及未来

王超

中国北方车辆研究所一、历史及发展概况

结构拓扑优化是近20年来从结构优化研究中派生出来的新分支，它在计算结构力学中已经被认为是最富挑战性的一类研究工作。目前有关结构拓扑优化的工程应用研究还很不成熟，在国外处在发展的初期,尤其在国内尚属于起步阶段。1904 年Michell在桁架理论中首次提出了拓扑优化的概念。自1964 年Dorn等人提出基结构法,将数值方法引入拓扑优化领域，拓扑优化研究开始活跃。20 世纪80 年代初，程耿东和N. Olhoff在弹性板的最优厚度分布研究中首次将最优拓扑问题转化为尺寸优化问题，他们开创性的工作引起了众多学者的研究兴趣。1988年Bendsoe和Kikuchi发表的基于均匀化理论的结构拓扑优化设计，开创了连续体结构拓扑优化设计研究的新局面。1993年Xie.Y.M和Steven.G.P 提出了渐进结构优化法。1999年Bendsoe和Sigmund证实了变密度法物理意义的存在性。2002 年罗鹰等提出三角网格进化法，该方法在优化过程中实现了退化和进化的统一，提高了优化效率。

二、拓扑优化的工程背景及基本原理

通常把结构优化按设计变量的类型划分成三个层次：结构尺寸优化、形状优化和拓扑优化。尺寸优化和形状优化已得到充分的发展，但它们存在着不能变更结构拓扑的缺陷。在这样的背景下，人们开始研究拓扑优化。拓扑优化的基本思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求最优材料的分布问题。寻求一个最佳的拓扑结构形式有两种基本的原理：一种是退化原理，另一种是进化原理。退化原理的基本思想是在优化前将结构所有可能杆单元或所有材料都加上，然后构造适当的优化模型，通过一定的优化方法逐步删减那些不必要的结构元素,直至最终得到一个最优化的拓扑结构形式。进化原理的基本思想是把适者生存的生物进化论思想引入结构拓扑优化，它通过模拟适者生存、物竞天择、优胜劣汰等自然机理来获得最优的拓扑结构。

三、结构拓扑优化设计方法

目前常使用的拓扑优化设计方法可以分为两大类：退化法和进化法。

退化法即传统的拓扑优化方法，一般通过求目标函数导数的零点或一系列迭代计算过程求最优的拓扑结构。目前常用于拓扑优化的退化法有基结构方法、均匀化方法、变密度法、变厚度法等。

基结构方法（GSA）的思路是假定对于给定的桁架节点，在每两个节点之间用杆件连结起来得到的结构称为基结构。按照某种规则或约束,将一些不必要的杆件从基本结构中删除，认为最终剩下的构件决定了结构的最佳拓扑。基结构方法更适合于桁架和框架结构的拓扑优化。基结构法是在有限的子空间内寻优，容易丢失最优解，另外还存在组合爆炸、解的奇异性等问题。

均匀化方法（HA）引入微结构的单胞，通过优化计算确定其材料密度分布，并由此得出最优的拓扑结构。均匀化方法主要应用于连续体的拓扑优化设计，它不仅能用于应力约束和位移约束，也能用于频率约束。目前用均匀化方法来进行拓扑优化设计的有一般弹性问题、热传导问题、周期渐进可展曲面问题、非线性热弹性问题、振动问题和骨改造问题等。

变密度法是一种比较流行的力学建模方式，与采用尺寸变量相比，它更能反映拓

扑优化的本质特征。通常，单元密度与弹性模量之间的关系采用人为给出的幂函数规律，例如E E n ⋅=ρρ)(或者E C E r ⋅⋅=ρρ)(，亦有采用ρ的有理分式或者E 1和E 0的组合形式，例如：

E E ⋅-----=ρρρρ3/23/2)1(2)

1(1)(，0013/2010/)1)((1)()(E E E E E E E ρρρ--+-+=。

目前已有研究将变密度法应用于卧式千斤顶以及磁场等拓扑优化设计中，取得了一定的成效。另外，人们对变密度法的不足加以改进,研究出SIMP 法以及能量法则等方法，提高了计算能力及应用水平。

变厚度法采用满应力法，进行有限元分析得到的各单元在节点处的Misses 相当于应力，将围绕每一节点i 的所有单元在节点i 处应力的加权平均值作为接点i 的应力σi 。通过迭代不断改变各节点处的厚度使其应力趋近最大的允许值，达到满应力的设计目的。当节点处的厚度低于事先设定值时，节点即被删除。 进化法是一类全局寻优方法，目前常用于拓扑优化的进化法主要有遗传算法、模拟退火算法和渐进结构优化法等。遗传算法(GA)是一种新型的基于遗传进化机理的寻优技术，它通过选择、交叉、变异等过程可使群体性能趋于最佳，从而获得全局最优解。作为拓扑优化设计方法之一，GA 主要应用于建筑结构优化方面，如桁架结构优化设计、抗震结构智能优化设计等。近几年出现的改进GA 法，如自适应GA 的研究和复合型GA ，使GA 法得到进一步发展。

模拟退火算法（SA ）的思想源于固体退火过程：固体在加温时，其内部粒子随温升变为无序状，内能增大，而徐徐冷却时粒子渐趋有序，在每个温度都达到平衡状态,最后在常温时达到基态，内能减为最小。SA 法是最近十年来才应用于拓扑优化的，应用得比较广泛的是桁架的拓扑优化，最近几年将非线性问题转化为多目标SA 法的研究也取得一定的进展，SA 法适合解决比较复杂的问题。 渐进结构优化法（ESO 法）的基本概念很简单，即通过将无效的或低效的材料一步步去掉，结构也将逐步趋于优化。可用ESO 法优化的结构为桁架、刚架、板壳或三维连续体等，优化的约束条件包括应力、刚度、位移、频率、临界压力及动响应力等。目前ESO 法已被成功应用于热应力优化、热传导优化、触应力优化、可变荷载结构优化以及机床床身的优化等。

四、拓扑优化设计方法比较

1) 用于拓扑优化的优化算法比较

目前拓扑优化设计方法常采用的优化算法主要有两类：优化准则法(OC) 和数学规划法(MP) 。拓扑优化中常用的OC 法是满应力准则法，其最大的特点是收敛速度快，要求重分析的次数一般跟设计变量的数目关系不大。拓扑优化中常用MP 法是移动渐进算法(MMA 法)。对单约束拓扑优化问题，算法一般选用优化准则法；对模型中包含多约束情况或多工况的拓扑优化问题，优化准则法的推导过程复杂，计算和更新拉格朗日乘子比较麻烦，这时应优先选用MP 法。

2) 几种主要退化法比较

基结构方法中的应力的不连续性会引起奇异最优解现象。基结构方法更适合于桁架和框架结构的拓扑优化。均匀化方法不仅能用于应力约束、位移约束，也能用于频率约束，但其设计变量多，敏度计算复杂，优化后的结构常常含有多孔质材料。变密度法基于各向同性材料，程序实现简单，计算效率高。但是有一点需要说明：它是人为假设的，甚至是基于经验的。