《弹塑性力学》第十章弹性力学的能量原理.ppt

《《弹塑性力学》第十章弹性力学的能量原理.ppt》由会员分享，可在线阅读，更多相关《《弹塑性力学》第十章弹性力学的能量原理.ppt（130页珍藏版）》请在三一办公上搜索。

1、2023/8/31,1,第十章 弹性力学的能量原理,10-1 几个基本概念和术语,10-2 虚功方程,10-3 功的互等定理,10-4 虚位移原理和最小势能原理,10-5 虚应力原理和最小余能原理,10-6 基于能量原理的近似解法,2023/8/31,2,第十章 弹性力学的能量原理,弹性力学的解法之一为弹性力学边值问题求解体系静力法。在前面各章中就围绕平面问题、扭转问题和空间轴对称问题进行了具体分析和研究。,2023/8/31,3,第十章 弹性力学的能量原理,弹性力学问题的解法还有另一种解法：以能量形来建立弹性力学求解方程能量法（从数学意义上说也可认为变分法）。,本章主要介绍几个基本能量原理以

2、及基于能量原理的近似解法。在介绍能量原理以前，先介绍几个基本概念和术语。,2023/8/31,4,10-1 几个基本概念和术语,1.1应变能U和应变余能Uc：,应变能 U 在第四章中已定义过：,应变能密度,2023/8/31,5,弹性关系,如果将几何关系引入应变能,U、W 为位移的函数。,应变余能（类似应变能）定义,10-1 几个基本概念和术语,2023/8/31,6,应变余能密度,单位体积的应变余能,Wc 与积分路径无关，只与终止状态和初始状态有关。,Wc=ijij 为全微分,10-1 几个基本概念和术语,2023/8/31,7,逆弹性关系,且 W+Wc=ijij,10-1 几个基本概念和术

3、语,2023/8/31,8,但,10-1 几个基本概念和术语,材料为线弹性时,2023/8/31,9,各向同性线性材料的应力应变关系,10-1 几个基本概念和术语,将几何关系引入上式,U=U(ui)应变能是位移的函数,2023/8/31,10,代入Uc表达式,各向同性线性材料的应力应变关系,10-1 几个基本概念和术语,2023/8/31,11,10-1 几个基本概念和术语,应变能、应变余能的计算举例,图示等截面杆，承受轴向荷载 P 作用。杆截面面积为 A,材料应力应变关系分别为（1）=E，(2)=E 1/2.试计算外力功T、应变能U 和应变余能Uc。,解：,（1）=E,2023/8/31,1

4、2,10-1 几个基本概念和术语,T=U=Uc=P l/2,l=Pl/(EA),P=N=lEA/l,U=l 2EA/(2l),Uc=P 2 l/(2EA),（2）=E 1/2,T=U=,2023/8/31,13,10-1 几个基本概念和术语,T=U,2023/8/31,14,10-1 几个基本概念和术语,Uc=P l U=P l-U,2023/8/31,15,10-1 几个基本概念和术语,作业：图示结构各杆等截面杆，截面面积为A,结点C承受荷载P作用,材料应力应变关系分别为（1）=E，(2)=E 1/2。试计算结构的应变能U 和应变余能Uc。,2023/8/31,16,1.2可能位移 ui(k

5、)和可能应变 ij(k)：,可能应变ij(k)：由ui(k)通过几何方程导出的,可能位移ui(k)：在V内连续且可微，在 su上 满足：,10-1 几个基本概念和术语,2023/8/31,17,1.3可能应力 ij(k)：,ij,j(k)+fi=0(a),在s上满足(b),满足式(a)、（b)满足静力方程,可能应力 ij(k)：在V内满足,10-1 几个基本概念和术语,2023/8/31,18,两种可能位移ui(k1)和ui(k2)之差称为虚位移 ui，而由两种可能位移状态对应的可能应变 ij(k1)、ij(k2)之差称为虚应变ij。,1.4虚位移 ui和虚应变 ij：,ui=0 在su上齐次

6、位移边界条件。,ij=(ui,j+uj,i)/2 在V内,10-1 几个基本概念和术语,2023/8/31,19,1.5 虚应力 ij：,在s 上:njij=0；,ij=ij(k1)-ij(k2),10-1 几个基本概念和术语,在V内：ij,j=0,满足齐次静力方程。,2023/8/31,20,10-2 虚功方程,2.1虚功方程,在给定体力、面力和约束情况下，如果找到两种状态:,第一种状态：,在给定的体力 fi和面力，已知（找到）可能应力状态ij(k1),在V内：ij(k1)+fi=0;在s=s:,2023/8/31,21,第二种状态：,弹性体处于可能变形状态 ui(k2)、ij(k2),则第

7、一种状态外力在第二种状态可能位移作的外力虚功等于第一种状态可能应力在第二种状态可能应变上作的虚变形功。,虚功原理,10-2 虚功方程,在s=su:,2023/8/31,22,2.2虚功方程的证明：,10-2 虚功方程,2023/8/31,23,10-2 虚功方程,2023/8/31,24,代入虚功方程左端，得,并注意,虚功方程未涉及本构关系，所有在各种材料性质虚功方程成立。,10-2 虚功方程,则 We=Wi,2023/8/31,25,虚功方程虽然对两种不相干的可能状态成立，但一般应用是一种为真实状态，另一种为虚设可能状态（虚设状态）。,10-2 虚功方程,2023/8/31,26,10-3

8、功的互等定理,将虚功方程用于线弹性体可导出功的互等定理。同一弹性体处于两种真实状态。,第一种状态：,满足所有方程。,第二种状态：,满足所有方程。,2023/8/31,27,10-3 功的互等定理,根据虚功方程,第一种状态的外力在第二种状态的相应弹性位移上做功,第二种状态外力在第一种状态的相应弹性位移上做功,2023/8/31,28,10-3 功的互等定理,对于线弹性体本构关系,W12=W21,2023/8/31,29,10-3 功的互等定理,功的互等定理优点：可以避免求解物体内的应力、应变和位移场的复杂过程，而直接从整体变形的角度来处理问题。,第一种状态的外力在第二种状态的相应弹性位移上所做的

9、功等于第二种状态外力在第一种状态的相应弹性位移上所做的功。,第一状态：一对力P 作用在直杆的垂直方向，局部效应，在杆两端点伸长？,2023/8/31,30,10-3 功的互等定理,第二状态：让一对力Q 作用同一杆两端点，很易求得一对力Q引起杆横向缩短。,对两种状态应用功的互等定理 P=Q,Q第二状态引起的 易求:,2023/8/31,31,10-3 功的互等定理,2023/8/31,32,10-4 虚位移原理和最小势能原理,4.1虚位移原理,运用虚功原理，但一种状态为与真实外力平衡的状态，ij、fi、;而第二状态为可能变形状态，为真实状态位移的变分:,2023/8/31,33,10-4 虚位移

10、原理和最小势能原理,根据虚功方程，真实的外力与应力状态在虚设的齐次可能位移上做功,弹性体应力与外力处于平衡状态，对于任意虚设的齐次微小位移及应变，则外力在虚位移上做的虚功等于应力在虚应变上做的虚功 虚位移方程。,2023/8/31,34,10-4 虚位移原理和最小势能原理,代入原虚位移方程,将虚位移方程重新改写,2023/8/31,35,10-4 虚位移原理和最小势能原理,代入原虚位移方程,虚位移方程为平衡方程和力的边界条件的积分形式。,虚位移原理举例,2023/8/31,36,10-4 虚位移原理和最小势能原理,虚位移原理举例,图示受均布荷载q作用的等跨连续梁，EI为常数，中间支座为弹性支座

11、。试用虚位移原理写出梁的挠曲线方程和边界条件。,解：图示连续梁在荷载作用下，产生挠曲线w(x)、内力和支座反力。,设连续梁有虚位移 w，C点虚位移 wC。,2023/8/31,37,10-4 虚位移原理和最小势能原理,虚位移方程,利用对称性,在计算薄梁的内力虚功时，只考虑梁的正应力x 作的虚功。,2023/8/31,38,10-4 虚位移原理和最小势能原理,2023/8/31,39,10-4 虚位移原理和最小势能原理,分部积分，得,再积分，得,2023/8/31,40,10-4 虚位移原理和最小势能原理,虚位移方程为,2023/8/31,41,10-4 虚位移原理和最小势能原理,得,平衡微分方

12、程和力的边界条件（由虚位移方程得到）,2023/8/31,42,10-4 虚位移原理和最小势能原理,得,平衡微分方程和力的边界条件（由虚位移方程得到）,而,事先要求满足的位移边界条件,2023/8/31,43,10-4 虚位移原理和最小势能原理,4.2最小势能原理,1弹性体的总势能,(k)=U(k)(ui(k)+V(k)(ui(k)=(k)(ui(k),2023/8/31,44,10-4 虚位移原理和最小势能原理,(k)=U(k)(ui(k)+V(k)(ui(k)=(k)(ui(k),（1）和 给定;,（2）已将几何关系引入 ij=(ui,j+uj,i)/2;,（3）ui(k)为可能位移：在s

13、u上;,2023/8/31,45,10-4 虚位移原理和最小势能原理,2由(k)中寻求真实位移,ui(k)为可能位移,有无穷多。因此，与其对应的势能(k)也有无穷多。,要从(k)中找真实位移：,（1）=0,（2）引入本构关系,真实位移应满足的方程。,2023/8/31,46,10-4 虚位移原理和最小势能原理,取=0,得,引入本构关系,虚位移方程,2023/8/31,47,10-4 虚位移原理和最小势能原理,或,ui(k)为可能位移，同时满足本构方程。而=0，表明由ui(k)导出ij(k)满足静力方程，所以由=0 即为真解应满足的控制方程。,2023/8/31,48,10-4 虚位移原理和最小

14、势能原理,最小势能原理的表述：在位移满足几何方程和位移边界条件的前提下，如果由位移导出的相应应力还满足平衡微分方程和力的边界条件，则该位移必使势能 为驻值（极值）。如果可能位移使 的变分=0，则该位移相应应力必满足静力方程。=0等价与静力方程。,2023/8/31,49,10-4 虚位移原理和最小势能原理,作业：图示梁受荷载作用，试利用虚位移原理和最小势能原理导出梁的平衡微分方程和力的边界条件。,2023/8/31,50,10-4 虚位移原理和最小势能原理,最小势能原理举例,（1）已知图示桁架各杆EA相同，材料的弹性关系为=E，试用势能原理求各杆内力。,解：计算图示桁架的总势能,=U+V=(u

15、c、vc),应变能,2023/8/31,51,10-4 虚位移原理和最小势能原理,则，应变能为,2023/8/31,52,10-4 虚位移原理和最小势能原理,荷载势能为,=U+V=,由总势能 的变分=0，得,2023/8/31,53,10-4 虚位移原理和最小势能原理,解得,2023/8/31,54,10-4 虚位移原理和最小势能原理,各杆内力为,2023/8/31,55,10-4 虚位移原理和最小势能原理,(2)已知图示 a b 薄板（厚度 t=1）无体积力作用，试用最小势能原理求位移。,解：图示薄板受双向压缩作用，可猜应力和应变为常量。,所以，位移为x,y的线性式.,采用最小势能原理求位移

16、,2023/8/31,56,10-4 虚位移原理和最小势能原理,位移边界条件,x=0:u=0,y=0:v=0,考虑位移边界条件,取位移：,u=A1x,v=B1y,薄板的总势能,=U+V,2023/8/31,57,10-4 虚位移原理和最小势能原理,平面应力问题：,代入应变能表达式，得,2023/8/31,58,10-4 虚位移原理和最小势能原理,u=A1x,v=B1y,代入上式，得,2023/8/31,59,10-4 虚位移原理和最小势能原理,荷载势能V 为,=,2023/8/31,60,10-4 虚位移原理和最小势能原理,由总势能 的变分=0，得,解得,2023/8/31,61,10-4 虚

17、位移原理和最小势能原理,代回 u=A1 x,v=B1 y，得,求得应力为,x=-q1,y=-q2,xy=0,2023/8/31,62,10-5 虚应力原理和最小余能原理,5.1 虚应力原理,1虚应力方程,运用虚功原理，但第一种状态为真实变形状态，ui和 ij、fi、,第二状态为自平衡状态的可能应力（或真实应力的变分）ij；满足：,ij,j=0 在 V内,njij=0 在s上（在s 上无面力）,2023/8/31,63,10-5 虚应力原理和最小余能原理,ij在 su上产生 Xi（在su上有反力）,根据虚功方程,2虚应力方程表达,弹性体的应变与位移处于相容状态，对于任意虚设的齐次容许应力ij及位

18、移边界上的虚反力Xi，虚应力在应变上做的虚功等于虚反力在给定位移 上做的虚功。,2023/8/31,64,10-5 虚应力原理和最小余能原理,3虚应力原理举例,图示受均布荷载q作用的等跨连续梁，EI为常数。试用虚应力原理求梁的反力和弯矩。,解：连续梁在荷载q作用下，产生挠曲线、内力和支座反力。,根据平衡关系，知,2023/8/31,65,10-5 虚应力原理和最小余能原理,梁的弯矩为,设连续梁有虚反力和虚弯矩：,应用虚应力方程,2023/8/31,66,10-5 虚应力原理和最小余能原理,在 su 上无位移,2023/8/31,67,10-5 虚应力原理和最小余能原理,得,解得,2023/8/

19、31,68,10-5 虚应力原理和最小余能原理,5.2最小余能原理,已知变形体在体力 fi、面力 作用及在位移边界上有给定位移。,定义：由可能应力状态ij(k)表示,c(k)=U c(ij(k)+V c(Xi(k)=c(k)(ij(k)变形的总余能,1变形体的总余能c(k),2023/8/31,69,10-5 虚应力原理和最小余能原理,c(k)=U c(ij(k)+V c(Xi(k)=c(k)(ij(k)变形的总余能,而 Xi(k)=njij(k)在 su 上（位移边界上的反力）,应变余能,位移边界的余能,由ij(k)导出的在 su 边界上的反力,2023/8/31,70,10-5 虚应力原理

20、和最小余能原理,结构总余能 c(k)由可能应力ij(k)定义。,ij(k)满足 在V 内：ij(k)+fi=0,在 s=s:,（1）c=0,2由ij(k)定义的总余能中找出真实应力ij,（2）引入 关系式逆弹性关系,2023/8/31,71,10-5 虚应力原理和最小余能原理,c=0，即,虚应力方程,2023/8/31,72,10-5 虚应力原理和最小余能原理,由于ij 满足自平衡的应力状态:,ij,j=0 在 V 内,njij=0 在 s 上,c=0 为一个虚应力方程,则 ij(k)为可能变形状态，而ij(k)已满足静力方程，,由ij(k)导出的 ij(k)、ui(k)满足几何方程及位移边界

21、条件。,2023/8/31,73,10-5 虚应力原理和最小余能原理,因此，由 c=0 表明由ij(k)中找出真实应力ij。,3最小余能原理表述：,在应力满足平衡微分方程和应力边界条件的前提下，如果由应力导出的相应应变还满足相容条件，则该应力必使总余能 c 为极值；或可能应力使得总余能的变分 c=0,则由该应力导出相应位移必满足相容条件。,2023/8/31,74,10-5 虚应力原理和最小余能原理,最小余能原理举例,两跨连续梁受均布荷载 q作用的弯曲问题，试用最小余能原理求解。,解:,（1）确定梁的可能应力状态：,梁的弯曲问题可选取支座反力RB 为广义可能应力,梁的弯矩和其它支座反力可由RB

22、 表示。,2023/8/31,75,10-5 虚应力原理和最小余能原理,2确定结构的总余能 c(k),2023/8/31,76,10-5 虚应力原理和最小余能原理,c(k)=U c(R B)+V c(R B),3由余能的变分c=0 确定 RB,2023/8/31,77,10-5 虚应力原理和最小余能原理,即,积分，解得,2023/8/31,78,10-5 虚应力原理和最小余能原理,代入,可得梁的最后弯矩方程,2023/8/31,79,10-5 虚应力原理和最小余能原理,作业：利用虚应力原理和最小余能原理，求梁的反力和弯矩。,2023/8/31,80,10-5 虚应力原理和最小余能原理,4真实状

23、态总势能 总余能 c 关系：,虚功方程,所以真实状态=-c,对于真实状态的,2023/8/31,81,10-6 基于能量原理的近似解法,在前面几节介绍了几个最基本的能量原理，利用能量原理求问题的解，从理论上看是明确的步骤规范。如最小势能原理：（分两步）,1对于给定的外力和边界条件寻找满足几何方程和位移边界条件的ui(k)函数序列，并确定(k)。,2由=0寻求真解ui，即由ui(k)中找ui的控制方程，但由于问题求解域复杂性及约束的变化，利用能量原理求解析解也是无法实际进行的。但可由能量原理可以建立寻求问题近似解的有效途径。,2023/8/31,82,10-6 基于能量原理的近似解法,6.1基于

24、虚位移原理的近似解法,虚位移原理可以用来求问题的近似解法。,在给定的体积力、边界力和边界位移情况下，真实应力、应变和位移状态（它们满足所有方程）对于任意虚位移和虚应变满足虚位移方程。,虚位移原理近似解法的步骤,2023/8/31,83,10-6 基于能量原理的近似解法,（1）选取可能位移（近似解），包含若干待定系数；（满足位移边界方程）,（2）由可能位移求可能应变（满足几何方程）以及应力（应力一般不是可能应力），它们包含若干待定系数；,（3）设满足齐次边界条件的虚位移，并导出相应的虚应变；,2023/8/31,84,10-6 基于能量原理的近似解法,（4）将外力和包含若干待定系数的应力对虚位移

25、作功等于零虚位移方程（应力近似满足静力方程）；,（5）由虚位移方程得到确定若干待定系数的方程，并由方程解出待定系数，从而得位移的近似解。,举一个应用虚位移原理求问题的近似解的例题,2023/8/31,85,10-6 基于能量原理的近似解法,图示简支梁受均布荷载q作用,材料应力应变关系分别为=E，试确定梁挠曲线的近似解。,解：,（1）设梁的近似解为（包含若干待定系数）：,v=B1x(x-l)满足边界条件,2023/8/31,86,10-6 基于能量原理的近似解法,（2）由梁的近似解求可能应变以及应力（应力一般不是可能应力）：,（3）设满足齐次边界条件的虚位移，并导出对应的虚应变；,v=B1x(x

26、-l),=-v=-2 B1,2023/8/31,87,10-6 基于能量原理的近似解法,（4）应用虚位移方程：,（5）由虚位移方程解出待定系数B1，从而得位移的近似解。,2023/8/31,88,10-6 基于能量原理的近似解法,得,材料力学的精确解,2023/8/31,89,10-6 基于能量原理的近似解法,作业：利用虚位移原理的近似法求梁的弯矩。,2023/8/31,90,10-6 基于能量原理的近似解法,u(k)=u0+Amum,v(k)=v0+Bmvm,w(k)=w0+Cmwm,1选取可能位移(在给定的 条件下选可能位移),式中u0、um、v0、vm、w0、wm 均为已知连续可微函数，

27、,6.2基于最小势能原理的近似解法（Ritz 法）,2023/8/31,91,10-6 基于能量原理的近似解法,而 um=0、vm=0、wm=0、在 Su 上,可能位移中的 Am、Bm和Cm为待定系数。,在 Su 上,且u0、v0、w0满足Su 的位移边界条件,2023/8/31,92,10-6 基于能量原理的近似解法,2结构的总势能(k)及其变分(k),此时结构的总势能不是泛函了，而是Am、Bm、Cm的函数。,(k)=(k)(Am、Bm、Cm)=U(Am、Bm、Cm)+V(Am、Bm、Cm）,2023/8/31,93,10-6 基于能量原理的近似解法,3.利用=0求解方程,=0,由于Am、B

28、m、Cm的增量 Am、Bm、Cm的任意性，则=0,2023/8/31,94,10-6 基于能量原理的近似解法,需要求,对于线弹性体的应变能 U 为待定系数的二次式，荷载势能 V 为待定系数的一次式。,2023/8/31,95,10-6 基于能量原理的近似解法,2023/8/31,96,10-6 基于能量原理的近似解法,各向同性线弹性体有关Am、Bm和 Cm 的线性方程组（3m个方程）,2023/8/31,97,10-6 基于能量原理的近似解法,4Ritz 法在平面问题的应用,可能位移 u(k)=u0+Amum,v(k)=v0+Bmvm,u0、v0满足 Su 的位移边界条件,而 um=0、vm=

29、0 在Su 上。,平面应力和平面应变问题均不考虑位移分量w，而u、v为x、y的函数，体积力分量 fz=0，面力分量,2023/8/31,98,10-6 基于能量原理的近似解法,总势能 的变分=0，得,总势能=（Am,Bm）,2023/8/31,99,10-6 基于能量原理的近似解法,平面应力问题，取薄板厚度 t=1。其应变能为:,对于平面应变问题（取t=1）将上式中,替换,可得平面应变问题应变能 U 的表达式。,2023/8/31,100,10-6 基于能量原理的近似解法,Ritz 法在平面问题举行例,设有一无限长的薄板，上下两端固定，仅受竖向重力作用。求其位移解答。,2023/8/31,10

30、1,10-6 基于能量原理的近似解法,解：由于问题对y 轴对称，所以推论：,(1)选择可能位移：,设：,2023/8/31,102,10-6 基于能量原理的近似解法,满足上下两边的边界条件：,=U+V,2023/8/31,103,10-6 基于能量原理的近似解法,(2)计算应变能 U：,（取 y 轴两侧各1/2单位长度计算）,2023/8/31,104,10-6 基于能量原理的近似解法,(3)由=0 变分方程确定系数 Bk：,经推导得：,即,2023/8/31,105,10-6 基于能量原理的近似解法,或：,2023/8/31,106,10-6 基于能量原理的近似解法,得：,这是k个独立方程，

31、可求出k个待定系数 Bk。,解得：,2023/8/31,107,10-6 基于能量原理的近似解法,(4)位移解答,2023/8/31,108,10-6 基于能量原理的近似解法,值与级数项数的关系,作业：设位移的近似解为 u=0,v=B1 y(y-b)，求其位移解答。,2023/8/31,109,10-6 基于能量原理的近似解法,最小势能原理与 Ritz 法的比较,2023/8/31,110,10-6 基于能量原理的近似解法,在最小势能原理中，由可能位移ui(k)定义的总势能(k)，并由=0寻求真解ui。,也可改写为静力方程的积分形式,等价于静力方程。,6.3伽辽金法（1915年）,而=0 本身

32、表示为虚位移方程：,2023/8/31,111,10-6 基于能量原理的近似解法,伽辽金法步骤：,1设定满足强约束条件的可能位移ui(k),ui(k)需要满足的强约束条件：,在Su上,由u(k)、v(k)、w(k)导出的应力ij(k),在S 上：,可能位移ui(k)满足给定 和 条件，即满足所有边界条件强约束条件。,2023/8/31,112,10-6 基于能量原理的近似解法,2由结构总势能 的变分等于零导出求解方程,如果所设可能位移ui(k)的形式与Ritz法一样,u(k)=u0+Amum,v(k)=v0+Bmvm,w(k)=w0+Cmum,但可能位移ui(k)满足强约束条件。由=0 得,这

33、里 u=um Am,v=vmBm,w=wm Cm,2023/8/31,113,10-6 基于能量原理的近似解法,并注意Am、Bm、Cm的任意性以及,由 可导出三组方程,3m个方程,2023/8/31,114,10-6 基于最小势能原理上的近似解法,3伽辽金法在平面问题例题,与Ritz 法类似，不考虑w,而 u、v为x、y函数。伽辽金法在平面问题的求解方程为2m个。,对于平面应力问题为,2023/8/31,115,10-6 基于能量原理的近似解法,对于平面应变问题将上式中,例.已知：2a b薄板（厚度t=1），无体力作用。边界条件:在x=a和 y=0：u=0,v=0;在y=b：u=0,全部边界为

34、位移边界条件，无力的边界条件。,2023/8/31,116,10-6 基于能量原理的近似解法,解：不管采用Ritz法或伽辽金法，选取的近似位移场首先要为可能位移：,u(k)=u0+Amum,v(k)=v0+Bmvm,当取m=1时，u=u0+A1u1,v=v0+B1v1,2023/8/31,117,10-6 基于能量原理的近似解法,根据边界条件,可选 u0=0，,而u1 和v1 在所有边界上为零。,取,则,2023/8/31,118,10-6 基于能量原理的近似解法,由于全部边界均为位移边界，可能位移无需要求其它约束条件。可利用伽辽金法求A1、B1。,由于体力为零，则用伽辽金法的求解方程（2个）

35、：,求出A1 和 B1,2023/8/31,119,10-6 基于能量原理的近似解法,，,A1 和 B1 求出后可求应力的近似解。,2023/8/31,120,10-6 基于能量原理的近似解法,作业：1.试写出伽辽金法在梁弯曲问题的求解方程。2.利用伽辽金法求图示简支梁的近似解，设梁挠度的近似解为 v=B1 sinx/l。,2023/8/31,121,10-6 基于能量原理的近似法,6.4 基于最小余能原理的近似解法,c(k)=U c(ij(k)+V c(Xi(k)=c(k)(ij(k)变形的总余能,且 Xi(k)=njij(k)在 su 上（位移边界上的反力）,由可能应力 ij(k)定义的总

36、余能中找出真实应力ij,即要求 c=0,由可能应力ij(k)定义的余能,2023/8/31,122,10-6 基于能量原理的近似法,最小余能原理的近似法步骤与最小势能原理的 Ritz 法类似。,最小余能原理的近似法步骤：,1选取可能应力ij(k)为自变函数：,ij(k)=ij0+A ijmijm,其中 ij0 和 ijm 为已知函数且在域内分别满足平衡微分方程和齐次平衡微分方程。,2023/8/31,123,10-6 基于能量原理的近似法,ij0 在 s 上满足力的边界条件。ijm 在 s 上满足力的齐次边界条件。,而A ijm 为待定系数（共 6m个）。,2确定结构的总余能 c(k),c(k

37、)=U c(A ijm)+V c(A ijm)=c(k)(A ijm)总余能为多元函数,ij(k)=ij0+A ijmijm,2023/8/31,124,10-6 基于能量原理的近似法,3由余能的变分c=0 导出近似法的求解Aijm 方程,由于 A ijm 任意性,得,（共 6m 个方程）,求解方程。这6m个求解方程为线性方程组，可解出6m个待定系数Aijm。,2023/8/31,125,题10-1 图示结构各杆等截面杆，截面面积为A,结点C承受荷载P作用,材料应力应变关系分别为（1）=E，(2)=E 1/2。试计算结构的应变能U 和应变余能Uc。,能量法习题,2023/8/31,126,题1

38、0-2 分别利用虚位移原理、最小势能原理、虚应力原理和最小余能原理求解图示桁架的内力。已知桁架各杆 EA 相同，材料的弹性关系为=E。,2023/8/31,127,题10-4 利用最小余能原理求左图示梁的弯矩。,题10-3 左图示梁受荷载作用，试利用虚位移原理 或最小势能原理导出梁的平衡微分方程和力的边界条件。,2023/8/31,128,（1）悬臂梁受两个集中力 P 作用。,（2）简支梁受均布荷载 q 作用,设：v=B1x(x-l)+B2x2(x-l)。,题10-5 利用虚位移原理的近似法或Ritz 法求解图示梁的挠曲线。,2023/8/31,129,设位移的近似解为 u=0,v=B1 y(y-b)，求其位移解答。,题10-6 设有一无限长的薄板，上下两端固定，仅受竖向重力作用。利用Ritz 法求其位移解答。,2023/8/31,130,题10-7 1.试写出伽辽金法在梁弯曲问题的求解方程。2.利用伽辽金法求图示简支梁的近似解，设梁挠度的近似解为 v=B1 sin（x/l）。,