第一章理论力学第二册

1、第一章第一章分析力学基础分析力学基础 18 18世纪提出了处理多个约束的刚体系统动力学问题。世纪提出了处理多个约束的刚体系统动力学问题。 利用利用矢量力学矢量力学分析出现以下问题：分析出现以下问题： 对于复杂约束系统约束力的性质和分布是未知的；对于复杂约束系统约束力的性质和分布是未知的； 表述形式复杂。如球坐标系下的运动方程。表述形式复杂。如球坐标系下的运动方程。 质点系问题为大量方程的微分方程组。质点系问题为大量方程的微分方程组。 1788 1788年拉格朗日发表了年拉格朗日发表了分析力学分析力学一书，提出了一书，提出了解决动力学问题的新观点和新方法：解决动力学问题的新观点和新方法：采用功和

2、能量来描采用功和能量来描述物体的运动和相互作用力之间的关系。述物体的运动和相互作用力之间的关系。与矢量力学相比，与矢量力学相比，分析力学的特点：分析力学的特点： （3 3）追求一般理论和一般模型，对于具体问题，只要代入和展开）追求一般理论和一般模型，对于具体问题，只要代入和展开 的工作，处理问题规范化。的工作，处理问题规范化。（1 1）把约束看成对系统位置（速度）的限定，而不是看成一种力。）把约束看成对系统位置（速度）的限定，而不是看成一种力。 （2 2）使用广义坐标、功、能等标量研究系统运动，大量使用数学）使用广义坐标、功、能等标量研究系统运动，大量使用数学 分析方法，得到标量方程。分析方法

3、，得到标量方程。（4 4）不仅研究获得运动微分方程的方法，也研究其求解的一般方）不仅研究获得运动微分方程的方法，也研究其求解的一般方 法。法。 在完整约束的条件下，确定质点系位置的独立参数在完整约束的条件下，确定质点系位置的独立参数的数目，称为质点系的的数目，称为质点系的自由度数自由度数，简称，简称自由度自由度。 1 11 1 自由度和广义坐标自由度和广义坐标 例：确定一个质点在空间的位置需例：确定一个质点在空间的位置需3 3个独立的参量个独立的参量 自由质点为自由质点为3 3个自由度。个自由度。例：质点例：质点M 被限定只能在球面被限定只能在球面2222)()()(Rczbyax的上半部分运

4、动的上半部分运动由此解出由此解出222)()(byaxRcz这样该质点在空间中的位置就由这样该质点在空间中的位置就由x,y这两个独立参数所确定这两个独立参数所确定它的自由度数为它的自由度数为2。 n个质点组成的质点系个质点组成的质点系， 若受到若受到s个完整约束作用个完整约束作用自由度数为自由度数为 N=3n-s描述质点系在空间中的位置的独立参数称为描述质点系在空间中的位置的独立参数称为广义坐标。广义坐标。对于完整约束对于完整约束广义坐标的数目系统的自由度数广义坐标的数目系统的自由度数思考：思考：非完整约束，广义坐标数目和系统的自由度数目的关系？非完整约束，广义坐标数目和系统的自由度数目的关系

5、？12(,)0(1 2 3)knfrrrtks， ，12(3)NqqqNns， ，拉格朗日广义坐标拉格朗日广义坐标 )21()(21nitqqqrrNii， 约束方程为约束方程为 系统系统N个独立的坐标参量表示为个独立的坐标参量表示为系统的系统的n个坐标参量个坐标参量 设由设由n个质点组成的系统受个质点组成的系统受s个完整双侧约束个完整双侧约束 其中其中)21(Nkqk， 为广义坐标为广义坐标 的变分称为的变分称为广义虚位移广义虚位移。 kq1(12)Niikkkrrqinq， ， ，例：一单摆在空间摆动，摆长为例：一单摆在空间摆动，摆长为l。 Oxyz约束方程为约束方程为 2222( , ,

6、 )f x y zxyzl自由度数为自由度数为2 2。 x，y为独立变量为独立变量 ( , )xx x yx( , )yy x yy222( , )zz x ylxy, （单摆在（单摆在xy面上的投影与面上的投影与x轴夹角）为独立变量。轴夹角）为独立变量。 ( , )sincosxxl ( , )sinsinyyl ( , )coszzl 思考：思考：导弹在追踪飞机的情况下，广义坐标的数目和自由度导弹在追踪飞机的情况下，广义坐标的数目和自由度数目的关系如何？数目的关系如何？描述导弹的位置：描述导弹的位置：质心的位置质心的位置,CCxy导弹的纵轴和导弹的纵轴和x 轴的夹角轴的夹角独立的广义坐标数

7、目为独立的广义坐标数目为3 3 约束方程约束方程 导弹的速度方向要对准飞机的质心导弹的速度方向要对准飞机的质心 CPCCPCyyyxxx非完整约束非完整约束 CPCCPCyyyxxx独立的虚位移数目自由度数目独立的虚位移数目自由度数目2 2 设作用在第设作用在第i个质点上的主动力的合力个质点上的主动力的合力 iF在三个坐标轴上的投影分别为在三个坐标轴上的投影分别为 (,)xiyiziFFF，虚功方程虚功方程1111111()() 0nFFiinNNNiiixikyikzikikkkkkkNniiixiyizikkikkkWWxyzFqFqFqqqqxyzFFFqqqq 1 12 2 以广义坐标

8、表示的质点系平衡条件以广义坐标表示的质点系平衡条件 1.1.以广义坐标表示的质点系平衡条件以广义坐标表示的质点系平衡条件10NFkkkWQq称为与广义坐标称为与广义坐标 kq相对应的相对应的广义力广义力 。 由于广义坐标的独立性由于广义坐标的独立性 kq可以为任一值可以为任一值 021 NQQQ如令如令 )21()(1NkqzFqyFqxFQnikizikiyikixik， 质点系的平衡条件是系统所有的广义力都等于零。质点系的平衡条件是系统所有的广义力都等于零。 用广义坐标表示的质点系的平衡条件用广义坐标表示的质点系的平衡条件 求广义力的两种方法求广义力的两种方法 1.1.直接计算法（解析法）

9、直接计算法（解析法） 2.2.几何法几何法 令某一个令某一个 不等于零不等于零 kq而其他而其他N-1个广义虚位移都等于零个广义虚位移都等于零 FkkWQqFkkWQq)21()(1NkqzFqyFqxFQnikizikiyikixik， 利用广义虚位移的任意性，利用广义虚位移的任意性， 例例1 11 1已知：杆已知：杆OA和和AB以铰链相连，以铰链相连， O端悬挂于圆柱铰链上，端悬挂于圆柱铰链上， 杆长杆长OA=a AB=b，杆重和铰链的摩擦都忽略不计。杆重和铰链的摩擦都忽略不计。今在点今在点A和和B分别作用向下的铅锤力分别作用向下的铅锤力 和和AFBF又在点又在点B作用一水平力作用一水平力

10、F试求：平衡时试求：平衡时 与与 ， ， 之间的关系。之间的关系。 21，AFBFF系统有两个自由度。系统有两个自由度。 现选择现选择 和和 为系统的两个广义坐标为系统的两个广义坐标 12计算其对应的广义力计算其对应的广义力 和和1Q2Q用第一种方法计算广义力：用第一种方法计算广义力：22221111BBBAABBBAAxFyFyFQxFyFyFQ11212coscoscossinsinABByayabxab，)21()(1NkqzFqyFqxFQnikizikiyikixik， 解：解：故故111111cossinsinaxayayBBA，22222cossin0bxbyyBBA，系统平衡时

11、应有系统平衡时应有 0cossin0cossin)(222111FbbFQFaaFFQBBABBAFFFFF21tantan，用第二种方法计算：用第二种方法计算： 保持保持 不变，不变， 2只有只有 时时11111sincosABByyaxa 则对应于则对应于 的广义力为的广义力为 11111AABBBWFyFyF xQ可得一组虚位移可得一组虚位移 111cossin)(FaaFFQBA保持保持 不变，不变， 1只有只有 时时2可得另一组虚位移可得另一组虚位移 22220sincosABByybxb 对应于对应于 的广义力的广义力 2222222sincosAABBBBWFyFyF xQF b

12、Fb 120QQ例例 1 12 2已知：重物已知：重物A和和B分别连接在细绳两端，重物分别连接在细绳两端，重物A放置在粗放置在粗 糙的水平面上。重物糙的水平面上。重物B绕过定滑轮绕过定滑轮E铅直悬挂。铅直悬挂。 在动滑轮在动滑轮H的轴心上挂一重物的轴心上挂一重物C。 设重物设重物A重量为重量为P2重物重物B重量为重量为P，不计动滑轮不计动滑轮H的重量。的重量。试求：平衡时重物试求：平衡时重物C的重量的重量 ；CP以及重物以及重物A与水平面间的静滑动摩擦因数。与水平面间的静滑动摩擦因数。 系统具有两个自由度。系统具有两个自由度。,Ax广义坐标广义坐标： By首先令首先令 向右向右， Ax0By

13、主动力所做虚功的和为主动力所做虚功的和为 1()2AAACCACAWFxPyFPx 对应广义坐标对应广义坐标 的广义力为的广义力为 Ax12AxACAAWQPFx解：解：12ByBCBWQPPx 因为系统平衡时应有因为系统平衡时应有 0yBxAQQPPFPPCAC212 ，因此平衡时因此平衡时, ,要求物块与台面间静摩擦因数要求物块与台面间静摩擦因数 5 . 02PFfA再令再令 向下，向下， By0Ax 2BCyy ()2CFBCCBPWPyPyPy2.2.以广义坐标表示的保守系统的平衡条件及系统的稳定性以广义坐标表示的保守系统的平衡条件及系统的稳定性 如果作用在质点系上的主动力都是有势力，

14、势能为如果作用在质点系上的主动力都是有势力，势能为111()nnnVV xyzxyz， ， ， ， ， ，各力的投影为各力的投影为 iziiyiixizVFyVFxVF，( )( )FxiiyiiziiiiiiiiWFxFyFzVVVxyzVxyz 虚功为虚功为 虚位移原理的表达式成为虚位移原理的表达式成为 0V 在势力场中，具有理想约束的质点系的平衡条件为质点在势力场中，具有理想约束的质点系的平衡条件为质点 系的势能在平衡位置处一阶变分为零。系的势能在平衡位置处一阶变分为零。如果用广义坐标如果用广义坐标 表示质点系的位置。表示质点系的位置。 Nqqq，21则质点系的势能可以写成广义坐标的函数

15、则质点系的势能可以写成广义坐标的函数 )(21NqqqVV，()()(12)iiikxiyizikkkiiiikikikkxyzQFFFqqqxyzVVVxqyqzqVkNq ， ， ，由广义坐标表示的平衡条件可写成如下形式由广义坐标表示的平衡条件可写成如下形式 0(12)kkVQkNq ， ， ， 在势力场中具有理想约束的质点系的平衡条件是势能在势力场中具有理想约束的质点系的平衡条件是势能 对于每个广义坐标的偏导数分别等于零。对于每个广义坐标的偏导数分别等于零。不稳定平衡：不稳定平衡：在平衡位置上系统势能具有极大值。在平衡位置上系统势能具有极大值。随遇平衡：随遇平衡：系统在某位置附近其势能是

16、不变的。系统在某位置附近其势能是不变的。 稳定平衡：稳定平衡：在平衡位置处系统势能具有极小值。在平衡位置处系统势能具有极小值。对于一个自由度系统，对于一个自由度系统， 系统具有一个广义坐标系统具有一个广义坐标q，因此系统势能可以表示为因此系统势能可以表示为q的一元函数的一元函数即)(qVV 当系统平衡时，当系统平衡时，在平衡位置处有在平衡位置处有0ddqV如果系统处于稳定平衡状态，则在平衡位置处系统势能具如果系统处于稳定平衡状态，则在平衡位置处系统势能具有极小值。有极小值。即系统势能对广义坐标的二阶导数大于零即系统势能对广义坐标的二阶导数大于零 0dd22qV一个自由度系统平衡的稳定性判据一个

17、自由度系统平衡的稳定性判据例例 1 13 3已知：如图所示一倒置的摆，摆锤重量为已知：如图所示一倒置的摆，摆锤重量为 ，摆杆长度为，摆杆长度为l，在摆杆上的点在摆杆上的点A连有一刚度为连有一刚度为k的水平弹簧，摆在铅直的水平弹簧，摆在铅直位置时弹簧未变形。设位置时弹簧未变形。设OA=a摆杆重量不计。摆杆重量不计。P试求：摆杆的平衡位置及稳定平衡时所应满足的条件。试求：摆杆的平衡位置及稳定平衡时所应满足的条件。解：解：该系统是一个自由度系统，该系统是一个自由度系统， 选择摆角选择摆角 为广义坐标。为广义坐标。 摆的铅直位置为摆锤重力势能和弹簧弹性势能的零点。摆的铅直位置为摆锤重力势能和弹簧弹性势

18、能的零点。系统的总势能为系统的总势能为22222212sin221)cos1 (kaPlkaPlV由22sin22)(21PlkaV有有)(dd2PlkaV由由0ddV得到系统的平衡位置为得到系统的平衡位置为0kFPlkaV222dd对于稳定平衡对于稳定平衡要求要求0dd22V即即02 PlkakPla 1 13 3 动力学普遍方程动力学普遍方程n个质点组成的系统个质点组成的系统第第i个质点个质点 , ， ， ， 。imirir iFiFN惯性力为惯性力为iiirmF I理想约束作用理想约束作用0)()(11INiniiiiiniiiirrmFrFFF 在理想约束的条件下，质点系在任一瞬时所受

19、的主动力在理想约束的条件下，质点系在任一瞬时所受的主动力系和虚加的惯性力系在虚位移上所作的功的和等于零。系和虚加的惯性力系在虚位移上所作的功的和等于零。 写成解析表达式写成解析表达式1()()() )0nxiiiiyiiiiziiiiiFm xxFm yyFm zz动力学普遍方程动力学普遍方程特别适合于求解非自由质点系的动力学问题。特别适合于求解非自由质点系的动力学问题。例例 14已知：滑轮系统中，动滑轮上悬挂着质量为已知：滑轮系统中，动滑轮上悬挂着质量为 的的 重物，绳子绕过定滑轮后悬挂着质量为重物，绳子绕过定滑轮后悬挂着质量为 的重物。设滑轮和绳子的重量以及轮轴摩的重物。设滑轮和绳子的重量

20、以及轮轴摩 擦都忽略不计。擦都忽略不计。1m2m求：质量为求：质量为 的物体下降的加速度。的物体下降的加速度。2m1m2m解：解：取整个滑轮系统为研究对象。取整个滑轮系统为研究对象。1m g2m g1IF2IF1s2s222I111 IamFamF，由动力学普遍方程由动力学普遍方程222211 11()()0m gm asm gm as221122sasa，gmmmma12122424例例 1 15 5已知：两相同均质圆轮半径皆为已知：两相同均质圆轮半径皆为R，质量皆为质量皆为m。求：当细绳直线部分为铅垂时，轮求：当细绳直线部分为铅垂时，轮II中心中心C 的加速度。的加速度。解：解：研究整个系

21、统。研究整个系统。22121122IFmaMmRMmRI1I2，此系统具有两个自由度此系统具有两个自由度取转角取转角 为广义坐标为广义坐标21，令令1200，则点则点C 下降下降2hR动力学普遍方程动力学普遍方程120mg hF hMI2OC121IM2IMmgIF2102gaR（a）令令1200，则则1hR代入动力学普遍方程代入动力学普遍方程1110mg hFhMI或或1102gaR（b b）运动学关系运动学关系12aRR（c c）联立式（联立式（a a）（）（b b）（）（c c）解出）解出ga54 1 14 4 第二类拉格朗日方程第二类拉格朗日方程 设由设由n质点组成的系统受质点组成的系

22、统受s个完整约束作用，系统具有个完整约束作用，系统具有N=3n-s个自由度。个自由度。设设 为系统的一组广义坐标为系统的一组广义坐标Nqqq，211Niikkkrrqq11nNiikkikFrQ q11() () 0nniiiiii iiiiFFFrFmrrNI111() ()0nNniii iiki ikikikrFmrrQmrqq对于完整约束系统，其广义坐标是相互独立的。对于完整约束系统，其广义坐标是相互独立的。故故 是任意的，是任意的，(12)kq kN， ， ，为使上式恒成立，为使上式恒成立， 必须有必须有)21(01NkqrrmQnikiiik， 广义惯性力广义惯性力上式不便于直接应

23、用，上式不便于直接应用， 为此可作如下变换：为此可作如下变换：（1）kikiqrqr证明：证明：)21()(21nitqqqrrNii， 1ddsiiiikkkrrrrqtqt注意注意 和和 只是广义坐标和时间的函数只是广义坐标和时间的函数kiqrtri（2）kikiqrqrt)(dd证明证明：12()iiNkkrrqqqtqq， ， ，iikkrrqq)21()(21nitqqqrrNii， 对时间求微分对时间求微分kijNjkjikijkiNjjkiqtrqqqrqrtqqrqqrt2121)()()(dd而而tqrqqqrtrqqrqqrkijNjjkiijNjjikki2121)(若函

24、数若函数 的一阶和二阶偏导数连续的一阶和二阶偏导数连续)(21tqqqrrNii，kikiqrqrt)(dd111()()nnniiii iii iiiikkkrrrmrmrmrqtqtqdddd11nniiii iiikkrrmrmrtqqdd111()2nnii ii iiiikkrmrmr rtqqdd 221111()()22nniiiiiikkmmtqqdd()kkTTtqqdd得到得到)21(0)(ddNkQqTqTtkkk，第二类拉格朗日方程第二类拉格朗日方程 拉格朗日方程拉格朗日方程方程式的数目等于质点系的自由度数。方程式的数目等于质点系的自由度数。)21(01NkqrrmQn

25、ikiiik， 如果作用在质点系上的主动力都是有势力（保守力）如果作用在质点系上的主动力都是有势力（保守力）(12)kkVQkNq ， ， ，于是拉格朗日方程可以写成于是拉格朗日方程可以写成)21(0)(ddNkqVqTqTtkkk， 引入拉格朗日函数（又称为动势）引入拉格朗日函数（又称为动势）VTL则拉格朗日方程又可以写成则拉格朗日方程又可以写成)21(0)(ddNkqLqLtkk， 例例 1 16 6已知：轮已知：轮A沿水平面纯滚动，轮心以水平弹簧联于墙上。沿水平面纯滚动，轮心以水平弹簧联于墙上。A，B两轮皆为均质圆盘，两轮皆为均质圆盘，质量为质量为 的物块的物块C以细绳跨过定滑轮以细绳跨

26、过定滑轮B联于点联于点A1m半径为半径为R，质量为质量为2m弹簧刚度为弹簧刚度为k，质量不计。质量不计。试求：当弹簧较软，在细绳能始终保持张紧的条件下，试求：当弹簧较软，在细绳能始终保持张紧的条件下， 此系统的运动微分方程。此系统的运动微分方程。解解： 此系统具有一个自由度，此系统具有一个自由度， 以物块平衡位置为原点。以物块平衡位置为原点。取取x为广义坐标为广义坐标。以平衡位置为重力零势能点。以平衡位置为重力零势能点。取弹簧原长处为弹性力零势能点。取弹簧原长处为弹性力零势能点。系统在任意位置系统在任意位置x处的势能为处的势能为gxmxkV120)(21其中其中 为平衡位置处弹簧的伸长量为平衡

27、位置处弹簧的伸长量0此系统的动能为此系统的动能为222222122222111 111 1()()22 222 21()2xxTm xm Rm xm RRRmm x系统的动势为系统的动势为gxmxkxmmVTL120212)(21)21(代入拉格朗日方程代入拉格朗日方程0)(ddxLxLt得得0)2(1012gmkxkxmm注意到注意到gmk10则系统的运动微分方程为则系统的运动微分方程为0)2(12kxxmm 例例 1 17 7试求：此系统的运动微分方程。试求：此系统的运动微分方程。已知：运动系统中，已知：运动系统中，可沿光滑可沿光滑，两个物体两个物体重物重物 的质量为的质量为1M1m摆锤摆

28、锤 的质量为的质量为2M2m水平面移动。水平面移动。用无重杆连接，杆长为用无重杆连接，杆长为l。解：解：选选 和和 为广义坐标为广义坐标1x12120sincosyxxlyl（a a）将式（将式（a a）两端对时间求导数）两端对时间求导数sincos02121lylxxy，（b b）系统的动能系统的动能)cos2(2)(21)(212112221212222221xllmxmmyxmxmT)cos1 (2glmV则系统的势能为则系统的势能为选质点选质点 在最低处时的位置为系统的零势能位置在最低处时的位置为系统的零势能位置2M由此得由此得2212111cos)(0lmxmmxTxT，222121

29、1sincos)()(dd lmlmxmmxTt01xVQxcossin122212x lmlmTxlmT ，)sincos()(dd112 xxllmTtsin2glmVQ把以上结果代入拉格朗日方程中把以上结果代入拉格朗日方程中0sincos)(222121 lmlmxmmsin)sincos(2112glmxxllm 如果质点如果质点 摆动很小，摆动很小，2M可以近似地认为可以近似地认为1cossin，且可以忽略含且可以忽略含 和和 的高阶小量，的高阶小量，21x上式可改写为上式可改写为0)(2121 lmxmm（c c）gxl1 （d d）从以上两式中消去从以上两式中消去1x ，得到，得

30、到0121lgmmm （e e）这是自由振动的微分方程，其解为这是自由振动的微分方程，其解为)sin(0tA（f f）固有角频率为固有角频率为lgmmm1210摆动周期摆动周期101222mlTmm g（g）如果如果21mm 则质点则质点 的位移的位移 将很小将很小1M1x质点质点 的摆动周期将趋于普通单摆的周期的摆动周期将趋于普通单摆的周期2M1lim2mlTg若将式（若将式（e e）代入（）代入（d d） 得到得到gmmx121 （h）可见质点可见质点 沿沿x方向也作自由振动。方向也作自由振动。1M将式（将式（f f）代入，）代入，)21(0)(ddNkQqTqTtkkk， 对于保守系统，

31、在一定条件下，可以直接给出初积分的对于保守系统，在一定条件下，可以直接给出初积分的一般形式。一般形式。1.1.能量积分能量积分若系统所受到的约束均为定常约束，则若系统所受到的约束均为定常约束，则1Niiikikrvrqq12() (12)iiNrr qqqin， ， ， 1 15 5 拉格朗日方程的初积分拉格朗日方程的初积分1111111()() )2212nnNNiiiiiikliiklklNklklklrrTmv vmqqqqm q q ，其中其中likiniiklqrqrmm1广义质量广义质量TqqTkNkk21关于齐次函数的欧拉定理关于齐次函数的欧拉定理注意势能注意势能V不含不含 项项

32、iq 从而从而TqqTqqLkNkkkNkk211)21(0)(ddNkqLqLtkk， 0)2(dddddd2)()(dd)(dd)(dd1111LTttLtTqqLqqLqqLtqqLqqLqqLtqqLqqLtkkNkkkNkkkNkkkkkkkNkkkkk 2T-L=T+V=常数常数-保守系统的机械能守恒定律保守系统的机械能守恒定律-保守系统中拉格朗日方程的能量积分保守系统中拉格朗日方程的能量积分2.2.循环积分循环积分 如果拉格朗日函数如果拉格朗日函数L中不显含某广义坐标中不显含某广义坐标 ，则称该则称该坐标为坐标为循环坐标循环坐标 kq0)(dd0kkqLtqL，kqL常数常数-拉

33、格朗日方程的循环积分拉格朗日方程的循环积分如果引入广义动量如果引入广义动量kkqLpkp常数常数-广义动量守恒广义动量守恒例例 1 18 8 如图所示一个均质圆柱体，可绕其垂直中心轴自由转如图所示一个均质圆柱体，可绕其垂直中心轴自由转动，圆柱表面上刻有一倾角为动，圆柱表面上刻有一倾角为的螺旋槽，今在槽中放一的螺旋槽，今在槽中放一小球小球M，自静止开始沿槽下滑，同时使圆柱体绕轴线转动，自静止开始沿槽下滑，同时使圆柱体绕轴线转动，设小球质量为，设小球质量为 ，圆柱体的质量为，圆柱体的质量为 ，半径为，半径为R，不，不计摩擦计摩擦. .1m2m求：当小球下降的高度为求：当小球下降的高度为h h 时时

34、, ,小球相对于圆柱体的速小球相对于圆柱体的速 度，以及圆柱体的角速度度，以及圆柱体的角速度. .解：解： 小球与圆柱体组成的系统是具有两个自由度的系统小球与圆柱体组成的系统是具有两个自由度的系统取圆柱体的转角取圆柱体的转角 和沿螺旋槽方向的弧坐标和沿螺旋槽方向的弧坐标s为广义坐标为广义坐标取小球为动点，圆柱体为动系，则小球的动能取小球为动点，圆柱体为动系，则小球的动能圆柱体的动能为圆柱体的动能为2222222241)2(2121RmRmJT)cos2(2)cos(221212221re2r2e12111sRRsmvvvvmvmT 系统的动能为系统的动能为cos4)2(241122212121

35、sRmRmmsmTTT若选择小球起点为零势能点，若选择小球起点为零势能点，则系统势能则系统势能V可表示可表示sin1gsmV系统的拉格朗日函数系统的拉格朗日函数sincos4)2(24111222121gsmsRmRmmsmVTL由于由于L中不显含时间中不显含时间t 和广义坐标和广义坐标系统有能量积分和循环积分系统有能量积分和循环积分1CT2CVT11221cos22CsRmRmm211222121sincos4)2(241CgsmsRmRmmsm将初始条件将初始条件 t=0 时时000，ss代入上式代入上式得得021 CCcos)2(2211sRmmm令令sinsh ghsmmmm22sin

36、2221221由此得小球相对于圆柱体的速度为由此得小球相对于圆柱体的速度为12r212222sinmmvsghmm得圆柱体转动的角速度为得圆柱体转动的角速度为)sin2)(2(2cos2221211mmmmghRm 1 16 6 第一类拉格朗日方程第一类拉格朗日方程12(,)0(1 2 3)knfrrrtks， ，约束方程为约束方程为 设由设由n个质点组成的系统受个质点组成的系统受s个完整双侧约束个完整双侧约束两边取变分两边取变分)321(01skrrfiniik，其中其中kzfjyfixfrfikikikik引用拉格朗日乘子引用拉格朗日乘子)21(skk， 1111( )() 0snnskkkikikiikiiffrrrr 1() 0nii iiiFmrr11() 0nskii ikiikifFmrrr 在在3n个质点坐标中，独立坐标有个质点坐