灵敏度过滤技术在拓扑优化中的应用研究.pdfVIP

灵敏度过滤技术在拓扑优化中的应用研究 1 2 赵丹宁 ，丁繁繁 1 河海大学工程力学系，江苏南京（210098 ） 2 苏州二建集团钢结构分公司，江苏苏州 (215000) E-mail：zd163.com 摘 要：棋盘格现象是结构拓扑优化中的普遍现象。本文介绍了连续体结构拓扑优化中的棋 盘格现象以及产生原因，并对目前解决这一问题的各种方法进行了分析研究。本文借鉴灵敏 度再分配抑制棋盘格的过滤方法，给出基于一阶灵敏度再分配的算法流程及程序框图，编制 相应的计算机程序，减弱结构拓扑优化中的棋盘格现象，并通过算例进行比较分析，证明此 方法的可行性。 关键词：棋盘格，拓扑优化，灵敏度过滤技术，连续体 1. 引 言 结构拓扑是当前结构优化领域内最具有挑战性的课题之一，也是当前结构优化设计领域 的研究热点。由于其难度大，连续体结构拓扑优化理论的研究直到近十多年才取得飞速发展。 结构拓扑优化主要是指在特定的设计区域内，满足给定的约束条件，寻求目标函数最小或最 大的最优的材料分布形式，从而确定结构中空洞的数量、位置及其连接形式[1][2] 。 棋盘格式是连续体结构拓扑优化中常见的一种现象。对于棋盘格产生的原因，目前尚没 有本质上的定论。在拓扑优化中，棋盘格式的出现与所采用的材质设计变量无关，即不论是 采用均匀化设计方法还是采用密度法，均会出现棋盘格式。在早期结构拓扑优化研究中，有 些学者误人为这事一种较为较优的理想的结构。但从生产角度来讲，不但加工困难，成本也 很高，不符合实际应用。从材料本身讲，也不符合材料的最优分布原则，棋盘格使结构形状 的提取和制造变得困难，因为设计者必须在元素与孔洞之间作取舍。因而这并不是我们需要 的最优解。 近代一些学者针对此现象提出了自己的看法，Jog 和Haber 认为是密度变量和唯一变量 的混合变分问题，并证明在特定的组合下是可以避免棋盘格的出现，其中Haber 提出的周长 约束法，此法是通过限制结构的周长来抑制棋盘格式的出现但周长的约束值事先难以确定， 只能通过试验方法得到，给实际应用带来困难[3][4] ；Kikuchi 认为棋盘格式的出现是由于相邻 单元间设计变量值变化过大引起的，基于此理论Kikuchi 提出了一种超级单元法，但这种方 法并不能彻底消除棋盘格式[5] 。 解决棋盘格式最简单的方法是采用高阶有限单元代替低阶有限单元，增加单元自由度， 如采用多节点单元。但在单元数目相同的情况下，分析高节点单元结构的时间要比分析低节 点单元结构的长得多，尤其在三维结构上更加明显。 上述这些方法都能不同程度的解决棋盘格现象，但是不能完全解决，而是在某些程度上 减弱了棋盘格现象。本文中使用的方法是参考Li 、Qing 和Steven 等人在二维结构上减弱棋 盘格单元的应用[6] 。该方法是在拓扑优化过程中对作为材料删除判据的灵敏度进行再分配 （如本文作为材料删除判据的灵敏度就为单元应力），有效地改善由于数值不稳定造成的灵 敏度在单元之间的高低分布现象，抑制材料删除后的拓扑构形中出现孔洞单元和实体单元呈 周期性分布，从而有效地抑制棋盘格式，此方法同样适用于三维结构。 - 1 - 2. 灵敏度再分配的抑制法及其计算原理 2.1 二维结构单元权重的一阶再分配计算原理 灵敏度再分配抑制棋盘格单元的方法是重新计算每个单元的应力值，使得当邻近单元 应力值较高或较低时，此单元的应力值也会随之提高或减少。 这种方法首先按式(1)计算连接到第k 个节点的单元应力平均值α , k