Chapter 12 矩阵位移法.ppt

2. 总刚的带宽与存储方式 结构总刚度矩阵有许多零元素，大型结构的总刚矩阵尤其明显。 稀疏矩阵：具有大量零元素的矩阵 总刚中的非零元素通常集中在主对角线附近的斜带形区域内，为 带状矩阵 在带状矩阵中，每行从主对角线元素起到该行最外一个非零元素止所包含的元素个数，称为该行的带宽。 某行带宽＝该行结点位移分量号－最小相关结点位移分量号＋1 所有各行带宽中的最大值称为矩阵的最大带宽 最大带宽＝相关结点位移分量号的最大差值+1 满阵存储：编程时，可以把总刚的全部元素存储在计算机中起来 编程简单。 等带宽存储：只存储下半带（上半带）在最大带宽范围内的元素 最大带宽愈大，存储量愈大。编程有些复杂。 结点编号时，应该尽量使相关结点编号的最大差值为最小。 变带宽存储：对于对称带状矩阵，每行只存储下半带在该行带宽内的元素。编程比较复杂。 3. 支承条件的引入 前面所述引入支承条件的办法，是删去零位移相应的行和列，修改原始刚度矩阵。 刚度矩阵的阶数降低，但总刚原来的行列编号也要改变，对于编程不便。 实际编程中，常采用“乘大数法”、“置大数法”和“划零置一法”来施加支承条件（边界条件） a)置大数法 原始刚度方程 设结点位移分量?j等于已知值Cj，将总刚中的主元素Kjj置换为一个充分大的数N（如N=1020）,把结点荷载向量中对应的分量Pj置换为NCj 相对充分小，计算时可略去 b)划零置一法 将总刚中的主元素Kjj换为1，j行和j列的其他元素均改为零；结点荷载向量中Pj改为Cj，其余分量Pi改为Pi-KijCj 设 得 4. 铰结点的处理 方法1：刚架中有铰结点时，处理方法如位移法，不把铰结点的转角作为基本未知量。 采用具有铰接端的单元刚度矩阵。 方法2：将各铰结端的转角均作为基本未知量求解。 优点：单元类型统一，编程简单，通用性强。 在铰接点处，各杆端线位移一样，但转角不同，铰结点处转角未知量不只一个，在结点位移向量中要把这些转角未知量都加进去 第十二章 矩阵位移法 §12-1 概述 力法和位移法— 传统结构分析方法，依靠手算，可解决的结构简单。 计算机提供了强大的分析计算手段，使包括结构力学在内的众多学科发生了巨大变化。 许多过去靠手算难以解决，甚至是不可想象的结构问题，由于以计算机为基础的计算分析技术的出现，都迎刃而解。 在结构力学中，适合电算的计算方法主要为结构矩阵分析。 其特点：使用矩阵数学工具，适合编程及机器运算。 杆件结构的矩阵分析，也称杆件有限元法。 一、结构离散化 把结构分解为有限个单元。对于杆件结构，一般一个杆件可作为一个单元。对单元进行分析：考察单元的的内力和位移之间关系。 二、组装单元 化整为零 把各单元集合成原来结构，要求各单元满足平衡条件和几何条件：如支撑条件和结点变形连续条件。对结构整体分析 化零为整 根据所选未知量的不同，结构矩阵分析有矩阵位移法（刚度法）和矩阵力法（柔度法）。 杆件有限元法的主要内容： 矩阵位移法的程序简单，通用性强，运用最广。 矩阵位移法很早提出，对手算很笨，很长时间没有进展。 计算机出现后，由于它特别适合电算，一段时间发展非常迅速，现在很成熟。 §12-1 概述 单元分析－建立杆端力和杆端位移的关系。 除了转角位移方程外，还考虑轴向位移，以矩阵形式表示。 局部坐标系 根据胡克定律和表10－1，得出某杆端位移分量为1，其余分量为零时，各杆端力分量，然后叠加 杆端力为 杆端位移为 杆端力和位移正负号： 轴力 以 正向为正；剪力 以 正向为正； 弯距逆时针为正 单元刚度方程 简写为 其中 单元刚度矩阵为 单元刚度矩阵的系数－刚度系数的物理意义： 当其所在列对应的杆端位移分量＝1，其余分量＝0时，所引起的其所在行对应的杆端力分量的数值。 单元刚度矩阵性质： 1）对称性 反力互等定理 2）奇异性 行列式＝0，不存在逆矩阵 杆端位移－杆端力； 由杆端力不能反求杆端位移 物理意义：自由单元，两端无约束，有刚体位移 l i j e 特例：平面桁架中的杆件 只有轴力，无剪力、弯距 单元刚度方程为 §12-3 单元刚度矩阵的坐标转换 ? ? ? 单元刚度矩阵，在局部坐标系下建立。 对整体结构，各杆方位各不相同，局部坐标系也各不相同 需设定一个统一的坐标系：整体坐标系 ? 在整体坐标系中，杆端力和杆端位移向量为 杆端力在局部坐标系和整体坐标系中的关系为 或者 同样 坐标转换矩阵为 性质：正交矩阵 其中，整体坐标系中的单元刚度矩阵 按单元的始末结点i,j进行分块 单元刚度矩阵的具体形式见式12-23 平面桁架杆件，两端只承受轴力，在整体坐标系下的杆端力和杆端位移向量为 ? 坐标转换矩阵 §