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

封面图：【荷】伦勃朗 (1606–1669)，《天使拉斐尔离开托比特家族》，1637年，收藏于巴黎卢浮宫

本文先讨论开口系的动量守恒定律，重点介绍有强间断特性的冲击波波阵面上的动量守恒，然后过渡至连续波情形。波阵面上动量守恒的物理意义可参考 波阵面上动量守恒的物理意义。

1. 开口系的动量守恒定律

开口系的动量守恒定律表述为：

在任一时刻，开口系体系的动量变化率 = 该时刻体系所受外力矢量和 + 外界向体系的动量纯流入率

考虑初始面积为$A_0$、初始质量为${\rho }_0$的一维杆中的纵波传播，在杆中取一个处于两截面之间的体系，这两个截面的拉格朗日坐标分别为$X_1(t)$ 和 $X_2(t)$，只要$X_1(t)$和$X_2(t)$随时间变化，此体系就是开口系；否则是闭口系。

考虑一个开口系，记为$(X_1(t),X_2(t))$，此开口系的动量守恒写为：

$\frac{\mathrm{d}}{\mathrm{d}}{\int }_{X_1(t)}^{X_2(t)}{\rho }_0{A}_{0}v(X,t)\mathrm{d}X=A_0({\sigma }_2-{\sigma }_1)+{\rho }_0{A}_0\frac{\mathrm{d}{X}_2}{\mathrm{d}t}{v}_2-{\rho }_0{A}_0\frac{\mathrm{d}{X}_2}{\mathrm{d}t}{v}_1(1)$

当取$X_1(t)=X_2(t)=X_{\ast }(t)$时，恰好为冲击波的波阵面，当此开口系为附着在波阵面上的一个无限薄的薄层，有：

$\frac{\mathrm{d}{X}_1(t)}{\mathrm{d}t}=\frac{\mathrm{d}{X}_2(t)}{\mathrm{d}t}=\frac{\mathrm{d}{X}_{\ast }(t)}{\mathrm{d}t}=C(2)$

式(2)恰好为波阵面的拉格朗日波速，记为 $C$，且动量守恒方程式(1)左端为零。在贴近波阵面的紧前方和紧后方，介质应力和质点速度的跳跃量分别记为 ${\sigma }_1-{\sigma }_2=[\sigma ]$ 和 $v_1-v_2-[v]$，于是式(1)写为：

$[\sigma ]=-{\rho }_{0}C[v](3)$

式(3)成立的基础是附着在波阵面上的无限薄层动量守恒，也称为波阵面上的动量守恒条件，它将跨越波阵面时的质点速度跳跃量和应力跳跃量联系起来。类似地，以上结论也可由闭口系的动量守恒方程导出，此处略。

2. 冲击波和连续波的分类

为了继续讨论，这里需要区分冲击波和连续波的概念。在数学上，波阵面是一个奇异面（singular surface），当跨过这个奇异面时介质中的某些物理量发生间断。
把受扰动部分和未受扰动部分的分界面（或新被扰动的部分和旧的被扰动的部分），称为波阵面（wave front）。

在连续介质力学中，我们考虑无裂纹的情况，即保证位移是连续的。区分如下情形：
（1）把跨过波阵面位移连续、但一阶导数发生间断的波阵面称为一阶奇异面，也称为强间断波（strong discontinuity wave）、冲击波（shock wave），此处位移的一阶导数主要指质点速度、应变（与之相连的应力）；

（2）把跨过波阵面位移、位移的一阶导数连续，但是位移的二阶导数间断的波阵面称为二阶奇异面，也称为弱间断波（weak discontinuity wave）、连续波（continuity wave）、加速度波，此处位移的二阶导数主要指加速度、应变梯度、应变率。

3. 一个注记：应力波波速

应力波在介质中的传播速度称为波速（wave speed），特别地，波速是指应力扰动信号的传播速度，这个概念与波扰动引起的介质本身质点的质点速度是完全不同的。一般后者远小于前者，且二者方向可能不同（李永池，2018）。

4. 应力波波速的定义

为了对右/左行波的波速表述形式统一，将$C$定义为波速的绝对值，对于左/右行波均有以下重新定义的波速（基于式(3)）：

$[\sigma ]={\rho }_0{C}^2[\epsilon ](4)$

式中 $\epsilon$ 为工程应变，表示为 $\epsilon =(\partial u/\partial X){|}_t$，下标表示保持时间不变求导。于是可反解$C$：

$C=\sqrt{\frac{[\sigma ]}{{\rho }_0[\epsilon ]}}(5)$

5. 连续波的波阵面动量守恒条件

以上讨论是针对冲击波而言的，即强间断波，对于连续波（弱间断波），可将其视为由无穷多个无限小的增量波依次构成的，则跨过每一个增量波的扰动量可由其微分表达，于是式(3)(4)写为：

$\mathrm{d}\sigma =\pm {\rho }_{0}C\mathrm{d}v(6\mathrm{a})$
$C=\sqrt{\frac{\mathrm{d}\sigma }{{\rho }_0\mathrm{d}\epsilon }}(6\mathrm{b})$

式中正、负号分别表示左、右行波。

参考资料

• 波阵面上动量守恒的物理意义
• Achenbach, J., 1973. Wave propagation in elastic solids.
• North-Holland Publishing Company.
• 李永池, 2018. 波动力学（第二版）. 中国科学技术大学出版社.
• 王礼立. 2023. 应力波基础（第三版）. 国防工业出版社.