应力与位移约束下基于HCA算法的结构拓扑优化研究.pdfVIP

应力和位移约束下基于 HCA 算法的 结构拓扑优化研究 1) 周珍珍 陈建桥 (华中科技大学力学系，武汉 430074) 摘要：混合元胞自动机(HCA)算法[1]是一种将元胞自动机原理与有限元分析结合起来的方法。利用元胞自动机离散性、 自适应、自进化的特征，与有限元分析方法相结合，可以解决优化过程中复杂结构的力学分析问题，为有效进行结 构的拓扑优化设计提供新的分析思路和技术手段。该方法不仅极大地减少了计算量，提高了求解效率，而且克服了 数值不稳定问题，如棋盘格式问题和网格依赖性问题；不仅适用于二维连续体结构拓扑优化，而且可以直接应用于 三维连续体结构拓扑优化中。 关于HCA算法进行结构拓扑优化设计，国外有些初步的研究成果，但技术环节的实现、复杂问题的求解有待进 一步发展。本文基于此HCA算法，更深入研究了考虑应力约束和位移约束时的改进方法(引入罚函数方法) ，算例结果 较直观地证明了，提出基于HCA算法的改进方案，对解决结构拓扑优化设计问题是有效的。对于其他更复杂的应用 条件，如何建立适合问题的局部控制规则，算法实现技术等问题，有待进一步的研究和发展。 关键词： 元胞自动机法，有限元法，结构拓扑优化，局部控制规则，应力约束，位移约束 引言  目前，最优化设计理论和方法在机械结构设计中得到了深入的研究和广泛的应用。根据设计变 [2,3] 量的类型，结构优化设计可以分为三个层次：(1)尺寸优化；(2)形状优化；(3)拓扑与布局优化 。 其中，拓扑优化是结构设计领域的关键和难题，因为结构拓扑形态在很大意义上，基本决定了结构 的功能、载荷、约束和材料配置等的适用范围。即用现有结构、材料、力学解析法(有限元法等)和 结构优化方法，由于表达结构拓扑设计的自由度庞大，只能在相对结构初始拓扑形态变化较小的范 围内，实现结构的拓扑形态的再设计和强度、应力分布等的再校核[4] 。为此，近年来基于元胞自动 机(CA) [5, 6]方法的结构拓扑优化设计的研究，对解决这一设计难题，开创了一个新的方向，具有十 [4] 分重要的学术价值和应用前景 。 元胞自动机(Cellular Automaton，简称CA)，是定义在一个由具有离散、有限状态的元胞组成的 元胞空间上，并按照一定局部规则，在离散的时间维上演化的动力学系统[7] 。CA 以其组成单元的简 单规则性，单元之间作用的局部性和信息处理的高度并行性，并表现出复杂的全局性等特点而备受 [8] 关注，成为探索复杂系统的一种有效工具 。 [1] HCA算法 简单易懂、收敛性好、计算效率高、克服了数值不稳定问题，但目前只有些初步的研究， 对于理论模型的构建、复杂问题的求解、技术环节的实现等方面还有待进一步发展。总的来说，HCA 算法比较实用且有较大的改进空间，因此本文其将作为研究的基础，结合罚函数法等对其加以改进 和完善，希望寻找到一种更精确实用的拓扑优化设计方法。 1 混合元胞自动机算法(HCA) 1.1 HCA 算法优化模型 变密度法是一种常用的拓扑优化方法[9,10] 。变密度法中常用的插值模型是固体各向同性惩罚微 结构模型(Solid isotropic microstructure with penalization，简称SIMP模型)[11]，设材料模型为： 1) E-mail: zhou_zhenzh@163.com E (x ) xp E (p≥1) (1) i i i 0