完整约束保守系的拉格朗日方程.pdf

四、完整约束保守系的拉格朗日方程：上次课我们导出了在完整约束下的第二类拉格朗日方程：),2, 1(sQqTqTdtd,并用它解了一些题目考虑到如果我们要研究的系统所在的内外力场均是保守力场，或者其它作用于系统的力均不作虚功在这种情况下， 上面这条完整约束下的第二类拉格朗日方程还可以进一步简化这次课准备要讲的内容就是，先由这条拉格朗日方程推出完整约束下保守系的拉格朗日方程，并举例应用， 然后再讨论完整约束保守系的拉氏方程的一次积分我们由前面学过的知识可以知道，如果系统处在保守力场中，保守力系必有与其对应的势能V ，此势能是系统中各个质点的位置函数，即：V=V(nrrr21)，且有VFi1它的三个分量表达式为：iiz iiy iixzVFyVFxVF，，如果将ir用广义坐标表示：),(tqrrii则势能也就是广义坐标及时间t 的函数： V=V(q ，t)，由此我们很容易求得在保守力场中广义力Q的表达式由广义力的定义得：iiiiiiii i ii iqzzVyyyVqxxVqrVqrFQ[ 根据复合函数的微分规则可知其结果为 qV]将此结果代入第二类拉格朗日方程就可将它写成为：0)(qVTqTdtdqVqTqTdtd∵0 qV∴左边的式子又可写成为：0qVTqVTdtd在这里就定义：VTL， L 称作为拉格朗日函数，简称为拉氏 函 数 ， 它 就 等 于 系 统 的 动 能 与 势 能 之 差 。

那 么 上 式 就 可 写 成 为 ：sqLqLdtd,2, 10这个方程就是完整约束保守系的拉格朗日方程有时也叫它为拉格朗日方程或拉氏方程由前面的推导可知这个方程适用的条件是：完整约束， 保守力系或者除了保守力系之外的其它力均不作虚功，T 和 V 即 L 都是相对惯性系的量下面我们就举个例子用它来求解应用保守系拉氏方程解题的步骤基本上和用第二类拉氏方程解题的步骤相同，只是将上次课中的求广义力这一步改成为求势能V 及拉氏函数L 就可在这里要注意， 势能必须包含系统内力的势能和外力的势能还有由于保守系的拉氏方程中包含有 L 对 ,的偏导数，所以拉氏函数L 也一定要表达成广义坐标q和广义速度q的函数：tLL,.例 1：有两个质点，质量分别为1m和2m中间连一弹簧，其弹性系数为k，在光滑水平面上作直线运动求两质点的相对运动规律解：∵题目要我们求的是两质点的相对运动规律，那么，如果我们取m1、2m和这根相连的弹簧为研究系统的话， 整个系统在运动过程中所受的外力只有垂直水平面方向上的力，而这些垂直方向上的力对系统是不作功的，系统的内力是弹性力，我们知道弹性力是保守力，因此内力是有势的，所以该系统是可以应用保守系的拉格朗日方程的。

根据这样一分析，那么在这个问题中我们就取m1，m2和弹簧组成的整个系统为研究对象那么对于我们所取的这个系统应该有几个自由度？由于整个系统被约束在光滑水平面上作直线运动，所以系统中每一个质点只有一个自由度， 而弹簧是不能看作一个约束的，因此可知系统共有两个自由度有两个自由度，就得取两个广义坐标，这里的两个广义坐标应该怎么取呢？如图所示，先取任一固定点o 为坐标原点，并沿系统运动直线向右的方向为x 轴的正方向在此坐标系中两个质点的位置坐标分别为1x和2x由于在这个问题中讨论的是两质点组成的质点组的相对运动规律，所以我们在这里不取1x和2x为广义坐标， 而取质心坐标xc和两质点的相对距离x 分别为两个广义坐标，q1=xc，q2=x从图上分明看出：x=2x-1x研究系统的质心位置坐标：xc=212211 mmxmxm广义坐标确定了，那么我们就可以着手求出系统的动能、势能及拉格朗日函数系统的动能：T= 21m1x2 1+21m2x2 2，还需将它化成广义坐标xc,x 和广义速度xc, x的函数，为此我们将对上面两式求导将上面两个式子对时间t 求一次导数可得：12xxx212211 mmxmxmxc由 这 两 个 方 程 可 以 解 得 ：cxxmmmx212 1cxxmmmx211 2将这两个式子代入动能表达式中去：2 222 112121xmxmT22 212 212211 22212 1212121x mmmmxxmmmmxxmmmmcc2 21221212 212 2 212 12)(2121)(2121 ccxmmxmmmmxmmx mmmm这就得到了用广义速度表示的动能表达式，刚才已经分析过， 系统的内力即弹性力是有势的，其势能就是弹性势能V，故有系统的总势能：clxkV2)(21，c 为常数，l为弹簧的自然长度。

系统的拉氏函数：22 2122121)(21)(21 21lxkxmmxmmmmVTLc于是我们对应广义坐标cxq1就可列出一条保守系的拉氏方程为：0ccxLxLdtd，将得到的拉氏函数代进去则有：0)(21cxmmdtd，0cx质心加速度等于零同样道理 对 广 义 坐 标xq2的 拉 氏 方 程 为 ：0 xLxLdtd，拉 氏 函 数 代 进 去 求 偏 导 得0)()(0)(0)]([)(212121212121ymmmmklxlxkxmmmmlxkxmmmmdtd再令21212)( mmmmk则02yy这正是一个简谐振动所满足的微分方程，很容易得到它的解为：)sin(tAy即：ltAxtAlx)sin()sin(，其中的A 和是两个积分常数由初始条件决定，而就等于：2121)(mmmmk由最后所得到的结果可见两质点的相对运动是以角频率为的简谐振动到现在为止我们可以看出用拉格朗日方程求运动规律是很方便的，而且应用拉格朗日方程解题的步骤总是相同的最后我再强调一下，应用拉氏方程解题必须注意的几点：1.所选的广义坐标的数目必须与自由度相等2.系统的动能T 必须按惯性系计算3.计算广义力或势能时，每一质点所受的力iF必须计算一切内外力，但光滑约束的约束力及钢体的内力可不计算，换句话说，一切不做功的力都可以不用计算，4. 动能、广义力及势能均要换成t、、的函数才能代入到拉氏方程中。

要注意求偏导与全导数的区别，不要搞错了§5 完整保守系的拉氏方程的一次积分接下去讲完整约束保守系的拉氏方程的第一次积分由前面的讨论知道完整约束保守系的拉格朗日方程为：0 qLqLdtd（s21、）它是 S 个二阶常微分方程组，在某些特殊情况下可以得到它的第一次积分，通过一次积分可以使此运动微分方程降阶常见的拉氏方程的第一次积分有两种，一种是循环积分，另一种是能量积分下面我们先介绍循环积分一、循环积分：在怎样的特殊情况下才有循环积分呢？只有在拉氏函数中存在循环坐标时，才存在循环积分，所以在讨论循环积分之前得先介绍循环坐标的概念1、循环坐标 ：如果在拉氏函数中不显含广义速度对应的某个广义坐标，比如是iq，那么这 个 坐 标 就 被 称 为 循 环 坐 标 或 者 称 可 遗 坐 标 例 如 有 一 个 质 点 的 拉 氏 函 数 为L=21m(2x+2y) – mgy ，在此拉氏函数中不存在与x对应的广义坐标x，这个坐标x就是一个循环坐标 对于循环坐标iq来说，则有iqL=0， 此时拉氏函数就变成为：iqLdtd=0 于是得到拉氏方程的一次积分：iqL=ib(常量 ) 此积分就称它为2、循环积分：由此可见对于任一循环坐标就有一对应的循环积分。

在分析力学中定义系统的动能T对广义速度iq的偏导数为广义动量，通常用符号ip表示广义动量， 即ip=iqT 一般来讲势能V 与广义速度iq无关，所以拉氏函数L=T-V 对iq的偏导数也就是广义动量iP因此，上式也可写成为：iP（q,q,t）=ib这正好说明广义动量是守恒的，所以循环积分又称为广义动量积分如果在拉氏函数L 中出现 m 个循环坐标时，就可相应地得到m 个一次积分例如：在质点的抛射运动中，(不考虑空气阻力)用直角坐标表示的拉氏函数为：L= 21m(2x+2y+2z) –mgz ，显而易见，在此拉氏函数中出现几个循环坐标？有两个循环坐标L 中不含x 和 y所 以x 、 y是 循 环 坐 标 ， 在 此 对 应 的 两 个 循 环 积 分 为 ：xp=xL=mx=const, yp=yL=my=const, ，此结果表明抛射体在水平面内两个方向上的动量是守恒的又例如：质点在有心力场中运动时，在极坐标下的拉氏函数是：L= 21m(2r+2r2)+rm k2 ，显然在拉氏函数中不含有，极角 θ 是循环坐标与此对应的一个循环积分p=L=m2r= 常数，这就是我们早已熟悉的有心力场中的动量矩守恒式。

下面我们来讨论在什么条件下可由拉氏方程推出能量积分二、能量积分：也就是T+V= 常数1、积分的推导：能量积分的推导条件有三个：①完整约束，②保守力场，③L=L （q,q）拉氏函数不显含t 由①②两个条件可以用 qLdtd— qL=0 这条拉氏方程然而由第三个条件可得拉氏函数对时间t 的 全导数： dtdL=LL如果拉氏函数显含时间t 的话，这里还要多一项 tL现在根据第三个条件由于L 不是 t 的显含函数， 所以就没有 tL这一项了 由于在①②两个条件下可用拉方程： qLdtd--qL=0 由这个拉氏方程可得： qL= qLdtd将此代到前面的式子中去则有：dtdL= dtdL）（ qLdtdq= dtd（L）= dtd（q qL）将它移到等式的左边，则有： dtd(L— L)=0 ，显然对它进行一次积分可得：L—L =常量，在分析力学中就定义左边的量为哈密顿函数，并用符号H 表示，即HL—L 可见由上面三个条件利用拉氏方程推得的一次积分是H=E ，而上式的H 显然是 q 和q的函数： H=H(q,q) ，那么这个H 包含什么物理意义？这当然是我们所要关心的问题由前面的推导知道H=E 是在所给的三个条件下得到的。

对于在什么情况下才能有①②两个条件这是清楚的，至于在什么条件下才有L=L(q,q) 这个条件，现在我们还不知道的因此，我们得先来讨论在什么条件下拉氏函数L=L(q,q) 不显含 t2、什么条件下L=L(q,q)由第一节对约束概念的讨论知道：对一般的完整不稳定约束：f(1r⋯⋯nr,t)= 0 来说，质点的位矢ir不仅是广义坐标而且还是时间的函数，即ir=ir(q, t ) 这是不稳定的情况由于在不稳定约束的情况下ir=ir(q, t ) 在这种不稳定情况下位矢对时间的一阶导数：ir=qriq+ tri则系统的总动能T=i21 im2 ir=i21 im（ri+ tri）2 = 21s1s1ni 1imqriqriq q+s1ni 1imqritriq+21ni 1im（tri ）2由此可见按广义速度q来分，不稳定约束的情况下系统的动能T 是广义速度q的非齐次的二次函数我们可以将它简写成为：T=T2+T1+T0其中 T2是广义速度q的二次函数， T1是广义速度q的一次函数，T0则是广义速度的零次函数这是在非稳定情况下得到的结果 如果系统是处在稳定约束的情况下，结果又会怎么样呢？如果是稳定约束在稳定情况下， 系统中质点的位矢ir只是广义坐标的函数：ir=ir（q） 。

不显含时间t，此时，系统的动能： T= T2它仅仅是q广义速度的二次齐次函数，而且所选的广义坐标与固定坐标是没有相对运动的，显然在此情况下系统的动能是广义坐标q 和广义速度的函数，而不显合时间 t 即 T=T(q,q)因为 L=T-V ，而 V 一般是位置ir的函数， 所以在稳定的情况下ir=ir(q) 势能 V 也只是广义坐标的函数：)(qVV，也不是显含时间t，故有：),(LL——不显含时间t， 通过前面的讨论说明了在稳定约束条件下拉氏函数不显含时间t， 有 L=L(q, q)这样的形式，下面补充一些关于齐次函数的概念和欧拉齐次函数定理，即：3．欧拉齐次函数定理：(1)齐次函数的定义：设：mxxxfy2, 1是mxxx2, 1这 m 个变量的n 次齐次。