习题9.1是非判断题
（1） 矩阵位移法既可计算超静定结构，又可以计算静定结构。（ ）
（2） 矩阵位移法基本未知量的数目与位移法基本未知量的数目总是相等的。（ ）
（3） 单元刚度矩阵都具有对称性和奇异性。（ ）
（4） 在矩阵位移法中，整体分析的实质是建立各结点的平衡方程。（
）
（5） 结构刚度矩阵与单元的编号方式有关。（
）
（6） 原荷载与对应的等效结点荷载使结构产生相同的内力和变形。（
）
【解】（1）正确。
（2） 错误。位移法中某些不的杆端位移不计入基本未知量。 （3） 错误。不计结点线位移的连续梁单元的单刚不具奇异性。 （4） 正确。
（5） 错误。结点位移分量统一编码会影响结构刚度矩阵，但单元或结点编码则不会。 （6） 错误。二者只产生相同的结点位移。
习题9.2填空题
（1） 矩阵位移法分析包含三个基本环节，其一是结构的,其二是 分析，其三
是 分析。
（2） 已知某单元的定位向量为[3 5 6 7 8 9]七则单元刚度系数炫应叠加到结构刚度
矩阵的元素 中去。
（3） 将非结点荷载转换为等效结点荷载，等效的原则是。
（4） 矩阵位移法中，在求解结点位移之前，主要工作是形成 矩阵和
_________________ 列阵。
（5） 用矩阵位移法求得某结构结点2的位移为4=[. V2 ft]T=[0.8 0.3 0.5]T,单元①的始、
末端结点码为3、2,单元定位向量为尸＞=[0 0 0 3 4 5]T ,设单元与x轴之间的夹角为a =买，则 尹＞ =O
（6 ）用矩阵位移法求得平面刚架某单元在单元坐标系中的杆端力为 Fe =[7.5 -48 -70.9 -7.5 48 -121.09]T ,则该单元的轴力 心=kN。
【解】（1）离散化，单元，整体； （2） 灯8；
（3） 结点位移相等；
（4） 结构刚度，综合结点荷载； （5） [0 0 0 0.3 -0.8 0.5]。 （6） -7.5o
2
离、空的值以及K⑴中元素妍、愚、姒的值。
习题9.3图
【解】各刚度系数的物理意义如习题解9.3图所示。因此，各刚度系数的值为
=EA/l ,
C =12EI/l3,嫩=0,襟=0。
=6EI/l2,恐）=-6£7/尸；
(c)右？的物理意义
J-
(d)妃〉的物理意义
(e)馅？的物理意义
(f)以？的物理意义
习题解9.3图
习题9.4根据结构刚度矩阵元素的物理意义，直接求出习题9.4图所示刚架结构刚度矩阵中的 元素如、如、心的值。各杆E、A, I相同。
3EI
4F
习题9.4图
【解】各刚度系数的物理意义如习题解9.4图所示。因此，各刚度系数的值为
7 12EI EA 7 八， 佑］=—^―+ 容，幻1=0, k32
力
1=1
kn
(a)知和心的物理意义 (b\"32的物理意义
习题解9.4图
习题9.5用简图表示习题9.5图所示刚架的单元刚度矩阵左⑴中元素史°, K⑵中元素待)的物
理意义。
习题9.5图
【解】各刚度系数的物理意义如习题解9.5图所示。
2
(a)说)的物理意义
习题解9.5图
习题9.6习题9.6图所示刚架各单元杆长为I, EA. 意义，写出单元刚度矩阵K⑴、K⑵的第3列和第5列元素。
政为常数。根据单元刚度矩阵元素的物理
习题9.6图
【解】各列刚度系数的物理意义如习题解9.6图所示。因而
双中第3列元素:
砂＞ 中第5列元素:
6EI 4EI 6EI 2EI
0 — -------------- 0 z——
I2 I I2 I
n 12EI 6EI n 12EI 6EI
Z2
K②中第3列元素：0 — 一哮0 —
I2 I i1 I
L
_ 「 EA EA ~|T 中第5列元素：0 --- 0 0 — 0
习题解9.6图
习题9.7用先处理法，对习题9.7图所示结构进行单元编号、结点编号和结点位移分量编码, 并写出各单元的定位向量。
习题9.7图
【解】离散化结果如习题解9.7图所示。
因而，各单元定位向量为
舟=[1 0 0 2 3 4]T, 人⑵二任3 4 2。)=[5 6 7 0 0 9]T 人”)=[5 6 8
5 6 7]T, 1(100) 0 0 0]To 2(2,3.4)
X
1 V
②④
3(5,6,7) 0 ■
4(5,6,8) 5(0,0,9) <
I
|6(0,0,0) 1
习题解9.7图
本题可有多种离散化方法，因此上述答案不是唯一的正确答案。 习题9.8用先处理法形成习题9.8图所示结构的综合结点荷载列阵。
5
4m
【解】离散化如习题解9.8图所示。
1(1,0,2)
_顷_|
习题9.8图
3(6,7,8)
心
方①2&4,5)②③4息产
L
5(0,0,0)
④
习题解9.8图
非结点荷载引起的单元固端力为
邦2)=[0 -12 -8 0 -12 8]T,
各单元的等效结点荷载列阵为
=[0 -9 -4.5 0 -9 4.5]T
人⑵-> 3
=-尸1尺2)=一尺2)=[0
4 5 6 12 8 0
7 8 12 - 9 0 9 -4.5]
/J 6 7 8 0 0
理)
=_尸拜3)=_拜3)=[0
=
9 4.5
集成为结构的等效结点荷载列阵
化=[0 0 0 12 8 0 21 -3.5 9]「
直接结点荷载列阵为
f；=[0 -5 0 4 0 0 0 0 0]T
综合结点荷载列阵为
P = Pi+PE =[o -5 0 16 8 0 21 -3.5
习题9.9用先处理法求习题9.9图所示连续梁的结构刚度矩阵和结构的综合结点荷载列阵。已 知：EZ=2.4xl04kNm2
8kN 8kN
5kN-m
.2
6kN/m
EI 3' 习题9.9图
EI
'4
【解】离散化如习题解9.9图所示。本题无需坐标转换。
妲？匚 2,2) H 3(3)「言(4)
® X _®X ®
习题解9.9图
先求结构刚度矩阵。各单元的单刚为
1 2
2 3 ■ 1 1/2- 1 _1/2
1 2,
_
3 4 ~2/3 1/3- 2 =EI
_1/3 2/3_ 3, 0 1/3 22/15 2/5
■4/5 2/5-
K& = EI _2/5 4/5_
集成即可得到结构刚度矩阵
K = EI
1 1/2
5/3
称
对
0 - 「2.4 1.2 0.0 0 4.0 0.8
4 =10
对
0.0 -
0.0 3.52 0.96 1.92
4/5
称
再求综合结点荷载列阵。非结点荷载作用单元的等效结点荷载列阵为
2 3 3 4
=[10.67 -10.67]T, P^} =[12.5 -12.5]T
集成为结构的等效结点荷载列阵
《=[0 10.67 1.83 12.5]T
综合结点荷载列阵为
尸=当+弓=[5 0 0 0] +《=[5 10.67 1.83 12.5]T
习题9.10用先处理法求习题9.10图所示结构刚度矩阵。忽略杆件的轴向变形。各杆
£/=5xl05kNm2o
习题 9.10图
【解】离散化如习题解9.10图所示。因为不计各杆轴向变形，所以本题只涉及转角位移未知量,
无需坐标转换。
各单元的单刚为
1 2
-2
4/5 2/5 1
3
■4/5 2/5 2 _2/5 4/5_ 3,
-
2
0
-
3 0
~ 1 1/2- Kw= EI _1/2 1
2, 2/5 一4/5_
Km = EI
■ 1 1/2 2 = EI _1/2 1 °,
集成即可得到结构刚度矩阵
4/5 2/5
K = EI 2/5 13/5
-
0 - 0一2
5 4 13 一2/5 =10 2 2 9/5
2(2)
0
1(1)
2/5 0 2
3(3)
9
①
1
②
③
4(0)
X
④
bo
)
习题9.11用先处理法建立习题9.11图所示结构的矩阵位移法方程。已知：各杆EA=4xlO*N,
习题解9.10图
£/=5xl04kN-m2。
【解】1)离散化如习题解9.11(a)图所示。
1(0,1,0) 2(2,3,4)
①
3(5,0,6)
4(0.0.0)
(a)离散化
习题解9.11图
2) 计算结构刚度矩阵
各单元单刚分别为：单元①
0
13.33 0 0 -13.33 0
1 0 2.222 3.333
0 0 3.333 6.667
2 -13.33
K⑴=欧)=］。4
0 0 -2.222 -3.333 3.333
3.333 :
0
4 0 0
0 -2.222 3.333 1
3.333 0 ° -3.333
2 13.33 0 0
0 2.222 -3.333 3 0 -3.333 6.667 4
3
0
单元②
2 ■ 10.00 0
=Km =104
0 -10.00 0 0
3 0 0.9375 1.875
4 0 1.875 5.000
5 -10.00 0
0
6
0 一 2
-0.9375 1.875 3 -1.875 2.500
4
0 0 5 0.9375 -1.875 0 -1.875 5.000
6
0 0
0 0 -0.9375 -1.875 1.875
2.500
10.00 0 0
单元③
_2 0.9375 0 4 0 -1.875 ； -0.9375
10.00 0 0 0 0 -10.00 0
5.000 1.875 0 2.500
1.875 0.9375 0 1.875
3 0 0 0 -10.00 0 0 10.00 0
0
-1.875- 2 0
3
2.500 4 1.875 0 0 0 5.000
0
Km =rTK(3)r = io4
-1.875 -0.9375 0 -1.875
集成为总刚
_
2.222 0 0 24.27 0
-2.222 3.333 0 0 -1.875 -10.00 13.16
-1.458 16.67
0
0 0 10.00 0
0 0 1.875 2.500 0 5.000
K=104
-2.222 3.333 0 0
-1.875 -1.458 -10.00 0 0
1.875 2.500
2）计算综合结点荷载列阵
除可以按照习题9.8的方法计算外，还可以直接根据其物理意义形成综合结点荷载列阵。具体 做法如下： 将原结构上各结点位移未知量利用附加约束住后，施以原结构所受荷载。这一过程可理解 成在矩阵位移
=因琅生％4 鼻5 4j =[-8 0 -18 -12
0 12「
因此，综合结点荷载列阵为
P =-外=[8 0 18 12 0 -12]
T
3）列出结构刚度方程瓦仁F
一 2.222
0 24.27
-2.222 3.333 0
-1.875 -10.00
13.16 -1.458
16.67
0 0 0 10.00
0 - 0 1.875 2.500 0
_
u2
8 - 0
104 对
—
18 12 0 -12
称
“3
5.000 _ A.
9.12图所示刚架的结构刚度矩阵。已知：EA=3.2xlO5kN ,
法（先处理法）的基本结构上，作用外荷载，形成如习题解9.11（b）图所示的矩阵位移 法基本体系。由此，可得各附加约束上的反力为
玫=4.8 x \"kN 5。
习题9.12用先处理法计算习题
习题9.12图
【解】离散化如习题解9.12图所示。各单元单刚分别为
2(2,3,4)
①
N8 x ©
3(0,0,0)
I
习题解9.12图
\\y
单元①
2 ■ 6.400
3 0 0.4608 1.152
4
0 0 ； -6.400
0
1 0 -0.4608 -1.152
0 0
2 0 0
1.152 1.152 3 1.920 4
K
{1>
=10
4
-6.400 0
0 -0.4608 -
0
1.152
3.840 ； 0
1
0 1 6.400 -0
1.152 ； 0 1.920 |
0
0.4608 -1.152
0 0 -1.152 1 3.84 0
单元②
2
_
3 4 0 0 0 -8.000
0 -1.800] 2
0 3
0.9000
0
0 -1.800 ]-0.9000 8.000 0 o 4.800 ；1.800
i 0 1.800 i 0.9000 -8.000 0 0
-1.800
0 0 2.400 4
-0.9000 0 -1.800
0 1.800 0 8.000 0 0 0
4.800 J 0
0 2.400 1.800
集成为总刚
'0.4608
对
K=104
0 -0.4608
7.300 0
称
-1.152 一
-1.800 1.152 8.0
8.461
习题9.13用先处理法计算习题9.13
EA=3.2 xl05kN , £7=4.8 xlO4kN • m2
图所示组合结构的刚度矩阵
链杆单元的£A=2.4xl05kN。
K.已知：梁杆单元的
习题9.13图
【解】离散化如习题解9.13图所示。这里利用一般单元来计算链杆单元③，令其刃为零，则
该单元的杆端转角为无意义的杆端位移，可为任意值。单元③的杆端位移编码如习题解9.13图所示, 其杆端转角在结点4处为“0”，表示无杆端转角；在结点2处为“3”, 转角相同。
点位移分量统一编码应给为“0”，再令该单元反为零。
表示与单元①和②在该端的
习题解9.13图
各单元单刚分别为
单元①和②
2(1) T 0
疔 T 1
0 2 0 0.900 1.800
0 3 12 0
4 0
3 0 0
A(1 2(2) J 1 0 2 0 3 0
_
8.000 0 0
0 i -8.000 1.800 | 4.800 !
0 0
0
=Ka>=Km =104
-0.900 1.800 -1.800 2.400
-8.000 0 0
1
0 0 i 8.000 -0.900 -1.800 1 1.800
2.400 ； 0
—
0 0 0 1 0.900 -1.800 4 2
0 3
-1.800 4.800 _
单元③
0 '3.072 0
0 0 1 2 3 2.304 0 ； -3.072 -2.304 0一 2.304
0
0
= &(3>
r = io
4
1.728 0 ； -2.304 -1.728 0 0
0 ! 0 0 0 0
-3.072 -2.304 0 1 3.072
-2.304 -1.728 0 ； 2.304
0 0 0 0
2.304
3.528
2.304 0 1
2 1.728 0
0 0 3 0 一 -0.900 -1.800 0.900
集成为总刚
-19.07
对
X=104
0
0 9.600
称
习题9.14若用先处理法计算习题9.14图所示结构，则在结构刚度矩阵K中零元素的个数至少 有多少个?
习题9.14图
【解】离散化如习题解9.14图所示，则各单元定位向量为
5(8,9,10)
2(1,2,3) @ ® _
3(4,5,6)
4(4,5,7)
6(0,0,11)“，
①
1(0,0,0)
④
7(0,0,0)
习题解9.14图
2(1) =[1 2 3 0 0 O]T, 2,2) =[1 2
3 4 5 6]丁，万3)=[8 9 10 4 5 7]T
2(4) =[4 5 7 0 0 0]T, 2(5) =[8 9 10 0 0 11]T
根据单元定位向量，判定各结点位移分量间的相关性。这里参考【例10.2］的方法，具体为: 位移分量1~3、6均与位移分量7~11无关，得到无关分量20对；位移分量4、5、7与位移分量11 无关，得到无关分量3对；合计无关分量共23对。说明K上半三角中，至少有23个元素为零，因 此整个K中至少应有46个零元素。
习题9.15试用矩阵位移法计算习题9.15图所示连续梁，并画出弯矩图。各杆常数。
6kN/m
20kN|
疝口 IIIII
\\A
4m
1 JC L
4m I 2m I 2m
习题9.15图
【解】1)离散化如习题解9.15(a)图所示。连续梁无需坐标转换。
1(0) 2(1) _ 3(2) _ 4(0)
学
-® X © X ®\"\"1
(a)离散化
习题解9.15图
2)计算总刚
各单元刚度矩阵为
0 1 1 2
2 EI
1 1/2- 0 = EI
1 1/2- 1
-EI
1 1/2-
_1/2
1 _ 1,
_l/2
1 2,
_l/2 1
集成为总刚
K = EI
2 1/2 1/2 2
3）计算综合结点荷载列阵
按照习题9.11中综合结点荷载列阵的解法,
在2、3两结点上附加刚臂, 易求得
外=[型 %]『=[8 -W]T
因此，综合结点荷载列阵为
尸=—耳=[-8 10]T
4）解结构刚度方程K4=P,
—[-5.6 EIl 6.4]T
5）求单元杆端力
根据 Fe = Ke8e + F；,得
A
（1
）
Fw
=
= K
{i）Sa>
+ Ff = El
1
1/2- 1
x —一
_l/2 1
EI o 0 --8- -10.8'
-5.6 1 + _8_ —
_ 2.4
F
（2）
= K% ⑵+ 玲少=1 1/2-
1 ~-5.6九⑶ 1 ~-2A~
Ei
_=EI
_l/2 1
1/21
- x — _l/2 1
x 1 ——EI ~6A6.4 2 2 + ■-0 —3.6
--3.6EI ~0 0 +
io10 - 13.2
一
6）绘弯矩图，如习题解9.15（b）图所示。
习题9.16用先处理法计算习题9.16图所示刚架的内力，并绘内力图。已知: £=3xlO7kN/m2, A = 0.16m2, Z= 0.002m4 „
习题9.16图
各杆
【解】1)离散化如习题解9.16(a)图所示。
(c) FQ 图(kN)
习题解9.16图
2)计算总刚
单元①无需坐标转换，其单元刚度矩阵为
0 120 0 0 -120 0 0
0 0 1.125 2.25 0 -1.125 2.25
0 0 2.25 6 0 -2.25 3
1 -120 0 0 120 0 0
2 0 -1.125 -2.25
3
0 一 0 2.25 0 3
0 1 2 3
Km=K^=104
0 0 1.125 -2.25 -2.25
6
单元②的坐标转换矩阵为
'0.8 -0.6 0.6 0.8 0 0 0
0 0 0 0
0 1 0 0 0 1 1。
0
0 0
o _ 0 0 0 0 1
0 r o.8 -0.6 0 1 0.6 0.8 0
0
0
0
则其刚度矩阵为
_
0 0 61.65 -45.8
34.93
0 0.8 1.152 4.8 -0.8 -1.152 2.4
1 -61.65 45.8 61.65 -45.8
2 3
45.8 0.8 - 0 —34.93 1.152
2.4
0 0
-45.8
K⑵
-0.8 1.152
= TK'T^ 1()4
-61.65 45.8
45.8 -34.93
T
2
-0.8 -1.152
-45.8 -0.8 1 34.93 -1.152 2
4.8
3
0.8 1.152 -0.8 -1.152
集成为总刚
181.6 -45.8 -0.8- K =\" -45.8 36.05 -3.402
-0.8 -3.402 10.8
_
3）计算综合结点荷载列阵
各单元的等效结点荷载列阵为
0 0 0 1 2 3
集成得结构的等效结点荷载列阵
理)=-留=[
0 16 10.67 0 16 -10.67]
弓=[0 16 -10.67]T
综合结点荷载列阵为
P = PE+P1=PE+[O 5 0]T =[0 21 -10.67]T
4）解结构刚度方程K4=P，得
T
4 = 1。5 [1.8513 7.4810 -7.3719]
5）求单元杆端力
根据Fe =T\\Ke5e) + F；,得
120 0 0 -120 0
0 1.125 2.25 0 -1.125 2.25
0 2.25 6 0 -2.25 3
-120 0 0 120 0 0
0 -1.125 -2.25 0 -2.25
0 2.25 3 0 6
殆, =京站。）+玲＞ =\"
1.125 -2.25
单元①定位向量 0
0 0 0 0 -16
0 0
xlO-5
0
-
-22.22- -18.50 + -10.67 — -14.56 1.8513 1 0 22.22 7.4810 2 -16 -13.50 -7.3719 3 10.67 4.56 _
0.8 -0.6 0.6 0.8
0 | 0 0 1 1 0
0
0 0 0
0 0 0 0 0 0 0 0 0 0 0 F ⑵=]⑵(K(2)於)=
0 1 0.8 -0.6 0 1 0.6 0.8
0
0
0
0
0 1 单元②定位向量
_
61.65 -45.8 0.8 -61.65 45.8 0.8 _ _ o _ o
-45.8 34.93 1.152 45.8 -34.93 1.152 0 '0 in4
0.8 1.152 4.8 -0.8 -1.152 2.4 0 0 -61.65 45.8 -0.8
61.65
-45.8
-0.8 1.8513 1 45.8 -34.93
-1.152 -45.8 34.93 -1.152 7.4810 2
0.8
1.152
2.4 -0.8 -1.152
4.8
-7.3719 3 J =[28.87 -1.47 -2.79 -28.87 1.47
-4.56 f
6)绘内力图，如习题解9.16(b)~(d)图所示。
习题9.17用矩阵位移法计算习题9.17图所示平面桁架的内力。已知：E=3xl07kN/m2
A = 0.1m2 o
【解】1)离散化如习题解9.17(a)图所示。
2)计算总刚
单元①和②无需坐标转换，其结构坐标系中的单刚分别为
0 0 1 2 0 0 3 4
■ 7.5 0 1 -7.5 o _ 0
■ 7.5 0 l -7.5 o - 0 W=km=
=
0 0 ；0 0 0 0 0 0 K105 -7.5 0 1 1 7.5 0 1， Km =
= 105
0 -7.5 0 1 0 1 7.5 0 3
0
0 1 o
0 2
0
0 1 o
0 4
单元③的a = »d,结构坐标系中的单刚为
各杆
12
3
一0 0 T
⑶=TK ⑶T = 1()6
0
0
4 0 0 1 0 0 0 -1 0
0- 1 -1 2 0 3 1
4 单元④的a = 0.35 rad，结构坐标系中的单刚为
0 ■ 3.84 2.88
0 2.88 2.16
3 -3.84 3.84 2.88
4
K(4) =rTK(4)r = io5
一0 2.88- 0 -2.88 -2.16
2.88 3 2.16 4
-3.84 -2.88 -2.88 -2.16
单元⑤的a = -0.35 rad ,结构坐标系中的单刚为
0 0 1 ■ 3.84 -2.88 -3.84
=TKT = 10
T
2 2.88 _ 0 -2.16 0 -2.88 1 2.16 2
&5
-2.88 2.16 2.88
集成总刚为
-3.84 2.88 3.84 2.88 -2.16 -2.88
11.34 -2.88 0 -2.88 12.16 0
K=105
0 0 11.34
0
0
-10 2.88 -10 2.88 12.16
—o
3) 计算综合结点荷载列阵
P = P} =[8 0
4) 解结构刚度方程K4=P，得
』=1。5 [2.5574 7.2919
0 20]T
-2.09 8.1304]T
5) 求单元杆端力
根据 Fe =T\\KeSe),
可求得各单元的杆端轴力。
这里以单元⑤为例，其杆端轴力为 [单元⑤定位向量 0.8 —0.6 0
Q8
0
2.88 -2.16 _ Lb”\")'。。
_06j
105 '3.84 -2.88 I -3.84 -2.88 2.16 1 2.88 1 -3.84 2.88 2.88 1 3.84 xlO-5
-2.88 -2.16 ； -2.88 2.16 p;7.2919 _ 13.98 _ -13.98 6)各杆轴力如习题解9.17(b)图所示。
因篇幅问题不能全部显示，请点此查看更多更全内容