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

1、.羀羈芄螁蚀膃膀螀螂羆蒈蝿袅膂蒄螈肇羅莀螇螇芀芆螆衿肃薅螅羁芈蒁螅肃肁莇袄螃芇芃蒀袅聿腿葿羈芅薇蒈螇肈蒃蒈袀莃荿蒇羂膆芅蒆肄罿薄蒅螄膄蒀薄袆羇莆薃罿膃节薂蚈羅膈薂袁膁薆薁羃肄蒂薀肅艿莈蕿螅肂芄薈袇芈膀蚇罿肀葿蚆虿芆莅蚆螁聿芁蚅羄莄芇蚄肆膇薆蚃螆羀蒂蚂袈膅莈蚁羀羈芄螁蚀膃膀螀螂羆蒈蝿袅膂蒄螈肇羅莀螇螇芀芆螆衿肃薅螅羁芈蒁螅肃肁莇袄螃芇芃蒀袅聿腿葿羈芅薇蒈螇肈蒃蒈袀莃荿蒇羂膆芅蒆肄罿薄蒅螄膄蒀薄袆羇莆薃罿膃节薂蚈羅膈薂袁膁薆薁羃肄蒂薀肅艿莈蕿螅肂芄薈袇芈膀蚇罿肀葿蚆虿芆莅蚆螁聿芁蚅羄莄芇蚄肆膇薆蚃螆羀蒂蚂袈膅莈蚁羀羈芄螁蚀膃膀螀螂羆蒈蝿袅膂蒄螈肇羅莀螇螇芀芆螆衿肃薅螅羁芈蒁螅肃肁莇袄螃芇芃蒀袅聿腿

2、葿羈芅薇蒈螇肈蒃蒈袀莃荿蒇羂膆芅蒆肄罿薄蒅螄膄蒀薄袆羇莆薃罿膃节薂蚈羅膈薂袁膁薆薁羃肄蒂薀肅艿莈蕿螅肂芄薈袇芈膀蚇罿肀葿蚆虿芆莅蚆螁聿芁蚅羄莄芇蚄肆膇薆蚃螆羀蒂蚂袈膅莈蚁羀羈芄螁蚀膃膀螀螂羆蒈蝿袅膂蒄螈肇羅莀螇螇芀芆螆衿肃薅螅羁芈蒁螅肃肁莇 第八章 平面问题的有限元分析及三角形单元的应用第一节 概述分析弹性力学平面问题时，最简单的单元式由三个结点组成的三角形单元。当用以分析平面应力问题时，可将其视为三角板；当用以分析平面应变问题时，则可式为三棱柱。各单元在结点处为铰结。图8-1所示位移悬臂梁离散为三角形单元的组合体以矩阵形式列出弹性力学平面问题的基本量和基本方程。谈形体所受体力分量可表示为 (

3、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）式中：，表示虚位移；，表示与虚位移相对应的虚应变

4、。为了便于计算，作用于弹性体上的体力和面力替换为作用在结点上的集中力，即等效结点荷载。设作用于各个结点上的外力分量用如下列阵来表示与这些结点外力分量相对应得结点虚位移分量列阵为则外力在虚位移上做的虚功为如平面弹性体的厚度为t，该虚功除以t，即可得出单位厚度薄板上的外力虚功。于是，式（8-11）所示虚功方程可写成 （8-11）虚功方程不仅仅应用于弹性力学，也可用于塑性力学。其应用条件是：只要变形体的全部外力和应力满足平衡方程；位移是微小的，并满足边界条件，位移与应变满足几何方程。所以，通常称为变形体虚功方程。第二节 单元分析图8-2所示为一个三角形单元。三个结点按逆时针顺序编号分别为i、j、m，

5、结点坐标分别为。图8-2由于每个结点有两个位移分量，单元共有六个结点位移分量：，如图8-2a）所示，因此三角形单元的结点位移分量e可表示为 （8-13）与这六个结点位移分量相对应得结点力也有六个分量，如图8-2b)所示 （8-14）在每个单元上，都可以把结点力用结点位移来表示，即建立如下关系式 （8-15）式中ke称为单元刚度矩阵。寻求ke的过程称为单元分析。单元分析按如下步骤ke结点力应力结点位移内部各点位移应变 一、 位移函数为了求单元内任一点（x，y）的位移，设该点的位移u、v为其坐标x、y的某种函数，单元有六个结点位移分量，在位移函数中取六个任意参数i（i=1，2，6），并将位移函数取

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的另两组公式。

7、将式（c）和(d)代入（8-16）并展开，得到以结点位移表示的位移函数 （8-20）式中，反映了单元的位移形态，故称为单元位移的形态函数或形函数。矩阵N称为形函数矩阵。选取得位移函数是否合理，要看随着单元网格的逐步细分，有限元解是否逼近于精确解。为了保证收敛型所选择的单元位移函数应满足以下条件：（1） 包含单元的刚体位移；（2） 包含单元的常量应变；（3） 保证相邻单元在公共边界处位移的连续性。二、 单元的应变和应力选择了位移函数并以结点位移表示单元内点的位移后，重新写出平面问题的几何方程 （f）由式（8-20）得 （g）将式（g）代入式（f），并利用下式 （i、j、m） (h)得单元应变 （

8、8-24）或简写成 （8-25）式中 （8-26）式（8-25）就是由结点位移求应变的转换式，其转换矩阵B称为几何矩阵。将（8-25）代入平面问题的物理方程式（8-8）有 （8-27）或写成 （8-28）S=DB (8-29)称为应力矩阵。三、 单元刚度矩阵在有限单元法中，常利用虚功方程代替平衡方程。图8-3a)所示为三角形单元的实际力系，其结点力为Fe，应力为；图8-3b)所示为单元虚位移状态，其结点位移为*e,应变为*。利用式（8-12）可得图8-3 （8-30）式中：由式（8-25）可知因此 将此公式代入式（8-30），由于中的元素是常量，公式右边的可以提到积分号的前面，得 由于虚位移是

9、任意的，则 因为B和都是常量矩阵，并且积分，所以 （8-31）利用（8-27），可得 （8-32）令 （8-33）则式（8-32）就变成式（8-15），即 单元刚度矩阵为一个66矩阵，它时单元结点位移与单元结点力之间的转换矩阵，具有以下性质：（1）示对称矩阵，其元素；（2）是奇异矩阵，由它的元素组成的行列式等于零，即它不存在逆矩阵；（3）具有分快性质。第三节 等效结点荷载 为简化各单元得受力情况，便于分析计算，应将单元所受各种载荷向结点移置，化为结点荷载，荷载的移置应安静力等效原则进行。静力等效是令原来的荷载与移置后的荷载在任意虚位移上的虚功相等。一、 集中荷载 图8-5a）所式为单元内的任意

10、一点M受到集中荷载P的作用，沿x、y方向的分量分别为Px、Py，用矩阵表示为。设移置到该单元结点上的等效荷载列阵为图8-5所示单元发生虚位移，其中单元内部任意点的虚位移为 设单元各结点虚位移为由式（8-22），则 （a）图8-5如果设Re为单元得等效结点荷载，则Re在结点虚位移上所做的虚功应与原来集中荷载在其作用点的虚位移上做的虚功相等，即将式(a)代如上式，得由于虚位移可以式任意的，因此 （8-36）或写成 （i、j、m）二、 分布体力 设单元受分布体力的作用将微分体积tdxdy上的体力当作集中荷载P，利用式（8-36）的积分得出 （8-38）如果单元上作用的分布体力为自重，为材料的重度，即

11、那么利用式（8-38）得到等效结点荷载列阵为 （8-39）三、分布面力设单元在边界ij上首有分布面力的作用 将微分面积tds上的面力当作集中荷载P，利用式（3-36）的积分可以得出 （8-40）第四节 整体刚度矩阵有限元求解弹性力学问题也要用结点平衡方程求解作为基本为质量的整体结点位移列阵。在求出各单元的单元刚度矩阵ke和解点荷载列阵后，就可以用几何的方法建立其平面弹性体的整体刚度矩阵和整体平衡方程。如果弹性体划分为m各单元，n各结点，则有 （8-44）式中：整体刚度矩阵K为2n2n阶方程；、R分别为整体结点位移列阵和整体结点荷载列阵，都是2n1阶列阵。第六节 单元网格的划分和计算成果的整理一

12、、 单元网格的划分 划分单元网格时，要综合考虑单元的数目和划分的合理性，注意以下方面：1 合理安排单元网格的疏密分布结构的不同部位疏密不同，边界比较曲折的部位，网格可以密一些，比较平滑的部位，可以疏一些，对于应力应变相对小的部位疏一些。在保证计算精度的前提下，减少单元划分数目。2 对称性的利用 当对称结构上有对称荷载或反对成荷载作用是，可以利用对称性取半结构或1/4结构。3 单元形状的合理性三角形单元的三个边长较为接近最好，可以避免由于周围应力长分布不均而引起较大的计算误差。4 不同材料截面处单元划分当计算对象为不同材料组成时，应以材料性质发生变化的不同材料界面作为单元的边界。 虿薈肂肈芅蚁袅羄芅袃膀莃芄薃羃艿芃蚅腿膅节螇羁肁芁袀螄荿莀蕿羀芅荿蚂螂膁荿螄羈肇莈薄螁肃莇蚆肆莂莆螈衿芈莅袀肅膄莄薀袇肀蒄蚂肃羆蒃螅袆芄蒂蒄肁芀蒁蚇袄膆蒀蝿腿肂葿袁羂莁蒈薁螅芇蒈蚃羁膃薇螆螃聿薆蒅罿羅薅薈螂莃薄螀羇艿薃袂袀膅薂薂肅肁薂蚄袈莀薁螆肄芆蚀衿袆膂虿薈肂肈芅蚁袅羄芅袃膀莃芄薃羃艿芃蚅腿膅节螇羁肁芁袀螄荿莀蕿羀芅荿蚂螂膁荿螄羈肇莈薄螁肃莇蚆肆莂莆螈衿芈莅袀肅膄莄薀袇肀蒄蚂肃羆蒃螅袆芄蒂蒄肁芀蒁蚇袄膆蒀蝿腿肂葿袁羂莁蒈薁螅芇蒈蚃羁膃