应力波波阵面上的能量守恒条件

本文介绍波阵面上的能量守恒条件，并导出动能和内能的表达式。可以根据文章应力波波阵面上的动量守恒条件
中关于动量守恒方程的推导过程，类似地导出开放系波阵面的能量守恒方程。此处采用另一种方法，即闭口系的能量守恒。这种处理可以加深对物理意义的理解。

1. 系统的能量守恒定律

闭口系的能量守恒方程为：
$\mathrm{d}U+\mathrm{d}K=\mathrm{d}{W}^{\ast }+\mathrm{d}Q(1)$

式中$\mathrm{d}U$ 和 $\mathrm{d}K$分别表示任意时间间隔$\mathrm{d}t$内的闭口体系内能增量和动能增量；$\mathrm{d}{W}^{\ast }$ 和 $\mathrm{d}Q$分别表示外部在时间间隔$\mathrm{d}t$内对该体系所做的功和纯供热。由于波动过程极快，因此通常认为波动过程是绝热的，热效应$\mathrm{d}Q$来不及影响波动。

对于不太剧烈的连续波波动过程，可视为绝热等熵（可逆的）；对于剧烈的冲击波波动过程，可视为绝热熵增过程（不可逆的），也就是绝热冲击波。外部供热为零，则式(1)：
$\mathrm{d}U+\mathrm{d}K=\mathrm{d}{W}^{\ast }(2)$

2. 波阵面上的能量守恒条件

初始面积为$A_0$、初始质量为${\rho }_0$的一维杆中有一个拉格朗日波速为$C$的冲击波，拉格朗日波阵面在$t$ 和 $t+\mathrm{d}t$时刻分别到达拉格朗日坐标为$X$ 和 $X+C\mathrm{d}t$ 的位置。

将$\mathrm{d}t$时间段内冲击波扫过的一段杆介质视为闭口系，其质量为${\rho }_0{A}_{0}C\mathrm{d}t$。冲击波的波阵面紧前方和紧后方的物理量记为 ${\sigma }^{+},{\epsilon }^{+},v^{+},e^{+}$ 和 ${\sigma }^{-},{\epsilon }^{-},v^{-},e^{-}$，分别表示工程应力、工程应变、质点速度、比内能。

则该闭口体系的绝热条件下能量守恒条件为：

${\rho }_0{A}_{0}C\mathrm{d}t\left[\frac{(v^{-}{)}^2}{2}+e^{-}-\frac{(v^{+}{)}^2}{2}-e^{+}\right]=A_0{\sigma }^{+}{v}^{+}\mathrm{d}t-A_0{\sigma }^{-}{v}^{-}\mathrm{d}t(3\mathrm{a})$

约去$A_0$ 和 $\mathrm{d}t$：

${\rho }_{0}C\left[\frac{(v^{-}{)}^2}{2}+e^{-}-\frac{(v^{+}{)}^2}{2}-e^{+}\right]={\sigma }^{+}{v}^{+}-{\sigma }^{-}{v}^{-}$

整理得：

$\frac{1}{2}{\rho }_{0}C[v^2]+{\rho }_{0}C[e]=-[\sigma v](3\mathrm{b})$

式(3b)可视为极限意义下得单位时间内单位面积波阵面扫过得闭口体系能量守恒条件。符号 $[\phi ]={\phi }^{-}-{\phi }^{+}$ 表示由前方至后方跨过波阵面的物理量 $\phi$ 跳跃量。

而 $-[\sigma v]$ 恰好表示外力（面力）对此闭口系的功率，它可进一步分解为：

$-[\sigma v]=-\frac{1}{2}(v^{+}+v^{-})[\sigma ]-\frac{1}{2}({\sigma }^{+}+{\sigma }^{-})[v](4)$

式(4)的物理意义为，外力（面力）对闭口系的功率可分解为两项，$\frac{1}{2}(v^{+}+v^{-})$ 和 $\frac{1}{2}({\sigma }^{+}+{\sigma }^{-})$ 可分别理解为平均速度和平均应力。

3. 波阵面上的动能定理和内能

借助波阵面的动量守恒条件 $[\sigma ]=-{\rho }_{0}C[v]$ （参考 应力波波阵面上的动量守恒条件中的式(3)），结合式(3b)和(4)，可得：

$-\frac{1}{2}(v^{+}+v^{-})[\sigma ]=\frac{1}{2}{\rho }_{0}C[v^2](5\mathrm{a})$
$-\frac{1}{2}({\sigma }^{+}+{\sigma }^{-})[v]={\rho }_{0}C[e](5\mathrm{b})$

式(5a)和(5b)分别表示动能定理和内能，物理意义为：

（1）闭口系的动能变化率等于波阵面前后方的不均衡应力在波阵面前后方平均速度上的功率；

（2）闭口系的内能变化率等于波阵面前后方的平均应力在波阵面前后方速度差上的功率。即介质的内能就是应变能。

式(5a)和(5b)实际上是能量守恒式(3b)的两个分解式。

参考资料

• 应力波波阵面上的动量守恒条件
• 波阵面上动量守恒的物理意义