分析力学中不同解题方法的比较.docxVIP

分析力学中不同解题方法的比较 分析力学中不同解题方法的比较 第18卷第1期 沈阳师范学院学报(自然科学版) 2000年1月 Jour nal of Shenyang N or mal Univ ersity (N atural Science) 文章编号:1008-374X(2000) 01-0026-04 Vol. 18, N o. 1Jan. 2000 分析力学中不同解题方法的比较 (沈阳师范学院物理系, 辽宁沈阳 110031) 哈密顿原理及哈密顿正则方程求出同一力学体系的运动规律, 从摘 要:用拉格朗日方程、 而比较了三种方法的异同. 关键词:拉氏函数; 哈密顿函数 中图分类号:O316 文献标识码:A 对于完整保守力学体系, 可以分别采用拉格朗日方程、哈密顿原理及哈密顿正则方程求出体系的运动规律. 下面选用同一例题, 来比较它们的异同 . 例题:如图1所示, 轴为竖直且顶点在 下的抛物线形金属丝, 以匀角速度 绕轴转动. 一质量为m 的小环, 套在此金属丝上, 并可沿着金属丝滑动, 金属丝是光滑的. 已知抛物线方程为x =4ay , 式中a 为常数. 求小环在x 方向的运动微分方程. 1 用拉格朗日方程求解 小环在重力及约束力作用下运动, 而且金属丝是光滑的, 故小环为完整保守力学体系. 由约束方程x 2=4ay , 小环自由度为1, 选取小环坐标x 为广义坐标. 由y =, y =4a 2a , 小环的动能为 收稿日期:1999-08-31 作者简介:姜玉明(1943-) , 女, 吉林集安人, 沈阳师范学院副教授. 2222222 T =2m (x +y + x ) =2m (x +2+ x ) 小环的势能为V =mgy =mg 4a L =T -V =m (x ++ x ) -mg 24a 4a 2 完整保守系的拉格朗日方程为 -) =0 ( =1, 2, , S ) d t ( 对此题来说, q =x , q , S =1, 计算 =x 2() =m (2x ! -+2) d t 22a 2a (5) (6) =m (+2 2x ) -mg 22a 2a 将(5) 、(6) 代入(4) 中, 得 22m (1+! +m ) x -m x +mg 2a =04a 4a 上式就是应用拉格朗日方程求出的小环沿x 方向运动微分方程. 2 用哈密顿原理求解 完整保守系的哈密顿原理为 ?L d t = 将(3) 代入上式, 得 [m (x +]d t =02+ x ) -mg t 24a 4a 1 进行计算, 得 2[-m x ! x -x -x + x +m 2x x -mg x ]d t =02 2a 4a 当t =t 1, t =t 2时, x =0, 则上式为 22[-m x ! -! -] x d t =0x +m x -mg t 12a 4a 4a 又x 是独立的, x 是任意的, 使上式为零, 只有积分函数为零, 即 222m (1+! +=0) x -m x +mg 2a 4a 4a 上式就是应用哈密顿原理求出的小环沿x 方向运动微分方程, 与(7) 相同. 3 用哈密顿正则方程求解 完整保守系的哈密顿正则方程为 ( =1, 2, , S ) (13) p =- 其中p =是广义动量, H 是哈密顿函数, H = p q -L , 或者在稳定约束下, H =T +V , 在非稳定约束下, H =T 2-T 0+V. 此题的约束为非稳定约束, 从(1) 中可知 22T 2=2m (x +) T 0=m 2x 2 (14) (15) H =2m (x +2) -2m x +mg 4a 4a 由(1) 求出 2p x ==m x + x 4a 解出 将上式代入(16) 中, 得出应求出的H 表达式为p x 2222 H =() -m x +mg 2m 4a 2+x 224a 将上式代入(13) 中的p , 计算得x =-4a 2p 2x x 2p x =222+m x -mg 2 (18) 与(20) 为用正则方程所求两个一阶微分方程式. 由此题要求, 将(17) 对t 求导数等于(20) , 得 4a 2p 2x x 22m (1+) x