[理学]10矩阵位移法.pptVIP

（3）结点力矩 平衡条件 （4）将杆端弯矩代入上式，并用矩阵表示 ——位移法方程（平衡方程）。即为结构（整体）刚度方程。可简记为： 其中，k——为结构（整体）刚度矩阵。 F2 1 F1 M12 2 M21 M23 F3 M32 2 §10-4 连续梁的整体刚度矩阵 两单元贡献矩阵集成 2、单元集成法（直接刚度法）步骤 1）单刚“换码” 单元①：k①= ①= λ① 1 2 1 2 §10-4 连续梁的整体刚度矩阵 —单元杆端位移编码与整体结点位移编码关系(若结点位移为零，其总码编为零), 即可得到单元定位向量 ，并分别标在 的上方和左（右）边。 1 3 2 F1 F2 F3 ② ① i1 i2 Δ1 Δ2 Δ3 2）总刚“集成” 1 2 3 1 2 3 单元②：k②= ②= λ② 2 3 2 3 §10-4 连续梁的整体刚度矩阵 ——按单刚各元素的行码和列码（单元定位向量）分别将其置于结构刚度矩阵相应的位置，即“对号入座”，集成总刚。 1 3 2 F1 F2 F3 ② ① i1 i2 Δ1 Δ2 Δ3 例题：试求图10-10所示连续梁的整体刚度矩阵k。 解：1）结点编码、单元划分。 结点位移编码—— 均示于图10-10中（凡 是给定为零值的结点 位移分量，其总码均编为零。） 2）列出各单元定位向量λe λ①=[1 2]T λ②=[2 3]T λ③=[3 0]T §10-4 连续梁的整体刚度矩阵 图10-10 i1 i2 i3 1 2 0 3 ② ① ③ Δ1 Δ2 Δ3 Δ0=0 4）“对号入座”， 集成总刚。 k①= 1 2 1 2 §10-4 连续梁的整体刚度矩阵 i1 i2 i3 1 2 0 3 ② ① ③ Δ1 Δ2 Δ3 Δ0=0 k②= 2 3 2 3 k ③= 3 0 3 0 1 2 3 1 2 3 k = 本节讨论用单元集成法求图10-11所示平面刚架的整体刚 度矩阵k。原理、方法同连续梁，但具体求解过程要复杂些。 1、求解方法与步骤 1）编码 （1）单元：①、②； （2）结点：1、2、3； （3）结点位移——位移编码顺序为： 即： 。 图10-11 §10-5 刚架的整体刚度矩阵 ② ① (1、2、3) (0、0、4) (0、0、0) 1 2 3 ——对已知为 零的结点位移 分量，其总码 为零。 2）选取结构坐标系和单元坐标系（确定单元的始端和末端）如左图 单元坐标系 x y 结构坐标系 ② ① (1、2、3) (0、0、4) (0、0、0) 1 2 3 ② ① (1、2、3) (0、0、4) (0、0、0) 1 2 3 §10-5 刚架的整体刚度矩阵 3）各单元定位向量λe λ①=(1 2 3 0 0 4)T λ②=(1 2 3 0 0 0)T 4）单元刚度矩阵 ∴单元①： ② ① (1、2、3) (0、0、4) (0、0、0) 1 2 3 单元坐标系 x y 结构坐标系 k①= ①＝104× 1 2 3 0 0 4 1 2 3 0 0 4 图10-11 §10-5 刚架的整体刚度矩阵 单元②： ② ① = T= k②=TT T = ② 1 2 3 0 0 0 1 2 3 0 0 0 §10-5 刚架的整体刚度矩阵 5）对号入座，集成总刚 ——即为所求平面刚架的整体刚度矩阵。 1 2 3 4 1 2 3 4 §10-5 刚架的整体刚度矩阵 2、整体刚度矩阵特点 1）整体刚度系数（k i j）的意义 ——表示当第j个结点位移分量Δ1=1（其它结点位移分量为零）时所产生的第i个结点力F i ； 2）k是对称矩阵（反力互等定理）； 3）k是满秩非奇异矩阵（先处理法，已考虑约束条件）； 4）k是稀疏、带状矩阵（即:k有许多零元素，且非零元素都分布在以主对角线为中心的倾斜带状区城内）。 §10-5 刚架的整体刚度矩阵 3、先处理法处理有中间铰刚架 图10-12所示为具有铰结点的刚架。现说明具有铰结点刚架的