弹性力学—第六—用有限单元法解平面问题.pptVIP

第六章 用有限元法解平面问题 胡 衡 武汉大学土木建筑工程学院 有限元法的概念 有限元法的基本量 有限元法的基本量 有限元法的量 有限元法的基本方程 有限元法的基本方程 有限元法的求解步骤(1) 有限元法的求解步骤(2) 有限元法的求解步骤(3) 有限元法的求解步骤(4) 有限元法的求解步骤(5) 位移模式 单元的应变列阵和应力列阵（1） 单元的应变列阵和应力列阵（2） 单元的应变列阵和应力列阵（3） 单元的应变列阵和应力列阵（4） 单元的结点力列阵与劲度矩阵（1） 单元的结点力列阵与劲度矩阵（2） 单元的结点力列阵与劲度矩阵（3） 单元的结点力列阵与劲度矩阵（4）- 举例 单元的结点力列阵与劲度矩阵（5） 载荷向结点移置（1） 载荷向结点移置（2） 载荷向结点移置（3） 载荷向结点移置（4） 载荷向结点移置（5） 载荷向结点移置（6） 载荷向结点移置（7） 结构整体分析（1） 结构整体分析（2） 结构整体分析（3） 结构整体分析（4） 有限元计算简例(1) 有限元计算简例(2)— 单元刚度矩阵的意义 有限元计算简例(3) 结构整体分析（4） 结构整体分析（5） 结构整体分析（6） 结构整体分析（7）- 计算单元刚度矩阵 结构整体分析（8）- 得到整体刚度矩阵 结构整体分析（9）- 在整体刚度矩阵中引入边界条件 结构整体分析（10）- 结点载荷 结构整体分析（11） 结构整体分析（12）- 求解 结构整体分析（13）- 整理结果求应力 结构整体分析（14）- 整理结果求应力 结构整体分析（14）- 整理结果求应力 结构整体分析（15） 有限单元法的单元划分（1） 有限单元法的单元划分（2） 计算结果的整理（1） 计算结果的整理（2） 根据虚功原理，我们也推导了结点力与结点位移的关系： 对于 i 点, 一个单元上的结点力为: i 点的力平衡要求围绕 i 点的各单元产生的结点力与各单元分配到 i 点的结点载荷相等。 综上所述，i 点的平衡方程为。 将所有结点的平衡方程写在一起，就得到了整个结构的结点平衡方程： 这里 称为整体劲度矩阵， 是整体结点载荷列阵， 是整体结点位移列阵。 其中： n为结点数, 下标1, 2, 3…是结点在整体结构中的编号。注：i, j, m 是结点在单元中的局部编号。 整体刚度矩阵的一个元素Krs 就是由按整体编码的，同下标的单元劲度矩阵元素的叠加得到的： x y 2m 2m 2m 2m 2N/m 2N/m 1 2 3 4 5 6 i--3 I j--1 m--2 1 2 3 5 6 I II III IV 4 kij 表示由于j结点沿x或y方向上的位移而引起的在i结点的沿x或y方向结点力。 i I II IV j m III i j m 1 2 3 5 6 I II III IV 4 I II III IV 整体刚度矩阵中的一个元素Krs表示，由于s结点的位移而引起的在r结点的结点力，下面举例说明： i I II IV j m III i j m 1 2 3 5 6 I II III IV 4 整体刚度矩阵中的K23, 表示3号结点的位移由于结构劲度在2号结点上引起的结点力，因此，它将由两部分组成，一部分来自I单元，一部分来自III单元。观察以下表格，可知 。 I III 整体劲度矩阵的获得过程：首先将整体刚度全部充零，然后逐个单元的建立单元劲度矩阵，最后根据结点的局部编码与整体编码的关系，将所有单元劲度矩阵的每一个子矩阵都叠加到整体劲度矩阵中去。 E E 1 2 3 5 6 I II III IV 4 需求解的结点还剩： 因此关于这六个零分量的六个平衡方程不用建立，须将整体刚度矩阵的第1，3，7，8，10，12以及同序列的各列去掉。最后得到： 1 2 3 5 6 I II III IV 4 i I II IV j m III i j m 不考虑位移边界条件的情况下： 1N/m 1 2 3 5 6 I II III IV 4 i I II IV j m III i j m 1N/m 考虑边界条件，去掉零位移对应的结点载荷： 至此，结构的整体结构刚度矩阵以及整体载荷列阵都已经得到，下面根据平衡方程 求解结点位移。 求应力所需关系式如下： 1 2 3 5 6 I II III IV 4 i I II IV j m III i j m I, II，IV II 1 2 3 5 6 I II III IV 4 对于单元III： III i j m 有限元法的求解步骤： 划分有限元， 利用已知的结点坐标以及结构的物理特性写出单元劲度矩阵， 利用整体编码与局部编码的关系写