hzf第6章——用有限单元法解平面问题2.ppt

三、集中力的移置： 如集中力G做用于其一边界上如图：先将G分解为Px，Py 后，分别按线段的比例把它们分别移置到i，j两点上。即： 即按静力平衡方法分配。 0 x i j m y §6-7 整体刚度矩阵 讨论了单元的力学特性之后，就可转入结构的整体分析。假设弹性体被划分为N个单元和n个节点，对每个单元按前述方法进行分析计算，便可得到N组形如（6-33）式的方程。将这些方程集合起来，就可得到表征整个弹性体的平衡关系式。为此，我们先引入整个弹性体的节点位移列阵 {?}2n×1 ，它是由各节点位移按节点号码以从小到大的顺序排列组成，即 其中子矩阵 (j) (i =1,2, …, n ) (k) 是节点i的位移分量。 继而再引入整个弹性体的载荷列阵{R}2n×1 ，它是移置到节点上的等效节点载荷依节点号码从小到大的顺序排列组成，即 (l) 其中子矩阵 (i =1,2, …, n ) (m) 是节点i上的等效节点载荷。 现将各单元的节点力列阵{R}e6×1 加以扩充，使之成为2n×1阶列阵 其中，子矩阵 (n) (i, j, m 轮换) (o) 是单元节点i上的等效节点力。 (n)式中的省略号处的元素均为零，矩阵号上面的i, j, m 表示在分块矩阵意义下Ri 所占的列的位置。此处假定了i, j, m 的次序也是从小到大排列的、并且与节点号码的排序一致。 各单元的节点力列阵经过这样的扩充之后就可以进行相加，把全部单元的节点力列阵叠加在一起，便可得到 (l)式所表示的弹性体的载荷列阵，即 这是由于相邻单元公共边内力引起的等效节点力，在叠加过程中必然会全部相互抵消，所以只剩下载荷所引起的等效节点力。 同样，将（6-35）式的六阶方阵[k]加以扩充，使之成为2n阶的方阵 (p) (q) 不难看出，（6-35）式中的2×2阶子矩阵[ki j ]将处于上式中的第i双行、第j双列中。 考虑到[k]扩充以后，除了对应的i, j, m 双行和双列上的九个子矩阵之外，其余元素均为零，故（6-33）式中的单元位移列阵{?}e2n×1 便可用整体的位移列阵{?}2n×1 来替代。这样，（6-33）式可改写为 把上式对N个单元进行求和叠加，得 (r) 上式左边就是弹性体所有单元刚度矩阵的总和，称为弹性体的整体刚度矩阵（或简称为总刚），记为[K]。注意到（6-28）式，有 (6-38) (6-39) 若写成分块矩阵的形式，则 显然，其中的子矩阵为 它是单元刚度矩阵扩充到2n×2n 阶之后，在同一位置上的子矩阵之和。由于(q)式中许多位置上的子矩阵都是零，所以（6-36）式不必对全部单元求和，只有当[krs ]的下标r = s或者属于同一个单元的节点号码时，[krs ]才可能不等于零，否则均为零。 (6-40) * §6-4 应变和应力列阵 有了单元的位移模式，就可以利用平面问题的几何方程 求得应变分量。将 (e) 、(f) 两式代入上式，即得： (g) 可简写成 ?其中 [B] 矩阵叫做单元应变矩阵，可写成分块形式 而子矩阵 由于?和bi 、bj 、bm 、ci 、cj 、cm 等都是常量，所以矩阵[B]中的诸元素都是常量，因而单元中各点的应变分量也都是常量，通常称这种单元为常应变单元。 （i , j , m轮换） (6-15) (6-14) (6-13) 求得应变之后，再将（6-13）式代入物理方程 ，便可推导出以节点位移表示的应力。即 (6-16) (h) (6-17) 令 则 其中 [S]叫做应力矩阵，若写成分块形式，有 (6-18) 对于平面应力问题，弹性矩阵[D]为 (i) 所以，[S]的子矩阵可记为 （i , j , m轮换） (6-19) 对于平面应变问题，只要将 (i) 式中的E换成E/1-?2 ，?换成 ? /1-?，即得到其弹性矩阵 (j) （i , j , m轮换）(6-20) 注意到（6-7）式，则有 (6-21) 由（6-19）、（6-20）式不难看出，[S]中的诸元素都是常量，所以每个单元中的应力分量也是常量。 可见，对于常应变单元，由于所选取的位移模式是线性的，因而其相邻单元将具有不同的应力和应变，即在单元的公共边界上应力和应变的值将会有突变，但位移却是连续的。 一