弹性直梁问题的变分原理及有限元素法.docx

第二章 弹性直梁问题的变分原理及有限元素法 讨论的问题：一变剖面的梁，一端 x=0固支，另一端 x=：l简支。承受轴向拉 力N，分布横向载荷q x以及端点弯矩M|的作用。 控制微分方程及边界条件（以梁的挠度 w表示） d2dP.2Ejd_wEJ 2dx.2 d2dP .2 Ejd_w EJ 2dx .2 d w _N 牙二q dx d2dx2 .2 d w . EJ 2 Nw dx2 在x = 0处: 在x =1处: （广义固支） -EJ d 2w dx2 自然边界条件（M M|^0 称谓：把满足方程及全部边界条件的挠度叫真实挠度，精确解；把满足基本边界条 件但不满足微分方程和自然边界条件的挠度叫（变形）可能挠度。 i） 最小势能原理（变分原理） 把载荷看作是不变的已知函数，把挠度看作是可变的自变函数。 整个系统的势能包括三部分: （1）梁的应变能： （2）轴向应变能:ioEJ1 Ydw （2）轴向应变能: i oEJ 1 Ydwf 2oN 1 V列2 2 dx （3）横向载荷势能: l —二 l — 二 p = - 0qwdx M lw l 后项取加号，是为着能够得到自然边界条件的结果 ⑷系统总势能二： N -dw I — qwSx + MMI) idx丿 *除w为可变外，其余变量假定为已知的不变量。 最小势能原理：在所有变形可能的挠度中， 精确解使系统的总势能取最小值 。 由于二w是w的二次函数，不用变分法而用较初等的方法也能作出数学证明。 证明过程: 设wx是精确解，它满足微分方程及所有边界条件。 设wk x是某一变形可能的挠度，仅知道它满足: dw ... w = w0, 0 dx 在X =1处: 在X =1处: 令：.w = w 令：.w = wk -w= wk = w lw --w是一有限量，而不是变 分 当w- 0时，才是w的变分。 由上式关系知 w满足: 在 x =0处，lw 在 x =0处，lw = 0j_w = 0 与Wk相应的总势能: .2 \ d wk dx2 ； 2 -. -qwk dx M lwk l J =1「wk 二二 w =w 二二 w ］亠 2〔」11 w,二w 厂丨 = w 其中: 01 01 dx dx2 詈 7 wdx Ml w l 二2 w「鳥 EJd LW-N2■d^w 二2 w「鳥 EJ d LW -N 2 ■d^w 丫 dx dx Note: ：w,w满足基本边界条件，可得 （功的互等定理） 冲wdx-M.w I J EJd2 Wd2w -dWdW dx dx dx dx dx 功的互等定理：第一组力在第二组位移上所做的功等于第二组力在第一组位移上所做 的功。（真实力状态关于虚拟位移作功；虚拟位移产生的虚拟力关于真 实位移作功）此式表明： ——11 n （w,心w ）= 0 二；Wk =二 w j亠丨.上 W 2 2 由二：w的算式，当N 一0时，二 ：w -0 或 当N :: 0 （轴受压）未达到临界压力时， 二2丄w _0。 所以： 二 wk「：i.[w 式中的等号只有在 w为刚体位移时才能成立（即弹性能仅为零的情况） 。 以上即是证明的最小势能原理。 上述的证明可普遍适用于其他类型的边界条件；也适用于其他复杂弹性力学体系。 精确解既然是使总势能取最小值，那么必有： 二 w = 0 即： ]~d2 f d2 、 d2 _ 6n = [ |—- EJ —耳-N —弓 一4 弘dx + Ejw(l )+ M i 5w(l )= 0 0 ^dx\ dx2 丿 dx2 .2 、 EJdF-N EJw”I M| =0 （正好补足了可能挠度尚未满足的边界平衡条件） 所以，最小势能原理与平衡条件（边界及内部）完全等价。 说明：尽管变分法与原问题的微分方程系统等价，但具体求解时，变分法涉及的导 数阶次要低（这是能量法的优点，求得的解可能满足微分方程的连续性要求， 也可能不满足） 能量逼近解（不是真解）不能满足力的自然边界条件（选取可能位移时一般很难取得满足力的边界条件，故一般情况有限元不能获得满足力边界条件的解） 。从这点上看，只 能是近似解。当然，如果在满足力的边界条件的那些函数集合中选， 则解的精度要大大提高 （即更为逼近）。 强调三点： 上述证明的是真实挠度使 系统势能取最小值。 又通过变分法证明了系统能量的极值曲线满足梁的微分控制方程（平衡方程） 。需 要深入认识的是：系统能量泛函中要求的挠度为可能挠度 （即满足连续性与位移边界条件的 曲线，但不满足微分方程）；变分的结果恰使极值曲线应满足微分方程及自然边界条件，补 足了真实解的所需条件（理论上的等价性） 。 微分方程解与能量泛函解对函数连续性的要求是不同的， 故能量方法对函数的连续 性要求较放松。以后在广