分析力学给出了力学系统在完全一般性的广义坐标下具有不变性是的动力学方程组,并突出了能量函数的意义。分析力学甚至概括了比牛顿力学广泛得多的系统,例如在电气系统、控制系统中的应用就是一个明显的例子。分析力学的数学形式有着极好的性质,它不仅提供了解决天体力学及一系列动力学问题的较佳途径,同时给量子力学的发展提供了启示,最适宜于成为引向现代物理的跳板。其最小作用量原理提供了建立相对论量子力学最简练而富有概括性的出发点。——梁昆淼
学习时间 | 大二上 |
周课时 | 3 |
本人成绩 | 83 |
课程教材 | 梁昆淼,理论力学(第四版),高等教育出版社,2009 |
个人建议参考教材 | 同上 |
先修课程 | 微积分 力学 电磁学 |
其他推荐教材:
- 鞠国兴,理论力学学习指导与习题解析(第二版),?,?
- 朗道,力学(第五版),高等教育出版社,2007
- 鞠国兴,朗道《力学》解读,高等教育出版社,2014
本人学习成果评论:本人的理论力学学得简直就是个笑话。脱离了牛顿力学的直观,所谓“虚”“广义”“拉格朗日”“哈密顿”“勒让德”“泊松”等字样令我十分痛苦。在课堂上我自我评价 80% 的时间不知道老师在说什么,所幸课下反复阅读、思考、抄作业答案,终于有了些进步。感谢老师不挂之恩!
我们为什么还要学习经典力学这种有局限性的结论?
任何逻辑上封闭的的、其真实性再一定精确程度上被实验所证实的理论永远不会失去价值,二是作为特殊情况下成立的近似结果被所有后来更为精确的理论所包含。当然,这不包括那些存在内在矛盾的理论,它们仅在理论物理发展的某一个阶段具有价值。——朗道
注意:鞠老师的习题解与课本习题在序号上并不是对应的,还请看具体指出是课本习题序号还是习题解序号。
课外阅读:理论力学与机器学习——Lagrangian Neural Networks, Hamiltonian Neural Networks. 我自个儿还没看,有机会看看。
虚功原理 达朗贝尔原理
虚功原理 甩掉约束力
概念:几何约束、运动约束、(非)完整约束、(非)定常约束、约束力、主动力。
完整约束:可积分的运动约束与几何约束在物理实质上没有区别,合称为完整约束,e.g., 无滑动的滚动。
带有约束的力学问题,按“牛顿方式”来求解时,未知数增多,方程个数相应地增多。
自由度
- 总结:\(n\) 个质点(位矢为 \(\vec r_i\))的力学系统,若存在着 \(m\) 个完整约束和 \(k\) 个非完整约束,那么,质点的直角坐标数 \(N=3n\),广义坐标 \(q_i(t)\) 个数等于 \(s=N-m\),自由度等于 \(N-m-k\)。对于只存在完整约束的系统(\(k=0\)),广义坐标的个数就是自由度。如果存在非完整约束(\(k>0\)),广义坐标的个数大于自由度。运用广义坐标后,不再需要考虑完整约束,但非完整约束仍需考虑,并应将它用相应的广义速度表示。
- 基础课处理重点:完整系统——约束全是完整约束的系统。
虚位移
- 一切想象的符合约束条件的无限小可能位移定义为虚位移(时间被冻结)
理想约束(以下处理的重点)
- 满足系统约束力的虚功为 \(0\)(约束力针对系统所有约束力)。\(\delta W=\displaystyle\sum_{i=1}^n\vec N_i\cdot\delta \vec r_i=0\).
- 举例:光滑性、刚性、不可伸长性。
虚功原理
- 当一个只有理想约束的力学系统处于平衡状态时,作用于该系统所有的主动力的虚功之总和为零。\(\delta W=\displaystyle\sum_{i=1}^n\vec F_i\cdot\delta \vec r_i=0\)
- 虚位移有很大任意性 -> 甩掉虚位移 -> 虚位移不独立 换成广义坐标 -> 广义坐标的虚位移 \(\delta\vec r_i=\displaystyle\sum_{\alpha=1}^s\frac{\partial\vec r_i}{\partial q_\alpha}\delta q_\alpha\) -> 广义坐标下的虚功原理 -> 广义力(广义力对应广义坐标,而非单个主动力) \(Q_\alpha=\displaystyle\sum_{i=1}^n\vec F_i\cdot\frac{\partial\vec r_i}{\partial q_\alpha}=0\) -> 完整约束广义虚位移相互独立 \(Q_\alpha=0\).
- 解题步骤:判断是否满足理想约束,识别主动力,确定自由度和广义坐标。列所有主动力的虚功为 \(0\)(真实坐标微分),用广义坐标表示真实坐标,算出用广义坐标微分表示的真实坐标微分。若广义坐标相互独立,则令每个广义坐标微分的系数为 \(0\),解方程。
约束力的求解——拉格朗日乘子法(虚功原理中约束力被甩掉了,如果要求怎么办?)
- 一开始可以不承认相应的约束条件。换句话说,在选取广义坐标时,可以暂不考虑相应的约束条件。这样选取的广义坐标实际上并不独立,他们之间有相应的约束条件相联系。具体见相关笔记。实操起来计算量巨大,我没有一次能独立算对的(逃
达朗贝尔原理
- 以上研究的是静力学问题,现在转到动力学问题。代入牛顿第二定律,
- 在理想约束作用下,主动力和达朗贝尔力的虚功之总和为零。\(\displaystyle\sum_{i=1}^n(\vec F_i-m_i\ddot{\vec r}_i)\cdot\delta \vec r_i=0\).
拉格朗日动力学
上文把约束力甩掉了,但是我们还想把虚位移甩掉。力学系统在广义坐标描述下所具有的最一般形式的动力学方程——拉格朗日动力学。拉格朗日动力学中,一个力学系统的性质可用单个函数——拉格朗日函数——来概括。
拉格朗日方程 甩掉虚位移
对于静力学问题,虚位移中的位矢不独立,我们引入与广义坐标相关的广义力。在动力学问题中,达朗贝尔原理带来的 \(\ddot{\vec{r}}_i\) 怎么处理?
拉格朗日关系:对位矢 \(\vec r_i=\vec r_i(q_1,q_2,...,q_s,t)=\vec r_i(q,t)\) 求导,得
\[\frac{\partial\dot{\vec r}_i}{\partial\dot q_\beta}=\frac{\partial\vec{r}_i}{\partial q_\beta}, \\ \frac{\partial}{\partial q_\beta}(\frac{d}{dt}\vec{r}_i)=\frac{d}{dt}(\frac{\partial\dot{\vec r}_i}{\partial q_\beta}).\]坐标变换带入达朗贝尔原理,再用拉格朗日关系做变换。对于完整系统,广义虚位移独立,每个广义坐标对应的方程为零,得到拉格朗日方程 \(\displaystyle\frac{d}{dt}\frac{\partial T}{\partial\dot q_\alpha}-\frac{\partial T}{\partial q_\alpha}=Q_\alpha\),其中 \(T=\displaystyle\sum_{i=1}^n \frac{1}{2}m_i\vert\dot{\vec r_i}\vert^2\) 是动能。
把广义坐标、广义速度和 \(t\) 当作独立变量,因为考虑的不是实际运动,而是虚幻的所有可能的运动。
拉格朗日函数
如果主动力全是保守力,存在势能函数 \(V\),与广义力的关系为 \(\displaystyle Q_\alpha=-\frac{\partial V}{\partial q_\alpha}\)。定义拉格朗日函数(不可观测量)\(L=T-V\)。得到拉格朗日方程 \(\displaystyle\frac{d}{dt}\frac{\partial L}{\partial\dot q_\alpha}-\frac{\partial L}{\partial q_\alpha}=0\)(您要是实在学不会,把它放进大脑 RAM 上考场,应该能混个总评及格)。
这种类型的拉格朗日方程还可以通过哈密顿原理得到,它是哈密顿作用量 \(S=\displaystyle\int_{t_1}^{t_2}Ldt\) 取极值时的欧拉 - 拉格朗日方程。详见后文变分法。
广义动量 \(\displaystyle p_\alpha=\frac{\partial L}{\partial\dot q_\alpha}\)
基础课程内主要考虑 \(V\) 与 \(\dot q_\alpha\) 无关的情况,于是 \(p_\alpha=\displaystyle\frac{\partial T}{\partial\dot q_\alpha}\))。这时,类比牛顿力学,拉格朗日方程可读作:广义动量的时间变化率等于广义力。
广义动量包含了通常的动量、角动量等物理量,它的量纲取决于广义坐标的量纲。
例:拉格朗日动力学。
诺特定理
运动积分
在系统运动过程中,存在 \(q_\alpha\) 和 \(\dot q_\alpha\) 的某些函数,它们不随时间而变,这些函数成为系统的运动积分(守恒律概念的推广)。
广义动量积分
\(\displaystyle\frac{\partial L}{\partial q_\beta}=0\)(\(q_\beta\) 是可遗坐标/循环坐标)-> \(\displaystyle\frac{d}{dt}\frac{\partial L}{\partial \dot q_\beta}=0\) -> 广义动量守恒 \(p_\beta=\text{const}\) 成为广义动量积分
拉格朗日系统并不以牛顿第三定律为先决条件,而是根据某广义坐标是否为可遗坐标,也可以处理电磁场问题(其中质点机械运动动量并不守恒,在拉格朗日函数中计入广义势之后,广义动量包括电磁场的动量,广义动量守恒)。
广义能量积分 \(H=\displaystyle\sum_{\alpha=1}^sp_\alpha\dot q_\alpha-L\)
- \(\displaystyle\frac{dH}{dt}=\frac{\partial L}{\partial t}=0\).
- 若 \(\displaystyle\frac{\partial\vec r_i}{\partial t}=0\)(定常约束)-> 动能是广义速度的二次齐次式 -> \(H\) 就是机械能 (齐次函数的欧拉定理参见微积分)
- 关于矢量力学和拉格朗日力学关于机械能守恒条件的区别,见教材 73-74 页。
诺特定理
- 德语小知识:th 中的 h 不发音。Amalie Emmy Noether, German mathematician.
入门级材料阅读 David Morin “Introduction to Classical Music” Chapter 5 5.6 Noether’s Theorem
守恒量 \(I(q,\dot q)=\displaystyle\sum_{k=1}^n\frac{\partial L}{\partial\dot q_k}\frac{d}{d\epsilon}\Phi(q,\epsilon)\vert_{\epsilon=0}\).
教材原话对于一直到大四的我来说都是相当晦涩和不可接近的,不如直接看 tangible 的例子:鞠老师《习题解》的讲义 3.1.2 部分和习题 3.30 3.31(两质点间有势力相互作用,时间平移、空间平移、空间转动拉格朗日函数守恒)。注意诺特定理在变换后与原来不相等的推广形式。其中习题 3.31 就是诺特定理的推广形式,时间平移变换所对应的守恒量守恒量就是机械能。
非完整系统的动力学在 3.3 节讲解,我懒得看,喜欢花样滑冰的同学可以看一看。
拉格朗日力学的推广
例:带电粒子在电磁场中受洛伦兹力 \(\vec F=q(\vec E+\vec v\times\vec B)\)。在洛伦兹规范下,\(\vec E+\displaystyle\frac{\partial\vec A}{\partial t}=-\nabla V\), \(\vec B=\nabla\times\vec A\)。于是 \(\vec F=q[\displaystyle-\nabla V-\vec v\times(\nabla\times\vec A)]\)…… -> \(\displaystyle\frac{d}{dt}(\vec p+q\vec A)=\nabla[q(V-\vec v\cdot\vec A)]\), \(\displaystyle\frac{d}{dt}(T+qV)=\frac{\partial}{\partial t}[q(V-\vec v\cdot\vec A)]\)。详见 Griffiths 10.1.4.
从而定义广义势能 \(U(q,\dot q,t)=e\varphi-e\vec A\cdot\vec v\),拉格朗日函数成为 \(L=T-U\),原来的拉格朗日方程仍成立。此时,研究的守恒量,如广义动量守恒,是不以牛顿第三定律为基础的。
通过我每篇表头,可以知道我们先学理论力学后学电动力学。因此,在理论力学教材看到第二段的时候,我很蒙蔽,想着学了电动力学可能就会了。然而学了电动力学,我的矢量分析还是毫无长进,也不敢回过头来看能不能懂这个推广的拉格朗日函数。当然,在看 Griffiths 不厌其烦的细致讲解之前,我是不能懂的。在此,我代表全体物理学院笨蛋,强烈支持循循善诱的讲解,而非点到为止。
有心力 散射问题
关于有心力作用下运动的性质,见力学 - 有心力。
我们关心两体问题,因为三体问题是邪恶的(见下文哈密顿力学-泊松括号-泊松定理)。如何在拉格朗日动力学的框架下简化两题问题?重点:理解《力学》中讨论过的质心系。
采用直角坐标系研究质心的运动,采用极坐标系研究第二个质点相对于第一个质点的相对运动……得到拉格朗日函数 \(\displaystyle L=\frac{1}{2}m_0(\dot x_0^2+\dot y_0^2+\dot z_0^2)+\frac{1}{2}m^\prime(\dot\rho^2+\rho^2\dot\varphi^2)-V(\rho)\). 其中 \(m_0\) 是质心质量,\(m^\prime\) 是约化质量。
观察拉格朗日函数,显然质心动量能量守恒,甩掉。相对运动具有
- \(\varphi\) 是可遗坐标,质点相对于极点的角动量守恒 \(h=\rho^2\dot\varphi\).
- \(L\) 不显含时间,机械能守恒 \(\displaystyle E^\prime=\frac{1}{2}m^\prime\dot\rho^2+\frac{1}{2}m^\prime\frac{h^2}{\rho^2}+V(\rho)\). 通常称 \(\displaystyle V_{\text{eff}}(\rho)=\frac{1}{2}m^\prime\frac{h^2}{\rho^2}+V(\rho)\) 为有效势能,用于研究相对径向的运动。
与可遗坐标对应的广义动量守恒,给解算带来不少方便,可遗坐标的个数越多,自然就越方便。同一个问题,由于广义坐标选取方法不同,可遗坐标的有无或个数随之就不同。例如,研究有心力问题时,如果选用直角坐标,那就没有一个可遗坐标,但如果变化到平面极坐标系,则极角 \(\varphi\) 就是可遗坐标。(9.1.1 正则变换的条件)
有了这两个关系,参见 4.1 节,不难求出轨道方程,但是我们也不必从有效势能出发,可以用比内公式做个转化,直接从相互作用力的形式出发,求出轨道方程。
平方反比力给出的轨道方程是圆锥曲线。
引力 开普勒行星运动定律
著名平方反比引力:万有引力。引力情况离心率大于还是小于 1 对应相对运动机械能是正还是负。轨道方程 \(\displaystyle \rho=\frac{p}{1+\epsilon\text{cos}(\sigma+C)}\),其中离心率 \(\displaystyle \epsilon=\sqrt{1+\frac{2E^\prime h^2m^\prime}{m^2k^4}}\)
平方反比引力下的椭圆轨道:半长轴 \(\displaystyle a=\frac{mk^2}{2\vert E^\prime\vert}\)、半短轴 \(\displaystyle b=h\sqrt{\frac{m^\prime}{2\vert E^\prime\vert}}\)、能量 \(\displaystyle E^\prime=-\frac{mk^2}{2a}\)。
引力下的圆轨道:稳定性 \(\displaystyle \frac{F^\prime(r_0)}{F(r_0)}+\frac{3}{r_0}>0\) -> 稳定。星际航行怎么行。
斥力 散射问题
斥力情况(由于在两题问题中规定距离无穷远相互作用势能为 0)相对运动机械能必为正,离心率必大于 1,也就对应双曲线轨道。
原子的结构? -> 卢瑟福 \(α\) 粒子散射实验 -> 核式模型
卢瑟福 α 粒子散射实验:先用套路在质心系求解相对运动散射角 \(\theta\),然后转化到静止(实验室)坐标系求散射角 \(\Theta\)。关于散射后的例子关于散射角 \(\theta\) 的分布,导出散射截面的卢瑟福公式 \(\displaystyle \sigma(\theta)=\frac{1}{4}(\frac{ZZ^\prime e^2}{2eE^\prime})\frac{1}{\text{sin}^4(\frac{\theta}{2})}\)。然后问题来了,怎么转化成实验室坐标系下的散射截面 \(\sigma(\Theta)\)?太复杂了!不过别担心,卢瑟福的实验里金原子核和 \(\alpha\) 粒子的质量差别很大,来个近似。
小振动
把拉格朗日函数在平衡位置附近展开,保留二阶近似(现在不是很明白教材里“拉格朗日方程缺只是一阶近似”是啥意思)。
概念:简正坐标、简正频率、模式。无需赘述。
小振动的一般理论 - 矩阵表述。其中教材 5.3.5 节写道,\(\mathbf M\) 与 \(\mathbf K\) 是两个对称矩阵,而且 \(\mathbf M\) 是正定的,于是按照矩阵理论,一定可以找到一个实的非奇异满秩代换使 \(\mathbf M\) 和 \(\mathbf K\) 同时对角化。请问这个矩阵理论在哪?参见 Theorem 7.6.1, Horn, Roger A.; Johnson, Charles R. (2013). Matrix Analysis, second edition. Cambridge University Press。
刚体力学
力学中我们已经学过,刚体的定点转动的三个欧拉角可以确定刚体的运动,于是可以选取欧拉角为广义坐标。如果主动力全是保守力,则可以写出拉格朗日函数。例:陀螺,参见梁老师 6.5 对称重刚体的定点运动(拉格朗日 - 泊松情况)。重力矩与进动的抗衡。
哈密顿力学
拉格朗日动力学用广义坐标和广义速度描写力学系统的运动,哈密顿动力学则用广义坐标和广义动量描写力学系统的运动,后者与广义坐标的量纲无关,总是作用量(能量 \(\times\) 时间的量纲)。而量子论的量子化条件正是把作用量加以量子化。可见哈密顿动力学可作为从经典力学到量子力学的“跳板”。
哈密顿正则方程
此部分下文均讨论只有完整约束,并且主动力都是具有势能或广义势能的,即拉格朗日方程为 \(\displaystyle\frac{d}{dt}\frac{\partial L}{\partial\dot q_\alpha}-\frac{\partial L}{\partial q_\alpha}=0\)。
勒让德变换:,\(L(q,\dot q,t)\) -> \(\bar L(q,p,t)\),运动规律微分方程(原为拉格朗日方程)如何改变?构造 \(H=\displaystyle\sum_{\beta=1}^sp_\beta\dot q_\beta-L\),这个 \(H\) 我们曾在“运动积分”见过的,然而独立变量的选择发生了改变。运动规律成为看起来很好看的哈密顿正则方程。哈密顿力学中,以广义动量、广义坐标、时间为变量的系统中,运动规律由哈密顿正则方程给出。例:哈密顿力学。
\[\begin{cases} \displaystyle\frac{\partial H}{\partial q_\alpha}=-\dot p_\alpha \\ \displaystyle\frac{\partial H}{\partial p_\alpha}=\dot q_\alpha \quad(\alpha=1,2,...,s). \end{cases}\]- P.S. “正则”的意思是形式简单而对称。对称到打字靠复制粘贴再改改都很方便。
刘维尔定理
在经典力学问题求解中,如果用哈密顿力学求解,总是逃不开拉格朗日函数的求解。那么我们为什么要用哈密顿力学呢?对于这个问题,梁昆淼老师在其《理论力学》教材第 7-9 章深入探讨。从一切可能运动的拉格朗日函数 \(L(q,\dot q,t)\) 出发(拉格朗日动力学)或从一切可能运动的哈密顿函数 \(H(q,p,t)\) 出发(哈密顿动力学),此时 \(q\) 和 \(\dot q\) 相互独立,\(q\) 和 \(p\) 也相互独立。但是进入虚位移求实际运动时 \(\dot q\) 就不再独立于 \(q\)(我的理解是”关联“指的就是拉格朗日方程?),然而 \(p\) 却独立于 \(q\)(这句话详见下文哈密顿原理部分)。于是拉格朗日动力学在 \(s\) 维位形空间研究,哈密顿力学在 \(2s\) 维位形空间中研究。
统计物理 - 统计系综理论(刘维尔定理)。把广义坐标和广义动量当作直角坐标而构成 \(2s\) 维的空间叫做相空间。哈密顿正则方程的重要性:刘维尔定理是 \(2s\) 维的相空间中的定理,在普通空间或 \(s\) 维的位形空间中并不存在类似的定理。
位力定理(Virial Theorem)
The word virial for the right-hand side of the equation derives from vis, the Latin word for “force” or “energy”, and was given its technical definition by Rudolf Clausius in 1870.[1] —Wikipedia
梁老师书中把这个定理放在哈密顿力学刘维尔定理之后,但我不是很明白和哈密顿动力学有什么关联,我个人把它放在这里。
\(S=\displaystyle\sum_{i=1}^n\vec K_i\cdot\vec r_i\) 有界(\(\vec K_i\) 指动量),取 \(\displaystyle\lim_{\tau\rightarrow\infty}\langle\frac{dS}{dt}\rangle\),得 \(\displaystyle\langle T\rangle=-\frac{1}{2}\langle\sum_{i=1}^n\vec F_i\cdot\vec r_i\rangle\)。
【GRAD-UPDATE】研究随机过程统计量的系综平均用到的数学手段也很像,详见 (En) hilfinger2011separating。
例:两粒子有心力 \(r^n\) 作用,势能满足 \(V=kr^{n+1}\),位力定理给出 \(\displaystyle\langle T\rangle=\frac{n+1}{2}\langle V\rangle\)。对于平方反比有心力,\(\displaystyle T=-\frac{1}{2}\langle V\rangle\)。解答了我高中时代的疑惑。
例:推导理想气体压强方程 \(pV=nkT\)。
泊松括号
把任意力学量用广义坐标和广义动量表示,\(\varphi(p,q,t)\)。两个力学量的泊松括号为 \([\varphi,\psi]=\displaystyle\sum_{\alpha=1}^s(\frac{\partial\varphi}{\partial q_\alpha}\frac{\partial\psi}{\partial p_\alpha}-\frac{\partial\varphi}{\partial p_\alpha}\frac{\partial\psi}{\partial q_\alpha})\)。
于是,力学量随时间变化率 \(\displaystyle\frac{d\varphi}{dt}=\frac{\partial\varphi}{\partial t}+[\varphi,H]\)。
基本泊松括号 \([q_\alpha,q_\beta]=0\),\([p_\alpha,p_\beta]=0\),\([q_\alpha,p_\beta]=\delta_{\alpha\beta}\)。
泊松定理 可积系统
运动积分 \(\varphi\) 满足 \(\displaystyle\frac{d\varphi}{dt}=0\) 即 \(\displaystyle[\varphi,H]=-\frac{\partial\varphi}{\partial t}\)。
雅可比恒等式 \([\varphi,[\psi,\theta]]+[\psi,[\theta,\varphi]]+[\theta,[\varphi,\psi]]=0\)。
泊松定理:若 \(\varphi\) 和 \(\psi\) 都是运动积分,将它们和 \(H\) 代入雅可比恒等式,得到 \([\varphi,\psi]\) 也是运动积分。
刘维尔定理:对于有 \(s\) 个自由度的力学系统,如果可以找到 \(s\) 个相互独立的运动积分 \(\phi_1,\phi_2,...,\phi_s\),且满足 \([\phi_\alpha,\phi_\beta]=0(\alpha,\beta=1,2,...,s)\),这样的系统称为可积系统,此时称 \(\phi_1,\phi_2,...,\phi_s\) 是相互对合的。
系统可积时,其哈密顿函数属于运动积分之一或者可以用它们表示出来,而系统的运动方程可以通过积分的方法求解。
具体证明书上没写我也没自己推导,相关词条:Integrable system、Liouville–Arnold theorem。
可积系统举例:对称重刚体的定点运动、各向同性有心力作用下的质点、保守完整的小振动系统……
不可积系统举例:绝大多数多体问题,e.g., 著名的三体问题。
可积系统的运动是确定论性的,即给定初始条件就可以知道系统过去和未来的运动情况;同时,如果初始条件有微笑得到变化,则系统的运动情况在变化前后的差别也是微小的,也就是运动轨道是稳定的;可积系统的运动限制在相空间的 \(s\) 维环面上做规则运动。——10.5
力学变分原理
Significance: 见梁老师本章引言。
哈密顿原理
哈密顿原理:力学系统的动力学归结为一个变分原理:力学系统从时刻 \(t_1\) 到时刻 \(t_2\) 的一切可能运动之中,使哈密顿作用量 \(S=\displaystyle\int_{t_1}^{t_2}L(t,q,\dot q)dt\) 取极值的运动才是实际发生的运动。
位形空间的哈密顿原理:\(\delta \displaystyle\int_{t_1}^{t_2}L(q,\dot q,t)dt=0\) <-> \(\displaystyle\int_{t_1}^{t_2}\displaystyle\sum_{\alpha=1}^s(\frac{d}{dt}\frac{\partial L}{\partial\dot q_\alpha}-\frac{\partial L}{\partial q_\alpha})\delta q_\alpha=0\)(拉格朗日方程).
相空间的哈密顿原理:\(\delta \displaystyle\int_{t_1}^{t_2}[\displaystyle\sum_{\alpha=1}^s p_\alpha\dot q_\alpha-H(p,q,t)]dt=0\) <-> \(\displaystyle\int_{t_1}^{t_2}\displaystyle\sum_{\alpha=1}^s [(\dot q_\alpha-\frac{\partial H}{\partial p_\alpha})\delta p_\alpha-(\dot p_\alpha-\frac{\partial H}{\partial q_\alpha})\delta q_\alpha]dt=0\)(哈密顿正则方程).
变分原理的等价:如果上面两个哈密顿原理里的被积函数加上某个函数 \(U\) 对时间的全导数 \(\displaystyle \frac{dU}{dt}\),如果代表点的初始和终末位形是给定的,则 \(\delta U\vert_{t_1}=\delta U\vert_{t_2}=0\),不改变动力学方程。(后面哈密顿 - 雅可比方程要用)
位形世界的哈密顿原理:将时间 \(t\) 看作和广义坐标 \(q\) 地位平等的变量,\(t\) 相当于第 \(s+1\) 个广义坐标,这 \(s+1\) 个广义坐标构成了推广的位形世界。……参数 \(\tau\),\(\displaystyle\dot q_\alpha=\frac{q_\alpha^\prime}{t^\prime}\)……修改的拉格朗日函数 \(\displaystyle\Lambda=L(q,t,\frac{q^\prime}{t^\prime})t=\Lambda(q,t,q^\prime,t^\prime)\),于是哈密顿原理成为 \(\delta\int\Lambda d\tau=0\)。
最小作用量原理
实际上,还有其他形式的最小作用量原理和哈密顿原理相似,最小作用量的名称并不意味着真是运动的作用量永远是最小,而只说明,真实运动的作用量稳定值。只有当积分长度是一个相当小的有限量时,真实运动的作用量才取极小值。
如果拉格朗日函数 \(L\) 不显含时间(保守系统),修改后的拉格朗日函数 \(\Lambda\) 中 \(t\) 是可遗坐标……-> \(\displaystyle\Delta\int\sum_{\alpha=1}^sp_\alpha\dot q_\alpha dt=0\) 雅可比最小作用量原理
在 \(\Delta\) 的变分中,考虑的时哈密顿函数 \(H\) 守恒的一切可能的运动(但它们的 \(H\) 各不相同)。这样,\(t\) 不再可以任意给定,时间 \(t\) 的变分 \(\Delta t\) 不为零,特别是在两端点。
对于一般的力学系统 \(\Delta\int vds=\Delta\int\sqrt{2(E-V)}ds=0\)。
力学以外的变分原理:几何光学中也有条变分原理,即费马原理 \(\displaystyle\delta\int\frac{1}{u}ds=\frac{1}{c}\delta\int nds=0\)。几何光学是波动光学的短波长极限,这就启发我们也可以把经典力学当作某种波动力学的短波长极限。
正则变换 哈密顿-雅可比方程
揭示了向量子力学的过渡……
哈密顿方程对于范围更广的变换(正则变换)具有不变性。因此可以借助于这种变换使哈密顿函数变得极简单(变为常数或零),以使求解变换后的哈密顿方程成为极简单的事。
此章例子见正则变换与谐振子.
正则变换
哈密顿原理 \(\delta \displaystyle\int_{t_1}^{t_2}L(q,\dot q,t)dt=0\)(位形空间中)或者 \(\delta \displaystyle\int_{t_1}^{t_2}[\displaystyle\sum_{\alpha=1}^s p_\alpha\dot q_\alpha-H(p,q,t)]dt=0\)(相空间中),被积函数如果加上初末值给定的一个时间全微分 \(U\),则变分为 \(\delta\displaystyle\int_{t_1}^{t_2}\frac{dU}{dt}dt=0\),并不改变动力学方程。
我们希望对哈密顿函数里的广义坐标和广义动量做变换,
\[\begin{cases} P_\alpha=P_\alpha(p,q,t)\\ Q_\alpha=Q_\alpha(p,q,t)\quad(\alpha=1,2,...,s). \end{cases}\]得到一个与之等价的函数 \(K(P,Q,t)\),变换后满足的变分原理与原来等价,
\[\begin{cases} \delta \displaystyle\int_{t_1}^{t_2}[\displaystyle\sum_{\alpha=1}^s p_\alpha\dot q_\alpha-H(p,q,t)]dt=0\\ \delta \displaystyle\int_{t_1}^{t_2}[\displaystyle\sum_{\alpha=1}^s P_\alpha\dot Q_\alpha-K(P,Q,t)]dt=0. \end{cases}\]如何实现?按照上面的说法,\((\displaystyle\sum_{\alpha=1}^s p_\alpha\dot q_\alpha-H)-(\displaystyle\sum_{\alpha=1}^s P_\alpha\dot Q_\alpha-K)=\frac{dU}{dt}\),即 \(\displaystyle\sum_{\alpha=1}^s p_\alpha dq_\alpha-\displaystyle\sum_{\alpha=1}^s P_\alpha dQ_\alpha+(K-H)=dU\) 即可,这就是正则变换。我们希望经过这个变换,能变换出来尽可能多的可遗坐标和相应守恒的“广义动量”。
怎么变?应用勒让德变换把宗量给变了。见本章例子。
泊松括号不变性:在正则变换下,\([\varphi,\psi]_{pq}=[\varphi,\psi]_{PQ}\)。
无限小正则变换:选母函数 \(\displaystyle U_3(q,P,t)=\sum_{\alpha=1}^sq_\alpha P_\alpha+\epsilon G(q,P)\),正则方程变为
\[\begin{cases} \displaystyle dq_\alpha=\epsilon\frac{\partial G}{\partial p_\alpha}\\ \displaystyle dp_\alpha=-\epsilon\frac{\partial G}{\partial q_\alpha}\quad(\alpha=1,2,...,s). \end{cases}\]- 取 \(G=H\),\(\epsilon=dt\),则上面正则方程就是哈密顿正则方程。
- 若 \(H\) 不显含时间,\(dH=\epsilon[H,G]\). 根据上文泊松括号中力学量随时间的变化率,若 \([H,G]=0\) -> \(dH=0\) -> \(G\) 是运动积分。
- 例子:有心力作用下的哈密顿量 \(H\) 取绕某轴转动的角动量为母函数 \(G\),算得 \([H,G]=0\)。
哈密顿 - 雅可比方程
如何找到上文中能把正则变量都变成常数(初始值)的正则变换?我们希望 \(K(P,Q,t)=0\),正则方程就保证了正则变量是常数作为变换后的“动量”。在哈密顿函数代入 \(K=0\) 的条件,把原来的一阶 linear system of ODEs 变换为一阶 a nonlinear PDE,这就是哈密顿-雅可比方程,方程的解叫做哈密顿主函数 \(S\),其实就是拉格朗日函数 \(S=\displaystyle\int_{t_1}^{t_2}Ldt\)。
当 \(H\) 不显含时间,可以把哈密顿-雅可比方程的空间变量与 \(t\) 分离,令 \(S(q,P,t)=W(q,P)+f(t)\),得到 \(\displaystyle H(q,\frac{\partial W}{\partial t})=-f^\prime(t)=E\)。也称前面这个方程为哈密顿-雅可比方程,\(W(q,P)\) 为哈密顿特征函数。
- 对于自由度多于 \(1\) 的系统,一种特殊的可积系统(泊松定理):完全可分离的系统,哈密顿函数是 \(s\) 个独立部分的和, \(H(q_1,...,q_s,p_1,...,p_s)=\displaystyle\sum_{\alpha=1}^s H_\alpha(q_\alpha,p_\alpha)\)…… 哈密顿特征函数可写为 \(\displaystyle W=\sum_{\alpha=1}^sW_\alpha(q_\alpha,E,C_2,\ldots,C_s)\),哈密顿 - 雅可比方程可以分离为 \(s\) 个 \(\displaystyle\Phi_\alpha(q_\alpha,\frac{dW_\alpha}{dq_\alpha},E,C_2,\ldots,C_s)=0\) 积分求解以得到特征函数。
书中给出了一些用哈密顿-雅可比方程对力学系统进行正则变换求解然后再变回去的例子。其中谐振子的问题很有意思。变换后 \(\displaystyle K(P,X)=\frac{1}{2\pi}\sqrt{\frac{k}{m}}P\)。在这个例子里我们发现,能量和时间其实就是对系统做正则变换后的一对共轭变量,其实就是变换后的哈密顿函数的一对“动量”(能量守恒)和“坐标”(\(K\) 中不显含 \(X\),\(X\) 是可遗坐标)。
在平方反比有心吸引力的例子里,提到哈密顿 - 雅可比方法总是直接给出轨道。
正则微扰理论
在力学中,有许许多多的问题是无法精确求解的。但是经常遇见这样的情况,这些无法精确求解的问题与能够严格求解的问题之间,仅有微小的差别,因而可以在后者严格解的基础上近似求解前者,这就是微扰理论。建立在正则变换基础上的正则微扰理论就是微扰理论的一种。
\(H=H_0+\epsilon H^\prime\)。蒙上眼睛不看微扰,解 \(H_0\) 假定可以按哈密顿-雅可比方程求解出主函数 \(\displaystyle S(q_1,q_2,...,q_s;\frac{\partial S}{\partial q_1},\frac{\partial S}{\partial q_2},\frac{\partial S}{\partial q_s};t)\),解出 \(S(q_1,q_2,...,q_s;\alpha_1,\alpha_2,...,\alpha_s;t)=0\) 一堆积分常数 \(\alpha_\alpha,\beta_\alpha\) 满足 \(\displaystyle\frac{\partial S}{\partial \alpha_\alpha}=\beta_\alpha\)。现在考虑微扰,这些常数就不再是常数,但是以 \(S\) 为母函数,正则变换的关系不变,原来求出来的常量是正则变量,满足正则方程
\[\begin{cases} \displaystyle\dot\alpha_\alpha=-\frac{\partial(\epsilon H^\prime)}{\partial\beta_\alpha},\\ \displaystyle\dot\beta_\alpha=-\frac{\partial(\epsilon H^\prime)}{\partial\alpha_\alpha}. \end{cases}\]满足以微扰函数 \(\epsilon H^\prime\) 为哈密顿函数的正则方程。这个方程一般称为微扰方程。
例:非线性振动一阶微扰修正(见本章附件)。我们看到非线性谐振ß子的振幅修正为零,频率修正后为 \(\displaystyle\omega_0(1+\frac{3b}{8k}A^2)\)。由于谐振子能量为 \(\displaystyle\frac{1}{2}m\omega^2A^2\),不难看出能量修正后为 \(\displaystyle E_0+\frac{3}{8}A^4b+o(b)\)。而非微扰解 \(x=A\sin(\omega_0(t+\beta))\) 在微扰 \(bx^4\) 下的期望 \(bA^4\langle\sin^4(\omega_0(t+\beta))\rangle\) 就是 \(\displaystyle\frac{3}{8}bA^4\)。
作用量变量与角变量
物理学中的周期性运动中,我们关注其频率。分为天平动和转动(回忆单摆相图,天平动是在底下来回摆,转动是哇哇绕大圈)。现只讨论哈密顿特征函数 \(W\) 可分离为 \(\displaystyle\sum_{\alpha=1}^s W_\alpha(q_\alpha;E,C_2,...,C_s)\) 的情况。变换后的“动量”常数取为对一个周期积分的作用量 \(J_\alpha=\oint p_\alpha dq_\alpha\),\(J_\alpha\) 只是积分常数 \(E,C_2,...,C_s\) 的函数,于是 \(W=\displaystyle\sum_{\alpha=1}^s W_\alpha(q_\alpha;J_1,J_2,...,J_s)\)。把 \(J_\alpha\) 的共轭变量 \(\omega_\alpha\) 叫做角变量满足 \(\displaystyle\dot\omega_\alpha=\frac{\partial E}{\partial J_\alpha}\),是 \(q_\alpha\) 周期运动的频率。
例:一维谐振子。\(\displaystyle E=(\frac{1}{2\pi}\sqrt{\frac{k}{m}})J\)。量子力学过渡:量子化作用量 \(J=nh(n=0,1,2,...)\),则 \(\displaystyle E=n\frac{h}{2\pi}\omega\)。不过,氢原子量子力学解并不需要这种认为量子化,量子化是薛定谔方程的结果。为啥关注作用量?下面会讲,浸渐过程下作用量是浸渐不变量。
例:平方反比有心力吸引下的质点周期运动频率。延申:各个坐标周期运动如果有不同的运动频率……
浸渐不变量与哈内角
如果哈密顿函数含有一个随时间变化的参数 \(\lambda\),且这个参数随时间的变化非常缓慢,“缓慢”是指在系统运动周期 \(\tau\) 的时间内,\(\lambda\) 的相对变化很小,即 \(\displaystyle\tau\frac{d\lambda}{dt}\ll\lambda\),如有某力学量于此条件下不随时间改变,这个力学量就称为浸渐不变量。
“浸渐”译自 adiabatic,这里强调的是上式所表明的缓慢性。这同一个词在热学中则译为“绝热”。热学中的绝热过程,其变化相对于系统趋于平衡态的弛豫过程是缓慢的,以至于可以认为系统时时处于平衡态,但另一方面,其变化相对于热传导又是迅速的,以致可以认为热量来不及进出系统。在热学中使用“绝热”一词往往强调其迅速的一面,其实它还有缓慢的一面。
例:摆长缓慢变化的小振幅单摆能量乘以周期为浸渐不变量。高中解法见 Elementary examples of adiabatic invariance, II. Simple Pendulum,其实用的就是功能原理。正则变换求解见本章附件。
讨论了 \(J\) 的不变性,与其共轭的变量 \(\omega\) 怎么变?如果参数缓慢变化,然后又变回去了,系统高频周期运动的相位还和原来一样吗?并不,累积的相位差就是哈内角。
从“几何力学”到“波动力学”
先自己看看吧,我有的东西需要再看看再写。
文档信息
- 本文作者:L Shi
- 本文链接:https://SHI200005.github.io/2022/04/30/Advanced-Classical-Mechanics/
- 版权声明:自由转载-非商用-非衍生-保持署名(创意共享3.0许可证)