有限元法基础.ppt

约束处理过程—乘大数法 关联约束方程 将关联约束加到有限元平衡方程对应自由度行，第3行或第4行，这里取第4行 约束处理过程—乘大数法 如果关联约束方程可以消除有限元平衡方程的奇异性，约束后的方程组就存在唯一的一组解。 约束后的方程组的系数矩阵是非对称的。 利用初等变换方法将系数矩阵变换成对称的 约束处理过程—乘大数法 关联约束方程 再乘以系数-k 加到约束后方程组的第3行 约束处理过程—乘大数法 约束处理过程—乘大数法 强制位移边界条件 约束后的方程组简化为 约束处理过程—乘大数法 固定位移边界条件 k = 0，C = 0 约束后的方程组简化为 约束处理过程—乘大数法 固定位移和强制位移边界条件的乘大数约束处理相对比较简单，而且它们的系数矩阵约束后是相同的，只是简单地将方程组系数矩阵中要约束自由度的对角线元素加上一个相对大数A即可 乘大数法的叫法并不十分准确，应该叫加大数法更贴切 乘大数和加大数的效果是一样的 约束处理过程—两种方法比较 赋0赋1法在约束处理过程中是严格精确的，而乘大数法是一种近似约束处理方法，它的精度取决于所乘大数A值 两种方法都可以消除有限元平衡方程的奇异性，得到符合实际边界条件的唯一一组解。但两种方法还是有很大的区别 约束处理过程—两种方法比较 采用乘大数法约束处理后的有限元平衡方程在求解时可能造成解的失真，大数A值越大可能解的偏差会越大，而赋0赋1法就不会出现类似的问题，它在约束过程和求解过程都是精确的 乘大数法相对于赋0赋1法在约束处理过程上简单一些 约束处理过程—两种方法比较 赋0赋1法实际上是将关联位移约束方程代入到有限元平衡方程中的，是代入法。而乘大数是将占绝对优势的关联位移约束方程合并到有限元平衡方程中的，是罚方法，计算误差来自于合并过程，计算精度取决于关联位移约束方程的优势大小 商业软件中，位移边界条件的约束处理都采用赋0赋1法，乘大数很少被采用主要原因是它是一种近似方法，而且大数的大小也不好确定，有时还会造成求解失败 单元平衡方程组装过程 2 1 ① U3 U4 U2 U1 1’ u1 u2 u3 u4 2’ 单元平衡方程组装过程 3 ② U6 1’ 1 U2 U1 2’ u1 u2 U5 u3 u4 单元平衡方程组装过程 单元①矩阵扩阶 2 1 ① U3 U4 U2 U1 1’ u1 u2 u3 u4 2’ 单元平衡方程组装过程 3 ② U6 1’ 1 U2 U1 2’ u1 u2 U5 u3 u4 单元②矩阵扩阶 单元平衡方程组装过程 组装总刚 单元平衡方程组装过程 总体平衡方程 单元平衡方程组装过程 组装单元 ① 单元平衡方程组装过程 再组装单元 ② 总体刚度方程 K 称为总体刚度矩阵 U 称为位移向量 F 称为载荷向量 总体刚度矩阵K的特性 对称性 奇异性 稀疏性 非零元素带状分布 约束处理过程 为什么要约束处理 ? 总体平衡方程组是奇异的 消除无限制的刚体运动 使总体平衡方程组存在唯一一组解 约束处理过程—边界条件 边界条件分类 力(载荷)边界条件 位移边界条件 集中载荷力 表面分布力 自重力 热交换引起的温度载荷 固定位移约束 强制位移约束 关联位移约束 约束处理过程—模型简化 y x 约束处理过程—模型简化 y x x y 约束处理过程—约束方程 1 2 3 4 5 6 7 8 9 10 11 12 y x 约束处理过程—约束处理方法 位移约束处理方法 赋0赋1法 乘大数法 约束处理过程—赋0赋1法 位移约束处理例子 关联约束方程 约束处理过程—赋0赋1法 赋0赋1法 约束处理过程—赋0赋1法 有6个方程，5个未知数，如果约束方程可以消除有限元平衡方程组的奇异性，则取任意5个方程联立求解，都会得到方程组的唯一一组解。 系数矩阵由原来的对称的变成了非对称的，这对于大规模有限元方程组求解是十分不利的，采用相同的求解方法，在求解时间和矩阵存贮容量方面都增加了一倍。 约束处理过程—赋0赋1法 如何解决 ? 可以发现, 约束处理过程实际上是在系数矩阵上做了一次列初等变换。 为了保持平衡方程组系数矩阵的对称性，对变换后的系数矩阵再做一次相同的行初等变换。 约束处理过程—赋0赋1法 具体做法： 第4行乘以系数k加到第3行，并去掉第4行。 为了保持系数矩阵的阶数，将第4行的所有元素赋0，在其对角线元素位置赋1。 即所谓赋0赋1法。 约束处理过程—赋0赋1法 经过初等变换，方程组的系数矩阵仍然保持对称性 初等变换不会改变方程组的解 约束后的方程组可以求得5个未知数 通过关联约束方程回代求解U4 约束处理过程—赋0赋1法