13 矩阵位移法

2020/6/14,1,13-1概述,13-2单元刚度矩阵(局部座标),第十三章矩阵位移法,进入,进入,13-3单元刚度矩阵(整体座标系),13-4连续梁的整体刚度矩阵,13-5刚架的整体刚度矩阵,13-7计算步骤和算例,进入,进入,进入,进入,13-6等效结点荷载,进入,13-8桁架及组合结构的整体分析,13-8忽略轴向变形时矩形刚架的整体分析,进入,进入,2020/6/14,2,矩阵代数复习,1、矩阵定义,一组元素按行、列次序排列成的矩形阵列称为矩阵。若矩阵的元素排列为m行和n列，称为mn阶矩阵。,2、方阵,3、行矩阵和列矩阵,一个单独的行组成的矩阵称为行矩阵，如：,由单列组成的矩阵称为列矩阵，如：,2020/6/14,3,4、纯量,仅由一个单独的元素所组成的11阶矩阵称为纯量。,5、矩阵乘法,两个规则：,（1）两个矩阵仅当他们是共形时才能相乘，即,（2）不具有交换律，即,2020/6/14,4,6、转置矩阵,将一个阶矩阵的行和列依次互换，所得的阶矩阵称之为原矩阵的转置矩阵，如：,其转置矩阵为,当连乘矩阵的乘积被转置时，等于倒转了顺序的各矩阵的转置矩阵之乘积。若,A=BCD,7、零矩阵,元素全部为零的矩阵称为零矩阵，用0表示。,2020/6/14,5,任意矩阵与单位矩阵相乘仍等于原矩阵，即,2020/6/14,6,10、逆矩阵在矩阵运算中，没有矩阵的直接除法，除法运算由矩阵求逆来完成。例如，若,一个矩阵的逆矩阵由以下关系式定义：,矩阵求逆时必须满足两个条件：,（1）矩阵是一个方阵。,（2）矩阵的行列式不为零，即矩阵是非奇异矩阵（行列式为零的矩阵称为奇异矩阵）。,正交矩阵的逆矩阵等于其转置矩阵，即,2020/6/14,7,13-1概述,矩阵位移法的理论基础是传统的位移法，只是它的表达形式采用矩阵代数，而这种数学算法便于编制计算机程序，实现计算过程的程序化。一、矩阵位移法的基本思路矩阵位移法又可以称为杆件结构的有限元法；矩阵位移法的两个基本步骤是（1）结构的离散化；（2）单元分析；（3）整体分析，,2020/6/14,8,指杆件除有弯曲变形外，还有轴向变形和剪切变形的单元，杆件两端各有三个位移分量，这是平面结构杆件单元的一般情况。,符号规则：图(a)表示单元编号、杆端编号和局部座标，局部座标的,座标与杆轴重合；,1,2,EAI,l,(a),图(b)表示的杆端位移均为正方向。,单元编号杆端编号局部座标,1,2,(b),杆端位移编号,1,2,杆端力编号,(c),二、杆端位移、杆端力的正负号规定,一般单元：,2020/6/14,9,1,2,1,2,（1）单元杆端位移向量,（2）单元杆端力向量,凡是符号上面带了一横杠的就表示是基于局部座标系而言的。,2020/6/14,10,现在讨论单元刚度方程。单元刚度方程是指由单元杆端位移求单元杆端力时的一组方程，可以用“”表示，由位移求力称为正问题。,在单元两端加上人为控制的附加约束，使基本杆单元的两端产生任意指定的六个位移，然后根据这六个杆端位移来推导相应的六个杆端力。,1,2,我们忽略轴向受力状态和弯曲受力状态之间的相互影响，分别推导轴向变形和弯曲变形的刚度方程。,13-2单元刚度矩阵(局部座标系),进行单元分析，推导单元刚度方程和单元刚度矩阵。,一、一般单元,2020/6/14,11,分别推导轴向变形和弯曲变形的刚度方程。,首先，由两个杆端轴向位移,可推算出相应的杆端轴向力,1,2,其次，由杆端横向位移,可以用角变位移方程推导出相应的杆端,2020/6/14,12,将上面六个方程合并，写成矩阵形式：,2020/6/14,13,EAl,6EIl2,6EIl2,EAl,12EIl3,12EIl3,4EIl,2EIl,上面的式子可以用矩阵符号记为,这就是局部座标系中的单元刚度方程。,可求单元杆端力,0,0,0,0,0,0,6EIl2,0,6EIl2,0,-EAl,-6EIl2,-6EIl2,EAl,-12EIl3,12EIl3,2EIl,4EIl,0,0,0,0,0,0,-6EIl2,0,6EIl2,0,只与杆件本身性质有关而与外荷载无关,通过这个式子由单元杆端位移,局部座标系的单元刚度矩阵,2020/6/14,14,二、单元刚度矩阵的性质,（1）单元刚度系数的意义,代表单元杆端第j个位移分量等于1时所引起的第i个杆端力分量。,例如,代表单元杆端第2个位移分量时所引起的第5个杆端力分量的数值。,（2）单元刚度矩阵是对称矩阵，,（3）一般单元的刚度矩阵是奇异矩阵；,因此它的逆矩阵不存在,从力学上的理解是，根据单元刚度方程,由,有一组力的解答(唯一的)，即正问题。,由,如果,不是一组平衡力系则无解；若是一组平衡力系，则解答不是唯一的，即反问题。,2020/6/14,15,三、特殊单元若单元六个杆端位移中有某一个或几个已知为零，则该单元称为特殊单元，其刚度方程是一般单元刚度方程的特例。,以连续梁为例：,1,2,2020/6/14,16,为了程序的标准化和通用性，不采用特殊单元，只用一般单元，如果结构有特殊单元，可以通过程序由一般单元来形成。,2020/6/14,17,13-3单元刚度矩阵(整体座标系),座标转换矩阵,单元杆端力的转换式、单刚的转换式,一、单元座标转换矩阵,2020/6/14,18,正交矩阵,T-1=TT,或TTT=TTT=I,于是可以有,同理可以有,2020/6/14,19,在局部座标系中杆端力与杆端位移的关系式表达为：,在整体座标系中杆端力与杆端位移的关系式可以表达为：,二、整体座标系中的单元刚度矩阵,（a）式可转换为：,两边前乘TT,比较式(b)和(d)可得：,2020/6/14,20,三、单元刚度矩阵的性质,整体座标系中的单元刚度矩阵,2020/6/14,21,例1.试求图示刚架中各单元在整体座标系中的刚度矩阵k。设和杆的杆长和截面尺寸相同。,l=5m，bh=0.5m1m,解:,(1)局部座标系中的单元刚度矩阵,(2)整体座标系中的单元刚度矩阵,单元1：=0，T=I,单元2：=90，单元座标转换矩阵为,2020/6/14,22,单元2：=90，单元座标转换矩阵为,2020/6/14,23,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的贡献大小进行叠加而计算所得。,传统位移法,2020/6/14,24,一、单元集成法的力学模型和基本概念,分别考虑每个单元对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的贡献,略去其它单元的贡献。,2020/6/14,25,k=,4i2,2i2,4i2,2i2,F2,=,4i1,2i1,4i1,2i1,0,0,0,0,0,1,2,3,K,F=,设i1=0，则,单元的贡献矩阵,单元对结点力F的贡献,略去单元的贡献。,2020/6/14,26,k,K,K,F=F,+F,=（K,+K,）,1,2,F=K,整体刚度矩阵为：,单元集成法求整体刚度矩阵步骤：,根据单元和单元分别对结点力F的贡献，可得整体刚度方程：,2020/6/14,27,2i2,2i2,4i2,K=,4i1,2i1,4(i1+i2),2i1,0,2i2,0,2i2,4i2,4i1+4i2,整体刚度矩阵:,2020/6/14,28,二、按照单元定位向量由,(1)在整体分析中按结构的结点位移统一编码，称为总码。,(2)在单元分析中按单元两端结点位移单独编码，称为局部码。,以连续梁为例,位移统一编码，总码,单元,对应关系,局部码总码,单元定位向量,(1)1,(2)2,=,(1)2,(2)3,=,位移单独编码局部码,由单元的结点位移总码组成的向量,2020/6/14,29,单元,单元,(1),(2),(1),(2),4i1,2i1,2i1,4i1,k=,(1),(2),(1),(2),4i2,2i2,4i2,2i2,单元定位向量描述了单元两种编码（总码、局部码）之间的对应关系。,单元定位向量定义了整体坐标系下的单元刚度矩阵中的元素在整体刚度矩阵中的具体位置，故也称为“单元换码向量”。,单元贡献矩阵是单元刚度矩阵，利用“单元定位向量”进行“换码重排位”。,2020/6/14,30,三、单元集成法的实施,（定位累加）,K,（1）将K置零，得K=0；,（2）将k的元素在K中按定位并进行累加，得K=K；,（3）将k的元素在K中按定位并进行累加，得K=K+K；,按此作法对所有单元循环一遍，最后即得整体刚度矩阵K。,2020/6/14,31,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,例.求连续梁的整体刚度矩阵。,2020/6/14,32,四、整体刚度矩阵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,2020/6/14,33,13-5刚架的整体刚度矩阵,思路要点：（1）设各单元已形成了整体座标系下的单元刚度矩阵；,与连续梁相比：(1)各单元考虑轴向变形；(2)每个刚结点有三个位移；(3)要采用整体座标；(4)要处理非刚结点的特殊情况。,一、结点位移分量的统一编码总码,=1234T,规定：对于已知为零的结点位移分量，其总码均编为零。,=uAvAACT,整体结构的结点位移向量为：,相应地结点力向量为：,=XAYAMAMCT,F=F1F2F3F4T,2020/6/14,34,二、单元定位向量,单元,单元,局部码总码,局部码总码,（1）1（2）2（3）3（4）0（5）0（6）4,（1）1（2）2（3）3（4）0（5）0（6）0,三、单元集成过程,2020/6/14,35,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,2020/6/14,36,四、铰结点的处理,K,求单元常数,T,单元刚度矩阵,程序设计框图（局部：集成整体刚度矩阵）,刚结点：变形连续，截面1和截面2具有相同的结点位移。,铰结点：部分变形连续，截面1和截面2具有相同的结点线位移；而其角位移不相等。,2020/6/14,37,C1,C2,结点位移分量总码,结点C1,456,结点C2,457,单元定位向量,2020/6/14,38,K=,2020/6/14,39,13-6等效结点荷载,F=K(1),结构体系刚度方程：,一、位移法基本方程,K+FP=0.(2),F+FP=0.(3),将(1)式代入(2)式：,表示结点位移和结点力F之间的关系，反映了结构的刚度性质，而不涉及原结构上作用的实际荷载，并不是原结构的位移法基本方程。,基本体系在荷载单独作用下产生的结点约束力。,基本体系在结点位移单独作用下产生的结点约束力。,2020/6/14,40,二、等效结点荷载的概念,结点结束力FP,结点结束力FP,等效结点荷载P,原荷载,显然P=FP解决了计算等效结点荷载的问题,等效原则是两种荷载在基本体系中产生相同的结点约束力,K=F,FP,+,=,2020/6/14,41,三、按单元集成法求整体结构的等效结点荷载P,(1)局部座标单元的等效结点荷载P,(2)整体座标单元的等效结点荷载P,(3)结构的等效结点荷载P,2020/6/14,42,单元1:,单元2:,12,10,-10,+4,+0,-5,2020/6/14,43,K,求单元常数,T,P,原始数据、局部码、总码,解方程K=P求出结点位移,开始,结束,13-7计算步骤和算例,2020/6/14,44,例.求图示刚架的内力。设各杆为矩形截面，横梁b2h2=0.5m1.26m，立柱b1h1=0.5m1m。,（1）原始数据、局部码、总码（设E=1）,0,0,柱,梁,2020/6/14,45,2020/6/14,46,（2）形成局部座标系中的单元刚度矩阵,单元1和3,=10-3,=10-3,单元2,（3）计算整体座标系中的单元刚度矩阵,2020/6/14,47,单元1和3的座标转换矩阵（=900）,单元2（=0）,2020/6/14,48,(4)用单元集成法形成整体刚度矩阵K,2020/6/14,49,（5）求等效结点荷载P,单元固端约束力,单元1（=90）,按单元定位向量,2020/6/14,50,（6）解基本方程,求得结点位移:,（7）求各单元杆端力,2020/6/14,51,=10-3,同样可得出：,2020/6/14,52,（8）绘制内力图,A,B,C,D,8.49,2.09,3.04,4.38,M图(kNm),4.76,1.24,0.43,1.24,Q图(kN),N图(kN),0.43,0.43,1.24,2020/6/14,53,k3,k2,k1,13-8忽略轴向