平面三角形单元

第八章平面问题的有限元分析及三角形单元的应用第一节概述分析弹性力学平面问题时，最简单的单元式由三个结点组成的三角形单元。当用以分析平面应力问题时，可将其视为三角板；当用以分析平面应变问题时，则可式为三棱柱。各单元在结点处为铰结。图8-1所示位移悬臂梁离散为三角形单元的组合体以矩阵形式列出弹性力学平面问题的基本量和基本方程。谈形体所受体力分量可表示为Tyxyxppppp(8-1)所受面力分量可表示为Tyxyxppppp（8-2）体内任一点应力分量可表示为Txyyx（8-3）任一点的应变分量可表示为Txyyx（8-4）任一点的位移分量可表示为Tvu（8-5）弹性力学平面问题的几何方程的矩阵表达式为xuyvyvxuxyyx（8-6）平面应力问题的物理方程的矩阵表达式为xyyxxyyxE2100010112（8-7）或简写成为D（8-8）式中2100010112ED（8-9）称为弹性矩阵。平面应变问题的物理方程也可写成式（8-8），但须将式（8-9）中的E换成21E，μ换成21，因此得出)1(22100011011)21)(1()1(22ED（8-10）平衡微分方程及边界条件也可以用矩阵表示，但弹性力学有限元位移法中，通常用虚功方程代替平衡微分方程和应力边界条件。虚功方程的矩阵表达式为tdxdytdspfptdxdyfTT（8-11）式中：Tvuf，表示虚位移；Txyxx，表示与虚位移相对应的虚应变。为了便于计算，作用于弹性体上的体力和面力替换为作用在结点上的集中力，即等效结点荷载。设作用于各个结点上的外力分量用如下列阵来表示TnnVUVUVUF2211与这些结点外力分量相对应得结点虚位移分量列阵为Tnnvuvuvu2211则外力在虚位移上做的虚功为FvVuUvVuUvVuUTnnnn22221111如平面弹性体的厚度为t，该虚功除以t，即可得出单位厚度薄板上的外力虚功。于是，式（8-11）所示虚功方程可写成tdxdyFTT（8-11）虚功方程不仅仅应用于弹性力学，也可用于塑性力学。其应用条件是：只要变形体的全部外力和应力满足平衡方程；位移是微小的，并满足边界条件，位移与应变满足几何方程。所以，通常称为变形体虚功方程。第二节单元分析图8-2所示为一个三角形单元。三个结点按逆时针顺序编号分别为i、j、m，结点坐标分别为),(),(),(mmjjiiyxmyxjyxi、、。图8-2由于每个结点有两个位移分量，单元共有六个结点位移分量：mmjjiivuvuvu、、、、、，如图8-2a）所示，因此三角形单元的结点位移分量δe可表示为Tmmjjiievuvuvu（8-13）与这六个结点位移分量相对应得结点力也有六个分量，如图8-2b)所示TmmjjiieVUVUVUF（8-14）在每个单元上，都可以把结点力用结点位移来表示，即建立如下关系式eeekF（8-15）式中ke称为单元刚度矩阵。寻求ke的过程称为单元分析。单元分析按如下步骤一、位移函数为了求单元内任一点（x，y）的位移，设该点的位移u、v为其坐标x、y的某种函数，单元有六个结点位移分量，在位移函数中取六个任意参数αi（i=1，2，…，6），并将位移函数取为线性函数，即yxyxvyxyxu654321),(),(（8-16）一般情况下，一个弹性变形体在外界作用下，内部点的位移变化比较复杂，不能用简单结点位移内部各点位移应变应力结点力ke的线性函数描述。但是，当把弹性体离散为许多微笑单元时，在每一个单元内部有限小的局部内，各点位移可以用线性函数描述。式（8-16）可写成矩阵形式65432110000001yxyxvuf（8-17）为了求出内部结点位移f与结点位移δe之间的关系，需求出δe与α间的关系。降格结点坐标和位移代入式（8-16），可得321111mmjjiimjiyxyxyxuuu（a）654111mmjjiimjiyxyxyxvvv（b）三角形单元的面积为mmjjiiyxyxyxA11121（8-18）求解方程组(a)得mjimjimjimjiuuucccbbbaaaA21321（c）求解方程组(b)得mjimjimjimjivvvcccbbbaaaA21654(d)式中，mmmiiicbacba、、、、、、由下式计算mjimjijmmjixxcyybyxyxa（i、j、m）上式中的（i、j、m）表示脚标依次轮换，可写出计算aj、bj、cj以及am、bm、cm的另两组公式。将式（c）和(d)代入（8-16）并展开，得到以结点位移表示的位移函数jijijimjimjivuvuvuyxNyxNyxNyxNyxNyxNyxvyxu),(0),(0),(00),(0),(0),(),(),(（8-20）式中，mjiNNN、、反映了单元的位移形态，故称为单元位移的形态函数或形函数。矩阵N称为形函数矩阵。选取得位移函数是否合理，要看随着单元网格的逐步细分，有限元解是否逼近于精确解。为了保证收敛型所选择的单元位移函数应满足以下条件：（1）包含单元的刚体位移；（2）包含单元的常量应变；（3）保证相邻单元在公共边界处位移的连续性。二、单元的应变和应力选择了位移函数并以结点位移表示单元内点的位移后，重新写出平面问题的几何方程vuxyyxxyyx00（f）由式（8-20）得mmjjiimmijiivNvNvNvuNuNuNu（g）将式（g）代入式（f），并利用下式AcyNAbxNiiii22（i、j、m）(h)得单元应变mmjjiimmjjiimjimjixyyxvuvuvubcbcbccccbbbA00000021（8-24）或简写成eBδε（8-25）式中mmjjiimjimjibcbcbccccbbbAB00000021（8-26）式（8-25）就是由结点位移求应变的转换式，其转换矩阵B称为几何矩阵。将（8-25）代入平面问题的物理方程式（8-8）有eeDBD（8-27）或写成eeS（8-28）S=DB(8-29)称为应力矩阵。三、单元刚度矩阵在有限单元法中，常利用虚功方程代替平衡方程。图8-3a)所示为三角形单元的实际力系，其结点力为Fe，应力为σ；图8-3b)所示为单元虚位移状态，其结点位移为δ*e,应变为ε*。利用式（8-12）可得图8-3tdxdyFTeeT（8-30）式中：为单元虚应变。为单元结点虚位移；e由式（8-25）可知eB因此TeTTB将此公式代入式（8-30），由于e中的元素是常量，公式右边的eT可以提到积分号的前面，得tdxdyBFTeTeeT由于虚位移e是任意的，则tdxdyBFe因为B和σ都是常量矩阵，并且积分Adxdy，所以tABFTe（8-31）利用（8-27），可得tADBBFeTe（8-32）令DBtABkT（8-33）则式（8-32）就变成式（8-15），即eeekF单元刚度矩阵ek为一个6×6矩阵，它时单元结点位移与单元结点力之间的转换矩阵，具有以下性质：（1）ek示对称矩阵，其元素jiijkk；（2）ek是奇异矩阵，由它的元素组成的行列式等于零，即它不存在逆矩阵；（3）ek具有分快性质。第三节等效结点荷载为简化各单元得受力情况，便于分析计算，应将单元所受各种载荷向结点移置，化为结点荷载，荷载的移置应安静力等效原则进行。静力等效是令原来的荷载与移置后的荷载在任意虚位移上的虚功相等。一、集中荷载图8-5a）所式为单元内的任意一点M受到集中荷载P的作用，沿x、y方向的分量分别为Px、Py，用矩阵表示为TyxPPP。设移置到该单元结点上的等效荷载列阵为TmmjjiieYXYXYXR图8-5所示单元发生虚位移，其中单元内部任意点的虚位移为Tvuf设单元各结点虚位移为Tmmjjiievuvuvu由式（8-22），则eNf（a）图8-5如果设Re为单元得等效结点荷载，则Re在结点虚位移上所做的虚功应与原来集中荷载在其作用点的虚位移上做的虚功相等，即PfRTeeT将式(a)代如上式，得PNRTeTeeT由于虚位移可以式任意的，因此PNRTe（8-36）或写成yiixiiPNYPNX（i、j、m）二、分布体力设单元受分布体力的作用Tyxppp将微分体积tdxdy上的体力当作集中荷载P，利用式（8-36）的积分得出ptdxdyNRTe（8-38）如果单元上作用的分布体力为自重，γ为材料的重度，即Tp0那么利用式（8-38）得到等效结点荷载列阵为TeAtR1010103（8-39）三、分布面力设单元在边界ij上首有分布面力的作用Tyxppp将微分面积tds上的面力当作集中荷载P，利用式（3-36）的积分可以得出tdspNRTe（8-40）第四节整体刚度矩阵有限元求解弹性力学问题也要用结点平衡方程求解作为基本为质量的整体结点位移列阵。在求出各单元的单元刚度矩阵ke和解点荷载列阵后，就可以用几何的方法建立其平面弹性体的整体刚度矩阵和整体平衡方程。如果弹性体划分为m各单元，n各结点，则有121222nnnnRK（8-44）式中：整体刚度矩阵K为2n×2n阶方程；δ、R分别为整体结点位移列阵和整体结点荷载列阵，都是2n×1阶列阵。第六节单元网格的划分和计算成果的整理一、单元网格的划分划分单元网格时，要综合考虑单元的数目和划分的合理性，注意以下方面：1．合理安排单元网格的疏密分布结构的不同部位疏密不同，边界比较曲折的部位，网格可以密一些，比较平滑的部位，可以疏一些，对于应力应变相对小的部位疏一些。在保证计算精度的前提下，减少单元划分数目。2．对称性的利用当对称结构上有对称荷载或反对成荷载作用是，可以利用对称性取半结构或1/4结构。3．单元形状的合理性三角形单元的三个边长较为接近最好，可以避免由于周围应力长分布不均而引起较大的计算误差。4．不同材料截面处单元划分当计算对象为不同材料组成时，应以材料性质发生变化的不同材料界面作为单元的边界。