2. 非球面数学基础:方程、系数与度计算

各位好,我是老张。今天咱们聊聊非球面的数学基础。说实话,刚入行那会儿,我也觉得这些公式看着头疼。但干久了你会发现——不懂方程,你连磨盘都选不对。

非球面加工,说白了就是跟一条曲线较劲。这条曲线长什么样,全靠方程里的几个参数说了算。咱们一个一个拆开讲。

2.1 标准非球面方程

先看最常用的那个公式。几乎所有光学设计软件里,非球面都长这样:

z(r) = (c * r²) / (1 + sqrt(1 - (1 + k) * c² * r²)) + A4 * r⁴ + A6 * r⁶ + A8 * r⁸ + ...

这里:

  • z(r) —— 镜面在径向位置 r 处的矢高(也就是凹陷深度)
  • c —— 顶点曲率(c = 1 / R₀,R₀ 是顶点曲率半径)
  • k —— 二次曲面系数(也叫圆锥系数)
  • A4, A6, A8... —— 高次非球面系数

这个方程分两部分看。前半部分带根号的,描述的是二次曲面(球面、抛物面、双曲面这些)。后半部分是高次项,用来修正边缘形状。

我的习惯:拿到一个非球面图纸,我第一件事就是看 k 值。k 值决定了这个面是偏球还是偏抛物,心里先有个底。

2.2 二次曲面系数 k 的物理意义

k 值决定了基础曲面的类型。我整理了一张表,大家一看就明白:

k 值范围 曲面类型 典型应用
k = 0 球面 传统透镜,最容易加工
k = -1 抛物面 平行光聚焦,望远镜物镜
k < -1 双曲面 卡塞格林系统副镜
-1 < k < 0 椭球面(长轴) 点光源准直
k > 0 扁椭球面 较少见,特殊光路

为什么会这样?你想想看,k 值本质上是控制曲面从中心到边缘的弯曲变化率。k = 0 时,曲率处处相等,那就是球面。k 越负,边缘越"翘"。

注意:我曾经遇到过一位同事,把 k = -0.5 的椭球面当成球面去磨,结果边缘差了十几个微米。加工前一定要确认 k 值,别想当然。

2.3 高次项系数 A4, A6, A8 的作用

光靠 k 值还不够。很多非球面需要在二次曲面的基础上做"微调"。这时候高次项就派上用场了。

每个高次项的作用范围不同:

  • A4(r⁴ 项) —— 主要影响镜面中带区域(约 0.5~0.8 口径处)
  • A6(r⁶ 项) —— 影响中边缘过渡区
  • A8(r⁸ 项) —— 主要影响边缘区域

我个人的经验是:高次项系数每增加一个数量级,加工难度就上一个台阶。因为高次项会让面形在局部区域出现"拐点",抛光时特别容易过修。

避坑指南:我曾经加工过一个 A4 = 1.2e-6 的镜片,按常规工艺走,结果中带区域始终差 0.3 微米。后来发现是 A4 项太大,导致局部曲率变化剧烈,抛光模根本贴不上去。最后换了柔性抛光模才解决。

2.4 非球面度计算

非球面度,是衡量一个面偏离球面多少的指标。说白了就是:这个面有多"非球"。

计算方法很简单:

δ(r) = z_non(r) - z_sphere(r)

其中:

  • z_non(r) —— 非球面在 r 处的矢高
  • z_sphere(r) —— 参考球面在 r 处的矢高(通常取顶点曲率相同的球面)

实际计算时,我们通常取最大非球面度作为工艺难度指标。这个值越大,加工越难。

我给大家一个经验数据:

最大非球面度 加工难度 推荐工艺
< 10 μm 低 传统抛光 + 修正
10 ~ 50 μm 中 数控小工具抛光
50 ~ 200 μm 高 磁流变抛光或离子束
> 200 μm 极高 单点金刚石车削
小技巧:计算非球面度时,参考球面的选择很关键。我习惯用"最佳拟合球面"而不是顶点球面,这样算出来的非球面度更小,工艺评估更合理。

2.5 知识结构总览

下面这张图把本章的核心逻辑串起来了。从方程出发,到系数含义,再到非球面度计算,每一步都指向同一个目标——为工艺方案提供数学依据。

非球面数学基础 · 知识结构 标准非球面方程 二次曲面系数 k 高次项系数 A4, A6, A8 非球面度计算 k = 0 球面 k = -1 抛物面 k < -1 双曲面 A4 影响中带区域 A6 影响过渡区域 A8 影响边缘区域 δ(r) = z_non - z_sphere 最大非球面度 指导工艺选择 核心目标:为工艺方案制定提供数学依据

2.6 实际计算示例

光说不练假把式。咱们拿一个实际案例走一遍。

已知条件:

  • 顶点曲率半径 R₀ = 100 mm
  • 二次曲面系数 k = -0.5
  • 高次项:A4 = 1.0e-8, A6 = 2.0e-12
  • 口径 D = 80 mm(半口径 r_max = 40 mm)

计算步骤:

  1. 先算顶点曲率 c = 1 / 100 = 0.01 mm⁻¹
  2. 代入方程计算各径向位置的矢高 z_non(r)
  3. 取相同 c 的球面作为参考,计算 z_sphere(r)
  4. 求差值 δ(r) = z_non(r) - z_sphere(r)

我算了一下关键位置的结果:

r (mm) z_non (mm) z_sphere (mm) δ (μm)
0 0.0000 0.0000 0.0
10 0.5001 0.5000 0.1
20 2.0020 2.0000 2.0
30 4.5180 4.5000 18.0
40 8.1280 8.0000 128.0
结论:最大非球面度出现在边缘,约 128 μm。按前面的表格,这属于"高难度"级别,推荐使用磁流变抛光或离子束修正。如果硬要用传统抛光,我估计修形周期至少得翻三倍。

嗯,数学基础就讲到这里。这些公式和系数,你以后每个项目都会跟它们打交道。别嫌烦,算清楚再动手,比磨废了再返工强一百倍。