第一章题目为 “Newton’s Equations in a Rotating Coordinate System”。这一章内容比较简单。这里只对一些书上没有交待,但是容易混淆的知识点做一些澄清。
一个位置变量 A ( t ) A(t) A(t) 表示的是空间的一个点。这个点在不同的坐标系下表现形式是不同的。但是我们写 A ( t ) A(t) A(t) 时不用说明 A ( t ) A(t) A(t) 是在哪个坐标系下的。因为无论在哪个坐标系下都表示的是同一个点。
A ( t ) = A ′ ( t ) = A 1 e ^ 1 + A 2 e ^ 2 + A 3 e ^ 3 = A 1 ′ e ^ 1 ′ + A 2 ′ e ^ 2 ′ + A 3 ′ e ^ 3 ′ \mathbf A(t) = \mathbf A'(t)= A_1 \mathbf{\hat e}_1 + A_2 \mathbf{\hat e}_2 + A_3 \mathbf {\hat e}_3 = A_1' \mathbf{\hat e_1'} + A_2' \mathbf{\hat e_2'} + {A_3'} \mathbf{\hat e_3'} A(t)=A′(t)=A1e^1+A2e^2+A3e^3=A1′e^1′+A2′e^2′+A3′e^3′
一个速度变量 v ( t ) = d A ( t ) / d t \mathbf v(t) = d \mathbf A(t)/dt v(t)=dA(t)/dt 是与坐标系直接关联的。不同坐标系下速度是不同的。所以我们写速度时需要指明所在的坐标系,例如 d A ( t ) d t ∣ M \left . \frac{dA(t)}{dt} \right \vert_M dtdA(t) M 或者 d A ( t ) d t ∣ L \left . \frac{dA(t)}{dt} \right \vert_L dtdA(t) L。 前面式子中 L L L 表示实验室坐标系(静止坐标系), M M M 表示运动坐标系。
d A ( t ) d t ∣ L = d A ( t ) d t ∣ M + A 1 ′ d e ^ 1 ′ d t + A 2 ′ d e ^ 2 ′ d t + A 3 ′ d e ^ 3 ′ d t \left . \frac{d \mathbf A(t)}{dt} \right \vert_L = \left . \frac{d \mathbf A(t)}{dt} \right \vert_M + A_1' \frac{d \mathbf{\hat e_1'}}{dt} + A_2' \frac{d \mathbf{\hat e_2'}}{dt} + A_3' \frac{d \mathbf{\hat e_3'}}{dt} dtdA(t) L=dtdA(t) M+A1′dtde^1′+A2′dtde^2′+A3′dtde^3′
单位坐标向量求导后与原坐标向量还是有联系的:
( e ^ ˙ 1 ′ e ^ ˙ 2 ′ e ^ ˙ 3 ′ ) = ( 0 a 12 a 13 a 21 0 a 23 a 31 a 32 0 ) ( e ^ 1 ′ e ^ 2 ′ e ^ 3 ′ ) \begin{pmatrix} \mathbf{\dot{\hat e}_1'} \\ \mathbf{\dot{\hat e}_2'} \\ \mathbf{\dot{\hat e}_3'} \end{pmatrix} = \begin{pmatrix} 0 & a_{12} &a_{13}\\ a_{21} & 0 & a_{23}\\ a_{31} & a_{32} & 0 \end{pmatrix} \begin{pmatrix} \mathbf{\hat e}_1' \\ \mathbf{\hat e}_2' \\ \mathbf{\hat e}_3' \end{pmatrix} e^˙1′e^˙2′e^˙3′ = 0a21a31a120a32a13a230 e^1′e^2′e^3′
上面的变换矩阵是反对称矩阵,因为单位向量互相有正交关系:
e ^ m ′ ⋅ e ^ n ′ = 0 , m ≠ n d ( e ^ m ′ ⋅ e ^ n ′ ) d t = e ^ ˙ m ′ e ^ n ′ + e ^ m ′ e ^ ˙ n ′ = 0 \mathbf{\hat e}_m' \cdot \mathbf{\hat e}_n' = 0 , m \neq n\\ \frac{d (\mathbf{\hat e}_m' \cdot \mathbf{\hat e}_n')}{dt} = \mathbf{\dot{\hat e}}_m' \mathbf{\hat e}_n' + \mathbf{\hat e}_m' \mathbf{ \dot {\hat e}}_n' = 0 e^m′⋅e^n′=0,m=ndtd(e^m′⋅e^n′)=e^˙m′e^n′+e^m′e^˙n′=0
所有有如下关系:
d A ( t ) d t ∣ L = d A ( t ) d t ∣ M + A 1 ′ d e ^ 1 ′ d t + A 2 ′ d e ^ 2 ′ d t + A 3 ′ d e ^ 3 ′ d t = d A ( t ) d t ∣ M + A 1 ′ ( a 12 e ^ 2 ′ + a 13 e ^ 3 ′ ) + A 2 ′ ( a 21 e ^ 1 ′ + a 23 e ^ 3 ′ ) + A 3 ′ ( a 31 e ^ 1 ′ + a 32 e ^ 2 ′ ) = d A ( t ) d t ∣ M + ( A 2 ′ a 21 + A 3 ′ a 31 ) e ^ 1 ′ + ( A 1 ′ a 12 + A 3 ′ a 32 ) e ^ 2 ′ + ( A 1 ′ a 13 + A 2 ′ a 23 ) e ^ 3 ′ = d A ( t ) d t ∣ M + ( − A 2 ′ a 12 − A 3 ′ a 13 ) e ^ 1 ′ + ( A 1 ′ a 12 − A 3 ′ a 23 ) e ^ 2 ′ + ( A 1 ′ a 13 + A 2 ′ a 23 ) e ^ 3 ′ = d A ( t ) d t ∣ M + ( e ^ 1 ′ e ^ 2 ′ e ^ 3 ′ a 23 − a 13 a 12 A 1 ′ A 2 ′ A 3 ′ ) = d A ( t ) d t ∣ M + ω × A \begin{equation} \begin{aligned} \left . \frac{d \mathbf A(t)}{dt} \right \vert_L &= \left . \frac{d \mathbf A(t)}{dt} \right \vert_M + A_1' \frac{d \mathbf{\hat e_1'}}{dt} + A_2' \frac{d \mathbf{\hat e_2'}}{dt} + A_3' \frac{d \mathbf{\hat e_3'}}{dt} \\ &= \left . \frac{d \mathbf A(t)}{dt} \right \vert_M + A_1' (a_{12} \mathbf{\hat e_2'} + a_{13} \mathbf{\hat e_3'}) + A_2'(a_{21} \mathbf{\hat e_1'} + a_{23} \mathbf{\hat e_3'}) + A_3' (a_{31} \mathbf{\hat e_1'} + a_{32} \mathbf{\hat e_2'}) \\ &= \left . \frac{d \mathbf A(t)}{dt} \right \vert_M + (A_2' a_{21} + A_3' a_{31}) \mathbf{\hat e_1'} + (A_1' a_{12} + A_3' a_{32}) \mathbf{\hat e_2'} + (A_1' a_{13} + A_2' a_{23}) \mathbf{\hat e_3'} \\ &= \left . \frac{d \mathbf A(t)}{dt} \right \vert_M+ (- A_2' a_{12} - A_3' a_{13}) \mathbf{\hat e_1'} + (A_1' a_{12} - A_3' a_{23}) \mathbf{\hat e_2'} + (A_1' a_{13} + A_2' a_{23}) \mathbf{\hat e_3'} \\ &= \left . \frac{d \mathbf A(t)}{dt} \right \vert_M + \begin{pmatrix} \mathbf{\hat e_1'} & \mathbf{\hat e_2'} & \mathbf{\hat e_3'} \\ a_{23} & -a_{13} & a_{12}\\ A_1' & A_2' & A_3' \end{pmatrix} \\ &= \left . \frac{d \mathbf A(t)}{dt} \right \vert_M + \mathbf {\omega} \times \mathbf A \end{aligned} \end{equation} dtdA(t) L=dtdA(t) M+A1′dtde^1′+A2′dtde^2′+A3′dtde^3′=dtdA(t) M+A1′(a12e^2′+a13e^3′)+A2′(a21e^1′+a23e^3′)+A3′(a31e^1′+a32e^2′)=dtdA(t) M+(A2′a21+A3′a31)e^1′+(A1′a12+A3′a32)e^2′+(A1′a13+A2′a23)e^3′=dtdA(t) M+(−A2′a12−A3′a13)e^1′+(A1′a12−A3′a23)e^2′+(A1′a13+A2′a23)e^3′=dtdA(t) M+ e^1′a23A1′e^2′−a13A2′e^3′a12A3′ =dtdA(t) M+ω×A
上面推导过程中定义了一个新的量: ω = ( a 23 − a 13 a 12 ) \omega = \begin{pmatrix} a_{23} & -a_{13} & a_{12}\end{pmatrix} ω=(a23−a13a12) 。 我们研究一下这个 ω \omega ω 的物理意义。
a 23 a_{23} a23 对应的是单位矢量 e ^ 2 ′ \mathbf{\hat e}_2' e^2′ 运动时速度矢量在 e ^ 3 ′ \mathbf{\hat e}_3' e^3′ 方向上的分量。以我们常见的 x , y , z x,y,z x,y,z 坐标系为例,这个就是 y y y 轴上单位矢量向上运动的速度分量。这个速度其实就是沿 x x x 轴逆时针旋转角速度。所以 a 23 a_{23} a23 就是 x x x 方向的角速度的值。
a 13 a_{13} a13 对应的是单位矢量 e ^ 1 ′ \mathbf{\hat e}_1' e^1′ 运动时速度矢量在 e ^ 3 ′ \mathbf{\hat e}_3' e^3′ 方向上的分量。这个就是 x x x 轴上单位矢量向上运动的分量。这个运动其实就是沿 y y y 轴顺时针旋转角速度。也就是 y y y 方向角速度的负值。
类似的分析过程 a 12 a_{12} a12 就是 z z z 方向角动量值。 所以 ω \omega ω 就是角速度矢量 :
ω = ( a 23 − a 13 a 12 ) = ( ω x ω y ω z ) \omega=\begin{pmatrix} a_{23} & -a_{13} & a_{12}\end{pmatrix} = \begin{pmatrix} \omega_x & \omega_y & \omega_z \end{pmatrix} ω=(a23−a13a12)=(ωxωyωz)
其实我们也可以把 ω = ( a 23 − a 13 a 12 ) \omega = \begin{pmatrix} a_{23} & -a_{13} & a_{12}\end{pmatrix} ω=(a23−a13a12) 作为角速度矢量的定义式。当然这样定义后需要我们去验证这个矢量是否满足矢量的运算律,也就是是否满足矢量加法的平行四边形法则。
这里还可以多说几句,角速度是矢量这个结论不是显而易见的。角位移不是矢量,不满足加法的交换律。无限小的角位移可以定义为矢量,因为无限小的角位移满足矢量加法定律。角速度是与无限小角位移相关联的,所以也是矢量。
Introduction of the operator D ^ \widehat D D
定义两个算符,分别表示在两个坐标系下的求导操作:
D ^ L = ∂ ∂ t ∣ L D ^ M = ∂ ∂ t ∣ M \widehat D_L = \left .\frac{\partial}{\partial t} \right |_L \\ \widehat D_M = \left .\frac{\partial}{\partial t} \right |_M D L=∂t∂ LD M=∂t∂ M
则有:
D ^ L A = D ^ M A + ω × A \widehat D_L \mathbf A = \widehat D_M \mathbf A + \omega \times \mathbf A D LA=D MA+ω×A
注意这里的 A \mathbf A A 是任意的矢量。如果 r \mathbf r r 表示位置矢量。那么有:
a ∣ L = r ¨ ∣ L = D ^ L r ˙ = D ^ L ( D ^ L r ) = ( D ^ M + ω × ) ( D ^ M r + ω × r ) = D ^ M 2 r + D ^ M ( ω × r ) + ω × ( D ^ M r ) + ω × ω × r = D ^ M 2 r + ( D ^ M ω ) × r + 2 ω × ( D ^ M r ) + ω × ω × r \mathbf a |_L = \ddot{ \mathbf{r}} |_L= \widehat D_L \dot{\mathbf{r}}=\widehat D_L{(\widehat D_L \mathbf r)} = (\widehat D_M+\omega \times) {(\widehat D_M \mathbf r + \omega \times \mathbf r)} \\ = \widehat D_M^2 \mathbf r + \widehat D_M (\omega \times \mathbf r) + \omega \times (\widehat D_M \mathbf r)+ \omega \times \omega \times \mathbf r \\ = \widehat D_M^2 \mathbf r + (\widehat D_M \omega) \times \mathbf r +2 \omega \times (\widehat D_M \mathbf r)+ \omega \times \omega \times \mathbf r a∣L=r¨∣L=D Lr˙=D L(D Lr)=(D M+ω×)(D Mr+ω×r)=D M2r+D M(ω×r)+ω×(D Mr)+ω×ω×r=D M2r+(D Mω)×r+2ω×(D Mr)+ω×ω×r
Formulation of Newton’s equation in the rotating coordinate system
牛顿第二定律可以表示为:
F = m a ∣ L = m r ¨ ∣ M + ω ˙ ∣ M × r + 2 ω × r ˙ ∣ M + ω × ω × r m r ¨ ∣ M = F − ω ˙ ∣ M × r − 2 ω × r ˙ ∣ M − ω × ω × r \mathbf F = m \ \mathbf a|_L = m \ \ddot{\mathbf r}|_M + \dot{\omega}|_M \times \mathbf r + 2 \omega \times \dot{\mathbf r}|_M + \omega \times \omega \times \mathbf r \\ m \ \ddot{\mathbf r}|_M = \mathbf F - \dot{\omega}|_M \times \mathbf r - 2 \omega \times \dot{\mathbf r}|_M - \omega \times \omega \times \mathbf r F=m a∣L=m r¨∣M+ω˙∣M×r+2ω×r˙∣M+ω×ω×rm r¨∣M=F−ω˙∣M×r−2ω×r˙∣M−ω×ω×r
需要注意的是所有的速度和加速度都要标明是在哪个坐标系,但是位置矢量是与坐标系无关的。
Newton’s equations in systems with arbitrary relative motion
上面我们说过一个矢量 A \mathbf A A 在不同坐标系下表现形式不同,但是都是同一个矢量。其实这个说法是不完全准确的。只有当两个坐标系的原点是重合的时候上述结论是正确的。之所以会这样就是因为单位矢量本身不含有这个矢量的起点所在位置的任何信息。
设 M 系的原点在 L系的坐标为 R \mathbf R R,这里的 M 系的原点指的是 M 系的旋转中心。那么 M 系的任意点 r ’ \mathbf r’ r’ 在 L 系中的位置矢量为 r = R + r ′ \mathbf r = \mathbf R + \mathbf r' r=R+r′ 。 注意式子中的矢量都是位置矢量,与坐标系是无关的。带入牛顿第二定律。
F = m r ¨ ∣ L = m R ¨ ∣ L + m r ′ ¨ ∣ L = m R ¨ ∣ L + m r ′ ¨ ∣ M + ω ˙ ∣ M × r ′ + 2 ω × r ′ ˙ ∣ M + ω × ω × r ′ \begin{aligned} \mathbf F &= m \ \ddot{\mathbf r}|_L \\ &= m \ \ddot{\mathbf R}|_L + m \ \ddot{\mathbf r'}|_L \\ &= m \ \ddot{\mathbf R}|_L + m \ \ddot{\mathbf r'}|_M + \dot{\omega}|_M \times \mathbf r' + 2 \omega \times \dot{\mathbf r'}|_M + \omega \times \omega \times \mathbf r' \end{aligned} F=m r¨∣L=m R¨∣L+m r′¨∣L=m R¨∣L+m r′¨∣M+ω˙∣M×r′+2ω×r′˙∣M+ω×ω×r′
整理一下,可以得到:
m r ′ ¨ ∣ M = F − m R ¨ ∣ L − ω ˙ ∣ M × r ′ − 2 ω × r ′ ˙ ∣ M − ω × ω × r ′ m \ \ddot{\mathbf r'}|_M = \mathbf F - m \ \ddot{\mathbf R}|_L - \dot{\omega}|_M \times \mathbf r' - 2 \omega \times \dot{\mathbf r'}|_M - \omega \times \omega \times \mathbf r' m r′¨∣M=F−m R¨∣L−ω˙∣M×r′−2ω×r′˙∣M−ω×ω×r′