4. 吸湿等温线与扩散模型

吸湿等温线描述的是在恒定温度下,材料平衡吸湿量随环境相对湿度(RH)或水蒸气活度 aw 的变化关系。它是判断高分子材料吸湿机理、评估环境可靠性的基础曲线。同一组数据用不同吸附模型拟合,得到的参数含义完全不同:Henry 模型反映理想溶解,Langmuir 模型反映有限位点的单层吸附,BET 模型描述多层吸附,Flory-Huggins 则从聚合物—溶剂相互作用出发描述非理想混合。理解这些模型的适用边界,比记住公式更重要。

4.1 四类吸附模型的物理图像与判据

Henry 模型最简单:C = S·aw,吸湿量与活度成正比,等温线为过原点的直线。它适用于水在疏水聚合物(如聚乙烯、聚四氟乙烯)中的稀溶液行为,此时水分子彼此孤立、无相互作用。若实测等温线在低活度区呈直线、高活度区上翘,说明出现了团聚或微孔填充,Henry 模型失效。

Langmuir 模型假设吸附位点数量有限且各点等价、单层吸附:C = Cm·b·aw / (1 + b·aw)。曲线在低活度区近似线性,高活度区趋于饱和平台 Cm。亲水性基团(酰胺、羟基、羧基)与水形成氢键的位点吸附,常可用 Langmuir 描述。判断依据是等温线是否出现明显饱和平台。

BET 模型把 Langmuir 扩展到多层吸附:C = Cm·K·aw / [(1 − aw)(1 + (K−1)aw)]。它适用于 0.05 < aw < 0.5 区间,可给出单层饱和量 Cm 和 BET 常数 K。当 aw 接近 1 时 BET 曲线发散,因此不能外推到高湿区。

Flory-Huggins 模型从混合自由能出发:ln aw = ln φw + (1 − φw) + χ(1 − φw)2,其中 φw 为水体积分数,χ 为相互作用参数。它适合描述高吸湿聚合物(如尼龙、环氧)在整个湿度范围内的行为,尤其能反映高活度区的强上翘。若 χ 随浓度变化,说明体系偏离简单 Flory-Huggins 假设,需要引入浓度依赖的 χ(φ)。

模型选择经验:低吸湿、疏水体系优先 Henry;有明显饱和平台用 Langmuir;中低湿区多层吸附用 BET;高吸湿、全湿度范围用 Flory-Huggins。拟合时至少报告 R² 和参数物理合理性,不能只看 R²。

4.2 Fick 扩散与非 Fick 行为

吸湿动力学常用 Fick 第二定律描述:∂C/∂t = D·∂²C/∂x²。对厚度 2l 的平板,在恒定表面浓度、初始浓度均匀的条件下,短时间吸湿量 Mt 与平衡吸湿量 M∞ 之比近似为:

M_t / M_∞ ≈ (4/√π)·√(D·t / l²)   (M_t/M_∞ < 0.6)

以 M_t/M_∞ 对 √t 作图,若为过原点直线,则符合 Fick 行为,
由斜率 k 可得:D = π·(k·l/4)²

实际操作中,把试样置于恒温恒湿箱,定期取出称重,记录 Mt,直到连续三次称重变化小于 0.1% 视为平衡。用上述公式计算 D。若 Mt/M∞–√t 曲线出现明显非线性、S 形或长时间不收敛,则属于非 Fick 行为。

非 Fick 行为的常见来源包括:聚合物链段松弛时间与扩散时间可比(Case II 扩散,吸湿前沿以恒定速度推进)、材料内部存在微孔或界面、吸湿引起塑化使 D 随浓度变化、以及表面结皮或氧化层阻碍。判断方法:用 Mt/M∞ 对 t 作图,若呈直线则为 Case II;若曲线呈两段斜率,则可能存在双阶段扩散。

4.3 扩散系数的测定与温度依赖性

扩散系数测定常用三种方法:重量法(上述称重法,设备简单,适合板材)、动态蒸汽吸附(DVS,可连续记录,适合粉末和薄膜)、以及渗透池法(适合薄膜,测稳态渗透率)。重量法要注意试样厚度均匀、边缘密封,避免侧面吸湿干扰一维假设。

温度对扩散系数的影响通常符合 Arrhenius 关系:D = D0·exp(−Ea/RT)。以 ln D 对 1/T 作图,斜率即 −Ea/R。例如某环氧试样在 30、40、50、60 ℃ 下测得 D 分别为 1.2×10−12、2.0×10−12、3.3×10−12、5.4×10−12 m²/s(示例数据),拟合得 Ea ≈ 42 kJ/mol。若 ln D–1/T 曲线出现折点,说明发生了玻璃化转变或次级松弛,此时不能简单外推到使用温度。

需要特别注意的是,吸湿等温线给出的是平衡态信息,扩散系数给出的是动力学信息,两者必须配套使用。工程评估中,先用等温线确定使用湿度下的 M∞,再用 D 和 Arrhenius 外推估算达到该吸湿量所需时间,才能判断器件在寿命期内是否会发生吸湿失效。

吸湿等温线模型与 Fick 扩散判据 等温线模型(C 对 a_w) 水活度 a_w 吸湿量 C Henry Langmuir 饱和平台 C_m BET Flory-Huggins 吸湿动力学:M_t/M_∞ 对 √t √t M_t / M_∞ Fick 线性段 斜率 k → D = π(k·l/4)² 非 Fick(S 形) Case II(对 t 线性) M_∞ 等温线定平衡吸湿量 M_∞,扩散曲线定 D 与机理;两者配套才能评估寿命期吸湿风险

4.4 工程应用中的参数配套

在环境可靠性评估中,建议按以下顺序组织数据:第一步,在目标温度下测 3~5 个湿度点的平衡吸湿量,绘制等温线并选择合适模型拟合,得到使用湿度下的 M∞;第二步,在至少三个温度下测吸湿动力学曲线,用短时间线性段计算 D,并做 Arrhenius 拟合;第三步,用 D 和 M∞ 估算达到 90% 平衡吸湿所需时间 t0.9 ≈ 0.85·l²/D(平板近似)。若 t0.9 小于器件预期服役时间,则必须考虑吸湿引起的尺寸变化、介电常数漂移或界面分层风险。

最后提醒:所有模型参数都应在报告中标明测试温度、湿度范围、试样厚度和拟合区间。BET 参数不能外推到 aw > 0.5,Arrhenius 外推不应跨越玻璃化转变温度。只有参数边界清晰,吸湿等温线与扩散模型才真正具备工程判断价值。