13矩阵位移法1.pptVIP

9-1 结构的离散化与杆端位移、杆端力的符号 9-3 整体坐标下的单元刚度矩阵 9-4 连续梁的整体刚度矩阵 1. 整体刚度矩阵的集成 2. 整体刚度矩阵的性质 整体刚度矩阵置零 集成单元①的刚度矩阵 集成单元②的刚度矩阵 集成单元③的刚度矩阵 (3)换码重排座 ① ② ③ * * 第九章 矩阵位移法 矩阵代数复习 1、矩阵定义 一组元素按行、列次序排列成的矩形阵列称为矩阵。若矩阵的元素排列为m 行和n列，称为m?n 阶矩阵。 A = a a a a a a a a a n n m m mn 11 12 1 21 22 2 1 2 L L M O M L é ? ê ê ê ê ê ù ? ú ú ú ú ú 2、方阵 一个具有相同的行数和列数的矩阵，即 m = n 时，称为 n 阶方阵。 3、行矩阵和列矩阵 一个单独的行组成的矩阵称为行矩阵，如： A = [ ] a a a a n 11 12 13 1 ? ? ? 由单列组成的矩阵称为列矩阵，如： 4、纯量 仅由一个单独的元素所组成的1?1阶矩阵称为纯量。 5、矩阵乘法 两个规则： （1）两个矩阵仅当他们是共形时才能相乘，即 A B C p l m p l n m n ′ ′ ′ = = 当 时才能相乘 A B = a a a a b b 11 12 21 22 11 21 é ? ù ? ú é ? ù ? ú 共形 2 × 2 2 × 1 B A= b b a a a a 11 21 11 12 21 22 é ? ù ? ú é ? ù ? ú 非 共形 2 × 1 2 × 2 （2）不具有交换律，即 AB 1 BA 6、转置矩阵 将一个阶矩阵的行和列依次互换，所得的阶矩阵称之为原矩阵的转置矩阵，如： A= a a a a a a 11 12 21 22 31 32 é ? ê ê ê ù ? ú ú ú 其转置矩阵为 A T = é ? ê ù ? ú a a a a a a 11 21 31 12 22 32 当连乘矩阵的乘积被转置时，等于倒转了顺序的各矩阵的转置矩阵之乘积。若 A=B C D 则 A T =D T C T B T 7、零矩阵 元素全部为零的矩阵称为零矩阵，用0表示。 若 AB=0 ， 但不一定 A=0 或 B=0。 任意矩阵与单位矩阵相乘仍等于原矩阵，即 AI = A IA = A 10、逆矩阵 在矩阵运算中，没有矩阵的直接除法， 除法运算由矩阵求逆来完成。例如，若 AB = C 则 B = A - 1 C 此处 A-1 称为矩阵 A 的逆矩阵。 一个矩阵的逆矩阵由以下关系式定义： A A - 1 = A - 1 A = I 矩阵求逆时必须满足两个条件： （1）矩阵是一个方阵。 （2）矩阵的行列式不为零，即矩阵是非奇异矩阵（行列式为零的矩阵称为奇异矩阵）。 11、正交矩阵 若一方阵A 每一行（列）的各个元素平方之和等于1，而 所有的两个不同行（列）的对应元素乘积之和均为零，则称该矩阵为正交矩阵，则 A = cos sin sin cos a a a a - é ? ê ù ? ú 正交矩阵的逆矩阵等于其转置矩阵，即 A - 1 = A T 结构的离散化：等截面直杆单元 ① ② ③ ④ ⑤ A B C E D F ① ② ③ ④ A B D E C 杆端力和杆端位移的符号 ■弯矩、转角：绕杆端顺时针为正； ■其它：与坐标轴同向为正。 i E，I，A，l j ■ ：顺时针为正 e e 杆端位移 杆端力 局部坐标系 e 9-2 局部坐标系中自由单元的单元刚度矩阵 … … … … e e e 1.一般单元的刚度方程和刚度矩阵 用 e 表示 e 用 e 表示 e e e 局部坐标下的单元刚度方程 e e e e 局部坐标下的自由单元的单元刚度矩阵 2. 单元刚度矩阵的性质 （1）单元刚度系数的意义 单位杆端位移引起的杆端力 （2）单元刚度矩阵是对称矩阵 反力互等定理 （3）自由单元刚度矩阵是奇异矩阵 矩阵行列式等于零，逆阵不存在。 解不唯一 ★由杆端力只能求出变形，不能求杆端总的位移 （刚体位移+变形）。 解唯一 e e e e e e 3. 特殊单元 若单元六个杆端位移中有某一个或几个已知为零，则该单元称为特殊单元，其刚度方程是一般单元刚度方程的特例。 e 1 2 e （1）连续梁单元的刚度方程 单元两端只有转角位移 e e e e e 1 2 e e e e （2）桁架单元刚度方程 （3）刚架中忽略轴向变形的梁单元刚度方程 1. 单元坐标转换矩阵 局部坐标系下的杆端力 整体坐标系下的杆端力 F(1) x y x y α e i j F(