结构拓扑优化设计的理论与方法.docxVIP

结构拓扑优化设计的理论与方法 结构拓扑优化设计的理论与方法 2021年第1期 2021年1月 化学工程与装备 Chemical Engineering Equipment 结构拓扑优化设计的理论与方法 张德恒，鹿晓阳 （山东建筑大学 力学研究所，山东 济南 250101） 摘 要：从离散结构拓扑优化和连续体结构拓扑优化两个方面,分别对结构拓扑优化的理论研究、优化方法、实际应用、拓展研究等方面的研究工作进行了介绍。并对拓扑优化中应用的主要数学求解方法进行了简要介绍，指出了它们的适用范围和优缺点。 关键词：拓扑优化；离散结构；连续体结构；数学规划法；最优准则法 结构优化设计是用系统的、目标定向的过程与方法代替传统设计,其目的在于寻求既经济又适用的结构形式,以最少的材料、最低的造价实现结构的最佳性能。结构形式包含了关于尺寸、形状和拓扑等信息,相应的结构优化设计分为三个层次:(1)尺寸优化:优化变量为杆件的横截面尺寸,或板壳的厚度分布；(2)形状优化:优化变量为杆系结构的节点坐标或表示连续体结构外形的变量；(3)拓扑优化:优化变量为杆系结构的节点布局、节点间的连接关系,或连续体结构的开孔数量和位置等拓扑信息。结构的拓扑构形选择恰当与否,决定了产品设计的主要性能,因而,在复杂结构的选型和轻量化设计中,拓扑优化比形状和尺寸优化更有价值。 结构拓扑优化包括离散结构的拓扑优化和连续体结构的拓扑优化。其主要思想是将寻求结构的最优拓扑问题转化为在给定的设计区域内寻求材料的最优分布问题。虽然结构拓扑优化的概念已经提出100多年了, 但直到近几十年才得到快速的发展,在建筑、机械、航空、航天、海洋工程及船舶制造等领域有着广阔的应用前景。 1 离散结构拓扑优化 结构拓扑优化设计研究最早是从桁架类离散体结构开始的, 1904年Michell用解析方法研究了单载荷作用下应力约束的结构设计,提出了桁架结构设计的Michell准则，并将符合Michell 准则的桁架称为Michell 桁架, 这被认为是结构拓扑优化设计理论研究的一个里程碑。但是, Michell提出的桁架理论只能用于单工况并依赖于选择适当的应变场,并且解析方法所涉及的复杂数学推导限制了它在实际工程中的应用。直到1964 年, Dorn、Gomory、Greenberg 等人提出基结构法( ground structure approach),将数值方法引入到结构拓扑优化领域, 此后拓扑优化的研究重新活跃了起来。所谓的基结构就是一个由结构节点、荷载作用点和支承点组成的节点集合, 集合中所有节点之间用杆件相 连所形成的结构。该方法的基本思想是: 从基结构的模型出发, 应用优化算法( 数学规划法或优化准则法) , 按照某种规则或约束, 将一些不必要的杆件从基结构中删除, 保留的杆件决定了结构的最优拓扑。从基结构出发的拓扑优化方法,由于单工况、应力约束下使桁架结构重量最轻的最优拓扑,必定是一个静定结构,因此早期研究者常忽略变形协调条件,以杆件内力作为设计变量、节点平衡方程作为约束条件构造成线性规划问题来处理。这种方法虽然计算效率高,但却无法推广到多工况和考虑位移约束的情况,因为此时结构的最优拓扑往往是超静定的,必须涉及变形协调条件,并采用非线性规划法来求解。Dobbs 和Fetton 使用最速下降法求解多工况应力约束下桁架结构的拓扑优化。Shen 和Schmidt 采用分枝定界法求解在应力和位移两类约束下桁架结构在多工况作用下的最优拓扑。王光远等提出了结构拓扑优化的两相法。Kirsch 针对离散结构的拓扑优化问题,提出了一种两阶段算法, 即首先以赘余内力和杆件截面积为设计变量, 忽略变形协调条件和位移约束, 将问题简化成容易求解的线性规划问题, 求出解的下界；第二阶段考虑全部约束, 在已得到的拓扑下解非线性规划, 得到杆件的截面面积。 针对大型结构的拓扑优化问题, Zhou和Rozvany发展了一种优化准则类算法,,即DCOC 算法, 认为采用这种算法可使准则法求解拓扑优化的能力大为提高。Bendsoe等对于桁架结构质量约束下的最小柔顺性(或最大刚度)设计及其对偶问题进行了较深入的研究,将它表达为仅含结构状态变量(位移、内力)的极大极小问题,设计了高效的处理非光滑优化问题的算法,大大提高了这类优化问题的求解效率。 近年来,一些适合于并行计算且对函数性态要求较低的全局搜索算法, 如遗传算法、神经元网络算法和模拟退火算法等开始被应用于拓扑优化上, 但目前这些方法仅能解决 张德恒：结构拓扑优化设计的理论与方法 较小规模的问题。这方面代表性的工作有: 许素强和夏人伟、Grierson