单侧接触约束下桁架拓扑优化的若干结果.pdfVIP

单侧接触约束下桁架拓扑优化的若干结果 ! 杨德庆 刘正兴 宣兆成 （上海交通大学建筑工程与力学学院，上海， 200030 ） （大连理工大学机械工程学院，大连， 116024 ） 摘 要 对考虑应力和无摩擦单侧接触位移约束的桁架拓扑优化问题进行了研究 . 采用解 析方法，通过对典型算例不同初始接触间隙下的拓扑优化解的讨论，发现这一优化问题中存在 解不唯一现象 . 研究也表明，接触状态对结构拓扑有一定影响，不同初始接触间隙可能导致不同 的拓扑结构 . 对利用数值方法求解该问题进行了研究 . 为克服结构拓扑优化中奇异最优解问题， 采用 ! 松弛方法建立桁架拓扑优化模型列式，并利用双层优化方法求解本拓扑优化问题 第一 层，通过求解与原接触问题等价的二次规划问题来进行结构接触分析；第二层，利用第一层得到 的位移及应力，采用 ! 松弛方法进行结构拓扑优化 . 关键词 结构，最优化，拓扑，接触，非光滑 ! 引言 桁架结构拓扑优化的研究历史可追溯到 1904 年提出的 MicheII 桁架理论 . 然而真正使 结构拓扑优化具有可操作性的研究工作是 Dorn 等人提出的基结构方法［1］，以及近 40 年来 发展起来的有限元理论 . 在基结构法基础上，目前发展出两类不同的桁架结构拓扑优化模 型，即连续变量拓扑优化模型和离散变量拓扑优化模型 . 对于连续变量拓扑优化模型，现有 方法多以单元截面积、内力或位移等为拓扑设计变量，构造线性规划或非线性规划并求 解［2 ~ 11］ . 在这种模型列式下， Sved 和 Ginos （ 1968 ）首先发现了结构拓扑优化中的奇异最优解 现象 . Kirsch （ 1990 ）和 Cheng （ 1992 ）等对此进行了研究， Cheng 揭示了产生奇异最优解的根本 原因， Cheng （ 1997 ）提出 ! ! 松弛算法解决奇异最优解问题 . 对于离散变量拓扑优化模型，单元 的拓扑（有和无）用整型量 1 和 0 描述，同时要优化单元的尺寸参数 . 这类模型列式中不存在 奇异最优解现象，多采用随机组合优化方法求解，如分枝定界法、模拟退火法和遗传算法 等［12 ~ 15］ . 基于上述两类模型，现有的工作总体上研究了单、多工况下，考虑应力、位移、频率 及局部稳定性等约束的桁架拓扑优化问题 . 接触问题属于边界待定问题，它与结构拓扑优化问题均是力学领域中较难处理的问题 . 考虑接触的拓扑优化问题更是非光滑优化问题，这方面的研究工作目前不多见［17］ . 本文基 于连续变量拓扑优化模型，对考虑应力及无摩擦单侧接触位移约束的桁架拓扑优化问题进 行研究 . 建立了多工况考虑接触约束的桁架拓扑优化模型，并采用双层优化方法求解；对不 同初始接触间隙下结构拓扑的变化也进行研究 . 第 22 卷 第 1 期 2001 年 3 月 固 体 力 学 学 报 ACTA MECHANICA SOLIDA SINICA VoI .22 No.1 March 2001 ! 上海市博士后科研资助计划项目（ H9902006 ） . 1999-09-14 收到第 1 稿， 2000-04-05 收到修改稿 . !应力及无摩擦单侧接触约束下桁架拓扑最优解的不唯一性 图 1 三杆桁架 我们首先研究图 1 所示著名的三杆桁架［8］ . 材料及几何参 数为：弹性模量 E = 1 .0 ，密度 ! = 1 .0 ， I = 1 .0 . 杆 1 和 3 的许用 应力 \ \ = 5 ；杆 2 的许用应力 \ \ = 20 . 三种荷载工况下的荷载 为 P = 40 ， 30 ， 20 ， # = 45 ， 90 ， 130 ， = 1 ， 2 ， 3 . 各杆截面积分 别为 A 1 ， A 2 ， A 3 ； A 点为可能接触节点 . 1968 年 Sved 和 Ginos 正是通过研究这一算例的拓扑优化， 首次发现了多工况应力约束下结构拓扑优化存在奇异最优解 现象 . 原问题共有四种结构拓扑形式 . 研究表明，若不采取特殊措施，用常规数学规划方法是 无法得到原问题的全局最优解的，因此它也成为衡量拓扑优化新算法能否克服奇异最优解 现象的算例之一 . 无接触约束时，原问题的最优解为（ 8 . 0 ， 1 . 5 ， 0 . 0 ）T，最优重量 W ! = 12.814 . 采用位移法对该结构进行分析，得结构刚度阵［ ! ］为 ［ ! ］ = 2 4 A 1 + A 3 A 1 - A 3 A 1 - A 3 A 1 + 2 2A 2 + A [ ] 3 （ 1 ） 外荷载节点力向量 分两种情况： 无接触时， = P cos # P sin # { } ； 接触产生时， = P cos # P