

是任意的(除了端点条件),积分要为零的唯一方式是被积函数为零。此外,
选择上的自由度使我们能够推断出被积函数中乘以
的因子必须恒为零,从而推导出拉格朗日方程。
现在
的选择并非完全自由。我们仍然可以从
的任意性推断出被积函数必须为零,88 但我们不再能推断出乘以
的因子为零(只能推断该因子在可接受变化上的投影为零)。因此我们有

受到约束。
路径 q 满足约束,如果
[q] =
o
[q] = 0。约束甚至在变化后的路径上也必须满足,因此我们只允许那些约束的变化为零的变化
:

与约束曲面
相切的条件为



?拉格朗日方程的残差正交89 于任何与约束曲面法线正交的
。一个与所有与给定矢量正交的矢量都正交的矢量平行于该给定矢量。因此,拉格朗日方程的残差平行于约束曲面的法线;两者必须成比例:

是时间的某个函数,可能依赖于所考察的路径。这些方程与约束方程
o
[q] = 0 一起构成控制方程。这些方程足以确定路径 q 并消去未知函数
。
视为一个坐标:

相关的拉格朗日方程就是约束方程。(注意
不出现在增广拉格朗日量中。)因此,这个增广拉格朗日量的拉格朗日方程完全封装了由添加显式坐标约束所施加的对拉格朗日方程的修正,代价是引入了额外的自由度。注意这个拉格朗日量的形式与我们在刚性系统推导 L = T - V 时使用的拉格朗日量(方程1.89)相同(第1.6.2节)。
,或者拉格朗日量(1.182)中引入的额外自由度纯粹是形式上的?
如果
可以写成解状态路径的函数,那么显然它由状态决定从而可以被消去。假设
可以写成一个状态依赖函数与路径的复合:
=
o
[q]。考虑拉格朗日量

的乘积产生的附加项。将欧拉-拉格朗日算符 E(参见第1.9节)应用于这个拉格朗日量得到90

[q] 复合得到路径 q 的拉格朗日方程。利用约束在路径上满足
o
[q] = 0 且因此 Dt
o
[q] = 0 的事实,我们得到

=
o
[q]。如果我们现在利用我们只处理坐标约束这一事实,
2
= 0,那么

=
o
[q] 由未增广的状态决定。这等价于说
可以被消去。
仅从增广拉格朗日量拉格朗日方程的形式有效性考虑,我们不能推断
可以写成一个状态依赖函数
与
[q] 的复合。从增广拉格朗日量导出的显式拉格朗日方程依赖于加速度 D2q 以及
,因此我们不能单独推断出两者中任一个是状态依赖函数与
[q] 的复合。然而,现在我们看到
就是这样一个复合。这使我们得以推断出 D2q 也是一个与路径复合的状态依赖函数。系统的演化由动力学状态决定。







,第二个方程乘以 sin
并相加,我们得到

得到的摆的拉格朗日方程相同。为完整起见,我们可以用其他变量求出
:

确实是状态函数与状态路径的复合。注意 2 l
是一个力——它是重力的向外分量与离心力之和。使用这种解释,在两个坐标运动方程中,涉及
的项是必须施加于无约束粒子以使其按照约束要求在圆上运动的力。等价地,我们可以将 2 l
视为摆杆中保持质量的张力。92
的弹簧和质量。我们可以分别为每个部分构造拉格朗日量。然后我们选择复合系统的拉格朗日量为两个分量拉格朗日量之和,并加上约束
= X1 + x1 来实现耦合。





是将系统维系在一起的约束力。我们现在可以消去“胶水”坐标
和
,得到坐标 x1 和 x2 中的运动方程:

练习1.34. 组合拉格朗日量 a. 制作另一个与本节中描述的弹簧-质量结构兼容的原始组件。例如,制作一个可以连接到弹簧-质量系统的摆。构建一个组合并推导运动方程。注意,如果坐标选择不当,代数运算将非常复杂。 b. 作为一个小项目,构建一族兼容的机械零件,每个零件由适当的拉格朗日量刻画,可以通过多种方式组合以构成有趣的机构。记住,在好的语言中,组合部分的结果本身应该是一个相同种类的部分,可以进一步与其他部分组合。
练习1.35. 三轴曲面上的珠子 再次考虑受约束在三轴曲面上运动的珠子(练习1.18)。使用直角坐标作为广义坐标,并附加以珠子保持在曲面上的显式约束,重新表述该问题。找到拉格朗日量并证明拉格朗日方程与练习1.18中找到的等价。
练习1.36. 微小高尔夫球的运动 将高尔夫球理想化为一个质点,在均匀引力场(加速度为 g)中,受约束在高度变化的 h(x, y) 无摩擦光滑表面上运动。 a. 为该系统寻找一个增广拉格朗日量,并推导支配质点在 x 和 y 中运动的方程。 b. 在什么条件下这可以近似为势函数 V(x, y) = mgh(x, y)? c. 假设 h(x, y) 关于 x = y = 0 轴对称。你能找到这样的 h 使得运动具有闭合轨道吗?
= 0。
是全时间导数意味着存在一个与速度无关的函数
使得

与速度无关意味着
2
= 0。作为状态函数,
与
之间的关系为

,我们可以通过求解这个线性偏微分方程找到
。解可确定到相差一个常数,因此
= 0 意味着
= K,其中 K 为某个常数。另一方面,如果我们知道
= K,则
= 0 随之成立。因此速度相关约束
= 0 等价于速度无关约束
= K,并且我们知道如何找到此类系统的拉格朗日方程。
如果 L 是无约束问题的拉格朗日量,则具有约束
= K 的拉格朗日方程为

是一个时间函数,将在求解过程中被消去。常数 K 不影响拉格朗日方程。函数
与速度无关,
2
= 0,因此拉格朗日方程变为


= 0 的拉格朗日方程为

写出拉格朗日方程,而无需构造积分
。但这些拉格朗日方程的有效性依赖于积分
的存在。
结果表明,增广拉格朗日量技巧在这里也适用。如果我们将拉格朗日量加上约束
乘以一个时间函数
',就得到这些拉格朗日方程:


= - D
',则与拉格朗日方程(1.212)相同。
有时,一个问题可以自然地用速度相关约束来表述。我们发展的形式体系可以处理任何可以写为坐标约束的导数的速度相关约束。这样的约束称为可积约束。任何约束可以置于坐标约束形式(或已经在该形式中)的系统称为完整系统。
练习1.37. 证明增广拉格朗日量(1.213)确实导致拉格朗日方程(1.214),考虑到
是
的全时间导数。

,圆环上任意点从任意参考方向转过的角度;以及 x,沿斜面向下的线性位移。约束是圆环不滑动。因此
的变化精确反映在 x 的变化中;约束函数为

+ c。我们可以使用积分后的约束或其导数来构造增广拉格朗日量。
动能有两部分:圆环的转动能量和其质心运动的能量。94 圆环的势能随高度降低而减少。因此我们可以写出增广拉格朗日量:




可以解释为实施约束所涉及的摩擦力。约束的摩擦力为


具有形式

不是一个全时间导数。如果 L 是无约束系统的拉格朗日量,则运动方程被断言为

= 0,系统封闭,系统的演化被确定。注意这些方程与
是全时间导数情况下的拉格朗日方程(1.212)相同,但这里这些方程的推导不再有效。
对于坐标约束
= 0(其中
2
= 0)的拉格朗日方程推导的关键步骤是注意必须满足两个条件:


[q] 与
正交,而
被约束为与
1
o
[q] 正交,所以两者在每一时刻必须平行:

与速度相关约束函数
相容,它必须满足(见方程1.179)

,因为
不再与
1
o
[q] 正交,而且由于
按假设不可积,我们不能将约束改写为坐标约束。
以下是 Arnold 等人 [6] 中对非完整方程的推导,转换为我们所用的记号。定义一个“虚速度”
为任何满足下式的速度


成立。因为
除了要求与
2
o
[q] 正交外是任意的,并且任何这样的
与 E[L] o
[q] 正交,所以
2
o
[q] 必须平行于 E[L] o
[q]。因此

替换为
。然而,为从方程(1.229)得到方程(1.230),我们必须令
= 0 并将 D
替换为
。所有非完整方程的“推导”都有类似的等同处理。归根结底:非完整方程并非来自作用量原理。它们是别的东西。它们是否正确取决于它们是否与实验相符。
对于坐标约束或导数约束系统,我们发现拉格朗日方程可以从用约束增广的拉格朗日量推导出来。然而,如果约束不可积,增广拉格朗日量的拉格朗日方程与非完整系统的方程(1.225)不同。96 令 L' 为带有不可积约束
的增广拉格朗日量:


相关的拉格朗日方程就是约束方程

和 D
。因此,通常的状态变量 q 和 Dq 连同约束一起,不足以确定推导出的拉格朗日方程的一组完整初始条件;我们还需要指定
的初始值。
一般而言,对于任何特定的物理系统,方程(1.225)和(1.234)并不相同,实际上它们具有不同的解。两组方程中哪一组能准确建模物理系统并不显然。第一种处理非完整系统的方法不能通过推广完整情形的论证来证明其合理性,而另一种方法则不是完全确定的。这可能表明这些模型是不充分的,需要指定更多关于约束如何维持的细节。

91 这个约束与第1.6.2节中证明 L = T - V 可用于刚性系统时使用的约束具有相同的形式。这里它是一个更一般约束集的具体例子。 92 事实上,如果像在牛顿约束力的讨论中那样对约束方程进行缩放,我们本可以将
等同于约束力 F 的大小。然而,尽管
通常与约束力有关,但它并非其中之一。我们选择保持缩放的天然形式,而不是为了使结果人为地好看而进行改动。
93 这个例子出现在 [20] 中,第49-51页。
94 我们将在第2章看到如何计算转动动能,但目前答案是 (1/2) M R2
2。
95 关于非完整系统的一些处理方法,参见例如 Whittaker [46]、Goldstein [20]、Gantmakher [19] 或 Arnold 等人 [6]。
96 Arnold 等人 [6] 将约束添加到拉格朗日量中的变分力学称为 Vakonomic 力学。