结构力学之矩阵位移法

《结构力学之矩阵位移法》由会员分享，可在线阅读，更多相关《结构力学之矩阵位移法（18页珍藏版）》请在金锄头文库上搜索。

1、第十二章 矩阵位移法 【例12-1】 图 a所示 连 续 梁 ，EI=常数，只 考 虑 杆 件 的 弯 曲 变 形 。分别用位移法和矩阵位移法计算。 图12-1解：（1）位移法解基本未知量和基本结构的确定用位移法解的基本结构如图c所示。这里我们将结点1处的转角也作为基本未知数，这样本题仅一种基本单元，即两端固定梁。位移法基本方程的建立将上式写成矩阵形式系数项和自由项计算（须绘出单位弯矩图和荷载弯矩图）由图d，结点力矩平衡条件，得，由图e，结点力矩平衡条件，得，由图f，结点力矩平衡条件，得，由图g，结点力矩平衡条件，得，将系数项和自由项代入位移法基本方程，得解方程，得由叠加法绘弯矩图，如图h所示

2、。（2）矩阵位移法解对单元和结点编号（图a）本题只考虑弯曲变形的影响，故连续梁每个结点只有一个角位移未知数。若用后处理法原始结构刚度阵为阶；用先处理法结构刚度阵为阶（已知角位移）。下面采用先处理法来说明矩阵位移法计算过程。单元标准形式为（图b） 求局部坐标系下的单元刚度矩阵求整体坐标下的单元刚度矩阵，因连续梁的局部坐标和整体坐标是一致的，所以有，得（注：本题用先处理法换码） ， ， 按“对号入座”规则集成总刚，得 形成荷载列阵(1) 计算单元固端列阵， （2）将单元固端列阵反号，并按“对号入座”规则送入荷载列阵 （本题结点荷载为零）=将结构刚度矩阵及荷载列阵代入矩阵位移法方程，得解方程，得计算

3、杆端弯矩=得各单元杆端弯矩后，再叠加上一相应简支弯矩图即得各单元弯矩图。将各单元弯矩图组合在一起，得整个结构的弯矩图（图h）。小结：通过本题的计算可看到：（1）基本未知量和基本结构。位移法与矩阵位移法二者都是以结点位移为基本未知量，以单根杆件（单元）为计算对象。位移法为方便计算，有三类杆件；而矩阵位移法只有一类杆件，即两端固定等截面梁。（2）刚度矩阵与荷载列阵的形成。位移法是用单位弯矩图和荷载弯矩图并由结点的平衡条件计算系数项和自由项的，而后形成刚度矩阵与荷载列阵的；而矩阵位移法是以单元杆端刚度元素、单元杆端荷载元素，按“对号入座”规则形成刚度矩阵与荷载列阵的。 矩阵位移法基本方程的建立，归结

4、为两个问题：一是根据结构的几何和弹性性质建立整体刚度矩阵，二是根据受载情况形成整体荷载列阵。（3）有（1）、（2）可知，二者的关系是：“原理同源，作法有别”。因此矩阵位移法不是一个新方法，它是新的计算工具（电子计算机）与传统力学原理（位移法）相结合的产物。【例12-2】试求图a所示结构原始刚度矩阵中的子块 ，已知单元 的整体坐标的单元刚度矩阵如图c所示。 图12-2解：本题每个结点有两个基本位知量（竖向线位移和角位移），如图b所示。单元刚度矩阵为阶（图c）。由图d所示子块形式，的元素应为单元的j端元素（图c右下角子块）与单元i端元素（图c左上角子块乘以2）之和，即 【例12-3】只计弯曲变形时

5、，用先处理法写出结构刚度矩阵 。(设 EI = 1) 图12-3 解：由图d及先处理法结点位移编号图c写出各单元刚度矩阵，并按“对号入座”规则集成整体刚度矩阵。, ,【例12-4】用先处理法写出图a所示结构刚度矩阵,E=常数。不计轴向变形影响。 图12-4解：本题虽然是刚架，但不计轴向变形影响，即每一个结点只有一个角位移未知量。根据图b所示结点位移编号，则整体刚度矩阵为阶。由于每个单元杆端只有角位移未知量，故单元刚度矩阵为阶的连续梁单刚形式。=，=，=，=【例12-5】图示连续梁 ，不计轴向变形 ，EI =常数 ，已知结点位移 。试求单元的杆端力列阵 。 图12-5解：根据图a的约束条件和图b

6、的结点位移编号，已知给出的结点位移是：有，。单元的杆端力列阵为 【例12-6】用矩阵位移法求图a所示桁架各杆内力。单元、的截面面积为A，单元的截面面积为2A，各杆E相同。图12-6解：桁架每个结点两个线位移未知量（图b）。局部坐标系下的单元刚度矩阵为阶，即=，=整体坐标系下的单元刚度矩阵为由图b可知，单元 ，。单元 ，。单元 ，。=，=。 整体刚度矩阵及荷载列阵 =，=矩阵位移法方程= 解方程，得计算各杆轴力=（拉）=（拉）=（压）【例12-7】已知图示桁架的自由结点位移列阵 ，求杆12在局部坐标系中的杆端 力 。设 ，杆12 的横截面积。图12-7 解：，。= 【例12-8】 用位移法和矩阵

7、位移法计算图a所示结构。各杆材料及截面均相同，。要求：（1）不考虑轴向变形影响的位移法解。 （2）考虑轴向变形影响的位移法解。 （3）用矩阵位移法(采用先处理法)解。 图12-8解：（1）不考虑轴向变形影响的位移法求解不考虑轴向变形影响下，仅有结点1处的角位移未知量。位移法的基本方程为 系数和自由项由图b、c得，将系数和自由项由代入位移法的基本方程，并解得 弧度。由叠加法作弯矩图，即。整个结构的弯矩图如图d所示。 （2）考虑轴向变形影响的位移法求解 基本结构如图e所示。位移法的基本方程为系数和自由项计算由图f：，由图g：， 由图h：由图c：，将系数和自由项由代入位移法的基本方程，并解得，弧度考

8、虑轴向变形影响的结构弯矩图如图i所示（剪力图和轴力图未画出）。（3）用矩阵位移法(采用先处理法)解用矩阵位移法求解时，单元和结点编号如图j所示。采用先处理法时其整体刚度矩阵为阶。两单元对应的整体编码如下图所示。 按“对号入座”规则集成结构刚度矩阵=注：（1）单元局部坐标与整体坐标一致，所以有。（2）单元局部坐标与整体坐标的夹角，须进行坐标变换，即。运算的结果是将中相关元素作行列交换。另外当局部坐标与整体坐标的夹角时，我们也可直接在整体坐标系下进行对换，如图k所示。按先后再转角的次序，则可直接在局部坐标的单元上标注相应的整体编码，本题就是采用这一方法。注意到坐标进行了轴交换，变号，故副系数须反号

9、。见本题中单元中送入结构刚度矩阵的元宵和。 荷载列阵的集成。方法一是按及进行。另一作法是，由,于是有=将结构刚度矩阵和荷载列阵基本方程，得与前位移法解得的相同结果，即=同样得结构弯矩图如图i所示（剪力图和轴力图未画出）。【例12-9】试 求 用 矩 阵 位 移 法 求解 图 a所 示 结 构 时 ，结 点 2 的 综 合 结 点 荷 载 列 阵 。解：刚架每个结点有三个基本未知量（），同时也有三个方向结点荷载项。 图12-9（1）结点2的直接结点荷载： （2）结点2的等效结点荷载涉及到单元、及的2端的固端力（见图c、d、e）。按式（）应首先应计算局部坐标系下的固端反力，而后进行坐标变换得整体坐

10、标系下单元固端反力，再“按对号如座”规则反其符号集成。这里我们直接根据图c、d、e求出整体坐标系下的单元固端反力。由图b及d、e、c得=，=，=结 点 2 的 综 合 结 点 荷 载 列 阵为。【例12-10】 试 用 先 处 理 法 写 出 图 a所 示 结 构 刚 度 矩 阵 。各 杆 杆 长 均 为 ，EI = 常 数 ，自 由 结 点 位 移 分 量 的 编 号 如 图 示 。图12-10 解：单元与整体坐标一致。而单元、按图b所示整体坐标系下来进行换码（注意到坐标进行了轴交换，变号，故副系数须反号），而后按下图“对号入座”规则集成总刚。=【例12-11】 用 先 处 理 法 求 图

11、a所 示 刚 架 的 结 构 刚 度 矩 阵 ，略去轴向 变 形 影响。图12-11解：由图b的位移编号可知，横梁各结点仅有一个向的水平位移，其变形如图c的所示（这就是“手算”），按“对号入座”规则集成总刚（这就是“机算”）=用经典位移法解时，其系数。【例12-12】按先处理法计算图a所示结构的刚度矩阵。各杆长度为 ，EA、EI 均为相同 。 图12-12 解：单元、结点及位移编号入图b所示。作为理解画出了结点位移的变形图，如图c、d及e所示（这就是“手算”）。按下图“对号入座”规则集成总刚（这就是“机算”）。= 【例12-13】图示刚架只考虑弯曲变形 ，按先处理法求在荷载和支座位移共同作 用

12、下的结点荷载列阵。已知各杆 。 图12-13解：图b为结点、单元编号，单元固端反力如图c所示 ，是由支座位移产生的。=，=，=+=【例12-14】 图 示 刚 架 各 杆 ，结点6有支 座 的 水 平 位 移 ，竖 向 位 移 ，忽 略 轴 向 变 形 ，已 求 得 结 点 位 移 为： 。求 单 元 的 杆 端 内 力 。图12-14解：本题有两各特点：（1） 不计轴向变形影响，单元刚度为阶，如图b所示，不需坐标变换。（2） 结点6的支座移动只有对单元有影响，将它作为单 元 杆端位移值，则有所以 【例12-15】对图示刚架的结点和单元进行编号,并以子块形式写出结构的原始刚度矩阵。图12-15

13、 解：所谓子块是按单元的始末端点( 结点号)进行分块的。在形式上类似于连续梁的的单元刚度矩阵形式，但对于刚架来说，则每一子块又是阶的。分块单元刚度矩阵形式为：对本例有5个结点，故分块总刚应是的，如图b所示（即将一个结点视为“一个位移子块”）。实际上本题 以结点位移未知量考虑按后处理法，则原始刚度矩阵为阶；先处理法整体刚度矩阵为阶的。本题小结：（1）同交于一个结点的各杆件称为该结点的相关单元（例如结点1的相关单元为、，结点3的相关单元为、）；而两个结点之间有杆件直接联结者称为相关结点（例如1、2；3、4和3、5）。（2）总刚的主子块（对角线上的子块）是由结点的各相关单元的主子块叠加求得，即，如、所示。（3）总刚的副子块（非主角线上的子块），当、为相关结点时即为联结它们的单元的相应副子块，即，如、等；当、为非相关结点时即为零子块，如、等。【例12-16】试 用 直 接 刚 度 法 求 图 示 结 构 的 刚 度 矩 阵 。各 杆 长 度