第4章 插值与基函数(下).pptVIP

在等参数单元中，最常用的平面单元是四边形四节 点、四边形八节点，前面已谈过即使是任意四边形 ，也不能过份“任意”，在8节点元中还要注意，由于 标准矩形中的 点都是在边中点的，故单元棱边上 的点 也尽量在中点。 否则对应变精度有影响。右图情况中点移到 处， 角点1的应变 （ 当然 （五）三维情形 在实际工程中并非所有问题都可简化为平面问题的， 如拱坝就是一个典型的空间问题，因此要引进三维 有限元。 最简单的三维插值是四面体线性插值，基函数为 ，在第二章“有限元方法基本方程 断裂力学中用此来解裂缝问题，另当别论），曲边单元也要注意。 j m i l 的建立”中已谈过。 实际上，这4个基函数即为体积坐标 ，与 面积坐标有相似的性质和定义，且有算式 （4.60） 为单元体积。 但这种单元精度差，工作量大（占内存多），部分 易出错，用途不大，略过,重点介绍下面两种单元。 1．8节点等参单元 这是个六面体等参元 在三维空间中取一个标准区域 A1 A2 A3 A4 A5 A7 A6 A8 O ξ ζ η 比平面问题多一个ζ坐标，问题就由平面扩展为三线 性插值，基函数表为 （4.61） 空间任意给定八个对应节点的整体坐标值 ，由基函数可决定一个8节点等参单元，相应坐标变换为： （4.62） Jocobi矩阵 （4.63） （4.64） 此式将整体坐标的积分转到局部坐标的积分。 由于 得到 （4.65） 单元刚度阵 8结点q=8 其中 单刚中的元素 （4.66） 由于 很难写出解析式，所以都要用数值积分。 8点元六面一般说不是平面，而是两族直线交织成的 直信面。但剖分时不需知道这些面的准确形状，剖分 时只要有大致形状即可。 2．20结点等参单元 内部插值有更高精确度， 边界上对曲面边界的逼近 更好，在每边中点给一个 插值条件。 不完全的三二次多项式插 值，所谓三点是指三个方 1 2 3 4 5 7 6 8 O ξ ζ η 15 11 18 12 19 20 10 13 14 9 16 17 向均为不完全二次插值。完全的三二次插值有 点，故现在要加7个限制条件。 去掉 公项， 插值多项式仍唯一。 在ξ=1,η=1, ζ=1三平面上，除1，9，13，17四点外， 其余16点都已包括，这样再过9，13，17三点作一平 面，其方程 故 类似方法得全部基函数（注意：按上页图示编号情况下） （4.67） 需要说明的问题： （1）每条棱边均为二次曲线，侧面为两族二次曲线交 织成的曲面，不必知这些曲面、曲线的准确形状，只要 知20点坐标即可。（即 )； （2）坐标变换式，Jacobi矩阵，基函数对整体坐标的 微商与对局部坐标微商的关系均与8点等参单元相似， 只要把“8”改成“20”即可。 （3）上述两种元形状也不可过份“任意”，也要考虑 （4）8点三线性和20点三二次插值所缺七项都对局部 坐标而言，对整体坐标并非如此。 （六）数值积分和单元特性的计算 求“刚度矩阵”和“荷载向量”时，并计算 ， 之类积分，用高次插值尤其等参元时，被积函数复杂，只能数值积分。因此数值积分占很重要地位。 问题。 1．一维高斯积分 不失一般性，可设积分区间为[-1,1]，数值积分公式一 般表为 （4.68） 称积分点， 称权系数。 衡量积分公式好坏的标准称代数精确度。 代数精确度m——即该公式对m阶多项式精确成立，如几个积分点给定，则由几个点可唯一确定一个n-1次多项式，它就是Lagrange插值多项式。 o x f(x) -1 +1 f(x1) f(x2) f(x3) f(x4) x2 x3 x1 x4 若证 则Lagange插值多项式为 若n取2，则为Lagrange线性插值 （与原来一样的） （4.69） 如Lagrange插值多项式的积分代替f