第 1–2 章 · 牛顿回顾与拉格朗日方程的诞生
Goldstein 的开局看似复习牛顿力学,实际在埋一个大转折:牛顿的矢量力学(力、加速度、逐个坐标列方程)处理约束系统极其笨拙,而拉格朗日力学用标量能量就能绕开约束力。这两章的任务是把转折的必要性讲透。
约束与广义坐标:扔掉多余变量
单摆其实是一维问题(角度 $\theta$),却被牛顿力学当作二维平面问题 + 一根绳子的约束力来解。广义坐标 $q_i$ 直接按自由度数参数化系统,约束"内建"在坐标选择里。约束分类:完整约束(可写成坐标关系式 $f(\mathbf{q},t)=0$)与非完整(速度不等式/不可积条件,如滚动圆盘)——拉格朗日方法对后者要格外小心。理想约束(约束力不做虚功)是使用拉格朗日方程的合法性条件。
虚功与达朗贝尔原理
虚位移 $\delta\mathbf{r}$:瞬时、想象、满足约束的位移。达朗贝尔原理:
理想约束的约束力在虚功中自动消失——这就是标量化方法的魔法来源。配上广义坐标,牛顿第二定律变形为:
拉格朗日方程
使用流程固定三步:选广义坐标 → 写 $T$ 与 $V$(都化成 $q, \dot q$ 的函数)→ 逐坐标列方程。守恒律自动显形:$L$ 不含某坐标(循环坐标)⟹ 对应的 $\partial L/\partial\dot q$(广义动量)守恒——对称性 ⟹ 守恒律的第一次正面亮相。
老师的提醒
- 列 $L$ 的头号坑:$T$ 必须用广义速度表达后再求偏导——把笛卡尔速度直接代入会漏掉坐标耦合项(极坐标 $T = \frac{m}{2}(\dot r^2 + r^2\dot\theta^2)$ 的 $r^2\dot\theta^2$ 项是试金石)。
- 拉格朗日方程不比牛顿"更准",它换的是效率与系统性:约束越多,优势越大。
自检
- 完整约束与非完整约束的区别?各举一例。
- 用拉格朗日方程解单摆并指出守恒量。
- 循环坐标与守恒律的关系?
更深一层:拉格朗日方法的"记账效率"账本
拉格朗日方法相对牛顿矢量法的优势可以精确量化。自由度匹配:单摆用牛顿法要列二维方程再处理绳子张力(未知约束力),拉格朗日方法一个坐标一步到位;$n$ 体系统有 $c$ 个约束时,牛顿法需要 $3n$ 个方程 + $c$ 个约束方程 + $c$ 个未知约束力,拉格朗日方法直接 $3n - c$ 个广义坐标方程——约束越多省得越多。约束力从不出现是方法的最大红利:理想约束的约束力在虚功中自动消失,而它们恰恰是牛顿法里最难算的部分。
广义坐标的选择自由还带来"按需定制":极坐标适合中心力、角度适合摆、弧长适合约束轨道——坐标选得好,方程少一半。这个选择自由在数学上由"广义坐标变换下拉格朗日方程形式不变"担保——协变性的第一次登场(广义相对论的坐标无关性是它的终极版)。
推导细读:d'Alembert 原理把动力学"化装"成静力学
达朗贝尔原理 $ \sum (\mathbf{F}_i - \dot{\mathbf{p}}_i)\cdot\delta\mathbf{r}_i = 0 $ 的巧思是把 $-\dot{\mathbf{p}}_i$ 当作"惯性力",于是动力学问题变成了"包含惯性力在内的一切力虚功为零"的静力学问题。虚功原理在静力学中说"理想约束下平衡的充要条件是主动力虚功为零"——把牛顿第二定律移项后套进同一个模子。之后的两步操作:选广义坐标展开虚位移 $\delta\mathbf{r}_i = \sum\frac{\partial \mathbf{r}_i}{\partial q_j}\delta q_j$,利用 $\delta q_j$ 的任意性让每个系数为零——拉格朗日方程的雏形。"化动为静"不是文字游戏:它把约束力从方程里合法地开除了,这是矢量力学做不到的。
在真实世界里:拉格朗日方法的工程报表
- 机器人学:$n$ 关节机械臂的动力学方程标配推导法——拉格朗日量逐关节写出,$\frac{d}{dt}\frac{\partial L}{\partial\dot q_j} - \frac{\partial L}{\partial q_j} = \tau_j$ 直接产出关节力矩方程(M 个关节耦合的质量矩阵 + 科氏/离心项);
- 车辆动力学:半车模型(车身俯仰 + 两个悬挂)用广义坐标列 $L$,悬挂设计参数扫描的基础方程;
- 航天器姿态:带内部动量轮的卫星姿态控制——循环坐标(轮转角)对应的广义动量守恒给出"轮转则体反转"的姿控原理。
三处的共同点:约束复杂、需要力矩/力的显式表达、参数要做扫描——正是拉格朗日方法的主场。
专题深入:广义坐标的"按需定制"哲学
广义坐标的选择自由值得展开成"按需定制"的方法论:极坐标适合中心力($r, \theta$ 天然分离径向与角向)、角度适合摆($\theta$ 直接是广义坐标)、弧长适合约束轨道(把二维问题降成一维)。坐标选得好,方程少一半——这个选择自由在数学上由"广义坐标变换下拉格朗日方程形式不变"担保(协变性的第一次登场)。协变性的深层含义:物理定律不依赖坐标选择——这是广义相对论"广义协变性"的力学预演。
约束力从不出现是方法的最大红利:理想约束的约束力在虚功中自动消失,而它们恰恰是牛顿法里最难算的部分。$n$ 体系统有 $c$ 个约束时,牛顿法需要 $3n$ 个方程 + $c$ 个约束方程 + $c$ 个未知约束力,拉格朗日方法直接 $3n - c$ 个广义坐标方程——约束越多省得越多。
数学细节展开:d'Alembert 原理的"化动为静"完整版
达朗贝尔原理 $\sum (\mathbf{F}_i - \dot{\mathbf{p}}_i)\cdot\delta\mathbf{r}_i = 0$ 的巧思值得展开:把 $-\dot{\mathbf{p}}_i$ 当作"惯性力",于是动力学问题变成了"包含惯性力在内的一切力虚功为零"的静力学问题。虚功原理在静力学中说"理想约束下平衡的充要条件是主动力虚功为零"——把牛顿第二定律移项后套进同一个模子。"化动为静"不是文字游戏:它把约束力从方程里合法地开除了,这是矢量力学做不到的。
之后的两步操作:选广义坐标展开虚位移 $\delta\mathbf{r}_i = \sum\frac{\partial\mathbf{r}_i}{\partial q_j}\delta q_j$,利用 $\delta q_j$ 的任意性让每个系数为零——拉格朗日方程的雏形。变分法($S = \int L\,dt$ 的极值条件)是它的"全局版"——真实轨道使作用量取极值,变分即得拉格朗日方程。
工程案例细节:拉格朗日方法的工程报表(工程版)
- 机器人学:$n$ 关节机械臂的动力学方程标配推导法——拉格朗日量逐关节写出,$\frac{d}{dt}\frac{\partial L}{\partial\dot q_j} - \frac{\partial L}{\partial q_j} = \tau_j$ 直接产出关节力矩方程(质量矩阵 + 科氏/离心项);
- 车辆动力学:半车模型(车身俯仰 + 两个悬挂)用广义坐标列 $L$,悬挂设计参数扫描的基础方程;
- 航天器姿态:带内部动量轮的卫星姿态控制——循环坐标(轮转角)对应的广义动量守恒给出"轮转则体反转"的姿控原理。
三处的共同点:约束复杂、需要力矩/力的显式表达、参数要做扫描——正是拉格朗日方法的主场。
例题
用拉格朗日方程求单摆(摆长 $l$、质量 $m$)的运动方程,并指出循环坐标(若有)与守恒量。
解答
广义坐标取 $\theta$:
小角度 $\sin\theta \approx \theta$:$\ddot\theta + \omega^2\theta = 0$,$\omega = \sqrt{g/l}$。
- $L$ 不显含 $t$ ⟹ 能量 $E = \frac12 ml^2\dot\theta^2 + mgl(1-\cos\theta)$ 守恒(这就是循环坐标思想的"时间版")。