结构力学10矩阵法2.ppt

1、1,一、矩阵位移法的基本思路 矩阵位移法的两个基本步骤是 （1）结构的离散化；（2）单元分析；（3）整体分析，,2,指杆件除有弯曲变形外，还有轴向变形和剪切变形的单元，杆件两端各有三个位移分量，,符号规则：图(a)表示单元编号、杆端编号和局部座标，局部座标的,座标与杆轴重合；,1,2,E A I,l,(a),图(b)表示的杆端位移均为正方向。,单元编号 杆端编号 局部座标,1,2,(b),杆端位移编号,1,2,杆端力编号,(c),二、杆端位移、杆端力的正负号规定,一般单元：,3,局部座标系中的单元刚度方程,只与杆件本身性质有关而与外荷载无关,局部座标系的单元刚度矩阵,4,10-3 单元刚度矩阵

2、(整体座标系),座标转换矩阵,一、单元座标转换矩阵,正交矩阵,T-1 =TT,5,三、单元刚度矩阵的性质,二、整体座标系中的单元刚度矩阵,6,10-4 连续梁的整体刚度矩阵,按传统的位移法,1,2,3,每个结点位 移对F的单 独贡献,F1,F2,F3,4i1,2i1,0,2i1,4i1+4i2,2i2,0,2i2,4i2,1,2,3,=,F=K,根据每个结点位移对附加约束上的约束力F的贡献大小进行叠加而计算所得。,传统位移法,7,一、 单元集成法的力学模型和基本概念,分别考虑每个单元对F的单独贡献，整体刚度矩阵由单元直接集成,令 i2 =0，则,k =,4i1,2i1,4i1,2i1,F2,=

3、,1,2,（a）,（b）,F2,=,4i1,2i1,4i1,2i1,0,0,0,0,0,1,2,3,单元 1 的贡献矩阵,单元 1 对结点力F的贡献,略去其它单元的贡献。,8,k =,4i2,2i2,4i2,2i2,F2,=,4i1,2i1,4i1,2i1,0,0,0,0,0,1,2,3,K ,F =,设 i1 =0，则,单元 的贡献矩阵,单元对结点力F的贡献,略去单元的贡献。,9,k,K,K,F=F,+F,=（K,+K,）,1,2,F=K,整体刚度矩阵为：,单元集成法求整体刚度矩阵步骤：,根据单元和单元分别对结点力F的贡献，可得整体刚度方程：,10,2i2,2i2,4i2,K=,4i1,2i

4、1,4(i1+i2),2i1,0,2i2,0,2i2,4i2,4i1+4i2,整体刚度矩阵:,11,二、按照单元定位向量由,(1)在整体分析中按结构的结点位移统一编码，称为总码。,(2)在单元分析中按单元两端结点位移单独编码，称为局部码。,以连续梁为例,位移统一编码，总码,单元,对应关系,局部码总码,单元定位向量,(1)1,(2)2,=,(1)2,(2)3,=,位移单独编码 局部码,由单元的结点 位移总码组成 的向量,12,单元,单元,(1),(2),(1),(2),4i1,2i1,2i1,4i1,k =,(1),(2),(1),(2),4i2,2i2,4i2,2i2,单元定位向量描述了单元两

5、种编码（总码、局部码）之间的对应关系。,单元定位向量定义了整体坐标系下的单元刚度矩阵中的元素在整体刚度矩阵中的具体位置，故也称为“单元换码向量”。,单元贡献矩阵是单元刚度矩阵，利用“单元定位向量”进行“换码重排位”。,13,三、 单元集成法的实施,（定位 累加）,K,（1）将K置零，得K=0；,（2）将k的元素在K中按定位并进行累加，得K=K；,（3）将k的元素在K中按定位并进行累加，得K=K+K；,按此作法对所有单元循环一遍，最后即得整体刚度矩阵K。,14,1,2,3,0= 0,（1）结点位移分量总码,（2）单元定位向量,=,=,=,（3）单元集成过程,1,2,2,1,k =,2,3,3,2

6、,k =,0,3,3,0,K =,4i1,2i1,2i1,2i2,2i2,4i2,4i1,4i2+4i3,4i1+4i2,例.求连续梁的整 体刚度矩阵。,15,四、整体刚度矩阵 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,16,10-5 刚架的整体刚度矩阵,思路要点：（1）设各单元已形成了整体座标系下的单元刚度

7、矩阵；,与连续梁相比： (1)各单元考虑轴向变形；(2)每个刚结点有三个位移； (3)要采用整体座标；(4)要处理非刚结点的特殊情况。,一、结点位移分量的统一编码总码, =1 2 3 4 T,规定： 对于已知为零的结点位移分量，其总码均编为零。,=uA vA A C T,整体结构的结点位移向量为：,相应地结点力向量为：,= XA YA MA MC T,F = F1 F2 F3 F4 T,17,二、单元定位向量,单元,单元,局部码总码,局部码总码,（1） 1 （2） 2 （3） 3 （4） 0 （5） 0 （6） 4,（1） 1 （2） 2 （3） 3 （4） 0 （5） 0 （6） 0,三、单元集成过程,18,K=,11,12,13,21,22,23,31,32,33,61,62,63,66,16,26,36,11,12,13,21,22,23,31,32,33,19,四、铰结点的处理,K,求单元常数,T,单元刚度矩阵,程序设计框图（局部：集成整体刚度矩阵）,刚结点：变形连