q第十七章普遍原理与拉氏方程

1、 如何研究较复杂的非自由质点系的动力学问题？ 引言12摆长不定，如何确定其摆动规律？K混沌摆问题多杆摆问题达朗伯（1717-1785）通过引入惯性力的概念，建立了著名的达朗伯原理（用静力学建立平衡方程的方法处理动力学问题）；约翰伯努利（1667-1748）于1717年精确表述了虚位移原理（建立虚位移、虚功的概念，用动力学的方法研究静力学中的平衡问题）；拉格朗日应用达朗伯原理，把虚位移原理推广到非自由质点系的动力学问题中，建立了动力学普遍方程，进一步导出了拉格朗日方程。 引言17-1 动力学普遍方程 动力学普遍方程是虚位移原理与达朗伯原理简单 结合的产物。设质点系由n个质点组成，第i个质点质量为

2、mi，受力加速度为ai虚加上其惯性力Fgi=miai则根据达朗伯原理， Fi 、FNi 与Fgi应组成形式上的平衡力系，即Fi + FNi +Fgi= 0若质点系受理想约束作用，应用虚位移原理，有或(171)FiFNiFgiaiMFNiFNiMMFgiaiFgiaiFiFi约束反力FNi主动力Fi在理想约束条件下，质点系的各个质点在任意瞬时所受的主动力和虚加的惯性力在任一虚位移上所做虚功之和等于零。即则D-L方程的坐标分解式为若动力学普遍方程或达朗伯-拉格朗日方程（D-L方程）(172)例17-1 两均质轮质量皆为m1，半径皆为r，对轮心的转动惯量为J；中心用质量为m2的连杆连接，在倾角为的斜

3、面上纯滚动。求连杆的加速度。 见后续例17-1续已知轮m1，r，J，纯滚；杆m2，求杆a。解：研究整个系统，进行受力分析；m2gm1gm1gN1N2F1F2Fg1Fg2Fg1MgMg虚加各刚体的惯性力。设杆的加速度为a，则Fg1= m1a，Fg2= m2a，a给连杆以平行于斜面向下的虚位移s，则相应地两轮有转角虚位移，且根据动力学普遍方程，得s于是解得m2gm1gm1gN1N2F1F2Fg1Fg2Fg1MgMgasm2gm1gm1gN1N2F1F2Fg1Fg2Fg1MgMgm2gm1gm1gN1N2F1F2Fg1Fg2Fg1MgMgas解毕。解毕。17-2 拉格朗日方程将动力学普遍方程用独立的

4、广义坐标表示，即可推导出 第二类拉格朗日方程。拉格朗日拉格朗日 （Lagrange 1736 1813）法籍意大利人，数学家、力学家、天文学家，十九岁成为数学教授，与欧拉共同创立变分法，是十八世纪继欧拉后伟大的数学家。 设质点系由n个质点组成， 用q1，q2，qN表示系统的广义坐标，D-L方程可写成对上式求变分得 ri= ri(q1，q2，qN，t)(173)上式中(174)第i个质点质量为mi，矢径为ri。则具有s个完整理想约束，则有N=3n-s个自由度（广义坐标）。广义力 根据虚位移原理中广义力与广义虚位移的表示形式，有 因为系统为完整约束，广义坐标相互独立，所以广标的变分qk是任意的，为

5、使上式恒成立，须有(175)广义惯性力以广义坐标表示的达朗伯原理广义力广义虚位移(k =1，2，N)对式(176)中广义惯性力进行变换：将下列两个恒等式（有关证明请参阅教材P282）代入得(177)(178)(1710)所以将式（17-10）代入式（17-5）得 这就是第二类拉格朗日方程，(1711) 若作用于质点系的主动力均为有势力（保守力），于是，对保守系统，拉格朗日方程可写成(1712)是一个方程组，该方程组的数目等于质点系的自由度数，各方程均为二阶常微分方程。揭示了系统动能的变化与广义力之间的关系。则广义力Qk可写成质点系势能表达的形式用函数L表示系统的动能T与势能V之差，即L = T

6、V(1713)L称为拉格朗日函数或动势。势能V仅仅是广义坐标qk的函数，而与广义速度无关，有 于是，在保守系统中，用动势表示的拉格朗日方程的形式为(1714)对保守系统，拉格朗日方程可写成 拉格朗日方程的形式 在保守系统中，用动势表示的拉格朗日方程的形式为 拉格朗日方程是解决具有完整约束的质点系动力学问题的普遍方程，是分析力学中重要的方程。 拉格朗日方程形式简洁，应用时只需计算系统的动能和广义力；对于保守系统，只需计算系统的动能和势能。 拉格朗日方程是标量方程，以动能为方程的基本变量，是用广义坐标表示的质点系运动微分方程。对受完整约束的多自由度质点系的平衡问题，根据虚位移原理，采用广义坐标，得

7、到与自由度相同的一组独立平衡方程。这种用分析方法建立的平衡条件，避开了未知的约束反力，使非自由质点系的平衡问题的求解变得简单。同理，对受完整约束的多自由度质点系的动力学问题，可以根据能量原理，采用广义坐标，推导出与自由度相同的一组独立的运动微分方程。这种用广义坐标表示的动力学普遍方程，称为拉格朗日第二类方程，或简称拉格朗日方程。拉格朗日方程概述拉格朗日方程概述拉格朗日方程是着眼于整个系统，避开约束反力，用分析方法给出了系统动力学问题的统一表述，为处理受约束的复杂的系统动力学问题开辟了新的捷径。这里从能量原理中的功率方程出发，采用广义坐标，推导拉格朗日方程。由于拉格朗日方程是用广义坐标且从能量的

8、观点研究系统的动力学问题，而能量是自然界各种不同物理形态的物质运动的统一度量，因此，拉格朗日方程的应用就具有较大的普遍性，它不仅适用于机械系统，也适用于电学系统和机电系统的动力学问题。拉格朗日方程概述续设质点系由n个质点组成，各质点系的速度则为jq 所以vi是广义速度的齐次线性函数。其中是广义坐标的函数，设质点系中第i个质点的质量为mi，位置矢径为ri，系统受s个定常的完整的理想约束， 因此有N =3n-s个自由度。用k个广义坐标q1、q2、qn表示质点的位置。由于约束是定常的，所以矢径ri 仅是广义坐标的函数用功率方程推导拉格朗日方程用功率方程推导拉格朗日方程又系统的动能T为其中 称为广义质

9、量，是广义坐标的函数。现令系统的动能为T，功率为P，则功率方程为所以，在定常约束中，动能T不是时间t的函数，仅是广义坐标和广义速度的函数。将上式对时间求全微商，有又由于动能T仅含的二次项，代入上式可得根据欧拉齐次函数定理将系统的动能T写成一般式是广义速度 的齐二次函数。再看功率方程的右边，功率为根据虚位移原理中广义力Qk的表达式可知所以将(3) (4)式代入(2)式并移项可得将(5)式改写成对受定常约束的系统，微小的实位移必为虚位移中之一，即dqk必为qk之一，但因各广义坐标qk是彼此独立的，为使上式在各qk具有任何值都能满足，必须即上式即为以广义坐标表达的质点系运动微分方程，称为拉格朗日方程

10、。为了使系统在任何初始条件下可能发生的位移都能满足方程（6），就必须用拉格朗日方程解题的步骤： 确定系统的自由度数（广义坐标数）； 选广义坐标； 计算系统的动能T，且用广义速度来表示动能； 计算广义力（对保守系统可计算势能）； 代入拉格朗日方程即可得质点系运动微分方程。例17-2 位于水平面内的行星轮机构中，质量为m1的均质系杆OA，可绕O轴转动，另一端装有质量为m2、半径为r的均质小齿轮，小齿轮沿半径为R的固定大齿轮纯滚动。当系杆受力偶M的作用时，求系杆的角加速度。 见后续rRMMO AM例17-2续 已知杆m1；轮m2，r；R；力偶M。求OA。解：研究整个系统，选广义坐标，则行星轮瞬心为P

11、，系统的动能为T = TOA+ T轮ARMrO 见后续P P P角速度为vAvAvA例17-2续已求得用广义坐标表示的系统动能为O AvARMr又关于广义坐标的广义力为代入Lagrange方程：于是得解毕。解毕。O例17-3 质量为m的质点悬在不计质量的软线上，线的另一端绕在半径为R的固定圆柱上。设在平衡位置时，线的下垂部分长度为l。求此摆的运动微分方程。 见后续Rml l lRlO例17-3续1m已知m；R ； l；求摆的运动微分方程。m解：此摆为单自由度保守系统，选广义坐标，系统的动能为选=0处为系统势能的零势点，则V = mg（l+Rsin）（lR）cos系统的动势为 见后续例17-3续

12、2已求得将式、代入保守系统的拉氏方程得摆的运动微分方程解毕。O O O例17-4 已知质量为m1的三棱柱放在光滑水平面上，质量为m2的均质圆柱体O由静止沿三棱柱的斜面向下纯滚动；求三棱柱的加速度。 见后续OO解(设圆柱O的半径为r)选x1、x2为广义坐标，x1x2O则三棱柱速度为圆柱中心的速度为 圆柱的角速度为vO系统具有两个自由度，o1o2所以，系统的动能为 见后续加速度为例18-4续1 已知三棱柱m1；圆柱m2，纯滚动；求三棱柱的加速度。x1x2vOx1x2vOx1x2vOx2x2x2x2系统关于广义坐标x1 、x2的广义力分别为：系统受力如图，例17-4续2已求得系统的动能为x1x2Oo

13、1o2m1gFNm2gx1联立解得：代入L程：m1gFNm2gx1m1gFNm2gx1m1gFNm2gx1例17-5 图示均质杆AB质量为m1，长为3l，B端铰接一质量为m2，半径为r的均质圆盘。杆AB在O处为铰支，两弹簧的刚性系数均为k；杆在水平位置平衡。求系统的微幅振动的固有频率。 见后续Okkllll2rACB例17-5续1 已知杆m1，长3l；盘m2，半径r，弹簧k。求：系统微幅振动的固有频率。解系统具有两个自由度，且为保守系统。选1、2为广义坐标，OkkllllACB2r12则杆的角速度为圆盘的角速度为 所以，系统的动能为 见后续1212系统的势能为OkACkllll2rB2m1gm

14、1gm1gm1gm2gm2gm2gm2g1l1l1F1F2l1l1l1l1重力与振动方向相同，系统受力如图，例17-5续2 已求得系统的动势为 见后续F1F2F1F2F1F2取平衡位置处为零势点，弹性力变形从平衡位置处计算，可以不计重力势能！例18-5续3已求得代入保守系统的拉氏方程可见，圆盘的角加速度为零！圆盘作平动！得所以系统的固有频率为圆盘作平动！圆盘作平动！解毕。解毕。OkACkllll2rB2m1gm2g1l1l1F1F2*例17-6 杆OA与AB以铰链相连，且OA=a，AB=b，O悬挂于圆柱铰链上， A、B处质点质量分别为 m1和m2，各处摩擦及两杆质量均不计，求系统微幅摆动的微分

15、方程。 见后续m1bam2OABvAvAvAvA例17-6续1 已知OA=a，AB=b； m1、m2，求系统微摆方程。baOAB12解选1、2为广义坐标，则系统具有两个自由度，系统动能为 见后续vAvB2 1vBA12vAvB2 1vBA12vAvB2 1vBA12vAvB2 1vBA系统作微幅摆动，cos(21)11221例17-6续2已求得求系统关于广义坐标1的广义力：系统受力如图。m2g求系统关于广义坐标2的广义力：112XOYOm2gm1g2m1gm2gXOYOm1gYOm2gXOm1gm2gYOm1gXO1122b2b2b2a1a1a1a1a1a1 见后续给1，则给2，则例17-6续3已求得代入Lagrange方程：化简得解毕。解毕。 本章小结 动力学普遍方程是将虚位移原理与达朗伯原理结合起来，形成如下的方程