结构力学之矩阵位移法

第十二章 矩阵位移法
【例12-1】 图 a 所示 连 续 梁 ，EI=常数，只 考 虑 杆 件 的 弯 曲 变 形 。

分别用位移法和矩阵位移法计算。

图12-1
解：（1）位移法解
∙基本未知量和基本结构的确定
用位移法解的基本结构如图c 所示。

这里我们将结点1处的转角也作为基本未知数，这样本题仅一种基本单元，即两端固定梁。

∙位移法基本方程的建立
⎪⎭
⎪
⎬⎫
=+θ+θ+θ=+θ+θ+θ=+θ+θ+θ000333323213123232221211313212111P P P R K K K R K K K R K K K 将上式写成矩阵形式
⎪⎭
⎪⎬⎫
⎪⎩⎪⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧θθθ⎥⎥
⎥⎦⎤⎢⎢⎢⎣⎡0003213213332
31
2322
21131211P 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 ji
e ij
e 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(⎥⎦⎤⎢⎣⎡=l
EI
， )
2(k 324224)
2(⎥⎦
⎤⎢⎣⎡=l EI ，)
3(k 0
3
4224)
3(⎥⎦
⎤⎢⎣⎡=l EI ∙按“对号入座”规则集成总刚，得
=K 3
21820282024⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡l EI
∙形成荷载列阵P
(1) 计算单元固端列阵
=)1(F F 218181⎭⎬⎫⎩⎨⎧-Pl ，=)2(F F 324141⎭⎬⎫⎩⎨⎧-Pl ，=)
3(F
F 0
34141⎭⎬⎫⎩⎨⎧-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 21
00
45.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 40
30667.2333.1333.1333.1333.1889.0333.1889.0333.1333.1667.2333.1333.1889.0333.1889.0⎥
⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡------=,K 4
3
2
1667.2333.100333.1667.415.10160
05
.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 =3
21
8404202024⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡ 【例12-5】图示连续梁 ，不计轴向变形 ，EI =常数 ，已知结点位移
∆T
43
812⎥⎦
⎤⎢⎣⎡-
-
=EI ql EI ql 。

试求单元②的杆端力列阵 。

图12-5
解：根据图a 的约束条件和图b 的结点位移编号，已知给出的结点位移是：
{}⎭
⎬⎫
⎩⎨⎧=⎭⎬⎫⎩⎨
⎧∆∆=∆3221v θ 有03211=θ==θ=v v ，EI ql 1232-=θ，EI ql v 843-=。

单元②的杆端力列阵为
)
2(F ⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧-=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩
⎪⎪⎪⎨⎧--⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤
⎢⎢⎢
⎢⎢⎢⎢⎢⎢⎣
⎡---=224
32322323127 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 =
⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢
⎣⎡---00
00
010********l EA ，T =⎥⎥⎥⎥
⎦
⎤⎢⎢⎢⎢⎣⎡αα-αααα-α
αcos sin 00sin cos 0000cos sin 00sin cos
∙整体坐标系下的单元刚度矩阵为
T k T k e T e )()(=
由图b 可知，单元① 0
30=α，23sin =α，21cos =α。

单元② 0
45=α，
22sin =α，22cos =α。

单元③ 090=α，1sin =α，0cos =α。

)1(k =210013133333131333338⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--------l EA ，)2(k =2
1002222222222222222222212222222228⎥⎥
⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡--------l EA )
3(k =
2
1
01601600000160160000
08⎥
⎥⎥⎥⎦⎤⎢
⎢⎢⎢
⎣⎡--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 =⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢
⎣⎡--23100
21230000232
100
2123)1(k ⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-38497.067381.100EA Pl =⎪⎪⎭⎪
⎪⎬⎫
⎪⎪⎩⎪⎪⎨⎧-06285.006285.0P （拉） )2(F =)2()2(∆Tk =⎥⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎢
⎣
⎡--
22
2
2002222000022220
222
2)2(k ⎪⎪⎭⎪⎪⎬⎫
⎪⎪⎩⎪⎪⎨⎧-38497.067381.100EA Pl =⎪⎪
⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-06442.006442.0P （拉） )3(F =)3()3(∆Tk =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢
⎣⎡--01
00
10000001
0010
)3(k ⎪⎪⎭
⎪⎪⎬⎫
⎪⎪⎩⎪⎪⎨⎧-38497.067381.100EA Pl =⎪⎪
⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-07699.007699.0P （压）
【例12-7】已知图示桁架的自由结点位移列阵∆ ，求杆12在局部坐标系中的杆端 力 。

设2kN/cm 3000=E ，杆12 的横截面积2
cm 18=A 。

∆m 10834.341803.613 6
-⨯⎭
⎬⎫⎩⎨
⎧-=
图12-7
解：0
16.53=α，8.0sin =α，6.0cos =α。

cm kN cm
cm cm kN l EA /906001830002
2=⨯=
)
()()
(e e e T k F
∆
==kN
⎪⎪⎭
⎪⎪⎬⎫
⎪⎪⎩⎪⎪⎨⎧-=⨯⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧-⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---032.85032.851000834.341803.6136.08.0008.06.000006.08.0008.06.000
00
0900900000090090
4
【例12-8】 用位移法和矩阵位移法计算图a 所示结构。

各杆材料及截面均相同，
28/100.2m kN E ⨯=，451032m I -⨯=，22101m A -⨯=。

要求：（1）不考虑轴向变形影响的位移法解。

（2）考虑轴向变形影响的位移法解。

（3）用矩阵位移法(采用先处理法)解。

图12-8
解：（1）不考虑轴向变形影响的位移法求解
不考虑轴向变形影响下，仅有结点1处的角位移未知量1Z 。

位移法的基本方程为 01111=+P R Z K 系数和自由项由图b 、c 得
m kN l EI K .128000811==，m kN ql R P .3401221-=-=
将系数和自由项由代入位移法的基本方程，并解得 4110042.1-⨯=Z 弧度。

由叠加法作弯矩图，即11Z M M M P +=。

整个结构的弯矩图如图d 所示。

（2）考虑轴向变形影响的位移法求解
基本结构如图e 所示。

位移法的基本方程为
⎪⎭
⎪
⎬⎫
=+++=+++=+++000333323213123232221211313212111P P P R Z K Z K Z K R Z K Z K Z K R Z K Z K Z K 系数和自由项计算
由图f ：m kN l EA l EI K 5
5
4
3
111012.510.5102.112⨯=⨯+⨯=+=
021=K ，kN EI K 4231104.26⨯=-=
由图g ：m kN l EA l EI K K 5
3
11221012.512⨯=+==，
kN l EI K 4232104.26⨯=-=
由图h ：m kN EI K .128000833==
由图c ：01=P R ，kN ql R P 2022==，m kN ql R P .401223-=-= 将系数和自由项由代入位移法的基本方程，并解得
m Z 6110621.4-⨯=，m Z 5210444.3-⨯-=，5310858.9-⨯=Z 弧度
考虑轴向变形影响的结构弯矩图如图i 所示（剪力图和轴力图未画出）。

（3）用矩阵位移法(采用先处理法)解
用矩阵位移法求解时，单元和结点编号如图j 所示。

采用先处理法时其整体刚度矩阵为33⨯阶。

两单元对应的整体编码如下图所示。

按“对号入座”规则集成结构刚度矩阵
K =⎥⎥
⎥⎥⎥⎥⎦⎤
⎢⎢⎢
⎢
⎢⎢⎣⎡+-
--+-
+l EI l EI l EI l EI l EI l EA
l EI l EI l
EA
l EI 4466612060122222
注：（1）单元①局部坐标与整体坐标一致，所以有)1()
1(k k
=。

（2）单元②局部坐标与整体坐标的夹角0
90=α，须进行坐标变换，即
T k T k e T e )()(=。

T k T e T )(运算的结果是将)1(k 中相关元素作行列交换。

另外当局部坐
标与整体坐标的夹角0
90=α时，我们也可直接在整体坐标系下进行对换，如图k 所示。

按先x 后y 再转角的次序，则可直接在局部坐标的单元上标注相应的整体编码，本题就是采用这一方法。

注意到坐标进行了y x ,轴交换，αsin 变号，故副系数须反号。

见本题中单元②中送入结构刚度矩阵的元宵13K 和31K 。

荷载列阵的集成。

方法一是按∑=-
=n
i i F
E F
P 1
及)()
(e F T e F
F T F =进行。

另一作法是，由
R P -=,于是有
P =321340200⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧-
将结构刚度矩阵K 和荷载列阵P 基本方程，得与前位移法解得的相同结果，即
∆=θ⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧⨯⨯-⨯---v u
55610858.910444.310621.4
同样得结构弯矩图如图i 所示（剪力图和轴力图未画出）。

【例12-9】试 求 用 矩 阵 位 移 法 求解 图 a 所 示 结 构 时 ，结 点 2 的 综 合 结 点 荷 载 列 阵 2P 。

解：刚架每个结点有三个基本未知量（θ,,v u ），同时也有三个方向结点荷载项。

图12-9
（1）结点2的直接结点荷载：D P 23218002
⎪⎪⎭
⎪
⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧=ql
（2）结点2的等效结点荷载涉及到单元①、②及③的2端的固端力（见图c 、d 、e ）。

按式∑=-
=n
i i
F
E F
P 1
（)()
(e F T e F
F T F =）应首先应计算局部坐标系下的固端反力)(e F F ，而后
进行坐标变换得整体坐标系下单元固端反力)
(e F F ，再“按对号如座”规则反其符号集成。

这里我们直接根据图c 、d 、e 求出整体坐标系下的单元固端反力)(e F F 。

由图b 及d 、e 、c 得
)1(F F =3210001202120222⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧---ql ql ql ql ，)2(F F =98732182082022⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧-ql ql ql ql ，)3(F F =6
54321
1202120222
⎪⎪⎪⎭⎪⎪
⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧---ql ql ql ql
E
P 2321
82120282012022)
3(2)
2(2)
1(2⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧--=⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧-+⎪⎭⎪⎬⎫
⎪⎩⎪⎨⎧--+⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧=ql ql ql ql ql ql ql ql ql 结 点 2 的 综 合 结 点 荷 载 列 阵为
E
D P P P 222+=3210
28280022⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧-=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧--+⎪⎭⎪⎬⎫⎪⎩
⎪⎨⎧=ql
ql ql ql ql ql 。

【例12-10】 试 用 先 处 理 法 写 出 图 a 所 示 结 构 刚 度 矩 阵 K 。

各 杆 杆 长 均 为 l ，EI =
常 数 ，自 由 结 点 位 移 分 量 的 编 号 如 图 示 。

图12-10
解：单元①与整体坐标一致。

而单元②、③按图b 所示整体坐标系下来进行换码（注意到坐标进行了y x ,轴交换，αsin 变号，故副系数须反号），而后按下图“对号入座”规则集成总刚。

K =⎥⎥⎥⎥⎥⎥⎥⎦
⎤⎢⎢
⎢
⎢⎢⎢⎢⎣
⎡++--++++)3()2()1()1(2
)
1(2
)
3()2()
1(3
)
3(3)2(3)1(444606120
1212l EI
l EI l EI l
EI l
EI l EA l EA l EI
l EI l EI l EA
【例12-11】 用 先 处 理 法 求 图 a 所 示 刚 架 的 结 构 刚 度 矩 阵 K ，略去轴向 变 形 影响。

图12-11
解：由图b 的位移编号可知，横梁各结点仅有一个x 向的水平位移，其变形如图c 的所示（这就是“手算”），按“对号入座”规则集成总刚（这就是“机算”）
K ==
++)3()2()1(k k k 333336121212l
EI
l EI l EI l EI =++
用经典位移法解时，其系数3
1136l
EI
K =。

【例12-12】按先处理法计算图a 所示结构的刚度矩阵K 。

各杆长度为 l ，EA 、EI 均为相同 。

图12-12
解：单元、结点及位移编号入图b 所示。

作为理解画出了结点位移的变形图，如图c 、d 及e 所示（这就是“手算”）。

按下图“对号入座”规则集成总刚（这就是“机算”）。

K =⎥⎥
⎥⎥⎥
⎥⎥⎦⎤⎢⎢
⎢
⎢⎢
⎢⎢⎣
⎡⨯++++++l EI l
EI l EI l
EA l
EA
l EI l EI l EA l EA 440
0012120
1212)
2(3
)
1(3
)
4()
3()
4(3)3(3)2()1(
【例12-13】图示刚架只考虑弯曲变形 ，按先处理法求在荷载和支座位移共同作 用下的结点荷载列阵P 。

已知各杆2
2
m kN 102⋅⨯=EI 。

图12-13
解：图b 为结点、单元编号，单元①固端反力如图c 所示 ，是由支座位移产生的。

D P =⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧0005kN ，
E P =⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧0.205.1m kN kN ，P =D P +E P =⎪⎪⎭
⎪
⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧0.205.6m kN kN 【例12-14】 图 示 刚 架 各 杆 2
m kN 64⋅=EI ，结点6有支 座 的 水 平 位 移 m 01.05=∆，竖 向 位 移 m 01.06-=∆，忽 略 轴 向 变 形 ，已 求 得 结 点 位 移 为：
{}[]T 140.000547
- 0.001719- 0.000547- 005208.0=∆⨯ 。

求 单 元 ③ 的 杆 端 内 力 。

图12-14
解：本题有两各特点：
（1） 不计轴向变形影响，单元刚度为44⨯阶，如图b 所示，不需坐标变换。

（2） 结点6的支座移动只有5∆对单元③有影响，将它作为单 元 ③ 杆端位移值，则有
)3(δ[][]T
T
j j
i
i v v 000547.0005208
.0001.0---==θθ
所以 )
3()3()
3(δ=k F
CF CF
FC
FC
M Q M Q
1500.00706.0 1325.00706.0000547.0005208.00 01.0 64 24 32 24 24 12 24 1232 24 64 24 24 12 24 12 ⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧---=⎪⎪⎭⎪⎪⎬⎫⎪
⎪⎩⎪⎪⎨⎧---⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------=
【例12-15】对图示刚架的结点和单元进行编号,并以子块形式写出结构的原始刚度矩阵。

图12-15
解：所谓子块是按单元的始末端点( 结点号)j i ,进行分块的。

在形式上类似于连续梁的22⨯的单元刚度矩阵形式，但对于刚架来说，则每一子块又是33⨯阶的。

分块单元刚度矩阵形式为：
对本例有5个结点，故分块总刚应是55⨯的，如图b 所示（即将一个结点视为“一个位移子块”）。

实际上本题 以结点位移未知量考虑按后处理法，则原始刚度矩阵为1515⨯阶；先处理法整体刚度矩阵为66⨯阶的。

本题小结：
（1）同交于一个结点的各杆件称为该结点的相关单元（例如结点1的相关单元为①、②，结点3的相关单元为③、④）；而两个结点之间有杆件直接联结者称为相关结点（例如1、2；3、4和3、5）。

（2）总刚的主子块（对角线上的子块）ii K 是由结点i 的各相关单元的主子块叠加求得，即∑=
)
(e ii ii k K ，如11K 、33K 所示。

（3）总刚的副子块（非主角线上的子块）im K ，当i 、m 为相关结点时即为联结它们
的单元的相应副子块，即)
(e im im k K =，如12K 、13K 等；当i 、m 为非相关结点时即为零子
块，如14K 、15K 等。

【例12-16】试 用 直 接 刚 度 法 求 图 示 结 构 的 刚 度 矩 阵 K 。

各 杆 长 度 为 l 。

图12-16
解：本题特点：
（1） 在结点D 具有半铰的情况；
（2） 结构中有两类单元（桁架单元DC 和刚架单元AD 、DB ）。

为了统一每一个结点均为三个基本未知量。

但桁架单元的杆端是以线位移为基本未知量的，故转角方向为无效未知量，其变号为零，如图b 所示。

要注意到单元①、②、③在结点D 的x 、y 的位移是一致的，所以结点4、5的u 、v 编号均为3、4编号。

因单元③垂直且为桁架单元，故只能将轴向刚度l EA 送如主元44K 中（注：
3对应的项不送）。

按“对号入座”规则集成结构刚度矩阵K 。