一、同一批接触角,换个模型差多少
表面自由能不是直接测出来的量,它由接触角数据推算得到。同一组接触角交给不同的计算模型,结果会有差异。为了把模型差异与数据波动分开,下面用同一个已知表面做对比:先设定表面的色散分量与极性分量,用 OWRK 模型反推出各类液体在该表面上的接触角,再把同一组接触角分别交给不同模型求解。
需要说明的是,这组数据由 OWRK 模型生成,所以 OWRK 的还原结果精确属于必然。这样设计的目的不是判断哪个模型更准,而是量化模型之间的差异量级,以及这个差异在不同表面上是否一致。
| 计算模型 | 中等极性聚合物 (真值 44.0 mN/m) | 亲水处理后的玻璃 (真值 68.0 mN/m) | 低表面能涂层 (真值 26.0 mN/m) |
|---|---|---|---|
| OWRK | 43.98(−0.02) | 68.01(+0.01) | 26.00(0.00) |
| Wu | 46.75(+2.75) | 72.49(+4.49) | 29.77(+3.77) |
| van Oss & Good | 42.96(−1.04) | 数据不足 | 25.75(−0.25) |
| Fowkes(仅求色散分量) | 35.98 | 41.98 | 24.00 |
| Zisman(临界表面张力) | 数据不足 | 数据不足 | 23.84 |
括号内是与真值的差。三点结论值得记下来。
模型之间的差异量级为 3 至 4.5 mN/m。这一量级换算成相对偏差,在亲水玻璃上是 6.6%,在低表面能涂层上是 14.5%。同样是换个模型,低表面能材料上的相对偏差更大,因为总表面能的基数小。
Wu 的结果系统性偏高。三种表面上分别高出 2.75、4.49、3.77 mN/m,方向一致。Wu 采用调和平均的写法,与 OWRK 的几何平均不同,对低能表面给出更高估值是这个模型的已知特性,不是计算错误。
极性分量的差异比总量更大。中等极性表面上,真值极性分量为 8.0 mN/m,Wu 给出 12.04,高出 50%;低表面能涂层上,真值 2.0 mN/m,Wu 给出 4.07,翻了一倍。如果分析目的是判断表面极性、匹配涂料配方,那么分量层面的差异比总量差异更需要注意。
还有一层区别:Fowkes 与 Zisman 求的不是同一个物理量。Fowkes 只给色散分量,Zisman 给临界表面张力 γc。低表面能涂层上 Zisman 给出 23.84 mN/m,与总表面能 26.0 mN/m 相差 2.16 mN/m,两者接近但不等价。把 γc 当成表面自由能写进报告,是一个容易被忽略的口径错误。

二、液体搭配:中心值都对,可信度差几倍
选定模型之后,下一个变量是测试液体的选择。用同一组表面(色散 36.0、极性 8.0 mN/m),分别用不同液体组合生成接触角并求解 OWRK,结果如下。
| 液体组合 | 条件数 | 还原的 γs | ±2° 误差下的结果区间 | 区间宽度 |
|---|---|---|---|---|
| 水 + 二碘甲烷 | 1.90 | 43.97 | 42.27 ~ 45.68 | 3.41 |
| 水 + 1-溴萘 | 1.95 | 43.98 | 42.45 ~ 45.51 | 3.06 |
| 水 + 二碘甲烷 + 乙二醇 | 2.32 | 43.98 | 42.59 ~ 45.36 | 2.77 |
| 水 + 二碘甲烷 + 乙二醇 + 甲酰胺 | 2.75 | 43.98 | 42.66 ~ 45.28 | 2.62 |
| 水 + 乙二醇 | 6.52 | 44.02 | 41.56 ~ 47.27 | 5.71 |
| 甘油 + 甲酰胺 | 13.81 | 43.94 | 41.63 ~ 50.29 | 8.66 |
| 乙二醇 + 甘油 | 27.42 | 43.90 | 41.41 ~ 57.69 | 16.29 |
表里有一处需要单独指出:还原的 γs 一列,从 43.90 到 44.02,几乎所有组合都给出相近的结果。单看这一列,会得出液体怎么选都一样的结论。
差异藏在不确定度里。假设接触角存在 ±2° 的测量误差,把误差传递到结果上,水配二碘甲烷给出 42.27 至 45.68 的区间,宽度 3.41 mN/m;换成乙二醇配甘油,区间变成 41.41 至 57.69,宽度 16.29 mN/m,是前者的 4.8 倍。
也就是说,液体搭配不佳时,算出的数字看上去正常,但它对接触角误差的放大倍数明显更高。同样 2° 的测量误差,在良态组合下不影响结论,在病态组合下足以让结果跨过判断阈值。
条件数是量化这一项的指标:在 5 以内说明结果对接触角误差不敏感,5 到 20 之间可接受,超过 20 提示液体搭配一般,建议换用分量结构差异更大的组合。上表里乙二醇配甘油的条件数 27.42,已经落在提示区间。药剂选择上,水配二碘甲烷这类极性差异大的组合条件数较低,原因是二碘甲烷没有极性分量,方程组的两个方向被拉开,稳定性随之提高。

三、接触角数据本身的两类问题
液体搭配之后,还有一类问题出在接触角读数本身。
1. 接近完全铺展的读数不携带信息
高能表面上液滴会完全铺展开,接触角接近 0°。以水 12°、二碘甲烷 30° 这样的读数为例,仍能算出一个结果,但这组数据对表面能的约束很弱:接触角已经接近下限,继续减小不再反映表面能的真实变化。
同一类现象在另一侧也有:高能表面上用低表面张力液体去测,液滴直接铺展,cosθ 触及物理上限 1。这类数据点在计算中会被截断,参与拟合反而会把结果往低处拉。测量时如果发现某种液体在高能表面上铺展,换用表面张力更高的液体更合适。
2. 全是非极性液体时求不出极性分量
OWRK 与 vOCG 需要把表面能分解成两个以上分量,方程组里的未知数多于方程数时无法求解。实际会遇到的两种情况:
液体都不含极性分量,例如只用烷烃类液体。此时极性分量的系数全为零,方程无法解出该分量,计算会被拒绝。用户容易踩到的是手上只有水和乙二醇两种液体,以为数量够了,实际两种都带极性,同样解不出稳定结果。
液体只有一种。一个方程解不出两个未知量,这种情况会直接提示液体数量不足。
反过来,如果液体里没有任何非极性液体,数学上有时仍能求解,但条件数会明显抬高。上表里水配乙二醇的条件数 6.52,远高于水配二碘甲烷的 1.90,就是这个原因。这一项在工具里会单独提示,提示的是缺非极性液体,而不是数量不够。
四、测量误差怎么传到结果上
接触角读数的误差按什么比例传到表面能,下面是同一组基准数据上施加不同扰动的结果。
| 扰动 | γs | γd | γp | 相对基准偏差 |
|---|---|---|---|---|
| 水 +2° | 43.04 | 35.98 | 7.06 | −2.11% |
| 水 −2° | 44.94 | 35.98 | 8.96 | +2.21% |
| 二碘甲烷 +2° | 43.22 | 34.89 | 8.33 | −1.71% |
| 二碘甲烷 −2° | 44.72 | 37.06 | 7.66 | +1.71% |
| 两者各 +2° | 42.27 | 34.89 | 7.38 | −3.87% |
| 两者反向各 2° | 43.82 | 37.06 | 6.75 | −0.35% |
基准接触角为水 71.1°、二碘甲烷 46.9°,基准结果 43.97 mN/m。
两条规律可以直接用。单一液体 2° 的误差带来约 2% 的结果偏差;两种液体的误差同向时叠加,达到 −3.87%,反向时相互抵消,只剩 −0.35%。
误差抵消这一条有实用含义。如果两种液体读数的误差方向一致,比如同一批数据都由同一位操作员读取、读数习惯偏大,抵消机制不起作用,结果偏差反而更大。跨实验室比对时,这类系统性偏差比随机误差更难发现。
从分量上看还有更细的规律。水接触角的变化主要影响极性分量,从 7.99 变为 7.06 或 8.96;二碘甲烷的变化主要影响色散分量,从 35.98 变为 34.89 或 37.06。这与两种液体的分量结构对应:水的极性分量高,二碘甲烷是非极性液体。想改善某一分量的精度,应当从对应类型的液体入手。

五、从接触角到粘附功
接触角数据落到工艺上,常用判据是表面处理效果。下面是一组不同处理强度下的读数与计算结果。
| 处理状态 | γs | γd | γp | 极性占比 | 水的粘附功 |
|---|---|---|---|---|---|
| 未处理(水 95°、二碘甲烷 62°) | 28.93 | 27.42 | 1.51 | 5.2% | 66.46 mN/m |
| 轻度处理(水 72°、二碘甲烷 45°) | 44.27 | 37.01 | 7.26 | 16.4% | 95.30 mN/m |
| 中度处理(水 45°、二碘甲烷 35°) | 61.94 | 42.03 | 19.92 | 32.2% | 124.28 mN/m |
| 充分处理(水 20°、二碘甲烷 28°) | 75.27 | 45.03 | 30.24 | 40.2% | 141.21 mN/m |
处理带来的变化分两部分看。总表面能从 28.93 升到 75.27 mN/m,增幅 160%;极性占比从 5.2% 升到 40.2%。色散分量也从 27.42 升到 45.03,但幅度小得多。这说明处理主要改变的是极性分量,与等离子、电晕处理引入极性基团的作用机理一致。
极性占比这一项在实际判断时比总量更灵敏。表面能从 44.27 到 61.94 的跨度里,总量增加 40%,极性占比从 16.4% 升到 32.2%,接近翻倍。如果处理的目标是改善涂覆或粘接,关注极性占比的变化比关注总量更贴合用途。
粘附功一项可以直接对应工艺判断。水的粘附功按 γL(1 + cosθ) 计算,未处理时为 66.46 mN/m,低于水的表面张力 72.8 mN/m,说明水在表面不铺展,涂料与胶黏剂同样缺少附着基础;处理后升到 141.21 mN/m,接近水表面张力的两倍,铺展与附着的条件明显改善。
判断阈值从粘附功入手比从表面能绝对值入手更直观:粘附功大于液体表面张力时,液体能自发铺展。这条判据与液体种类直接挂钩,不需要为每种工艺另定一个表面能阈值。
六、让结果具备可比性的四个条件
把上面的变量收起来,让表面自由能数据可以跨批次、跨实验室比对,需要固定四件事。
计算模型。OWRK、Wu、vOCG 之间的差异可达 3 至 4.5 mN/m,本量级下比测量误差带来的偏差还大。报告里写明所用模型是基本要求,跨部门或与供应商比对数据时,先对齐模型再比数值。
液体组合与各自的表面张力参数。同一批接触角数据,换成另一套液体参数来源,结果也会有差异。文献里同一种液体的分量取值存在多种版本,二碘甲烷的色散分量在 47.4 至 50.8 mN/m 之间都有引用。工具里的液体参数可以改写,用来对齐自家仪器软件的内置值,这样算出的结果才能和设备读数对得上。
接触角类型。前进角、静止角、后退角在同一个表面上读数不同,对低表面能材料的差异更明显。用 van Oss & Good 做酸碱分量分析时,取前进角更合适,静止角会引入接触角滞后的误差。
读数的有效性。检查有没有混入接近完全铺展的数据点,这类读数的接触角接近 0°,对表面能的约束很弱。计算前先看每个液体接触角的分布,比事后分析异常值更省事。
七、延伸阅读与工具
表面自由能数据通常与前处理效果评价一起使用,需要看测量配置与设备选型时,可以参考表面清洁度与活化效果的测量方案;想了解接触角测量仪的功能与常见场景,可以看另一篇入门介绍。