整体分析及总体刚度矩阵的性质.pptVIP

整体分析及总体刚度矩阵的性质 整体分析 整体分析 整体分析 整体刚度矩阵的形式 整体刚度矩阵的形式 整体刚度矩阵的形式 2、整体刚度矩阵的集成方法 具体集成方法是：先对每个单元求出单元刚度矩阵 ，然后将其中的每个子块 送到结构刚度矩阵中的对应位置上去，进行迭加之后即得出结构刚度矩阵[K]的子块，从而得出结构刚度矩阵[K]。 关键是如何找出 中的子块在[K]中的对应位置。这需要了解单元中的节点编码与结构中的节点编码之间的对应关系。 整体刚度矩阵的形式 整体刚度矩阵的形式 整体刚度矩阵的形式 用同样的方法可得出其他单元的扩大的单元刚度方程: 据节点力平衡，各个单元相应节点力叠加： 整理可得，整体平衡方程： 整体刚度矩阵的形式 2. 整体刚度矩阵的特点 整体刚度矩阵的特点 整体刚度矩阵的特点 整体刚度矩阵的特点 整体刚度矩阵的特点 整体刚度矩阵的特点 * 图示结构的网格共有四个单元和六个节点。在节点1、4、6共有四个支杆支承。结构的载荷已经转移为结点载荷。 整体分析的四个步骤： 1、建立整体刚度矩阵； 2、根据支承条件修改整体刚度矩阵； 3、解方程组，求节点位移； 4、根据节点位移求出应力。 单元分析得出单元刚度矩阵，下面，将各单元组合成结构，进行整体分析。 1、建立整体刚度矩阵(也叫作结构刚度矩阵) 上图中的结构有六个节点，共有12个节点位移分量和12个节点力分量。由结构的节点位移向量求结构的节点力向量时，转换关系为： 分块形式为： 其中子向量 和 都是二阶向量，子矩阵 是二行二列矩阵。整体刚度矩阵[K]是12*12阶矩阵。 2、根据支承条件修改整体刚度矩阵。 建立整体刚度矩阵时，每个节点的位移当作未知量看待，没有考虑具体的支承情况，因此进行整体分析时还要针对支承条件加以处理。 在上图的结构中，支承条件共有四个，即在节点1、4、6的四个支杆处相应位移已知为零： 建立节点平衡方程时，应根据上述边界条件进行处理。 3、解方程组，求出节点位移。 通常采用消元法和迭代法两种方法。 4、根据节点位移求出应力。 整体刚度矩阵 是单元刚度矩阵 的集成。 1、刚度集成法的物理概念： 刚度矩阵中的元素，即由节点作单位位移时引起的节点力。在单元刚阵 中， 表示j节点单位位移，其他节点位移为零时，单元e在i节点引起的节点力；类似，在整体刚阵中， 表示j节点单位位移，其他节点位移为零时，整体结构在i节点引起的节点力（由于结构已被离散为一系列单元，即所有与i、j节点相关的单元在i节点引起的节点力之和）。 如上图结构，计算 时，与节点2和3相关的单元有单元①和③，当节点3发生单位位移时，相关单元①和③同时在节点2引起节点力，将相关单元在节点2的节点力相加，就得出结构在节点2的节点力 。由此看出，结构的刚度系数是相关单元的刚度系数的集成，结构刚度矩阵中的子块是相关单元的对应子块的集成。 2、刚度矩阵的集成规则： 1）在整体离散结构变形后，应保证各单元在节点处仍然协调地相互连接，即在该节点处所有单元在该节点上有相同位移， 2）整体离散结构各节点应满足平衡条件。即环绕每个节点的所有单元作用其上的节点力之和应等于作用于该节点上的节点载荷Ri， 1 2 i 3 4 1 2 i Ri 3 4 结构中的节点编码称为节点的总码，各个单元的三个节点又按逆时针方向编为i,j,m,称为节点的局部码。 单元刚度矩阵中的子块是按节点的局部码排列的，而结构刚度矩阵中的子块是按节点的总码排列的。因此，在单元刚度矩阵中，把节点的局部码换成总码，并把其中的子块按照总码次序重新排列。 以单元②为例，局部码i,j,m对应于总码5,2,4，因此 子块按照总码重新排列后，得出扩大矩阵 为： 而相应的单元刚度方程为（或节点力表达式）： 整体平衡方程： 1）其中[K]为将各单元的扩大矩阵迭加所得出的结构刚度矩阵： 集成包含搬家和迭加两个环节： A、将单元刚度矩阵 中的子块搬家，得出单元的扩大刚度矩阵 。 B、将各单元的扩大刚度矩阵 迭加，得出结构刚度矩阵[K]。 2） 为节点载荷向量， 为节点位移向量。 在有限元法中，整体刚度矩阵的阶数通常是很高的，在解算时常遇到矩阵阶数高和存贮容