平面杆系结构的单元分析局部坐标系下的单元刚.ppt

首先对各单元及结点进行编号。用①、②、…表示单元编号；用1、2、…表示结点编号，这里支座也视为结点。其次，选取结构坐标系和各单元的局部坐标系如图9-6(b)所示。各单元的始末两端i、j的结点号码如表9-1表示，则按式(9-23)表示的各单元刚度矩阵的四个子块应为： （a） 表9-1 各单元始末端的结点编号 1 2 2 3 3 4 ① ② ③ 始末端结点号 i j 单元 在平面刚架中，每个结点有两个线位移和一个角位移。此刚架有四个结点，共有12个结点位移分量，按一定顺序排列成一列阵，称为结点位移列向量，即 式中 其中这里，[Δi]表示结点i的位移列向量，、和分别为结点i沿结构坐标系x、y轴的线位移和角位移，它们分别以沿坐标轴的正向和顺时针方向为正。 设刚架只受到结点荷载(非结点荷载可等效为结点荷载，见§9-6)，则与结点位移列向量相应的结点外力(包括荷载和反力)列向量为 ，式中， 这里，[Fi]代表结点i的外力列向量，、和分别为作用于结点i的沿x、y方向的外力和外力偶，它们的正负号规定与相应的结点位移相同。在结点2、3处，结点外力[F2]、[F3]就是结点荷载，它们通常是给定的。而在结点1、4上，当没有给定结点荷载时，结点外力{F1}、[F4]就是支座反力；当支座处还有给定的荷载作用时，则应为结点荷载与支座反力的代数和。下面考虑结构的平衡条件和变形条件。各 单元和各结点的隔离体图如图9-6(c)所示。图中各单元上的杆端力均是沿着结构坐标系的正向作用的。在前面单元分析中，已经保证了各单元本身的平衡和变形连续，因此只需考虑各单元联结处，即各结点处的平衡条件和变形连续条件。现以结点2为例，由平衡条件ΣFx=0、ΣFy=0和ΣM=0，可得 写成矩阵形式 上式左边为结点2的荷载列向量[F2]，右边两列分别为单元①和单元②在2端的杆端力列向量[F2]①和[F2]②，故上式可简写为 [F2]=[F2]①+[F2]② (b) 根据式(9-22)，上述杆端力列向量可用杆端位移列向量来表示 (c) 再根据结点处的变形连续条件，应该有 (d) 将式(c)和(d)代入(b)，则得到以结点位移表示的结点2的平衡方程 (e) 同理，对结点1、3、4都可建立类似的平衡方程。将所有四个结点的方程汇集在一起， 就有 (9-27) 写成矩阵形式则为 (9-28) 上式便是用结点位移表示的所有结点的平衡方程，它表明了结点外力与结点位移之间的关系，通常称为结构的原始刚度方程，“原始”之意是指尚未引入支承条件。上式可简写为 [F]=〔K〕[Δ] (12-29) 式中 称为结构的原始刚度矩阵，也称结构的总刚度矩阵 (简称总刚)。它的每个子块都是3×3阶方阵，故 为12×12阶方程。其中的每一个元素的物理意义就是当其所在列对应的结点位移分量等于1(其余结点位移分量均为零)时，其所在行对应的结点外力分量的数值。 结构的原始刚度矩阵 具有如下性质：  (1) 对称性。原始刚度矩阵 是一个对称方阵，这可由反力互等定理得知。 (2) 奇异性。原始刚度矩阵是奇异的，其逆阵不存在。这是由于建立方程式(9-29)时，没有考虑结构的约束条件，结构还可以有任意刚体位移，结点位移的解答不是唯一的。故还不能由式(9-29)来求结点位移，只能在引入支承条件，对结构的原始刚度方程进行修改后，才能求解未知的结点位移，这将在下一节中讨论。 下面讨论如何由单元刚度矩阵直接形成结构原始刚度矩阵。 考察前面式(a)及式(9-28)可知，结构的原始刚度矩阵是由每个单元刚度矩阵的四个子块按其两个下标号码送入结构刚度矩阵中的相应位置上而形成的。也就是将各单元子块“对号入座”即形成总刚。以单元③的四个子块为例，其入座位置如图9-7所示。一般而言，某单刚子块 就应被送入总刚(以子块形式 形式表示)中第i行第j列的位置上去。这种利用结构坐标系中的单元刚度矩阵子块对号入座直接形成总刚的方法，称为直接刚度法。 图9-7 在对号入座时，具有相同下标的各单刚子块，在总刚中被送到同一位置上，各单刚子块要进行叠加，而没有单刚子块送入的位置上则为零子块。位于主对角线上的子块称为主子块；其余子块称为副子块。同交于一个结点上的各杆件称为该结点的相关单元；两个结点之间有杆件直接相联者称为相关结点。则由单刚子块对号入座形成总刚