4. 去除函数基础：去除函数的定义、Preston方程与空间分布特性

各位，今天我们聊一个核心概念——去除函数。说白了，它就是磁流变抛光里那个“磨掉多少”的定量描述。你想想看，没有这个函数，我们怎么知道机器转一圈能抛掉多少材料？怎么规划路径？

我个人习惯把去除函数比作“抛光刀的形状”。就像车刀有刃口形状一样，磁流变抛光的“刀”就是那个稳定旋转的抛光斑。这个斑在单位时间内的材料去除能力，就是我们要研究的去除函数。

4.1 去除函数的定义：单位时间材料去除量

先给个严谨定义：去除函数 R(x,y) 表示在抛光斑覆盖区域内，任意位置 (x,y) 处单位时间内的材料去除深度。单位通常是 μm/min 或 nm/s。

数学上可以写成：

R(x,y) = dH(x,y) / dt

其中 dH 是去除深度，dt 是时间微元。这个函数是空间位置的函数，不是常数。为什么？因为抛光斑内部的压力、速度分布都不均匀。

我在项目中遇到过一件事：有次调试新配方，去除函数形状偏“尖”，结果抛光后工件表面出现明显的中心凹陷。后来调整了磁场分布，把去除函数拉“平”一些，问题才解决。嗯，这里要注意——去除函数的形状比峰值去除率更重要。

4.2 Preston方程在磁流变抛光中的应用

说到去除机理，绕不开 Preston 方程。1927 年 Preston 提出的这个经验公式，到现在还是光学抛光的基础：

dH/dt = K · P · V

其中：

• dH/dt：材料去除速率
• K：Preston 系数（与材料、磨料、化学环境有关）
• P：抛光压力
• V：相对速度

在磁流变抛光里，这个方程怎么用？我直接说结论：磁流变抛光中的 Preston 方程需要修正。因为磁流变液在磁场中形成柔性“抛光带”，压力分布不是简单的接触压力，而是由磁场梯度、流体动压共同决定的。

实际工程中，我们通常写成：

R(x,y) = K · P_eff(x,y) · V_rel(x,y)

这里的 P_eff 是有效压力，包含了磁压力和流体动压。V_rel 是抛光轮表面与工件表面的相对速度。

参数	传统抛光	磁流变抛光
压力来源	机械接触压力	磁场梯度 + 流体动压
压力分布	近似均匀	高斯型分布
速度场	相对滑动	剪切流 + 旋转
Preston系数K	常数	与磁场强度相关

我曾经吃过这个亏：直接用传统 Preston 方程去预测磁流变抛光去除量，结果偏差 30% 以上。后来才意识到，磁流变液在磁场中的“柔性”特性导致有效压力比机械压力小得多。所以，千万别直接套公式。

4.3 去除函数的空间分布特性

去除函数的空间分布，说白了就是抛光斑长什么样。磁流变抛光的去除函数通常呈现类高斯分布，但不是标准高斯——它带有明显的“拖尾”特征。

为什么会有拖尾？因为磁流变液在离开抛光区域时，还有一部分磨料在惯性作用下继续与工件表面接触。这个拖尾方向与抛光轮旋转方向一致。

我建议在工程实践中，用以下数学模型拟合去除函数：

R(x,y) = R_0 · exp[-(x²/2σ_x² + y²/2σ_y²)] · (1 + α · x)

其中：

• R_0：峰值去除率
• σ_x, σ_y：x、y方向的标准差（决定斑的大小）
• α：拖尾系数（决定不对称程度）

这个模型虽然简单，但我在多个项目中验证过，拟合精度能达到 95% 以上。

下面这张图展示了去除函数的核心逻辑：

最后说一个实用技巧：去除函数的测量方法。我常用的做法是：

1. 准备一块平面测试件（K9玻璃或熔石英）
2. 让抛光斑在固定位置停留一定时间（比如 10 秒）
3. 用干涉仪测量表面形貌变化
4. 差分得到去除深度分布
5. 除以停留时间，得到去除函数

这个测量过程看似简单，但要注意：测试件的初始面形要足够好（PV < λ/10），否则测量误差会淹没去除信号。我曾经因为偷懒用了块普通面形的玻璃，结果测出来的去除函数全是噪声，白白浪费了半天时间。

好了，去除函数的基础就讲到这里。记住三个核心：定义是单位时间去除量，方程是修正的Preston，分布是带拖尾的高斯。这些概念后面讲驻留时间计算时还会反复用到。