8.1 基本方程及其简单应用 基本方程 理想约束的完整有势系统 L=T-V d dt =0, =1, 2,, s q q 存在非有势力的理想约束的完整系统 d dt = Q q q, =1, 2,, s 注 : 通常约束力不出现在动力学方程中 ; 方程组的数目等于 自由度数 ; 每一个方程都是二阶常

Similar documents
高等数学A

untitled

例15


untitled

第七章 分析力学基础

PowerPoint Presentation

(p.29). (a) F Qq r 2 ()() N (b) Q 2 r 2 F ( 2 )() Q 0 5 C 2. (a) F (b) F 3. 7 (p.42). (a) T (b) F (c) T 2. (a) A (b) (c) 4. (a) 4 (b) (

微积分 授课讲义


幻灯片 1

Microsoft Word doc

一、选择:1~8小题,每小题4分,共32分.下列每题给出的四个选项中,只有一个选项是符合题目要求的.

试卷

untitled


第一章三角函数 1.3 三角函数的诱导公式 A 组 ( ) 一 选择题 : 共 6 小题 1 ( 易诱导公式 ) 若 A B C 分别为 ABC 的内角, 则下列关系中正确的是 A. sin( A B) sin C C. tan( A B) tan C 2 ( 中诱导公式 ) ( ) B. cos(

Homework9 答案 1. 粒子在一维势场 V (x) = k x ν 中运动,k, ν > 0. 利用测不准关系估算基态能级和 x. 解 : 粒子在某个态下的能量期望值可以写成 : p 2 E = 2m + k x ν (1) 考虑 x ν x 2 ν/2 V (x) = k x ν 是偶函数

湖北文都考研官网 : 考研数学二考试真题 ( 完整版 ) 来源 : 文都教育 一 选择题 1~8 小题, 每小题 4 分, 共 32 分, 下列每题给出的四个选项中, 只有一个选项是符合题目要求 的. k 1. 当 x 0 时, x tan x与 x 同阶



08-01.indd

untitled

untitled

8 9 < ; ; = < ; : < ;! 8 9 % ; ϑ 8 9 <; < 8 9 <! 89! Ε Χ ϑ! ϑ! ϑ < ϑ 8 9 : ϑ ϑ 89 9 ϑ ϑ! ϑ! < ϑ < = 8 9 Χ ϑ!! <! 8 9 ΧΧ ϑ! < < < < = 8 9 <! = 8 9 <! <

(r s) {φ r1, φ r2,, φ rn } {φ s1, φ s2,, φ sn } u r (t) u s (t). F st ι u st u st k 1 ι φ i q st i (6) r β u r β u r u r(t) max u st r φ

! Ν! Ν Ν & ] # Α. 7 Α ) Σ ),, Σ 87 ) Ψ ) +Ε 1)Ε Τ 7 4, <) < Ε : ), > 8 7

动 Prepared by Jiang Xiao

,..,.,,,,,.,,.,., ,.,,.,,.,, 1,,, ; 2,,,,.,,,,.,,.,,,.,.,.,,.,.,,,.,,,.,,,,.,.,,,, i

# # 4 + % ( ) ( /! 3 (0 0 (012 0 # (,!./ %

2016考研数学三线性代数题目及试题答案

07-3.indd


6.3 正定二次型

第9章 排队论

PowerPoint 演示文稿

&! +! # ## % & #( ) % % % () ) ( %

untitled

Ρ Τ Π Υ 8 ). /0+ 1, 234) ς Ω! Ω! # Ω Ξ %& Π 8 Δ, + 8 ),. Ψ4) (. / 0+ 1, > + 1, / : ( 2 : / < Α : / %& %& Ζ Θ Π Π 4 Π Τ > [ [ Ζ ] ] %& Τ Τ Ζ Ζ Π

Homework7 答案 1. 均匀磁场中的中性自旋 1/2 粒子, 磁场方向 x, 强度 B, n 为任意方向 (θ, φ) 的单位矢量,σ n = n σ 为泡利算符在该方向的投影, 如果初态是 σ n = 1 的本征态, 请解出态的时间演化, 并分别计算 S x, S y, S z 的测量值与


参考文献:

!! )!!! +,./ 0 1 +, 2 3 4, # 8,2 6, 2 6,,2 6, 2 6 3,2 6 5, 2 6 3, 2 6 9!, , 2 6 9, 2 3 9, 2 6 9,

幻灯片 1

untitled

. 微积分课程 微积分 2 复习 2019 年 5 月 2 日 暨南大学数学系 吕荐瑞 (lvjr.bitbucket.io)

幻灯片 1

Microsoft PowerPoint - Eng-math-lecture14.ppt [Compatibility Mode]

4= 8 4 < 4 ϑ = 4 ϑ ; 4 4= = 8 : 4 < : 4 < Κ : 4 ϑ ; : = 4 4 : ;

PowerPoint Presentation

第四章 数值积分与数值微分

(CIP) ().: ISBN O4-44 CIP (2004)


., /,, 0!, + & )!. + + (, &, & 1 & ) ) 2 2 ) 1! 2 2

7. 下列矩阵中, 与矩阵 相似的为. A.. C.. B.. D. 8. 设 AB, 为 n 阶矩阵, 记 rx ( ) 为矩阵 X 的秩,( XY?) 表示分块矩阵, 则 A. r( A? AB) r( A). B. r( A? BA) r( A). C. r A B r A r B (? )

<4D F736F F D BFBCD1D0CAFDD1A7B6FED5E6CCE2BCB0B4F0B0B8BDE2CEF65FCDEAD5FBBEABD7BCB0E65F>

(8) 设 A = ( α α α α) 是 阶矩阵 A 为 A 的伴随矩阵 若 ( ) T 是方程组 A = 的一个基础解系 则 A= 的基础解系可为 ( ) (A) α α (B) α α (C) α α α (D) α α α 二 填空题 (9~ 小题 每小题 分 共 分 请将答案写在答题纸

WL100014ZW.PDF

lim f(x) lim g(x) 0, lim f(x) g(x),


, ( 6 7 8! 9! (, 4 : : ; 0.<. = (>!? Α% ), Β 0< Χ 0< Χ 2 Δ Ε Φ( 7 Γ Β Δ Η7 (7 Ι + ) ϑ!, 4 0 / / 2 / / < 5 02

! /. /. /> /. / Ε Χ /. 2 5 /. /. / /. 5 / Φ0 5 7 Γ Η Ε 9 5 /

1

# # # #!! % &! # % 6 & () ) &+ & ( & +, () + 0. / & / &1 / &1, & ( ( & +. 4 / &1 5,

<3935BCC6A5D2C1CDB6D52E747066>

8 9 8 Δ 9 = 1 Η Ι4 ϑ< Κ Λ 3ϑ 3 >1Ε Μ Ε 8 > = 8 9 =

/ Ν #, Ο / ( = Π 2Θ Ε2 Ρ Σ Π 2 Θ Ε Θ Ρ Π 2Θ ϑ2 Ρ Π 2 Θ ϑ2 Ρ Π 23 8 Ρ Π 2 Θϑ 2 Ρ Σ Σ Μ Π 2 Θ 3 Θ Ρ Κ2 Σ Π 2 Θ 3 Θ Ρ Κ Η Σ Π 2 ϑ Η 2 Ρ Π Ρ Π 2 ϑ Θ Κ Ρ Π

untitled

Microsoft PowerPoint - 概率统计Ch02.ppt [Compatibility Mode]

Β 8 Α ) ; %! #?! > 8 8 Χ Δ Ε ΦΦ Ε Γ Δ Ε Η Η Ι Ε ϑ 8 9 :! 9 9 & ϑ Κ & ϑ Λ &! &!! 4!! Μ Α!! ϑ Β & Ν Λ Κ Λ Ο Λ 8! % & Π Θ Φ & Ρ Θ & Θ & Σ ΠΕ # & Θ Θ Σ Ε

! Β Β? Β ( >?? >? %? Γ Β? %? % % %? Χ Η Ιϑ Κ 5 8 Λ 9. Μ Ν Ο Χ? Π Β # % Χ Χ Θ Ρ% Ρ% Θ!??? % < & Θ

2019 考研数学三考试真题及答案详解 来源 : 文都教育 一 选择题 :1~8 小题, 每小题 4 分, 共 32 分, 下列每题给出的四个选项中, 只有一个选项是符合题目要 求的. k 1. 当 x 0 时, 若 x - tan x 与 x 是同阶无穷小, 则 k = A. 1. B. 2. C

untitled

) Μ <Κ 1 > < # % & ( ) % > Χ < > Δ Χ < > < > / 7 ϑ Ν < Δ 7 ϑ Ν > < 8 ) %2 ): > < Ο Ε 4 Π : 2 Θ >? / Γ Ι) = =? Γ Α Ι Ρ ;2 < 7 Σ6 )> Ι= Η < Λ 2 % & 1 &

数字带通 带阻 高通滤波器的设计 把一个归一化原型模拟低通滤波器变换成另一个所需类型的模拟滤波器, 再将其数字化 直接从模拟滤波器通过一定的频率变换关系完成所需类型数字滤波器的设计 先设计低通型的数字滤波器, 再用数字频率变化方法将其转换成所需类型数字滤波器


第十一章 颈部疾病病人的护理

untitled

幻灯片 1

!! # % & ( )!!! # + %!!! &!!, # ( + #. ) % )/ # & /.


Α 3 Α 2Η # # > # 8 6 5# Ι + ϑ Κ Ι Ι Ι Η Β Β Β Β Β Β ΔΕ Β Β Γ 8 < Φ Α Α # >, 0 Η Λ Μ Ν Ο Β 8 1 Β Π Θ 1 Π Β 0 Λ Μ 1 Ρ 0 Μ ϑ Σ ϑ Τ Ο Λ 8 ϑ

> # ) Β Χ Χ 7 Δ Ε Φ Γ 5 Η Γ + Ι + ϑ Κ 7 # + 7 Φ 0 Ε Φ # Ε + Φ, Κ + ( Λ # Γ Κ Γ # Κ Μ 0 Ν Ο Κ Ι Π, Ι Π Θ Κ Ι Π ; 4 # Ι Π Η Κ Ι Π. Ο Κ Ι ;. Ο Κ Ι Π 2 Η

314 动力学与控制学报 2016 年第 14 卷 x 2 =-x 2 +A(x 1 -x 2 )+B( x 1 - x 2 )+ε(1-x 2 2) x 2 =g 2 (ε 0 )+εg 2 (ε 1 ) (2b) 与系统 (2) 相应的未受微扰作用系统的运动微分方程可表示成 x 1 =-x 1


WL100079ZW.PDF

9!!!! #!! : ;!! <! #! # & # (! )! & ( # # #+

PowerPoint 演示文稿

. /!Ι Γ 3 ϑκ, / Ι Ι Ι Λ, Λ +Ι Λ +Ι

避孕篇


<4D F736F F D2036A1A BFBCD1D0CAFDD1A7B6FED5E6CCE2BCB0B4F0B0B8BDE2CEF6A3A8CEC4B6BCB0E6A3A9>

2 2 Λ ϑ Δ Χ Δ Ι> 5 Λ Λ Χ Δ 5 Β. Δ Ι > Ε!!Χ ϑ : Χ Ε ϑ! ϑ Β Β Β ϑ Χ Β! Β Χ 5 ϑ Λ ϑ % < Μ / 4 Ν < 7 :. /. Ο 9 4 < / = Π 7 4 Η 7 4 =

2009ÄêÈ«¹ú˶ʿÑо¿ÉúÈëѧͳһ¿¼ÊÔÊýѧ¶þÊÔÌâ

, & % # & # # & % & + # & # # # & # % #,

! # % & # % & ( ) % % %# # %+ %% % & + %, ( % % &, & #!.,/, % &, ) ) ( % %/ ) %# / + & + (! ) &, & % & ( ) % % (% 2 & % ( & 3 % /, 4 ) %+ %( %!

2014 年全国硕士研究生入学统一考试 数学三试题 一 选择题 :1~8 小题, 每小题 4 分, 共 32 分, 下列每小题给出的四个选项中, 只有一项符合题目要求 的, 请将所选项前的字母填在答题纸... 指定位置上. (1) 设 lim a = a, 且 a 0, 则当 n 充分大时有 ( )

从分析力学到量子力学

幻灯片 1

= Υ Ξ & 9 = ) %. Ο) Δ Υ Ψ &Ο. 05 3; Ι Ι + 4) &Υ ϑ% Ο ) Χ Υ &! 7) &Ξ) Ζ) 9 [ )!! Τ 9 = Δ Υ Δ Υ Ψ (

Transcription:

第八章 拉格朗日动力学

8.1 基本方程及其简单应用 基本方程 理想约束的完整有势系统 L=T-V d dt =0, =1, 2,, s q q 存在非有势力的理想约束的完整系统 d dt = Q q q, =1, 2,, s 注 : 通常约束力不出现在动力学方程中 ; 方程组的数目等于 自由度数 ; 每一个方程都是二阶常微分方程

简单应用 处理问题的基本步骤 ( 套路固定 ): (1) 判断是否为理想约束的完整有势系统 (2) 判断系统自由度并选择合适的广义坐标 (3) 将 L=T-V 表示成只含广义坐标 广义速度和时间的函数 * (4) 对于有势系统, 将 L 代入拉格朗日方程得到系统的运动微分方程 (5) 对于非有势系统, 还通过定义要求出非有势力对应的广义力, 连同 L 一起代入拉格朗日方程得到系统的运动微分方程注 : 从处理问题角度来看, 拉格朗日方法比较规范, 不需要太多的技巧 不像牛顿方法对同一个问题的处理可以采用多种方案.

例题 1 质量为 m 的质点, 被约束在半顶角为 α 的光滑固定圆锥面上运动, 试通过拉格朗日方程, 写出质点的运动微分方程. 解 : 此为理想约束完整有势系. 建立图示本征系 Oxyz 以及取柱坐标 ρ,θ,z 面内运动自由度为 2, 可选 ρ,θ 为广义坐标. v= e e ż e z =z tan T = m 2 2 csc 2 2 2 V =mgz=mg cot L=T V = m 2 2 csc 2 2 2 mg cot 代入拉格朗日方程得到 csc 2 2 g cot =0 d 2 dt =0 2 =const.

例题 2 求弹簧摆的振动方程. 已知质量为 m 的摆锤挂在轻弹簧上, 弹簧一端固定于点 O 并能在 xy 面内自由旋转. 弹簧的自然长度为 l, 劲度系数为 k. ρ 解 : 理想约束完整有势系. 自由度为 2, 选极坐标 ρ,θ 为广义坐标. T = m 2 2 2 2, V = mg cos k 2 l 2 将 L=T-V 代人拉格朗日方程可得... mg 2 g cos k m l =0 2 g sin =0 大角度摆 对于在小角度平衡位置附近的微振动, 上述方程可化为... g l 0 =0, k m =0, 其中 l 0 =l mg k, = l 0

对称性与守恒量的关系 8.2 守恒定律 定义 对称性: 力学系统在坐标或时间的变换下的不变性 反映为在该变换下系统的拉格朗日函数不变. 例如 : 若,L x, y, z, ẋ, ẏ,ż,t =L x, y, z, ẋ, ẏ, ż,t 我们说系统在 x 方向的平移 ( 变换 ) 下不变, 即具有对称性 容易证明, 这等价于说 / x=0. 定义 连续变换 : 变换过程中系统的坐标或时间 可以表示成随某个参数连续变化的形式. 例如 : q Q, Q 0 =q t, 0 =t, 且, 且 lim 0 lim 0 Q =Q =

诺特定理 对于系统相应于连续变换的对称性, 总有一个系统的守恒量与之对应. 证明 : 首先就空间对称性部分证明. 即考虑变换 q Q, Q 0 =q 而不对时间进行变换, 即时间与参数 无关 我们无非是选了另外一套广义坐标而已 记号 : Q dq d t, Q ' lim 0 Q Q 0, Q ' lim 0 Q Q 0 L ϵ =L(Q α ϵ, Q α ϵ, t), L 0 =L(q α, q α,t) 对称性意味着变换前后系统的拉格朗日函数不变, 即 L( )=L(0) 0=lim ϵ 0 ( L ϵ L 0 ) ϵ 拉格朗日方程 = = q Q ' q = d dt Q ' d dt Q ' Q ' =0 =const.

下面考虑对时间的变换, 为了用上面证明已有的结果, 我们可以 把时间看成某个假想的变量 χ 的函数 : t=t(χ) 记 t' =dt/dχ, q =q t, 则 S= t 1 t 2 2 L q, q,t dt= 1 d q d = d q dt L q, q,t t ' d 2 1 t '= q t ' q = 1 t ' L q, d q d d q d,t, t ', d 于是我们得到了一个扩展的力学系统, 它以 q 和 t 为广义坐标, 新的拉格朗日函数为 L. 考虑对时间的变换 L t ' ' =const. L=L q, q,t t ' q = 1 t ' t, 0 =t 其中 L t ' = ' lim 0 d q d t ' = 1 d q t ' 2 d = q t ', 上页的证明思路仍然有效, 可得 0 t ' t ' L L q '=const. 证毕.

广义动量和广义动量积分 定义 广义动量 p = 注 : 可证明 p = n m n v n r n q 定义 广义坐标平移变换 :q Q =q, q q 推论 若力学系统具有某广义坐标的平移不变性, 则 该广义坐标相应的广义动量守恒. 证明 : 若系统在某广义坐标 q α 的变换下保持不变 根据诺特定理有 因为 Q ' lim 0 Q Q 0 Q ' =const. =1, Q ' =0 q Q =q p = =const. 证毕. 广义动量积分 注 : 也可直接用拉格朗日方程证明...( 可遗坐标 循环坐标 )

例题 3 自由质点在重力场中运动, 试分析其广义动量积分 解 : 建立 Oxyz 本征系, z 轴竖直向上 广义坐标可取 x,y,z, 则 L= m 2 ẋ2 ẏ 2 ż 2 mgz 由于 L 不含 x 和 y, 所以沿 x 和 y 具有平移不变性, 故 p x = ẋ =m ẋ=const. p 如果选取球坐标 (r, θ, φ), 则 y = =m ẏ=const. ẏ 上两式分别表示质点在 x 和 y 方向的动量守恒. L= m 2 ṙ 2 r 2 2 r 2 sin 2 2 mgr cos 由于 L 不含 φ, 所以沿 φ 具有平移不变性, 故 p = =mr 2 sin 2 =const. 此式表示对 z 轴角动量守恒

可以看到 : (1) 广义动量积分存在与否和数量多少, 与广义坐标的选取有关, 故选取适当的广义坐标可以找到较多的广义动量积分. (2) 广义坐标不同, 对应的广义动量的物理意义也不同 如 x 对应的 广义动量即 x 方向的动量, 而 φ 对应的广义动量为对质点 z 轴的 角动量. 也可能选择某种广义坐标, 其对应的广义动量无简单的 物理意义. 推论 若力学系统具有沿空间固定方向平移不变性, 则 系统沿该固定方向的动量守恒. 证明 : 假定该固定方向单位向量为 l, 系统沿 l 平移无穷小量 之后 q Q 系统平移前位矢 r n q, t 系统平移后位矢 r n Q,t =r n q, t l r n q Q q = l r n q Q ' =l

诺特定理 Q ' =const. n ṙ n ṙ n Q ' =const. r 注意到 ṙ n = n 是 q, q, t的函数 q r n q t 只是 q, t 的函数 ṙ n = r n q, n ṙ n r n q Q ' =const. r n q Q ' =l n ṙ n l=const. m L= n n 2 ṙ 2 n V r n, t =m ṙ n ṙ n n p l = n m n ṙ n l= n m n ṙ n l =const. 此即沿 l 方向的动量守恒. ( 证毕 )

推论 若力学系统具有绕空间固定轴的转动不变性, 则 系统对该固定轴的角动量守恒. 证明 : 设该固定轴对应单位矢量 l. 在轴上取 l O 为原点, 把转动前第 n 质点位矢记为 r n 用广义坐标表示为 r n q, t 假定转动无限小角度, 广义坐标变化 q Q 位矢变化 r n q,t r n Q,t =r n q, t l r n q,t r(q α,t) r n Q, t r n q Q q = l r n q,t r n q Q ' =l r n q,t 诺特定理 Q ' =const. n ṙ n ṙ n Q ' =const. n ṙ n [l r n q,t ]=const. L l =l n r n m n ṙ n =const.

广义能量和广义能量积分 定义 广义能量 H = q L 定义 时间平移变换 :t =t 推论 若力学系统具有时间平移不变性, 则系统 广义能量守恒. 证明 : 在时间平移变换下 诺特定理 L H = ' lim 0 q '=const. 0 =1 q L=const. ( 证毕 ) 广义能量积分 注 : 也可直接用拉格朗日方程证明 dh/dt=- L/ t, 然后得到结论.

问 : 广义能量是否就是牛顿力学中质点系的机械能呢? r n q, t v n = d r n dt = r n q T = 1 2 n m nv n v n = 1 2 n m n r n q q r n t q r n t r n q q r n t = 1 2 n m n[ r n q r n q 2 r n t r n q r n t 2] q q q = 1 2 n m n r n r n q q q q 广义速度二次齐次式 T 2 n m n r n t r n q 1 2 n m n r n t 2 q 广义速度一次齐次式 T 1 广义速度零次齐次式 T 0

引理 齐次函数欧拉定理 : 若 f(λx,λy,λz)=λ n f(x,y,z), 则 证明 : 令 χ=λx, η=λy, ζ=λz. 把定义式两边同时对 λ 求偏导可得 然后令 λ=1 即可得证. f x x f y y f z z=n f f x f y f z=n n 1 f 定理 系统广义能量满足 : H=T 2 -T 0 +V 证明 : L=T-V=T 2 +T 1 +T 0 -V T 2 q =2T 2 上一引理 => T 1 q =T 1 T 0 q =0 H = q L=T 2 T 0 V ( 证毕 )

也就是说, 写出 T 后, 要看 T 是不是广义速度的二次齐次式 如果 T=T 2,( 即 T 1 =0,T 0 =0 ), 此时 H=T 2 +V=T+V, 即广义能量等于系统 的机械能. 一般情况下, 广义能量不是系统的机械能. 例题 4 长 2a 质量为 m 的匀质直杆 AB, A 端与光滑水平面接触, 在重力作用下从竖直位置被自由释放倒下. 求杆落地瞬间的角速度. 解 : 由于重力和水平面支持力在竖直面内, y B C 由对称性可知杆一直在竖直面内运动. N φ mg 自由度为 2, 受理想约束的完整系. O A x 建立图示本征系 Oxyz. 选 A 端在 x 轴上的坐标 x 以及杆与 x 轴夹角 φ 为广义坐标. 下面关键计算杆的拉格朗日函数, 主动力是有势的, 以 O 为势能零点, 势能可写为 V =mg a sin x

根据柯尼希定理, 杆的动能可以写为 v C =v A AC T = m 2 v 2 C I C 2 2 而质心速度可以用刚体速度公式求得 v A =ẋ x = z AC=a cos x sin y O y x A φ C v C = ẋ a sin x a cos y v C 2 =ẋ 2 2 a ẋ sin a 2 2 B x T = m 2 ẋ2 2 a ẋ sin a 2 2 ma2 6 2 = m 2 ẋ2 2a ẋ sin 2ma2 3 2 L=T V = m 2 ẋ2 2 a ẋ sin 2 ma2 3 2 mg a sin 首先看看有没有守恒量. L 不含 x, 故有守恒量 p x = ẋ =m ẋ a sin =c 1 ( 水平方向动量守恒 ) 利用初始条件 t=0 时 : ẋ=0, =0 c 1 =0 ẋ=a sin

L 不显含 t, 且 T=T 2, 故机械能守恒 T V = m 2 ẋ2 2a ẋ sin 2 ma 2 3 2 mga sin =E 利用初始条件 t=0 时 : = /2, ẋ=0, =0 E =mga ẋ=a sin 2 6 g 1 sin = a 4 3sin 2 = [ 6 g 1 sin a 4 3 sin ]1/2 2 y O x A φ C B x ( 想想为什么不取 + 号?) 杆落地瞬间, =0 = 3g /2a. 注 : 利用分析力学解题也会优先考虑守恒定律. 另外, 本题用分析 力学并不比牛顿力学有优势, 实际上牛顿力学可直接写出守恒定律

例题 5 质量为 m 的小环 P 被限制在一半径为 R 的光滑大圆环上, 大圆环绕过环心的铅垂轴以角速度 ω 匀速转动. 初始时小环在大环的最高点, 且相对大环静止, 然后无初速地滑下. 试通过存在的第一积分建立小环相对大环的运动微分方程. 解 : 以小环为研究对象, 它是受理想约束的完整系统. 取球坐标, 有 2 个约束方程 r=r, = t 0 因此自由度为 1. 可以取图示角度 θ T = m 2 ṙ2 r 2 2 r 2 sin 2 2 = m 2 R2 2 R 2 2 sin 2 =T 2 T 0 以 O 为势能零点, 势能可表示为 V =mg R cos 由于 L=T-V 不显含时间, 所以广义能量守恒 : H =T 2 T 0 V = m 2 R2 2 m 2 R2 2 sin 2 mgrcos =const.

利用初始条件 t=0 时 : =0, =0 H =mgr H =T 2 T 0 V = m 2 R2 2 m 2 R2 2 sin 2 mgrcos =mgr 2 2 sin 2 2 g / R cos 1 =0 此即小环的运动微分方程. 注 : (1) 本题中小环受到的约束 = t 0 是非定常约束 (2) 本题中广义能量不是相对地面惯性系的机械能 (T T 2 ), 而是相对于固连于大环的非惯性系中小环的总能量. H =T 2 T 0 V = m 2 R2 2 m 2 R2 2 sin 2 mgrcos 相对动能 惯性离心能 ( 想想为什么?) 重力势能

思考题质量为 m 1 三角形楔置于光滑水平面上, 质量为 m 2 半径 为 r 匀质圆柱可沿斜面无滑滚动 (1) 试用牛顿力学的质点系动力学方法寻找系统的守恒量 (2) 试用分析力学的拉格朗日方法寻找系统的守恒量 特别注意 : 柱体质心相对于楔的速度沿斜面方向! 无滑条件 ( 接触点相对速度为零 )

8.3 多自由度系统的微振动 问题的背景 世界是由原子和分子组成的 H H 在真实世界中原子总是在平衡位置附近作小幅振动 V(r) 微振动 微振动 r0? r 弹簧质点系统 ( 多自由度 )

简化模型 ( 透视多自由度振动特征 ) 考虑两质量为 m 的质点, 被限制在光滑水平线上运动, 3 个轻弹簧沿该水平线连接两质点, 中间弹簧的劲度系数为 k 1, 两侧弹簧的劲度系数为 k 2, 它的外端被固定, 两质点静止时各弹簧无伸长. 回顾 : k m x 运动微分方程 m ẍ k x=0 固有频率方程 2 m 1 k=0 自然的问题 : 在多自由度振动中, 上述方程是否有对应的形式?

简正频率 ( 固有频率 ) 和简正模式 如图所示, 选取 x 1 和 x 2 作为广义坐标, 分别表示两质点相对自身平衡位置的位移. ( 推导运动微分方程...) 引入记号 [ M ]=[ m 0 0 m], [ K ]= [ k 1 k 2 k 1 k 1 k 1 k 2], = [ x 1 x 2] 系统质量矩阵 ( 对称正定 ) 系统刚度矩阵 ( 对称正定 ) 定理 系统的运动微分方程为 [M ] [K ] =0 与单自由度比较 : m ẍ k x=0

证明思路 : 系统拉格朗日函数为 L=T V = m 2 ẋ1 2 m 2 ẋ2 2 k 2 2 x 2 1 k 2 2 x 2 2 k 1 2 x x 2 1 2 特点 : 这里动能是广义速度的二次型, 势能是广义坐标的二次型. 代入拉格朗日方程可得 写成矩阵形式即可. ( 证毕 ) m ẍ 1 k 1 k 2 x 1 k 1 x 2 =0 m ẍ 2 k 1 x 1 k 1 k 2 x 2 =0

( 推导频率满足的方程...) 推论 系统的运动微分方程有 ψ=acos(ωt+α) 形式的解, 其中 ω 满足特征方程 证明思路 : det 2 I [ M ] 1 [ K ] =0 2 m 1 k=0 或 det [ K ] 2 [ M ] =0 ( 单自由度对应 ) [M ] [K ] =0 [ M ] 1 [K ] =0 这个方程类似于谐振子方程, 故可猜测有 ψ=acos(ωt+α) 形式的解. 将 ψ=acos(ωt+α) 代入 [ M ] 1 [ K ] =0 可得 [M ] 1 [K ] A= 2 A 可见 A 是矩阵 [M ] 1 [K ] 对应于特征值 2 的特征向量. 故有 det [ M ] 1 [K ] 2 I =0 det [ K ] 2 [ M ] =0( 证毕 )

注 : 将矩阵 [ M ]=[ m 0 0 m], [ K ]= [ k 1 k 2 k 1 k 1 k 1 k 2] 代入特征方程 det [ K ] 2 [ M ] =0 det[ k 1 k 2 2 m k 1 k 1 k 1 k 2 2 m] =0 k 1 k 2 m 2 2 k 1 2 =0 = k 2 /m 或 2 k 1 k 2 /m 定义 简正频率 : 从特征方程解出来的频率 注 : 它只依赖于系统的质量矩阵和刚度矩阵, 是系统固有的. 定义 简正模式: 简正频率对应的振动模式, 即矩阵 [M ] 1 [K ] 对应于简正频率平方的特征向量决定的振动模式.

例题 6 分析系统的简正模式. 解 : 将 代入 [M ] 1 [K ] A= 2 1 = k 2 /m 1 A 或 [ K ] 1 2 [M ] A=0 可以得到 A [ 1 1] 1= [ x 1 2] x = A cos 1t 1 [ 1 1] cos 1t 1 此时二质点的步调完全一致, 称为对称模式. [ K ] 2 2 [M ] A=0 将 2 = 2 k 1 k 2 /m 代入可得 A [ 1 1] 2= [ x 1 x 2] = A cos 2 t 2 [ 1 1] cos 2 t 2 此时二质点的步调完全相反, 称为反对称模式. 注 : 系统微分方程的通解是上面两种模式的线性叠加, 可写为 =c 1 1 c 2 2 =c 1[ 1 1] cos 1 t 1 c 2[ 1 1] cos 2 t 2 这里 4 个待定实常数 c 1,c 2,α 1,α 2 由初始条件来确定.

简正坐标 定义 简正坐标 : 一组特殊的广义坐标 {ξ α }, 使得系统的 运动微分方程表现为一系列简单的振动 ( 2 =0). 这意味着系统的拉格朗日函数可表示为如下形式 L= m 2 2 2 2. why? 例题 7 寻找图示系统简正坐标. 解 : 我们已经知道得到 [M ] 1 2 2 [K ] 对应特征值 1 和 的特征向量分别为 [ 1 1] 2 和 [ 1. 1 经过归一化后得到和 1] 2[ 1 1] 2[ 1 1] 1. 根据线性代数知, 存在矩阵 2[ [U ]= 1 1 1 1 1] 使得

[U ] 1 [ M ] 1 [ K ][U ]=[ 2] 2 1 0 0 2 [U ] 1 [ M ] 1 [ K ]=[ 2 1 0 0 2 2] [U ] 1 系统运动微分方程 =[ [U ] 1 [ M ] 1 [ K ] 2] 2 1 0 0 2 [U ] 1 [ M ] 1 [ K ] =0 [U ] [ 1 2 1 0 2 2] 0 [U ] 1 =0 只需令 = [ 1 2] =[U ] 1 即可得到 2 =0, =1,2. 1 = x 1 x 2 / 2, 2 = x 1 x 2 / 2 此即简正坐标 x 1 = 1 2 / 2, x 2 = 1 2 / 2 L= m 2 ẋ1 2 m 2 ẋ2 2 k 2 2 x 2 1 k 2 2 x 2 2 k 1 2 x 2 x 1 2 = = m 2 1 2 1 2 1 2 m 2 2 2 2 2 2 2

例题 8 分析双摆简正频率和简正坐标. 设线长均为 l 且不可伸长, 小球半径可以忽略, 质量均为 m. 假定摆角很小. 解 : 理想约束的完整系统, 自由度为 2, 可选摆角 θ 1 和 θ 2 为 θ 1 广义坐标 当摆角很小时, 两个球的速度分别可以写为 v 1 =l 1, v 2 v 1 l 2 =l 1 2 θ 2 动能 T = 1 2 ml2 1 2 1 2 ml2 1 2 2 以最高悬挂点为势能零点, 则势能 V = mgl cos 1 mgl cos 1 cos 2 mgl 1 2 2 2 /2 3 将 L=T-V 代入拉格朗日方程有... 2 1 2 2 g/l 1 =0 1 2 g /l 2 =0 [M ]=[ 2 1 1 1], [K ]= g l [ 2 0 0 1], 代入特征方程 det [K ] 2 [M ] =0求得... = 2± 2 g /l

将 [ K ] 2 1 = 2 2 g /l 代入 1 [M ] A=0 可得 A [ 1 2] 1 = [ 1 2] = Acos 1 t 1 [ 1 2] cos 1 t 1 对称模式 将 [ K ] 2 2 = 2 2 g/l 代入 2 [ M ] A=0 可得 A [ 1 2] 2 = [ 1 2] = A cos 2 t 2 [ 1 2] cos 2 t 2 反对称模式 3[ [U ]= 1 1 2] 1 2 = [ 1 2] =[U 2 =[ 3 ] 1 3 2 3 ] 2 2 3 2 2 [ 1 2] 简正坐标

总结 ( 多自由度振动 ) 系统的运动微分方程 : [M ] [K ] =0 简正频率 ω 满足特征方程 : det 2 I [ M ] 1 [ K ] =0 或 det [ K ] 2 [ M ] =0 简正模式 : 简正频率对应的振动模式 [ M ] 1 [ K ] A= 2 A 简正坐标 : =[U ] 1 2 =0

8.4 非惯性系中的运动 哈密顿原理以及拉格朗日方程本身不受参考系的限制. 选择不同的参考系在某种意义上来看相当于选取了不同的广义坐标来描述一个力学系统而以. 因此, 力学系统的拉格朗日函数在不同参考系中都是同一个. 最多相差一个坐标和时间函数的全导数. 由于采用了不同的广义坐标, 所以表现形式可能很不同. 平动参考系 S' S 0 : 惯性系 S': 平动系 S' 系以已知速度 u(t) 相对 S 0 系做平动. 记 :w= d u dt 考虑质点的运动, 在惯性系 S 0 中的拉格朗日函数可写为 L 0 = mv 0 2 2 V r 0,t m d v 0 dt = V r 0

质点在惯性系 S 0 和平动参考系 S' ( 相对 S 0 以 u 运动 ) 之间满足关系 r 0 =r' u t, v 0 =v' u u 事先给定, 非匀速 L 0 = mv 2 0 2 V r 0, t = m v' 2 2 2 mu m v' u V r ' u t, t 2 m v' u=m d r' dt u = d dt mu r ' m d u dt r ' 含时间和坐标的函数 的全导数, 可略. 拉格朗日方程 d dt ' v' ' r ' w =0 m d v' dt 仅含时间, 可写成含时间 函数的全导数, 故可略. 故可以将拉格朗日函数用广义 坐标 r' 和广义速度 v' 表示为 L '= m v ' 2 m w r ' V r ' u t, t 2 相当于惯性力势能 = V r ' m w

任意非惯性系 S' 系以已知速度 u(t) 相对惯性系 S 0 平动. 而 S 与 S' 原点重合, 且以角速度 ω 相对 S' 系转动. 记 w=du/dt. r 0 r' ut S 0 : 惯性系 S: 转动系 S': 平动系 质点在 S 和 S' 中的位矢是重合的, 故 r=r' 质点相对 S 的速度 v= d' r => 质点相对 S' 的速度 v '= r v dt L '= mv' 2 mw r ' V r ' u t, t 2 = m v2 2 m r v m 2 r 2 m w r V r u t, t L r,v, t 代入拉格朗日方程 d ' dt v r =0 可得 m d ' v dt = V r m w m( ω r) m ω (ω r) 2 m(ω v) 牵连惯性力科氏力

例题 5' 质量为 m 的小环 P 被限制在一半径为 R 的光滑大圆环上, 大圆环绕过环心的铅垂轴以角速度 ω 匀速转动. 初始时小环在大环的最高点, 且相对大环静止, 然后无初速地滑下. 试通过存在的第一积分建立小环相对大环的运动微分方程. 解 : 以大环作为参考系, 并建立固连坐标系 Oxyz. 理想约束的完整系统, 自由度为 1, 选 θ 作为 广义坐标. v=r cos y sin x v = z, r=r sin y cos z L= mv2 2 m r v m 2 r 2 mw r V r u t,t = m R2 2 2 2 sin 2 2 sin 2 mgrcos L 不含时间, 广义能量守恒 H = m R2 2 2 2 sin 2 mgrcos =const. 考虑初始条件 2 2 sin 2 2 g/ R cos 1 =0

拉格朗日方法的特点和意义 拉格朗日方程是完整系最普遍的动力学方程. ( 1 ) 不过问约束力, 方程个数减到最低限度. ( 2 ) 高度的概括性和统一性 : 惯性系和非惯性系 ; 任何广义坐标 ; 质点 质点系 刚体, 无限自由度的连续介质系统 其他物理领域 ( 如电磁场 ). ( 3 ) 程序统一. ( 4 ) 较少采用几何 矢量的方法, 较多采用能量和数学分析的方法 利于推广到物理学的其他领域 只要找到拉格朗日函数, 可用同一形式的拉格朗日方程求出系统的运动方程 因此, 理论和实用性上, 都显示出很大的优越性.