弹性力学有限元第九章FEM部分之三一维单元及二维单元.pptVIP

河海大学机电工程学院力学教研室第九章一维单元和二维单元

第九章一维单元和二维单元谨向徐芝纶院士表示深深的敬意

第九章一维单元和二维单元§9-1线性单元形函数的概念散热片一维热传导问题温度沿单元的线性近似

第九章一维单元和二维单元§9-1线性单元形函数的概念典型单元的线性温度分布：按照单元节点的温度：计算出c1c2并回代入温度方程前一章中的典型柱体单元位移用形函数表示

第九章一维单元和二维单元§9-1线性单元形函数的性质形函数固有的性质之一是在相应节点上的值为1，相邻节点上为0。

第九章一维单元和二维单元§9-1线性单元形函数的性质形函数固有的性质之二是形函数的和为1。线性形函数对于X的导数之和为0。

第九章一维单元和二维单元§9-2二次单元通过增加分析中使用的线性单元的数量，或者使用高阶的插值函数，能够增加有限元计算的精度。二次函数要求至少使用3个节点来定义一个单元。

第九章一维单元和二维单元§9-2二次单元一般而言，二次单元由节点值表示的任何可变参数可以写为二次形函数的性质相应节点值为1，相邻节点值为0；形函数之和为1；二次形函数关于X的导数之和一般不为0。

第九章一维单元和二维单元§9-3三次及以上单元三次函数需要用4个节点来定义一个单元。当插值函数的阶数增加时，用Lagrange插值函数代替以上方法得到形函数。设节点i的整体坐标为Xi则对应的形函数Si满足从而求出参数a1

第九章一维单元和二维单元`§9-3三次及以上单元三次以上的单元，可以从Lagrange插值多项式得到N-1阶多项式的形函数：形函数的性质就是Lagrange插值基函数的性质

第九章一维单元和二维单元§9-4整体坐标，局部坐标和自然坐标整体坐标：表示节点的位置和单元的方向；施加边界条件和荷载；表示计算的结果；比如节点的位移。局部坐标和自然坐标：构造几何关系，计算积分。

第九章一维单元和二维单元§9-4整体坐标，局部坐标和自然坐标自然坐标是局部坐标的无量纲形式。在计算单元刚度矩阵或者热传导矩阵时，需要用到从上限1到下限－1之间的积分(典型的Gauss-Legendre积分)。一维线性自然坐标

第九章一维单元和二维单元§9-4整体坐标，局部坐标和自然坐标自然坐标和局部坐标的关系一维线性自然坐标x是局部坐标局部坐标下的形函数：自然坐标下的形函数：一维散热片单元的温度分布：

第九章一维单元和二维单元§9-5等参单元等参(isoparametric)的概念用自然形函数表示位移：用同样的自然形函数表示几何关系(局部坐标和全局坐标)：所谓等参，就是对单元几何形状和单元内的参变量函数采用相同数目的节点参数和相同的形函数进行变换。用这种方式表示的单元称为等参单元。

第九章一维单元和二维单元§9-5等参单元问题求X=8cm处的温度X=8cm处位于单元③注意到，X=8cm在单元③的局部坐标单元③l=5

第九章一维单元和二维单元§9-5等参单元用一维自然坐标表示二次和三次形函数自然坐标下二次形函数：自然坐标下三次形函数：

第九章一维单元和二维单元§9-6Gauss-Legendre积分自然坐标是局部坐标的无量纲形式有限元需要对单元进行数值积分数值积分的方法主要有：梯形法，Simpson法?计算离散数据的积分Gauss-Legendre积分?计算离散数据的积分?计算不等间距点上已知函数的积分Gauss-Legender积分的主要目标是用特定权系数的乘积之和，以及一些选定点上的函数值来计算积分。

第九章一维单元和二维单元§9-6Gauss-Legendre积分数值积分的形式Guass-Legendre积分的形式注意到Guass-Legendre积分的上下限为－1到1(自然坐标)Guass-Legendre积分的权因子和样本点点数权因子wi样本点li21.000000001.00000000-0.5773502690.57735026930.555555560.888888890.555555560.77459666900.774596669

第九章一维单元和二维单元§9-7矩形单元一维问题用线段近似描述，而二维问题用平面近似描述。

第九章一维单元和二维单元§9-7矩形单元温度沿单元的分布是关于坐标X和Y的函数。任意单元在局部坐标下的温度分布上式描述的单元具有■在单元边缘上线性变化；■在单元内部非线性变化；双线性单元：BilinearElement考虑节点的温度，必须满足条件：应用这些条件，可求出各参数，并得到各形函数。

第九章一维单元和二维单