第三章 平面问题的有限单元法.pptVIP

　　如果任意三角形单元ijk的重心c上受有自重 则按刚体静力等效原理可把W直接移置到i，j，m三个结点上而组成 ： 等效结点力、载荷列阵小结 例:图示三角形单元的i j m的jm边上作用有如图所示线性均匀分布面载荷，三角形单元厚度为t，求结点载荷向量。 解：运用叠加原理 三角形单元jm边长度为： X m (5,6) j (6,3) i (2,2) q1 q2 O Y 习 题 1、 设有三角形单元ijm受自重作用，容重为Pg。 求出其等效结点载荷。 2．设有：三角形单元ijm的ij边受线性分布的法向压力作用，ij两点的压力分别为Pi和Pj试求其等效结点载荷。 3.4 整体刚度矩阵 讨论了单元的力学特性之后，就可转入结构的整体分析。假设弹性体被划分为N个单元和n个结点，对每个单元按前述方法进行分析计算，便可得到N组形如（3-21）式的方程。将这些方程集合起来，就可得到表征整个弹性体的平衡关系式。为此，我们先引入整个弹性体的结点位移列阵 {?}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)式所表示的弹性体的载荷列阵，即 这是由于相邻单元公共边内力引起的等效结点力，在叠 加过程中必然会全部相互抵消，所以只剩下载荷所引起 的等效结点力。 同样，将单元刚度矩阵的六阶方阵[k]加以扩充，使之 成为2n阶的方阵 (p) 返回 (q) 不难看出，其中的2×2阶子矩阵[ki j ]将处于上式中的第i双行、第j双列中。 考虑到[k]扩充以后，除了对应的i, j, m 双行和双列 上的九个子矩阵之外，其余元素均为零，故（3-33）式 中的单元位移列阵{?}e2n×1 便可用整体的位移列阵 {?}2n×1 来替代。这样，（3-33）式可改写为 把上式对N个单元进行求和叠加，得 (r) 上式左边就是弹性体所有单元刚度矩阵的总和， 称为弹性体的整体刚度矩阵（或简称为总刚），记为 [K]。注意到（3-28）式，有 (3-38) (3-39) 若写成分块矩阵的形式，则 返回 显然，其中的子矩阵为 它是单元刚度矩阵扩充到2n×2n 阶之后，在同一 位置上的子矩阵之和。由于(q)式中许多位置上的子矩 阵都是零，所以（3-36）式不必对全部单元求和，只 有当[krs ]的下标r = s或者属于同一个单元的结点号码 时，[krs ]才可能不等于零，否则均为零。 将（3-34）式和 (p) 式代入 (r) 式，便可得到关于 结点位移的所有2n个线性方程，即 [K]{?}={R} (3-41) (3-40) 1 2 3 4 2 1 q 图 3-5 组装总刚[k]的一般规则： 1. 当[krs]中r=s时，该点被哪几个单元所共有，则总刚子矩阵[krs]就是这几个单元的刚度矩阵子矩阵[krs]e的相加。 2. 当[krs]中r s时，若rs边是组合体的内边，则总体刚度矩阵[krs]就是共用该边的两相邻单元单刚子矩阵[krs]e的相加。 3. 当[krs]中r和s不同属于任何单元时，则总体刚度矩阵 [krs]=[0]。 下面，我们考查一个组装总刚的实例： 1. 整体刚度矩阵及载荷列阵的 组集 根据叠加原理，整体结构的