矩阵法结构力学.ppt

图示四单元的l,EI,EA 相同，它们n局部坐标系下的单元刚度矩阵各不相同。 一般单元的单元刚度矩阵一定是奇异矩阵，而特殊单元的单元刚度矩阵一定是非奇异矩阵。 如特殊单元是几何不变体系，其单元刚度矩阵一定是非奇异矩阵。 　 由一般单元的单元刚度方程，任给一组满足平衡的杆端力，可唯一地确定杆端位移。 由一般单元的单元刚度方程，任给一组可能的杆端位移，可唯一地确定杆端力。 l l P EA= c 1 2 3 4 5 A B C D x y 1 2 3 4 5 6 7 8 例：试用后处理法计算图示桁架各杆内力，设各杆EA为常数。 （1）单元和结点编码，准备基本数据。 （2）建立结点位移向量和结点力向量： {?}={?1 ?2 ?3 ?4 ?5 ?6 ?7 ?8}T {P}={F1 F2 F3 F4 0 P 0 0}T （3）建立整体座标系单刚 e [k] 对称 对称 1 2 5 6 5 6 7 8 对称 5 6 7 8 1 2 5 6 3 4 7 8 3 4 7 8 1 2 3 4 5 A B C D 1 2 3 4 5 6 7 8 （3）建立整体座标系单刚 e [k] 对称 1 2 5 6 5 6 7 8 对称 5 6 7 8 1 2 5 6 对称 3 4 7 8 3 4 7 8 对称 3 4 5 6 3 4 5 6 1 2 7 8 1 2 7 8 对称 4.形成原始总矩阵位移法方程 5.引进支座位边界条件 ?1 = ?2= ?3 = ?4=0 划去1 → 4行(列) 对称 1 2 3 4 5 6 7 8 1 2 3 4 5 6 7 8 6.解矩阵位移法方程 对称 2 组合结构 ★区分梁式杆和桁架杆. ★对梁式杆采用梁单元刚度方程和相应的计算公式; 对桁架杆采用桁架单元刚度方程和相应的计算公式. 例题 试求图示组合结构的内力。已知横梁的EA=2EI。 吊杆的E1A1=EI/20。 10kN/m 15m 20m 20m 20m ① ② ③ ④ ⑤ ① ② ③ ④ ⑤ 1 2 3 4 5 6 1 确定坐标系、单元和结点位移编码 解： ① ③ ② ⑤ ④ 2 形成刚度矩阵 (1)形成局部坐标系下的单元刚度矩阵 1、结构刚度方程矩阵形式为： ，它是整个结构所应满足的变形条件。 2、在直接刚度法中，定位向量的物理意义是变形连续条件和位移边界条件。 3、等效结点荷载数值等于汇交于该结点所有固端力的代数和。 4、矩阵位移法中，等效结点荷载的“等效原则”是指与非结点荷载的结点位移相等。 5、矩阵位移法既能计算超静定结构，也能计算静定结构。 k3 k2 k1 §9-8 忽略轴向变形时矩形刚架的整体分析 1 2 3 A B D x y 0 0 0 1 0 2 1 0 3 C1 C2 1 0 4 0 0 0 1 2 3 单元定位向量 1 2 3 4 5 6 1 2 3 4 5 6 1 0 2 1 0 3 1 0 2 1 0 3 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 0 2 0 0 0 1 0 2 0 0 0 1 0 4 0 0 0 1 0 4 0 0 0 1 2 3 4 [K]= 1 2 3 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 k3 k2 k1 1 2 3 4 5 6 1 2 3 4 5 6 1 0 2 1 0 3 1 0 2 1 0 3 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 0 2 0 0 0 1 0 2 0 0 0 1 0 4 0 0 0 1 0 4 0 0 0 例题 做图示刚架的内力图 12m 1kN/m 6m A B C D ① ② 2 1 A B C D ③ 1 3 解 1 确定整体和局部坐标系、单元和结点位移编码 2 形成刚度矩阵 (1)形成局部坐标系下的单元刚度矩阵。 ① ③ ② ① ③ (2)形成整体坐标系下的单元刚度矩阵。 ① ③ ① ① ① ② ② (3) “换码重排座”，形成整体刚度矩阵 集成单元① ① ② 2 1 A B C D ③ 1 3 ① 集成单元② ① ② 2 1 A B C D ③ 1 3 ② 集成单元③ ① ② 2 1 A B C D ③ 1 3 ③ 整体刚度矩阵 3 形成等效结点荷载 (1)形成局部坐标系下的单元固端约束力 ① (3) “换码重排座”，形成整体坐标系下的等效结点荷载 (2)形成整体坐标系下的单元等效结点荷载 ① ① ① 4 解基本方程，求结点位移 5 求各单元的杆端内力 (1)由结点位移，确