平面問题的有限元分析及三角形单元的应用.doc

平面问题的有限元分析及三角形单元的应用 第一节 概述 分析弹性力学平面问题时，最简单的单元式由三个结点组成的三角形单元。当用以分析平面应力问题时，可将其视为三角板；当用以分析平面应变问题时，则可式为三棱柱。各单元在结点处为铰结。图8-1所示位移悬臂梁离散为三角形单元的组合体 以矩阵形式列出弹性力学平面问题的基本量和基本方程。 谈形体所受体力分量可表示为 (8-1) 所受面力分量可表示为 （8-2） 体内任一点应力分量可表示为 （8-3） 任一点的应变分量可表示为 （8-4） 任一点的位移分量可表示为 （8-5） 弹性力学平面问题的几何方程的矩阵表达式为 （8-6） 平面应力问题的物理方程的矩阵表达式为 （8-7） 或简写成为 （8-8） 式中 （8-9） 称为弹性矩阵。 平面应变问题的物理方程也可写成式（8-8），但须将式（8-9）中的E换成，μ换成，因此得出 （8-10） 平衡微分方程及边界条件也可以用矩阵表示，但弹性力学有限元位移法中，通常用虚功方程代替平衡微分方程和应力边界条件。虚功方程的矩阵表达式为 （8-11） 式中：，表示虚位移； ，表示与虚位移相对应的虚应变。 为了便于计算，作用于弹性体上的体力和面力替换为作用在结点上的集中力，即等效结点荷载。设作用于各个结点上的外力分量用如下列阵来表示 与这些结点外力分量相对应得结点虚位移分量列阵为 则外力在虚位移上做的虚功为 如平面弹性体的厚度为t，该虚功除以t，即可得出单位厚度薄板上的外力虚功。于是，式（8-11）所示虚功方程可写成 （8-11） 虚功方程不仅仅应用于弹性力学，也可用于塑性力学。其应用条件是：只要变形体的全部外力和应力满足平衡方程；位移是微小的，并满足边界条件，位移与应变满足几何方程。所以，通常称为变形体虚功方程。 第二节 单元分析 图8-2所示为一个三角形单元。三个结点按逆时针顺序编号分别为i、j、m，结点坐标分别为。 图8-2 由于每个结点有两个位移分量，单元共有六个结点位移分量：，如图8-2a）所示，因此三角形单元的结点位移分量δe可表示为 （8-13） 与这六个结点位移分量相对应得结点力也有六个分量，如图8-2b)所示 （8-14） 在每个单元上，都可以把结点力用结点位移来表示，即建立如下关系式 （8-15） 式中ke称为单元刚度矩阵。寻求ke的过程称为单元分析。单元分析按如下步骤 位移函数 为了求单元内任一点（x，y）的位移，设该点的位移u、v为其坐标x、y的某种函数，单元有六个结点位移分量，在位移函数中取六个任意参数αi（i=1，2，…，6），并将位移函数取为线性函数，即 （8-16） 一般情况下，一个弹性变形体在外界作用下，内部点的位移变化比较复杂，不能用简单的线性函数描述。但是，当把弹性体离散为许多微笑单元时，在每一个单元内部有限小的局部内，各点位移可以用线性函数描述。式（8-16）可写成矩阵形式 （8-17） 为了求出内部结点位移f与结点位移δe之间的关系，需求出δe与α间的关系。降格结点坐标和位移代入式（8-16），可得 （a） （b） 三角形单元的面积为 （8-18） 求解方程组(a)得 （c） 求解方程组(b)得 (d) 式中，由下式计算 （i、j、m） 上式中的（i、j、m）表示脚标依次轮换，可写出计算aj、bj、cj以及am、bm、cm的另两组公式。 将式（c）和(d)代入（8-16）并展开，得到以结点位移表示的位移函数 （8-20） 式中，反映了单元的位移形态，故称为单元位移的形态函数或形函数。矩阵N称为形函数矩阵。 选取得位移函数是否合理，要看随着单元网格的逐步细分，有限元解是否逼近于精确解。为了保证收敛型所选择的单元位移函数应满足以下条件： 包含单元的刚体位移； 包含单元的常量应变； 保证相邻单元在公共边界处位移的连续性。 单元的应变和应力 选择了位移函数并以结点位移表示单元内点的位移后，重新写出平面问题的几何方程 （f） 由式（8-20）得 （g） 将式（g）代入式（f），并利用下式 （i、j、m） (h) 得单元应