结构力学之能量原理

例1 用瑞利—里兹法求图示简支梁的挠度和弯矩。
该题材料力学已有精确 解，在梁中点挠度
Fp x
vmax4Fp8 E l3 I0.020F 8E pl33I
l/2
l/2
y
中点弯矩 MmaxF4pl 0.25Fpl
Fp x
l/2
l/2
y
解：设该简支梁的挠曲线（几何可能位移）为
ix
v ai
1,3,5
sin
e
e
内力（弯矩）在单独在荷载作用下基本结构的变形（曲率）
上所做的内力虚功。
而当Zi=1时基本结构的外力（r1i、 r2i …… rni ）在单独在 荷载作用下基本结构的变形上所做的外力虚功为0。所以
根据虚功方程得
E v i I v p d xM ip d U x p iW p i0
0
81 4 EI 2l 3
a3
Fp
0
解之得
a1
2Fpl 3
4EI
,
a3
2Fpl3
81 4EI
则 v82F 1p 4lE 3 I8s1inlxsi3 nlx
M (x) EvI2 9 F p 2 l 9silx nsi3 n lx
vmax0.0207F8EplI3 ，比精确解少0.24%
该原理说明，如果内力满足全部的静力平衡条件， 而且还能使结构的余能为驻值，则与此内力相应的变形必 然满足变形协调条件，即余能驻值条件与变形协调条件是 等价的。
可以证明：超静定结构中，在同时满足静力平衡方 程、几何方程和物理方程的解具有唯一性的情况下，结构 的真实内力不仅使余能为驻值，而且该驻值一定为极小值。 这就是最小余能原理。
e
e
或
n
j 1 e
E v iI v jd Z x je
E v iI v p d xF p i 0
p
当Zi=1时的基本结构外力，在基本结构单独在荷载作用下 的变形上所做得的虚功为0，而在基本结构单独在荷载作
用下的外力在Zi=1时的基本结构的变形上所做的虚功为
Rip1 Fpi
p
根据功的互等定理有 Rip Fpi 0
n
v vi Zi vp i1
上式中，v i 为基本结构由于Zi=1时引起的各杆任一截面的
位移方程。vp为基本结构在荷载作用下任一截面的位移方 程。
与广义荷载Fp相应的广义位移也可表示为
n
i Zi p i1
n
iwenku.baidu.comZi p i1
上式中， i 为基本结构由于Zi=1时引起的与广义荷载相应
的广义位移。△p为基本结构在荷载作用下引起的与广义
p
即 Rip Fpi
p n
由上述讨论可得 rij Zj Rip 0 j1
这就是杆系结构的位移法典型方程。
多提意见与建议
谢谢！
作业：
3.2瑞利—里兹法(Rayleigh-Ritz Method)
建立在能量原理基础之上的解题方法是一种精确方法， 但在精确解难以求得或不能求得的许多工程实际问题中， 能量原理又能为我们提供一种求近似解的有效途径。瑞 利—里兹法就是其中之一。
即应力 应变曲线中OAB所围的面积。
1.2 应变余能密度
应变余能密度定义 :单位体积内 dFN 的应变余能称为应变余能密度。
仍以简单拉伸杆件为例，应变余能为
U*W* dN F （12-3）
FN dx
图(a)简单拉伸曲线
根据应变余能密度的定义，则应变余能密度u*N为
u* NU V*A1dxdN Fd（12-4）
E*p Fp
（12-18）
p
上式中p为荷载的序号，为Fp方向上的位移。
2 势能驻值原理 势能驻值原理：在所有几何可能的位移状态中，真实的 位移应使结构势能为驻值。
这一能量原理说明，如果位移满足全部的变形协调条件， 而且还能使势能为驻值，则与此位移相应的内力必然满 足全部的静力平衡条件。即说明势能驻值条件与平衡条 件是等价的。
U* e
1 2
F EN2AkG FQ 2A M E2Idx（12-20）
E*C为结构的支座位移余能，或称给定边界位移余能， 即在支座位移c上相应支座反力R所做的虚功总和的负值。
EC* Rc
c
2 余能驻值原理
（12-21）
超静定杆件结构的余能驻值原理可表述如下：在所有静力 可能内力中，真实的内力应使结构的余能为驻值。
显然有
u1
EA
FN
u 1 GA k
F
Q
u 1 EA M
（12-8）
设：杆截面形心的轴向位移为u，横向位移为v，截面的转
角为。则几何方程为
du dx
dv dx
d dx
将上式代入（12-7）式得
（12-9）
u 11 2E d d A u 2 x1 2G k d d A x v 21 2E Id d x 2（12-10）
(2)仅取级数的前两项 ，则 va1sinlxa3sin 3lx
上式中没有取
a2
sin
2x
l
项，是因为在Fp的作用下，
内力和变形都是对称的，而此项在中点处v=0，变形
是反对称的。
U 0 lE 2v I 2d x E 22I0 l a l2 2a 1silx n 3 2 l22si3 ln x 2dx
§12.2势能原理
1 势能的定义
杆件结构的势能Ep定义为
Ep UE*p
（12-16）
上式中，U为杆件结构的应变能，对于刚架而言，通常仅
考虑弯曲应变能，则
U e
1 2
EIdd2xv2 2dx
（12-17）
上式中e为结构中杆件的排序号。E*p为结构的荷载势能，
通常以结构未变形前的荷载位置为起始位置，则
结构力学
第12章 能量原理
主要内容
1 杆件的应变能及应变余能计算 2 结构势能定义及势能原理 3 结构余能定义及余能原理
能量的概念大家早已了解，在第六章分析静定结构 的位移计算中，曾介绍了虚功方程的两种应用：虚设单 位力求位移和虚设单位位移求未知力。在本章中将介绍 基于能量原理基础上的解题方法。
§12.1 杆件的应变能及应变余能计算
几何可能位移
如果变形体的应变、、与位移u、v、满足几何方
程，而且在结点处满足位移连接条件，在边界上能与约束 几何相容。则此种位移称为几何可能位移。
在变形体上，这种几何可能位移有无穷组，但只有同 时能满足静力平衡条件的那一组才是真实的解答。
结构的总势能是一个泛函，对于稳定的平衡问题而言， 按位移法求解时，就归结为求泛函的极值问题。瑞利—里 兹法就是建立在泛函求极值基础之上的一种求近似解的方 法。下面举例说明。
∵ l sin2ixdxl (i1,2)
0
l
2
0lsin ilxsinjlxdx0 (ij)
∴ U0 l E 2Iv2dx 44lE 3 Ia128a13 2
E * p F pvx l F pa 1 a 3 2
由势能驻值条件 Ep 0i 1,3得
ai
4 EI 2l 3
a1
Fp
当杆件同时承受拉伸、剪切和弯曲时，其杆件的应
变余能密度为
FN
FQ
M
u1* dF N dF Q dM（12-13）
0
0
0
对于线弹性材料，用类似的方法，可以得
u1 *2E 1F A N 22G k F A Q 221 EM I2 （12-14）
一根杆的应变余能为
U *u 1 * d x 1 2 E 1F A N 2G kF Q 2 A E 1M I2 d x （12-15）
在介绍瑞利—里兹法之前，先介绍两个基本概念：
静力可能内力 对于变形体而言，如果它的内力与外力满足全部
的静力平衡条件，即满足杆件的平衡微分方程，而且 在边界上和结点处满足力的平衡条件，则此种内力称 为静力可能内力。
对于静定结构而言，静力可能的内力是唯一的， 而对于超静定结构而言，静力可能的内力不是唯一的。
l
这个函数不仅满足简支梁的两端的位移边界条件，而
且满足两端力的边界条件：
v(0)v(l)0; v(0)v(l)0
(1)仅取级数的首项，则
v
a1
sinx
l
∵ U 0 lE 2vI 2d x E 22I0 l a l2 2 2s2 in lxd x4 4 lE 3a 1 2 I
E* pFpvxl Fpa1 2
用u1表示。
则当杆件同时承受拉伸、剪切和弯曲时，其杆件的应 变能密度为
u1dAdAddA
0
0
A0
即 u1 FNd FQdMd （12-6）
0
0
0
对于线弹性材料，有FN=EA. ， FQ=GA. /k（k为截面形状 系数），M=EI . 。则
u11 2EA 21 2G kA 21 2EI2 （12-7）
1 应变能密度和应变余能密度
1.1 应变能密度
应变能密度定义 :单位体积内的应变能称为应变能密度 例如简单拉伸杆件，取出dx微段，其拉伸曲线如图(a)
所示
FN
dx
应变能U为
UW FNd
（12-1）
C
A
d 图(a)简单拉伸曲线
O d
B
根据应变能密度的定义，则应变能密度为 图(b) 应力 应变曲线
uNU VA 1 dx F Ndd （12-2）
∴ EpUE* p44lE 3 aI1 2Fpa1
由势能驻值条件 E p 0 得
24lE3 Ia1
Fp
a1
0
即
a1
2Fpl 3
4 EI
则
v2 F 4p E l3s Iiln x0.0
2F 0 pl35 s3 in x
EI l
vmax0.0205F3EplI3 ，比精确解少1.44%，
Mmax0.202Fp6l，比精确解少19%。
荷载相应的广义位移。
则结构的势能为
E p e
1 2 E i n 1 Iv iZ i v p 2 d x pF p i n 1 iZ i p
根据势能驻值条件得
Ep 0 (i1,2n) Zi
n
即 e E Ij 1vjZjvp vid xpF p i 0
或
3 余能驻值原理应用
3.1 利用余能驻值原理推导力法典型方程
设力法基本未知量向量为{X}={X1 X2 … … Xn}T ，在力法 基本结构中，各杆任一截面的内力可表示为
F N
n
FN
i
X
i
F Np
i1
n
F Q F Q i X i F Qp
Mma x 0.22F5pl，比精确解少10%
误差仍较大，但位移和弯矩的精度均有所提高，随 着级数项数增加，位移和弯矩都将趋于精确解。
§12.3余能原理
1 余能的定义 ：杆件结构的余能EC定义为
EC U*EC *
（12-19）
上式中，U*为杆件结构的应变余能，对于线弹性材料而 言，杆件结构的应变余能为
一根杆的应变能为
U u 1d x 1 2 E d d A 2 u xG k d d A x v 2E Id d 2 x d（1x2-11）
当忽略较小的剪切变形后,
0,
dv dx
,
则
Uu1d x1 2EA d du x2EId d2v 2x2dx（12-12）
定义：单位杆长上的应变余能为杆件的应变余能密度， 用u*1表示。
即应力 应变曲线中OAC所围的面积。
C
A
d
O
B
图(b) 应力 应变曲线
对于线弹性材料，=E.有，则
uN0d12E212E2 u*N
（12-5）
2 杆件的应变能和应变余能
象纯拉伸一样，当杆件处于纯剪切和纯弯曲时，其应 变能密度分别为
uQ d ; uM d
0
0
定义：单位杆长上的应变能为杆件的应变能密度，
根据虚功方程U ij =W ij得
rij EIvivjdx e
或
n
j 1 e
E v iI v jd Z x je
E v iI v p d xF p i 0
p
又因为 vp p 为单独在荷载作用下的基本结构的变形
（曲率）。
E v iI v p d x M ip d U xp代i 表了当Zi=1时基本结构的
n
j 1 e
E v iI v jd Z x je
E v iI v p d xF p i 0
p
因为，EIviMi为Zi=1时的基本结构的内力（弯矩），
vj j 为Zj=1时的基本结构变形（曲率）。则
EIvivjdxUij 为基本结构Zi=1时的内力（弯矩）
在eZj=1时的变形（曲率）上所做的内力虚功（虚应变能）。 而当Zi=1时基本结构的外力（r1i、 r2i …… rni ）在Zj=1时的 位移上所做的外力虚功为W ij=rij1=rij。
可以证明，在小变形、线弹性的稳定平衡问题中， 满足几何方程、物理方程和静力平衡方程的解是唯一的。 此时真实的位移不仅使势能取得极值，而且该极值为极 小值。这就是最小势能原理。
3势能驻值原理应用 3.1 利用势能驻值原理推导位移法典型方程 设：位移法的基本未知量向量为{Z} ={Z1 Z2 …… Zn}T 在位移法基本结构中，各杆任一截面的位移方程可表示为