第三章弹性力学有限元法§33 单元分析.ppt

有限元模型 §3.4 整体分析 显然，单元分析得出的仅仅是局部的信息，各个单元靠节点连接起来组成整体，因而必须从全局进行分析。就是说，将各个单元的方程(单元刚度矩阵)按照保持节点处位移连续性的方式组合起来，就可得到整个物体的平衡方程(整体刚度矩阵)，再按照给定的位移边界条件修改这些方程，使平衡方程组有解。 第三章 弹性力学有限元法 §3.4 整体分析 1．整体分析的基本思路 “积零为整”是由单元方程建立整体方程的集合过程，应遵循以下原则：根据节点的平衡条件，相应的刚度和载荷必须相加。这个过程包括由各个单元刚度矩阵 集合成整体的刚度矩陈[K]，以及由各单元节点力向量 集合成总的载荷向量 。本节运用的集合法是最常用的直接刚度法。集合所依据的原理是：各单元在公共节点处相互连接，因而要求在公共节点处连接的单元在该接点上的位移相同。集合而成的代数方程组反映了整体的平衡关系 第三章 弹性力学有限元法 §3.4 整体分析 1．整体分析的基本思路 其中整体刚度矩阵[K]仍是奇异的。为使上述方程组有唯一解，必须利用几何边界条件——在物体边界处给定某些位移约束条件，给上述方程组以适当修正。 现以下图所示的平面问题为例，进行整体分析。 第三章 弹性力学有限元法 §3.4 整体分析 2．整体刚度矩阵的集合 在平面问题中，每个节点的位移和节点力的分量，都取为沿其x轴和y轴的两个分量，对于每个子单元而言，均含有三个节点,因而单元的位移向量 和节点力向量 都包含3×2=6个分量，所以，单元刚度矩阵 是6×6阶矩阵．对此例的整体而言，共含6个节点,因而单元的位移向量 和节点力 都包括6×2=12个分量，所以，整体刚度矩阵[K]是12×12阶矩阵。 为了将单元刚度矩阵叠加到整体刚度矩阵中,需要将单元刚度矩阵中的局部节点编号转换成结构的统一编号,这样,各单元刚度阵中的各子块的下标就表示出它在整体矩阵中的位置。 第三章 弹性力学有限元法 §3.4 整体分析 对上图所示的图形的第一单元,取节点号为l=1,m=2,n=3,则单元刚度阵为: 第二单元取节点号为l=2,m=5,n=3,则单元刚度阵为 第三章 弹性力学有限元法 §3.4 整体分析 第三、四单元l,m,n,分别取2,4,5与3,5,6,则单元刚度阵为: 由于单元刚阵都是在统一的直角坐标系中建立的,故可以直接迭加。首先将 按其元素的下标分块加到整体刚度矩阵中,得 第三章 弹性力学有限元法 §3.4 整体分析 再将 加入,得 第三章 弹性力学有限元法 §3.4 整体分析 用同样的方法将 及 加入,得整体刚度矩阵 第三章 弹性力学有限元法 §3.4 整体分析 3.整体载何向量的建立 (l)节点处载荷向量的处理：可以按照上述建立整体刚度矩阵的集合方法，去建立整体载荷向量 。但是，实际上往往根据已知条件直接给出各节点的载荷。一般而言，如果节点在某方向上有支承，那么在这方向上的节点力为未知量；反之，如果节点在某方向上没有约束及外力，则在这方向上的节点力为零，例如，图(b)中的整体载荷向量为 第三章 弹性力学有限元法 §3.4 整体分析 3.整体载何向量的建立 (2)对非节点载荷的处理 如果结构上受有不直接作用于节点上的载荷，则要根据静力等效的原则，向节点移置。静力等效原则是指原载荷和移置到节点的载荷在任何虚位移上所作的虚功相等。如下图所示的单元123的边12长为l,在边12上受有沿x 方向的载荷F，作用点距节 点l和2的距离分别为l1和l2， 可将载荷F移置到节点l和 2上，得 第三章 弹性力学有限元法 §3.4 整体分析 3.整体载何向量的建立 (2)对非节点载荷的处理 又如，设该单元的边12上受有沿x方向r的均布载荷q,如图b所示，可将分布载荷转换成节点l和2上的集中载荷，得 第三章 弹性力学有限元法 §3.5 约束条件的处理 前面已经讨论了引人支承条件并对方程组 进行修正的必要性，即不进行这样的修正，方程组便无法得到唯一解。 (1)当节点n带有水平支承，使水平位移 ，则应对总体平衡方程作如下修正：在刚度矩阵[K]中