We re-examine the conservation of energy and momentum we learned in middle school from a purely mathematical perspective, and then we derive the conservation of angular momentum from the same perspective.
我们从纯数学的角度来重新审视中学学过的能量守恒和动量守恒,然后我们再通过同样的角度得出角动量守恒。
守恒定律是什么⌗ 能量守恒,动量守恒和机械能守恒都是运动积分。 为了明白什么是运动积分,需要弄清楚下面两个问题。
动能定理 和 动量定理 与牛顿第二定律 的关系动能定理 和动量守恒 以及 动量定理 和动量守恒 的关系。前置知识:⌗ Nabla算子 ∇
\nabla
∇
标量场Ψ ( r ⃗ ) \Psi(\vec{r}) Ψ ( r ) 的梯度 是向量 :
∇ Ψ = ∂ Ψ ∂ x i ⃗ + ∂ Ψ ∂ y j ⃗ + ∂ Ψ ∂ z k ⃗ \nabla\Psi = \frac{\partial{\Psi}}{\partial{x}}\vec{i} + \frac{\partial{\Psi}}{\partial{y}}\vec{j}+ \frac{\partial{\Psi}}{\partial{z}}\vec{k} ∇ Ψ = ∂ x ∂ Ψ i + ∂ y ∂ Ψ j + ∂ z ∂ Ψ k
矢量场A ⃗ ( r ⃗ ) \vec{A}(\vec{r}) A ( r ) 的散度 是标量 :
散度:
∇ ⋅ A ⃗ = ∂ A x ∂ x + ∂ A y ∂ y + ∂ A z ∂ z
\nabla \cdot \vec{A}=\frac{\partial{A_x}}{\partial{x}} + \frac{\partial{A_y}}{\partial{y}}+ \frac{\partial{A_z}}{\partial{z}}
∇ ⋅ A = ∂ x ∂ A x + ∂ y ∂ A y + ∂ z ∂ A z
旋度:
∇ × A ⃗
\nabla \times \vec{A}
∇ × A
结论:
∇ × ( ∇ Ψ ) = 0
\nabla\times\left(\nabla\Psi\right) = 0
∇ × ( ∇ Ψ ) = 0
标量场的梯度的旋度为0
∇ ⋅ ( ∇ × A ⃗ ) = 0
\nabla\cdot\left(\nabla\times\vec{A}\right) = 0
∇ ⋅ ( ∇ × A ) = 0
矢量场的旋度的散度为0
这两个反命题也成立,也就是说,一个旋度为0的(向量)场,可以表示为一个标量场的梯度.
一个散度为0的(向量)场,可以表示为一个(向量)场的旋度。
例子:
对于一个无旋场 ,我们可以用一个标量场 来表示,这个标量场一般称为 “势” 。比如静电场(标量场)的电势场(只是和位置有关的一个数)和 重力场(标量场)的重力势(只是一个和位置有关的标量) 。 同时,我们把这种旋度为0 的力场称为保守力 。
能量守恒⌗ F ⃗ = m r ⃗ ¨ ⇒ F ⃗ = m d v ⃗ d t
\vec{F} = m\ddot{\vec{r}} \Rightarrow \vec{F} = m \frac{d\vec{v}}{dt}
F = m r ¨ ⇒ F = m d t d v
如果我们定义一个量
T = 1 2 m v 2 T = \frac{1}{2}mv^2 T = 2 1 m v 2
那么就有
F ⃗ d r ⃗ = m d v ⃗ d t d r ⃗ = d ( 1 2 m v 2 )
\vec{F}d\vec{r} = m\frac{d\vec{v}}{dt}d\vec{r} = d(\frac{1}{2}mv^2)
F d r = m d t d v d r = d ( 2 1 m v 2 )
得到动能定理
F ⃗ d r ⃗ = d T \vec{F}d\vec{r} = dT F d r = d T
其积分形式
T 2 − T 1 = ∫ r 1 ⃗ r 2 ⃗ F ⃗ d r ⃗
T_2-T_1 = \int_{\vec{r_1}}^{\vec{r_2}}\vec{F}d\vec{r}
T 2 − T 1 = ∫ r 1 r 2 F d r
如果F ⃗ = 0 \vec{F} = 0 F = 0 ,运动前后动能不变 .
动量守恒⌗ F ⃗ = m r ⃗ ¨ ⇒ F ⃗ = m d v ⃗ d t
\vec{F} = m\ddot{\vec{r}} \Rightarrow \vec{F} = m \frac{d\vec{v}}{dt}
F = m r ¨ ⇒ F = m d t d v
如果我们定义一个量p ⃗ = m v ⃗ \vec{p} = m\vec{v} p = m v
就得到
F ⃗ = d p ⃗ d t
\vec{F} = \frac{d\vec{p}}{dt}
F = d t d p
首先我们需要明白,动量定力和牛顿第二定律F ⃗ = m r ⃗ ¨ \vec{F} = m\ddot{\vec{r}} F = m r ¨ 都是二阶微分方程,也就是他们的次数是一样的,从形式上来讲他们是等价的。实际上,牛顿最开始提出来的第二定律其实就是这个我们今天所熟悉的动量定理 F ⃗ = d p ⃗ d t \vec{F} = \frac{d\vec{p}}{dt} F = d t d p
然后我们把这个动量定理对时间 积分,得到
p 2 ⃗ − p 1 ⃗ = ∫ t 1 t 2 F ⃗ d t
\vec{p_2} - \vec{p_1} = \int_{t_1}^{t_2}\vec{F}dt
p 2 − p 1 = ∫ t 1 t 2 F d t
如果F ⃗ = 0 \vec{F} = 0 F = 0 ,运动前后动量不变 , 因为和牛顿第二定律等价,所以也可以从牛顿第二定律的角度来解释,如果合外力为0,那么速度不变 。
当然以上都是显而易见的结论,但重点在下面:
也就是说
F ⃗ = 0 ⇒ p ⃗ = C ⃗
\vec{F} = 0 \Rightarrow \vec{p} = \vec{C}
F = 0 ⇒ p = C
其中的这个 p ⃗ = C ⃗ \vec{p}= \vec{C} p = C 是一个守恒量 。
我们知道,p ⃗ \vec{p} p 里面的速度是r ⃗ \vec{r} r 的一阶导数,因此,如果我们以计算出系统的位置r ⃗ \vec{r} r 为目标,那么这个守恒量p ⃗ = C ⃗ \vec{p} = \vec{C} p = C 就是一个一阶微分方程,相比于牛顿第二定律或者说动量定理,已经降了一阶了。
动量矩守恒(角动量守恒)⌗ 位置的矢量和矢量的叉积
位置矢量和力
M ⃗ = r ⃗ × F ⃗
\vec{M} = \vec{r}\times\vec{F}
M = r × F
动量矩(角动量)⌗ 位置矢量和动量
J ⃗ = r ⃗ × p ⃗
\vec{J} = \vec{r}\times\vec{p}
J = r × p
所谓的动量矩定理,就是
M ⃗ = d J ⃗ d t
\vec{M} = \frac{\mathrm{d}{\vec{J}}}{\mathrm{d}{t}}
M = d t d J
也就是动量矩对时间的导数是的力矩,求导 + 牛顿第二定律 就可以证明。
所以,所以动量矩定理和牛顿第二定律等价。
我们类比之前的动量守恒,如果这里的力矩M ⃗ = 0 \vec{M} = 0 M = 0 ,也可以得到积分后的量是一个不变量。这就是动量矩(角动量)守恒 。
动量守恒或动量矩守恒发生前提(外力为0)是在是有点苛刻。然而我们解决的问题基本上不可能不受力的,所以,对于这两种守恒,我们一般使用其分量形式来解决问题。即在某一固定方向(对动量定理来说)或对于某一固定点(对于角动量守恒)分别使用这两个守恒。在某一固定方向上应用动量定理的典型问题是一个滑块在静止的斜坯上滑下,问这两个物体的运动规律。这当然是一个在水平方向上动量守恒的例子。对于角动量守恒的经典问题是只受有心力的物体叫动量守恒,比如天体的轨道运动规律。
机械能守恒⌗ 我们看到,前面的动能定理和动量定理本质上是在某种特殊的条件下(F=0)去得到降过阶的微分方程 。 使我们不用再去解F ⃗ = m r ⃗ ¨ \vec{F}=m\ddot{\vec{r}} F = m r ¨ 这种复杂的二阶微分方程了。有了这种想法,我们自然更多的要去挖掘更多的运动积分 。 然而,我们从牛顿第二定律出发,似乎已经不能得到新的运动积分了。但是,如果我们继续施加某种条件 ,可能又会发现一些守恒量。幸运的是,我们发现,如果只有保守力 为主动力,我们又可以的到机械能守恒 这个非常重要的运动积分。
之前的前置知识里面说了,保守力可以写作一个标量场的形式,就是势能。我们无时无刻不在的保守力重力场中,有:
∇ × F ⃗ = 0 ⇒ F ⃗ = − ∇ V ( r ⃗ )
\nabla\times\vec{F} = 0 \Rightarrow \vec{F} = -\nabla V(\vec{r})
∇ × F = 0 ⇒ F = − ∇ V ( r )
再通过之前得到动能定理的方式得到:
F ⃗ ⋅ d r ⃗ = − ∇ V ( r ⃗ ) ⋅ d r ⃗ = − d V ( r ⃗ )
\vec{F}\cdot\mathrm{d}{\vec{r}} = -\nabla V({\vec{r}})\cdot\mathrm{d}{\vec{r}} = -\mathrm{d}{V(\vec{r})}
F ⋅ d r = − ∇ V ( r ) ⋅ d r = − d V ( r )
中间那一步关键的过程直接在迪卡尔坐标系下写完全就能把Nabla算子转换成一个三元函数的全微分.
因为之前我们得到的动能定理也是
F ⃗ ⋅ d r ⃗ = d T \vec{F}\cdot d\vec{r} = \mathrm{d}T F ⋅ d r = d T
因此
d T = − d V ( r ⃗ )
\mathrm{d}T = -\mathrm{d}{V(\vec{r})}
d T = − d V ( r )
即保守力的势能 和动能 相互转换。
积分的到新的运动积分为
T + V = C
T + V = C
T + V = C
这个就是机械能守恒 。
一维势能曲线⌗ 在保守力 的前提下,有一个这样的曲线。