Processing math: 100%

两周速通朗道《力学》(3)

(3/14)

Posted by Whitney on November 25, 2024

早上十点慵懒的起床,状态不错。今天下雨,体育课暂停一次,有一整个下午和晚上来学习。

运动方程的积分

一维运动

  如果取笛卡尔坐标,那么只有一个自由度的系统的运动的拉格朗日函数就可以写为L=1/2m˙xU(x)。由于对于一个封闭系统,能量是守恒的,那么可以直接写出其能量表达式E=1/2m˙x+U(x)。这是一个一阶微分方程,可以直接分离变量再积分得到运动方程的解:

dxdt=2m[EU(x)]t=m2dxEU(x)+const

由于EU(x)需要恒正,那么实际的运动只能发生在U(x)<E的区域内,而U(x)=E给出了这个运动的边界,这个方程给出的点速度为零,被称作转折点。由两个转折点给出的运动区域被称作有界运动,由一个或者没有转折点的运动被称作无界运动。一维有界运动的运动区域被称作势阱,其运动被称作振动。在牛顿力学的体系下,时间是可逆的,因此两个转折点x1,x2之间往返的时间应该是一致的。定义x1x2x1的运动所花的时间为周期T,那么

T=2mx2x1dxEU(x)

注意这里积分的上下限(即两个转折点)仍需给定力学系统的总能量才可以得到。

例子:平面单摆

  平面单摆的总能量为

E=mglcosφ0=ml2˙φ22mglcosφ

其中φ0是最大摆角。由于摆动过程的对称性,周期可以被写作0φ0过程所需时间的四倍,因此

dφdt=2ml2(E+mglcosφ)t=ml22φ00dφE+mglcosφT=4t=4l2gφ00dφcosφ0+cosφ=2lgφ00dφsin2(φ0/2)sin2(φ/2).

  这个积分怎么算?考虑第一类完全椭圆积分:

K(k)=π/20dθ1k2sin2θ=π2n=0[(2n)!22nn!2]2k2n

那么我们只需在周期积分中令sin(φ/2)/sin(φ0/2)=sinξ,那么就会有

dξ=12cos(φ/2)sin(φ0/2)11sin2(φ/2)sin2(φ0/2)dφdφ=2sin(φ0/2)cos(φ/2)1sin2(φ/2)sin2(φ0/2)dξT=2lgφ00dφsin2(φ0/2)sin2(φ/2)=4lgπ/20sin(φ0/2)cos(φ/2)1sin2(φ/2)sin2(φ0/2)1sin2(φ0/2)sin2(φ/2)dξ=4lgπ/201cos(φ/2)dξ=4lgπ/2011sin2(φ/2)dξ=4lgπ/2011sin2(φ0/2)sin2ξdξ=4lgK(sinφ02)

当最大摆角很小时,sinφ0/2φ0/2 按照幂级数展开可得

T=2πlg(1+116φ20+)

舍去高阶小量,余下的一项便是我们熟知的单摆小振动的周期表达式。

根据周期确定势能

  我们能够通过势能确定周期。现在我们试图已知周期求势能。这其实就是求解周期的积分方程。对于有界的运动,x=x(U)是双值的,因此我们将积分化为两个部分:

T(E)=2mE0(dx2dUdx1dU)dUEUα0T(E)αEdE=2mα0E0(dx2dUdx1dU)dUdE(EU)(αE)α0T(E)αEdE=2mα0(dx2dUdx1dU)dUE0dE(EU)(αE)α0T(E)αEdE=π2m[x2(α)x1(α)]x2(U)x1(U)=1π2mα0T(E)αEdE

也就是说给定周期,只能确定到两个坐标的差值,然而坐标本身无法确定,这也反映了势能可以任意加减一个常数的事实。但是如果势能曲线关于U轴对称,那么U(x)就可以唯一确定。

约化质量

  仅由两个质点组成的力学系统非常重要,我们在此单独研究。两个质点间的势能仅依赖于他们之间的距离,因此拉格朗日函数可以写为

L=12m1˙r1+12m2˙r2U(|r1r2|)

我们可以引入位矢差r=r1r2,并取质心系,那么在质心系中就有m1r1+m2r2=0,反解出

r1=m2m1+m2r,r2=m1m1+m2r

代入拉格朗日函数得

L=12m1m2m1+m2˙r2U(r)

因此一个二体问题可以简化为一个质点在外场U(r)中的运动,其中质点的质量为

m=m1m2m1+m2

我们称其为约化质量

有心力场中的运动

朗道对于一些符号的使用有些奇怪,在此使用更熟悉、更常见的符号。

  由角动量的定义J=r×p,我们可以发现,如果一个力学系统的角动量不变,那么径矢总是保持在一个平面内的,因此轨道都在一个平面内。在那个平面内研究力学系统,引入极坐标r,φ,拉氏量就可以写作

L=12m(˙r2+r2˙φ2)U(r)

由于上式不显含坐标φ(这一点也说明了有心力场下空间的旋转对称性,或者说角动量守恒的系统空间旋转对称),那么我们就可以对坐标φ写运动方程

ddtL˙φ=Lφ=0

这也表明广义动量L/˙φ也是一个守恒量,p=mr2˙φ,这即是角动量守恒。

  接下来我们来推导开普勒第二定律。考虑转过一小角形成的扇形面积,ds=1/2r2dφ,那么角动量守恒就可以写作2m˙s=const,这正是开普勒第二定律,相等时间间隔内质点径矢扫过相同的面积,掠面速度不变。

  为了将能量表达式换成r˙r的函数,将J=mr2˙φ代入能量表达式,我们有

E=12m˙r2+J2mr2+U(r)

用相似的方法分离变量积分能得到两个坐标关于t的关系

t=dr2m[EU(r)]J2m2r2φ=(J/r2)dr2m[EU(r)]J2/r2

能量的表达式又可以看作是 E=T+Ueff,是一个一维运动,有效势能是 Ueff=U(r)+J2/2mr2,我们称后一项J2/2mr2离心势能Ueff=E给出的转折点是运动过程的最小和最大径矢。然而在这个情况下,转折点处仍有角向的速度˙φ。如果这个方程给出两个转折点,那么轨道就全部在两个转折点确定的环形区域内,然而不一定是封闭曲线:只有两个转折点之间角向的积分

Δφ=rmaxrmin(J/r2)dr2m[EU(r)]J2/r2

2π的整数倍时,轨道才能够封闭。轨道是否封闭完全取决于势能的形式;只有当势能与1/rr2成正比时,其中一切有界运动的轨道都是封闭的。

  当r0时,离心势能应当足够快的趋近于质点才可能坠入场的中心。离心势能像J2/2mr2一样趋近于+,因此U(r)必须像α/r2(α>J2/2m)或者像k/rn这样趋向于

开普勒问题

  我们来研究引力场下的运动。引力场的势能U(r)=α/rUeff=α/r+J2/2mr2。这个势能在r0时趋近于+,在r时趋近于0,最小值为U(J2/mα)=mα2/2J2。当E>0时运动无界,E<0时运动有界。

  通过φr的积分可以得到极坐标系下的轨道方程

φ=arccosJ/rmα/J2mE+m2α2/J2

p=J2/mαe=1+2EJ2/mα2,轨道方程可以重写为

r=p1+ecosφ

这正是我们熟悉的圆锥曲线方程。当E<0时,运动有界,对应e<1,此时轨道是椭圆。当能量取最小值时,e=0,椭圆变成圆。椭圆轨道的周期仅与总能量有关。用面积积分公式能得到

T=παm2E3.

E0时,运动无界,此时轨道可能取双曲线或是抛物线,对应的偏心率e1。当一个静止的质点从无穷远跌向场心,那么对应的E=0,轨迹变为抛物线。

拉普拉斯-龙格-楞次矢量

  在形如U(r)=α/r的有心力场中,有一个特别的守恒量,称作拉普拉斯-龙格-楞次矢量,简称LRL矢量:

p×J+αrr=const.

下午花了四个小时写了这么多,也算是速通了。现在是北京时间 11.25 16:03,济南小雨。可能回宿舍还需要一些时间把习题写一下?今日p24-41。