应力流驱动结构拓扑优化.pdf

第 27卷 第 6期 南 京 理 工 大 学 学 报 Vo1．27 No．6 2003年 l2月 JournalofNan．1ingUniversityofScienceandTechnology Dec．2003 应力流驱动 的结构拓扑优化 汤文成 易 红 幸 研 (东南大学机械系 ，南京 210096) 摘要 以应力流作为结构优化的控制 函数，通过建立 以图形和有限元 网格模型相 对应的结构描述 ，采用删除低应力区域结构 的 自适应方式 ，实现结构的拓扑变化和 几何模型生成 ，所提方法具有扩展性和应用价值，完成结构形式遍历，为结构优化 设计作有益的尝试 。 关键词 优化设计 ，结构分析 ．拓扑 分类号 TP315 学科代码 130．15 结构优化设计经历了截面尺寸优化 、几何形状优化 2个阶段 。在结构优化方法 中最为 成熟的是准则法，设计过程 中的结构改进是通过迭代地使用根据该准则而构造 的函数来实 现的，函数 中包含着 由经验决定的参数 ，控制这些参数就可调整收敛速度 。虽然优化设计 和有限元分析各 自形成了完整 、成熟的理论体系，但 由于 目前在结构形状优化过程 中普遍采 用边界节点控制方法 ，从某种程度上讲 ，所得结果不一定是最佳结构 。‘本文 以初始结构的 应力流作为结构优化的控制 函数，通过几何描述和拓扑描述的分离，实现结构的拓扑生成 ， 为结构优化设计作 了有益的尝试 。 1 结构模型描述与拓扑数据处理策略 结构拓扑优化的数学模型(以板壳单元为例)是求 E和 t，使 rrfinⅢ = piA，(E)t，成立， 且满足： ≤ ≤ 。其 中，，z为结构总质量 ； 为第 个单元的材料密度 ；A (E)为以E边 界函数定义的板壳面积；t为板壳厚度-7l为单元数； 和 分别为许用压力的上、下限。 1．1 结构特征表达 定义以下变量或函数作为结构及其有限元模型的特征表达 ： ct 区界函数厂cn = ／ ．。 1 嚣 上的点集。 (2)加权 图。以V，表示单元的节点，E(V ，V，)表示有限元模型中对应于节点 V 至V 的边界(当为其它单元形式时 ，维数可 以扩展 。例如为块体单元 时，可 以扩展到面 函数 收稿 日期：2002—04—25 * 汤文成 男 45岁 教授 第 133期 汤文成 易 红 幸 研 应力流驱动的结构拓扑优化 729 F(E ，E．，E ，… )，那 么对应于 图 G(E Vk)的每个 节点或每 个边 都 有一 个实 数集 {W，(V)或 {W，[E(V，，V )]：其中，、，一分别是模型的点权图和边界权图序列。 (3)节点 V，的度q(V，)：它表示与节点 V，相联系的边界集合，当为板壳单元形式时，这 个集合 中只有边 E(V，，V )；如果包含块体单元 ，则节点的度是该节点与面 F(E，，E ，E ， … )、边 E(V ，V)相关的并集 ： (4)边界 E(V，，V)的度 qE(V，，V，)0基本定义与节点的度概念一样 ： (5)面 S(E，，E)(，J= l，2，…)的度 q；S(E ，E，)：同样 ，它是边界的度定义在三维 空间的推广 ，表示与面 S(E ，E )相关的体单元的集合： 几何上可 以证 明，对于无 向边表述的结构，节点 V 的总度 (或度)等于投射到这个节点 v 的边界的数 目，也即节点 V 的度等于节点V相联系的边界数 。将节点的度 的概念应用到 节点集合 中，设 S= (V，E， )是一有 向结构 ，且 V V是一节点集合 ，那么，这时始节点 属于子集 V ，而终节点不属于该子集 V 的边界数为节点子集V 的出度 ；同样 ，终节点属于 子集 V ，而始节点不属于该子集 V 的边界数称为V 的入度： 1．2 变异结构的有 限元模型处理 本