结构力学 矩阵位移法.pptx

《结构力学 矩阵位移法.pptx》由会员分享，可在线阅读，更多相关《结构力学 矩阵位移法.pptx（192页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、结构矩阵分析方法是电子计算机进入结构力学领域而产生的一种方法。第1页/共192页结构力学传统方法与结构矩阵分析方法，二者同源而有别：在原理上同源，在作法上有别。简单地说，前者在“手算”的年代形成，后者则着眼于“电算”，计算手段的不同，引起计算方法的差异。第2页/共192页与传统的力法、位移法相对应，在结构矩阵分析中也有矩阵力法和矩阵位移法，或称柔度法与刚度法。矩阵位移法由于具有易于实现计算过程程序化的优点而广为流传，本章只对矩阵位移法进行讨论。第3页/共192页矩阵位移法是有限元法的雏形，因此结构矩阵分析有时也称为杆件结构的有限元法。在本章中将使用有限元法中的一些术语和提法。第4页/共192页

2、先把整体拆开，分解成若干个单元（在杆件结构中，一般把每个杆件取作一个单元），这个过程称作离散化。然后再将这些单元按一定条件集合成整体。在一分一合，先拆后搭的过程中，把复杂结构的计算问题转化为简单单元的分析和集合问题。有限元法的要点：第5页/共192页有限元法包含两个基本环节：1.单元分析 2.整体分析 第6页/共192页单元分析的任务:建立单元刚度方程，形成单元刚度矩阵 第7页/共192页整体分析的主要任务:将单元集合成整体，由单元刚度矩阵按照刚度集成规则形成整体刚度矩阵，建立整体结构的位移法基本方程，从而求出解答。第8页/共192页本节和下一节对平面结构的杆件单元进行单元分析，得出单元刚度方

3、程和单元刚度矩阵。第9页/共192页位移法中给出的转角位移方程实际上就是梁单元的刚度方程。梁单元是杆件单元的特例。本节推导单元刚度方程时有几点新的考虑：重新规定正负号规则，讨论杆件单元的一般情况，采用矩阵表示形式。第10页/共192页1.一般单元 图11-1所示为平面刚架中的一个等截面直杆单元 图111第11页/共192页设杆件除弯曲变形外，还有轴向变形。左右两端各有三个位移分量（两个移动、一个转动），杆件共有六个杆端位移分量，这是平面结构杆件单元的一般情况。由端点1到端点2的方向规定为杆轴的正方向，在图中用箭头表明。第12页/共192页图中采用坐标系，轴与杆轴重合。这个坐标系称为单元坐标系或

4、局部坐标系。上面都划上一横，作为局部坐标系的标志。第13页/共192页在局部坐标系中，一般单元的每端各有三个位移分量 和对应的三个力分量 第14页/共192页图11-2中所示的位移、力分量方向为正方向。图112 第15页/共192页单元的六个杆端位移分量和六个杆端力分量按一定顺序排列，形成单元杆端位移向量 和单元杆端力向量 如下：（111）第16页/共192页向量中的六个元素的序码记为（1），（2），（6）。由于它们是在每个单元中各子编码的（不是在刚架所有单元中统一编码的），因此称为局部码杆端位移分量（或杆端力分量）的局部码。数码（1），（2），都加上括号，作为局部码的标志。第17页/共192

5、页单元刚度方程是指由单元杆端位移求单元杆端力时所建立的方程记为 第18页/共192页为了建立单元刚度方程，我们按照位移法基本体系的作法，在杆件两端加上人为控制的附加约束，使基本体系在两端发生任意指定的位移 如下图11-3所示。然后根据 推算相应的杆端力 第19页/共192页图11-3忽略轴向受力状态和弯曲受力状态之间的相互影响，分别推导轴向变形和弯曲变形的刚度方程。第20页/共192页首先，由杆端轴向位移 可推算出相应的杆端轴向力（112）第21页/共192页其次，由杆端横向位移 和转角 可推算出相应的杆端横向力 和杆端力矩 第22页/共192页根据转角位移方程（85）和（86），并改用本章的

6、记号和正负号，即得（85）（86）（113）第23页/共192页上面六个刚度方程（112）和（113）实际上在位移法中已经推导过。现在将它们合在一起，写成矩阵形式如下：（114）第24页/共192页上式可记为：（115）（116）其中 第25页/共192页式（115）即为所求的 称为在局部坐标系中的单元刚度方程。矩阵 称为局部坐标系中的单元刚度矩阵。它是 方阵。第26页/共192页2单元刚度矩阵的性质 1）单元刚度系数的意义 中的每个元素称为单元刚度系数，代表由于单位杆端位移所引起的杆端力。例如，第（6）行第（3）列元素（即元素 ）代表当第（3）个杆端位移时引起的第（6）个杆端力分量 第27页

7、/共192页一般来说，第（i）行第（j）列元素 代表当第（j）个杆端位移分量 等于1（其他位移分量为零）时所引起的第（i）个杆端力分量 的值。第28页/共192页 中某一列的六个元素分别表示当某个杆端位移分量等于1时所引起的六个杆端力分量。例如，第1列对应于单位位移 所引起的杆端力。为了帮助理解，在式（116）中，在 每一列的上方都标明了对应的单位位移分量。第29页/共192页2）是对称矩阵 的对称性是指其元素有如下关系：（117）这实际上是根据反力互等定理得出的结论。第30页/共192页3）一般单元的是奇异矩阵 的奇异性是指其行列式等于零，即（118）直接计算式（116）的矩阵行列式，便可验

8、证上述结论。第31页/共192页由此可知，不存在逆矩阵。也就是说，根据单元刚度方程（115），可以由杆端位移 推算出杆端力 且 的解是唯一解；但不能由杆端力 反推出杆端位移 ，可能无解，如有解，则为非唯一解。第32页/共192页为了避免混淆，我们把正反两个问题再从数学提法、力法模型、解的性质等方面作一对比。见下表：第33页/共192页正问题反问题数学提法为任意指定值，为待求量。为任意指定值，为待求量。力学模型把单元按“两端有六个人工控制的附加约束的杆件”（位移法基本体系）来分析由控制附加约束而加以指定。把单元按“两端自由的杆件”来分析直接加在自由端作为指定的杆端力。解的性质为任意值时，都有解，

9、且为唯一解。总是一个平衡力系，不可能是不平衡力系。为不平衡力系时，没有解。为平衡力系时，有解，但为非唯一解（因为自由杆件除本身变形外还可有任意刚体位移）。不存在。第34页/共192页总之，正反两个问题的力学模型是截然不同的，不能把单元笼统地统称为“自由单元”。逆矩阵 的性质是跟据反问题确定的，这里的反问题是按“自由单元”分析，故得出 不存在的结论。第35页/共192页3.特殊单元 式（114）是一般单元的刚度方程，其中六个杆端位移可指定为任意值。在结构中还有一些特殊单元，单元的某个或某些杆端位移的值已知为零，而不能任意指定。各种特殊单元的刚度方程无需另行推导，只需对一般的单元刚度方程（114）

10、作一些特殊处理便可自动得到。第36页/共192页举例来说，计算连续梁时，我们通常忽略轴向变形。如取每跨梁作为单元（图114），则只有两个杆端位移分量 可指定为任意值，而其余四个分量均已知为零：图114 返回第37页/共192页（a）将式（a）代入式（114），即自动得出此特殊单元的刚度方程如下：（119）第38页/共192页此时单元刚度矩阵为（1110）返回第39页/共192页在结构矩阵分析中，我们着眼于计算过程的程序化、标准化和自动化。因此只采用一种标准化形式一般单元的刚度矩阵（116），关于单元刚度矩阵的各种特殊形式将由计算机程序去自动形成。第40页/共192页某些特殊单元的刚度矩阵是可逆

11、的。例如式（1110）中的 ，其逆矩阵存在。对于图114所示特殊单元来说，正问题的力学模型如图115a所示，每端有两个支杆和一个控制转角的附加约束，可指定为任意值。和 第41页/共192页图115返回返回第42页/共192页反问题的力学模型如图115b所示，每端有两个支杆，杆端力矩 为任意值。由于反问题的力学模型是一个几何不变体系，因此，当 为任意值时，杆端转角 有解，且为唯一解。由此得出 存在的结论。第43页/共192页选用局部坐标系的目的是希望导出的单元刚度矩阵具有最简单的形式。为了便于进行整体分析，必须选用一个统一的公共坐标系，称为整体坐标系。为了区别，用 表示局部坐标，用 表示整体坐标

12、。第44页/共192页单元坐标转换矩阵（1）图116a所示为一单元e，局部坐标系中的杆端力分量用 表示。整体坐标系中则用 表示，如图116b所示。第45页/共192页图116返回第46页/共192页显然，二者有下列关系：（1111）第47页/共192页将式（1011）写成矩阵形式：（1112）或简写成（1113）返回第48页/共192页式中T 称为单元坐标转换矩阵（1114）式（1113）是两种坐标系中单元杆端力的转换式。第49页/共192页T1T T（1115）或 T TTT T TI（1116）式（1113）的逆转换式为（1117）设局部坐标系中单元杆端位移列阵为 ，整体坐标系中单元杆端位

13、移列阵为 ，则返回第50页/共192页（1118）（1119）整体坐标系中的单元刚度矩阵（2）单元杆端力与杆端位移在整体坐标系中的关系式可写为（1120）返回第51页/共192页 单元e在局部坐标系中的刚度方程为（a）将式（1113）和（1118）代入式（a），得到等式两边各前乘 ，并引入式（1116），得（b）第52页/共192页表较式（b）与（1120），可知（1121）整体坐标系中的单元刚度矩阵与同阶，具有类似的性质：元素 表示在整体坐标系中第（j）个杆端位移分量等于1时引起的第(i)个杆端力分量。（1）是对称矩阵。（2）一般单元的是奇异矩阵。（3）第53页/共192页例 111 试求图

14、117所示刚架中各单元在整体坐标系中的刚度矩阵 。设各杆的杆长和截面尺寸相同。图117 第54页/共192页解:局部坐标系中的单元刚度矩阵由式（116）得=（1）第55页/共192页整体坐标系中的单元刚度矩阵（2）单元：=单元：单元坐标转换矩阵为 第56页/共192页=T=第57页/共192页11-4连续梁的整体刚度矩阵 前两节进行了单元分析，建立了单元刚度方程，推导了单元刚度矩阵。从本节起转到整体分析，建立整体刚度方程，导出整体刚度矩阵。本节以连续梁为例，下节讨论刚架的一般情况，并考虑杆件轴向变形的影响。第58页/共192页整体刚度方程是按位移法建立的,具体做法有两种：1.传统位移法 2.单

15、元集成法（也称为刚度集成法或直接刚度法）单元集成法的优点是便于实现计算过程的程序化。第59页/共192页对于图11-8a所示的连续梁，位移法基本体系如图118b所示。图118 第60页/共192页位移法的基本未知量为节点转角他们可指定为任意值，在基本体系中用控制附加约束加以指定。他们组成整体结构的节点位移向量：（）T与 对应的力是附加约束的力偶 它们组成整体结构的结点力向量F第61页/共192页在传统作法中，分别考虑每个结点转角 独自引起的节点力偶，如图11-9a bc所示。第62页/共192页图119第63页/共192页叠加得结点力偶；如下（1122）记为：（1124）（1123）式（11-

16、22）或(11-23)称为整体刚度方程，K称为整体刚度矩阵。第64页/共192页1.单元集成法的力学模型和基本概念 传统位移法求结构的结点F时，分别考虑每个结点位移对F的单独贡献（采用图11-9中的力学模型），然后进行叠加。单元集成法求F时，分别考虑每个单元对F的单独贡献，然后进行叠加其特点就是“由单元直接集成”。第65页/共192页首先，考虑单元的贡献，力学模型见图11-10。整个结构的结点力是由单元单独产生的，记为 F=(F1 F2 F3)T图1110第66页/共192页F1表示单元对结构结点力F的贡献。F1和F2可由单元的单元刚度矩阵k算出。已知F30 （a）k（1025）得：（b）第6

17、7页/共192页由（a）和（b）得：（1026）记为F3K（1027）其中 K（1028）K表示单元对刚度矩阵提供的贡献，称为单元的贡献矩阵。第68页/共192页其次，考虑单元的贡献。力学模型见图11-11所示。图1111已知 k 故得 记为F=K（1129）（1130）（1131）第69页/共192页其中 K（1132）K称为单元的贡献矩阵。将式（11-27）和式（11-31）叠加，得：F=F+F=(k+k)（1133）由此得出整体刚度矩阵K为K=K+K=（1134）第70页/共192页单元集成法求整体刚度矩阵的步骤可表示为 其中：为单元刚度矩阵，单元贡献矩阵，K为整体刚度矩阵第71页/共1

18、92页2.按照单元定位向量由 求 注意以下3点：1)结点位移（或结点力）有两种编码：在整体分析中，结点位移在结构中统一进行编码，称为总码。在单元分析中，每个单元的两个结点位移各自编码为（1）和（2），称为局部编码。（见下图11-12）第72页/共192页（b）（a）图1112第73页/共192页2)注意每个单元的结点位移分量两种编码之间的对应关系，具体见下表：单元对应关系单元定位向量局部码总码（1）1（2）2（1）2（2）3第74页/共192页3)注意单元刚度矩阵 和单元贡献矩阵 中元素的排列方式,见下表 在单元刚度 矩阵中在单元贡献矩阵中换码元素的原行码(i)原列码（j）换成新行码新列码（i

19、）（j）重排座原排在(i)行（j）的元素改在行列第75页/共192页总之，由 求 的问题实质上就是 中的元素在 中如何定位的问题。定位规则是：（1136）参见下表:第76页/共192页单元单元刚度矩阵单元定位向量单元贡献矩阵（1）（2）（1）4i12i1（2）2i14i1(1)(2)123（1）14i12i10（2）22i14i103000（1）（2）（1）4i22i2（2）2i24i2(1)(2)123（1）1000（2）204i22i2302i24i2第77页/共192页3.单元集成法的实施方案 单元集成法形成K的过程：1)先将K置零，这时K=0 2)将k的元素在K中按定位 并进行累加，这

20、时，K=K 3)将k的元素在K中按定位 并进行累加，这时，K=K+K按此作法对所有单元循环一遍，最后得到第78页/共192页现以图11-8a所示连续梁为例，说明上过程：将k集成后，得到：在此基础上将k集成得最终结果：第79页/共192页例11-2 试求图11-13a所示连续梁得整体刚度矩阵K 解（1）结点位移分量总码（见图11-13a）图1113第80页/共192页（2）各单元得定位向量 第81页/共192页（3）单元集成过程 单元单元刚度矩阵按单元定位向量换码 集成过程中得阶段结果 （1）（2）（1）4i1 2i1（2）2i1 4i1（1）1（2）2 1 2 3（1）1 4i1 2i1 0（

21、2）2 2i1 4i1 0 3 0 0 0 （1）（2）（1）4i2 2i2（2）2i2 4i2(1)2(2)3 1 2 3（1）1 4i1 2i1 0（2）2 2i1 4i1+(4i2)2i2 3 0 2i2 4i2 （1）（2）（1）4i3 2i3 （2）2i3 4i3(1)3(2)0(2)0 1 2 3（1）1 4i1 2i1 0（2）2 2i1 4i1+(4i2)2i2 3 0 2i2 4i2+4i 第82页/共192页4.整体刚度矩阵的性质（1）整体刚度系数的意义 K中的元素 称为整体刚度系数。它表示当第j个结点位移分量 （其他结点位移分量为零）时所产生的第i个结点力(2)K是对称矩

22、阵(3)按本节方法计算连续梁时，K时可逆矩阵。第83页/共192页图11-8a所示为下图11-14为例的反问题力学模型。当F为指定值时，均可得 的唯一解，故 是存在的。图1114（4）K是稀疏矩阵和带状矩阵。第84页/共192页对下图11-15，可导出其整体刚度矩阵：图1115（1137）第85页/共192页11-5 刚架的整体刚度矩阵 本节讨论用单元集成法求平面刚架的整体刚度矩阵K。思路的要点：K由 直接集成；集成包括将 的元素在K中定位和累加两个环节；定位是依据单元定位向量 进行的。第86页/共192页情况的复杂性表现在下列几个方面：1）在一般情况下要考虑刚架中各杆的轴向变形，而忽略杆件轴

23、向变形的情况则作为特例来处理；2）刚架中每个结点的位移分量要增加到三个：角位移和两个方向的线位移；3）刚架中各杆方向不尽相同，在整体分析中需采用整体坐标；4）刚架中除刚结点外，还要考虑铰结点等其他情况。第87页/共192页1.结点位移分量的统一编码总码 图11-16所示刚架整体结构的结点位移向量 ：（）T（）T相应结点力向量为F=(F1F2F3F4)T第88页/共192页图1116第89页/共192页2.单元定位向量 单元单元局部编码总码单元定位向量局部编码总码单元定位向量（1）1（2）2（3）3（4）0（5）0（6）4（1）1（2）2（3）3（4）0（5）0（6）0第90页/共192页3.单

24、元集成过程 首先，考虑单元；（1）（2）（3）（4）（5）（6）（1）300kN/m00-300kN/m00（2）012kN/m30kN0-12kN/m30kN（3）4（4）-300kN/m00300kN/m00（5）0-12kN/m-30kN012kN/m-30kN（6）（1138）第91页/共192页K的阶段结果(1)(2)(3)(6)(1)1300kN/m000(2)2012kN/m30kN30kN104（1139）第92页/共192页其次，考虑单元k=（1）（2）（3）（4）（5）（6）（1）12kN/m0-30kN-12kN/m0-30kN（2）0300kN/m00-300kN/m0

25、（3）4（4）-12kN/m030kN12kN/m030kN（5）0-300kN/m00300kN/m0（6）（1140）第93页/共192页K=(1)(2)(3)1234(1)1300+(12)kN/m0+(0)0+(-30)kN0(2)20+(0)12+(300)kN/m30+(0)kN30kN104（1141）第94页/共192页4铰结点的处理图11-17所示为具有铰结点的刚架。图1117第95页/共192页1.考虑结点位移分量的统一编码（图中已标出）2.考虑单元定位向量各单元定位向量如下；分析过程：（123456）T（123000）T（457000）T第96页/共192页第一阶段结果，

26、见式（1143）。（1）（2）（3）（4）（5）（6）（1）1（2）2（3）3（4）4（5）5(6)6730000-300000012300-123000301000-30500-300003000000-12-30012-300030500-30100000000001234567（1143）第97页/共192页（1）（2）（3）1234567（1）1（2）24（3）3300+（12）0+（0）0+（30）-3000000+（0）12+（300）30+（0）0-123000+（30）30+（0）100+（100）0-30500-300003000000-12-30012-300030500-

27、3010000000000567（1144）第98页/共192页在式（1140）中已给出。将其中的元素按在K中定位并与前阶段结果累加，即得K的第二阶段结果，见式（1144）。其次，考虑单元：最后，考虑单元：与 相同。由 即得K最后结果，见式（1145）。第99页/共192页（1）（2）（3）（1）（2）（3）123456712345673120-30-3000000312300-12300-30302000-30500-30000300+（12）0+（0）00+（-30）0-12-300+（0）12+（300）-300+（0）030500-3010000000+（30）0+（0）00+（100

28、）（1145）第100页/共192页以上式（1143）、（1144）、（1145）中各物理量是有单位的，这里只是表示单元集成的过程，故式中未标单位。第101页/共192页1.1.位移法基本方程 116 等 效 结 点 荷 载第102页/共192页前两节讨论了结构的整体刚度矩阵K，建立了整体刚度方程（1146）整体刚度方程（1146）是根据原结构的位移法基本体系建立的，它表示由结点位移推算结点力（即在基本体系的附加约束中引起的约束力）F的关系式。它只反映结构的刚度性质，而不涉及到原结构上作用的实际荷载。它并不是用以分析原结构的位移法基本方程。第103页/共192页为了建立位移法基本方程，我们回顾

29、一下本书（）85中的推导方法，分别考虑位移法基本体系的两种状态：(1)设荷载单独作用（结点位移 设为零）此时在基本结构中引起的结点约束力，记为 。(2)设结点位移 单独作用（荷载设为零）此时在基本结构中引起的结点约束力为 。第104页/共192页位移法基本方程为（1147）2.等效结点荷载的概念 等效的原则是要求这两种荷载在基本结构中产生相同的结点约束力。第105页/共192页如果原来荷载在基本结构中引起的结点约束力记为，则等效结点荷载在基本结构中引起的结点约束力也应为。由此即可得出如下结论：（1148）将式（1148）代入式（1147），则位移法基本方程可写为（1149）第106页/共192

30、页由式（1146）和式（1149）可知，如果把刚度方程（1146）中的结点约束力F换成等效节点荷载P，即得到位移法基本方程（1149）。3按单元集成法求整体结构的等效结点荷载(1)单元的等效结点荷载 先考虑局部坐标系。第107页/共192页在单元两端加上六个附加约束，使两端固定。在给定荷载作用下，可求出六个固端约束力，它们组成固端约束力向量 ：（1150）在表111中给出了几种典型荷载所引起的固端约束力。将固端约束力 反号，即得到单元等效结点荷载 （局部坐标系）：第108页/共192页（1151）(2)单元的等效结点荷载（整体坐标系）现考虑整体坐标系。由坐标转换公式（1117），得（1052）

31、(3)整体结构的等效结点荷载P依次将每个 中的元素按单元定位向量 在P中进行定位并累加，最后即得到P。第109页/共192页表 111单元固端约束力（局部坐标系）荷载简图始 端 1末 端 212第110页/共192页34第111页/共192页56第112页/共192页7第113页/共192页例113试求图1116a所示刚架在图1118给定荷载下的等效结点荷载向量P。图 1018第114页/共192页解 （1）求局部坐标系中的固端约束力 单元：由表111第1行，得第115页/共192页单元：由表111第2行，得第116页/共192页因此(2)各单元在整体坐标系中的等效结点荷载 单元、的倾角分别为

32、 第117页/共192页由式（1151）和（1152）得 第118页/共192页(3)求刚架的等效结点荷载 两个单元的结点局部和总码见图1116。总码在图1118中用虚线重新示出。单元定位向量已知为第119页/共192页将 中得元素，按 在 中进行定位并累加即可得出 。首先，考虑单元：的阶段结果为（4）、（5）行元素在 中无座位（1）1（2）2（3）3(6)4第120页/共192页其次,考虑单元（1）1（2）2（3）34117计算步骤和算例 用矩阵位移法计算平面刚架的步骤如下：第121页/共192页(1)整理原始数据，对单元和刚架进行局部编码和总体编码。(2)形成局部坐标系中的单元刚度矩阵 ，

33、用式（116）。(3)形成整体坐标系中的单元刚度矩阵 ，用式（1121）。(4)用单元集成法形成整体刚度矩阵 ，参看式（1135）。第122页/共192页(5)求局部坐标系的单元等效结点荷载 ，转换成整体结构的等效结点荷载 ，用式（1151）和式（1152）；用单元集成法形成整体结构的等效结点荷载 。(6)解方程 ，求出结点位移 。(7)求各杆的杆端内力 ，用下面的式（1153）。各杆的杆端内力是由两部分组成：第123页/共192页一部分是在结点位移被约束住的条件下的杆端内力，即各杆的固端约束力。另一部分是刚架在等效结点荷载作用下的杆端内力，可由式（115）求出。将两部分内力叠加，即得（115

34、2）例 114试求图1119a所示刚架得内力。设各杆为矩形截面，横梁 立柱 第124页/共192页解 （1）原始数据及编码 原始数据计算如下（为了计算上得方便，设 ）。图1119第125页/共192页柱：梁：第126页/共192页(2)形成局部坐标系中的单元刚度矩阵 单元和：第127页/共192页单元：(3)计算整体坐标系中的单元刚度矩阵 第128页/共192页单元和的坐标转换矩阵为第129页/共192页单元：(4)用单元集成法形成整体刚度矩阵K由图1119b中单元局部编码与结点位移统一编码的关系，各杆的单元定位向量可写出如下：第130页/共192页按照单元定位向量 ，依次将各单元 中的元素在

35、K中定位并累加，最后得到K如下：第131页/共192页(5)求等效结点荷载P首先，求单元固端约束力 ：只有单元有 ，第132页/共192页其次，求单元在整体坐标系中的等效结点荷载 ：单元的倾角 ，由式（1151）和（1152）第133页/共192页按单元定位向量，将中的元素在中定位，得第134页/共192页(6)解基本方程 第135页/共192页求得 第136页/共192页(7)求各杆杆端力 第137页/共192页第138页/共192页单元：第139页/共192页单元：第140页/共192页第141页/共192页(8)根据杆端力绘制内力图，如图1120所示。图1120第142页/共192页11

36、8忽略轴向变形时矩形刚架得整体分析 以图1121所示矩形刚架为例进行说明。图1021第143页/共192页首先，考虑结点位移分量统一编码：在固定端B和D处，三个位移分量都为零，总码编为（000）。在刚结点A处，在铰结点 和 处，竖向位移分量都为零，故其编码也应为零码。结点的总码编为（102）结点的总码编为（103），结点的总码编为（104）第144页/共192页其次，考虑单元定位向量：在图1121中，单元、的 轴正方向用箭头标明。各杆的单元定位向量可由图直接写出如下：（a）第145页/共192页最后，按单元、的次序进行单元集成。设刚架中各杆尺寸相同，均采用例111中各杆的尺寸 首先，考虑单元：

37、在式（1138）中已给出。根据式（a）中的 进行定位，即得 的第一阶段结果如下：第146页/共192页（1）、（4）1（3）2（6）34（1）、（4）（3）（6）300-300-300+3000+00+000100500050100000001234其次，考虑单元：第147页/共192页 在式（1040）中已给出。将其中的元素按 在 中定位并与前阶段结果累加，即得 的第二阶段结果如下：（1）1（3）234（1）（3）12340+（12）（-30）00（-30）100+（100）50005010000000第148页/共192页最后，考虑单元:与相同。12+12-300-30-302005000

38、501000-30001001234（1）（3）1234（1）（3）即得k的最后结果如下：由 第149页/共192页例11-5试求图11-19a所示刚架的内力，截面尺寸同例11-4。忽略轴向变形的影响。解（1）原始数据及编码图11-22第150页/共192页除结点C和D为固定端，用0编码外，结点A和B的竖向位移均为零，用0编码。结点A和B的水平位移相同，编码同为1。（2）形成局部坐标系中的单元刚度矩阵（3）计算整体坐标系中的单元刚度矩阵由图11-22中单元局部编码与结点位移的统一编码的对应关系，各杆的单元定位向量可写出如下：第151页/共192页按照单元定位向量，依次将各单元中的元素在K中定位

39、并累加，最后得到K如下：（5）求等效结点荷载P单元在整体坐标系中的等效结点荷载同11-4按单元定位向量，将中的元素在P中定位得第152页/共192页（6）解基本方程求得（7）求各杆杆端力单元：第153页/共192页第154页/共192页单元：单元：第155页/共192页（8）根据杆端弯矩和剪力绘制M图和图M图和图如图11-23a、b所示。与图11-20相比，可见轴向变形的影响不大。由于假设杆件的轴向变形为零，因此根据刚度方程求出的杆端轴力为零。图11-23c中的图是根据平衡条件由得出的。第156页/共192页图11-23第157页/共192页119桁架及组合结构的整体分析1.桁架桁架单元的刚度

40、方程（局部坐标系，图11-24）已在式（11-2）中给出，其矩阵形式为（a）图11-24第158页/共192页对于斜杆单元，其轴力和轴向位移在整体坐标系中将有沿x轴和y轴的两个分量。因此，整体坐标系中的杆端力向量和杆端位移向量为（图11-25）图11-25第159页/共192页，可以从图11-25导出局部坐标系和整体坐标系中杆端力向量和杆端位移向量间的转换关系。11-3中对一般单元导出的转换关系具有一般性。桁架单元是一般单元的一个特殊情况，只需按照桁架单元的特点，把11-3中坐标系转换矩阵T修改一下，则11-3中的其他转换公式均可应用。第160页/共192页为了便于利用以前的坐标转换关系，我们

41、先将局部坐标系中的单元刚度方程（a）扩大为四阶的形式：这里，在 和 中引入了 、和 、，在 中添上了相应的零元素。式（11-54）与式（a）是等价的。（11-54）第161页/共192页式（11-12）是一般单元杆端力的转换式。对于桁架单元，由于 ，所以，删去坐标转换矩阵T式（11-14）中相应的行和列，便得到桁架单元的坐标转换矩阵T 如下：（11-55）单元集成法求整体刚度矩阵的步骤仍为式（11-35）。编码时应注意，桁架单元的结点转角不是基本未知量。第162页/共192页例11-6 试求图11-26a所示桁架的内力。各杆EA 相同。解（1）单元和结点位移分量的统一编码如图11-26b中。单

42、元的局部坐标用箭头方向表示，示于图11-26b中。整体坐标示于图11-26b中。图11-26第163页/共192页(2)形成局部坐标系中的单元刚度矩阵按四阶方阵的形式式（11-54）形成各单元在局部坐标系中的单元刚度矩阵 。第164页/共192页(3)形成整体坐标系中的单元刚度矩阵单元和单元：。由式（11-55）得单元和单元：得第165页/共192页单元：，由式（11-55）第166页/共192页单元：，由式（11-55）得第167页/共192页(4)用单元集成法形成整体刚度矩阵K 由图11-26b，各杆的单元定位向量可写出如下：按照单元定位向量 ，将各单元 中的元素在K中定位，并与前阶段结果

43、累加，最后得到K如下：12341234第168页/共192页（5）结点荷载P 由图11-26a和b，P 直接写出如下：（6）解基本方程第169页/共192页（7）求各杆杆端力单元：单元：第170页/共192页单元：单元：第171页/共192页单元：第172页/共192页单元：各杆内力值标在图11-26c中桁架各杆旁边。第173页/共192页设横梁截面抗拉和抗弯刚度分别为EA和EI，且。又吊杆截面抗拉刚度2.组合结构计算组合结构时，先区分梁式杆和桁杆。对梁式杆，采用一般单元的单元刚度方程及相应的计算公式。对桁杆，采用桁架单元的单元刚度方程及相应的计算公式。例11-7试求图11-27所示组合结构的

44、内力。图11-27第174页/共192页解（1）单元和结点位移分量的统一编码如图 11-27中所示。横梁固定端的三个位移分量都为零，也用0编码。拉杆和在支座处的两个线位移分量均为零，也用0编码。单元和间刚结点A，编码为（1 2 3）。单元和间刚结点B，编码为（4 5 6），单元、和横梁的铰接点处，线位移不独立，因此应采用同码，分别（1 2）和（4 5）。所取的整体坐标系也示于图11-27中。（2）形成局部坐标系中的单元刚度矩阵第175页/共192页单元、为梁式杆，按一般单元的式（11-6）形成单元刚度矩阵单元、为桁架，按桁架单元的（11-54)形成单元刚度矩阵第176页/共192页(3)形成整

45、体坐标系中的单元刚度矩阵单元、：，所以单元：，由式（11-55）。第177页/共192页单元：由式（11-55），（1）（2）（3）（4）（1）（2）（3）（4）第178页/共192页（1）（2）（3）（4）（1）（2）（3）（4）(4)用单元集成法形成整体刚度矩阵K 图11-27中，各单元的 轴正方向用箭头标明。各杆单元定位向量可由图直接写出如下：第179页/共192页按照单元定位向量 ，将各单元 中的元素在K 中定位并与前阶段结果累加，最后得到K如下（集成时用 的关系）：123456123456（5）求等效结点荷载P首先，求单元固端约束力第180页/共192页只有单元有。因为无需转换坐标，

46、所以（1）（2）（3）（4）（5）（6）第181页/共192页123456按照单元定位向量 ，将 中的元素在P中定位，得（6）解基本方程第182页/共192页得（7）求各杆杆端力第183页/共192页第184页/共192页（8）作内力图内力图如图11-28所示。图11-28第185页/共192页1110小结矩阵位移法是新的计算工具（电子计算机）与传统力学原理（位移法）相结合的产物。矩阵位移法要与传统位移法对照起来学习，注意它们之间“原理上同源、作法上有别”的关系。位移法最便于实现计算过程的程序化。矩阵位移法（有限元位移法）是结构矩阵分析（有限元法）中占主导地位的方法。由式（10-49）可知位移

47、法基本方程的矩阵形式为（a）第186页/共192页这是一个非常普遍、非常简洁的方程。它既可用于分析梁、桁架等平面和空间结构，又可用于分析板、壳和弹性力学问题，具有普遍性。这样丰富的内容凝聚在三个符号、一个方程之中，形式简洁而具有高度的概括。矩阵位移法基本方程的建立，归结为两个问题：一是根据结构的几何和弹性性质建立整体刚度矩阵K，二是根据结构的受载情况形成整体荷载向量P。第187页/共192页()按 集成推导整体刚度矩阵K时，采用单元集成法，其推导过程为 ()坐标转换(a)定位(b)累加(b)第步，进行坐标转换，由局部坐标系的单元刚度矩阵 导出整体坐标系的单元刚度矩阵 ，为第步作好准备。第188

48、页/共192页第步，根据单元定位向量 ，依次由各单元的刚度矩阵 进行集成，得出整体刚度矩阵 。“集成”实际上包括将 的元素在 中“定位”以及将定在同一座位上的诸元素进行“累加”两个环节。为了便于理解，不妨对“定位”和“累加”分别说明其含义。为此，在 与 之间再引入单元贡献矩阵 ，如式（b）中实线所示。第189页/共192页在进行整体分析时，有“先处理”和“后处理”两种作法。本章采用的是先处理法，即在形成整体刚度矩阵时事先已根据结构的支承条件进行了处理。后处理法的特点是，先不考虑支承条件，按单元集成法得出原始的整体刚度矩阵，然后再引入支承条件，进行处理，得出整体刚度矩阵。关于后处理法可参阅有关书籍。第190页/共192页为了加深对矩阵位移法的理解，应当同时对计算程序有所了解并上机进行实践。本书在附录A中给出平面刚架源程序（FORTRAN77和Fortran90）,可供学习和上机实习时参考。第191页/共192页感谢您的观看！第192页/共192页