有限元分析第三章.docVIP

最小势能原理和分片插值 　　最小势能原理和分片插值是有限单元方法的核心内容之一，在这一章中我们将介绍最小势能原理的具体应用，同一力学问题的几种不同的表达方式及它们之间的联系；Ritz方法在单元内的应用；两种最常用的插值形式（Lagrange型和 Hermite型）。协调的位移型单元的收敛条件。 §3-1最小势能原理 　　一个平衡问题，可以至少用以下叁种不同的方式加以描述： 　　　　（i）平衡方程 　　　　（ii）虚位移原理 　　　　（iii）总势能取驻值（函数的极值问题）　 有限自由度系统　质点系 图3-1(a)为两个重分别为PA, PB的小球，由 不计重量，弹性系数为k的弹簧相连，放置在光滑的曲面F(x,y)=0上。该系统的平衡问题可由以下三种方法来描述： （１）平衡方程　（在平衡位置描述平衡） 分别取Ａ、Ｂ为研究对象，建立一直角坐标系，如图3-1 (b)所示。由Ａ、Ｂ两球组成的系统的平衡条件为： 其中，弹簧力以压力为正，即： （２）虚位移原理 (在平衡位置附近描述平衡) 系统的平衡位置可描述为 rA，rB　或　 xＡ、yＡ，xＢ、yＢ　如图3-1 (c) 所示。　　　　　　　　　 两质点的虚位移分别为： δrＡ, δrＢ 或　δxＡ,δyＡ,δxＢ,δyＢ　　　　　　 由虚位移原理可得 上式也称为微分型变分原理。进一步改写为 (3) 总势能取驻值 （在全范围内描述平衡（寻找平衡位置））引入总势能函数 由总势能函数取驻值，即： 寻找平衡位置。以上三种描述系统平衡的方法，对于有限自由度系统而言是等价的。然而，三种方法的特殊性是：平衡方程是在平衡位置上描述平衡问题；虚位移原理是在平衡位置附近描述平衡问题；而最小势能原理是在全局范围上描述（寻找）平衡问题。 最小势能原理可以表述为：在所有满足给定位移边界条件和协调条件的位移中，满足平衡条件的位移使总势能取驻值，若驻值是最小值，则平衡是稳定的。 作为一条基本原理，它的正确性应由事实加以检验，不必从理论上给以证明。但为了便于理解“位移边界条件“和“协调条件“的含义，下面通过三个例子对最小势能原理和平衡方程的等价关系加以验证。 ２、无限自由度系统　弹性体　 轴向受拉的直杆 设杆长为Ｌ，截面积为Ａ，轴向分布载荷f(x)。x=0端固定，x=L端受端点集中力P。 设位移u(x)满足： （i）u(0)=0 (位移边界条件) （ii）u(x) 在[O, L]上连续（协调条件） (iii) 使 取最小值。 若u(x)+δu(x)为不同于u(x) 的另外一种位移分布函数，也满足上述的位移边界条件和协调条件，由 可推出 而 不难求得： 因πP(u) 取最小值，即 的充分必要条件是：对任意满足(3-1-2)的δu(x) 有 式(3-1-3)即势能驻值条件。 　　若补充假定u’’(x) 存在、连续，则对（3-1-3）分部积分一次，并利用（3-1-2）, 可得到 (3-1-4) 式对任意δu(x) 都成立的充分必要条件是： 由上述过程不难看出： 　　①　由势能取驻值可以推出平衡方程。反之也对，说明两种描述方法在力学上等价。 　　②　两种描述方法对边界条件的要求不同。用微分方程描述时，u必须满足： 用最小势能原理描述时，要求函数满足位移边界条件 　u(x) 使πP精确取驻值时，平衡方程和力边界条件将精确地得到满足；当u(x) 使πP近似取驻值时，平衡方程和力边界条件只能近似的得到满足。在这种意义下，力边界条件又称为自然边界条件(Natural Boundary Condition )，而位移边界条件由称为强制边界条件（Essential Boundary Condition）。 　　③两种描述方法对函数的光滑程度（即可微性）要求不同。用微分方程描述时要求u(x)　有连续的二阶导数（记作u∈C2(0, L)）。而用最小势能原理描述时，为了保证变形能存在，要求u’(x)平方可积（记作u∈H１(0, L)） 平面应力问题 正方形区域边长为a，厚度为t，受到体积力 (fx, fy), 边界AB固定。边界BC、CD自由。边界AD的法向应力为q(x)，切向应力为p(x)。势能πP和强制边界条件为： u∣AB=v∣AB=0 利用 以及格林公式 （其中LX，LY为区域Ω之边界Г的外法线n的方向余弦）。势能πP的驻值条件为： 注意到沿边界Г，外法线n的方向余弦为 AB BC CD DA LX -1 0 1 0 LY 0 -1 0 1 以及沿线AB：δu=δv=0 则(3-1-5)对任意δu,δv都成立的充分必要条件为 沿BC：σy=τxy=0 沿CD：σx=τxy=0 沿AD：σy=q,τxy=p 关于边界条件的讨论以及u, v可微性的讨论与上例相似，这里从略。 梁的平面弯曲 对于图3-4