粘弹性材料的本构模型简介
粘弹性材料兼具弹性体的瞬时恢复能力和黏性流体的时间依赖特性，在工程结构、地质介质、生物组织等领域都有广泛应用。对这类材料进行力学分析时，必须选取合适的本构模型来描述其应力–应变关系及其随时间的演化。本文系统梳理粘弹性材料的基本概念、常用本构模型及其参数辨识方法，并通过简化实例展示模型的实际使用。
基本概念
- 应力 (\sigma)、应变 (\varepsilon)：宏观力学变量，分别描述内部力和形变。
- 松弛函数 (G(t))：在保持恒定应变的条件下，材料随时间衰减的应力响应。
[
\sigma(t)=G(t),\varepsilon_0
] - 蠕变函数 (J(t))：在保持恒定应力的条件下，材料随时间增长的应变响应。
[
\varepsilon(t)=J(t),\sigma_0
] - 线性粘弹性：满足叠加原理，即应力和应变均为时间卷积形式。
[
\sigma(t)=\int_{0}^{t} G(t-\tau),\dot{\varepsilon}(\tau),d\tau
]
线性粘弹性模型
1. Maxwell 模型
由弹簧（弹性模量 (E)）与黏性阻尼器（粘度 (\eta)）串联组成。其松弛函数为
[
G(t)=E,e^{-\frac{E}{\eta}t}
]
特点：能够描述应力的快速松弛，但在长期载荷下应变趋于无界（流动）。
2. Kelvin‑Voigt（Kelvin）模型
弹簧与阻尼器并联。蠕变函数为
[
J(t)=\frac{1}{E}\left(1-e^{-\frac{E}{\eta}t}\right)
]
特点：能够捕捉瞬时弹性响应和随时间的迟滞，但不具备应力松弛能力。
3. 标准线性固体（Standard Linear Solid, SLS）模型
在 Maxwell 支路上并联一个弹簧，常用的表达式为
[
G(t)=E_{\infty}+\left(E_{0}-E_{\infty}\right)e^{-\frac{t}{\tau}},\qquad \tau=\frac{\eta}{E_{1}}
]
其中 (E_{0}) 为瞬时模量，(E_{\infty}) 为长期模量。该模型兼具松弛与蠕变特性，适用于多数工程材料的近似描述。
常用本构模型
| 模型 | 结构 | 关键参数 | 适用范围 |
|---|---|---|---|
| Maxwell | 串联弹簧+阻尼器 | (E,\eta) | 高温塑性、流体化材料 |
| Kelvin‑Voigt | 并联弹簧+阻尼器 | (E,\eta) | 低频振动、粘性固体 |
| 标准线性固体 | Maxwell + 并联弹簧 | (E_{0},E_{\infty},\eta) | 大多数聚合物、土体 |
| 广义 Maxwell（多支路） | 多个 Maxwell 支路并联 | ({E_i,\eta_i}) | 需要精确拟合宽频谱松弛行为 |
| 广义 Kelvin（多支路） | 多个 Kelvin 支路并联 | ({E_i,\eta_i}) | 需要精确描述蠕变曲线 |
| 分数阶模型 | 采用分数阶导数 | (E,\eta,\alpha)（(0<\alpha<1)） | 长时记忆效应显著的材料 |
1. 广义 Maxwell 模型的松弛函数
[
G(t)=\sum_{i=1}^{N} E_i,e^{-\frac{E_i}{\eta_i}t}
]
通过增加支路数 (N)，可在宽频带内逼近实验测得的松弛模量曲线。
2. 分数阶模型（常用的 Scott‑Blair 模型）
[
\sigma(t)=E,\frac{d^{\alpha}\varepsilon(t)}{dt^{\alpha}},\qquad 0<\alpha<1
]
该模型在数学上等价于一条连续的松弛谱，能够简洁地描述“幂律”衰减行为。
参数辨识与实验方法
应力松弛实验
- 步进应变 (\varepsilon_0) 保持不变，记录应力随时间的衰减 (\sigma(t))。
- 采用最小二乘法拟合松弛函数 (G(t))，得到模型参数。
蠕变实验
- 步进应力 (\sigma_0) 保持不变，记录应变随时间的增长 (\varepsilon(t))。
- 通过逆卷积或直接拟合蠕变函数 (J(t)) 获得参数。
动态力学分析（DMA）
- 在不同频率 (\omega) 下施加小幅正弦应变，测得储能模量 (E'(\omega)) 与损耗模量 (E''(\omega))。
- 将实验数据转换为时域松弛谱（如采用 Prony 展开），得到一组 ({E_i,\eta_i})。
数值优化
- 使用 Levenberg‑Marquardt、遗传算法等全局优化手段，避免局部最小导致的参数偏差。
示例计算
假设对一种聚合物进行应力松弛实验，得到以下离散数据（单位：MPa）：
| 时间 (t) (s) | 应力 (\sigma(t)) |
|---|---|
| 0 | 10.0 |
| 1 | 7.2 |
| 5 | 4.5 |
| 20 | 2.1 |
| 100 | 0.8 |
采用 两支路标准线性固体模型（即广义 Maxwell，(N=2)）：
[
G(t)=E_1 e^{-t/\tau_1}+E_2 e^{-t/\tau_2}
]
使用最小二乘拟合得到：
- (E_1 = 6.5\ \text{MPa},\ \tau_1 = 2.3\ \text{s})
- (E_2 = 2.8\ \text{MPa},\ \tau_2 = 35\ \text{s})
验证：在 (t=5) s 时，预测应力
[
\sigma(5)=6.5,e^{-5/2.3}+2.8,e^{-5/35}\approx4.48\ \text{MPa}
]
与实验值 4.5 MPa 误差 < 1%。该模型即可用于后续的有限元时域分析或寿命预测。
小结
- 粘弹性材料的本构模型本质上是 应力–应变时域卷积 的数学表述，核心在于松弛函数或蠕变函数的选取。
- 线性模型（Maxwell、Kelvin‑Voigt、标准线性固体）适用于小应变范围，结构简洁、参数易辨识。
- 广义模型（多支路 Maxwell/Kelvin）和 分数阶模型 能够更精确地捕捉宽频带或幂律衰减行为，但参数辨识相对复杂。
- 实验手段（松弛、蠕变、DMA）配合数值优化是获取模型参数的常规路径。
- 在工程实践中，往往在 精度需求 与 计算成本 之间做权衡，选取最合适的本构模型，以实现可靠的结构分析和材料设计。
通过本文的介绍，读者应能够快速判断不同粘弹性本构模型的适用场景，掌握基本的参数辨识流程，并在实际工程中进行模型的实现与验证。