结构力学PPT备课课件版 13矩阵法1.ppt

1,矩阵位移法,第十三章,2,矩阵代数复习,1、矩阵定义,一组元素按行、列次序排列成的矩形阵列称为矩阵若矩阵的元素排列为m 行和n列，称为mn 阶矩阵2、方阵,3、行矩阵和列矩阵,一个单独的行组成的矩阵称为行矩阵，如：,由单列组成的矩阵称为列矩阵，如：,3,4、纯量,仅由一个单独的元素所组成的11阶矩阵称为纯量5、矩阵乘法,两个规则：,（1）两个矩阵仅当他们是共形时才能相乘，即,（2）不具有交换律，即,4,6、转置矩阵,将一个阶矩阵的行和列依次互换，所得的阶矩阵称之为原矩阵的转置矩阵，如：,其转置矩阵为,当连乘矩阵的乘积被转置时，等于倒转了顺序的各矩阵的转置矩阵之乘积若,A=B C D,7、零矩阵,元素全部为零的矩阵称为零矩阵，用0表示5,任意矩阵与单位矩阵相乘仍等于原矩阵，即,6,10、逆矩阵 在矩阵运算中，没有矩阵的直接除法， 除法运算由矩阵求逆来完成例如，若,一个矩阵的逆矩阵由以下关系式定义：,矩阵求逆时必须满足两个条件：,（1）矩阵是一个方阵2）矩阵的行列式不为零，即矩阵是非奇异矩阵（行列式为零的矩阵称为奇异矩阵）正交矩阵的逆矩阵等于其转置矩阵，即,7,13-1 概 述,矩阵位移法的理论基础是传统的位移法，只是它的表达形式采用矩阵代数，而这种数学算法便于编制计算机程序，实现计算过程的程序化。

一、矩阵位移法的基本思路 矩阵位移法又可以称为杆件结构的有限元法； 矩阵位移法的两个基本步骤是 （1）结构的离散化；（2）单元分析；（3）整体分析，,8,指杆件除有弯曲变形外，还有轴向变形和剪切变形的单元，杆件两端各有三个位移分量，这是平面结构杆件单元的一般情况符号规则：图(a)表示单元编号、杆端编号和局部座标，局部座标的,座标与杆轴重合；,1,2,E A I,l,(a),图(b)表示的杆端位移均为正方向单元编号杆端编号局部座标,1,2,(b),杆端位移编号,1,2,杆端力编号,(c),二、杆端位移、杆端力的正负号规定,一般单元：,9,1,2,1,2,（1）单元杆端位移向量,（2）单元杆端力向量,凡是符号上面带了一横杠的就表示是基于局部座标系而言的10,现在讨论单元刚度方程单元刚度方程是指由单元杆端位移求单元杆端力时的一组方程，可以用“ ”表示，由位移求力称为正问题在单元两端加上人为控制的附加约束，使基本杆单元的两端产生任意指定的六个位移，然后根据这六个杆端位移来推导相应的六个杆端力1,2,我们忽略轴向受力状态和弯曲受力状态之间的相互影响，分别推导轴向变形和弯曲变形的刚度方程13-2 单元刚度矩阵(局部座标系),进行单元分析，推导单元刚度方程和单元刚度矩阵。

一、一般单元,11,分别推导轴向变形和弯曲变形的刚度方程首先，由两个杆端轴向位移,可推算出相应的杆端轴向力,1,2,其次，由杆端横向位移,可以用角变位移方程推导出相应的杆端,12,将上面六个方程合并，写成矩阵形式：,13,EA l,6EI l2,6EI l2,EA l,12EI l3,12EI l3,4EI l,2EI l,上面的式子可以用矩阵符号记为,这就是局部座标系中的单元刚度方程可求单元杆端力,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,只与杆件本身性质有关而与外荷载无关,通过这个式子由单元杆端位移,局部座标系的单元刚度矩阵,14,二、单元刚度矩阵的性质,（1）单元刚度系数的意义,代表单元杆端第j个位移分量等于1时所引起的第i个杆端力分量例如,代表单元杆端第2个位移分量 时所引起的第5个杆端力分量 的数值2）单元刚度矩阵 是对称矩阵，,（3）一般单元的刚度矩阵 是奇异矩阵；,因此它的逆矩阵不存在,从力学上的理解是，根据单元刚度方程,由,有一组力的解答(唯一的)，即正问题。

由,如果,不是一组平衡力系则无解；若是一组平衡力系，则解答不是唯一的，即反问题15,三、特殊单元 若单元六个杆端位移中有某一个或几个已知为零，则该单元称为特殊单元，其刚度方程是一般单元刚度方程的特例以连续梁为例：,1,2,16,为了程序的标准化和通用性，不采用特殊单元，只用一般单元，如果结构有特殊单元，可以通过程序由一般单元来形成17,13-3 单元刚度矩阵(整体座标系),座标转换矩阵,单元杆端力的转换式、单刚的转换式,一、单元座标转换矩阵,18,正交矩阵,T-1 =TT,或 TTT=TT T =I,于是可以有,同理可以有,19,在局部座标系中杆端力与杆端位移的关系式表达为：,在整体座标系中杆端力与杆端位移的关系式可以表达为：,二、整体座标系中的单元刚度矩阵,（a）式可转换为：,两边前乘TT,比较式(b)和(d)可得：,20,例1. 试求图示刚架中各单元在整体座标系中的刚度矩阵k 设 和 杆的杆长和截面尺寸相同l=5m，bh=0.5m 1m,解:,(1) 局部座标系中的单元刚度矩阵,(2) 整体座标系中的单元刚度矩阵,单元 1 ： = 0，T =I,单元 2 ： = 90，单元 座标转换矩阵为,21,单元 2 ： = 90，单元座标转换矩阵为,22,13-4 连续梁的整体刚度矩阵,按传统的位移法,1,2,3,每个结点位移对F的单独贡献,F1,F2,F3,4i1,2i1,0,2i1,4i1+4i2,2i2,0,2i2,4i2,1,2,3,=,F=K,根据每个结点位移对附加约束上的约束力F的贡献大小进行叠加而计算所得。

传统位移法,23,一、 单元集成法的力学模型和基本概念,分别考虑每个单元对F的单独贡献，整体刚度矩阵由单元直接集成,令 i2 =0，则,k =,4i1,2i1,4i1,2i1,F2,=,1,2,（a）,（b）,F2,=,4i1,2i1,4i1,2i1,0,0,0,0,0,1,2,3,单元 1 的贡献矩阵,单元 1 对结点力F的贡献,略去其它单元的贡献24,k =,4i2,2i2,4i2,2i2,F2,=,4i1,2i1,4i1,2i1,0,0,0,0,0,1,2,3,K ,F =,设 i1 =0，则,单元 的贡献矩阵,单元对结点力F的贡献,略去单元的贡献25,k,K,K,F=F,+F,=（K,+K,）,1,2,F=K,整体刚度矩阵为：,单元集成法求整体刚度矩阵步骤：,根据单元和单元分别对结点力F的贡献，可得整体刚度方程：,26,2i2,2i2,4i2,K=,4i1,2i1,4(i1+i2),2i1,0,2i2,0,2i2,4i2,4i1+4i2,27,二、按照单元定位向量由,(1)在整体分析中按结构的结点位移统一编码，称为总码2)在单元分析中按单元两端结点位移单独编码，称为局部码以连续梁为例,位移统一编码，总码,单元,对应关系,局部码总码,单元定位向量,(1)1,(2)2,=,(1)2,(2)3,=,位移单独编码局部码,由单元的结点位移总码组成的向量,28,单元,单元,(1),(2),(1),(2),4i1,2i1,2i1,4i1,k =,(1),(2),(1),(2),4i2,2i2,4i2,2i2,单元定位向量描述了单元两种编码（总码、局部码）之间的对应关系。

单元定位向量定义了整体坐标系下的单元刚度矩阵中的元素在整体刚度矩阵中的具体位置，故也称为“单元换码向量”单元贡献矩阵是单元刚度矩阵，利用“单元定位向量”进行“换码重排位”29,三、 单元集成法的实施,（定位 累加）,K,（1）将K置零，得K=0；,（2）将k的元素在K中按定位并进行累加，得K=K；,（3）将k的元素在K中按定位并进行累加，得K=K+K；,按此作法对所有单元循环一遍，最后即得整体刚度矩阵K30,1,2,3,0= 0,（1）结点位移分量总码,（2）单元定位向量,=,=,=,（3）单元集成过程,1,2,2,1,k =,2,3,3,2,k =,0,3,3,0,K =,4i1,2i1,2i1,2i2,2i2,4i2,4i1,4i2+4i3,4i1+4i2,例.求连续梁的整 体刚度矩阵31,四、整体刚度矩阵 K 的性质,（1）整体刚度系数的意义: Kijj=1 (其余=0)时产生的结点力Fi,（2）K是对称矩阵,（3）对几何不变体系，K是可逆矩阵，如连续梁,F=K,=K-1F,（4）K是稀疏矩阵和带状矩阵，如连续梁,1,2,3,n,n+1,4i1,2i1,2i1,2i2,2i2,4i2+4i3,4i1+4i2,4in,2i3,2in,。