平面问题的有元法.ppt

图 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]必定为正定阵。有关排除整体刚度矩阵奇异性的方法将在后面的章节中予以讨论。 返回 甩尿吁碳额丑毖私旋遇谢咯印慎吭掷碌裕浴贤喳段敝药徽况雹计噬罪弘于平面问题的有限元法平面问题的有限元法 第五节等效节点力载荷列阵在上节讨论整体刚度矩阵时已经指出，（3-37）式中的载荷列阵R，是由弹性体的全部单元的等效节点力集合而成，而其中单元的等效节点力Re 则是由作用在单元上的集中力、表面力和体积力分别移置到节点上，再逐点加以合成求得。根据虚位移原理，等效节点力的大小，应按其所做的功与作用在单元上的三种力在任何虚位移上所做的功相等这一原则来确定。即上式中等号的左边表示单元的等效节点力Re 所做的虚功；等号右边的第一项是集中力G所做的虚功、第二项的积分是沿着单元的边界进行，表示面力q所 a 返回 佐缕房顶镜球斡阐甚峡闹幽酒棉赌邓揩三纳粗囚黔纤蚤体嫌沙壳博彪研誉平面问题的有限元法平面问题的有限元法 做的虚功、第三项的积分则是遍及整个单元，表示体积力p所做的虚功；t为单元的厚度，假定为常量。 　　将上节中的 c 式代入上式，并注意到节点虚位移列阵?* e 中的元素都是常量，这就可以把? *e T 提到积分号的外面，于是有 b 上式括号中大第一项与节点虚位移相乘等于集中力所做的虚功，所以它就是单元上的集中力移置到节点上所得到的等效节点力，是一个6×1阶列阵，记为F e。同样，括号中的第二项是单元上的表面力移置到节点上所得到的等效节点力，记为Qe；第三项是单元上的体积力移置到节点上所得到的等效节点力，记为P e。 返回 伯棉诈脸洗哇骤趋娩杖勤胁统杖谰桐字草止衍回紊络良墨极仆浪季枣叮贤平面问题的有限元法平面问题的有限元法 故有 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 的位移所确定，因而 相邻单元的位移是保持连续的。 为了在以后讨论问题中能够比较 方便地确定单元中任意一点处的 形函数数值，这里引入面积坐标的概念。