ch10矩阵位移法[李].ppt

由式(10-40)，并将单元①的，单元②、③的代入计算，可得各单元在整体坐标系中的固端力为 由式(10-41)可求出结点2、3上的等效结点荷载为 再由式(10-42)求得综合结点荷载为 于是可得结点外力列向量为 这里需要指出，对于支座结点1，4同样可按式(10-41)和式(10-42)算出等效结点荷载和综合结点荷载。但由于[F1]、[F4]是综合结点荷载与支座反力的代数和，而其中的反力仍是未知的，又由于在引入支承条件时，未知力[F1]、[F4]将被删去或被修改，故在此可不必计算支座结点1、4的等效结点荷载及综合结点荷载。 (5) 引入支承条件，修改原始刚度方程。结构的原始刚度方程为 见p.256 结点1和4为固定端，故已知 在原始刚度矩阵中删去与上述零位移相应的行和列，同时在结点位移列向量和结点外力列向量中删去相应的行，得结构的刚度方程为 见p.257 (6) 解方程。求得未知结点位移为 见p.257 (7) 计算各单元杆端力。按式(10-45)计算。 单元①：见p.258 单元②：见p.258～259 单元③：见p.259 (8) 绘内力图。绘刚架的内力图时，内力正负号规定仍与第三章中的规定相同，而 上述计算出的杆端内力正负号是按本章中的规定求得的。故画 内力图时应特别注意。  刚架的弯矩图如图10-11所示。 例10-3 试用矩阵位移法计算图10-12所示桁架的内力。各杆EA相同。 图10-12 解：（1）对各结点和单元进行编号，并选结构坐标系如图10-12所示。在确定各单元始、末端结点i、j的号码时，我们约定，这样就确定了各单元的局部坐标系。 （2） 计算各单元刚度矩阵的各子块 用矩阵位移法计算桁架的步骤和计算刚架的完全一样。在桁架中，任一结点i只有两个位移分量u i、vi和两个结点外力分量Fxi、Fyi。单刚按式(10-26)计算，根据表10-4中数据，可求得单元②、④、⑦、⑨的刚度矩阵的各子块为 单元③、⑧的刚度矩阵的各子块为 单元①、⑤的刚度矩阵的各子块为 单元⑥、⑩的刚度矩阵的各子块为 表10-4 各单元几何数据 (3) 将各单刚子块对号入座即形成总刚，结构的原始刚度方程为见p.262 (4) 引入支承条件，u1=v1=u2=v2=0，在原始刚度方程中删去与这些零位移相应的行和列，并将已知的结点荷载代入，得结构刚度方程为 见p.263 (5) 解方程，得未知结点位移 见p.264 (6) 计算杆端内力，由式(10-35)及式(10-36)计算。例如单元⑨，有 见p.264 又如单元⑩ 见p.264 同理可求出其他各杆内力，各杆内力值标在图10-12中杆件旁边。 §10-8 几点补充说明 为了便于电算，对实际的计算问题要做以下补充说明。 1．结点位移分量的编号，单元定位向量 单刚子块对号入座形成总刚，必须落实 到每个元素对号入座。单刚子块的两个下 标号码是由单元两端的结点编号确定的， 而每个元素的两个下标号码则由单元两端 的结点位移分量的编号确定。因此，我们 图10-13 不仅要对结点进行编号，而且还要对结点 位移的每个分量进行编号。对图10-13所示刚架，单元、结点和 结点位移分量编号如图所示，它们的对应关系如表10-5所示。 表10-5 各单元始末端结点及结点位移分量编号 对结点位移分量的编号，同时也就是对结点外力分量的编号，两者是一一对应的。 有了结点位移分量的编号，单刚中的每个元素便可按其两个下标号码送到总刚中相应的下标号码送到总刚中相应的行列位置上去 。图10-14表示单元②的单刚元素 的入座位置。一个 平面刚架的一般单元有6个杆端结点位移，依靠这6个号码，其单刚的36个元素才能确定在总刚中的位置，因此这6个号码称为单元定位向量。 图10-13所示，单元②的定位向量为： 图10-14 2．总刚的带宽及存储方式 结构总刚度矩阵为稀疏矩阵。具有大量的零元素，而非零元素通常集中在主对角线附近的斜带形区域内，称为带状矩阵。在带状矩阵中，从主对角元素起至该行（列）最外一个非零元素止所包含的元素个数，称为该行（列）的带宽。由总刚的形成规律可知: 某行（列）的带宽=该行（列）结点位移分量号-最小相关结点位移分量号+1 （10-47） 所有各行（列）带宽中的最大值称为矩阵的最大带