1. 为什么系列走到这里必须聊分析力学写这个系列之前五篇的时候我一直把重点放在牛顿力学的框架里受力分析、动量定理、角动量定理、能量守恒、刚体的平面运动、转动惯量的计算……这些都是“力”的语言。到了第六篇该换一套叙述方式了。这一篇的核心关键词是广义坐标、拉格朗日方程、哈密顿正则方程、循环坐标、守恒量。说得再直白一点从牛顿力学进入分析力学。很多同学学理论力学读到分析力学这一章会有一种“突然从物理课变成数学课”的错觉——到处都是偏导数、变分、泛函拉着拉着就不知道自己在算什么物理问题。其实分析力学并没有引入新的物理规律它只是换了一副框架来描述同一个世界。这个框架的好处在于对于约束多、自由度多、系统复杂的动力学问题分析力学往往能比牛顿力学更机械、更省脑地列出运动方程。这篇文章不打算重复教材里的全部推导而是想把分析力学里最容易让人卡壳的几个点拆开讲约束和广义坐标怎么选、虚位移到底“虚”在哪、拉格朗日方程为什么长那个样子、循环坐标和守恒量怎么找、哈密顿方程和拉格朗日方程的切换逻辑是什么。最后还会聊聊考试和做题时常见的坑。老规矩凡是教材上写得很详细的标准推导我尽量用“为什么是这样”的思路带你过一遍凡是教材上含糊带过、但实际做题时总出问题的地方我多花点笔墨。2. 约束、自由度与广义坐标先选好“坐标”再做力学2.1 约束方程怎么数自由度怎么算分析力学和新牛顿力学最根本的区别从“选坐标”这一步就开始了。牛顿力学里我们习惯用笛卡尔坐标系比如一个质点的位置是(x, y, z)三个坐标决定了它在空间中的位置。但很多系统里的质点并不能任意动比如单摆的小球被一根轻绳拉着它只能在一个圆弧上运动再比如刚体上的一个点它和刚体上其他点的距离保持不变。这类对系统位置的限制在分析力学里叫约束。一个系统里如果有N个质点每个质点在空间中都要3个坐标来确定位置所以一共有3N个坐标。但约束方程会建立坐标之间的关系把独立的坐标数量减下来。系统的自由度 3N - 独立约束方程的个数。自由度的含义就是要唯一确定一个系统的位形最少需要多少个独立的坐标参数。举个例子一根杆子在平面里运动杆上有两个端点每个端点是2个坐标平面运动x和y所以总坐标数是4杆长固定是一个约束方程所以自由度是4 - 1 3。事实也确实如此杆子质心的位置有2个坐标再加一个杆的转角就能完全描述杆的位形——这就是广义坐标的一个具体例子。这里有个常见误区约束方程有没有写独立。比如平面上一个质点被限制在一个圆周上运动如果只写x² y² R²这是一个约束方程自由度是2 - 1 1。有的同学会再写一个θ arctan(y/x)结果把自由度算成0——但这个式子其实是把同一个约束换了个写法它并没有新的独立性不能重复计入。2.2 广义坐标的选取原则自由度确定之后就要选广义坐标。广义坐标可以是长度、角度、弧长甚至可以是一个比例系数只要能唯一地、方便地确定系统位形就行。选广义坐标的基本原则有三条坐标数等于自由度坐标之间必须独立不能再有约束关系夹在里面同一个系统可以有多种广义坐标的选法但方程的形式和质量会不一样。我见过太多同学在选广义坐标上翻车不是选少了漏掉一个转动自由度就是选多了把约束方程又写出来当坐标用最后拉格朗日方程根本列不干净。给个经典例子一个圆盘在斜面上做纯滚动。如果只选圆盘质心沿斜面的位移x作为广义坐标你会损失圆盘的转角θ因为纯滚动条件x Rθ其实是一个约束。这里正确的自由度数应该是2质心位移加转角减去1纯滚动约束等于1所以选x作为广义坐标就足够了转角通过纯滚动条件表达出来不用再单独设变量。但如果你一开始选了两个坐标x和θ又同时引用了纯滚动条件再列方程时就要注意方程的数量和约束力的处理——这就是非完整约束与完整约束的差异。2.3 完整约束与非完整约束约束分两大类完整约束和非完整约束。完整约束是指约束方程可以写成坐标和时间的函数关系比如f(x, y, z, t) 0它的特点是约束把系统的位形空间直接限制到某个低维曲面上。前面说的绳长固定、杆长固定、纯滚动位置约束形式都属于完整约束。非完整约束是指约束方程是速度之间的关系比如某些滚动系统的约束只能写成dx - R sinθ dφ 0这种微分形式无法积分成有限的坐标关系式。这类约束在分析力学里处理起来更麻烦第一类拉格朗日方程乘子法和第二类拉格朗日方程最小坐标法的区别也正是针对这两种约束问题。对于非完整约束达朗伯原理结合拉格朗日乘子是比较常见的处理方式。对应到教材上就是“带约束乘子的拉格朗日方程”。很多同学看到书上出现好几种形式的拉格朗日方程就懵了。我的建议是先搞清楚系统究竟是完整约束还是非完整约束再决定用哪一套公式。90%的课程题目都是完整约束别自己把问题想复杂。提示做题之前先问自己一句话——“这个系统的约束方程能不能写成坐标与时间的关系不出现速度”能就是完整约束不能才需要动用乘子法。3. 虚位移原理与达朗贝尔原理分析力学的两块基石3.1 虚位移不是真的位移从牛顿力学到分析力学最让人不适应的一步就是“虚位移”这个概念。一个系统在某一瞬间所有约束都“冻结”了不随时间变化在这种情况下系统允许发生的、无限小的、满足所有约束条件的位移叫虚位移。关键不是“位移”而是“虚”——它不是真实发生的位移而是想象中可能的位移。真实位移是动力学方程组解出来的结果有确定的大小和方向虚位移是“每个瞬间都重新假想一遍”的试探位移。怎么理解想象一个单摆小球被绳拉着真实位移是沿着圆弧的一小段弧长虚位移是“假设那一瞬间绳长不变、所有外力都不变小球在约束允许的方向上动那么一丁点”的位移。因为单摆的约束是完整约束虚位移方向和切向速度方向重合很多场景下虚位移个数就等于自由度个数。虚位移原理说的是一个受约束的静止系统处于平衡态的充要条件是——所有主动力在任意虚位移上做的虚功之和为零。写作∑ Fᵢ · δrᵢ 0这里Fᵢ只包含主动力不含约束力。为什么约束力可以不算因为理想约束光滑接触、光滑铰链、不可伸长的轻绳等的约束力方向总是和虚位移方向垂直或者做功为零约束力的虚功恒等于零。这是整个分析力学链条里非常关键的一步把约束力的效应从动力学方程中“剥离”出去。3.2 达朗贝尔原理把动力学问题变成静力学问题虚位移原理处理的是平衡问题而达朗贝尔原理把动力学问题转化为“平衡”问题来处理思路非常巧妙。对牛顿第二定律 F ma把ma移到等号另一边写成F - ma 0。如果我们把-ma定义为一个虚构的“惯性力”那原来的动力学方程在形式上就变成了一个“平衡”方程主动力 约束力 惯性力 0。然后把虚位移原理扩展到动力学问题上在任意虚位移下主动力、约束力和惯性力的虚功之和为零。因为理想约束力虚功为零就得到达朗贝尔原理的积分形式∑ (Fᵢ - mᵢaᵢ) · δrᵢ 0这里的“惯性力”是个数学技巧不是真实存在的力但它给了我们一个威力无穷的工具——可以用虚功的方法来处理动力学问题绕开解约束反力。在实际应用中达朗贝尔原理对某些系统非常高效比如有几个刚体通过铰链连接的机构每个刚体的质心有加速度你要列整个系统的运动方程直接用牛顿第二定律每个刚体都列会牵扯大量铰链反力用达朗贝尔原理把这些反力通过虚功消掉整个方程一下就简洁很多。3.3 广义力在广义坐标下的“受力”既然有了广义坐标力的概念也得跟着升级。原来笛卡尔坐标下的力是矢量有方向有大小在广义坐标qⱼ里我们用广义力Qⱼ来描述“某个广义坐标方向上的有效作用”。广义力的定义式是Qⱼ ∑ Fᵢ · (∂rᵢ / ∂qⱼ)如果主动力是保守力即力可以写成某个势能函数V(q)的负梯度那么广义力就是Qⱼ -∂V / ∂qⱼ这里注意V要写成广义坐标和时间的形式不能还停留在笛卡尔坐标的V(x, y, z)。很多同学第一步就卡在这里——动能T很容易换到广义坐标但势能V忘了把约束关系代入。比如一个质点在抛物线上运动重力势能是mgy如果y是x的二次函数你就要把y用广义坐标x表示出来。4. 拉格朗日方程从虚功到运动方程的标准路线4.1 方程的推导逻辑把达朗贝尔原理和广义坐标结合起来经过一套标准的变分推导就能得到完整约束系统下的第二类拉格朗日方程d/dt (∂L/∂q̇ⱼ) - ∂L/∂qⱼ 0其中L T - V叫拉格朗日函数q̇ⱼ是广义速度。这个方程对每个独立的广义坐标qⱼ都成立。很多同学背下了这个式子但不清楚它的物理含义。我拆解一下∂L/∂q̇ⱼ叫作广义动量pⱼ。如果qⱼ是笛卡尔坐标pⱼ就等于m·v如果qⱼ是角度pⱼ就可能是角动量。广义动量是“对应坐标方向上的动量”。∂L/∂qⱼ叫作广义力。它反映了势能和其他因素在qⱼ方向上的变化。整个方程的意思是广义动量的时间变化率等于广义力。看到没有这其实还是牛顿第二定律的广义版本力 动量的变化率。只是我们把坐标系从笛卡尔换成了广义坐标把力换成了广义力。4.2 具体计算中的操作要点实际操作中列拉格朗日方程只需要五个步骤确定自由度选择广义坐标qⱼ写出系统的动能T用广义坐标和广义速度表达写出系统的势能V用广义坐标表达注意代入约束构造拉格朗日函数L T - V对每个广义坐标分别求∂L/∂q̇ⱼ、d/dt(∂L/∂q̇ⱼ)和∂L/∂qⱼ代入方程。但真正做题时最难的不是套公式而是“如何把T写出来”。举一个经典例子平面双摆两个单摆的端点铰接在一起。系统自由度是2选两个摆角θ₁、θ₂为广义坐标。第一个质点的速度是l₁θ̇₁的切向速度第二质点的速度却包含两部分θ̇₁引起的牵连速度和θ̇₂引起的相对速度两个速度方向不同必须用矢量合成。很多同学直接写成v₂ l₁θ̇₁ l₂θ̇₂这是代数加和错了——速度是矢量两者方向不同要用几何关系比如用极坐标或投影来合成。正确的做法是把第二个质点的位置坐标显式写出来然后对时间求导。x₂ l₁sinθ₁ l₂sinθ₂y₂ -l₁cosθ₁ - l₂cosθ₂求导后得到速度分量再平方相加自然得到含cos(θ₁ - θ₂)的交叉项。这个交叉项是双摆动能的标志性特征也是许多同学算错的地方。4.3 耗散系统怎么处理实际工程系统很少没有摩擦。如果系统里有摩擦力且摩擦力可以写成瑞利耗散函数D (1/2)∑ cᵢ q̇ᵢ²的形式线性阻尼那么拉格朗日方程修正为d/dt (∂L/∂q̇ⱼ) - ∂L/∂qⱼ ∂D/∂q̇ⱼ 0这里D的下标是广义速度物理意义是阻尼力在广义坐标方向上做的功率的二分之一。日常做题里最常用到的就是弹簧-质量-阻尼系统L包含弹簧势能和质量的动能D包含阻尼器的耗散项。按这个修正方程列式比硬把阻尼力当主动力去算广义力省事得多。5. 循环坐标、守恒量与首次积分找对称性比硬解方程更重要5.1 循环坐标就是“看不见的坐标”如果拉格朗日函数L中不显含某个广义坐标qⱼ只含它的广义速度q̇ⱼ那么这个坐标就叫循环坐标对应的广义动量pⱼ ∂L/∂q̇ⱼ守恒。直接用拉格朗日方程看d/dt (∂L/∂q̇ⱼ) ∂L/∂qⱼ如果∂L/∂qⱼ 0那么d/dt(∂L/∂q̇ⱼ) 0说明∂L/∂q̇ⱼ是常数。这就得到了一个首次积分——一个包含广义速度但不含广义加速度的守恒式。很多系统都有明显的循环坐标。比如一个不受外力作用的自由质点它的L (1/2)m(ẋ² ẏ² ż²) - V(x, y, z)如果势能V与某个坐标无关比如在均匀重力场中V mgzx和y就是循环坐标那么x和y方向的动量守恒。这正是“对称性导致守恒律”在拉格朗日框架下的体现——坐标平移不变性对应动量守恒角度旋转不变性对应角动量守恒时间平移不变性对应能量守恒。这组对应关系教材上叫诺特定理。了解它有几个很实际的价值做题时先看有没有循环坐标如果有可以直接降阶方程的求解难度减少一个量级。5.2 能量守恒和广义能量除了循环坐标对应的广义动量守恒还有一个常见的守恒量是能量也就是哈密顿函数H。一般地H ∑ q̇ⱼ ∂L/∂q̇ⱼ - L当L不显含时间t的时候H就是常数而且如果动能T是广义速度的二次齐次函数约束不显含时间时满足则H T V就是系统的总机械能。这里有一个很多同学容易犯迷糊的点L不显含时间不代表系统没有外力做功而是说明系统的约束不随时间变化且外力和势能也不依赖时间。此时系统的总能量才是守恒量。如果约束是随时间变化的比如一个长度随时间变化的单摆L里会出现显含t的项这时候H虽然依然可以计算但它不再是总机械能可能多出来一些“广义能量”的修正项。做题时我建议先判断L是否显含t如果显含别指望能量守恒老老实实解方程或者找别的守恒量。5.3 怎样寻找首次积分寻找首次积分是理论力学考试的经典套路也是很多同学觉得“看答案好简单自己做不出”的痛点。这里给一套可复用的检索顺序检查L中是否缺少某个坐标——如果有立即写出对应广义动量守恒。检查L是否显含t——如果否立即写出能量守恒或广义能量守恒。如果系统是中心力场势能只依赖r角动量守恒。如果系统存在对称轴比如势能不依赖某个角度θ则对应的角动量分量守恒。找到首次积分后拉格朗日方程就可以降阶。比如一个二维中心力场问题用极坐标L不显含θ那么角动量p_θ守恒方程立刻变成关于r的一维问题再用能量守恒就只剩一个一阶方程——求解难度大幅下降。6. 哈密顿正则方程从二阶微分方程到一阶方程组6.1 勒让德变换的逻辑拉格朗日方程是关于广义坐标的二阶微分方程组求解大系统时计算量很大。哈密顿想了个办法把广义速度换成广义动量把二阶方程降为一阶方程。做法是通过勒让德变换定义哈密顿函数H(q, p, t) ∑ pⱼ q̇ⱼ - L(q, q̇, t)形式上只是一个函数替换意义在于拉格朗日描述用的是(q, q̇)作为状态变量哈密顿描述用的是(q, p)作为状态变量。状态从N维升到2N维但方程从N个二阶方程变成2N个一阶方程。不要小看这个变换它是通往量子力学、统计力学和很多现代物理理论的桥梁。推导出来的哈密顿正则方程是q̇ⱼ ∂H/∂pⱼ ṗⱼ -∂H/∂qⱼ这个方程组高度对称而且形式极其优美广义速度等于哈密顿函数对广义动量的偏导广义动量变化率等于负的哈密顿函数对广义坐标的偏导。6.2 哈密顿函数的物理意义哈密顿函数在大多数物理背景下就代表系统的总能量H T V。但要注意这只是“通常情况”。严格来说当约束不显含时间且动能是速度的二次齐次函数时H才等于总机械能。对于显含时间的系统H同样可以计算但它代表的是“广义能量”和总机械能之间可能差一些与约束运动相关的项。举个例子一个圆环竖直放置一个小珠在环上滑动同时环自己绕竖直轴以恒定角速度ω旋转。取广义坐标为小珠与环心的连线跟竖直方向的夹角θ那么系统有两个贡献小珠相对环的滑动动能和环转动带来的牵连动能。你会得到一个显含时间或者显含广义速度二次项系数的L而H算出来可能不等于小珠的总机械能因为环的驱动不断向系统注入能量。不用担心哈密顿框架依然成立只是不要机械地认定H T V。6.3 什么时候用哈密顿方程很多同学会有疑问既然拉格朗日方程已经能解题为什么还要学哈密顿方程我的观点是对于基础的理论力学考试题拉格朗日方程往往已经够用哈密顿方程的优势更多体现在两个方向——数值计算上一阶方程组比二阶方程组更容易做数值积分标准的四阶龙格-库塔法直接作用在一阶形式上分析上哈密顿力学是相空间理论的起点李雅普诺夫稳定性、混沌、辛几何这些高阶内容都建立在哈密顿框架之上。如果只是求解一个双摆的小振动两套方法差别不大。但如果系统自由度到了几十上百比如分子模拟哈密顿形式才是标准选择。所以我的建议是考试时优先用你熟悉的框架拿分但练习题阶段一定把哈密顿方程也练熟后面学振动理论或分析力学进阶时需要它。7. 典型系统实操从单摆到双摆全步骤7.1 单摆最简单的验证单摆的约束是绳长固定自由度1选θ为广义坐标。位置坐标x l sinθy -l cosθ。速度v² l² θ̇²。动能T (1/2)ml² θ̇²。势能V -mgl cosθ选最低点为零势能点。所以L (1/2)ml² θ̇² mgl cosθ。代入拉格朗日方程∂L/∂θ̇ ml²θ̇d/dt ml²θ̈∂L/∂θ -mgl sinθ。得到ml²θ̈ mgl sinθ 0即θ̈ (g/l)sinθ 0。再看能量L不显含tH T V是常数这正是能量守恒。单摆是最简单的案例但所有方法都能在这里得到验证——如果你在更复杂的系统里算出来的结果不能退化回单摆形式那大概率是哪里出错了。7.2 平面双摆交叉项是分水岭双摆是理论力学课程里“劝退”很多人一个经典系统。两根轻杆长度l₁、l₂两个质点质量m₁、m₂广义坐标为θ₁、θ₂。质心1的位置和速度很简单重点是质心2的速度x₂ l₁sinθ₁ l₂sinθ₂ y₂ -l₁cosθ₁ - l₂cosθ₂求导ẋ₂ l₁θ̇₁cosθ₁ l₂θ̇₂cosθ₂ ẏ₂ l₁θ̇₁sinθ₁ l₂θ̇₂sinθ₂平方相加v₂² l₁²θ̇₁² l₂²θ̇₂² 2l₁l₂θ̇₁θ̇₂cos(θ₁ - θ₂)注意那个最后的交叉项cos(θ₁ - θ₂)——这是速度矢量合成的结果两个速度方向夹角恰好是(θ₁ - θ₂)。很多同学算错就是漏了它或把它算成了cos(θ₁ θ₂)后者完全错误。有了v₁² l₁²θ̇₁²和v₂²动能就是T (1/2)m₁v₁² (1/2)m₂v₂²。势能V -m₁gl₁cosθ₁ - m₂g[(l₁cosθ₁ l₂cosθ₂)]。现在可以写L了。拉格朗日方程严格代入后会得到两个互相耦合的二阶非线性方程(m₁m₂)l₁²θ̈₁ m₂l₁l₂θ̈₂cos(θ₁-θ₂) m₂l₁l₂θ̇₂²sin(θ₁-θ₂) (m₁m₂)gl₁sinθ₁ 0m₂l₂²θ̈₂ m₂l₁l₂θ̈₁cos(θ₁-θ₂) - m₂l₁l₂θ̇₁²sin(θ₁-θ₂) m₂gl₂sinθ₂ 0看着复杂但每一项都有清晰来源第一项是本杆的惯性项第二项是另一杆运动带来的牵连惯性项第三项是科氏/离心耦合项第四项是重力项。方程写成矩阵形式会看到质量矩阵里也有cos(θ₁-θ₂)的耦合这才导致双摆运动有极其丰富的动力学行为。7.3 带电粒子在磁场中运动一场速度依赖势的实战分析力学里还有一个很有趣的案例带电粒子在磁场中运动。磁场力不是保守力不能直接用U(q)来写势能但通过矢势A可以构造广义势U q(φ - v·A)拉格朗日函数可以写成L (1/2)mv² - qφ qv·A。这里真正有用的经验是当力与速度有关常规牛顿力学里“力是速度的函数”很难处理时分析力学依然可以给出漂亮的拉格朗日函数。比如一个在均匀磁场B B₀ẑ中运动的带电粒子取矢势A (1/2)(-B₀y, B₀x, 0)L里会出现与速度线性耦合的项而计算广义动量时粒子的机械动量会多出一项qA——这正是磁场中“正则动量不等于机械动量”的来源。这一块在很多教材里讲得比较快初学时可以放一放但如果你以后要学电动力学或量子力学一定要回来把这一节彻底弄懂。正则动量、规范变换、相空间这些概念在量子力学里会反复碰到。8. 做题易错点与判断技巧批改作业十年里最常见的几类坑8.1 广义速度求动能时忘记约束代入动能T里既有广义速度也可能隐含广义坐标。比如一个质点在半径为R的球面上运动很多人从笛卡尔坐标直接写v² ẋ² ẏ² ż²然后忘了用球面约束x² y² z² R²把其中一个坐标消掉导致T里含了三个速度分量结果自由度变成3方程多了一条且互相矛盾。判断办法列完T之后数一下T里出现了几个独立的“速度变量”如果大于自由度数一定还没代入约束。8.2 广义坐标选得不自然有的系统明明用角度作广义坐标最方便有人非要用某个长度再叠加热角函数表达最后方程里全是根号、反三角函数自己算了半天还总出错。我的建议是优先选择能直接反映约束对称性的坐标。转动的选角度平移的选位移滑块选沿导轨的位移不要为了“统一”强行用同一类坐标。8.3 拉格朗日函数里丢掉符号L T - V这个减号是分析力学里最容易出低级错误的地方。动能T是广义速度的二次型势能V是广义坐标的函数。如果势能V的零点和方向没设清楚符号错了后面整个运动方程的正负号全反。比如前面单摆例子如果把最低点势能写在V mgl(1 - cosθ)那L T - V方程里重力项符号依然正确但如果你写成V mglcosθ会造成方程里重力项符号整个反掉。做题时必须明确写出“零势能点在哪里正方向朝哪”。8.4 循环坐标看出来却不知道怎么用有同学找到循环坐标qₖ知道pₖ守恒但继续解方程时又把qₖ当活变量不降阶。正确做法是把qₖ的广义动量守恒式作为代数关系把q̇ₖ用其他广义坐标和常数表达出来然后代入剩余的方程中。这样方程里就不再出现qₖ或q̇ₖ自由度自动少一个。8.5 非完整约束当完整约束处理纯滚动问题是最典型的例子。一个圆盘在平面上纯滚动如果滚动方向是直线约束x Rθ是完整约束直接用最小坐标法就行。但如果圆盘沿平面做任意方向的运动像硬币在桌面上滚约束方程是速度分量之间的微分关系无法积分成有限形式这就是非完整约束。强行用第二类拉格朗日方程会得到错误的方程。看到“纯滚动且运动方向不固定”时务必警惕这个坑。注意考试时如果在题目里看到“纯滚动”“无滑动”先动手写约束关系再决定能不能用最小坐标法。多数课程考的都是直线纯滚但总有几个同学栽在上面的非完整处理上。8.6 拉格朗日方程的多余方程当广义坐标选多了比如选了3个坐标但自由度只有2时用第二类拉格朗日方程可能会得到3个方程但其中只有2个独立另外1个方程实际上包含了约束乘子。这时候如果你坚持3个方程都要满足解出来会互相矛盾。处理方式要么回到约束乘子法第一类拉格朗日方程要么重新调整坐标数目宁可多花一分钟检查也别在解题中途卡死。9. 一个完整的数值验证双摆两套方法对照只推导公式还不够我建议所有学习分析力学的同学都做一次数值验证——用同样的初值分别用牛顿力学笛卡尔坐标约束力和拉格朗日力学求解双摆对比结果。我自己带学生时特别推荐这么做因为“算出来两个结果对不上”是检验理解最好的方式。以一个m₁ m₂ 1、l₁ l₂ 1、初始θ₁ 0.2、θ₂ 0.1、初速为0的双摆为例。用拉格朗日方程做四阶龙格-库塔积分时间步长取0.005秒能明显看到两个摆角随时间呈不规则摆动而且对初值极其敏感。如果初值的小数点后第三位改变几十秒后轨迹就完全不同。牛顿力学角度怎么验证把两个质点的位置坐标设为(x₁,y₁)和(x₂,y₂)每根杆给出一个约束方程(x₂ - x₁)² (y₂ - y₁)² l₂²然后每根杆给一个未知的张力T₁、T₂对每个质点的x和y分别列牛顿第二定律加上约束方程的二阶导数形式一共5个方程解5个未知数。这种方法通用性极差推导复杂度高而且数值稳定性不如拉格朗日方法但它能作为结果校验。我实测的情况是两种方法在较短时间尺度内数值几乎一致误差在积分步长允许范围内。随着时间推移因为系统有混沌特性轨迹差异会逐渐变大——这不是方法错了而是初值敏感性在起作用。如果两个方法在早期就明显不一致那说明至少有一方列式有问题。10. 学习顺序和复习建议这一章到底怎么学才不吃力10.1 先打好牛顿力学基础再进入分析力学分析力学并不是牛顿力学的替代品而是更高维度的概括。如果你的牛顿力学本身就很薄弱受力分析画不清楚、动量定理和动能定理搞混那学分析力学只会雪上加霜。我见过不少同学直接抱着拉格朗日方程死磕结果连“虚功”都理解不了——因为虚位移的概念本质上是对约束几何的深刻理解而这来自对真实运动的分析积累。10.2 每做一个新系统先问三个问题我建议做题时先回答这三个问题再动笔这个系统的自由度是多少广义坐标选什么最自然系统的约束是完整的还是非完整的如果这三个问题都答得干脆恭喜后面就是套公式的体力活如果答不利索先去复习约束和广义坐标那部分别急着列方程。10.3 工程控制里的分析力学应用理论力学学到分析力学不是只为了考试。在很多工程方向分析力学才是真正的分析工具。比如机器人动力学多关节机械臂的每个关节就是一个自由度运动方程用拉格朗日方法列写比用牛顿-欧拉递推在某些场景下更直观尤其是需要分析动力学性质如能量、动量映射的时候。再比如多体系统仿真里约束方程和广义坐标的选取干脆就是算法的核心——你选得好不好直接决定仿真会不会发散。另外振动分析课程里大量使用拉格朗日方程推导多自由度系统的质量矩阵、刚度矩阵和阻尼矩阵。你会发现理论力学里学的广义坐标、动能矩阵、势能矩阵这些概念到了振动分析课程里直接对应到工程结构的有限元建模。矩阵对称性、质量矩阵的正定性这些都能在分析力学框架下得到统一理解。10.4 参考书和练习策略市面上的参考书很多我个人的建议是选一本你读起来最舒服的教材把它从头到尾啃完比同时翻五本教材但每本都读一半要好得多。练习题方面不要一开始就做复杂的双摆或陀螺进动先把单自由度系统练熟——弹簧振子、单摆、复摆然后用两到三个自由度系统验证自己对循环坐标和守恒量的理解最后再去碰那些需要完整推导的复杂系统。做题目时强烈建议坚持用“公式推导 文字说明”的方式写清每一步的理由。分析力学比牛顿力学更容易“算着算着就忘了自己在干什么”如果你能在每步旁边标一句话说明物理含义考试时就能少很多莫名其妙的符号错误。写在最后的经验回到这个系列第六篇的主题。从我给多届学生辅导的经验来看分析力学这一章最大的分水岭不是数学能力而是“物理图像是否清晰”。虚位移、广义坐标、循环坐标、正则动量这些概念看着抽象但它们每一个背后都有非常具体的物理直觉虚位移就是“符合约束的所有可能方向”循环坐标就是“拉格朗日函数里没出现的坐标”正则动量就是“广义坐标对应的动量分量”。把它们还原成物理图像公式就不再是一堆偏导数的机械操作。如果你读到这里觉得“框架清楚了”那建议立刻找三个系统——单摆、双摆、弹簧块组合——按本文的步骤各做一遍完整推导。做过这三道题分析力学就算真正入门了。后面再学哈密顿-雅可比方程、泊松括号、正则变换都是在今天这个框架上盖楼。理论力学中后期给人的感觉确实和前期完全不同前期是“看图说话受力分析”后期是“选坐标、找对称、写函数、列方程”。从牛顿到拉格朗日再到哈密顿本质上是一次比一次更“数学化”也是一次比一次更接近物理学的基本原理。如果你能把这一步跨过去以后面对分析力学的任何进阶内容都不会再觉得它是天书。