lec-13-等参数单元.pptVIP

位移插值函数取下列多项式形式 当 固定时为 的二次函数，当 固定时为 的 二次函数。因此这样的位移插值函数是称为双二次 函数，相应的插值称为双二次插值。 2、八结点等参数单元的形状函数 八结点四边形等参数单元在局部坐标系下的位移插值函数为 其中形函数 由下列条件唯一决定： （1）为双二次函数； （2） 在结点 的值为1，而在其余结点 的值为0，即 以 为例说明。 在结点1其值为1，在结点2-8其值为0。直线35、57与28通过2-8这7个结点，他们的方程分别为 因此，函数 在结点2-8的值为0，而且满足（1 ，是双二次函数，再利用 在结点1（-1，-1）的值为1的条件，可以写出 如下： 再来求 。它在结点2其值为1，在结点1和3~8其值为0。直线17 、35、与57通过这7个结点，他们的方程分别为 因此，函数 在这些结点的值为0，而且满足（1），是双二次函数，再利用在结点2（0，-1）的值为1的条件，可以写出 如下： 形函数的统一表达式： （5-41） 作为等参数变换，坐标变换式取插值函数（5-40） 同样的形式 （5-42） 3、等参数变换及其实现条件 给出整体坐标下8个点的坐标值 ， 由等参数变换的思想可写出局部坐标 到整体坐标的变换式 通过坐标变换式可以了解整体坐标下单元的形状。以局部坐标下单元345边为例： 通过这三个节点的直线方程为 ，将其代入（5-42），可得345边在整体坐标下的参数方程形式 消取 后，就得到二次抛物线方程，实际单元为四条抛物线围成的曲边四边形，如图5-9所示，它可以较好地模拟计算区域的曲线边界，使计算精度提高。 4、等参数单元的形态 为保证等参数坐标变换能顺利进行，仍需使 整体坐标下8个节点的位置配置和单元形状需要作一定限制。 （1）不能使单元太歪斜 （2）四边形的任意两条对边即使通过适当的延长也不能再单元上出现交点。 （3）每边的中间点应尽量在两角点的当中，若位于1/3分点，则会出现较大计算偏差；若位于1/4分点，则计算结果会完全不正确。 划分单元时应尽量做到：1）单元划分尽量接近正方形；2）中间节点尽量位于每边的1/2分点处。 图 5、等参数单元的特点 优点： 1）等参数单元的计算精度高，可以较好地拟合曲线边界的求解区域； 2）它所需要的数据量少。（例：对一个平面应力分析的问题，三角形单元用300个左右的单元，精度于应力集中处不理想；而8结点等参数单元只用50个，且精度高。） 缺点： 1）划分单元时有特殊要求。 2）要进行等参数变换，计算复杂，使程序编写复杂。 5-4空间问题等参数单元 1、八结点六面体等参数单元 先在局部坐标 下考察八结点的正六面体单元，边长为2，中心在原点，取八个顶点，如图5-10所示 取位移模式如下形式： 5-43 当一个自变量固定时，（5-43）为另两个 自变量的双线性函数。因此在立方体的每一侧面，位移模式完全由侧面上的四个结点的位移分量所唯一决定。 将位移模式写成如下形式 其中形函数 由下列两个条件唯一决定： （1） 为形如（5-43）的双二次函数； （2） 在结点 的值为1，而在其余结点 的值为0，即 各结点的局部坐标为： 下面具体分析形状函数 下面具体分析形状函数 以 为例来说明。它在结点1取值为1，在结点2~8其值为0。平面2376、3487与5678通过2~8这7个结点，他们的方程分别为 容易求得 （5-46） 类似可得到 至 的表达式，统一写成 按等参变换的思想，由局部坐标到整体坐标变换的表达式： （5-48） 1）经过变换局部坐标下的直棱边， 对应于整体坐标下也 具有有直棱边。 2）在局部坐标下的每一侧面经过变换 后在整体坐标下不一定表示为平面。 思考题 1、有限元分析中，采用等参数单元的主要优点是什么？（3点） 2、在等参变换时，规则（母）单元位移模式必须满足的条件是什么？（ 3点：完备性，连续性） 3、为保证等参数单元的精度，结点的位置配置和单元形状需要作哪些限制？ 计算题型： 已知规则单元的点和线，求等参单元的点的坐标和线的方程。 8.1.如图所示一个等参单元，其中各点的坐标为：1（1，1）， 2（3，1）， 3（2.5，2），4（1.25，1.75）， （1）试通过等参变换确定局部坐标中 的点映射到等参数单元的坐标。 （2） 求在母单元中， 的直线映射到等参单元 中的直线的方程。 8.2 试证明（a）局部坐标系下的八节点的正方形单元的781边（ ）在整体坐标下（b）的边为二次抛物线形状。 （a） （b） Lecture 13 等参数单元 1. 等参数单元的引入 在平面有限元分析中，采用简单的三结 点三角形单元和四个结点的矩形单元。 *三角形