2. 填料表面能测定与润湿性判据建立

为什么需要“表面能判据”而不是凭肉眼判断润湿

在例1的基准体系中，硅微粉在双酚A环氧树脂中的润湿状态往往靠“看是否挂壁、是否有干粉团”来判断，这种判断只能区分“明显不润湿”和“看起来还行”，无法回答两个关键问题：该填料与树脂体系在热力学上是否具备自发润湿条件，以及不同分散剂能把润湿性改善到什么程度。要建立量化判据，需要把填料表面能拆成色散分量 γ^d 与极性分量 γ^p，再与环氧树脂体系的表面能参数比较，计算润湿包络（wetting envelope）。

本节采用两条互补的测定路线：接触角法（用已知表面能参数的探针液体在压片填料表面测接触角，通过OWRK模型回归出填料 γ^d、γ^p）与毛细渗透法（Washburn法，用粉末柱毛细上升速率反推接触角，适合无法压成光滑片状的填料）。两者得到的接触角代入同一套判据，输出可比较的润湿性指标。

测定对象与探针液体选择

本节以两类典型填料为对象：硅微粉（以SiO₂为主，表面极性中等、易吸附水分）与重质碳酸钙（CaCO₃，表面能偏低、极性分量小，常经硬脂酸表面处理）。探针液体需覆盖“高色散+低极性”与“高极性”两端，才能让OWRK回归有足够的方程条件数。推荐组合如下（参数为文献常用示例值，实际测定前应以所用批次试剂的技术数据为准）：

接触角测定与OWRK回归

硅微粉压片：取干燥填料约0.5 g，在压片机上以约10 MPa保压60 s，得到直径约13 mm的圆片，表面用无水乙醇轻拭后于60 ℃干燥30 min。用座滴法每次滴加2 μL探针液体，滴后5 s内读数，每片测5个点取平均，剔除明显铺展不均的点。碳酸钙压片压力可略低（约8 MPa），避免片体开裂。

OWRK方程把固-液界面能拆成色散与极性两项的几何平均：

γ_L (1 + cosθ) = 2·√(γ_S^d · γ_L^d) + 2·√(γ_S^p · γ_L^p)

令 y = γ_L (1 + cosθ) / (2·√γ_L^d)
令 x = √(γ_L^p / γ_L^d)
则 y = √γ_S^d + √γ_S^p · x      → 线性回归
截距 = √γ_S^d，斜率 = √γ_S^p

下面是一段可直接运行的Python脚本，输入各探针液体的接触角，输出填料的 γ^d、γ^p 与总表面能，并计算与给定环氧树脂体系的润湿包络判据。示例接触角为教学假设值，用于演示计算路径，不代表实测结果。

import numpy as np

# 探针液体参数 (mN/m)：名称, γd, γp
liquids = [
    ("正己烷",   18.4,  0.0),
    ("二碘甲烷", 48.5,  2.3),
    ("乙二醇",   29.0, 19.0),
    ("水",       21.8, 51.0),
]

# 教学假设接触角（度）：硅微粉压片 / 碳酸钙压片
theta_silica = {"正己烷": 12.0, "二碘甲烷": 38.0, "乙二醇": 52.0, "水": 61.0}
theta_caco3  = {"正己烷": 20.0, "二碘甲烷": 47.0, "乙二醇": 66.0, "水": 78.0}

def owrk(theta_deg, liquids):
    xs, ys = [], []
    for name, gd, gp in liquids:
        th = np.radians(theta_deg[name])
        gl = gd + gp
        y = gl * (1 + np.cos(th)) / (2 * np.sqrt(gd))
        x = np.sqrt(gp / gd) if gd > 0 else 0.0
        xs.append(x); ys.append(y)
    xs, ys = np.array(xs), np.array(ys)
    # 最小二乘线性回归 y = a + b x
    b, a = np.polyfit(xs, ys, 1)
    gs_d = a**2
    gs_p = b**2
    return gs_d, gs_p, gs_d + gs_p

for label, th in [("硅微粉", theta_silica), ("碳酸钙", theta_caco3)]:
    gd, gp, gt = owrk(th, liquids)
    print(f"{label}: γd={gd:.1f}, γp={gp:.1f}, γ总={gt:.1f} mN/m")

# 环氧树脂体系表面能（示例：双酚A环氧 + 胺固化剂，未固化混合液）
ep_d, ep_p = 38.0, 8.0
ep_t = ep_d + ep_p

def wetting_envelope(gs_d, gs_p, gl_d, gl_p):
    # 固-液界面能（几何平均）
    gamma_sl = (gs_d + gs_p) + (gl_d + gl_p) \
               - 2*np.sqrt(gs_d*gl_d) - 2*np.sqrt(gs_p*gl_p)
    # 铺展系数 S = γ_S - γ_L - γ_SL
    S = (gs_d + gs_p) - (gl_d + gl_p) - gamma_sl
    # 粘附功 Wa = γ_S + γ_L - γ_SL
    Wa = (gs_d + gs_p) + (gl_d + gl_p) - gamma_sl
    return gamma_sl, S, Wa

for label, th in [("硅微粉", theta_silica), ("碳酸钙", theta_caco3)]:
    gd, gp, gt = owrk(th, liquids)
    gsl, S, Wa = wetting_envelope(gd, gp, ep_d, ep_p)
    print(f"{label} vs 环氧: γSL={gsl:.1f}, S={S:.1f}, Wa={Wa:.1f} mN/m")

判据的物理含义：铺展系数 S > 0 表示液体可在固体表面自发铺展（完全润湿）；S < 0 但粘附功 Wa 较大 表示虽不能完全铺展，但界面结合较强，属于“部分润湿、可分散”区间；Wa 偏小且 γSL 偏高 则提示界面能失配，需要分散剂介入降低界面张力。实际工程中更常用的是把 S 与 Wa 一起看，而不是单看接触角大小。

从判据到分散剂筛选的量化接口

测得填料表面能后，润湿包络给出的是“无分散剂时的热力学基线”。分散剂的作用可理解为在固-液界面引入吸附层，降低 γ_SL、提高 Wa。因此筛选时不必重新测每种分散剂下的填料表面能，而是固定填料参数，比较加分散剂后体系对填料的接触角下降幅度或毛细渗透速率提升倍数，再回代同一套 S、Wa 公式，即可得到可比较的改善量。这样例3、例4中的分散剂类型与用量对比就有了统一的量化标尺，而不是各说各话的“看起来更均匀”。

结果记录与判据卡

每次测定至少记录以下字段，保证不同批次、不同操作者之间可追溯、可复算：

完成本节后，你应得到每个填料批次的 γ^d、γ^p 与对应的 S、Wa 值，并明确其落在“自发铺展 / 部分润湿 / 界面失配”哪一区。这个分区结果就是后续分散剂类型与用量实验的起点判据：界面失配区优先选能显著降低 γ_SL 的分散剂，部分润湿区则重点比较润湿速度的提升幅度。