粘弹性材料的本构模型简介

粘弹性材料兼具弹性体的瞬时恢复能力和黏性流体的时间依赖特性，在工程结构、地质介质、生物组织等领域都有广泛应用。对这类材料进行力学分析时，必须选取合适的本构模型来描述其应力–应变关系及其随时间的演化。本文系统梳理粘弹性材料的基本概念、常用本构模型及其参数辨识方法，并通过简化实例展示模型的实际使用。

基本概念

• 应力 (\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}) 为长期模量。该模型兼具松弛与蠕变特性，适用于多数工程材料的近似描述。

常用本构模型

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
]

该模型在数学上等价于一条连续的松弛谱，能够简洁地描述“幂律”衰减行为。

参数辨识与实验方法

1. 应力松弛实验

   • 步进应变 (\varepsilon_0) 保持不变，记录应力随时间的衰减 (\sigma(t))。
   • 采用最小二乘法拟合松弛函数 (G(t))，得到模型参数。

2. 蠕变实验

   • 步进应力 (\sigma_0) 保持不变，记录应变随时间的增长 (\varepsilon(t))。
   • 通过逆卷积或直接拟合蠕变函数 (J(t)) 获得参数。

3. 动态力学分析（DMA）

   • 在不同频率 (\omega) 下施加小幅正弦应变，测得储能模量 (E'(\omega)) 与损耗模量 (E''(\omega))。
   • 将实验数据转换为时域松弛谱（如采用 Prony 展开），得到一组 ({E_i,\eta_i})。

4. 数值优化

   • 使用 Levenberg‑Marquardt、遗传算法等全局优化手段，避免局部最小导致的参数偏差。

示例计算

假设对一种聚合物进行应力松弛实验，得到以下离散数据（单位：MPa）：

采用 两支路标准线性固体模型（即广义 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）配合数值优化是获取模型参数的常规路径。
• 在工程实践中，往往在 精度需求 与 计算成本 之间做权衡，选取最合适的本构模型，以实现可靠的结构分析和材料设计。

通过本文的介绍，读者应能够快速判断不同粘弹性本构模型的适用场景，掌握基本的参数辨识流程，并在实际工程中进行模型的实现与验证。