4. 微分型本构模型

微分型本构模型用“弹簧”和“黏壶”两种基本元件的串并联组合，把应力、应变及其对时间的导数写进同一个方程。弹簧服从胡克定律 σ = Eε，代表瞬时弹性；黏壶服从牛顿黏性 σ = η·dε/dt，代表与速率相关的耗散。把二者组合，就能用常系数线性微分方程描述材料的松弛、蠕变和滞后。本节先建立三种基本模型的方程，再推导各自的松弛模量与蠕变柔量，最后讨论它们对真实材料行为的描述能力与局限。

一、Maxwell模型：串联组合

Maxwell模型把弹簧与黏壶串联。串联元件应力相同、应变相加，因此总应变率等于弹性应变率与黏性应变率之和：

dε/dt = (1/E)·dσ/dt + σ/η

这就是Maxwell模型的本构方程。它最典型的用途是描述应力松弛：在 t = 0 施加并保持恒定应变 ε₀，则 dε/dt = 0，方程退化为 (1/E)·dσ/dt + σ/η = 0，解得

σ(t) = σ₀·exp(−t/τ)，  τ = η/E

其中 τ 称为松弛时间，是应力衰减到初始值 1/e 所需的时间。松弛模量 E(t) = σ(t)/ε₀ = E·exp(−t/τ)。若在恒定应力 σ₀ 下考察蠕变，则 dσ/dt = 0，方程变为 dε/dt = σ₀/η，积分得 ε(t) = σ₀/E + σ₀t/η，即瞬时弹性应变叠加一个随时间线性增长的黏性流动。这说明Maxwell模型能描述瞬时弹性和稳态流动，但蠕变应变会无限增大，无法反映真实材料中应变趋于稳定的行为；同时它只有一个松弛时间，而真实材料的松弛往往跨越多个数量级。

二、Kelvin-Voigt模型：并联组合

Kelvin-Voigt模型把弹簧与黏壶并联。并联元件应变相同、应力相加，因此总应力等于弹性应力与黏性应力之和：

σ = E·ε + η·dε/dt

它最典型的用途是描述蠕变：在 t = 0 施加并保持恒定应力 σ₀，方程变为一阶线性常微分方程，结合初始条件 ε(0) = 0，解得

ε(t) = (σ₀/E)·[1 − exp(−t/τ)]，  τ = η/E

蠕变柔量 J(t) = ε(t)/σ₀ = (1/E)·[1 − exp(−t/τ)]。应变从零开始逐渐趋近于 σ₀/E，表现出延迟弹性。若考察应力松弛，令 ε = ε₀ 恒定，则 dε/dt = 0，方程给出 σ = Eε₀，即应力瞬间达到并保持不变，无法描述松弛。因此Kelvin-Voigt模型能描述延迟弹性蠕变，却不能描述瞬时弹性响应和应力松弛，这是它与Maxwell模型互补的地方。

三、标准线性固体：三元件组合

把Maxwell模型与一个弹簧并联，或把Kelvin-Voigt模型与一个弹簧串联，都得到三元件模型，统称标准线性固体（Standard Linear Solid, SLS）。以“弹簧E₁与Kelvin-Voigt单元（E₂、η）串联”为例，总应变等于弹簧应变与Kelvin-Voigt单元应变之和，而两者应力相同。消去中间变量后可得本构方程：

σ + (η/(E₁+E₂))·dσ/dt = (E₁E₂/(E₁+E₂))·ε + (E₁η/(E₁+E₂))·dε/dt

写成标准形式 σ + p₁·dσ/dt = q₀·ε + q₁·dε/dt，其中 p₁ = η/(E₁+E₂)，q₀ = E₁E₂/(E₁+E₂)，q₁ = E₁η/(E₁+E₂)。在恒定应变 ε₀ 下求解，得到松弛模量：

E(t) = E_∞ + (E₀ − E_∞)·exp(−t/τ_σ)

其中 E₀ 是瞬时模量，E_∞ 是长期平衡模量，τ_σ 是松弛时间。在恒定应力 σ₀ 下求解，得到蠕变柔量：

J(t) = J₀ + (J_∞ − J₀)·[1 − exp(−t/τ_ε)]

标准线性固体的关键改进是：它同时具有瞬时弹性响应、延迟弹性和有限的长期平衡态。应力松弛时应力不会衰减到零，而是趋于 E_∞ε₀；蠕变时应变不会无限增长，而是趋于 σ₀/E_∞。这比Maxwell和Kelvin-Voigt模型更接近真实黏弹性材料在中等时间尺度下的行为。

四、三种模型的描述能力与局限

把三种模型放在一起比较，可以清楚看到它们各自覆盖的时间尺度。Maxwell模型擅长应力松弛和稳态流动，但蠕变无界；Kelvin-Voigt模型擅长延迟弹性蠕变，但无瞬时弹性、不能松弛；标准线性固体同时具备瞬时弹性、延迟弹性和有限平衡态，是三者中唯一能同时给出有界松弛和有界蠕变的模型。不过，标准线性固体仍然只有一个松弛时间，而真实聚合物、生物组织或沥青的松弛谱通常跨越多个数量级，单一τ只能拟合某一时间窗口。

工程上处理这一局限的常用办法是把多个Maxwell单元并联，得到广义Maxwell模型（Prony级数）：

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

每个分支贡献一个松弛时间和一个模量权重，叠加后可以逼近宽频或长时间范围的松弛曲线。同理，多个Kelvin-Voigt单元串联可得到广义Kelvin-Voigt模型，用于拟合多阶段蠕变。参数 Eᵢ、τᵢ 通常由动态力学分析（DMA）或应力松弛实验数据拟合得到：先在对数时间轴上观察 E(t) 或 J(t) 的转折点，初步确定τᵢ的量级，再用最小二乘拟合权重。判断拟合是否可接受，可以看残差是否随机分布、拟合曲线是否覆盖实验时间窗口的两个数量级以上，以及外推到平衡态时 E_∞ 是否为正且物理合理。

实际选择模型时，可以先问三个问题：是否需要描述瞬时弹性？是否需要描述应力松弛？长期变形是否有界？如果三者都需要，标准线性固体是最低要求；如果松弛谱很宽，则用广义Maxwell或广义Kelvin-Voigt。模型不是越复杂越好，而是要在可辨识参数数量与实验数据覆盖范围之间取得平衡。