矩阵位移法习题课.pptVIP

* 矩阵位移法习题课 单元分析任务：建立单元两端位移与两端力之间的解析关系。 目标：生成局部坐标系下的单元刚度矩阵。 整体分析任务：建立结构体系结点位移与结点力之间的解析关系。 步骤：1. 通过坐标转换矩阵，形成整体坐标系下的单元刚度矩阵； 2. 根据单元定位向量，集成整体刚度矩阵； 3. 建立整体坐标系下的结构体系的静力平衡方程。 目标：计算外部作用下结构体系各结点位移值，进一步可计算各 杆端力，然后根据各杆件所受荷载特点计算整个结构体系 各截面的内力。 EA l 6EI l2 6EI l2 EA l 12EI l3 12EI l3 4EI l 2EI l 局部坐标系下的单元刚度矩阵 e 矩阵元素物理意义：单元杆端单位位移下所产生的杆端力。 e = (1) (2) (3) (4) (5) (6) (1) (2) (3) (4) (5) (6) 0 0 0 0 0 0 6EI l2 0 6EI l2 0 -EA l -6EI l2 6EI l2 EA l -12EI l3 12EI l3 2EI l 4EI l 0 0 0 0 0 0 -6EI l2 0 6EI l2 0 只与杆件本身性质有关而与外荷载无关 ? e x y X1 Y1 X2 Y2 e e e e e 座标转换矩阵 单元座标转换矩阵 目的：将局部坐标系中的杆端力 和杆端位移统一在整体坐 标系中。 正交矩阵 [T]-1 =[T]T 或 [T][T]T=[T]T [T] =[I] 存在： e e e e 在局部座标系中杆端力与杆端位移的关系式： 在整体座标系中杆端力与杆端位移的关系式： 整体坐标系下的单元刚度矩阵 单元定位向量 由单元结点位移总码组成的向量。 是单元刚度矩阵各元素在整体刚度矩阵中的位置指针向量。 注意：对结点编号时，应尽量使每一个单元的两端结点编号差值最小。 目的是使结构体系的整体刚度矩阵的带宽达到最小。 整体刚度矩阵 [K] 的性质 (1)整体刚度系数的物理意义： Kij：?j=1 (其余?=0)时产生的结点力Fi； (2)[K]是对称矩阵； (3)对几何不变体系，[K]是可逆矩阵，即存在： {F}=[K]{?} {?}=[K]-1{F} (4)[K]是稀疏矩阵和带状矩阵，如连续梁： ?1 ?2 ?3 F1 F2 F3 1 2 3 n ?n Fn ?n+1 Fn+1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4i1 2i1 2i1 2i2 2i2 4i2+4i3 4i1+4i2 4in 2i3 2in 例1 如图所示连续梁结构，试求在图示荷载作用下的各单元杆端力和支座反力。 已知 1. 未知量编号 连续梁可动结点自左向右连续编号，被约束的位移编号为0。 2. 单元刚度矩阵 对单元1、单元2有： (1) (1) (2) (2) (0) (0) (1) (1) x y 1 2 3 1 2 4 3 A B C D 2. 单元定位向量 单元1： (2) (3) (4) (2) (3) (4) 单元2： 单元3： 3. 集成总体刚度矩阵 （1） （1） （2） （2） （3） （3） （4） （4） 4. 等效荷载向量 单元2： 单元3： 荷载向量： 5. 解方程： 求结点位移 例2 试求图示门式刚架的总体刚度矩阵。已知各杆件截面特征如图所示，且不 考虑轴向变形的影响。 1. 体系基本未知量分析 2. 单元刚度矩阵 整体坐标系下的单刚矩阵 单元1： 单刚 元素 简记 (1) (1) (2) (2) (0) (0) (0) (0) 0 0 0 0 0 0 0 0 (3) (3) (1) (1) (0) (0) (0) (0) 单元2： 1 2 3 3 2 1 x 1