有限元第三章最小势能原理和分片插值论述.pptVIP

在这一章中我们将介绍： 最小势能原理及具体应用； 同一力学问题的几种不同的表达方式及它们之间的联系； Ritz方法在单元内的应用 两种最常用的插值形式（Lagrange型和 Hermite型）； 协调的位移型单元的收敛条件。 第三章 最小势能原理和分片插值 (有限单元方法的核心内容之一) A B 0 k j n m x, u y, v l i s 0 y,v 2b 2a x,u j i l k 0’ y’ x’ §3-1最小势能原理 平衡问题，可以至少用以下叁种不同的方式加以描述： 　（i）平衡方程 　（ii）虚位移原理 　（iii）总势能取驻值（函数的极值问题） 1. 有限自由度系统　质点系 图3-1(a)为两个重分别为PA, PB 的小球，由 不计重量，弹性系数为 k 的弹簧相连，放置在光滑的曲面F(x,y)=0上。该系统的平衡问题可由以下三种方法来描述： F(x,y)=0 A B A B y x PB PA T’ T rB rA o 　（c） 　图３－１ （a） A B T NA NB PB PA T’ x y o （b） (1) 平衡方程　 (2) 虚位移原理 (3) 总势能取驻值 在所有满足给定位移边界条件和协调条件的位移中，满足平衡条件的位移使总势能取驻值，若驻值是最小值，则平衡是稳定的。 最小势能原理和平衡方程是否等价？ ２. 无限自由度系统　弹性体 Ｏ Ｌ f(x) x, u Ｐ 图　3-2 (1) 轴向受拉的直杆。 设杆长为Ｌ，截面积为Ａ，弹性模量为 E 轴向分布载荷f(x)。x=0 端固定，x=L 端受端点集中力P。 设位移u(x)满足： (i) u(0)=0 (位移边界条件) (ii) u(x) 在[O, L]上连续（协调条件） (iii) 使总势能取最小值。 (3-1-1) u(x)即为该问题的解 最小势能原理： 总势能=变形能—外力之功 设：u(x)+δu(x)为不同于u(x) 的另外一种位移分布函数，也满足上述的位移边界条件和协调条件 ，则 (3-1-2) 将u(x)+δu(x)代入总势能函数 考察两总势能函数之差 因πP(u) 取最小值，即 的充分必要条件是： 对任意满足(3-1-2)的δu(x) 有： (3-1-3) 若假定 u’’(x) 存在、连续，则对（3-1-3）分部积分一次，并利用（3-1-2）, 可得到 (3-1-2) (3-1-3) (3-1-4) (3-1-4) 式对任意δu(x) 都成立的充分必要条件是： （平衡方程） （力边界条件） 　　①　由势能取驻值可以推出平衡方程。反之也对，说明两种描述方法在力学上等价。 （平衡方程） （力边界条件） （位移边界条件） 用最小势能原理描述时，要求函数满足位移边界条件 　 而力边界条件将作为势能取驻值的自然结果。 ③两种描述方法对函数的光滑程度（即可微性）要求不同。用微分方程描述时要求u(x)　有连续的二阶导数（记作u∈C2(0, L)）。而用最小势能原理描述时，为了保证变形能存在，要求u’(x)平方可积（记作u∈H１(0, L)） (Natural Boundary Condition ) （Essential Boundary Condition） 　　②　两种描述方法对边界条件的要求不同。用微分方程描述时，u必须满足： (2) 平面应力问题 正方形区域边长为a，厚度为t，受到体积力 (fx, fy), 边界AB固定。边界BC、CD自由。边界AD的法向力为q(x)，切向力为p(x)。 n n n x, u P(x) (fx, fy) B D C A O y, v 图　3-3 n q(x) u∣AB=v∣AB=0 （其中LX，LY 为区域Ω之边界Г 的外法线 n 的方向余弦）。 (能量泛函) 格林公式 总势能πP 的驻值条件为： (3-1-5) 注意到沿边界Г，外法线 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 略去了积分过程 (2) 梁的平面弯曲 Ｏ Ｌ q(x) x Q 图 ３－4 M v (3-1-6) 总势能和强制边界条件为 势能驻值条件 对上式分部积分两次，并注意到由于必须满足强制边界条件 δv(0)=δv/(0)=0 则有 (3-1-7) 使(3-1-7)对任意δv(x) 都成立的充分必要条