JGLXchap13矩阵位移法.pptVIP

§9-1 概 述 前面介绍的力法、位移法和力矩分配法都是建立在手算基础上的传统计算超静定结构的方法，因而只能解决一些比较简单的超静定结构的计算问题。随着电子计算机的广泛使用，给以计算为特征的结构力学在计算方法上带来了很大的变革。传统计算方法已不能适应新的计算技术的要求,于是适合于电算的结构矩阵分析便得到了迅速的发展。它是以传统的位移法作为理论基础,力学分析中采用数学中的矩阵形式,计算工具采用计算机。因为采用矩阵运算,不仅使得公式形式紧凑,而且运算规律性强，很适合电算的要求,便于编制简单而通用性强的计算程序,故该方法得到了广泛应用。 杆系结构矩阵分析又叫杆系结构的有限元法,分为矩阵力法和矩阵位移法,亦称为柔度法和刚度法。由于矩阵位移法比矩阵力法更容易实现计算过程程序化,因而应用很广泛,故本章只讨论矩阵位移法。 §9-2 单元刚度矩阵 1. 一般单元杆端力和杆端位移的表示方法 图9-1所示刚架中的一等截面直杆单元e。设杆件除弯曲变形外，还有轴向变形。杆件两端各有三个位移分量(两个移动、一个转动)，杆件共有六个杆端位移分量，这是平面杆系结构单元的一般情况，故称为一般单元。单元的两端采用局部编码i和j点。 §9-3 单元刚度矩阵的坐标变换 在上节中，单元刚度矩阵是建立在杆件的局部坐标系中的。其目的是推导出的单元刚度矩阵形式最简单。如果从整体分析的 角度来考虑，对于整个结构，由于各杆轴方向不尽相同，因而各单元的局部坐标也不尽相同，很不统一。为了便于整体分析，在考虑整个结构的几何条件和平衡条件时，必须选定一个统一的坐标系，称为结构坐标系(或整体坐标系)。为了与局部坐标相区分，结构坐标系用xoy表示。 为了推导结构坐标系下的单元刚度矩阵 ，可采用坐标变换的方法，即把局部坐标系中建立的单元刚度矩阵 转换为结构坐标系中的 ，为此，首先讨论两种坐标系中单元杆端力的转换式，得到单元坐标转换矩阵；其次再讨论两种坐标系中单元刚度矩阵的转换式。 §9-4 结构的原始刚度矩阵 从本节开始进行结构的整体分析，即在单元分析的基础上，考虑各结点的几何条件及平衡条件，建立结构的刚度方程和结构刚度矩阵。现以图9-6(a)所示刚架为例来说明。 §9-5 支承条件的引入 §9-6 非结点荷载的处理 §9-7 矩阵位移法的计算步骤及示例 §9-8 几点补充说明 作业： 9-2、9-4、9-5、9-6 解：(1) 将单元、结点编号，确定坐标系，如图所示。 (3) 求出各单元在整体坐标系中的单刚，见例9-1。 (3) 将各单刚子块对号入座，形成结构原始刚度矩阵，见例9-1。 (4) 计算固端力、等效结点荷载及综合结点荷载。根据表9-3可知，各单元在局部坐标系中的固端力为 由式(9-40)，并将单元①的，单元②、③的代入计算，可得各单元在整体坐标系中的固端力为 由式(9-41)可求出结点2、3上的等效结点荷载为 再由式(9-42)求得综合结点荷载为 于是可得结点外力列向量为 这里需要指出，对于支座结点1，4同样可按式(9-41)和式(9-42)算出等效结点荷载和综合结点荷载。但由于[F1]、[F4]是综合结点荷载与支座反力的代数和，而其中的反力仍是未知的，又由于在引入支承条件时，未知力[F1]、[F4]将被删去或被修改，故在此可不必计算支座结点1、4的等效结点荷载及综合结点荷载。 (5) 引入支承条件，修改原始刚度方程。结构的原始刚度方程为 见p.256 结点1和4为固定端，故已知 在原始刚度矩阵中删去与上述零位移相应的行和列，同时在结点位移列向量和结点外力列向量中删去相应的行，得结构的刚度方程为 (6) 解方程。求得未知结点位移为 (7) 计算各单元杆端力。按式(9-45)计算。 单元①： 单元②： 单元③： (8) 绘内力图。绘刚架的内力图时，内力正负号规定仍与第三章中的规定相同，而上述计算出的杆端内力正负号是按本章中的规定求得的。故画 结构坐标系下的单元刚度矩阵可按式(9-18)来计算，其四个子块为 （9-26）  图9-6（a） 图9-6 首先对各单元及结点进行编号。用①、②、…表示单元编号；用1、2、…表示结点编号，这里支座也视为结点。其次，选取结构坐标系和各单元的局部坐标系如图9-6(b)所示。各单元的始末两端i、j的结点号码如表9-1表示，则按式(9-23)表示的各单元刚度矩阵的四个子块应为： （a） 表9-1 各单元始末端的结点编号 1 2 2 3 3 4 ① ② ③ 始末端结点号 i j 单元 在平面刚架中，每个结点有两个