超弹性材料的本构理论

超弹性（Hyperelasticity）是指某些材料在经历大变形后，在卸载时能够完全恢复到初始状态的特性。这类材料最典型的代表是橡胶、弹性体以及某些生物软组织。与传统的线性弹性材料（遵循胡克定律）不同，超弹性材料的应力-应变关系呈现出强烈的非线性，且变形量通常可达数百个百分点。

在固体力学中，超弹性材料的本构关系并非直接通过实验拟合的应力-应变曲线给出，而是基于一个标量函数——应变能密度函数（Strain Energy Density Function, $W$）。通过对该函数求导，可以推导出材料在任意变形状态下的应力状态。

大变形运动学基础

由于超弹性材料涉及大变形，传统的微小应变理论不再适用，必须采用有限变形理论（Finite Deformation Theory）。

1. 变形梯度张量

描述物体从参考构型（$\mathbf{X}$）到当前构型（$\mathbf{x}$）映射的核心量是变形梯度张量 $\mathbf{F}$：
$$\mathbf{F} = \frac{\partial \mathbf{x}}{\partial \mathbf{X}}$$
其中 $\det(\mathbf{F}) = J$ 表示体积变化率。对于绝大多数橡胶类材料，通常假设其为不可压缩材料，即 $J = 1$。

2. 右柯西-格林变形张量

为了消除刚体旋转的影响，引入右柯西-格林变形张量 $\mathbf{C}$：
$$\mathbf{C} = \mathbf{F}^T \mathbf{F}$$
$\mathbf{C}$ 是一个对称正定张量，其不变量 $\text{I}_1, \text{I}_2, \text{I}_3$ 常用于构建本构模型：

• $\text{I}_1 = \text{tr}(\mathbf{C})$
• $\text{I}_2 = \frac{1}{2} [(\text{tr}\mathbf{C})^2 - \text{tr}(\mathbf{C}^2)]$
• $\text{I}_3 = \det(\mathbf{C}) = J^2$

超弹性本构理论核心

超弹性材料的核心假设是：材料的状态由一个势函数 $W$ 决定，该函数仅依赖于变形梯度 $\mathbf{F}$（或其相关的变形张量）。

1. 应变能密度函数 $W$

应变能密度 $W$ 表示单位参考体积内存储的弹性势能。对于各向同性材料，$W$ 可以表示为 $\mathbf{C}$ 的不变量的函数：
$$W = W(\text{I}_1, \text{I}_2, \text{I}_3)$$

2. 应力与能量的关系

根据热力学定律，第二类皮奥拉-基尔霍夫应力张量 $\mathbf{S}$（Second Piola-Kirchhoff stress）可以通过对 $W$ 关于 $\mathbf{C}$ 求偏导得到：
$$\mathbf{S} = 2 \frac{\partial W}{\partial \mathbf{C}}$$
而实际工程中更常用的柯西应力 $\boldsymbol{\sigma}$（Cauchy stress，即真实应力）则通过推拉关系转换：
$$\boldsymbol{\sigma} = \frac{1}{J} \mathbf{F} \mathbf{S} \mathbf{F}^T$$

对于不可压缩材料，由于 $J=1$ 且存在一个未知的静水压力 $p$，其应力表达式变为：
$$\boldsymbol{\sigma} = -p\mathbf{I} + 2 \mathbf{F} \frac{\partial W}{\partial \mathbf{C}} \mathbf{F}^T$$

常见的超弹性模型

根据对 $W$ 函数的不同数学假设，产生了多种经典的本构模型，适用于不同的应变范围。

1. Neo-Hookean 模型

这是最简单的超弹性模型，基于高分子链的统计力学理论。它假设 $W$ 仅与 $\text{I}_1$ 相关：
$$W = \frac{\mu}{2}(\text{I}_1 - 3)$$
其中 $\mu$ 为剪切模量。该模型在小变形（$<30%$）时精度较高，但在大变形时无法捕捉到橡胶的“硬化”效应。

2. Mooney-Rivlin 模型

为了提高拟合精度，Mooney-Rivlin 模型引入了 $\text{I}2$ 项：
$$W = C{10}(\text{I}1 - 3) + C{01}(\text{I}2 - 3)$$
其中 $C{10}$ 和 $C_{01}$ 是材料常数。该模型能更好地描述中等变形下的剪切行为。

3. Ogden 模型

Ogden 模型不使用不变量，而是直接基于主伸长率 $\lambda_1, \lambda_2, \lambda_3$ 构建：
$$W = \sum_{i=1}^N \frac{\mu_i}{\alpha_i} (\lambda_1^{\alpha_i} + \lambda_2^{\alpha_i} + \lambda_3^{\alpha_i} - 3)$$
通过调整参数 $\mu_i$ 和 $\alpha_i$ 以及项数 $N$，Ogden 模型可以极其精确地拟合极高应变（$>700%$）下的实验数据。

实例分析：单轴拉伸实验

为了理解上述理论，考虑一个不可压缩的 Neo-Hookean 材料块进行单轴拉伸。

• 变形状态：假设拉伸方向为 $X_1$，伸长率为 $\lambda$。由于不可压缩（$J=1$），横向伸长率为 $\lambda_2 = \lambda_3 = 1/\sqrt{\lambda}$。
• 不变量计算：$\text{I}_1 = \lambda^2 + \frac{2}{\lambda}$。
• 应力推导：
  根据 $\boldsymbol{\sigma} = -p\mathbf{I} + \mu \mathbf{B}$（其中 $\mathbf{B} = \mathbf{F}\mathbf{F}^T$），且边界条件 $\sigma_{22} = \sigma_{33} = 0$，可解得压力 $p = \mu/\lambda$。
• 最终应力-应变关系：
  $$\sigma_{11} = \mu (\lambda^2 - \frac{1}{\lambda})$$
  这个公式揭示了超弹性材料典型的非线性特征：应力随伸长率 $\lambda$ 的增加而呈非线性增长。

总结

超弹性本构理论为描述大变形弹性材料提供了严谨的数学框架。其核心逻辑链条为：
变形梯度 $\mathbf{F} \rightarrow$ 变形不变量 $\text{I}_n \rightarrow$ 应变能密度 $W \rightarrow$ 应力 $\boldsymbol{\sigma}$。

在实际工程应用（如轮胎设计、密封圈分析、生物组织模拟）中，选择合适的模型至关重要：

• 低应变/初步分析 $\rightarrow$ Neo-Hookean
• 中等应变/剪切主导 $\rightarrow$ Mooney-Rivlin
• 高应变/精确拟合 $\rightarrow$ Ogden