5平面分析整体刚度矩阵..pptVIP

对于同样的有限元单元网格，按照图（a）的结点编码，最大的半带宽为 按照图（b）的结点编码，最大的半带宽为B 按照图（c）的结点编码，最大的半带宽为B (a) (b) (c) * * 平面问题有限元分析-总刚 2×(10－6)＝10 5.2 整体刚度矩阵的特点 B = 2×(10－4)=14 2×(10－2)＝18 4. 是奇异矩阵，在排除刚体位移后，它是正定阵 无约束的弹性体（或结构物）的整体刚度矩阵是奇异的，不存在逆矩阵，即关于位移的解不唯一。这是因为弹性体在外力的作用下处于平衡，外力的分量应该满足三个静力平衡方程。这反映在整体刚度矩阵中就意味着存在三个线性相关的行或列，所以是个奇异阵，不存在逆矩阵。例如：设弹性体在外力的作用下处于平衡，这时相应的解为 ，然后在给予弹性体以刚体位移而相应的节点位移 ，这时， 仍是问题的解，因为刚体位移不会破坏平衡。 注：当排除刚体位移后，整体刚度矩阵是正定矩阵。 * * 平面问题有限元分析-总刚 5.2 整体刚度矩阵的特点 前文已经提到在排除刚体位移后，整体刚度矩阵是正定的，方程才可求得唯一解。排除刚体位移可以通过引入边界约束条件来实现，这里介绍两种比较简单的引入已知位移的方法。 1．代入法 2．乘大数法 * * 平面问题有限元分析-总刚 5.3 边界条件 1．代入法 该方法保持方程组仍为2n×2n系统，仅对整体刚度矩阵 和整体载荷列阵进行修正。下面以一个只有四个方程的简单例子加以说明，方程如下 假定系统中节点位移 、 ，则当引入这些节点的已知位移之后，方程就变成 * * 平面问题有限元分析-总刚 5.3 边界条件 若 ，则 然后，用这组维数不变的方程来求解所有的节点位移。显然，其解答仍为原方程的解答。 在手算时，可直接将零位移约束所对应的整体刚度矩阵中的行和列直接划去，使得整体刚阵的维数变小，更便于手算。 * * 平面问题有限元分析-总刚 5.3 边界条件 2．乘大数法 将 中与指定的节点位移相对应的主对角元素乘上一个大数，同时将 中的对应元素换成指定的节点位移值、该大数与节点位移相对应的主对角元素三者的乘积。若把此方法用于上面的例子，则方程就变成 该方程组的第一个方程为 解得 ，这种方法就是使 中相应行的修正项远大于非修正项。 * * 平面问题有限元分析-总刚 5.3 边界条件 1．代入法 2．乘大数法 在以上的两种方法中，代入法接近人工解法，虽然该方法比较直观，但该方法对刚度矩阵改变较多，程序效率不高。乘大数法对刚度矩阵改变较少，工作量较小，但相乘的“大数”若取得过大，求解时会发生“溢出”、若取得太小则会引起较大的误差。 * * 平面问题有限元分析-总刚 5.3 边界条件 对于三节点三角形单元，单元内各点的应力值相等，算出的应力一般作为单元形心处的应力。由于单元应力为常数，整个结构的应力场呈阶梯状，在单元之间不连续。而工程上往往更加关心边界和节点上的受力情况，因此，必须对所得到的应力再次进行处理，得到更加合理的应力场，并得到所需点上的应力值。 这里介绍两种简单的方法，一种方法称为节点平均法，即把环绕某一节点的各单元的应力加以平均作为该节点的应力值。例如图中节点3的应力为 * * 平面问题有限元分析-总刚 5.4 计算结果整理 为了使通过这样平均得来的应力比较接近实际情况，要求环绕节点的各单元尺寸不应相差太大。这种做法，对内节点比较好，对边界点则可能很差。 因此，边界节点处的应力不宜直接由单元应力平均来获得，而要根据内节点的应力构造插值函数推算出来。例如图中边界点1的应力，可以先用节点平均法求得节点2、3、4处的应力，在构造相应的插值函数推算边界点1的应力，如常用的抛物线插值公式如下 5.4 计算结果整理 * * 平面问题有限元分析-总刚 另一种方法称为单元平均法，即把两相邻单元的应力加以平均，用以表示公共边界中点的应力。为了由这样平均所得到的应力具有较好的精度，两相邻单元的面积不应相差太大。如图中单元②和③边界的中点处的应力为 在不同的有限元软件中均具有各自的