(优)平面问题的有限元法.pptVIP

图 3-6 a 返回 图 3-6 b 半带宽B=（相邻节点号的最大差值D+1）*2 返回 若第r双行的第一个非零元素子矩阵是[Krl ]，则从[Krl ] 到[Krr ]共有（r-l+1）个子矩阵，于是[K]的第2r行从第一个 非零元素到对角元共有2 (r-l+1) 个元素。显然，带状刚度矩 阵的带宽取决于单元网格中相邻节点号码的最大差值D。我 们把半个斜带形区域中各行所具有的非零元素的最大个数 叫做刚度矩阵的半带宽（包括主对角元），用B表示，即 B = 2 (D+1)。 通常的有限元程序，一般都利用刚度矩阵的对称和稀 疏带状的特点，在计算求解中，只存储上半带的元素，即 所谓的半带存储。因此，在划分完有限元网格进行节点编 号时，要采用合理的编码方式，使同一单元中相邻两节点 的号码差尽可能小，以便节省存储空间、提高计算效率。 返回 ⒌ 刚度矩阵[K]是一个奇异矩阵，在排除刚体位移后，它是正定阵。 弹性体在{R}的作用下处于平衡，{R}的分量应该满足三个静力平衡方程。这反映在整体刚度矩阵[K]中就意味着存在三个线性相关的行或列，所以[K]是个奇异阵，不存在逆矩阵。 因 代入（3-30）得 (u) 返回 上式左乘{?}T ，并注意到 (3-13) 式， 在集合过程中将[B]扩充到3×2n阶后，有[B]3×2n {?}e2n ×1 = [B]3 ×2n {?}2n ×1 ，故 (v) 由于弹性矩阵[D]是正定的，且t和?都是正的，所以只有当每个单元中都有{? }=0时，才有 否则 也就是说，当排除了弹性体的刚体位移{? }=0之后，若 {? }?0，则二次型{? }T [K]{? }恒大于零，于是[K]必定为正定阵。有关排除整体刚度矩阵奇异性的方法将在后面的章节中予以讨论。 返回 故有 (g) 另外，由（3-22）可以求得 (h) 利用形函数的这一性质可以证明，相邻单元的位移分别进行线性插值之后，在其公共边上将是连续的。 返回 例如，对图3-3所示的单元 jm和ijn ，具有公共边ij。 这样，不论按哪个单元来计算，根据（3-11）式，公共边ij上的位移均由下式表示 图3-3 由(3-23)式可知，在ij边上 式中 Ni , Nj 的表达形式如（3-23）式所示。 (i) 返回 由此可见，在公共边上的位移 u、v 将完全由公共边上的两个 节点i、j 的位移所确定，因而 相邻单元的位移是保持连续的。 为了在以后讨论问题中能够比较 方便地确定单元中任意一点处的 形函数数值，这里引入面积坐标的概念。 在图3-4所示的三角形单元ijm中， 任意一点P(x , y)的位置可 以用 以下三个比值来确定 图3-4 式中 ?为三角形单元ijm的面积，?i 、?j 、?m 分别是三角形Pjm、Pmi、Pij的面积。这三个比值就叫做P点的面积坐标。 (3-24) 返回 显然这三个面积坐标并不是完全独立的，由于 所以有： 而三角形pjm的面积为： 故有： 返回 类似地有 (3-25) (3-26) 由此可见，前述的三角形常应变单元中的形函数Ni 、Nj 、Nm 就是面积坐标Li 、Lj 、Lm 。 根据面积坐标的定义，我们不难发现，在平行jm边的直线上的所有各点，都有相同的坐标Li ，并且该坐标就等于“该直线至jm边的距离”与“节点i至jm边的距离”之比，图3-4中给出了Li 的一些等值线。 返回 容易看出，单元三个节点的面积坐标分别为 节点 i： Li =1 Lj =0 Lm =0 节点 j： Li =0 Lj =1 Lm =0 节点m： Li =0 Lj =0 Lm =1 不难验证，面积坐标与直角坐标之间存在以下变换关系： (3-27) 返回 当面积坐标的函数对直角坐标求导时，可利用下列公式： (3-28) 返回 求面积坐标的幂函数在三角形单元上的积分时，有 (3-29) 式中 ?、?、? 为整常数。若求面积坐标的幂函数在三角形某一边上的积分值时，则可用下式 (3-30) 式中 l为该边的长度。 返回 一. 单元刚度矩阵 第四节 刚度矩阵 为了推导单元的节点力和节点位移之间的关系，可应用虚位移原理对图3-2中的单元e进行分析。单元e是在等效节点力的作用下处于平衡的，而这种节点力可采用列阵表示为 (a) 假设在单元e中发生有虚