基于预条件Lanczos说你法的结构拓扑修改静态重分析方法.pdfVIP

基于预条件Lanczos 算法的结构拓扑修改静态重分析 方法1 徐涛，程飞，于澜，王永利，修航 吉林大学机械科学与工程学院，吉林长春 (130022) E-mail：chengfeifuture@ 摘 要：本论文提出了一种新的结构拓扑修改静态重分析方法。该方法易于实现并且适用于 拓扑修改中的三种情况：（a ）自由度不变化；（b ）自由度减少；（c ）自由度增加。根据初始 分析结果，对拓扑修改的三种情况分别选取了适当的预条件子，应用Lanczos 迭代算法得到 拓扑修改后结构响应。将此算法与直接Lanczos 算法分别应用于数值算例，比较了预条件前 后迭代收敛速度。应用本文提出的方法计算成本显著降低，计算结果显示了该算法的正确性 和有效性。 关键词：拓扑修改；结构静态重分析；Lanczos 算法；有限元系统 中图分类号：O0302 1．引言 [1, 2] 结构优化是近年来工程技术研究领域的重要课题，特别是拓扑优化方法的研究 。拓 扑优化的目的是在给定的约束条件下，寻找最优材料分布来达到目标函数（如质量，柔度） 的最小化。在拓扑优化中，随着结构的不断修改，需要求解一系列形如：K X P 的大规 i 模线性方程组，其中K 是在结构拓扑优化第i 步中结构刚度阵，X 是待求结构响应，P 是 i 右端载荷向量。在结构拓扑优化中求解一系列大规模线性方程组已成为计算的瓶颈。结构重 分析技术就是利用初始结构响应，设计高效重分析算法求解修改后结构的新响应而不对新结 构进行完整分析，因此重分析技术可以大大减少计算成本，摆脱了拓扑优化中的繁重计算， 从而加速优化过程。 重分析方法可分为精确方法和近似方法[3, 4] 。精确重分析一般基于Sherman-Morrison and Woodbury (SMW) [5] 公式，对于结构小修改给出修改结构的精确响应。近似重分析方法通常 适用于大中等规模系统的修改，给出结构修改后响应的近似解。国内外许多学者已进行了大 量的研究，文献[5]分析了几种结构修改静态精确重分析方法，并且根据SMW 公式将易于实 现的线性重分析技术扩展到一些非线性重分析问题中。文献[6]提出了一种基于序列矩阵逆 （SMI）的重分析方法,并将它用到优化设计中。文献[7, 8] 提出了著名的拓扑修改重分析的组 合近似方法（CA ），它将缩减基方法与一个级数展开的前若干项结合起来，计算是基于一 个给定点的精确分析结果，每次重分析求解一个小规模线性方程组。文献[9] 提出了一个处 理各种拓扑修改的结构重分析统一的CA 方法。文献[10, 11] 中改进了CA 方法，提出了一个处 理各种拓扑修改的统一的迭代组合近似方法（ICA ），该方法可以处理拓扑结构大修改。但 是 CA 方法和 ICA 方法都需要引进增量刚度矩阵，计算缩减刚度矩阵，存储基向量（特别 是对于大修改，为了得到适合的精度需要较多基向量），而这些操作都增加了计算成本。 在拓扑优化中，拓扑修改后的结构按自由度划分可以分为三种类型：（a ）自由度不变 1本课题得到高等学校博士学科点专项科研基金（项目编号：20050183018 ）和吉林大学“985 汽车工程’项 目资助 - 1 - 化；（b ）自由度减少；（c ）自由度增加。对于这三种拓扑修改，迭代法可以被用来求解结 构修改后系统[3, 12, 15] 。求解线性系统的迭代法中，Krylov 子空间法是应用最广泛的方法之 一[13] ，其