3. 线性黏弹性本构理论

线性黏弹性（linear viscoelasticity, LVE）是聚合物有限元建模中使用频率最高的材料模型族。它的核心假设只有两条：响应与载荷成正比（线性），以及响应具有时间记忆（黏弹性）。只要材料处于小应变、等温、未损伤区间，绝大多数热塑性塑料、热固性树脂、橡胶基体乃至部分生物组织都可以用 LVE 描述。本节把 Boltzmann 叠加原理、松弛模量、蠕变柔量、复模量之间的数学关系串成一条线，并说明如何从连续形式的本构方程走向有限元可用的离散递推格式。

一、Boltzmann 叠加原理：记忆的数学表达

设材料在 t=0 时受到一个阶跃应变 ε₀，其应力响应为 σ(t)=E(t)·ε₀，其中 E(t) 称为松弛模量（relaxation modulus），它随时间单调下降，反映应力松弛。若应变历史是任意可微的 ε(t)，Boltzmann 原理认为：每一段微小应变增量 dε(τ) 都独立产生一份应力响应，且各份响应可以线性叠加。于是

σ(t) = ∫₋∞ᵗ E(t−τ) · dε(τ)/dτ dτ

这就是遗传积分（hereditary integral）形式的本构方程。它说明当前应力不仅取决于当前应变，还取决于全部历史应变，权重由 E(t−τ) 给出。E(t) 是 LVE 的“记忆核”。

二、松弛模量与蠕变柔量：一对互逆的响应函数

与松弛模量对偶的是蠕变柔量（creep compliance）D(t)：在阶跃应力 σ₀ 下，应变响应为 ε(t)=D(t)·σ₀。E(t) 与 D(t) 并非简单倒数关系，而是在 Laplace 域中满足 s²·Ē(s)·D̄(s)=1。工程上常用 Prony 级数（广义 Maxwell 模型）表示松弛模量：

E(t) = E_∞ + Σᵢ Eᵢ · exp(−t/τᵢ)
D(t) = D₀ + Σⱼ Dⱼ · (1 − exp(−t/τⱼ)) + t/η₀

其中 E_∞ 是平衡模量，τᵢ 是松弛时间，η₀ 是稳态黏度。Prony 级数把连续谱离散为若干支 Maxwell 单元并联，这正是有限元实现的基础。

三、复模量：频域中的同一材料

对稳态简谐激励 ε(t)=ε₀·sin(ωt)，线性响应为 σ(t)=ε₀·[E′(ω)·sin(ωt)+E″(ω)·cos(ωt)]。定义复模量 E*(ω)=E′(ω)+i·E″(ω)，其中 E′ 为储能模量，E″ 为损耗模量，tanδ=E″/E′ 为损耗因子。E*(ω) 与 E(t) 通过 Fourier/Laplace 变换一一对应：

E*(ω) = iω · ∫₀^∞ E(t) · e^(−iωt) dt

因此 DMA 测得的 E′(ω)、E″(ω) 可以反演为时域 Prony 参数，反之亦然。这是“同一材料、两种视角”的典型体现。

四、从本构方程到有限元离散

有限元求解需要的是增量递推而非积分形式。以广义 Maxwell 模型为例，把总应力分解为平衡部分与各支非平衡部分之和：σ(t)=E_∞·ε(t)+Σᵢ hᵢ(t)，其中 hᵢ 是第 i 支 Maxwell 单元的内变量，满足一阶常微分方程 ḣᵢ + hᵢ/τᵢ = Eᵢ·ε̇。在时间步 [tₙ₋₁, tₙ] 内假设应变线性变化，可解析积分得到递推公式：

hᵢ,ₙ = exp(−Δt/τᵢ)·hᵢ,ₙ₋₁ + Eᵢ·(1−exp(−Δt/τᵢ))/(Δt/τᵢ) · Δε
σₙ = E_∞·εₙ + Σᵢ hᵢ,ₙ

该格式只需存储每个积分点的 hᵢ，无需保存全部历史，计算量与 Prony 项数成正比。切线刚度则由 ∂σₙ/∂Δε 给出，用于牛顿迭代。若使用 ABAQUS，可直接以 *VISCOELASTIC, TIME=PRONY 输入 Eᵢ、τᵢ；若自编 UMAT，则按上述递推实现。

五、工程示例：DMA 数据到 Prony 参数

假设某 PC 材料在 25 ℃ 下 DMA 频率扫描得到：E′(1 Hz)=2.4 GPa，E′(10 Hz)=2.7 GPa，E′(100 Hz)=2.9 GPa，tanδ 峰值约 0.03。用三支 Maxwell 拟合，可取 τ₁=0.01 s、τ₂=0.1 s、τ₃=1 s，对应 E₁=0.15 GPa、E₂=0.25 GPa、E₃=0.30 GPa，E_∞=2.0 GPa。将上述参数写入材料卡片后，做单轴松弛仿真：加载 1% 应变并保持 100 s，提取应力-时间曲线，应与 E(t)=2.0+0.15e^(−100t)+0.25e^(−10t)+0.30e^(−t) GPa 吻合。若偏差超过 5%，需检查时间步长是否小于最小 τᵢ 的 1/5，或增加 Prony 项数。

六、常见误区与检查要点

• 把 E(t) 与 D(t) 当倒数：仅在 t→0 或 t→∞ 的极限下近似成立，中间时段必须用卷积关系。
• 忽略平衡模量 E_∞：若 E_∞=0，长期载荷下应力会完全松弛，与热固性材料不符。
• 时间步长过大：Δt 超过最小 τᵢ 时，递推公式的指数项失真，导致应力振荡或漂移。
• 频域参数直接当时域用：E′(ω) 不等于 E(t)，必须经变换或拟合。