有限元第2章-杆件结构的有限元法_直接刚度法.ppt

类似，可以得到节点2处的力关系。 由此可得单元力在局部坐标系和总体坐标系件的 关系是 或 其中： ， 为矢量变换矩阵 节点位移也是矢量，其在局部坐标系和总体坐标 系中的变换关系与力的变换关系相同 根据力和位移的变换关系，就可以把局部坐标系 下的单元刚度矩阵，转换成总体坐标系下的单元 刚度矩阵。 （4）总体坐标系下的单元刚度矩阵 将 和 代入局部坐标 系下的单元刚度矩阵 中 两端同乘 ，并注意到 ，有 即 其中： 写成显式为 其中： （5）总体刚度矩阵 将单元刚度矩阵合成，得到 求解出节点位移后，即可计算杆的内力： 式中： ——杆单元的节点总体编号。 [例] 图示杆件系统。节点1、3固定，节点2处受 力 。所有材料相同，弹性模量E，截面 积A。求各杆受力。 解：按照有限元步骤求解。 （1）单元刚度矩阵 求解关键：确定各杆的方向。 单元1（节点1-2）： 单元2 （节点1-3）： 单元3 （节点2-3）： 由此可的单元刚度矩阵： 单元1： 单元2： 单元3： （2）合成 整个系统有3个节点，每个节点有两 个自由度，故总刚矩阵为 阶 将实际值代入 因为 ，所以最后只有两个方 程 改写为 由此可以求出： 各杆受力 如果 ，则 二 刚阵存储和节点排列 铰支杆系统每个节点有两个自由度，系统节点数 为n时，总体刚度矩阵的阶数为 。 如果节点数n=1000，那么总体刚度矩阵中包含的 元素为 个。 为了提高计算机的运行速度和节省计算机的内存 空间，有必要研究总体刚度矩阵中元素的分布情 况。 对于铰支杆系统，每一个单元只有两个节点，单 元刚度矩阵中的元素数为 个。 设单元节点的总体编号为i和j，将单元刚度矩阵 合成为总体刚度矩阵后，单元刚度矩阵中的16个 元素在总体刚度矩阵中的位置为 总体刚度矩阵的特点： （1）对称；（2）稀疏；（3）带状分布。 刚度矩阵中的元素以主对角线元素为对称，成带 状分布。带内有些是零元素，带外全部是零元素。 为了节约存储空间，可以只存带内一半的元素。 带的宽度与单元两节点的总体编号的差值有关。 刚度矩阵的最大半带宽 =节点自由度数*（单元中节点最大编号差+1） 差值越小，带越窄。为此，应尽量缩小同一单元 内节点编号差值的最大值。 为了做到这一点，节点编号的方式是有讲究的。 [例]两弹簧系统。 解：采用图示编号方案，总刚矩阵为 采用如下编号方案 总体刚度矩阵为 图示铰支杆系统，应该怎么样编号？最大半带宽 等于多少？需要存储多少个元素？ 结论： 为了减少半带宽，应顺着结构窄的方向编号，沿着结构长的方向运动。 大型问题的等带宽存储量，只有满阵存储的10%~20%。 小结 （1）本章介绍的方法称为直接刚度法。 （2）总刚矩阵对称、稀疏，主对角线上的元素为正，未引入边界条件前是奇异的。 （3）总体刚度矩阵由单元刚度矩阵叠加而成。叠加的原则： a.节点位移唯一；b.节点力的平衡条件。不与所关心的节点相联系的单元，对于该节点的力平衡没有贡献。 （4）可以利用对称性只存半带宽内的元素（含主对角线上的元素）。 （5）半带宽的尺寸与同一单元中节点编号差的最大值有关。注意节点编号的方式。 第二章 杆件结构的有限元法 2-1 引言 杆的特点：只能承受拉力或压力 变形和位移的关系： 对于杆件系统： 式中 刚度矩阵 节点力向量 节点位移向量 杆件系统存在一个整体刚度（矩阵）。知道了这个刚度就可以由系统外力求出各节点的位移。这就是本小结应掌握的一个重要概念。 如何求出杆系的刚度矩阵呢？ 2-2 弹簧系统的刚度矩阵 一 单个弹簧的刚度矩阵 弹簧系统是由节点（铰链）把一系列弹簧连接起来的系统。 如果把弹簧系统拆开，可以得到弹簧和铰链。 从其中取出任意一个弹簧，弹簧的两端会受到来自节点（铰链）的作用力。 单个弹簧两端节点的位移和力向量的关