结构力学习题集(下)-矩阵位移法习题及答案
结构力学课后习题解答：9矩阵位移法习题解答.docx
第9章矩阵位移法习题解答习题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 =买，则2 尹＞ =O（6 ）用矩阵位移法求得平面刚架某单元在单元坐标系中的杆端力为F e =[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离、空的值以及K ⑴中元素妍、愚、姒的值。
【解】各刚度系数的物理意义如习题解9.3图所示。
结构力学(9.14.1)--矩阵位移法习题2
5kN m
8m 8m
8m
三 . 整体分析
12. 试求图示结构 ( 不计轴变 ) 的荷载列阵 ( 先处理法 ).
1(1,0,2) 2(1,0ቤተ መጻሕፍቲ ባይዱ3) 3(1,0,3)
X1
X2
4(0,0,0)
P
X
1
0
X
2
0
四 . 求杆端力
1. 连续梁在一般荷载作用下 , 单元杆端力由下式计算 . 是否正确 ?
6
48
4
2
1(0,0,0)
12
1 6
k
6
48
4(1,0,3)
3
2(0,0,0)
3
1
2
3
例 . 不计轴变 , 作弯矩图
已知 : 各杆长均为 12m, 线刚度均为 12
P 10kN, q 5kN / m
P 10kN, q 5kN / m
解 : 1 6 1 6
k
1
6
1
48 6
6 1
24
6
6
24
6
48
3(1,0,2)
2
1
1 6 1 6 1 0
k
1
6 1
48 6
6 1
24
2
0
63 1
6 24
EI
EI
EA 2l
2 2
l
l
三 . 整体分析
4(1,0,0)
5(1,0,0)
矩阵位移法练习题一、判断题1-1、不计轴向变形，图示（a）（b）梁整体..
矩阵位移法练习题一、判断题1-1、不计轴向变形，图示（a）（b）梁整体刚度矩阵阶数相同，对应元素不同。
1-2、图示四单元的l,EI,EA相同，它们整体坐标系下的单元刚度矩阵各不相同。
1-3、矩阵位移法基本未知量的数目与位移法基本未知量的数目总是相等的。
1-4、一般单元的单元刚度矩阵一定是奇异矩阵，而特殊单元的单元刚度矩阵一定是非奇异矩阵。
1-5、如特殊单元是几何不变体系，其单元刚度矩阵一定是非奇异矩阵。
1-6、由一般单元的单元刚度方程： ,任给，并且为一平衡力系，有唯一解。
1-7、由一般单元的单元刚度方程： ,任给，有唯一解，并且为一平衡力系。
1-8、原荷载与对应的等效结点荷载产生相同的内力和变形。
1-9、在忽略轴向变形时，由单元刚度方程求出的杆端轴力为零。
应根据节点平衡由剪力求轴力。
1-10、如单元定位向量中的元素λi=0，说明该单元第 i 个杆端位移分量对应刚性支座。
二、单项选择题2-1 忽略轴向变形，用先处理法，单元①的定位向量是2-2、在图示约束情况下，单元①的单元刚度矩阵[k]=（）（ A ）2-3、图示结构单元固端弯矩列阵为，则等效结点荷载为（C ）2-4、将单元刚度矩阵分块，下列论述错误的是A 是对称矩阵B 不是对称矩阵C D2-5、在矩阵位移法中，基本未知量的确定与哪些因素无关？A 坐标系的选择B 单元如何划分C 是否考虑轴向变形D 如何编写计算机程序2-6、图示体系，忽略轴向变形，则矩阵位移法的基本未知量有几个?A 2B 3C 4D 72-7、不计轴向变形，图示（a ）（b ）梁整体刚度矩阵有何不同？A 阶数不同B 阶数相同，对应元素不同C 阶数相同，对应元素也相同D 阶数相同，仅元素k22不同2-8、不计轴向变形，图示（a ）（b ）梁整体刚度矩阵A 阶数相同，对应元素不同B 阶数相同，对应元素相同C 阶数不同，对应元素不同D 阶数不同，对应元素相同2-9、由一般单元的单元刚度方程： ,任给A 可唯一的求出，并且为一平衡力系。
《结构力学习题集》(含答案)
第三章 静定结构的位移计算一、判断题：1、虚位移原理等价于变形谐调条件，可用于求体系的位移。
2、按虚力原理所建立的虚功方程等价于几何方程。
3、在非荷载因素(支座移动、温度变化、材料收缩等)作用下，静定结构不产生内力，但会有位移且位移只与杆件相对刚度有关。
4、求图示梁铰C 左侧截面的转角时，其虚拟状态应取：A.;; B.D.C.M =15、功的互等、位移互等、反力互等和位移反力互等的四个定理仅适用于线性变形体系。
6、已知M p 、M k 图，用图乘法求位移的结果为：()/()ωω1122y y EI +。
M k M p 21y 1y 2**ωω( a )M 17、图a 、b 两种状态中，粱的转角ϕ与竖向位移δ间的关系为：δ=ϕ 。
8、图示桁架各杆E A 相同，结点A 和结点B 的竖向位移均为零。
a a9、图示桁架各杆EA =常数，由于荷载P 是反对称性质的，故结点B 的竖向位移等于零。
二、计算题：10、求图示结构铰A 两侧截面的相对转角ϕA ，EI = 常数。
q l l l /211、求图示静定梁D 端的竖向位移 ∆DV 。
EI = 常数 ，a = 2m 。
a a a 10kN/m12、求图示结构E 点的竖向位移。
EI = 常数 。
l l l l /3 2 /3/3q13、图示结构，EI=常数 ，M =⋅90kN m , P = 30kN 。
求D 点的竖向位移。
P 3m 3m 3m14、求图示刚架B 端的竖向位移。
q15、求图示刚架结点C 的转角和水平位移，EI = 常数 。
q16、求图示刚架中Ｄ点的竖向位移。
EI ＝常数。
l ll/217、求图示刚架横梁中Ｄ点的竖向位移。
EI＝常数。
18、求图示刚架中D点的竖向位移。
E I = 常数。
qll l/2219、求图示结构Ａ、Ｂ两截面的相对转角，EI＝常数。
l/23l/320、求图示结构A、B两点的相对水平位移，E I = 常数。
ll21、求图示结构B点的竖向位移，EI = 常数。
《结构力学习题集》-矩阵位移法习题及答案
第八章 矩阵位移法 – 老八校一、判断题：1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。
2、单元刚度矩阵均具有对称性和奇异性。
3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。
4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。
5、结构刚度方程矩阵形式为：[]{}{}K P ∆=，它是整个结构所应满足的变形条件。
6、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。
7、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。
8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。
9、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。
10、矩阵位移法既能计算超静定结构，也能计算静定结构。
11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是：(0,1,2)(0,0,0)(0,0,0)(0,1,3)(0,0,0)(1,2,0)(0,0,0)(0,0,3)(1,0,2)(0,0,0)(0,0,0)(1,0,3)(0,0,0)(0,1,2)(0,0,0)(0,3,4)A.B.C.D.2134123412341234( )二、计算题：12、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。
123ll4l5EI2EIEA(0,0,0)(0,0,1)(0,2,3)(0,0,0)(0,2,4)(0,0,0)EI13、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。
EI ，EA 均为常数。
l14、计算图示结构整体刚度矩阵的元素665544,,K K K 。
E 为常数。
ll1342A , I AA /222A I , 2A15、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵[][]K K 2224,。
[][]k k 1112 [][]k k 2122 []k =ii iii单刚分块形式为 ：16、已知平面桁架单元在整体坐标系中的单元刚度矩阵，计算图示桁架结构原始刚度矩阵[]K 中的元素,,7877K K EA ＝常数。
矩阵位移法例题
0
2 1 2
0
0
4 1 3
00 2 00 3
0
0
K③
41
3
0
0
0
00 3 000
5 集成总刚度矩阵
第8章矩阵位移法
4 2 2 2
0 1 8 4 0
K 2 2 4 2 4 1
21
2
4
12
2
0
2 1 4 1 4 1 3 0 2 8
1
2
3
6 形成荷载向量
P 60 190 62.5T
2 结点位移编号矩阵 3 各单元旳定位向量
0 0 0 C 0 0 1
0 0 2 0 0 0
2 3T
U1 0 0 0 0 0 1 U2 0 0 1 0 0 2 U3 0 0 2 0 0 0
-90 250
-250 187.5 -112.5
1
2
3
4
第8章矩阵位移法
4 各单元旳刚度矩阵
单元旳刚度矩阵与解法一相同
2 12i 2 BCx l2 Cy
12i (B l2 )CxC y
2 12i 2
BC Y
2 l
Cx
6i l Cy 6i l Cx
2 12i 2 BCx 2 C y
l 12i (B 2 )CxC y l
12i (B 2 )CxC y
l 2 12i 2 BCy 2 Cx
l
6i l Cy 6i l Cx
(e)
K
6i
4i
l Cy
6i l Cx
2i
2 12i 2 BCx 2 C y
l
12i (B 2 )CxC y
l
6i
(参考资料)结构力学答案(下册)
2l 2 -3l
-3l 6
l2
⎥ ⎥
-3l ⎥
⎥
⎢⎣3l l 2 -3l 2l2 ⎥⎦
（4）总刚度矩阵
12
3
4
1
2
3
⎡12 6l -12 6l 0 0
⎢⎢6l 4l 2 -6l 2l 2 0 0
⎢⎡k11① k12 ①
kθ
=
⎢k ⎢
① 21
⎢0
k 22 ① + k 22 ② k32 ②
0
0⎤
⎢-12
2m
解：（1）结构标识 y
②3 ③2
①
1x
单元 局部坐标系（ i → j ） 杆长
① 1→2
2
② 2→3
2
cosα
0
3 2
sin α
1
1 2
（2）建立结点位移向量，结点力向量
[ ] ∆ = µ2 ν 2 θ2 θ3 T
F = [20 0 - 30 0]T
（3）建立单元刚度矩阵（l=2m）
µ2
⎡12 EI
（3）计算单元刚度矩阵
1
2
⎡12 6l -12 6l ⎤
k①
= ⎢⎡k11①
⎢⎣k
① 21
k12 ① k 22 ①
⎤ ⎥ ⎥⎦
=
2EI l3
⎢⎢6l ⎢−12 ⎢
4l 2 -6l
-6l
2l
2
⎥ ⎥
12 -6l ⎥
⎥
⎢⎣6l 2l 2 -6l 4l 2 ⎥⎦
1
2
⎡6 3l -6 3l ⎤
k ② = ⎢⎡k22② ⎢⎣k32②
③
x
1
2
结构力学位移法题及答案
超静定结构计算——位移法一、判断题：1、判断下列结构用位移法计算时基本未知量的数目。
（1） （2） （3）（4） （5） （6）EIEIEIEI 2EI EI EIEIEAEA ab EI=EI=EI=244422、位移法求解结构内力时如果P M 图为零，则自由项1P R 一定为零。
3、位移法未知量的数目与结构的超静定次数有关。
4、位移法的基本结构可以是静定的，也可以是超静定的。
5、位移法典型方程的物理意义反映了原结构的位移协调条件。
二、计算题：12、用位移法计算图示结构并作M 图，横梁刚度EA →∞，两柱线刚度 i 相同。
213、用位移法计算图示结构并作M 图。
E I =常数。
—— 41 ——ll /2l /214、求对应的荷载集度q 。
图示结构横梁刚度无限大。
已知柱顶的水平位移为 ()5123/()EI →。
12m12m8mq15、用位移法计算图示结构并作M 图。
EI =常数。
ll ll16、用位移法计算图示结构，求出未知量，各杆EI 相同。
4m19、用位移法计算图示结构并作M 图。
qll20、用位移法计算图示结构并作M 图。
各杆EI =常数，q = 20kN/m 。
6m6m23、用位移法计算图示结构并作M 图。
EI =常数。
ll 224、用位移法计算图示结构并作M 图。
EI =常数。
q29、用位移法计算图示结构并作M 图。
设各杆的EI 相同。
qql l /2/232、用位移法作图示结构M 图。
E I =常数。
—— 43 ——qql l/2l /2l36、用位移法计算图示对称刚架并作M 图。
各杆EI =常数。
l l38、用位移法计算图示结构并作M 图。
EI =常数。
ql l l l42、用位移法计算图示结构并作M 图。
2m 2m43、用位移法计算图示结构并作M 图。
EI =常数。
lllql48、已知B 点的位移∆，求P 。
ll/2/2A∆51、用位移法计算图示结构并作M 图。
超静定结构计算——位移法（参考答案）1、（1）、4； （2）、4； （3）、9； （4）、5； （5）、7；（6）、7。
《结构力学习题集》(含答案)
第三章 静定结构的位移计算一、判断题：1、虚位移原理等价于变形谐调条件，可用于求体系的位移。
2、按虚力原理所建立的虚功方程等价于几何方程。
3、在非荷载因素(支座移动、温度变化、材料收缩等)作用下，静定结构不产生内力，但会有位移且位移只与杆件相对刚度有关。
4、求图示梁铰C 左侧截面的转角时，其虚拟状态应取：A.;;B.D.C.=1=15、功的互等、位移互等、反力互等和位移反力互等的四个定理仅适用于线性变形体系。
6、已知M p 、M k 图，用图乘法求位移的结果为：()/()ωω1122y y EI +。
M kM p21y 1y 2**ωω( a )M =17、图a 、b 两种状态中，粱的转角ϕ与竖向位移δ间的关系为：δ=ϕ 。
8、图示桁架各杆E A 相同，结点A 和结点B 的竖向位移均为零。
aa9、图示桁架各杆EA =常数，由于荷载P 是反对称性质的，故结点B 的竖向位移等于零。
二、计算题：10、求图示结构铰A 两侧截面的相对转角ϕA ，EI = 常数。
qlll /211、求图示静定梁D 端的竖向位移 ∆DV 。
EI = 常数 ，a = 2m 。
a a a10kN/m12、求图示结构E 点的竖向位移。
EI = 常数 。
ll l /32 /3/3q13、图示结构，EI=常数 ，M =⋅90kN m , P = 30kN 。
求D 点的竖向位移。
P 3m3m3m14、求图示刚架B 端的竖向位移。
q15、求图示刚架结点C 的转角和水平位移，EI = 常数 。
q16、求图示刚架中Ｄ点的竖向位移。
EI ＝常数。
l ll/217、求图示刚架横梁中Ｄ点的竖向位移。
EI＝常数。
18、求图示刚架中D点的竖向位移。
E I = 常数。
qll l/2219、求图示结构Ａ、Ｂ两截面的相对转角，EI＝常数。
l/23l/320、求图示结构A、B两点的相对水平位移，E I = 常数。
ll21、求图示结构B点的竖向位移，EI = 常数。
结构力学大作业(矩阵位移法)
矩阵位移法编程大作业姓名：学号：一、编程原理本程序的原理是基于结构力学矩阵位移法原理，以结构结点位移作基本未知量，将要分析的结构拆成已知节点力—结点力位移关系的单跨梁集合，通过强令结构发生待定的基本未知位移，在各个单跨梁受力分析结果的基础上通过保证结构平衡建立位移法的线性方程组，从而求得基本未知量。
二、程序说明本程序是计算10个节间距的悬索-拱组合体系主塔顶节点水平位移、主塔底截面弯矩、拱顶节点竖向位移、拱顶截面弯矩和轴力的程序。
首先将各杆件的交汇点作为结点，共有20个结点，51个位移，然后根据不同结构单元分别建立单元刚度矩阵，然后转换为整体坐标系下的刚度矩阵，然后将所有杆件的单元刚度矩阵整合成为总体刚度矩阵，在进行整合时连续运用for函数，最终形成51阶的总体刚度矩阵。
然后通过对荷载的分析确定出荷载矩阵，直接写进程序。
这样就可以把20个结点的51个位移求得，然后再利用各个单元的单元刚度矩阵和所得的位移求得单元杆件的内力。
三、算法流程建立各单位在局部结构离散化编号进行单元分析坐标系下的单位刚度方程确定各单位在总体将单元刚度矩阵集合确定综合结点坐标系下的单元矩阵方程成总体刚度矩阵点荷载矩阵建立方程利用杆件单元刚度矩阵输出结果求解位移和所求位移求内力结束四、源代码L=input('输入单节间L：');EIc=input('主塔的抗弯刚度EIc:');EAc=input('主塔的抗压刚度EAc:');EAb=input('悬索和斜索的抗拉刚度EAb:');EAt=input('吊杆的抗拉刚度EAt:');EIa=input('拱的抗弯刚度EIa:');EAa=input('拱的抗压刚度EAa:');q=input('拱上沿轴向均布荷载集度q:');T1=[0,1,0,0,0,0;-1,0,0,0,0,0;0,0,1,0,0,0;0,0,0,0,1,0;0,0,0,-1,0,0;0,0,0,0,0,1;];%主塔的转换矩阵h=(5*L)/2;KcO=[EAc/h,0,0,-EAc/h,0,0;0,12*EIc/(h*h*h),6*EIc/(h*h),0,-12*EIc/(h*h*h),6*EIc/(h*h);0,6*EIc/(h*h),4*EIc/h,0,-6*EIc/(h*h),2*EIc/h;-EAc/h,0,0,EAc/h,0,0;0,-12*EIc/(h*h*h),-6*EIc/(h*h),0,12*EIc/(h*h*h),-6*EIc/(h*h);0,6*EIc/(h*h),2*EIc/h,0,-6*EIc/(h*h),4*EIc/h;];%主塔的单元刚度矩阵x=atan(2*L/h);T2=[cos(x),sin(x),0,0;-sin(x),cos(x),0,0;0,0,cos(x),sin(x);0,0,-sin(x),cos(x);];y=-atan(2*L/h);T21=[cos(y),sin(y),0,0;-sin(y),cos(y),0,0;0,0,cos(y),sin(y);0,0,-sin(y),cos(y);];%斜索的转换矩阵s1=sqrt(2*L*2*L+h*h);KbO1=(EAb/s1)*[1 0 -1 0;0 0 0 0;-1 0 1 0;0 0 0 0;];%斜索的单元刚度矩阵f2(1)=5*L/2;f2(2)=58*L/25;f2(3)=109*L/50;f(4)=52*L/25;f2(5)=101*L/50;f2 (6)=2*L;f2(7)=101*L/50;f2(8)=52*L/25;f2(9)=109*L/50;f2(10)=58*L/25;f2(1 1)=5*L/2;y=zeros(10,1);for i=1:10y(i)=atan((f2(i+1)-f2(i))/L);endT3=zeros(4,40);for i=1:10T3(1:4,4*i-3:4*i)=[cos(y(i)),sin(y(i)),0,0;-sin(y(i)),cos(y(i)),0,0;0,0,cos(y(i)),sin(y(i));0,0,-sin(y(i)),cos(y(i));];end%悬索的转换矩阵s2=zeros(10,1);for i=1:10s2(i)=sqrt((f2(i+1)-f2(i))^2+L^2);endKbO2=zeros(4,40);KbO2(1:4,4*i-3:4*i)=(EAb/s2(i))*[1 0 -1 0;0 0 0 0;-1 0 1 0;0 0 0 0;];end%悬索的单元刚度矩阵f1(1)=0;f1(2)=9*L/20;f1(3)=4*L/5;f1(4)=21*L/20;f1(5)=6*L/5;f1(6)=5*L/4; f1(7)=6*L/5;f1(8)=21*L/20;f1(9)=4*L/5;f1(10)=9*L/20;f1(11)=0;z=zeros(10,1);for i=1:10z(i)=atan((f1(i+1)-f1(i))/L);endT4=zeros(6,60);for i=1:10T4(6*i-5:6*i,6*i-5:6*i)=[cos(z(i)),sin(z(i)),0,0,0,0;-sin(z(i)),cos(z(i)),0,0,0,0;0,0,1,0,0,0;0,0,0,cos(z(i)),sin(z(i)),0;0,0,0,-sin(z(i)),cos(z(i)),0;0,0,0,0,0,1;];end%拱的转换矩阵s3=zeros(10,1);for i=1:10s3(i)=sqrt((f1(i+1)-f1(i))^2+L^2);endKaO=zeros(6,60);for i=1:10KaO(1:6,6*i-5:6*i)=[EAa/s3(i) 0 0 -EAa/s3(i) 0 0;0 12*EIa/(s3(i)*s3(i)*s3(i)) 6*EIa/(s3(i)*s3(i)) 0-12*EIa/(s3(i)*s3(i)*s3(i)) 6*EIa/(s3(i)*s3(i));0 6*EIa/(s3(i)*s3(i)) 4*EIa/s3(i) 0 -6*EIa/(s3(i)*s3(i)) 2*EIa/s3(i);-EAa/s3(i) 0 0 EAa/s3(i) 0 0;0 -12*EIa/(s3(i)*s3(i)*s3(i)) -6*EIa/(s3(i)*s3(i)) 012*EIa/(s3(i)*s3(i)*s3(i)) -6*EIa/(s3(i)*s3(i));0 6*EIa/(s3(i)*s3(i)) 2*EIa/s3(i) 0 -6*EIa/(s3(i)*s3(i)) 4*EIa/s3(i);]; end%拱的单元刚度矩阵T5=[0 1 0 0;-1 0 0 0;0 0 0 1;0 0 -1 0;];%吊杆的转换矩阵s4=zeros(9,1);s4(i)=f2(i+1)-f1(i+1);endKtO=zeros(4,36);for i=1:9KtO(1:4,4*i-3:4*i)=(EAt/s4(i))*[1 0 -1 0;0 0 0 0;-1 0 1 0;0 0 0 0;];end%吊杆的单元刚度矩阵Kc=T1'*KcO*T1;%总体坐标下主塔的单元刚度矩阵Kb1=T2'*KbO1*T2;Kb11=T21'*KbO1*T21;%总体坐标下斜索的单元刚度矩阵Kb2=zeros(4,40);for i=1:10T3O=T3(1:4,4*i-3:4*i);Kb2(1:4,4*i-3:4*i)=T3O'*KbO2(1:4,4*i-3:4*i)*T3O;end%总体坐标下悬索的单元刚度矩阵Ka=zeros(6,60);for i=1:10T4O=T4(6*i-5:6*i,6*i-5:6*i);Ka(1:6,6*i-5:6*i)=T4O'*KaO(1:6,6*i-5:6*i)*T4O;end%总体坐标下拱的单元刚度矩阵Kt=zeros(4,36);for i=1:9KtOO=KtO(1:4,4*i-3:4*i);Kt(1:4,4*i-3:4*i)=T5'*KtOO*T5;end%总体坐标下吊杆的单元刚度矩阵%定义51阶0矩阵K1=zeros(51,51);K2=zeros(51,51);K3=zeros(51,51);K4=zeros(51,51);K5=zero s(51,51);X=zeros(51,51);Y=zeros(51,51);Z=zeros(51,51);%把主塔整合到整体刚度矩阵中：K1(1:3,1:3)=KcO(4:6,4:6);K1(22:24,22:24)=KcO(4:6,4:6);%把斜索整合到整体刚度矩阵中：K2(1:2,1:2)=Kb1(3:4,3:4);K2(22:23,22:23)=Kb11(1:2,1:2);%把悬索整合到整体刚度矩阵中:K3(1:2,1:2)=KbO2(1:2,1:2);K3(1:2,4:5)=KbO2(1:2,3:4);for i=2:10X(2*i:2*i+3,2*i:2*i+3)=KbO2(1:4,4*i-3:4*i);K3=K3+X;end%把拱整合到整体刚度矩阵中：K4(25:27,25:27)=KaO(4:6,4:6);K4(49:51,49:51)=KaO(1:3,55:57);for i=2:9Y(3*i+19:3*i+24,3*i+19:3*i+24)=KaO(1:6,6*i-5:6*i); K4=K4+Y;end%把吊杆整合到整体刚度矩阵中：for i=1:9Z(2*i+2:2*i+3,2*i+2:2*i+3)=KtO(1:2,1:2);Z(2*i+2:2*i+3,3*i+22:3*i+23)=KtO(1:2,3:4);Z(3*i+22:3*i+23,2*i+2:2*i+3)=KtO(3:4,1:2);Z(3*i+22:3*i+23,3*i+22:3*i+23)=KtO(3:4,3:4);K5=K5+Z;endK=K1+K2+K3+K4+K5;%荷载矩阵：P=zeros(51,1);P(26,1)=-q*L/(2*cos(s3(1)));P(27,1)=q*L*L/(12*cos(s3(1)));P(50,1)=-q*L/(2*cos(s3(10)));P(51,1)=-q*L*L/(12*cos(s3(10)));for i=2:9P0=zeros(51,1);P0(3*i+20,1)=-q*L/(2*cos(s3(i)));P0(3*i+21,1)=-q*L*L/(12*cos(s3(i)));P0(3*i+23,1)=-q*L/(2*cos(s3(i)));P0(3*i+24,1)=q*L*L/(12*cos(s3(i)));P=P+P0;endA=K\P;%结构的位移%主塔底截面的弯矩：Ac(4:6,1)=A(1:3,1);Bc=KcO*Ac;Mc=Bc(3,1);%拱顶截面的弯矩和轴力：Aa=A(34:39,1);KaO17=KaO(1:6,25:30);Ba=KaO17*Aa;Ma=Ba(6,1);Fa=Ba(4,1);%输出结果fprintf('主塔顶结点的水平位移%f\n',A(1,1)); fprintf('主塔底截面的弯矩%f\n',Mc);fprintf('拱顶结点的竖向位移%f\n',A(38,1)); fprintf('拱顶截面的弯矩%f\n',Ma);fprintf('拱顶截面的轴力%f\n',Fa);五、试算算例输入单节间L：1主塔的抗弯刚度EIc:1主塔的抗压刚度EAc:1悬索和斜索的抗拉刚度EAb:1吊杆的抗拉刚度EAt:1拱的抗弯刚度EIa:1拱的抗压刚度EAa:1拱上沿轴向均布荷载集度q:1主塔顶结点的水平位移NaN主塔底截面的弯矩NaN拱顶结点的竖向位移0.016046拱顶截面的弯矩3.791098拱顶截面的轴力0.000000。
《结构力学习题集》(下)-矩阵位移法习题及答案 (2)
第七章 矩阵位移法一、就是非题1、单元刚度矩阵反映了该单元杆端位移与杆端力之间得关系。
2、单元刚度矩阵均具有对称性与奇异性。
3、局部坐标系与整体坐标系之间得坐标变换矩阵T 就是正交矩阵。
4、结构刚度矩阵反映了结构结点位移与荷载之间得关系。
5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。
6、结 构 刚 度 矩 阵 就是 对 称 矩 阵 ,即 有K i j = K j i ,这 可 由 位 移 互 等 定 理 得 到 证 明 。
7、结构刚度方程矩阵形式为:,它就是整个结构所应满足得变形条件。
8、在直接刚度法得先处理法中,定位向量得物理意义就是变形连续条件与位移边界条件。
9、等效结点荷载数值等于汇交于该结点所有固端力得代数与。
10、矩阵位移法中,等效结点荷载得“等效原则”就是指与非结点荷载得结点位移相等。
11、矩阵位移法既能计算超静定结构,也能计算静定结构。
二、选择题1、已知图示刚架各杆EI = 常数,当只考虑弯曲变形,且各杆单元类型相同时,采用先处理法进行结点位移编号,其正确编号就是:(0,1,2)(0,0,0)(0,0,0)(0,1,3)(0,0,0)(1,2,0)(0,0,0)(0,0,3)(1,0,2)(0,0,0)(0,0,0)(1,0,3)(0,0,0)(0,1,2)(0,0,0)(0,3,4)A.B.C.D.21341234123412342、平面杆件结构一般情况下得单元刚度矩阵,就其性质而言,就是:A.非对称、奇异矩阵;B.对称、奇异矩阵;C.对称、非奇异矩阵;D.非对称、非奇异矩阵。
3、单元i j 在图示两种坐标系中得刚度矩阵相比:A.完全相同;B.第2、3、5、6行(列)等值异号;C.第2、5行(列)等值异号;D.第3、6行(列)等值异号。
4、矩阵位移法中,结构得原始刚度方程就是表示下列两组量值之间得相互关系:A.杆端力与结点位移;B.杆端力与结点力;C.结点力与结点位移;D.结点位移与杆端力。
结构力学习题集矩阵位移法习题及答案老八校
1文档收集于互联网，已整理，word 版本可编辑.第八章 矩阵位移法 – 老八校一、判断题：1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。
2、单元刚度矩阵均具有对称性和奇异性。
3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。
4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。
5、结构刚度方程矩阵形式为：[]{}{}K P ∆=，它是整个结构所应满足的变形条件。
6、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。
7、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。
8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。
9、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。
10、矩阵位移法既能计算超静定结构，也能计算静定结构。
11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是： 二、计算题：12、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。
13、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。
EI ，EA 均为常数。
14、计算图示结构整体刚度矩阵的元素665544,,K K K 。
E 为常数。
15、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵[][]K K 2224,。
16、已知平面桁架单元在整体坐标系中的单元刚度矩阵，计算图示桁架结构原始刚度矩阵[]K 中的元素,,7877K K EA ＝常数。
,cos α=C ,sin α=S ,C C A ⋅= S S D S C B ⋅=⋅=,，各杆EA 相同。
2文档收集于互联网，已整理，word 版本可编辑.17、计算图示刚架结构刚度矩阵中的元素8811,K K （只考虑弯曲变形）。
设各层高度为h ，各跨长度为l h l 5.0,=，各杆EI 为常数。
18、计算图示结构原始刚度矩阵的元素4544,K K 。
《结构力学习题集》(含答案)
第三章 静定结构的位移计算一、判断题：1、虚位移原理等价于变形谐调条件，可用于求体系的位移。
2、按虚力原理所建立的虚功方程等价于几何方程。
3、在非荷载因素(支座移动、温度变化、材料收缩等)作用下，静定结构不产生内力，但会有位移且位移只与杆件相对刚度有关。
4、求图示梁铰C 左侧截面的转角时，其虚拟状态应取：A.;;B.D.M C.=1=1=15、功的互等、位移互等、反力互等和位移反力互等的四个定理仅适用于线性变形体系。
6、已知M p 、M k 图，用图乘法求位移的结果为：()/()ωω1122y y EI +。
M kM p21y 1y 2**ωω( a )M =17、图a 、b 两种状态中，粱的转角ϕ与竖向位移δ间的关系为：δ=ϕ 。
8、图示桁架各杆E A 相同，结点A 和结点B 的竖向位移均为零。
aa9、图示桁架各杆EA =常数，由于荷载P 是反对称性质的，故结点B 的竖向位移等于零。
二、计算题：10、求图示结构铰A 两侧截面的相对转角ϕA ，EI = 常数。
qlll /211、求图示静定梁D 端的竖向位移 ∆DV 。
EI = 常数 ，a = 2m 。
a a a10kN/m12、求图示结构E 点的竖向位移。
EI = 常数 。
ll l l /32 /3/3q13、图示结构，EI=常数 ，M =⋅90kN m , P = 30kN 。
求D 点的竖向位移。
P 3m3m3m14、求图示刚架B 端的竖向位移。
q15、求图示刚架结点C 的转角和水平位移，EI = 常数 。
q16、求图示刚架中Ｄ点的竖向位移。
EI ＝ 常数 。
l/217、求图示刚架横梁中Ｄ点的竖向位移。
EI ＝ 常数 。
18、求图示刚架中D 点的竖向位移。
E I = 常数 。
qll/219、求图示结构Ａ、Ｂ两截面的相对转角，EI ＝ 常数 。
l/23l/320、求图示结构A 、B 两点的相对水平位移，E I = 常数。
ll21、求图示结构B 点的竖向位移，EI = 常数。
《结构力学习题集》(下)-矩阵位移法习题及答案
第八章 矩阵位移法1、(O)2、(X)3、(O)4、(X)5、(X)6、(O)7、(O)8、(X)9、(O) 10、(O) 11、(A)一、判断题：1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。
2、单元刚度矩阵均具有对称性和奇异性。
3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。
4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。
5、结构刚度方程矩阵形式为：[]{}{}K P ∆=，它是整个结构所应满足的变形条件。
6、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。
7、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。
8、等效结点荷载数值等于汇交于该结点所有固端力的代数和。
9、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。
10、矩阵位移法既能计算超静定结构，也能计算静定结构。
11、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是：(0,1,2)(0,0,0)(0,0,0)(0,1,3)(0,0,0)(1,2,0)(0,0,0)(0,0,3)(1,0,2)(0,0,0)(0,0,0)(1,0,3)(0,0,0)(0,1,2)(0,0,0)(0,3,4)A.B.C.D.2134123412341234xy M , θ( )二、计算题：12、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。
123ll4ll5EI2EIEA(0,0,0)(0,0,1)(0,2,3)(0,0,0)(0,2,4)(0,0,0)xyM , θEI13、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。
EI ，EA 均为常数。
l(0,0,1)(0,5,0)(2,3,4)l①②123xy M , θ14、计算图示结构整体刚度矩阵的元素665544,,K K K 。
习题课1 矩阵位移法(含答案作业)518706462
k44 k45 k46
k54 k55 k56
0 0 1 0 02
k64 k65 k66
0 k11 k12 k13
k16
0 k21 k22 k23
k26
[ ] = k 2 1 k31 k32 k33
k36
0
0
2 k61 k62 k63
k686
0 0 1 3 04
0 k11 k12 k13 k14 k15 k16 0 k21 k22 k23 k24 k25 k26
（0,0,0） （0,0,1） （0,0,2）
（0,0,0） （1,2,3） （0,0,4）
1
2
3
①
②
x③
y
（（53,,60,,85） ）5
4（3,0,4） （5,6,7）
④
（0,0,0） 6（0,0,0）
不考虑轴向变形 考虑轴向变形
7
（1) 不考虑轴向变形
0 0 0 0 01 0 0
0
[k ] 1= 0 0 1
(↑↓)
16
分别作上述两种情况下的弯矩图，如下图示。 据此容易得出含铰单元的刚度矩阵[K]。
3EIa2 a3 + b3 A
3EIab a3 + b3
B
3EIab a3 + b3 A
3EIb2 a3 + b3
B
3EIa
−3EIa
a3 + b3
a3 + b3
θ
e A
=1
a2
[k]e =
ab
3EIb
−3EIb
k k (2) (2) 63 66
0
0
00
k k k 5 (3) (3) (3) 41 42 43
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第七章 矩阵位移法一、是非题1、单元刚度矩阵反映了该单元杆端位移与杆端力之间的关系。
2、单元刚度矩阵均具有对称性和奇异性。
3、局部坐标系与整体坐标系之间的坐标变换矩阵T 是正交矩阵。
4、结构刚度矩阵反映了结构结点位移与荷载之间的关系。
5、用 矩 阵 位 移 法 计 算 连 续 梁 时 无 需 对 单 元 刚 度 矩 阵 作 坐 标 变 换。
6、结 构 刚 度 矩 阵 是 对 称 矩 阵 ，即 有K i j = K j i ，这 可 由 位 移 互 等 定 理 得 到 证 明 。
7、结构刚度程矩阵形式为：[]{}{}K P ∆=，它是整个结构所应满足的变形条件。
8、在直接刚度法的先处理法中，定位向量的物理意义是变形连续条件和位移边界条件。
9、等效结点荷载数值等于汇交于该结点所有固端力的代数和。
10、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。
11、矩阵位移法既能计算超静定结构，也能计算静定结构。
二、选择题1、已知图示刚架各杆EI = 常数，当只考虑弯曲变形，且各杆单元类型相同时，采用先处理法进行结点位移编号，其正确编号是：(0,1,2)(0,0,0)(0,0,0)(0,1,3)(0,0,0)(1,2,0)(0,0,0)(0,0,3)(1,0,2)(0,0,0)(0,0,0)(1,0,3)(0,0,0)(0,1,2)(0,0,0)(0,3,4)A.B.C.D.21341234123412342、平面杆件结构一般情况下的单元刚度矩阵[]k 66⨯，就其性质而言，是： A ．非对称、奇异矩阵； B ．对称、奇异矩阵；C ．对称、非奇异矩阵；D ．非对称、非奇异矩阵。
3、单元i j 在图示两种坐标系中的刚度矩阵相比：A ．完全相同；B ．第2、3、5、6行(列)等值异号；C ．第2、5行(列)等值异号；D ．第3、6行(列)等值异号。
xi4、矩阵位移法中，结构的原始刚度程是表示下列两组量值之间的相互关系： A ．杆端力与结点位移； B ．杆端力与结点力； C ．结点力与结点位移； D ．结点位移与杆端力 。
5、单 元 刚 度 矩 阵 中 元 素 k ij 的 物 理 意 义 是 ：A ．当 且 仅 当 δi =1 时 引 起 的 与 δj 相 应 的 杆 端 力 ；B ．当 且 仅 当 δj =1时 引 起 的 与 δi 相 应 的 杆 端 力 ；C ．当 δj =1时 引 起 的 δi 相 应 的 杆 端 力 ；D ．当 δi =1时 引 起 的 与 δj 相 应 的 杆 端 力。
三、填充题1、图示结构用矩阵位移法计算时（计轴向变形）未知量数目为8个。
2、图 示 刚 架 用 两 种 式 进 行 结 点 编 号 ，结 构 刚 度 矩 阵 最 大 带 宽 较 小 的 是 图。
35641271234567(a)(b)3、图 示 梁 结 构 刚 度 矩 阵 的 主 元 素 K K 1122== , 。
ll4、图 示 桁 架 结 构 刚 度 矩 阵 有 个 元 素 ，其 数 值 等 于 。
3m3m ABC DEAEAEA5、用 矩 阵 位 移 法 解 图 示 连 续 梁 时 ，结 构 的 综 合 结 点 荷 载 是l /2ll l /26、已知图示桁架杆件①的单元刚度矩阵为式(a)，又已知各结点位移为式(b)，则杆件①的轴力（注明拉力或压力）应为N①= 。
l[]k EA l u v u v u v u v Pl EA ①=--⎡⎣⎢⎢⎢⎢⎤⎦⎥⎥⎥⎥⎧⎨⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎫⎬⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪=-⎧⎨⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎫⎬⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪1(a) b)010*********005100230011223344ΛΛΛΛ(四、计算题1、用先处理法写出图示梁的整体刚度矩阵[]K 。
123llli 0123i i2、用先处理法写出图示梁的结构刚度矩阵[]K 。
123ll4lEI EI EI 233、计算图示结构的综合结点荷载列阵{}P 。
l /2l /2l /2l /2ll4、计算图示连续梁对应于自由结点位移的荷载列阵{}P 。
l /2l l /25、已 知 图 示 连 续 梁结 点 位 移 列 阵 {}θ如 下 所示 ，试 用 矩 阵 位 移 法 求 出 杆 件 23的 杆 端 弯 矩 并 画 出 连 续 梁 的 弯 矩 图 。
设 q = 20kN /m ，23 杆 的 i =⨯⋅10106.kN cm 。
{}θ=--⎧⎨⎪⎪⎩⎪⎪⎫⎬⎪⎪⎭⎪⎪⨯-365714572286104....rad6、已知图示梁结点转角列阵为{}[]∆=056516822-/ /Tql i ql i ，EI =常数。
计算B 支座的反力。
1m1m7、试 用 矩 阵 位 移 法 解 图 示 连 续 梁 ，绘 弯 矩 图 。
EI = 已 知 常 数 。
xθ8、试 求 结 构 原 始 刚 度 矩 阵 中 的 子 块 []K 22 ，已 知 单 元 ①的 整 体 坐 标 的 单 元 刚度 矩 阵 为 ：[]K ①=-⨯-⨯---⨯-⨯⎡⎣⎢⎢⎢⎢⎤⎦⎥⎥⎥⎥72360072360036003600723600360036001103600210442101107244ll9、用先处理法写出图示结构的结构刚度矩阵[]K 。
E =常数。
ll10、用先处理法计算图示连续梁的结点荷载列阵{}P 。
m4m4m411、计算图示连续梁对应于自由结点位移的荷载列阵{}P 。
m3m3m 4m 412、已 知 图 示 两 端 固 定 梁 跨 中 结 点 C 的 竖 向 位 移 为 ∆CV l EI =-5123() ，转 角 ϕC =0 ，l =5m ，EI =常 数 。
试 求 单 元 ① 、② 的 杆 端 力 列 阵 。
ll13、用先处理法计算图示结构刚度矩阵的元素133322,,K K K 。
123ll4l5EI2EIEA(0,0,0)(0,0,1)(0,2,3)(0,0,0)(0,2,4)(0,0,0)EI14、用先处理法计算图示刚架结构刚度矩阵的元素153422,,K K K 。
EI ，EA 均为常数。
l15、写出图示结构以子矩阵形式表达的结构原始刚度矩阵的子矩阵[][]K K 2224,。
[][]k k 1112 [][]k k 2122 []k =ii iii单刚分块形式为 ：16、已知图示结构在整体坐标系中的单元刚度矩阵。
用先处理法集成结构刚度矩阵[]K 。
（用子块形式写出）。
[][]k k 1112 [][]k k 2122 []k =ii iii单刚分块形式为 ：17、用先处理法写出图示刚架的结构刚度矩阵[]K ，只考虑弯曲变形。
EI EI EIEI=o ol ll18、用先处理法写出图示结构的结构刚度矩阵[]K 。
各杆长度为l ，EA 、EI 为常数。
ABCD19、用先处理法写出以子块表示的图示结构的结构刚度矩阵[]K 。
m12m20、用先处理法写出图示刚架结构刚度矩阵[]K 。
已知：[][][]k k k ①②③===⨯--------⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥10300003000001230012300301000305030000300000123001230030500301004x21、计算图示结构结点3的等效结点荷载列阵{}P 3E 。
2222、计算图示结构结点2的等效结点荷载列阵{}P 2E 。
l /2l /2q23、计算图示结构的综合结点荷载列阵元素431,,P P P 。
l24、用先处理法计算图示结构的综合结点荷载列阵{}P 。
l/2/225、计算图示结构结点荷载列阵中的元素654,,P P P 。
l /2l /2(0,7,8)326、计算图示结构综合结点荷载列阵中的元素431,,P P P 。
l l l27、计算图示结构综合结点荷载列阵{}P 中的元素9873,,,P P P P 。
ll228、计算图示刚架对应于自由结点位移的综合结点荷载列阵{}P 。
m3m3m29、计算图示刚架对应自由结点位移的综合结点荷载列阵{}P 。
各杆长度为 4m 。
30、计算图示结构结点2的综合结点荷载列阵{}P 2。
l /2l l /2l l31、计算图示刚架考虑弯曲、轴向变形时的综合结点荷载列阵{}P 。
32、若考虑弯曲、轴向变形，用先处理法写出图示结构综合结点荷载列阵{}P 。
l /2l/2ql33、考虑弯曲、轴向变形，计算图示结构综合结点荷载列阵{}P 。
m 4m2m2m334、考虑弯曲、轴向变形时，用先处理法计算图示结构综合结点荷载列阵{}P 。
8mm5m635、用先处理法计算图示结构的综合结点荷载列阵{}P 。
/2/2ll36、试 用 矩 阵 位 移 法 解 图 示 结构，绘 弯 矩 图 。
m 1mm1kN m.37、计算下图结构（a ）中杆34的杆端力列阵中的第3个元素和第6个元素。
不计杆件的轴向变形。
已知下图结构（a ）结点位移列阵为：{}[]T 0.66667 0.2 0.7556- 0 0.3667 0 0.3333 0.2 0.2- 0.1333 00.2- 0 0 0=∆。
1m1mll（a ） （b ） 38、计算上图结构（b ）单元③的杆端力列阵{}③F，已知各杆,cm 300 ,kN/cm 101.2424=⨯=I E ,cm 202=A cm l 100=，结点2位移列阵{}[][]T T 2rad cm cm 5313.04596.04730.01012222--⨯⨯==-θ∆v u 。
39、考虑杆件的轴向变形，计算图示结构中单元①的杆端力{}F ①。
已知：I =(/),124m4E =⨯3107kN /m 2， m 2A =05.。
结点1的位移列阵{}[]δ16110370022710151485=⨯⨯---...m m rad T。
5m40、计算图示刚架单元①在局部坐标下的杆端力{}F ①。
已知各杆E 、A 、I 、l 均为常数，不考虑杆件的轴向变形，{}[]∆=--ql EIl l 2100002727 0 5 19 0 0T。
lq41、已求得图示结构结点2、3的结点位移为式(a)、(b)并已知单元②的整体坐标的单元刚度矩阵。
计算单元②2端的弯矩。
（长度单位m ，力单位kN ，角度单位弧度）(b)10 , (a)1040-160-0.2=5-5-ΛΛΛΛ⨯⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧---=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⨯⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧108.1593.0333222φφv u v u[]510205.1105.1050005005.105.15.105.1105.1205.1050005005.105.15.105.1⨯⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-------=②k42、用先处理法写出图示桁架的结构刚度矩阵[]K 。