结构布局修改静力重分析的预条件共轭梯度法

【摘要】：在结构设计和优化过程中，有时需要对结构进行几十次甚至上百次的修改，而每作一次修改都执行一次完全的有限元分析其计算成本是巨大的。为了减少计算成本，以不直接求解结构修改后的隐式分析方程，而利用初始分析的相关信息设计高效算法来求解修改后结构响应的重分析方法，日益受到人们的重视并得到飞速发展。 传统的结构修改主要侧重于参数修改。随着结构设计与优化技术的不断发展，人们发现结构布局的设计与优化比布局固定的设计与优化更具有节约潜力。 然而，由于在布局修改中，结构的有限元模型发生变化，其自由度数量的改变将导致分析方程组的结构和阶数发生变化。目前大部分重分析方法都是针对结构自由度没有变化的情形。 发展适用于结构布局修改的重分析方法是具有巨大挑战性和经济价值的。 预条件技术是科学与工程计算中复杂问题高效求解器的最关键组成部分。 所谓预条件就是将一个难处理的问题转化成另一个解能够被迅速逼近的问题。 本文主要应用预条件技术解决结构布局修改的静力重分析问题。 论文主要由五部分组成。第一章绪论，主要介绍了重分析问题的工程背景和国内外的研究状况。 第二章首先回顾了解结构修改后自由度未发生变化的重分析问题的组合逼近方法和预条件共轭梯度法，然后从计算复杂度和收敛精度等方面将二者做了比较，指出了预条件共轭梯度法的优越性。第三、四、五章是作者的主要研究工作(即摘要的1－3部分)，分别建立了求解结构自由度减少、增加和一般布局修改的静力重分析问题的预条件共轭梯度法。 结点自由度减少情形的静力重分析方法 给定一个初始结点自由度为 的结构，相应的刚度矩阵记为 ，荷载向量记为 ，对应的位移 可由求解如下平衡方程 WP=93 (1) 来获得。由初始分析，刚度矩阵 可写成如下的Cholesky分解形式 (2) 式中 是上三角矩阵。 假设增加和删除一些单元并且从初始结构中删除一些结点后，结构有 个自由度。记 与 分别是结构修改后 阶刚度矩阵和荷载向量。则修改后的分析方程组可以写成如下形式： (3) 式中 是位移向量。我们选取 (4) 为预条件矩阵，再利用预条件共轭梯度法来求解。在方程(4)中， 是 阶投影矩阵，变换 通过除去m-维向量 中对应于被删除结点的分量而得到n-维向量 ； 是将n-维向量 转化为m-维向量 ， ， 在保留自由度上的分量与 的对应分量相同，而在被删除自由度上的分量为零。 本文建立了求解该问题的预条件共轭梯度法的有效辅助，研究了其计算复杂性。 结点自由度增加情形的静力重分析方法 在结点自由度为 的初始设计基础上，假设增加和删除一些单元并增加一些新的结点，这使得结构的自由度数目增加。修改后的分析方程可以写为如下形式： (5) WP=94 式中 表示在修改设计中新增加结点的自由度数目， 对应新增结点的子刚度矩阵， 为联系初始结点和新增结点的子刚度矩阵。方程(5)可以分解为以下两个线性系统： (6) (7) 由方程(7)解出 ，然后代入方程(6)可以得到如下方程： . (8) 式中 , , . (9) 这样就把一个 维的问题化为一个 维的问题。本文提出的第一个方法是以 为预条件矩阵应用预条件共轭梯度方法求解(8)。在求得 后代入(7)式即可求得 。 本文提出的第二个方法是取 (10) 为预条件矩阵，再利用预条件共轭梯度方法来求解方程(5)。 本文证明了 作为预条件矩阵的必要性质——对称正定性，建立了求解该问题的两种预条件共轭梯度法的有效辅助并研究了它们的计算复杂度。 结构一般布局修改的静力重分析方法 在结点自由度为 的初始设计基础上，假设增加和删除一些单元，删除了 个已有自由度并增加了 个新的自由度，修改后结构的分析方程为 (11) WP=95 式中 为 阶刚度矩阵。修改后的刚度矩阵可以分块的形式写为： (12) 下角标中小写字母表示子刚度矩阵的阶数, 其中 。 本文构造 (13) 为预条件矩阵，用预条件共轭梯度方法求解方程(11)。在方程(13)中， 是 阶投影矩阵，变换 通过除去m0-维向量 中对应于被删结点的分量而得到m-维向量 ； 将m-维向量 变换为m0-维向量 ， 中对应保留自由度上分量与 向量的对应分量相同，而在被删除自由度上的分量为零。 本文证明了 作为预条件矩阵的必要性质——对称正定性，建立了求解该问题的预条件共轭梯度法的有效辅助并研究了该算法的计算复杂度。 在本文提出的预条件方法中，预条件算子的构造主要是基于修改前的刚度矩阵，这样，预条件算子的构造是廉价的且预条件后的系统很容易求解。预条件方法的主要优点是计算量小，需要的存储空间小，收敛速度快，算法容易辅助并能自适应地监视逼近解的精度，因而该方法能够大大地节省计算费用。