最新第6章形函数.docxVIP

I j I j TOC \o 1-5 \h \z 第六章单元形函数的构造 1 6.1形函数构造的一般原理 1 \o Current Document 6.2形函数的性质 7 6.3用面积坐标表达的形函数 8 \o Current Document 6.4有限元的收敛准则 10 \o Current Document 6.5等效结点载荷列阵 11 \o Current Document 6.5.1 单元载荷的移置 11 \o Current Document 6.5.2结构整体载荷列阵的形成 11 \o Current Document 6.5.3载荷移置与静力等效关系 12 \o Current Document 习题 14 第六章 单元形函数的构造 在有限单元法的基本理论中， 形函数是一个十分重要的概念， 它不仅可以用作单元的内插函数, 把单元内任一点的位移用结点位移表示，而且可作为加权余量法中的加权函数，可以处理外载荷， 将分布力等效为结点上的集中力和力矩，此外，它可用于后续的等参数单元的坐标变换等。 根据形函数的思想，首先将单元的位移场函数表示为多项式的形式，然后利用结点条件将多项 式中的待定参数表示成场函数的结点值和单元几何参数的函数，从而将场函数表示成结点值插值形 式的表达式。在本节中，重点讨论几种典型单元的形函数插值函数的构造方式，它们具有一定的规 律。然后以平面三角形单元为例，讨论了形函数的性质，在此基础上分析了有限元的收敛准则。 6.1形函数构造的一般原理 单元的类型和形状决定于结构总体求解域的几何特点、问题类型和求解精度。根据单元形状 ，可 分为一维、二维、三维单元。单元插值形函数主要取决于单元的形状、结点类型和单元的结点数目。 结点的类型可以是只包含场函数的结点值，也可能还包含场函数导数的结点值。是否需要场函数导 数的结点值作为结点变量一般取决于单元边界上的连续性要求， 如果边界上只要求函数值保持连续， 称为CO型单元，若要求函数值及其一阶导数值都保持连续，则是 C1型单元。 在有限元中，单元插值形函数均采用不同阶次的幕函数多项式形式。 对于CO型单元，单元内的 未知场函数的线性变化仅用角（端）结点的参数来表示。结点参数只包含场函数的结点值。而对于 C1型单元，结点参数中包含场函数及其一阶导数的结点值。与此相对应，形函数可分为 Lagra nge 型（不需要函数在结点上的斜率或曲率）和 Hermite型（需要形函数在结点上的斜率或曲率）两大 类，而形函数的幕次则是指所采用的多项式的幕次，可能具有一次、二次、三次、或更高次等。 另外，有限元形函数［N］是坐标x、y、z的函数，而结点位移不是 x、y、z的函数，因此静力学 中的位移对坐标微分时，只对形函数 ［N］作用，而在动力学中位移对时间 t微分时，只对结点位移向 量作用。 （1） 一维一次两结点单元 Xi Xj 精选文库 精选文库 PAGE PAGE # 图3-8 一维一次两结点单元模型 设位移函数u(x)沿x轴呈线性变化,即u(x)% a2x写成向量形式为a1 u(x) 1 xa2设两个结点的坐标为 xi,xj 设位移函数u(x)沿x轴呈线性变化,即u(x) % a2x 写成向量形式为 a1 u(x) 1 x a2 设两个结点的坐标为 xi,xj ;两结点的位移分别为 1 1 Xi Ui 1 Xj Uj Ui，Uj，可以代入上式并解出 da，得 ai a2 位移函数U(X )记作形函数与结点参数乘积的形式 1 Xi Ui 得到形函数为 [N] 1 x u(x) Xj Uj Xj 1 1 Xj Xj xxx Nj Nj xj x Xj Xi X xi Xj Xi 在自然坐标系内进行定义，则可得到形函数的标准化形式 1_ 1_ 2 2 2丄 1 , L2 1 。 L [N] Ni Nj 其中，自然坐标的变换公式为 L (5.90) (5.91) (5.92) (5.93) (5.94) (5.95) X Xi L1 1 ■ L2 1 i 1 Xj 1 图3-9 一维一次两结点单元的局部坐标表达 二维一次三结点单元(平面三角形单元)在总体坐标系统下，任一点的某一方向的位移是 u(x,y) a1 a2X aay (5.96) 设三个结点的坐标是 xi, % : 「 Xj, yj : ，xk, yk ，ui ,Uj, Uk为三个结点在某方向上的位移， 具有如下关 a1 a 1 X 1 yi Ui 系 u 1 x y a2 a2 1Xj yj Uj (5.97) a3 a3 1 Xk yk Uk 得到形函数矩阵如下式 1Xi 1 yi [N] 1 x y 1Xj yj (5.98) 1 Xk yk  上