哈密顿正则方程例题
19哈密顿正则方程
——哈密顿正则方程
相应的广义动量, 坐标叫做正则变量, 它们组成的2s维 空间叫相空间, 一组数值对应相空间中一点, 叫相点.
一维弹簧振子的运动 • 哈密顿量 H=Ep+Ek
L T V
L pi qi
H Ek Ep mx2 / 2 kx2 / 2 p2 1 kx2 2m 2
哈密顿正则方程:
p
H
p2 cos mr 2 sin 3
p 0
r H pr pr m
H p
p mr 2
H
p
p mr 2 sin 2
m r m r2 m r sin 2 2 a
r2
pr m r p m r2θ
p m r2 sin 2 θ
C m r 2 sin 2
m r m r2 C 2 a 0 m r 3 sin 2 r 2
哈密顿动力学的优点之一是便于量子化.另一个优 点在变量的变换中比较自由:拉格朗日动力学采用的 变量广义坐标和广义动量并不对等, 只能对广义坐标进 行变换, 而广义速度也随之而变. 哈密顿动力学采用的 变量坐标和动量是完全对等的,不仅可以对广义坐标进 行变换,而且可以坐标和动量一起变换, 这个到下面正 则变换时进一步分析.
2
V
(r)
1 2
k(r
r0
)2
1 2
2
(r
r0
)2
L = T-V
1 2 (m1
m2 )(x2
y2
z2 )
1 2
(r2
r2
r2
sin
2
2
)
1 2
2
(r
r0
)2
p x (m1 m2 )x p y (m1 m2 ) y
5.5哈密顿正则方程作业
1哈哈密密顿顿正正则则方方程程作作业业1.如果哈密顿函数L qp H s1-=∑=ααα 能表示成坐标αq 和动量αp 的函数,则无论H 是(a)明显不包含时间变量t,还是(b)明显包含时间变量t,哈密顿正则方程均为: αα∂∂=p H q ;αα∂∂-=q H p )s 2,1( =α 1.证明: (a)H 明显不包含时间变量t 时:对L qp H s1-=∑=ααα 微分,得: ∑∑∑∑=ααα=ααα=ααα=ααα∂∂-∂∂-+=s 1s1s 1s 1q d q L dq q L dp q q d p dH 而αα=∂∂p q L ; ααα=⎪⎪⎭⎫ ⎝⎛∂∂=∂∂p qL dt d q L (利用拉氏方程) 于是:∑∑=ααα=ααα-=s 1s 1dq p dp qdH 又哈密顿函数L qp H s1-=∑=ααα 能表示成广义坐标αq 和广义动量αp 的函数,即:∑∑=ααα=ααα∂∂+∂∂=s 1s1dp p H dq q H dH 比较上面两式可得:αα∂∂=p H q ;αα∂∂-=q H p )s 2,1( =α (b)H 明显包含时间变量t 时:对L qp H s1-=∑=ααα 微分,得: dt t L q d q L dq q L dp q q d p dH s 1s1s 1s 1∂∂-∂∂-∂∂-+=∑∑∑∑=ααα=ααα=ααα=ααα2 而αα=∂∂p q L ; ααα=⎪⎪⎭⎫ ⎝⎛∂∂=∂∂p qL dt d q L (利用拉氏方程) 于是:dt t L dq p dp qdH s 1s 1∂∂--=∑∑=ααα=ααα 又哈密顿函数L q p H s1-=∑=ααα 能表示成广义坐标αq 、广义动量αp 和时间t 的函数,即:dt t H dp p H dq q H dH s 1s1∂∂+∂∂+∂∂=∑∑=ααα=ααα 比较上面两式可得:αα∂∂=p H q ;αα∂∂-=q H p )s 2,1( =α,且t L t H ∂∂-=∂∂。
Chapter5-分析力学05-哈密顿正则方程
主讲教师:邱晓燕
西南大学-物理科学与技术学院 理论力学-5.5 哈密顿正则方程
2.循环积分 H ( p , q , t )
H 中不显含
q , p
循环积分 (亦称 广义 动量积分)
1. 一个循环坐标,就有一个循环积分。从而产生一个对应的 广义动量的守恒量.
先考虑两个变量函数的勒让德变换 f f x, y
df udx vdy
f f u ,v x y
若选择u.y 做为独立变量
x xu, y , v vu, y
f u, y f xu, y , y
西南大学-物理科学与技术学院 理论力学-5.5 哈密顿正则方程 主讲教师:邱晓燕
p 0 p C
2 c a 2 m r mr 3 2 2 0 mr sin r 2 d c cos 2 mr dt mr 2 sin 3
电子作平面运动, 取 0的平面, c 0
a 2 mr mr 2 0 r d 0 mr 2 dt
H L pr r p p
L 2 2 p mr sin
L pr mr r L 2 p mr
H L pr r p p p 1 a 2 (p r 2 2 ) 2 2m r r sin r
西南大学-物理科学与技术学院 理论力学-5.5 哈密顿正则方程
主讲教师:邱晓燕
本节习题
P366: 5.20, 5.23
西南大学-物理科学与技术学院 理论力学-5.5 哈密顿正则方程
主讲教师:邱晓燕
哈密顿正则方程
=
V
(r0
)+
∂V ∂r
r0 (r − r0 )
+ 1 ∂2V 2 ∂r 2
r0 (r − r0 )2 + ⋅⋅⋅⋅⋅⋅
r0
r
V
(r)
=
1 2
k
r2
r = r − r0
L
=T
−V
=
1 2
m ( x& 2
+
y& 2
+
z&
2
)
−
1 2
k
r
2
+
1 2
μ (r&2
+
r
2θ& 2
+
r
2ϕ&
2
sin
2
θ
)
L
=T
=
1 ml2θ&2 (1+ cos2 θ ) + mgl sinθ
4
= E0
= mgl
θ& θ = 0 =
2g l
y m1
m1x1 + m2x2 = 0
坐标数 约束数
3 x1 = −x2 2
m2 θ
自由度数 1
x
取如图所示 θ 为广义坐标
yc
=
l 2
sin θ
y
y& c
=
l 2
θ&
cos
θ
yc
根据柯尼西定理
T
=
1 2
2my&c2
+
1 2
I cθ& 2
T = 1 ml2θ&2 (1+ cos2 θ )
哈密顿正则方程
哈密顿正则方程
哈密顿正则方程为(1)式中H称为哈密顿函数,是广义动量pi和广义坐标qi及时间t的函数。
H由式(2)确定。
括号外边的角标表示式中的妜i应该用N个方程pi= 解出N 个妜i为(E1,E2,…,EN;q1,q2,…,qN;t)的N 个函数,然后代入式(2)就得到哈密顿函数H。
对于直角坐标变换到广义坐标的变换式虽然显含时间t,但是动能的表示式不明显地包含t,此时H=T2-T0+V,
式中T2和T0可说明如下:用(E1,E2,…,EN;q1,q2,…,qN;t)表示的动能式T=T2+T1+T0,式中T2、T1和T0分别表示广义动量的二次齐次式、一次齐次式和不含广义动量的项。
如果直角坐标变换到广义坐标的变换式不显含t,势函数V也不显含t,则T=T2,H=T+V。
即对于保守系统,哈密顿函数是系统总机械能用广义动量表示的公式。
正则方程式(1)是2N个一阶微分方程组,而拉格朗日方程是N个二阶微分方程组,都只适用于完整系统(见约束)的动力学方程组。
由于式(1)的左边不再有变数q和p的导数,所以方程(1)成为如下形式的方程组保守系统的正则方程在天体力学和经典统计力学中有重要的应用。
在天体力学中从可解的二体问题出发,逐渐添加其他星球的引力,可以把所用的哈密顿函数H,从简单改变成较复杂的H┡。
这是天体力学中的摄动法,用来解决考虑太阳和各种行星、卫星的引力作用下的行星运动,由此可制定行星和月球的星历表,在统计力学中的刘维定理
就是应用正则方程推导出来的。
分析力学第七章正则方程
知 必须满足条件:
由此得出重要推论:
当不显含t时, 为运动常数的充要条件是:
3. 泊松定理
如果函数
和函数
分,则函数[f , g]也是正则方程的初积分。
证:由于是f和g正则方程的初积分,得
是正则方程的两个初积
由雅克比恒等式: 得 于是有 即得到:
因此[f,g]=C也是正则方程的初积分.
泊松定理指出: 由正则方程的两个已知的初积分, 可不断地求出新的初 积分.
那么有
;于是得到:
(即在该四种正则变换中哈密顿量保
持不变).
此时正则变换条件变为下列形式:
。
例1.寻求常数 ,使变换
解:由于此变换不显t,有
是正则变换。
即
, 由于q的任意性,得
因此有变换:
该变换被彭家莱应用于天体力学中
例2. 证明变换 关的四类母函数。 解:
是正则的,并求出与该变换相
因此该变换是正则的。其母函数为:
,其中
是n+1个任意常数。
另外,如果我们已知
,其中
是n+1个任意常数。同样可以得到哈密顿—雅克比偏微
分方程:
——这是哈密顿在当时推证所用的方法。 利用哈密顿—雅克比方程求出
---这样就能得到正则方程的全部积分。
由
及哈密顿正则方程
若力学体系的哈密顿函数H中不显函时间t,即 (h是积分常数)。
;则
当约束又是稳定的,则动能可表示为
2n个代数方程是相互独立的,所以可以解出逆变换为:
若通过变量的变换,使得正则方程的形式保持不变,即:
我们把这种变换叫做正则变换。 当取第二类母函数 则正则变换的条件: 变为:
令
ቤተ መጻሕፍቲ ባይዱ
5.5 哈密顿正则方程ppt课件
H
H ( p,q,t)
dH dt
s 1
H p
p
H q
q
H t
正则方程
s 1
H p
H q
H q
H p
H t
dH H dt t
H t
若H不显含t, 则H守恒。
8
5.5.3 能量积分与循环积分
循环积分
若H不显含某个广义坐标q ,则根据正则方程
可得:
p
H q
0
p
常数
即相应的广义动量守恒。
16
5.5 哈密顿正则方程
分析:四、最终的运动微分方程
例题 2
由(5')(8')可得
p mr2 h(常数)
(11)
即在垂直于运动平面方向上的动量矩守恒。
对(11)求导即得到横向运动微分方程。
由(4')(7')(11)可得
m(r
r2 )
r2
(12)
即径向运动微分方程
17
5.5 哈密顿正则方程
18
(1)(2)联立可得
x k x 0 m
即一维弹簧振子的 运动微分方程
11
5.5 哈密顿正则方程
例题 2
例 2 用哈密顿正则方程建立质点在有心力势场
V
(r)
r
中的运动微分方程.
解:采用球坐标 (r, ,) 描述 位矢 r rer
质点的速度
v
d dt
(rer
)
rer
re
r sine
拉氏函数 L T V 1 m(r2 r22 r2 sin2 2 )
注:这里的能量积分与循环积分,与从拉氏函数得到的 能量积分和循环积分,结果是一样的。
哈密顿正则方程
mx
x px
m
.
10
5.5 哈密顿正则方程
哈密顿函数
HpxxLxpx m
px
px 1mpx21kx2 m 2 m 2
px2 1 k x2 2m 2
.
正则方程
例题 1
x
H px
px m
(1)
p x
H x
kx
(2)
(1)(2)联立可得
x k x 0 m
即一维弹簧振子的 运动微分方程
由(4')(7')(11)可得 m(rr2)r2
(12)
即径向运动微分方程
.
17
5.5 哈密顿正则方程
.
18
例题 2
由(6)(9)可得 pm2rsi2nC(常)数
(10)
另一方面计算可得 J z (r m v )e z m 2sr 2 in
可见(10)代表沿z方向的动量矩守恒。 但是z轴的方向是任意选择的,故沿任何方向都 有动量矩守恒,从而系统动量矩守恒,即
J C(常矢量)
.
14
5.5 哈密顿正则方程
例题 2
(2) (3)
p r
H r
p
H
p
H
p2 mr 3
p2
mr 3 sin 2
p2 cos mr 2 sin 3
0
r2
.
(4) (5) (6)
r
H pr H p H p
pr mp mr 2ຫໍສະໝຸດ p mr 2 sin2
(7)
(8)
(9)
13
5.5 哈密顿正则方程
分析:一、动量矩守恒
广 义 动
pr
哈密顿正则方程
1
哈密顿正则方程
例1-13 半径为a 的光滑圆形金属圈,以匀角速ω绕铅直方向的轴z
转动,圈上套有一质量为m 的小环。
初始时,小环自圆圈的最高点无初
速地沿圆环下滑,试求当与圆圈中心的联线和铅直向上的z 轴成θ角时,
小环的运动方程。
解:这是一个非定常的完整系统,主动力有势。
应用哈密顿正则方程
求解,取广义坐标为θ。
系统的功能T 为
)sin (2
122222θωθa a m T += 系统的势能为
θcos mga V =
拉格朗日函数L 为 θθωθcos )sin (2122222mga a a m L -+= 为求哈密顿函数H ,先求广义动量p θ
θθ
θ 2ma L p =∂∂= 哈密顿函数H 为 θθωθθθθθcos sin 2122242222mga a a m p a m ma p p L p H +⎪⎪⎭⎫ ⎝⎛--=-=
正则方程为
θθθθp H p H -=∂∂=∂∂, 得到
θθ =p ma 21 (1) θθθθωp
mga ma -=--sin cos sin 22
(2)
2 由式(1)(2)得到
θ
θθθω 222sin cos sin a ga a -=-- 或者写成
θθθωθsin cos sin 2g a a += 哈密顿函数H 中不显含时间t ,有广义能量积分,即
h g a a '=+-θθωθcos sin 2
121222。
正则方程
则根据正则方程可得到 pα = pα (qα , pα , t ) = aα = 常量 这就是与可遗坐标 qα 对应的广义动量积分, 物理 意义与前相同. 如果存在可遗坐标 qα , 必然存在与它对应的 广 义动量守恒, pα = aα , 则哈密顿函数可 表示 为
H = ( q1 ,, qα −1 , qα +1 ,, qs , p1 ,, pα −1 , aα , pα +1 ,, ps , t ) .
α 看做与 qα 无关的独立变量 , 并 改用 X α 如果把 q
表示, 则
L = L(qα , X α , t )
拉格朗日方程就可变成以 qα 和 X α 为独立变量的一 阶微分方程组:
∂L d ∂L − =0 dt ∂X α ∂qα dqα = X α dt α = 1,2, , s
qα =
Cα 是积分常数.
∫
∂H dt + Cα ∂aα
虽 然 系统存在 广 义动量 积 分 , 但是 , 在拉格 朗日表述中 , 可 遗 坐标 的出现 不会 导 致 广 义 速 度
α 为常数 , 因此 拉格朗日函数中 仍包 含 s 个独立 q
变量,
1 ,, q s , t ) L = L(q1 ,, qα +1 ,qs , q
但 这 种 处 理方法并 没 有在理论 上带来实质性 的 好 处. 计算机数值计算即如此. 一、勒让德变换 两个变量的勒让德变换: 设函数 f = f ( x, y ) , 该函数的全微分为 ∂f ∂f df = dx + dy ∂x ∂y 令u =
∂f ∂x
df = udx +
∂f dy ∂y
由于 udx = d(ux ) − xdu 所以
用球坐标写出球摆的哈密顿正则方程 泊松括号
用球坐标写出球摆的哈密顿正则方程泊松括号哈密顿力学一般是指以哈密顿函数为主函数的力学体系。
哈密顿函数一个罗斯函数的一个显然的推广:将朗格朗日力学中的所有的qk,q˙k换成qk,pk即可,事实上,这也正是数学上定义的勒让德变换。
将L中显含的时间t视作参数不参与变换,我们有:dL(qk,q˙k)=∑k∂L∂qkdqk+∑k∂L∂q˙kdq˙k按照动量定义式pk=∂L∂q˙k和拉格朗日方程:p˙k=∂L∂qk带入上面的式子:dL(qk,q˙k)=∑kp˙kdqk+∑kpkdq˙k我们的目标是让上面的式子不含dq˙k,所以利用pkdq˙k=d(pkq˙k)−q˙kdpk带入上面的式子:dL(qk,q˙k)=∑kp˙kdqk+∑kd(pkq˙k)−∑kq˙kdpk简单的合并一下得到:d(∑kpkq˙k−L)=−∑kp˙kdqk+∑kq˙kdpk那么根据全微分的定义,被微分的函数是qk,pk的函数,那么定义H=∑kpkq˙k−L就是很自然的。
对于主函数为H(pk,qk,t),写出它的全微分,有:dH=−∑kp˙kdqk+∑kq˙kdpk所以:∂H∂qk=−p˙k,∂H∂pk=q˙k,左边的方程也称作哈密顿方程,他与前面的拉格朗日-欧拉方程的信息完全一样。
也就是,哈密顿力学给出了力学的另一个表述形式。
启动哈密顿力学有一些绝妙的好处,首先是其上可以定义一个代数(事实上,无歧义来说,代数是指的一个线性代数,而我们这里事实上是一个李代数,如果你没有学过高等代数,请忽略。
)考虑一个力学量,在(pk,qk)张成的空间中,哈密顿已知,为了探讨这个量,我们考虑做他的全微分:df(pk,qk,t)=∑k(∂f∂pkdpk+∂f∂qkdqk)+∂f∂tdt,由于系统的哈密顿已知,系统的全部信息都已经被获取,事实上任何力学量都只依赖于时间,带入哈密顿方程,消去dqk,dpk,得到:dfdt=∂f∂t+∑k(∂f∂qk∂H∂pk−∂f∂pk∂H∂qk),用一些带有启发性的标记标记括号中的内容:dfdt=∂f∂t+∑k∂(H,f)∂(pk,qk),引入泊松括号:{H,f}=∑k∂(H,f)∂(pk,qk)那么dfdt=∂f∂t+∑k∂(H,f)∂(p,q)=∂f∂t+∑k∂(H,f)∂(p,q)=∂f∂t+{H,f}p,q=[∂∂t+{H,⋅}]f那么抽调出“算符”的形式,有:ddt=∂∂t+{H,⋅}回想牛顿力学中,转动参照系的导数,即考虑一个矢量在原参考系的导数和转动参考系的导数,有:dG→dt|fixed=dG→dt|rot+ω→×G→;,抽调成算符:ddt|fixed=ddt|rot+ω→×我们自然的联想到,泊松括号这个二元算符(指有两个量参加运算,如:加法、减法、乘法等),和叉积的相似性。
哈密顿正则方程
/ 2m
m
x
2
由正则方程,得 x= H p p / m, p - H x m
2
x
由 上 两 式 削 去 p, 得 x 积分,得
2
x 0
x A c o s ( t ) p - m A s in ( t )
表明谐振子的相轨道为沿顺时针方向的 封闭椭圆。 物理与光电工程学院
第五章 分析力学
§5.5 哈密顿正则方程
知识回顾 • 基本形式的拉氏方程
d T d t q T q Q
• 保守系的拉氏方程
d L d t q L q 0
• 拉氏方程的应用
a. 确定自由度 b. 选取广义坐标 c. 写出体系的拉氏函数 d. 解拉氏方程并讨论
2 2
由 定 义 求 p , 并 进 而 求 哈 密 顿 函 数 H : L , p L m r 2 sin 2 pr m r, p m r r p 1 m ( r 2 r 2 2 r 2 sin 2 2 ) H L p r r p 2 r 1 2 m ( pr
g
二、正则方程
y
v,
g u
x
(正则形式)
前面提过,用 P代换L函数的 q有一定的优越性, q P 但只用 代换 而不改变函数的形式,则原函数 对新变量无正则形式,给计算带来麻烦.
下面把L函数: L ( q , q , t ) 和 f ( x , y ) 比较
a o
5 g sin ( R r ) 7
19
物理与光电工程学院
第七章哈密顿正则方程
H p q j dt 0 j t0 j1 q j
t1 k
对于完整系统,由于δqj 是相互独立的,且可取任何值, 则 H
j p
j
即得关于变量
q , p , t
j
q j
的Hamilton正则方程
t1
k t1 k k j H Qj q j dt L Qj q j dt p j q t0 t0 j j j 1
H H j p j p j q j q qj p j Qj q j dt t0 q j p j j 1
H j p Q j q j
j
1,2, ,k
其中Qj 为系统的非有势力对应于广义坐标 qj 的广义力。
例7-1 试用Hamilton正则方程求出水平弹簧质量振动 系统的运动微分方程 解:单自由度系统, x为广义坐标
L T V
1 2 1 2 1 2 kx L mx V kx 2 2 2 px L x mx 构造H函数 p x m x 1 2 1 2 L px x mx kx H Px x 2 2 px 2 1 2 kx H x, px 2m 2
t1 t1
对上式进行变分运算,得
H H p q q p p q dt 0 j j j j j j t0 p q j 1 j j
t1 k
将上式中的第一项改写成
d j p j q j p j q dt j 1 j 1
j H p j q j H q j p
正则方程
+ g (h f f h ) (h f f h) f p q p q p q p p q p q q
+ h (f g f g ) (f g f g ) f p q p q p q p p q p q q
=........
0
重要性质: 如果f, g是运动积分,则它们的泊松括号也是运动积分。 (注:表达成函数关系) 泊松定律
42. 泊松括号
已知函数
有
代入哈密顿方程 其中已定义
上面括号称为 泊松括号。 运动积分条件
不显含时间时,即要求
算符,对易,。。。
任意两个函数之间的泊松括号 泊松括号的重要性质
特殊情况,如果f, g之一是广义坐标或广义动量,则
雅克比恒等式
证明:
1:直接代入:麻烦计算可得
2:方便技巧方法
左边,对f ,第一项只包含f的一阶微分,第二、三项包含 f的二阶微分,现在来看第二、三项对f的二阶贡献
设 则
D1,D2的一般形式(不包含2次微分形式) 其中系数任意。由此
二、三项对f的2阶微分贡献为0 左边只有二阶微分贡献?
f {g, h} {g, h} f
p q
p q
g {h, f } {h, f } g
p q
p q
h { f , g} { f , g} h
p q
p q
= f (g h h g ) (g h h g ) f p q p q p q p p q p q q
① 分析体系约束类型,主动力性质;
② 确定自由度,选择适当的广义坐标;
③ 正确写出体系的L函数和H函数;
④ 将L或H代入相应的哈密顿理论的动力学方程,并进行运算,可
得出体系的运动微分方程;
哈密顿正则方程
§6.哈密顿正则方程引言:哈密顿正则方程是与拉氏方程:0=∂∂-∂∂a a q L qL dt d 等价的动力学方程。
s q L q L dt d a a ,....2,10==∂∂-∂α ,这组拉氏方程是s 个关于广义坐标a q 的二阶常微分方程。
在这组拉氏方程中的拉氏函数L 它是广义坐标q ,广义速度q以及时间t 的函数:),,(t qq L L =。
如果我们把拉氏函数中的广义速度a q 变换成→广义动量αp ,即),,(t p q L L =那么就可以将上面的s 个拉氏方程①化成2s 个一阶常微分方程,而且这2s 个一阶常微分方程还具有一定的很漂亮的对称性②具有一定的对称性。
要想把拉氏函数:),,(t qq L L =变成是广义坐标、广义动量P 及时间t 的函数→),,(t p g L L =,以及将s 个拉氏方程化成2s 个一阶常微分方程。
将会用到勒襄特变换这一数学工具。
∴得先介绍一下:一.勒襄特变换(只作了解,不作要求,大纲不要求讲这部分内容)现在先讨论两个变量的勒襄德变换,假设所给的函数是两个变量x 1 和x 2的函数,即:),(21x x f f =。
则由高等数学的知识可得此函数的全微分:2211dx x f dx x f df ∂∂+∂∂=在此我们令11x f u ∂∂=,22x f u ∂∂=,[ii x f u ∂∂=(i=1,2)]……①并以1u 和2u 为新的变量定义一个新函数g: ∑=-+=-≡212211i i if u x u x f u xg ……②如果我们从变换方程①解出i x ,使i x 是i u 的函数,即)(i i i u x x =,再代入上式②中去,那么,g 就是只含新变量i u 的函数了,即:),(21u u g g =。
我们先对②式两边进行微分,则得:∑∑∑∑=====⎥⎦⎤⎢⎣⎡⋅∂∂-+=∂∂-+=21212121)()(i i i i i i i i i i i i i i i i i du x dx x f u du x dx x f dx u du x dg 又∵将旧变量i x 换成新变量i u 之后,新函数g 就是新变量i u 的函数:),(21u u g g =那么对它微分就有:2211du u g du u g dg ∂∂+∂∂=……*′,将这个等式与上一等式进行比较就可得到变换关系:11u g x ∂∂=,222du u g x ∂∂=……③前面我们利用变换方程①把旧的变量x 1,x 2及旧的函数),(21x x f 变为新的变量21,u u 及新的函数),(21u u g g =的方法,就称为勒让德变换。
哈密顿正则方程例题
典型例题3 典型例题
经验告诉我们,用一枝铅笔的笔尖与水平桌面接触, 经验告诉我们,用一枝铅笔的笔尖与水平桌面接触,使之 竖直地稳定转动是很困难的,一长为 竖直地稳定转动是很困难的,一长为10cm、直径为 、直径为0.8cm的 的 铅笔,即使以角速度ω0=100rot/s高速转动,也不能稳定地 铅笔,即使以角速度 高速转动, 高速转动 竖直转动,试用分析力学方法解释 竖直转动, 分析: 分析: 铅笔是否能稳定地竖直转动
从拉格朗日函数的表达式知: ψ,ϕ为循环坐标,故: 从拉格朗日函数的表达式知: , 为循环坐标,
& & & pϕ = ∂L = J*ϕ sin 2 θ + J z cos θ (ϕ cos θ + ψ ) = C1 & ∂ϕ & & pψ = ∂L = J z (ϕ cos θ +ψ ) = C2 & ∂ψ
J z 2 ω0 2 (1 − cos θ) (1 − cos θ) cos θ J *&& = θ [ − 1] + mgrOC sin θ 2 J* sin θ sin θ & d J z 2 ω0 2 (1 − cos θ) 2 dθ & [ ] + mgrOC cos θ} = J *θ =− { d θ 2 J* sin θ dθ 1 & 2 J z 2 ω0 2 (1 − cos θ) 2 + mgrOC cos θ = E 积分有: 积分有: J *θ + 2 2 2 J * sin θ
哈密顿正则方程应用
mm & r 1 2 r & r=f m +m 1 2
s 1 µr 2 & T= 2
'
'
r r µ&& = f r
相对质心运动 相对质心 运动动能
θ
r
µ r
1 m x2 + y2 + z2) T= (& & & c 2 2 2 1
1 µ(r2 +r2θ2 +r2ϕ2 sin2 θ) & & & T= 2
& H = −L+θpθ
3 m R−r)2θ2 −mg(R−r)cosθ +θp & & =− ( θ 4
p θ H= −m (R−r)cosθ g 2 3m R−r) (
2
p θ H= −m (R−r)cosθ g 2 3m R−r) (
则 方 程
2
& = ∂H = 2pθ 根 θ 据 ∂pθ 3m(R−r)2 正 &θ = −∂H = −mg(R−r)sin θ p ∂θ & 2pθ 2g 2gsinθ &= & θ
的
正 则 方 程 解 题 步 骤 分析约束, 分析约束,确定自由度 选好广义坐标 写出系统的T,V,L 写出系统的 写出 H = −L+qα ∂L (α =12,3......) & , ∂qα H=H(q,p,t)
Attention:
广 义
& = ∂L = pα (q, q,t) α & ∂qα
& & qα = qα (q, p,t)
∂L = µr & pr = & ∂r ∂L = µr2θ & p = & θ ∂θ ∂L = µr2ϕsin 2 θ & p = ϕ & ∂ϕ
哈密顿正则方程
r1'
r
r2'
m2 rm1' 1
r2'
r2'
m1
m1 m2
r
2' 1'
m1
m1 m2
m2
m1 m2
S
r S'
m1 r1' c
r2m' 2
'2 2
m1 m1 m2
2
2
'2 1
m2 m1 m2
2
2
T 1 m1m2 2 1 2
2 m1 m2
z p z
m
H L p p p z z
H
1 2m
p
2
p2
2
p
2 z
V
(
,
,
z)
(3)在球面坐标系中
第8页/共17页
T 1 m(r2 r 2 2 r 2 2 sin 2 ) ,V=V(r,,)
2
L 1 m(r2 r 22 r 2 2 sin2 ) V(r,,)
2
pr
L r
mr,
p
L
mr 2, p
L
mr 2 sin 2
r pr , p ,
p
m
mr 2
mr 2 sin 2
H
1 2m
p
2 r
p2 r2
p2
r 2 sin 2
V(r,,)
第9页/共17页
[例10] 求弹性双原子分子的拉格朗日函数和哈密顿函数。设两 原子之间相互作用的弹性力为 F = -k(r-r0) 其中r为两原子间距离,r0为两原子处在平衡时的距离。
解: 为了求出拉格朗日函数,应先求分子的动能。 从寇尼格定理可知,分子动能
3-2哈密顿函数和正则方程
y
ω
x
mg
o x
V mgy
y
x
2
y
xx 2a
4a
L T V
1 2
2 2 2 2 2 m x ( 1 x / 4 a ) x mgx
2
/ 4a
H T2 To V
正则方程
p H 2 p ml H p mgl sin
g sin 0 消去 p 得: l
当 5 时 , sin
为简谐振动
g 0 l
例5:如图所示,匀质绳质量为 m ,长为 l ,一半放在光滑桌面上, 另一半垂挂于桌外,绳无初速下落,问绳全部离开桌面的瞬间,速 度多大? (用正则运动方程求解) 解:建立如图坐标系,取x为广义坐标,以桌面为重力势能零点。
j1
j
s
j
d q j q j dp j) dL
L
其中
L dL dq j1 q j
( p
j1 s j
L dq j dt q j t
L t dt
dq
j
p j d q j)
dH ( p j d q j q j dp j) dL , dL
L T V
由 px
动能 T
1 2
mx
2
mgx 2l px
m
2
yห้องสมุดไป่ตู้
x x
L x
1 2
m x 得: x
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1 2 力学系的动能为( 力学系的动能为(圆柱中心轴的转动惯量 J D = mr ): 2 1 1 1 & 1 1 1 & & & T = m ' y 2 + ms 2 + J D θ2 = ( m '+ m + m) y 2 2 2 2 2 2 4
2011年9月15日 年 月 日
理论物理导论
& & & 2 g (m '− m sin α) dy dy ydy = = 求得: y 求得: && = dt dy dy 2m '+ 3m
分离变量积分: 分离变量积分:
∫
& y
0
& & ydy = ∫
y+h y
2 g (m '− m sin α) dy 2m '+ 3m
2011年9月15日 年 月 日
理论物理导论
若把铅笔视为圆柱,则有: 若把铅笔视为圆柱,则有:
J z = 1 mr 2 = 0.08m, J x = J y = J* = 1 ml 2 = 100 m 2 3 3
铅笔定点转动时的拉格朗日函数为: 铅笔定点转动时的拉格朗日函数为:
L = T − V = 1 J* (ωx2 + ω y 2 ) + 1 J zωz 2 − V 2 2 & & & & = 1 J* (θ 2 + ϕ 2 sin 2 θ ) + 1 J z (ϕ cos θ + ψ )2 − mgrOC cos θ 2 2
理论物理导论 对拉格朗日方程化简有: 对拉格朗日方程化简有:
J z 2 ω0 2 (1 − cos θ) (1 − cos θ) cos θ J *&& = θ [ − 1] + mgrOC sin θ 2 J* sin θ sin θ & d J z 2 ω0 2 (1 − cos θ) 2 dθ & [ ] + mgrOC cos θ} = J *θ =− { d θ 2 J* sin θ dθ 1 & 2 J z 2 ω0 2 (1 − cos θ) 2 + mgrOC cos θ = E 积分有: 积分有: J *θ + 2 2 2 J * sin θ
广义坐标, 对θ广义坐标,拉格朗日方程为: 广义坐标 拉格朗日方程为: d ∂L ∂L & & & & ( &)− = J *&& − J *ϕ2 sin θ cos θ + J z ϕ sin θ(ϕ cos θ + ψ) θ dt ∂θ ∂θ −mgrOC sin θ = 0
2011年9月15日 年 月 日
& pθ = −
∂H ( m '+ m) sin θ cos θ 2 =− 2 pθ − mgl sin θ 2 2 ∂θ l ( m '+ m sin θ)
2
( ml cos θ) 2 & & & & θ + ml 2 θ 又 pθ = mlx cos θ + ml θ = − m '+ m 求导并代入上式化简可得关于θ的微分方程 的微分方程: 求导并代入上式化简可得关于 的微分方程:
注意到动能是广义速度的二次齐次函数, 注意到动能是广义速度的二次齐次函数,故:
p2 H = T +V = − (m '− m sin α) gy − mgl sin α 2m '+ 3m
代入哈密顿正则方程得: 代入哈密顿正则方程得:
∂H 2p 2m '+ 3m &= && = ⇒p y ∂p 2m '+ 3m 2 ∂H 2m '+ 3m & =− && p = ( m '− m sin α) g = y 2 ∂y & y=
2011年9月15日 年 月 日
理论物理导论
r r 代入初始条件 t = 0, θ = 0, ω = ω0 k ,ψ = ω0 k , ϕ = 0 & & r
求得 C1 = C2 = J zω0
& & & pϕ = ∂L = J*ϕ sin 2 θ + J z cos θ (ϕ cos θ + ψ ) = J zω0 即 & ∂ϕ & & pψ = ∂L = J z (ϕ cos θ +ψ ) = J zω0 & ∂ψ & cos θ + ψ , ϕ = J zω0 (1 − cos θ ) & ω & 化简有: 化简有: z = ω0 = ϕ J* sin 2 θ
α
ω
m
m'
r
s
x
哈密顿函数。 哈密顿函数。圆柱体作平动
y
和定轴转动,需两个坐标; 和定轴转动,需两个坐标;重物作一维运动需一个坐标 力系受轻绳约束,故只需一个坐标,不妨取重物坐标 力系受轻绳约束,故只需一个坐标,不妨取重物坐标y
2011年9月15日 年 月 日
理论物理导论
解: 设绳长为l,柱体轴心离坡顶距离为s,柱体滚动的角度 设绳长为 ,柱体轴心离坡顶距离为 柱体滚动的角度 大小为θ,下垂的重物到坡顶的距离为 ,约束关系为: 大小为 ,下垂的重物到坡顶的距离为y,约束关系为:
ζ
分析: 分析:
z
铅笔是否能稳定地竖直转动 关键是看在竖直方向(θ=0)时, 时 关键是看在竖直方向 铅笔的势能是否处于极小值。 铅笔的势能是否处于极小值。
ξ
2011年9月15日 年 月 日
C
y
θ
r mg
η ϕ
x
O
理论物理导论 解: 象对称陀螺一样取坐标系。那么 ,ψ,ϕ为广义坐标 象对称陀螺一样取坐标系。那么θ, , 显然, 显然,坐标轴皆为惯量主轴 初始时刻, 初始时刻,铅笔直立
理论物理导论
典型例题1 典型例题
一质量为m、半径为r的圆柱体置于坡角为α的斜面上;柱轴 一质量为 半径为 的圆柱体置于坡角为 的斜面上; 用轻绳过坡顶的光滑滑轮与一质量为m’的重物相连结, 用轻绳过坡顶的光滑滑轮与一质量为 的重物相连结,如图 的重物相连结 重物自静止开始下落,求重物的加速度和下落h距离后的速度 重物自静止开始下落,求重物的加速度和下落 距离后的速度 分析:利用哈密顿正则方程求解 分析: 需先求得力学系的动能和势能 从而求得拉氏函数, 从而求得拉氏函数,最终获得
注意到动能是广义速度的二次齐次函数, 注意到动能是广义速度的二次齐次函数,故:
H ( x, θ, px , pθ ) = T + V ml 2 px 2 − 2ml cos θpx pθ + (m '+ m) pθ 2 = − mgl cos θ 2 2 2ml (m '+ m sin θ)
2011年9月15日 年 月 日
& (m '+ m sin 2 θ)l&& + ml θ2 sin θ cos θ + (m '+ m) g sin θ = 0 θ
11年9月15日 年 月 日
理论物理导论
典型例题3 典型例题
经验告诉我们,用一枝铅笔的笔尖与水平桌面接触, 经验告诉我们,用一枝铅笔的笔尖与水平桌面接触,使之 竖直地稳定转动是很困难的,一长为 竖直地稳定转动是很困难的,一长为10cm、直径为 、直径为0.8cm的 的 铅笔,即使以角速度ω0=100rot/s高速转动,也不能稳定地 铅笔,即使以角速度 高速转动, 高速转动 竖直转动,试用分析力学方法解释 竖直转动,
理论物理导论 显然,哈密顿函数中不显含 ,故为循环坐标: 显然,哈密顿函数中不显含x,故为循环坐标:
& & px = ( m '+ m) x + ml θ cos θ = C |t =0, x =0,θ=0 = 0(x方向动量守恒) 方向动量守恒) 方向动量守恒 & ml θ cos θ & 即 x=− m '+ m 由关于θ的正则方程 的正则方程, 可得: 由关于 的正则方程,并注意到 px = 0 可得:
2011年9月15日 年 月 日
ζ
z
r r t = 0, θ = 0, ω = ω0 eζ = ω0 k r
任意时刻, 任意时刻,有:
r r r eζ = sin θ j + cos θ k
任意时刻,铅笔角速度为: 任意时刻,铅笔角速度为:
C
y
θ
r mg
η
ϕ
x
O
ξ r r r & r & eζ + ψ k & ω = θi +ϕ r r r & & & & = θ i + ϕ sin θ j + (ϕ cos θ + ψ ) k r r r = ωx i + ω y j + ω z k
力系的势能为: 力系的势能为: V = − mgl cos θ
2011年9月15日 年 月 日
理论物理导论 相应的广义动量为: 相应的广义动量为:
∂L ∂T & & px = = = (m '+ m) x + ml θ cos θ & & ∂x ∂x ∂L ∂T & & pθ = = = mlx cos θ + ml 2 θ & & ∂θ ∂θ 由以上两式可解得: 由以上两式可解得: pθ px − cos θ l & = (m '+ m) pθ − mlpx cos θ θ & x= 2 ml 2 ( m '+ m sin 2 θ) m '+ m sin θ