结构力学--矩阵位移法--110页.ppt

,第六章 矩阵位移法,6.1 概 述,矩阵位移法是以结构位移为基本未知量， 借助矩阵进行分析，并用计算机解决各种杆系 结构受力、变形等计算的方法。,理论基础：位移法 分析工具：矩阵 计算手段：计算机,基本思想:,化整为零,- 结构离散化,将结构拆成杆件,杆件称作单元.,单元的连接点称作结点.,单元分析,对单元和结点编码.,单元杆端力,集零为整,- 整体分析,单元杆端力,结点外力,单元杆端位移,结点外力,单元杆端位移,(杆端位移=结点位移),结点外力,结点位移,基本未知量:结点位移,6.2 矩阵位移法解连续梁,一.离散化,-整体编码,1,2,3 -结点编码,(1),(2),(3) -结点位移编码,结点位移逆时针为整, 结点力逆时针为整.,二.单元分析,建立单元杆端力和单元杆端位移的关系.,-单元杆端力,1,2-局部编码,单元分析的目的:,1,2,-单元杆端位移,单元杆端力和单元杆端位移 逆时针为正.,=,=,+,简记为,-单元刚度方程,其中 称作单元刚度矩阵(简称作单刚),二.单元分析,单元刚度矩阵中元素的物理意义,-发生 位移时在 i端所需加的杆端力.,单元刚度矩阵性质:对称矩阵,三.整体分析,整体分析的目的: 建立结点力与结点位移的关系.,简记为,-结构刚度方程,-结构刚度矩阵(总刚),单元刚度矩阵中元素的物理意义,-发生 其它结点位 移为零位移时在 i结点所需 加的结点力.,结构刚度矩阵性质:对称矩阵,总刚的形成方法 -“对号入座”,四.计算杆端力,计算结点位移,计算杆端力,例: 计算图示梁,作弯矩图,解: 1.离散化,2.计算总刚,总荷,3.解方程,求位移,4.求杆端力,6,7/2,1/2,3,m,0,6,3,五.(零位移)边界条件处理,方法:,先处理法,后处理法,后处理法:,置0置1法,乘大数法,(1)置0置1法,0,0,0,1,0,m,1/2,1,2,1,(2)乘大数法,若 ,则将总刚主对角 元素 乘以大数n.,第三个方程变为:,6.2 矩阵位移法解连续梁,六.非结点荷载,(1).等效结点荷载,-结构等效结点荷载,=,+,“等效”是指等效结点荷载引起的结点位移与非结点荷载引起的结点位移相同,(2).等效结点荷载的计算,(2).等效结点荷载的计算,“对号入座”形成结构的 等效结点荷载,单元固端力:荷载引起的固端弯矩. 逆时针为正.记作,单元固端力改变符号称为单元等效结点荷载,记作,由单元等效结点荷载 “对号入座”可形成结 构等效结点荷载,(3).结构综合结点荷载,-直接结点荷载,-等效结点荷载,-结构综合结点荷载 (总荷),(4).最终杆端力计算,-计算结点位移,-计算杆端力,=,+,七.例题 矩阵位移法解图示梁,作m图.,解:,1.离散化,2.求总刚,3.求总荷,七.例题 矩位移法解图示梁,作m图.,4.边界条件处理,5.解方程,6.求杆端力,七.例题 矩位移法解图示梁,作m图.,6.求杆端力,7.作m图,1.29,27.43,19.43,9.71,八.边界条件的先处理法,解:,1.离散化,2.求总刚,3,+8,4,4,8,+3,4.解方程,先处理法,后处理法,其它过程同后处理法,3.求总荷,九. 无结点线位移的刚 架计算,1(0),2(1),3(2),4(0),解:,1.离散化,2.求总刚,3.求总荷,4.解方程,5.求杆端力,6.3 矩阵位移法解平面刚架,一.离散化,将结构离散成单元的分割点称作结点.,结点的选择:转折点、汇交点、支承点、 刚度变化、荷载作用点等,整体编码：单元编码、结点编码、 结点位移编码。,（1,2,3）,（4,5,6）,（7,8,9）,（10,11,12）,（13,14,15）,（16,17,18）,坐标系:整体(结构)坐标系;,局部(单元)坐标系.,曲杆结构:以直代曲.,变截面杆结构:以等截面杆 代变截面杆,二.单元分析,建立单元杆端力和 单元杆端位移的关系.,单元杆 端力,单元分析的目的:,单元杆 端位移,单元杆端力和单元杆端位移 的方向与局部坐标系一致为正.,二.单元分析,则有:,若令:,二.单元分析,则有:,若令:,练习: 试求单刚第 三列元素.,二.单元分析,则有:,若令:,练习: 试求单刚第 三列元素.,二.单元分析,单刚的性质:,1.对称矩阵,2.奇异矩阵,单刚的分块矩阵表示:,三.坐标转换,1.问题的提出,2.整体坐标系下的杆端力与 局部坐标系下的杆端力之 间的关系,局部坐标系下的杆端力,整体坐标系下的杆端力,简记为:,其中,可直接验证坐标转焕矩阵是一个正交矩阵.,即,对于杆端位移有相同的关系:,3.整体坐标系下的单元刚度矩阵,-整体坐标系下的单元刚度方程,其中,-整体坐标系下的单元刚度矩阵,(整体单刚),4.整体单刚的计算,4.整体单刚的计算,已知:,求:各单元整体单刚,解:,利用物理意义求:,令:,则有:,作业:用此方法求2单 元整体单刚其它元素.,四.整体分析(后处理法),1.结点力与结点位移的关系,结点力,结点位移,由变形协调条件,有,同理,有,由结点1平衡条件:,由结点2平衡条件:,由结点3平衡条件:,-结构原始刚度矩阵,-结构原始刚度矩阵,2.结构原始总刚的形成,采用“对号入座”的方法,3.单刚子块在总刚中的分布规律,主子块-主对角线上的子块.,付子块-非主对角线上的子块.,相关结点-有单元相连的结点.,相关单元-与结点相连的单元称为 该结点的相关单元.,规律:,(1)若i,j为相关结点, 为连接 i,j结点的单元单刚的相应付 子块;若不是相关结点,(2)主子块 为 i 结点的相关单 元单刚主子块之和.,2.结构原始总刚的形成,3.单刚子块在总刚中的分布规律,4.总刚中元素的物理意义,则有:,若令:,练习:求总刚中第四列元素,则有:,若令:,(1).对称性,5.原始总刚的性质,(2).奇异性,(3).稀疏性,6.总刚的半带存贮,带宽,变带宽存贮,半变带宽存贮,等带宽存贮,半等带宽存贮,等半带宽计算公式:,d=(相关结点最大差值+1),半带宽与结点编码有关,19层,1,2,3,4,39,40,总刚占用存贮单元:,1,2,20,21,40,19,22,总刚占用存贮单元:,五.边界条件处理,已知某结点位移,1.乘大数法,做法:取大数n,总刚中元素 乘以n;并用 替换,2.置零置1法,做法:(1)用 中的第i列代替,(2)将总刚中第i行第i列的非主对角元素置0;,(3)将总刚中主对角元素 置为1,总荷中元素 置成c,经边界条件处理后的总刚称为结构刚度矩阵.,结构刚度矩阵是一个非奇异矩阵.,六.非结点荷载处理,-单元固端力,-单元等效结点荷载,例.计算图示结构的等效结点荷载,例.计算图示结构的等效结点荷载,例.计算图示结构的等效结点荷载,按物理意义求解:,1.划分单元、编码； 2.形成结构原始总刚； 3.形成综合结点荷载； 4.引入边界条件； 5.解方程； 6.求杆端力。,七.杆端力计算,八.计算步骤及程序流程,九.先处理法,2(4,5,6),1(1,2,3),3 (7,8,9),1(0,0,1),2(0,2,3),3 (0,0,0),4(4,5,6),后处理法,先处理法,总刚的形成:,4(10,11,13),九.先处理法,1(0,0,1),2(0,2,3),3 (0,0,0),4(4,5,6),-单元定位向量,九.先处理法,-2单元定位向量,4(4,5,6),九.先处理法,4(4,5,6),十.不计轴变刚架的计算,计轴变时的结点位移编码,不计轴变时的结点位移编码,梁单元的单刚,柱单元的单刚,局部单刚与梁相同.,解:,解:,6.4 矩阵位移法解平面桁架,一.离散化,二.局部坐标系下单元刚度方程,二.局部坐标系下单元刚度方程,局部坐标系单元刚度方程,局部单刚,三.整体坐标系下单元刚度方程,-整体部坐标系单元刚度方程,-整体单刚,四.整体分析,-整体单刚,整体分析及求杆端力与刚架类似.,例:矩阵位移法求图示桁架各杆轴力.,已知:ea=6 0 , p=100,解:,-整体单刚,-整体单刚,五.组合结构的计算,采用两种单元,其它过程与前类似.,例:矩阵位移法求图示桁架各杆轴力.,已知:ea=6 0 , p=100,解:,五.组合结构的计算,采用两种单元,其它过程与前类似.,例:矩阵位移法求图示桁架各杆轴力.,已知:ea=6 0 , p=100,解:,五.组合结构的计算,采用两种单元,其它过程与前类似.,例:矩阵位移法求图示桁架各杆轴力.,已知:ea=6 0 , p=100,解:,杆端力计算与前相同,杆端力计算与前相同,1,2,3,作业:1.对图示结构作结点位移编码(先处理), 并写出各单元的单元定位向量.,4,5,1,2,4,5,3,6,2. 215页 6-2,练,l/2,l/2,c,q,m,n,mp,mi,mp,矩阵位移法习题讨论,一.离散化,1.不计轴变时先处理法的结点位移编码,2. 计轴变时先处理法的 结点位移编码,二. 单元分析,1.单元刚度方程表示什么量之间的关系方程?,2.单元刚度矩阵(自由式单元)是什么样的矩阵?,3.单刚元素 的物理意义是什么?,4.坐标转换矩阵是一个什么样的矩阵?,5.局部坐标系下的杆端位移与整体坐标下的有何关系?,6.单元刚度矩阵均是奇异矩阵吗?,7.试写出自由式单元坐标转换矩阵.,二. 单元分析,8.求图示结构2单元的坐标转换矩阵中的元素,二. 单元分析,9.试写出桁架单元坐标转换矩阵中的第二行元素.,三. 整体分析,1.结构刚度方程 是整体结构所应满足的变形 协调条件吗?,2.总刚元素 的物理意义是什么?,3.试写出图示刚架2单元的单元定位向量.,4.图示结构2单元的整体单刚元素 应放在总刚的什么位置?,三. 整体分析,4.图示结构2单元的整体单刚元素 应放在总刚的什么位置?,第5行第6列,三. 整体分析,5.试求总刚元素 ea=常数,6.先处理法求图示结构总刚 (不计轴变),三. 整体分析,6.先处理法求图示结构总刚 (不计轴变),7.先处理法求图示结构总刚 (不计轴变),三. 整体分析,7.先处理法求图示结构总刚 (不计轴变),三. 整体分析,8.等效结点荷载的数值等于汇交于该结点的所有单元 固端力之和. 此结论对否?,9.试求图示结构的荷载列阵(先处理法).,三. 整体分析,10.试求图示结构的荷载列阵(先处理法).,11.试求图示结构(不计轴变)的荷载列阵(先处理法).,三