自学内容网 自学内容网

Greiner 经典力学(多体系统和哈密顿力学)第三章 学习笔记 (Foucault‘s Pendulum)

第三章 学习笔记 (Foucault’s Pendulum)

傅科摆是个经典的物理实验,这个实验证明了地球的自转。利用上一节获得的任意坐标系下的运动方程,可以很简洁的推导出傅科摆的各种性质。

首先,以傅科摆所在的位置建立一个局部坐标系。那么有:
m   r ′ ¨ ∣ M = F − m   R ¨ ∣ L − m ω ˙ ∣ M × r ′ − 2 m ω × r ′ ˙ ∣ M − m ω × ( ω × r ′ ) = F − m   ω × ( ω × R ) − 2 m ω × r ′ ˙ ∣ M − m ω × ( ω × r ′ ) m \ \ddot{\mathbf r'}|_M = \mathbf F - m \ \ddot{\mathbf R}|_L - m \dot{\omega}|_M \times \mathbf r' - 2 m\omega \times \dot{\mathbf r'}|_M - m\omega \times (\omega \times \mathbf r') \\ = \mathbf F - m \ \omega \times (\omega \times \mathbf R) - 2 m\omega \times \dot{\mathbf r'}|_M - m\omega \times (\omega \times \mathbf r') m r¨M=Fm R¨Lmω˙M×r2×r˙M×(ω×r)=Fm ω×(ω×R)2×r˙M×(ω×r)
上面推导中利用了地球是匀速旋转的,所以 ω ˙ = 0 \dot \omega = 0 ω˙=0 。傅科摆收到两个外力,分别是张力 T \mathbf T T 和万有引力 G \mathbf G G
m   r ′ ¨ ∣ M = T + G − m   ω × ( ω × R ) − 2 m ω × r ′ ˙ ∣ M − m ω × ( ω × r ′ ) m \ \ddot{\mathbf r'}|_M = \mathbf T + \mathbf G - m \ \omega \times (\omega \times \mathbf R) - 2 m\omega \times \dot{\mathbf r'}|_M - m\omega \times (\omega \times \mathbf r') m r¨M=T+Gm ω×(ω×R)2×r˙M×(ω×r)
这里万有引力和地球旋转产生的离心力共同的作用效果就是重力: G − m   ω × ( ω × R ) = m g \mathbf G - m \ \omega \times (\omega \times \mathbf R) = m \mathbf g Gm ω×(ω×R)=mg。因此上式可以进一步写为:
m   r ′ ¨ ∣ M = T + m g − 2 m ω × r ′ ˙ ∣ M − m ω × ( ω × r ′ ) m \ \ddot{\mathbf r'}|_M = \mathbf T + m \mathbf g- 2 m\omega \times \dot{\mathbf r'}|_M - m\omega \times (\omega \times \mathbf r') m r¨M=T+mg2×r˙M×(ω×r)
由于地球自转角速度很慢,所以等式右边第四项相比其他项可以忽略。因此上式进一步简化后得到:
m r ¨ = T + m g − 2 m ω × v m \ddot r = \mathbf T + m \mathbf g - 2m \omega \times \mathbf v mr¨=T+mg2×v

上面的式子全都是使用地面坐标系,所以我们省略了下标和上角标。

下面就是把矢量方程化为标量方程。先把张力 T \mathbf T T 分解为三个方向的分量。这个比较简单,参考下图很容易就可以写出分量。

在这里插入图片描述
T = − x l e ^ 1 + − y l e ^ 2 + l − z l e ^ 3 T = \frac{-x}{l} \mathbf {\hat e}_1 + \frac{-y}{l} \mathbf {\hat e}_2 + \frac{l-z}{l} \mathbf {\hat e}_3 T=lxe^1+lye^2+llze^3

实际上 x , y , z x,y,z x,y,z 并不是独立的,质点只能在一个球面上运动。

下面还需要把 $\omega $ 分解为三个方向的分量,这个可以参考下面的图:

在这里插入图片描述

有了 $\omega $ 的分量表示后就可以求 ω × v \omega \times \mathbf v ω×v
ω × v = ∣ e 1 e 2 e 3 − ω sin ⁡ λ 0 ω cos ⁡ λ x ˙ y ˙ z ˙ ∣ = − ω cos ⁡ λ y ˙ e ^ 1 + ω ( cos ⁡ λ x ˙ + sin ⁡ λ z ˙ ) e ^ 2 − ω sin ⁡ λ y ˙ e ^ 3 \omega \times \mathbf v = \begin{vmatrix} e_1 & e_2 & e_3\\ -\omega \sin \lambda & 0 & \omega \cos \lambda \\ \dot x & \dot y & \dot z \end{vmatrix} \\ = -\omega \cos \lambda \dot y \hat {\mathbf e}_1 +\omega (\cos \lambda \dot x + \sin \lambda \dot z) \hat {\mathbf e}_2 -\omega \sin \lambda \dot y \hat {\mathbf e}_3 ω×v= e1ωsinλx˙e20y˙e3ωcosλz˙ =ωcosλy˙e^1+ω(cosλx˙+sinλz˙)e^2ωsinλy˙e^3
到这里就可以写出运动方程的分量形式了。
m x ¨ = − x l T + 2 m ω cos ⁡ λ y ˙ m y ¨ = − y l T − 2 m ω ( cos ⁡ λ x ˙ + sin ⁡ λ z ˙ ) m z ¨ = l − z l T + 2 m ω sin ⁡ λ y ˙ m \ddot x = -\frac{x}{l} T + 2 m \omega \cos \lambda \dot y \\ m \ddot y = -\frac{y}{l} T - 2 m \omega (\cos \lambda \dot x + \sin \lambda \dot z) \\ m \ddot z = \frac{l- z} {l} T + 2m \omega \sin \lambda \dot y mx¨=lxT+2cosλy˙my¨=lyT2(cosλx˙+sinλz˙)mz¨=llzT+2sinλy˙
傅科摆很长,可以认为 z z z 方向式不动的。那么上面式子可以简化为:
m x ¨ = − x l T + 2 m ω cos ⁡ λ y ˙ m y ¨ = − y l T − 2 m ω cos ⁡ λ x ˙ m \ddot x = -\frac{x}{l} T + 2 m \omega \cos \lambda \dot y \\ m \ddot y = -\frac{y}{l} T - 2 m \omega \cos \lambda \dot x \\ mx¨=lxT+2cosλy˙my¨=lyT2cosλx˙
另外就是绳子的张力约等于重力: T ≈ m g \mathbf T \approx m\mathbf g Tmg。所以上式进一步可以消去 T T T m m m
x ¨ = − x l g + 2 ω cos ⁡ λ y ˙ y ¨ = − y l g − 2 ω cos ⁡ λ x ˙ \ddot x = -\frac{x}{l} g + 2 \omega \cos \lambda \dot y \\ \ddot y = -\frac{y}{l} g - 2 \omega \cos \lambda \dot x \\ x¨=lxg+2ωcosλy˙y¨=lyg2ωcosλx˙
我们设 g / l = k 2 , ω cos ⁡ λ = α 2 g/l = k^2, \omega \cos \lambda = \alpha^2 g/l=k2,ωcosλ=α2,那么上式可以化为:
x ¨ = − k 2 x + 2 α 2 y ˙ y ¨ = − k 2 y − 2 α 2 x ˙ \ddot x = -k^2 x + 2 \alpha^2 \dot y \\ \ddot y = -k^2 y - 2 \alpha^2 \dot x \\ x¨=k2x+2α2y˙y¨=k2y2α2x˙
下面的工作就是解这个微分方程组,就不详细写了。需要强调的是上面推导过程用到了大量的近似,这些近似都是有明确的物理意义的。解决物理问题一个很重要的能力就是学会把复杂的问题一步步简化为可以求解的简单问题。


原文地址:https://blog.csdn.net/liyuanbhu/article/details/142678008

免责声明:本站文章内容转载自网络资源,如本站内容侵犯了原著者的合法权益,可联系本站删除。更多内容请关注自学内容网(zxcms.com)!