拉格朗日方程的应用及举例08讲.doc

拉格朗日方程的应用及举例
拉格朗日方程有以下几个特点：（1）拉格朗日方程适用于完整系统，可以获得数目最少的运动微分方程，即可以建立与自由度数目相同的n个方程，是一个包含n个二阶常微分方程组，方程组的阶数为2n。求解这个方程组可得到以广义坐标描述的系统运动方程。（2）拉格朗日方程的形式具有不变性。对于任意坐标具有统一的形式，即不随坐标的选取而变化。特别是解题时有径直的程序可循，应用方便。（3）所有的理想约束的约束反力均不出现在运动微分方程中。系统的约束条件愈多，这个特点带来的便利越突出。（4）拉格朗日方程是以能量的观点建立起来的方程，只含有表征系统运动的动能和表征主动力作用的广义力，避开了力、速度、加速度等矢量的复杂运算。（5）拉格朗日方程不但可以建立相对惯性系的运动，还可以直接建立相对非惯性系的动力学方程，只要写出的动能是绝对运动的动能即可，至于方程所描述的运动是对什么参考系的运动，则取决于所选的广义坐标。
纵观拉格朗日方程，看出分析力学在牛顿力学的基础上，提出严密的分析方法，从描述系统的位形到建立微分方程都带有新的飞跃。我们还应看到，虽然拉格朗日方法在理论上和应用上都有重要的价值，但是，牛顿力学的价值并未降低，特别是它的几何直观性和规格化的方法使人乐于应用，由于计算机的广泛使用，牛顿一欧拉方法又有所发展。我们将会看到，用拉格朗日方程求解，在获得数量最少的运动微分方程时，其求导过程有时过于繁琐，并有较多的耦合项。
应用拉格朗日方程建立动力学方程时，应首先建立以广义坐标q和广义速度表示的动能函数和广义力Q。为此，首先讨论动能的计算和广义力的计算，在此基础上，再讨论拉格朗日方程的应用。
一、动能的计算
对于系统的动能，可以写出关于广义速度的齐次函数的表达式。在实际计算中，应用理论力学的有关知识就可以建立以广义坐标和广义速度所表达的动能函数。
例1-1 已知质量为m，半径为r的均质圆盘D，沿OAB直角曲杆的AB段只滚不滑。圆盘的盘面和曲杆均放置在水平面上。已知曲杆以匀角速度w1绕通过O点的铅直轴转动，试求圆盘的动能。
解：取广义坐标x和j，x为圆盘与曲杆接触点到曲杆A点的距离，j为曲杆OAB的转角，j = w1t。
应用柯尼希定理求圆盘的动能。为此，先求圆盘质心C的速度和相对于质心平动坐标系的角速度。若以曲杆OAB为动参考系，C为动点，

再应用刚体绕二平行轴转动的合成方法，圆盘的角速度为

于是圆盘的动能为

若将动能表达式展开，得到

可以看出，圆盘的动能包含广义速度的二次项，广义速度的一次项和它的零次项。

二、广义力的计算
概括地说，广义力有三种计算方法：
1）根据广义力的定义，有

我们可以按照这个公式来计算，但是，有时计算是繁冗的。
2）我们知道，作用在系统上的诸主动力对于任何虚位移元功之和等于诸广义力对于相应的广义坐标的虚位移元功之和，即

对于完整系统，广义坐标的变分dq1，dq2，…，dqn是彼此独立的。若给出某一广义坐标的变分为dqj，而令其它坐标变分均为零，即
dqj≠0，dq1 = dq2 = … = dqj－1 = dqj＋1 = … = dqn = 0
则上式为

于是
由于系统的主动力在给定的虚位移中元功之和的计算是我们熟悉的，则广义力Qj可较易地计算出。依次给出不同序数的坐标变分的同时，令其它坐标变分为零，则可依次计算出与广义坐标对应的广义力。这种方法是我们经常应用的。
3）若作用于系统上的主动力有势，则通过势能函数即可求出广义力。设势能函数为V，则可应用式

进行广义力的计算。
例1-3 均质杆OA和AB在A点铰链连接，并在O点用铰链支承。杆重分别为P1和P2，F1为作用于B点的水平力，试求对应于j和y的广义力。
解：系统具有两个自由度。依题意，取j和y为广义坐标，对应于j和y的广义力以Qj和Qy表示。于是，

当j获得变分dj，而y保持不变，即dy = 0时，

当y获得变分dy，而dj = 0时，

三、拉格朗日方程的应用
应用拉格朗日方程建立系统的动力学方程时，一般采用以下步骤：
1）分析系统的约束条件，判断系统的类型是否为完整系统，是定常还是非定常的，是保守的还是非保守的。
2）若系统为完整的，在确定其自由度数目后，选择恰当的广义坐标。
3）计算出以广义速度表达的动能T(q，，t)、势能V (q，t) 或广义力Q (q，t)，若主动力有势，计算出拉格朗日函数L(q，，t)。
4）列出拉格朗日方程。
例1-4 半径为R、质量为m的圆环挂在一半径为r的固定圆柱上。设圆环与圆柱间有足够大的摩擦力阻止相对滑动，试写出圆环的运动微分方程，并求微幅摆动的周期。
解：圆环具有一个自由度，是完整系统。取q为广义坐标，圆环的动能为

其中，瞬心为A，则

于是
主动力有势，系统的势能为
V =－mg (R－r) cosq

代入拉格朗日方程，得到系统的动力学方程：

即
考虑到微幅，有

周期为
由于主动力有势，可以写出拉格朗日函数：

代入式（1-25）中同样可以得到系统的动力学方程。
2. 已知摆线绕在固定圆柱上，尺寸如图；求此摆的运动微分方程。
解 这是单自由度保守系统，选q为广义坐标，选q = 0为系统的零势能位置，则

将T、V代入保守系统的拉格朗日方程

或将拉格朗日函数L = T - V代入如下形式的拉格朗日方程

皆可得运动微分方程

3. 已知三均质齿轮，半径皆为r，质量都是m，此机构位于水平面内，若无重系杆受矩为M的力偶作用；求系杆的角加速度a。

解 这是单自由度非保守系统，选系杆的转角j为广义坐标，则有关的角速度和速度为

该系统的广义力为 Qj = M
动能为
代入拉格朗日方程
得

例1-9 试求例1-1中圆盘的运动微分方程。又，若t = 0时，x = 10cm，= 0，求当x=20cm时，为多少？
例1-1 已知质量为m，半径为r的均质圆盘D，沿OAB直角曲杆的AB段只滚不滑。圆盘的盘面和曲杆均放置在水平面上。已知曲杆以匀角速度w1绕通过O点的铅直轴转动，试求圆盘的动能。
解：由例1-1已求得动能T为

水平台为零势面，则圆盘的势能为
V = 0
系统的拉格朗日函数L为

代入拉格朗日方程，有

由于系统是非定常的，虽然作用于圆盘上的主动力有势，但并不存在能量积分，由于拉格朗日函数L不显含时间t，系统有广义能量积分。由动能表达式得到

圆盘的广义能量积分为
T2－T0 + V＝常数.
于是得到

整理后，有

当t = 0时，x0 = 10cm，= 0，则

于是有
当x = 20cm时，
cm/s
例9 质量为m，半径为r的圆环O竖立在一粗糙平面上。圆环的边缘上刚连一质量为m的质点A。试写出系统的运动微分方程。

解：由圆环O和质点A组成的系统只能在地面上作纯滚动，自由度为1，取OA与铅垂线的夹角为广义坐标，以系统为研究对象 O点处水平面为零势能面，则系统的动能和势能分别为

于是有

代入拉格朗日方程，导出

例1-7 三角楔块A可沿水平光滑面作直线运动，楔块A的质量为m1，其上受有简谐力F＝Hsinwt的作用（H和w均为常量）。楔块斜边BD上有一质量为m2、半径为r的圆柱体，沿BD滚动而不滑动，二弹簧的刚体系数分别为k1和k2。试建立系统的运动微分方程。
解：系统具有二个自由度。取三角楔块的位移x和圆柱体相对于楔块的位移x为广义坐标，二者均以其静平衡位置为原点。
楔块A作平动，，圆柱体作平面运动，质心速度vC为

角速度w为

系统的动能T为

系统的势能V为

在平衡位置有关系式

于是势能V为

非有势力F相应的广义力分别为

又，

代入拉格朗日方程，得到系统的运动微分方程：

9