演示文稿有限元第三章最小势能原理和分片插值.pptVIP

演示文稿有限元第三章最小势能原理和分片插值目前一页\总数二十七页\编于二点（优选）有限元第三章最小势能原理和分片插值目前二页\总数二十七页\编于二点§3-1最小势能原理平衡问题，可以至少用以下叁种不同的方式加以描述：（i）平衡方程（ii）虚位移原理（iii）总势能取驻值（函数的极值问题）1.有限自由度系统质点系图3-1(a)为两个重分别为PA,PB的小球，由不计重量，弹性系数为k的弹簧相连，放置在光滑的曲面F(x,y)=0上。该系统的平衡问题可由以下三种方法来描述：F(x,y)=0ABAByxPBPAT’TrBrAo（c）图３－１（a）ABTNANBPBPAT’xyo（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)）(NaturalBoundaryCondition)（EssentialBoundaryCondition）②两种描述方法对边界条件的要求不同。用微分方程描述时，u必须满足：目前六页\总数二十七页\编于二点(2)平面应力问题正方形区域边长为a，厚度为t，受到体积力(fx,fy),边界AB固定。边界BC、CD自由。边界AD的法向力为q(x)，切向力为p(x)。nnnx,uP(x)(fx,fy)BDCAOy,v图3-3nq(x)u∣AB=v∣AB=0（其中LX，LY为区域Ω之边界Г的外法线n的方向余弦）。(能量泛函)格林公式目前七页\总数二十七页\编于二点总势能πP的驻值条件为：(3-1-5)注意到沿边界Г，外法线n的方向余弦为ABBCCDDALX-1010LY0-101以及沿边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)xQ图３－4Mv(3-1-6)总势能和强制边界条件为势能驻值条件对上式分部积分两次，并注意到由于必须满足强制边界条件δv(0)=δv/(0)=0则有(3-1-7)使(3-1-7)对任意δv(x)都成立的充分必要条件是：（平衡方程）（自然边界条件）目前九页\总数二十七页\编于二点微分方程的阶数为４。关于v’’、v’’’的边界条件为自然边界条件