[工学]9矩阵位移法.ppt

§9-1 概述 §9-2 单元刚度矩阵（局部坐标系） (element stiffness matrix) §9-3 单元刚度矩阵（整体坐标系） §9-4 连续梁的整体分析assembly analysis of continuous beam §9-5 刚架的整体刚度矩阵(assembly analysis of frame) §9-10 连续梁程序的框图设计 求图示刚架的内力。设各杆为矩形截面，横梁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 按传统的位移法 i1 i2 1 2 ?1 4i1?1 2i1?1 0 i1 i2 1 2 ?2 2i1?2 2i2?2 (4i1+4i2)?2 i1 i2 1 2 ?3 0 2i2?3 4i2?3 每个结点位 移对{F}的单 独贡献 F1 F2 F3 4i1 2i1 0 2i1 4i1+4i2 2i2 0 2i2 4i2 ?1 ?2 ?3 = {F}=[K]{?} 根据每个结点位移对附加约束上的约束力{F}的贡献大小进行叠加而计算所得。 传统位移法 一、 单元集成法的力学模型和基本概念 分别考虑每个单元对{F}的单独贡献，整体刚度矩阵由单元直接集成 i1 i2 1 2 ?1 ?2 ?3 F3 {F} 1 = [ F1 1 F2 1 1 ]T F1 1 F2 1 F3 1 令 i2 =0，则 F3 1 =0 [k] = 4i1 2i1 4i1 2i1 1 F1 1 F2 1 = 4i1 2i1 4i1 2i1 ?1 ?2 （a） （b） 单元 1 对结点力{F}的贡献 略去其它单元的贡献。 F1 1 F2 1 F3 1 = 4i1 2i1 4i1 2i1 0 0 0 0 0 ?1 ?2 ?3 1 [K] {?} {F} = 1 [K] = 1 4i1 2i1 4i1 2i1 0 0 0 0 0 单元 1 的贡献矩阵 i1 i2 1 2 ?1 ?2 ?3 F1 2 F2 2 F3 2 [k] = 4i2 2i2 4i2 2i2 2 F1 2 F2 2 F3 2 = 4i1 2i1 4i1 2i1 0 0 0 0 0 ?1 ?2 ?3 2 [K] {?} {F} = 2 设 i1 =0，则 F1 2 =0 [K] = 2 4i1 2i1 4i1 2i1 0 0 0 0 0 单元 ? 的贡献矩阵 F3 {F} 2 = [ F1 2 F2 2 2 ]T 单元?对结点力{F}的贡献 略去单元?的贡献。 1 [K] {?} {F} = 1 [K] = 1 4i1 2i1 4i1 2i1 0 0 0 0 0 2 [K] {?} {F} = 2 [K] = 2 4i1 2i1 4i1 2i1 0 0 0 0 0 i1 i2 1 2 1 2 1 2 [K]=（[K] +[K] ）= 1 2 e e [k] [K] [K] e e {F}={F} +{F} =（[K] +[K] ）{?} 1 2 {F}=[K]{?} 整体刚度矩阵为： 单元集成法求整体刚度矩阵步骤： 根据单元?和单元?分别对结点力{F}的贡献，可得整体刚度方程： [k] [K] [K] e e 1 2 [k] = 4i1 2i1 4i1 2i1 1 [K] = 1 4i1 2i1 4i1 2i1 0 0 0 0 0 [k] = 4i2 2i2 4i2 2i2 2 [K] = 2 4i2 2i2 4i2 2i2 0 0 0 0 0 1 2 1 4i1 2i1 4i1 2i1 0 0 0 0 0 2i2 2i2 4i2 [K]= 4i1 2i1 4(i1+i2) 2i1 0 2i2 0 2i2 4i2 4i1+4i2 二、按照单元定位向量由 [k] 求 e [K] e (1)在