4. 广义胡克定律与各向异性弹性

在前三章中，我们已经建立了应力张量 σ_ij 与应变张量 ε_kl 的描述方式。要把两者联系起来，就需要引入材料的本构关系。对于绝大多数工程材料在小变形、等温、线弹性范围内，这一关系就是广义胡克定律：应力张量的每一个分量都是应变张量分量的线性函数，反之亦然。写成张量形式为 σ_ij = C_ijkl ε_kl，其中 C_ijkl 是四阶弹性刚度张量。由于应力、应变各自对称，C_ijkl 对前两个指标和后两个指标分别对称，独立分量从 81 个降到 36 个，于是可以用 Voigt 记法把双指标对 (ij)、(kl) 缩并为 1～6 的单指标，得到工程上最常用的 6×6 刚度矩阵形式：

σ₁ = C₁₁ε₁ + C₁₂ε₂ + C₁₃ε₃ + C₁₄γ₄ + C₁₅γ₅ + C₁₆γ₆
σ₂ = C₂₁ε₁ + C₂₂ε₂ + C₂₃ε₃ + C₂₄γ₄ + C₂₅γ₅ + C₂₆γ₆
...
τ₆ = C₆₁ε₁ + C₆₂ε₂ + C₆₃ε₃ + C₆₄γ₄ + C₆₅γ₅ + C₆₆γ₆

这里 ε₁、ε₂、ε₃ 为正应变，γ₄、γ₅、γ₆ 为工程剪应变（γ₄ = 2ε₂₃，γ₅ = 2ε₁₃，γ₆ = 2ε₁₂）。注意 Voigt 记法中剪应变采用工程剪应变而非张量剪应变，这是后续与柔度矩阵、工程常数换算时最容易出错的地方。

全各向异性：36 个常数的完整矩阵

最一般的情况是全各向异性（三斜晶系），材料没有任何对称面，刚度矩阵 36 个分量全部独立。此时正应力会诱发剪应变，剪应力也会诱发正应变，即存在拉剪耦合。C₁₄、C₁₅、C₁₆ 等非对角项就是这种耦合的度量。工程中真正完全三斜的材料不多，但单向复合材料偏轴铺层、某些压电陶瓷、木材斜纹理试样在特定坐标系下都会表现出接近全各向异性的行为，因此理解 36 常数矩阵是理解后续所有简化模型的基础。

刚度矩阵 C 与柔度矩阵 S 互为逆矩阵：S = C⁻¹，C = S⁻¹。刚度矩阵描述“给定应变求应力”，柔度矩阵描述“给定应力求应变”。在实验测量中，我们更容易控制应力（如单轴拉伸、纯剪切），因此工程常数通常先通过柔度矩阵定义，再求逆得到刚度矩阵。

工程常数与柔度矩阵的关系

工程常数包括：三个方向的杨氏模量 E₁、E₂、E₃；三个泊松比 ν₁₂、ν₁₃、ν₂₃（ν_ij 表示 i 方向加载时 j 方向的收缩比）；三个剪切模量 G₂₃、G₁₃、G₁₂。对于正交各向异性材料，柔度矩阵可以完全用这 9 个工程常数写出：

S₁₁ = 1/E₁      S₂₂ = 1/E₂      S₃₃ = 1/E₃
S₁₂ = -ν₁₂/E₁ = -ν₂₁/E₂
S₁₃ = -ν₁₃/E₁ = -ν₃₁/E₃
S₂₃ = -ν₂₃/E₂ = -ν₃₂/E₃
S₄₄ = 1/G₂₃     S₅₅ = 1/G₁₃     S₆₆ = 1/G₁₂

其中泊松比满足互等关系 ν_ij/E_i = ν_ji/E_j，这是柔度矩阵对称性的直接结果，也是检查实验数据是否自洽的重要依据。若实测得到的 ν₁₂/E₁ 与 ν₂₁/E₂ 明显不等，说明试样制备、应变片粘贴方向或数据约化环节存在问题。

正交各向异性与横观各向同性

当材料具有三个相互正交的弹性对称面时，称为正交各向异性。此时拉剪耦合项全部为零，刚度矩阵和柔度矩阵都退化为对角块形式：正应力只与正应变耦合，剪应力只与对应剪应变耦合。独立常数从 36 个降到 9 个：E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₂₃、G₁₃、G₁₂。典型的正交各向异性材料包括单向纤维增强复合材料（沿纤维方向为 1 轴）、正交胶合木、冷轧钢板等。

若材料在垂直于某一轴的平面内各向同性，即存在一个各向同性面，则称为横观各向同性。设 2-3 平面为各向同性面，1 轴为对称轴，则独立常数进一步降到 5 个：E₁（轴向模量）、E₂ = E₃（面内模量）、ν₁₂ = ν₁₃（轴向泊松比）、ν₂₃（面内泊松比）、G₁₂ = G₁₃（轴向剪切模量），而面内剪切模量 G₂₃ = E₂/[2(1+ν₂₃)] 由其余常数导出，不再独立。单向复合材料、拉拔金属丝、沉积岩层在宏观上常按横观各向同性处理。

从实验数据到本构矩阵：一个可执行的换算流程

假设我们有一块单向碳纤维/环氧复合材料层板，已通过单轴拉伸和纯剪切试验测得以下工程常数（示例数据，用于说明换算过程）：E₁ = 135 GPa，E₂ = E₃ = 9.5 GPa，ν₁₂ = ν₁₃ = 0.30，ν₂₃ = 0.45，G₁₂ = G₁₃ = 5.0 GPa。由于是横观各向同性，G₂₃ 由面内关系导出：G₂₃ = E₂/[2(1+ν₂₃)] = 9.5/[2×1.45] ≈ 3.28 GPa。

接下来按以下步骤构建柔度矩阵并求逆得到刚度矩阵：

1. 填写柔度矩阵对角项：S₁₁ = 1/135 = 7.407×10⁻³ GPa⁻¹，S₂₂ = S₃₃ = 1/9.5 = 0.10526 GPa⁻¹，S₄₄ = 1/3.28 = 0.30488 GPa⁻¹，S₅₅ = S₆₆ = 1/5.0 = 0.2 GPa⁻¹。
2. 填写耦合项：S₁₂ = S₁₃ = -ν₁₂/E₁ = -0.30/135 = -2.222×10⁻³ GPa⁻¹，S₂₃ = -ν₂₃/E₂ = -0.45/9.5 = -0.04737 GPa⁻¹。
3. 检查对称性：ν₂₁ = ν₁₂·E₂/E₁ = 0.30×9.5/135 = 0.0211，验证 S₁₂ = -ν₂₁/E₂ = -0.0211/9.5 = -2.222×10⁻³ GPa⁻¹，与第 2 步一致。
4. 对 6×6 柔度矩阵求逆，得到刚度矩阵 C。由于正交各向异性矩阵为块对角形式，可分别对 3×3 正应力块和 3×3 剪应力块求逆，计算量远小于直接对 6×6 求逆。

求逆后可得 C₁₁ ≈ 136.5 GPa，C₁₂ = C₁₃ ≈ 3.05 GPa，C₂₂ = C₃₃ ≈ 10.2 GPa，C₂₃ ≈ 4.35 GPa，C₄₄ ≈ 3.28 GPa，C₅₅ = C₆₆ = 5.0 GPa。可以看到 C₁₁ 略大于 E₁，这是因为单轴拉伸时横向收缩受到泊松效应约束，刚度矩阵对角项并不等于杨氏模量。这一区别在有限元材料卡填写时非常关键：软件通常要求输入工程常数而非刚度矩阵分量，若误把 C₁₁ 当作 E₁ 填入，会导致轴向刚度偏高约 1%。

独立常数的物理意义与工程判断

理解独立常数数量，不只是记忆数字，而是判断“哪些实验必须做、哪些可以推算”。对于横观各向同性材料，5 个独立常数意味着至少需要 5 组独立试验：轴向拉伸（E₁、ν₁₂）、横向拉伸（E₂、ν₂₃）、面内剪切（G₁₂）。G₂₃ 不需要单独测量，可由 E₂ 和 ν₂₃ 导出。若材料实际是正交各向异性而非横观各向同性，则 G₂₃ 必须单独测量，否则面外剪切响应会严重失真。

在实际工程中，判断材料属于哪一类对称性，通常先观察微观结构或铺层方式：单向带材、拉拔棒材、沉积岩层通常按横观各向同性处理；正交编织布、胶合木、冷轧板则按正交各向异性处理；只有确实存在明显偏轴耦合或三斜对称时才使用全各向异性模型。过度简化会丢失关键耦合，过度复杂化则会引入大量难以测准的常数，反而降低模型可信度。