第11章-矩阵位移法.pptVIP

单元② 单元③ （3）利用单元定位向量 集成结构刚度矩阵K （4）计算综合结点荷载列阵P 已在例11.2中求得，为 （5）解结构刚度方程KD=P，求出结点位移列阵D 得结点位移D为 （6）计算单元杆端内力 按式（11.29），即 计算各单元杆端内力，分别为 单元①： a = 0° 单元定位向量 单元②： a = 0° 单元定位向量 单元②： a = 90° 单元定位向量 （7）绘制内力图 ③求综合结点荷载列阵P，即利用公式（11.28） 求出P。 5)解总刚方程 解结构刚度方程KD=P，求出结点位移列阵 D 6)计算单元杆端内力 如果计算单刚时，先计算了单元坐标系中的单刚 ，可用公式（11.29）求单元杆端内力，即 如果是直接使用式（11.24）和式（11.25）求出整体坐标系中的单刚 则用公式（11.32）来求单元杆端内力，即 例. 求图示刚架的内力。设各杆为矩形截面，横梁b2×h2=0.5m ×1.26m， 立柱b1 ×h1=0.5m ×1m。 （1）原始数据、局部码、总码（设E=1） 12m 6m A B C D q=1kN/m A B C D 1 2 3 x y 1 3 4 5 2 6 {0} {0} 柱 梁 e e （2）形成局部座标系中的单元刚度矩阵 k e 单元1和3 =10-3× =10-3× 单元2 （3）计算整体座标系中的单元刚度矩阵 e [k] [k] = [T]T k e [T] e 2 [k] – 1 [k] = – 3 [k] – 单元1和3 的座标转换 矩阵 （?=900） 1 [k] = =10-3× [k] = [T]T k 1 [T] 3 单元2 （?=0°） 2 [k] k 2 = =10-3× (4)用单元集成法形成整体刚度矩阵[K] A B C D 1 2 3 x y 1 3 4 5 2 6 {0} {0} 2 1 3 （5）求等效结点荷载{P} 12m 6m A B C D q=1kN/m A B C D 1 2 3 x y 1 3 4 5 2 6 {0} {0} 1 单元固端约束力 单元1 （?=90°） 1 1 1 按单元定位向量 1 （6）解基本方程 求得结点位移: （7）求各单元杆端力 e e e e e 1 1 1 1 1 1 单元1：先求 {F} 然后求 ① ① 1 1 1 1 =10-3× 1 1 2 3 同样可得出： （8）绘制内力图 1 2 3 1 2 e A B C D 8.49 2.09 3.04 4.38 M图(kN·m) 4.76 1.24 0.43 1.24 Q图(kN) N图(kN) 0.43 0.43 1.24 1）结点位移分量的统一编码——总码 在刚结点A铰结点C1和C2处， 竖向位移均为零，故其编码 也应为零，另外它们的水平 位移分量都相等，因此它们 的水平位移应采用同码。 0 0 0 1 0 2 2 1 A C1 B D 0 0 0 1 0 3 1 4 0 C2 3 2）单元定位向量： {λ} = 1 [1 0 2 1 0 3]T {λ} = 2 [1 0 2 0 0 0]T {λ} = 3 [1 0 4 0 0 0]T 3）按①②③次序进行单元集成： 忽略轴向变形时矩形刚架的整体分析 ×104 [k] = 1 1 0 2 1 0 3 1 2 3 4 1 2 3 4 1 0 2 1 0 3 300 0 －300 0 100 0 +0 50 －300 +0 +300 +0 0 50 +0 100 1 0 2 0 0 0 k 2 = ×104 1 0 2 0 0 0 0 0 0 0 100