早上十点慵懒的起床,状态不错。今天下雨,体育课暂停一次,有一整个下午和晚上来学习。
运动方程的积分
一维运动
如果取笛卡尔坐标,那么只有一个自由度的系统的运动的拉格朗日函数就可以写为L=1/2m˙x−U(x)。由于对于一个封闭系统,能量是守恒的,那么可以直接写出其能量表达式E=1/2m˙x+U(x)。这是一个一阶微分方程,可以直接分离变量再积分得到运动方程的解:
dxdt=√2m[E−U(x)]t=√m2∫dx√E−U(x)+const由于E−U(x)需要恒正,那么实际的运动只能发生在U(x)<E的区域内,而U(x)=E给出了这个运动的边界,这个方程给出的点速度为零,被称作转折点。由两个转折点给出的运动区域被称作有界运动,由一个或者没有转折点的运动被称作无界运动。一维有界运动的运动区域被称作势阱,其运动被称作振动。在牛顿力学的体系下,时间是可逆的,因此两个转折点x1,x2之间往返的时间应该是一致的。定义x1→x2→x1的运动所花的时间为周期T,那么
T=√2m∫x2x1dx√E−U(x)注意这里积分的上下限(即两个转折点)仍需给定力学系统的总能量才可以得到。
例子:平面单摆
平面单摆的总能量为
E=−mglcosφ0=ml2˙φ22−mglcosφ其中φ0是最大摆角。由于摆动过程的对称性,周期可以被写作0→φ0过程所需时间的四倍,因此
dφdt=√2ml2(E+mglcosφ)t=√ml22∫φ00dφ√E+mglcosφT=4t=4√l2g∫φ00dφ√−cosφ0+cosφ=2√lg∫φ00dφ√sin2(φ0/2)−sin2(φ/2).这个积分怎么算?考虑第一类完全椭圆积分:
K(k)=∫π/20dθ√1−k2sin2θ=π2∞∑n=0[(2n)!22nn!2]2k2n那么我们只需在周期积分中令sin(φ/2)/sin(φ0/2)=sinξ,那么就会有
dξ=12cos(φ/2)sin(φ0/2)1√1−sin2(φ/2)sin2(φ0/2)dφdφ=2sin(φ0/2)cos(φ/2)√1−sin2(φ/2)sin2(φ0/2)dξT=2√lg∫φ00dφ√sin2(φ0/2)−sin2(φ/2)=4√lg∫π/20sin(φ0/2)cos(φ/2)√1−sin2(φ/2)sin2(φ0/2)1√sin2(φ0/2)−sin2(φ/2)dξ=4√lg∫π/201cos(φ/2)dξ=4√lg∫π/201√1−sin2(φ/2)dξ=4√lg∫π/201√1−sin2(φ0/2)sin2ξdξ=4√lgK(sinφ02)当最大摆角很小时,sinφ0/2≈φ0/2 按照幂级数展开可得
T=2π√lg(1+116φ20+⋯)舍去高阶小量,余下的一项便是我们熟知的单摆小振动的周期表达式。
根据周期确定势能
我们能够通过势能确定周期。现在我们试图已知周期求势能。这其实就是求解周期的积分方程。对于有界的运动,x=x(U)是双值的,因此我们将积分化为两个部分:
T(E)=√2m∫E0(dx2dU−dx1dU)dU√E−U∫α0T(E)√α−EdE=√2m∫α0∫E0(dx2dU−dx1dU)dUdE√(E−U)(α−E)∫α0T(E)√α−EdE=√2m∫α0(dx2dU−dx1dU)dU∫E0dE√(E−U)(α−E)∫α0T(E)√α−EdE=π√2m[x2(α)−x1(α)]x2(U)−x1(U)=1π√2m∫α0T(E)√α−EdE也就是说给定周期,只能确定到两个坐标的差值,然而坐标本身无法确定,这也反映了势能可以任意加减一个常数的事实。但是如果势能曲线关于U轴对称,那么U(x)就可以唯一确定。
约化质量
仅由两个质点组成的力学系统非常重要,我们在此单独研究。两个质点间的势能仅依赖于他们之间的距离,因此拉格朗日函数可以写为
L=12m1˙r1+12m2˙r2−U(|r1−r2|)我们可以引入位矢差r=r1−r2,并取质心系,那么在质心系中就有m1r1+m2r2=0,反解出
r1=m2m1+m2r,r2=−m1m1+m2r代入拉格朗日函数得
L=12m1m2m1+m2˙r2−U(r)因此一个二体问题可以简化为一个质点在外场U(r)中的运动,其中质点的质量为
m=m1m2m1+m2我们称其为约化质量。
有心力场中的运动
朗道对于一些符号的使用有些奇怪,在此使用更熟悉、更常见的符号。
由角动量的定义J=r×p,我们可以发现,如果一个力学系统的角动量不变,那么径矢总是保持在一个平面内的,因此轨道都在一个平面内。在那个平面内研究力学系统,引入极坐标r,φ,拉氏量就可以写作
L=12m(˙r2+r2˙φ2)−U(r)由于上式不显含坐标φ(这一点也说明了有心力场下空间的旋转对称性,或者说角动量守恒的系统空间旋转对称),那么我们就可以对坐标φ写运动方程
ddt∂L∂˙φ=∂L∂φ=0这也表明广义动量∂L/∂˙φ也是一个守恒量,p=mr2˙φ,这即是角动量守恒。
接下来我们来推导开普勒第二定律。考虑转过一小角形成的扇形面积,ds=1/2r2dφ,那么角动量守恒就可以写作2m˙s=const,这正是开普勒第二定律,相等时间间隔内质点径矢扫过相同的面积,掠面速度不变。
为了将能量表达式换成r,˙r的函数,将J=mr2˙φ代入能量表达式,我们有
E=12m˙r2+J2mr2+U(r)用相似的方法分离变量积分能得到两个坐标关于t的关系
t=∫dr√2m[E−U(r)]−J2m2r2φ=∫(J/r2)dr√2m[E−U(r)]−J2/r2能量的表达式又可以看作是 E=T+Ueff,是一个一维运动,有效势能是 Ueff=U(r)+J2/2mr2,我们称后一项J2/2mr2为离心势能。Ueff=E给出的转折点是运动过程的最小和最大径矢。然而在这个情况下,转折点处仍有角向的速度˙φ。如果这个方程给出两个转折点,那么轨道就全部在两个转折点确定的环形区域内,然而不一定是封闭曲线:只有两个转折点之间角向的积分
Δφ=∫rmaxrmin(J/r2)dr√2m[E−U(r)]−J2/r2是2π的整数倍时,轨道才能够封闭。轨道是否封闭完全取决于势能的形式;只有当势能与1/r或r2成正比时,其中一切有界运动的轨道都是封闭的。
当r→0时,离心势能应当足够快的趋近于−∞质点才可能坠入场的中心。离心势能像J2/2mr2一样趋近于+∞,因此U(r)必须像−α/r2(α>J2/2m)或者像−k/rn这样趋向于−∞。
开普勒问题
我们来研究引力场下的运动。引力场的势能U(r)=−α/r,Ueff=−α/r+J2/2mr2。这个势能在r→0时趋近于+∞,在r→∞时趋近于0−,最小值为U(J2/mα)=−mα2/2J2。当E>0时运动无界,E<0时运动有界。
通过φ−r的积分可以得到极坐标系下的轨道方程
φ=arccosJ/r−mα/J√2mE+m2α2/J2设p=J2/mα,e=√1+2EJ2/mα2,轨道方程可以重写为
r=p1+ecosφ这正是我们熟悉的圆锥曲线方程。当E<0时,运动有界,对应e<1,此时轨道是椭圆。当能量取最小值时,e=0,椭圆变成圆。椭圆轨道的周期仅与总能量有关。用面积积分公式能得到
T=πα√m2E3.当E≥0时,运动无界,此时轨道可能取双曲线或是抛物线,对应的偏心率e≥1。当一个静止的质点从无穷远跌向场心,那么对应的E=0,轨迹变为抛物线。
拉普拉斯-龙格-楞次矢量
在形如U(r)=α/r的有心力场中,有一个特别的守恒量,称作拉普拉斯-龙格-楞次矢量,简称LRL矢量:
→p×→J+α→rr=const.下午花了四个小时写了这么多,也算是速通了。现在是北京时间 11.25 16:03,济南小雨。可能回宿舍还需要一些时间把习题写一下?今日p24-41。