下载文档原格式
第十二章 矩阵位移法【例12-1】 图 a 所示 连 续 梁 ，EI=常数，只 考 虑 杆 件 的 弯 曲 变 形 。
分别用位移法和矩阵位移法计算。
图12-1解：（1）位移法解•基本未知量和基本结构的确定用位移法解的基本结构如图c 所示。
这里我们将结点1处的转角也作为基本未知数，这样本题仅一种基本单元，即两端固定梁。
•位移法基本方程的建立⎪⎭⎪⎬⎫=+θ+θ+θ=+θ+θ+θ=+θ+θ+θ000333323213123232221211313212111P P P R K K K R K K K R K K K 将上式写成矩阵形式⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡000321321333231232221131211P P P R R R K K K K K K K K K•系数项和自由项 计算（须绘出单位弯矩图和荷载弯矩图）由图d ，结点力矩平衡条件∑=0M ，得 EI K 411=，l EI K 221=，031=K由图e ，结点力矩平衡条件∑=0M ，得l EI K 212=，l EI l EI l EI K 84422=+=，l EI K 232=由图f ，结点力矩平衡条件∑=0M ，得 013=K ，l EI K 223=，l EI EI EI K 84433=+=由图g ，结点力矩平衡条件∑=0M ，得81Pl R p -=，2Pl R P -=，03=P R将系数项和自由项代入位移法基本方程，得⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧--+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡0000118820282024321Pl l EI •解方程，得⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ14114162321EI Pl •由叠加法绘弯矩图，如图h 所示。
（2）矩阵位移法解•对单元和结点编号（图a ） 本题只考虑弯曲变形的影响，故连续梁每个结点只有一个角位移未知数。
《结构力学》第十章矩阵位移法矩阵位移法是结构力学中的一种重要分析方法，通过将结构的受力分析转化为矩阵运算，可以有效地求解复杂结构的位移和应力分布。
本文将分为四个部分来介绍矩阵位移法的基本原理和应用。
第一部分将介绍矩阵位移法的基本原理。
矩阵位移法基于结构的受力平衡方程和变形条件，建立了适用于不同类型结构的一般形式的位移函数。
通过对这些位移函数进行适当组合，可以得到一个较为简化的位移矩阵方程。
这个方程可以通过矩阵运算求解，从而得到结构的位移和应力分布。
第二部分将介绍矩阵位移法的应用。
矩阵位移法可以用于求解各种类型的结构，包括梁、柱、框架等。
具体应用时，首先需要确定结构的边界条件和受力情况，然后根据结构的几何形状和材料性质，建立相应的位移函数。
之后，将位移函数按照一定的规则组合起来，建立一个位移矩阵方程。
通过解这个方程，可以得到结构的位移和应力分布。
第三部分将介绍矩阵位移法的优点。
相比于传统的力方法，矩阵位移法具有计算简单、准确性高、适用范围广等优点。
这是因为矩阵位移法可以通过矩阵运算将结构的受力分析转化为代数运算，减少了繁琐的计算过程，并且可以应用于各种不规则结构。
第四部分将介绍矩阵位移法的局限性。
矩阵位移法虽然具有很多优点，但也有一些限制。
首先，矩阵位移法对结构的刚度矩阵的求取较为复杂，需要通过精确和谐振数法等途径进行求解。
其次，矩阵位移法不能用于解决非线性和动力问题。
总结起来，矩阵位移法是一种重要的结构力学分析方法，通过将结构的受力分析转化为矩阵运算，可以有效地求解复杂结构的位移和应力分布。
它具有计算简单、准确性高、适用范围广等优点，但也有一些局限性。
因此，在实际应用中需要根据具体情况选择合适的方法。
同时，矩阵位移法的进一步研究和发展也是一个非常重要的方向。
第十二章 矩阵位移法【例12-1】 图 a 所示 连 续 梁 ，EI=常数，只 考 虑 杆 件 的 弯 曲 变 形 。
分别用位移法和矩阵位移法计算。
图12-1解：（1）位移法解∙基本未知量和基本结构的确定用位移法解的基本结构如图c 所示。
这里我们将结点1处的转角也作为基本未知数，这样本题仅一种基本单元，即两端固定梁。
∙位移法基本方程的建立⎪⎭⎪⎬⎫=+θ+θ+θ=+θ+θ+θ=+θ+θ+θ000333323213123232221211313212111P P P R K K K R K K K R K K K 将上式写成矩阵形式⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡000321321333231232221131211P P P R R R K K K K K K K K K∙系数项和自由项 计算（须绘出单位弯矩图和荷载弯矩图）由图d ，结点力矩平衡条件∑=0M ，得 l EI K 411=，l EI K 221=，031=K由图e ，结点力矩平衡条件∑=0M ，得l EI K 212=，l EI l EI EI K 84422=+=，l EI K 232=由图f ，结点力矩平衡条件∑=0M ，得 013=K ，l EI K 223=，l EI l EI l EI K 84433=+=由图g ，结点力矩平衡条件∑=0M ，得1Pl R p -=，2Pl R P -=，03=P R将系数项和自由项代入位移法基本方程，得⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧--+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡0000118820282024321Pl l EI ∙解方程，得⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ14114162321EI Pl ∙由叠加法绘弯矩图，如图h 所示。
（2）矩阵位移法解∙对单元和结点编号（图a ） 本题只考虑弯曲变形的影响，故连续梁每个结点只有一个角位移未知数。
若用后处理法原始结构刚度阵为44⨯阶；用先处理法结构刚度阵为33⨯阶（已知角位移04=θ）。
下面采用先处理法来说明矩阵位移法计算过程。
单元标准形式为（图b ）)(e k ⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=)()()()()(4224e jj e jie ije ii e k k k k l EI l EI l EI l EI∙求局部坐标系下的单元刚度矩阵)(e k∙求整体坐标下的单元刚度矩阵T k T k e T e )()(=，因连续梁的局部坐标和整体坐标是一致的，所以有)()(e e k k=，得（注：本题用先处理法换码）)1(k 214224)1(⎥⎦⎤⎢⎣⎡=lEI， )2(k 324224)2(⎥⎦⎤⎢⎣⎡=l EI ，)3(k 034224)3(⎥⎦⎤⎢⎣⎡=l EI ∙按“对号入座”规则集成总刚，得=K 321820282024⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡l EI∙形成荷载列阵P(1) 计算单元固端列阵=)1(F F 218181⎭⎬⎫⎩⎨⎧-Pl ，=)2(F F 324141⎭⎬⎫⎩⎨⎧-Pl ，=)3(FF 034141⎭⎬⎫⎩⎨⎧-Pl （2）将单元固端列阵反号，并按“对号入座”规则送入荷载列阵P （本题结点荷载为零）P =E D P P +=321081814141418181000⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+-+-+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧Pl Pl∙将结构刚度矩阵及荷载列阵代入矩阵位移法方程P K =∆，得⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡011820282024321Pl l EI∙解方程，得⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ14114162321EI Pl ∙计算杆端弯矩)()()()()()()()()()(e e e F e e e F e e e F e T k F k F k F F ∆+=∆+=δ+=)1(F=⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧+⎭⎬⎫⎩⎨⎧-=⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡+⎭⎬⎫⎩⎨⎧-4502083852416525241641141642241812Pl Pl Pl EI Pl l EI Pl)2(F=⎭⎬⎫⎩⎨⎧-=⎭⎬⎫⎩⎨⎧+⎭⎬⎫⎩⎨⎧-=⎭⎬⎫⎩⎨⎧-⎥⎦⎤⎢⎣⎡+⎭⎬⎫⎩⎨⎧-544520841441610410441614416422441412Pl Pl Pl EI Pl l EI Pl )3(F=⎭⎬⎫⎩⎨⎧-=⎭⎬⎫⎩⎨⎧--+⎭⎬⎫⎩⎨⎧-=⎭⎬⎫⎩⎨⎧-⎥⎦⎤⎢⎣⎡+⎭⎬⎫⎩⎨⎧-51542082441610410441601416422441412Pl Pl Pl EI Pl l EI Pl 得各单元杆端弯矩后，再叠加上一相应简支弯矩图即得各单元弯矩图。
将各单元弯矩图组合在一起，得整个结构的弯矩图（图h ）。
小结：通过本题的计算可看到： （1）基本未知量和基本结构。
位移法与矩阵位移法二者都是以结点位移为基本未知量，以单根杆件（单元）为计算对象。
位移法为方便计算，有三类杆件；而矩阵位移法只有一类杆件，即两端固定等截面梁。
（2）刚度矩阵与荷载列阵的形成。
位移法是用单位弯矩图和荷载弯矩图并由结点的平衡条件计算系数项和自由项的，而后形成刚度矩阵与荷载列阵的；而矩阵位移法是以单元杆端刚度元素、单元杆端荷载元素，按“对号入座”规则形成刚度矩阵与荷载列阵的。
矩阵位移法基本方程的建立，归结为两个问题：一是根据结构的几何和弹性性质建立整体刚度矩阵K ，二是根据受载情况形成整体荷载列阵P 。
（3）有（1）、（2）可知，二者的关系是：“原理同源，作法有别”。
因此矩阵位移法不是一个新方法，它是新的计算工具（电子计算机）与传统力学原理（位移法）相结合的产物。
【例12-2】试求图a 所示结构原始刚度矩阵中的子块 22K ，已知单元 ①的整体坐标的单元刚度矩阵如图c 所示。
图12-2解：本题每个结点有两个基本位知量（竖向线位移和角位移），如图b 所示。
单元刚度矩阵为44⨯阶（图c ）。
由图d 所示子块形式，22K 的元素应为单元①的j 端元素（图c 右下角子块）与单元②i 端元素（图c 左上角子块乘以2）之和，即)2(22)1(22)2()1(22K K K K K ii jj +=+=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡--=60000360036002164000072007200144200003600360072【例12-3】只计弯曲变形时，用先处理法写出结构刚度矩阵K 。
(设 EI = 1)图12-3解：由图d 及先处理法结点位移编号图c 写出各单元刚度矩阵，并按“对号入座”规则集成整体刚度矩阵。
)1(k 210045.125.15.175.05.175.025.145.15.175.05.175.0⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=, )2(k 30210.25.10.15.15.15.15.15.10.15.10.25.15.15.15.15.1⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------= )3(k 4030667.2333.1333.1333.1333.1889.0333.1889.0333.1333.1667.2333.1333.1889.0333.1889.0⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=,K 4321667.2333.100333.1667.415.1016005.1025.2⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=【例12-4】用先处理法写出图a 所示结构刚度矩阵K ,E=常数。
不计轴向变形影响。
图12-4解：本题虽然是刚架，但不计轴向变形影响，即每一个结点只有一个角位移未知量。
根据图b 所示结点位移编号，则整体刚度矩阵为33⨯阶。
由于每个单元杆端只有角位移未知量，故单元刚度矩阵为22⨯阶的连续梁单刚形式。
)1(k =214224⎥⎦⎤⎢⎣⎡l EI ，)2(k =208448⎥⎦⎤⎢⎣⎡l EI ，)3(k =328448⎥⎦⎤⎢⎣⎡l EI ，K =3218404202024⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡ 【例12-5】图示连续梁 ，不计轴向变形 ，EI =常数 ，已知结点位移∆T43812⎥⎦⎤⎢⎣⎡--=EI ql EI ql 。
试求单元②的杆端力列阵 。
图12-5解：根据图a 的约束条件和图b 的结点位移编号，已知给出的结点位移是：{}⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧∆∆=∆3221v θ 有03211=θ==θ=v v ，EI ql 1232-=θ，EI ql v 843-=。
单元②的杆端力列阵为)2(F ⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧-=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧--⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡---=22432322323127 125 0 8120 4 6 12 2 6 4 6 12 6 21ql ql ql ql EI ql EI ql l EI l EI l EI l EI l EI l EI l EI l EI l EI l EI 称对【例12-6】用矩阵位移法求图a 所示桁架各杆内力。
单元①、②的截面面积为A ，单元③的截面面积为2A ，各杆E 相同。
图12-6解：桁架每个结点两个线位移未知量（图b ）。
∙局部坐标系下的单元刚度矩阵为44⨯阶，即)(e k =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---0000010********l EA ，T =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡αα-αααα-ααcos sin 00sin cos 0000cos sin 00sin cos∙整体坐标系下的单元刚度矩阵为T k T k e T e )()(=由图b 可知，单元① 030=α，23sin =α，21cos =α。
单元② 045=α，22sin =α，22cos =α。
单元③ 090=α，1sin =α，0cos =α。
)1(k =210013133333131333338⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--------l EA ，)2(k =21002222222222222222222212222222228⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--------l EA )3(k =210160160000016016000008⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--l EA 。
∙整体刚度矩阵及荷载列阵K =⎥⎦⎤⎢⎣⎡47855.257006.057006.072855.0l EA ，P =⎭⎬⎫⎩⎨⎧0P∙矩阵位移法方程⎥⎦⎤⎢⎣⎡47855.257006.057006.072855.0l EA ⎭⎬⎫⎩⎨⎧11v u =⎭⎬⎫⎩⎨⎧0P ∙解方程，得⎭⎬⎫⎩⎨⎧-=⎭⎬⎫⎩⎨⎧38497.067381.111EA Pl v u∙计算各杆轴力)1(F =)1()1(∆Tk =⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--23100212300002321002123)1(k ⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-38497.067381.100EA Pl =⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-06285.006285.0P （拉） )2(F =)2()2(∆Tk =⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--22220022220000222202222)2(k ⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-38497.067381.100EA Pl =⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-06442.006442.0P （拉） )3(F =)3()3(∆Tk =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--0100100000010010)3(k ⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-38497.067381.100EA Pl =⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-07699.007699.0P （压）【例12-7】已知图示桁架的自由结点位移列阵∆ ，求杆12在局部坐标系中的杆端 力 。
