hzf第6讲——用有限单元法解平面问题1.ppt

§ 6-3 三角形常应变单元的位移模式和解答的收敛性 一、离散化 在运用有限单元法分析弹性力学平面问题时，第一步就是要对弹性体进行离散化，把一个连续的弹性体变换为一个离散的结构物。对于平面问题，三角形单元是最简单、也是最常用的单元，在平面应力问题中，单元为三角形板，而在平面应变问题中，则是三棱柱。 假设采用三角形单元，把弹性体划分为有限个互不重叠的三角形。这些三角形在其顶点（即节点）处互相连接，组成一个单元集合体，以替代原来的弹性体。同时，将所有作用在单元上的载荷（包括集中载荷、表面载荷和体积载荷），都按虚功等效的原则移置到节点上，成为等效节点载荷。由此便得到了平面问题的有限元计算模型，如图6-1所示。 ?图6-1 弹性体和有限元计算模型 ?图6--2 平面三角形单元 二、位 移 首先，我们来分析一下三角形单元的力学特性，即建立以单元节点位移表示单元内各点位移的关系式。设单元e的节点编号为i、j、m，如图6-2所示。由弹性力学平面问题可知，每个节点在其单元平面内的位移可以有两个分量，所以整个三角形单元将有六个节点位移分量，即六个自由度。用列阵可表示为： 其中的子矩阵 （i，j，m 轮换） (a) 式中 ui、vi 是节点i在x轴和y轴方向的位移。 (6-7) 在有限单元法中，虽然是用离散化模型来代替原来的连续体，但每一个单元体仍是一个弹性体，其内部依然是符合弹性力学基本假设的，弹性力学的基本方程在每个单元内部都适用。 从弹性力学平面问题的解析解法中可知，如果弹性体内的位移分量函数已知，则应变分量和应力分量也就确定了。但是，如果只知道弹性体中某几个点的位移分量的值，那么就不能直接求得应变分量和应力分量。 　　因此，在进行有限元分析时，必须先假定一个位移模式。由于在弹性体内，各点的位移变化情况非常复杂，很难在整个弹性体内选取一个恰当的位移函数来表示位移的复杂变化，但是如果将整个区域分割成许多小单元，那么在每个单元的局部范围内就可以采用比较简单的函数来近似地表示单元的真实位移，将各单元的位移式连接起来，便可近似地表示整个区域的真实位移函数。 这种化繁为简、联合局部逼近整体的思想，正是有限单元法的绝妙之处。 基于上述思想，我们可以选择一个单元位移模式，单元内各点的位移可按此位移模式由单元节点位移通过插值而获得。线性函数是一种最简单的单元位移模式，故设 (b) 式中 ?1、?2、…?6是待定常数。因三角形单元共有六个自由度，且位移函数u、v在三个节点处的数值应该等于这些点处的位移分量的数值。假设节点i、j、m的坐标分别为（xi , yi ）、（xj , yj ）、（xm , ym ），代入 (b) 式，得： (c) 由 (c) 式左边的三个方程可以求得 (d) 其中 (6-8) 从解析几何可知，式中的 ? 就是三角形i、j、m的面积。为保证求得的面积为正值，节点i、j、m的编排次序必须是逆时针方向，如图6-2所示。 ?图6-2 平面三角形单元 将 (d) 式代入 (b) 式的第一式，经整理后得到 (e) 其中 同理可得 若令 这样，位移模式 (e) 和 (f) 就可以写为 （i , j , m轮换） (6-10) （i , j , m轮换） (6-9) (f) 式中 I是二阶单位矩阵；Ni 、Nj 、Nm 是坐标的函数，它们反映了单元的位移状态，所以一般称之为形状函数，简称形函数。矩阵 [N] 叫做形函数矩阵。 (6-11) 也可写成矩阵形式 (6-12) 形函数的性质 在前面的分析中，提出了形函数的概念，即 其中 （i , j , m轮换） 现在我们来讨论一下形函数所具有的一些性质。根据行列式的性质：行列式的任一行（或列）的元素与其相应的代数余子式的乘积之和等于行列式的值，而任一行（或列）的元素与其他行（或列）对应元素的代数余子式乘积之和为零，并注意到（6-9）式中的常数ai 、bi 、ci ，aj 、bj 、 cj 和am 、bm 、cm 分别是行列式2的第一行、第二行和第三行各元素的代数余子式，我们有 ⒈ 形函数在各单元节点上的值，具有“本点是1、它点为零”的性质，即 在节点i上， 在节点j、m上， (a) (b) (c) 类似地有 (d) ⒉ 在单元的任一节点上，三个形函数之和等于1，即 (e) 简记为 (6-22) 这说明，三个形函数中只有二个是独立的。 ?⒊三角形单元任意一条边上的形函数，仅与该边的两端节点坐标有关、而与其它节点坐标无关。例如，在i j 边上，有 (6-23) 事实上，因i j 边的直线方程方程为 (f) 代入（6-10）式