4. 常用黏弹模型及其工程选择

在聚合物黏弹有限元建模中，模型选择不是“越复杂越好”，而是要在材料真实响应、可获得的实验数据、参数可识别性以及求解代价之间取得平衡。本节从工程使用角度梳理五类常用描述方式：Maxwell、Kelvin-Voigt、标准线性固体（SLS）、广义Maxwell模型和Prony级数，重点说明它们的参数含义、能拟合什么、不能拟合什么，以及在有限元中如何落地。

一、从两个基本元件说起

弹簧服从胡克定律，应力与应变成正比，卸载后变形立即恢复，代表瞬时弹性；黏壶服从牛顿定律，应力与应变率成正比，代表不可逆的黏性流动。Maxwell模型是弹簧与黏壶串联：总应变等于两者应变之和，应力相同。其松弛模量为

E(t) = E · exp(-t/τ),  τ = η / E

其中 E 是瞬时弹性模量，η 是黏度，τ 是松弛时间。Maxwell能描述应力松弛和瞬时弹性，但恒定应力下应变会随时间线性增长（稳态黏性流动），因此不适合描述长期有界蠕变。Kelvin-Voigt模型是弹簧与黏壶并联：两者应变相同，应力相加，蠕变柔量为

J(t) = (1/E) · [1 - exp(-t/τ)],  τ = η / E

它能描述延迟弹性蠕变，但无法表现瞬时弹性响应，也不能描述应力松弛到零的过程。两者各偏一端，这正是工程中选择模型时的基本矛盾。

二、标准线性固体：兼顾瞬时弹性与有界蠕变

标准线性固体由一个弹簧与一个Maxwell元件并联（或等效地由一个弹簧与一个Kelvin元件串联）构成，包含三个参数：瞬时模量 E₀、长期模量 E∞ 和松弛时间 τ。其松弛模量可写为

E(t) = E∞ + (E₀ - E∞) · exp(-t/τ)

当 t=0 时 E(0)=E₀，给出瞬时弹性；当 t→∞ 时 E(t)→E∞，应力松弛到有限值而非零，蠕变也趋于有界。SLS是能同时描述瞬时弹性、应力松弛和有界蠕变的最小模型，参数少、物理意义清晰，适合单一松弛机制的粗略拟合。但真实聚合物往往在多个时间尺度上松弛，单松弛时间无法同时拟合宽频或长时程数据。

三、广义Maxwell模型与Prony级数

广义Maxwell模型把多个Maxwell元件并联，再与一个长期弹簧 E∞ 并联，其松弛模量为

E(t) = E∞ + Σᵢ Eᵢ · exp(-t/τᵢ)

每个支路贡献一个松弛时间 τᵢ 和模量权重 Eᵢ，整体构成离散松弛谱。理论上，只要支路足够多，它可以逼近任意线性黏弹响应。Prony级数正是这一形式的归一化写法，有限元软件通常要求输入无量纲模量比 gᵢ = Eᵢ/E₀ 和松弛时间 τᵢ，其中 E₀ = E∞ + ΣEᵢ 为瞬时模量。例如某聚碳酸酯在20 ℃下的剪切松弛数据，用5阶Prony拟合可写成：

g₁=0.32, τ₁=1e-3 s
g₂=0.25, τ₂=1e-2 s
g₃=0.18, τ₃=1e-1 s
g₄=0.12, τ₄=1e0  s
g₅=0.08, τ₅=1e1  s
g∞=0.05（长期模量比）

注意 Σgᵢ + g∞ = 1，这是Prony级数的归一化约束。τᵢ 的分布应覆盖实验时间窗口，通常按对数等间隔布置，如每十倍频程2～3个点。若实验只做到100 s，却外推到10⁶ s，则长期响应完全由 g∞ 和最大 τ 决定，外推可靠性很低。

四、拟合能力与数值代价的权衡

从Maxwell到广义Maxwell，拟合能力递增，但数值代价也同步上升。在有限元中，每个Maxwell支路引入一个内部状态变量（黏弹性应变或应力分量），显式积分时需存储历史变量，隐式积分时需在单元层面求解局部本构方程。支路数 N 增加，单元计算量和内存占用大致线性增长；若 N 过大，还可能因时间步与最小 τ 不匹配导致数值振荡或收敛困难。

工程选择可参考以下判断：若只关心短时松弛或单一机制，Maxwell或SLS足够；若需同时描述瞬时弹性、有界蠕变和宽频松弛，选广义Maxwell/Prony，阶数通常取3～8阶；若实验数据只覆盖1～2个数量级时间，强行用10阶以上拟合只会引入不可识别的参数，反而降低预测稳定性。此外，Prony级数只适用于线性黏弹范围，大应变、循环加载下的非线性效应需改用并行流变模型或非线性黏弹本构。

五、工程选择速查

实际建模时，先明确分析关心的物理量是松弛模量、蠕变柔量还是动态模量，再根据实验数据的时间/频率窗口选择模型阶数。模型不是越复杂越准确，而是要在可识别、可验证、可承受的数值代价内，复现对工程结论真正重要的那部分黏弹响应。