三维连续体结构形状优化设计ξ.pdfVIP

第 16 卷　第 2 期 应 用 力 学 学 报 . 16　 . 2 V o l N o 1999 年 6 月 CH INESE JOURNAL OF APPL IED M ECHAN ICS Jun. 1999 三维连续体结构形状优化设计 王　备　　　隋允康 ( ) ( ) 中国工程物理研究院　610003 　 大连理工大学　大连　116023 摘　　要 借助于一种名为“函数识别器”的工具, 采用“二级控制”法, 实现了用设计变量控 制网格变动的目的。除此之外, 又在单刚对设计变量的导数计算、目标函数的敏度计 算方面, 提出了新的解决方案。最后给出两个算例, 验证了以上方法是切实可行的。 关键词: 形状; 优化函数识别器; 网格变动控制; 敏度 1　引　言 连续体结构形状优化是通过改变结构的边界形状来优化结构的性能, 它以控制结构几何 形状的几何参数为设计变量, 这使得形状优化同截面优化相比在实现过程中会遇到更多的障 碍。概括起来主要包括以下三个方面问题: ) 1 　如何实现几何模型与有限元模型之间的关联, 即用设计变量控制网格变动。 2) 　有限元网格的自动生成与精化。 ) 3 　状态变量对设计变量的敏度分析。 近二十多年来, 国内外学者在上述三个方面进行了积极的研究工作, 并取得了一定的进 展[ 13 ] 。本文从形式语言的角度, 提出将函数关系表达式引入优化模型的思想, 研制了一个“函 [ 4 ] 数识别器” , 它可以接受任意带变量的函数表达式输入, 改变了以往程序只能接受数值输入 [ 4 ] 的状况。然后借助于这种工具, 提出了“二级控制”法 , 实现了几何模型与有限元模型之间的 关联。相应地, 该方法的采用, 又在单刚对设计变量的导数计算、目标函数的敏度计算方面提出 了新的解决方案。 2　几何模型与有限元模型之间的关联 在形状优化过程中, 设计变量同它控制的结构边界是不断变化的, 优化程序必须能根据设 ( ) 中国工程物理研究院科学技术基金资助课题 93040 来稿日期: 修回日期: © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 128 应 用 力 学 学 报 第 16 卷 计变量值自动调整有限元网格, 以保持网格与结构几何形状的一致性。这种制约关系要求在设 计变量与有限元模型之间建立起联系。本文首先用映射法[ 5, 6 ] 自动生成网格, 然后用二级控制 法建立起设计变量对有限元网格的控制关系。其原理如下: 2 1　映射法自动生成网格 先将结构人工划分为几个大单元, 然后再对这些大单元划分网格, 即在大单元映射到参数 空间中后所得到的立方体区域里划分网格, 生成节点和单元, 再由映射关系计算出节点在原空 间中的真实坐标。这里所说的大单元就称之为设计单元[ 5, 6 ]。设计单元上控制其形状的点称为 关键点。这些关键点不一定是有限元网格中的节点。关键点又可分为可动和不可动, 可动关键 点的位置是可变的, 形状优化的过程就是搜寻可动关键点的最佳位置。 2