动量守恒定律与运动方程
在固体力学中，研究物体在力作用下的变形与运动是核心目标。而所有力学分析的基石在于物理量的守恒定律。动量守恒定律（Law of Conservation of Momentum）不仅是经典力学的基本定律，更是推导连续介质运动方程的理论源头。通过将牛顿第二定律扩展到连续体，我们可以建立起描述固体内部应力分布与运动状态之间关系的数学模型。
动量守恒定律的基本概念
动量守恒定律在宏观物理层面定义为：一个系统在不受外力作用时，其总动量保持不变。对于一个质量为 $m$、速度为 $\mathbf{v}$ 的质点，其动量定义为 $\mathbf{P} = m\mathbf{v}$。
在固体力学研究的连续介质中，我们不再关注单个质点，而是关注一个特定的物质区域 $\Omega$。该区域的总动量 $\mathbf{P}$ 可表示为体积积分：
$$\mathbf{P} = \int_{\Omega} \rho \mathbf{v} , d\Omega$$
其中，$\rho$ 为材料密度，$\mathbf{v}$ 为速度矢量。根据牛顿第二定律，该区域动量的随时间变化率等于作用在物体上的所有外力的矢量和：
$$\frac{d\mathbf{P}}{dt} = \mathbf{F}_{ext}$$
运动方程的推导过程
为了将宏观的守恒定律转化为可用于计算的局部微分方程，我们需要分析作用在连续体上的外力组成。外力通常分为两类：
- 表面力 (Surface Forces)：作用在物体边界 $\Gamma$ 上的力，如压力、摩擦力。其单位面积上的力称为牵引力 $\mathbf{t}$。
- 体力 (Body Forces)：作用在物体内部每个微元上的力，如重力、电磁力。其单位体积的力记为 $\mathbf{b}$。
1. 力的积分表达式
总外力 $\mathbf{F}{ext}$ 可以表示为：
$$\mathbf{F}{ext} = \int_{\Gamma} \mathbf{t} , d\Gamma + \int_{\Omega} \rho \mathbf{b} , d\Omega$$
2. 引入柯西应力原理
根据柯西应力原理，表面牵引力 $\mathbf{t}$ 与内部应力张量 $\boldsymbol{\sigma}$ 的关系为：
$$\mathbf{t} = \boldsymbol{\sigma} \cdot \mathbf{n}$$
其中 $\mathbf{n}$ 是边界的单位外法矢量。将此代入表面力积分项，并利用高斯散度定理（Gauss's Divergence Theorem）将面积分转化为体积分：
$$\int_{\Gamma} \sigma_{ij} n_j , d\Gamma = \int_{\Omega} \frac{\partial \sigma_{ij}}{\partial x_j} , d\Omega$$
3. 建立局部平衡方程
将上述各项代入动量变化率方程，并假设控制体 $\Omega$ 为一个随物质运动的材质体，可得：
$$\int_{\Omega} \rho \frac{d\mathbf{v}}{dt} , d\Omega = \int_{\Omega} \left( \nabla \cdot \boldsymbol{\sigma} + \rho \mathbf{b} \right) , d\Omega$$
由于该等式对任意体积 $\Omega$ 均成立，因此被积函数必须相等。由此得到了连续介质的运动方程（Cauchy's Equation of Motion）：
$$\rho \frac{\partial v_i}{\partial t} + \rho v_j \frac{\partial v_i}{\partial x_j} = \frac{\partial \sigma_{ij}}{\partial x_j} + \rho b_i$$
在小变形且低速运动的假设下，对流项 $\rho v_j \frac{\partial v_i}{\partial x_j}$ 可忽略不计，方程简化为：
$$\rho \frac{\partial^2 u_i}{\partial t^2} = \sigma_{ij,j} + \rho b_i$$
其中 $u_i$ 为位移矢量，$\sigma_{ij,j}$ 表示应力张量对坐标的偏导数之和。
运动方程的物理意义
运动方程揭示了固体内部三种力的动态平衡关系：
- 惯性力项 ($\rho \ddot{u}_i$)：代表物体由于加速运动而产生的抗力。
- 内力项 ($\sigma_{ij,j}$)：代表应力梯度产生的合力。当物体内部应力分布不均时，会产生驱动物体运动的内力。
- 体力项 ($\rho b_i$)：代表外部场（如重力场）对物质微元的直接作用。
特殊情况分析
1. 静态平衡 (Static Equilibrium)
当物体处于静止状态或匀速直线运动状态时，加速度 $\frac{\partial^2 u_i}{\partial t^2} = 0$。此时运动方程退化为平衡方程：
$$\frac{\partial \sigma_{ij}}{\partial x_j} + \rho b_i = 0$$
这是结构力学和静力学分析的基础。
2. 波动传播 (Wave Propagation)
在动态分析中，若考虑弹性材料（$\sigma = C : \varepsilon$），运动方程将演变为波动方程。例如，在简单的 1D 弹性杆中，方程形式为：
$$\rho \frac{\partial^2 u}{\partial t^2} = E \frac{\partial^2 u}{\partial x^2}$$
这描述了应力波在固体中的传播特性。
实际应用示例：一维受压杆的动态响应
假设有一根长度为 $L$、截面积为 $A$ 的弹性杆，一端固定，另一端受到瞬时冲击力 $F(t)$。
- 建立模型：忽略重力（$\rho b_i = 0$），仅考虑轴向 $x$ 方向。
- 控制方程：$\rho \frac{\partial^2 u}{\partial t^2} = \frac{\partial \sigma_{xx}}{\partial x}$。
- 本构关系：$\sigma_{xx} = E \frac{\partial u}{\partial x}$。
- 组合方程：$\rho \frac{\partial^2 u}{\partial t^2} = E \frac{\partial^2 u}{\partial x^2}$。
通过求解该偏微分方程，我们可以计算出冲击波在杆中的传播速度 $c = \sqrt{E/\rho}$ 以及杆内各点随时间变化的位移和应力分布。
总结
动量守恒定律是固体力学从物理定律向数学模型跨越的桥梁。通过从全局积分形式推导出局部的柯西运动方程，我们能够将复杂的物体运动问题转化为对应力张量和位移场的偏微分方程求解。无论是简单的静力分析，还是复杂的冲击动力学，其核心逻辑均是对动量守恒定律的具体应用。