胡克定律（广义）的推导与应用

在经典力学中，胡克定律（Hooke's Law）通常以 $F = kx$ 的形式出现，描述了弹簧在小变形范围内受力与形变之间的线性比例关系。然而，在固体力学与连续介质力学的研究中，物体不再被视为质点，而是被视为一个连续的介质。在这种情况下，

在经典力学中，胡克定律（Hooke's Law）通常以 $F = kx$ 的形式出现，描述了弹簧在小变形范围内受力与形变之间的线性比例关系。然而，在固体力学与连续介质力学的研究中，物体不再被视为质点，而是被视为一个连续的介质。在这种情况下，受力不再是单一的力，而是分布在物体内部的应力（Stress），而变形也不再是单一的位移，而是应变（Strain）。

因此，我们需要将胡克定律从一维的标量形式推广到三维的张量形式，即“广义胡克定律”。它是线性弹性力学理论的核心，建立了应力张量与应变张量之间的本构关系（Constitutive Relation）。

基础概念：应力与应变张量

在推导广义胡克定律之前，必须首先明确两个核心物理量：应力张量和应变张量。

1. 应力张量 $\sigma_{ij}$

在连续介质内部，某一点的受力状态由柯西应力张量（Cauchy Stress Tensor）描述。由于应力是一个二阶张量，在三维空间中，它由 9 个分量组成：
$$\sigma = \begin{bmatrix} \sigma_{xx} & \tau_{xy} & \tau_{xz} \ \tau_{yx} & \sigma_{yy} & \tau_{yz} \ \tau_{zx} & \tau_{zy} & \sigma_{zz} \end{bmatrix}$$
其中，$\sigma_{ii}$ 代表正应力，$\tau_{ij}$ ($i \neq j$) 代表剪应力。

2. 应变张量 $\epsilon_{ij}$

为了描述物体微小变形后的几何变化，我们引入微小应变张量（Infinitesimal Strain Tensor）：
$$\epsilon_{ij} = \frac{1}{2} \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right)$$
其中 $u_i$ 是位移分量。应变张量同样是一个二阶对称张量，具有 6 个独立分量。

广义胡克定律的数学表达

广义胡克定律的核心思想是：应力张量是应变张量的线性函数。用张量符号表示为：
$$\sigma_{ij} = C_{ijkl} \epsilon_{kl}$$
这里的 $C_{ijkl}$ 被称为弹性刚度张量（Elastic Stiffness Tensor），它是一个四阶张量。

由于 $i, j, k, l$ 每个指标都有 3 种取值（$x, y, z$），从数学上看，这个四阶张量最初拥有 $3^4 = 81$ 个分量。然而，通过物理对称性的约束，这个数量可以大幅度减少。

刚度张量的对称性简化

在实际的物理系统中，刚度张量并非拥有 81 个独立的常数，而是通过以下对称性进行简化：

  1. 应力张量的对称性（小对称性）：由于角动量守恒，应力张量必须是对称的，即 $\sigma_{ij} = \sigma_{ji}$。这要求 $C_{ijkl} = C_{jikl}$。
  2. 应变张量的对称性（小对称性）：根据应变定义的公式，应变张量天然是对称的，即 $\epsilon_{kl} = \epsilon_{lk}$。这要求 $C_{ijkl} = C_{ijlk}$。
    • 结合上述两点，分量数从 81 降至 36。
  3. 应变能密度的对称性（大对称性）：根据热力学第二定律，弹性应变能密度 $W$ 是应变的二次型：$W = \frac{1}{2} C_{ijkl} \epsilon_{ij} \epsilon_{kl}$。为了保证能量函数的唯一性，必须满足 $C_{ijkl} = C_{klij}$。
    • 结合大对称性后，分量数从 36 进一步降至 21。

因此，对于一个完全各向异性的材料（Anisotropic Material），描述其力学特性的最小独立参数数量是 21 个。

特殊情况：各向同性材料

在工程实践中，大多数金属和常见的建筑材料可以近似视为各向同性材料（Isotropic Material），即材料在各个方向上的力学性质完全相同。对于各向同性材料，刚度张量 $C_{ijkl}$ 的复杂度会大幅降低，仅由两个独立的物理常数决定。

1. 拉梅常数（Lamé Parameters）表示法

在各向同性条件下，广义胡克定律可以写成如下形式：
$$\sigma_{ij} = \lambda \delta_{ij} \epsilon_{kk} + 2\mu \epsilon_{ij}$$
其中：

  • $\delta_{ij}$ 是 Kronecker delta 符号（当 $i=j$ 时为 1，否则为 0）。
  • $\epsilon_{kk} = \epsilon_{xx} + \epsilon_{yy} + \epsilon_{zz}$ 是应变张量的迹（即体应变）。
  • $\lambda$ 是第一个拉梅常数。
  • $\mu$ 是第二个拉梅常数，也称为剪切模量（Shear Modulus）。

2. 与工程常数的转换

在工程领域，我们更习惯使用杨氏模量（Young's Modulus, $E$）和泊松比（Poisson's Ratio, $\nu$）。它们与拉梅常数的关系如下：
$$\lambda = \frac{E\nu}{(1+\nu)(1-2\nu)}$$
$$\mu = \frac{E}{2(1+\nu)}$$

通过这种转换，我们可以将复杂的张量运算简化为熟悉的标量运算。例如，在三向应力状态下，各向同性材料的轴向应力公式为：
$$\epsilon_{x} = \frac{1}{E} [\sigma_{x} - \nu(\sigma_{y} + \sigma_{z})]$$

实际应用示例

示例：承受三向压力的立方体块

假设一个立方体材料块受到三向均匀压力 $\sigma_x = \sigma_y = \sigma_z = -P$（压力为负），且没有剪应力。求该材料产生的体积应变 $\epsilon_v$。

解题步骤：

  1. 确定应变分量：
    根据各向同性胡克定律：
    $$\epsilon_x = \frac{1}{E} [\sigma_x - \nu(\sigma_y + \sigma_z)] = \frac{1}{E} [-P - \nu(-P - P)] = -\frac{P}{E}(1 - 2\nu)$$
    同理，$\epsilon_y = \epsilon_z = -\frac{P}{E}(1 - 2\nu)$。

  2. 计算体积应变：
    体积应变 $\epsilon_v = \epsilon_x + \epsilon_y + \epsilon_z$：
    $$\epsilon_v = 3 \times \left[ -\frac{P}{E}(1 - 2\nu) \right] = -\frac{3P(1-2\nu)}{E}$$

这个结果展示了材料在压力作用下体积缩小的程度，直接体现了材料的压缩性。

总结

广义胡克定律是连接“几何变形”与“内力分布”的桥梁。从一维的简单比例，到三维的四阶张量关系，这一演变过程体现了力学描述从简化模型向真实物理世界的逼近。

  • 对于初学者：理解应力与应变的张量本质是掌握固体力学的关键。
  • 对于工程师：掌握各向同性简化形式是进行结构强度计算的基础。
  • 对于研究者：深入研究刚度张量的对称性与材料对称性（如正交各向异性、单斜对称性等）是进行先进复合材料设计的核心。

通过掌握广义胡克定律，我们可以利用有限元分析（FEA）等数值方法，模拟从微观晶体到宏观航空航天结构的复杂力学行为。