2. 连续介质力学基础回顾
在进入聚合物黏弹有限元建模之前，有必要把连续介质力学的四个基本方程重新梳理一遍。它们不是抽象数学，而是有限元软件内部真正求解的对象。理解每个方程在黏弹问题中扮演什么角色，以及小变形与有限变形假设的适用边界，直接决定了建模时该选哪种分析类型、该输入哪些材料参数、该关注哪些输出量。
一、四个基本方程及其在黏弹问题中的角色
连续介质力学描述一个变形体，需要同时回答四个问题：物体受什么力、怎么变形、变形如何协调、材料如何响应。对应四个方程：
- 平衡方程：描述应力场的空间变化与体力、惯性力的关系，即 ∇·σ + f = 0（静力）或 ∇·σ + f = ρa（动力）。它回答“力如何传递”。在黏弹问题中，平衡方程本身不随时间改变，但应力σ(t)是时间的函数，因此每个时刻都要满足平衡。
- 几何方程：描述应变与位移的关系，小变形下 ε = ½(∇u + ∇uᵀ)。它回答“位移如何转化为应变”。黏弹问题中位移随时间演化，几何方程把位移场映射为应变场，是时间步进计算的起点。
- 本构关系：描述应力与应变（及应变率）的关系。弹性材料是 σ = C:ε，黏弹材料则是 σ(t) = ∫₀ᵗ C(t−τ) : dε(τ)/dτ dτ，即卷积形式。它回答“材料如何响应”，是黏弹建模的核心，也是与弹性问题唯一本质不同的方程。
- 协调方程：保证应变场可由单值连续位移场导出，即应变分量之间满足相容条件。它回答“变形是否物理上可能”。在有限元中，协调性由形函数和单元连接自动保证，通常不显式求解，但它是位移法成立的前提。
这四个方程中，平衡方程和几何方程是普适的，与材料无关；本构关系是材料专属的；协调方程是运动学约束。黏弹有限元建模的难点几乎全部集中在本构关系的离散化和时间积分上。
二、本构关系：黏弹问题的唯一“非弹性”入口
弹性有限元中，本构关系退化为一个常数刚度矩阵 C，应力由当前应变直接算出。黏弹问题则不同：应力不仅取决于当前应变，还取决于整个加载历史。以最常见的广义Maxwell模型为例，本构关系可写为：
σ(t) = C∞ : ε(t) + Σᵢ Cᵢ : ∫₀ᵗ exp(−(t−τ)/τᵢ) dε(τ)/dτ dτ其中 C∞ 是长期弹性模量，Cᵢ 和 τᵢ 是第 i 个Maxwell分支的模量和松弛时间。有限元实现时，这个卷积积分被离散为递推公式，每个时间步只需更新内部变量，而不必存储全部历史。这正是黏弹有限元能够高效计算的关键。
实际建模中，用户通常输入的是松弛模量 G(t) 的Prony级数系数，软件内部自动转换为上述递推形式。例如某聚合物在25°C下的剪切松弛模量可拟合为：
| 分支 i | gᵢ = Gᵢ/G₀ | τᵢ (s) |
|---|---|---|
| 1 | 0.35 | 0.01 |
| 2 | 0.25 | 0.1 |
| 3 | 0.20 | 1.0 |
| 4 | 0.10 | 10.0 |
| ∞ | 0.10 | — |
这些参数直接决定时间步长的选择：最小松弛时间 τ₁ = 0.01 s 意味着时间步长应至少小于 0.01 s 才能捕捉快速松弛过程，否则本构积分会产生显著误差。
三、小变形与有限变形假设的适用边界
小变形假设的核心是：位移梯度远小于1，应变与位移呈线性关系，平衡方程可在未变形构型上建立。对于聚合物黏弹问题，这个假设在以下情况成立：
- 应变水平低于约5%，例如大多数动态力学分析（DMA）测试、小振幅振荡剪切（SAOS）；
- 结构刚度足够大，变形不引起几何构型显著改变；
- 不需要考虑大转动或大位移引起的几何非线性。
当应变超过5%～10%，或结构发生大位移、大转动时，必须切换到有限变形框架。此时需要：
- 使用Green-Lagrange应变或对数应变替代小应变；
- 使用第二类Piola-Kirchhoff应力或Kirchhoff应力；
- 本构关系必须满足客观性，即应力率需采用客观应力率（如Jaumann率或Green-Naghdi率）；
- 平衡方程建立在变形后构型上，需引入几何刚度矩阵。
常见误区：很多工程师在聚合物大变形问题中仍使用小变形黏弹模型，结果在应变超过10%后应力预测严重偏高。原因是小变形假设下几何非线性被忽略，而聚合物在大变形下往往伴随应变软化和几何软化，两者叠加会使实际应力远低于小变形预测值。
四、建模时的实际判断流程
面对一个具体的聚合物黏弹有限元问题，可按以下顺序判断：
- 估算最大应变水平。若 ε_max < 5%，优先选小变形黏弹分析；若 ε_max > 10%，必须选有限变形黏弹分析；5%～10%之间需做收敛性验证。
- 确认材料参数是否覆盖实际时间尺度。Prony级数的最小τ应小于最小时间步长，最大τ应大于总分析时间。
- 检查载荷类型。蠕变、应力松弛、动态振荡对应不同的本构积分策略，但都基于同一套Prony参数。
- 验证时间步长。建议初始步长取最小τ的1/5～1/10，并开启自动时间步进。
实用建议：在Abaqus中，小变形黏弹分析使用 *VISCOELASTIC 配合 *STATIC；有限变形黏弹分析则需在 *STEP 中设置 NLGEOM=ON，并确认材料定义中的 *VISCOELASTIC 与 *HYPERELASTIC 或 *ELASTIC 组合使用。切换NLGEOM后，必须重新检查时间步长和收敛容差，因为几何非线性会显著增加迭代次数。
总结来说，平衡方程和几何方程提供了黏弹有限元的“骨架”，本构关系提供了“材料灵魂”，协调方程保证了“运动学合法性”。小变形与有限变形的选择不是精度高低的问题，而是假设是否成立的问题。在聚合物黏弹建模中，先判断应变水平，再选择分析框架，最后配置材料参数和时间步长，这条逻辑链比任何软件操作都更重要。