基于人工材料密度的新型拓扑优化理论和算法研究-机械工程学报.PDFVIP

第 40 卷第 12 期 机 械 工 程 学 报 Vol.40 No.12 2004 年 12 月 CHINESE JOURNAL OF MECHANICAL ENGINEERING Dec. 2 0 0 4 基于人工材料密度的新型拓扑 优化理论和算法研究* 左孔天 陈立平 钟毅芳 王书亭 张云清 (华中科技大学机械工程学院 武汉 430074) 摘要：对基于人工材料密度的拓扑优化理论作了深入分析；推导出了一种基于人工材料密度的拓扑优化准则算法， 将该算法应用于拓扑优化的计算中，并给出了基于人工材料密度的拓扑优化迭代分析流程；提出了一种卷积因子 方法，用于消除拓扑优化计算结果中易出现的棋盘格式和多孔材料现象；通过数值计算验证了理论和算法的有效 性；分析讨论了不同优化参数对拓扑优化计算结果的影响。 关键词：人工材料密度 拓扑优化 优化准则 卷积因子法 数值计算 中图分类号：TH11 1 基于人工材料密度的拓扑优化理论 0 前言* 和算法 拓扑优化是结构优化中相对复杂的一类问题， 从 20 世纪 80 年代末期以来，产生了许多拓扑 近年来在工业、军事、航空航天等领域中得到了一 优化理论和方法，如均匀化理论、进化结构优化理 些初步应用。然而无论是理论上还是应用上，拓扑 论、人工材料密度优化理论、基因优化算法以及拓 优化还存在许多问题需要解决，是结构优化领域研 扑函数描述方法等，其中均匀化方法和人工材料密 究的热点和难点问题[1,2] 。 度方法最具有代表性。 均匀化方法[3~6]是一种经典的拓扑优化方法，它 1.1 均匀化拓扑优化理论和算法的不足之处 在数学和力学理论上最为严密，但其均匀化弹性 均匀化方法是 M.P.Bendsoe 和 N.Kikuchi 等[3] 张量的求解非常复杂，并且微单元的最佳形状和 于 1988 年提出的一种拓扑优化方法，该方法求解拓 方向难以确定。另外其计算结果中容易产生棋盘 扑优化问题分为三个步骤：①以复合材料力学为基 格式和多孔材料现象[7~9] ，优化结果的工程可制造 础，由最小势能原理出发并结合均匀化理论的微元 性差。 H 进化结构优化方法[10~12]符合工程直观性和直 体假设，求得结构的均匀化弹性张量Eijkl (x ) 。②采 用一定的优化方法更新设计变量，并保证计算过程 觉性，但其难以应用于复杂目标函数和多约束情况 的收敛性和稳定性。③用有限元方法求解结构位移 下的拓扑优化问题。也有学者用遗传算法对结构进 场和其他性能指标。其中尤以均匀化弹性张量的求 行拓扑优化，然而这种算法难以进行大规模的计算 解最为复杂。均匀化弹性张量的求