有限元中的子结构法.pdfVIP

J7 ，．寻纬枸獬 ，，寸 维普资讯 有 限 元 中 的 子 结 构 法 郑鑫元 他 、2f 在对汽轮机零部件实施有限元分析时， 排求解方程的先后次序，使得计算起来更为 经常遇到一个问题就是有限元形成的线代数 方便，这个就是子结构的摄原始想法。 方程阶数都很高，可以达l到几千阶方程组， 例_如对成组叶片组或进行强度分析或进行振 动分析，单只叶片已是近千阶方程，5片一组 鄢就要接近 5千阶，不要说微机上计算，就是 1 在大型的计算机 上实施也要化去相当机时。 就 目前我厂计算机设备主要是微型机，因此 计算机的r土J存容量和速度是找 们 求 解 问题 时，首先考虑的问题 。奉文阐述大型结构 问 题的静态和动态的乎结构辫法，提供了计算 圈 1 机上具体实施的步骤。 图I所示为一个任意结构。我们人为地 一 、 结构静力计算中子结构 划分丁一些区域④，③，………⑥，其中， 对于大型的结构有限元，由于所对短的 ①、②、③、④为子结构物 (包括界面 )， 方程组曲阶数通常商选数千阶，在求解时， ⑤、⑥、⑦、⑧称为界面。通常在子结构中 当然也可用波阵法或矩阵分块解法。还可采 再划分若干单元，同时在界面上具有相应的 用子结佝解法： ‘先对要计算分析的物体划 结点。倘若对未知量编号时，按区域①、② 成若干个子结构，遵一解出备子缔构，最终 ③………@的次序进行显然可看出，在求解 汇总得 一个较低阶次的蝗体方程，这个方 方程时，各子结构可分别独立进行消元，它 程 由于阶数较低，一般来说比较容易求解。 们只与界面有美。根据渡阵法原理，可 以先 我们知道，直接法解有 限元所形成的代数方 消元再组集。于是，对每个子 结构进 行消 程组，jf谐 是按存分块三角解法，还是艘l阵 元并盔策后，渡宽中就只剩下界面上的未知 法，它们都是 以有限元所形成的方程组是带 量。因此，最后要求解的整体方程的阶数就 状为先决条件，其核心是高斯消去法，它们 大大地缩减了，这就是子结构解题的基本想 两者之间的区别就在于方程组的未知数消去 法。因此，在某种意义上也可以称为广义的 的次序上。三角分解，按预先排定的结点序 波阵 。 号依次来i卜解消元的，丽波阵法是按单元序 对于一个物体，设所划分的子结构总数 号所对应构结点的先 焉次序消元。因此，它 为N，在未知量编号对，从第一 个 结构 开 们实质上都是变带宽高斯消去法，其特点是 始，在编完 N个子结构后再编界面，则组集 基于方程组的带状分布。从这个思想出发， 后的整体方程系数矩阵为一个加 边对 角块 对任何结拇，我们都可以人为地任意地去安 型 I — — 50 —— 维普资讯 、 K iI