结构力学之矩阵位移法.docVIP

第十二章 矩阵位移法 【例12-1】 图 a所示 连 续 梁 ，EI=常数，只 考 虑 杆 件 的 弯 曲 变 形 。分别用位移法和矩阵位移法计算。 图12-1 解：（1）位移法解 基本未知量和基本结构的确定 用位移法解的基本结构如图c所示。这里我们将结点1处的转角也作为基本未知数，这样本题仅一种基本单元，即两端固定梁。 位移法基本方程的建立 将上式写成矩阵形式 系数项和自由项 计算（须绘出单位弯矩图和荷载弯矩图） 由图d，结点力矩平衡条件，得 ，， 由图e，结点力矩平衡条件，得 ，， 由图f，结点力矩平衡条件，得 ，， 由图g，结点力矩平衡条件，得 ，， 将系数项和自由项代入位移法基本方程，得 解方程，得 由叠加法绘弯矩图，如图h所示。 （2）矩阵位移法解 对单元和结点编号（图a） 本题只考虑弯曲变形的影响，故连续梁每个结点只有一个角位移未知数。若用后处理法原始结构刚度阵为阶；用先处理法结构刚度阵为阶（已知角位移）。下面采用先处理法来说明矩阵位移法计算过程。 单元标准形式为（图b） 求局部坐标系下的单元刚度矩阵 求整体坐标下的单元刚度矩阵，因连续梁的局部坐标和整体坐标是一致的，所以有，得（注：本题用先处理法换码） ， ， 按“对号入座”规则集成总刚，得 形成荷载列阵 计算单元固端列阵 ，， （2）将单元固端列阵反号，并按“对号入座”规则送入荷载列阵 （本题结点荷载为零） == 将结构刚度矩阵及荷载列阵代入矩阵位移法方程，得 解方程，得 计算杆端弯矩 = = = 得各单元杆端弯矩后，再叠加上一相应简支弯矩图即得各单元弯矩图。将各单元弯矩图组合在一起，得整个结构的弯矩图（图h）。 小结：通过本题的计算可看到： （1）基本未知量和基本结构。位移法与矩阵位移法二者都是以结点位移为基本未知量，以单根杆件（单元）为计算对象。位移法为方便计算，有三类杆件；而矩阵位移法只有一类杆件，即两端固定等截面梁。 （2）刚度矩阵与荷载列阵的形成。位移法是用单位弯矩图和荷载弯矩图并由结点的平衡条件计算系数项和自由项的，而后形成刚度矩阵与荷载列阵的；而矩阵位移法是以单元杆端刚度元素、单元杆端荷载元素，按“对号入座”规则形成刚度矩阵与荷载列阵的。 矩阵位移法基本方程的建立，归结为两个问题：一是根据结构的几何和弹性性质建立整体刚度矩阵，二是根据受载情况形成整体荷载列阵。 （3）有（1）、（2）可知，二者的关系是：“原理同源，作法有别”。因此矩阵位移法不是一个新方法，它是新的计算工具（电子计算机）与传统力学原理（位移法）相结合的产物。 【例12-2】试求图a所示结构原始刚度矩阵中的子块 ，已知单元 ①的整体坐标的单元刚度矩阵如图c所示。 图12-2 解：本题每个结点有两个基本位知量（竖向线位移和角位移），如图b所示。单元刚度矩阵为阶（图c）。由图d所示子块形式，的元素应为单元①的j端元素（图c右下角子块）与单元②i端元素（图c左上角子块乘以2）之和，即 【例12-3】只计弯曲变形时， 。(设 EI = 1) 图12-3 解：由图d及先处理法结点位移编号图c写出各单元刚度矩阵，并按“对号入座”规则集成整体刚度矩阵。 , , 【例12-4】用先处理法写出图a所示结构刚度矩阵,E=常数。不计轴向变形影响。 图12-4 解：本题虽然是刚架，但不计轴向变形影响，即每一个结点只有一个角位移未知量。根据图b所示结点位移编号，则整体刚度矩阵为阶。由于每个单元杆端只有角位移未知量，故单元刚度矩阵为阶的连续梁单刚形式。 =，=，=，= 【例12-5】图示连续梁 ，不计轴向变形 ，EI =常数 ，已知结点位移 。试求单元②的杆端力列阵 。 图12-5 解：根据图a的约束条件和图b的结点位移编号，已知给出的结点位移是： 有，，。单元②的杆端力列阵为 【例12-6】用矩阵位移法求图a所示桁架各杆内力。单元①、②的截面面积为A，单元③的截面面积为2A，各杆E相同。 图12-6 解：桁架每个结点两个线位移未知量（图b）。 局部坐标系下的单元刚度矩阵为阶，即 =，= 整体坐标系下的单元刚度矩阵为 由图b可知，单元① ，，。单元② ，，。单元③ ，，。 =，= =。 整体刚度矩阵及荷载列阵 =，= 矩阵位移法方程 = 解