4、面型设计的基本原理：折射与反射定律、斯涅尔定律、能量守恒与边缘光线原理
目标照度分布反求面型（Inverse Design of Freeform Optics）的核心，是回答一个工程问题：给定光源的光强分布和靶面上的目标照度分布，如何确定光学元件（透镜或反射镜）的面型坐标？ 这一问题的数学本质是求解一个非线性偏微分方程（Monge-Ampère 型方程），而物理本质则完全由本章阐述的四条基本原理支配。任何自由曲面设计软件（如 LightTools、Zemax OpticStudio 或自研算法）的求解结果，都必须满足这些定律，否则设计出的面型在物理上不可实现。本章将逐一拆解这些原理，并说明它们如何协同作用，将“光线方向”与“能量分配”两个维度耦合在一起。
上图展示了四条基本原理在自由曲面设计中的耦合关系。光源发出的光线经过面型折射或反射后到达靶面，其中斯涅尔定律决定了光线在面型上的偏折方向，能量守恒决定了光线从光源到靶面的映射关系，边缘光线原理为求解提供了边界条件，而面型的梯度与曲率则直接决定了最终的照度分布形态。下面逐一展开。
4.1 斯涅尔定律与折射/反射定律：光线方向的“交通规则”
4.1.1 核心原理：从矢量形式到面型梯度
斯涅尔定律（Snell's Law）是几何光学中最基本的定律，它描述了光线在两种介质界面上的偏折行为。对于折射，其标量形式为：
n₁ · sin(θ₁) = n₂ · sin(θ₂)其中 n₁、n₂ 分别为入射介质和出射介质的折射率，θ₁ 为入射角，θ₂ 为折射角。在自由曲面设计中，我们更常使用其矢量形式：
n₂ · Ô = n₁ · Î + [n₂·cos(θ₂) - n₁·cos(θ₁)] · N̂式中 Î 为入射光线单位矢量，Ô 为出射光线单位矢量，N̂ 为面型法向单位矢量。对于反射面（n₁ = n₂ = 1），该式退化为：
Ô = Î - 2(ηN̂)N̂在自由曲面设计中，面型通常表示为 z = f(x, y)，其法向量 N̂ 与面型梯度 ∇f 的关系为：
N̂ = (-∂f/∂x, -∂f/∂y, 1) / √[(∂f/∂x)² + (∂f/∂y)² + 1]这意味着，只要确定了面型上每一点的梯度，就确定了该点的法向量，进而通过斯涅尔定律确定了出射光线的方向。因此，面型设计问题可以转化为求解一个关于梯度 ∇f 的偏微分方程。
4.1.2 判断方法：如何验证面型是否满足斯涅尔定律
在实际工程中，我们通常通过光线追迹来验证设计结果。具体步骤如下：
- 提取面型数据：从设计结果中导出面型点云或 NURBS 曲面参数。
- 计算法向量：对每个采样点，通过数值微分（如中心差分）计算 ∂f/∂x 和 ∂f/∂y，进而得到法向量 N̂。
- 计算入射方向：根据光源位置和面型点坐标，计算入射光线单位矢量 Î。
- 计算理论出射方向：代入斯涅尔定律矢量公式，得到理论出射方向 Ô_theory。
- 对比实际出射方向：在光学软件中追迹光线，得到实际出射方向 Ô_actual，计算两者夹角误差。
通常，当所有采样点的角度误差小于 0.1° 时，认为面型满足斯涅尔定律。若误差较大，则说明面型梯度计算有误或曲面光顺性不足。
4.2 能量守恒定律：光源与靶面之间的“账本”
4.2.1 核心原理：从光通量守恒到映射关系
能量守恒定律在光学设计中表现为：光源发出的总光通量必须等于靶面接收的总光通量（忽略光学元件的吸收、散射和反射损失）。数学表达式为：
∫∫_Ω I(θ,φ) dΩ = ∫∫_A E(x,y) dA其中 I(θ,φ) 为光源的发光强度分布（单位：cd），dΩ 为立体角微元；E(x,y) 为靶面上的照度分布（单位：lx），dA 为面积微元。
更重要的是，能量守恒还要求光源的每一个立体角微元 dΩ 与靶面的每一个面积微元 dA 之间存在一一对应的映射关系。这个映射关系是自由曲面设计的核心未知量。设光源出射方向为 (θ,φ)，靶面落点坐标为 (x,y)，则映射关系可表示为：
(x, y) = F(θ, φ)自由曲面设计的任务，就是找到一个面型 z = f(x,y)，使得通过该面型折射/反射后，光线按照映射 F 从光源到达靶面，同时满足斯涅尔定律。
4.2.2 操作流程：建立能量映射的数值方法
在实际求解中，我们通常采用以下步骤建立能量映射：
- 离散化光源：将光源的发光空间划分为 M 个立体角网格，每个网格的光通量为 ΔΦᵢ = I(θᵢ,φᵢ)·ΔΩᵢ。
- 离散化靶面：将靶面划分为 N 个面积网格，每个网格的目标光通量为 ΔΦⱼ = E(xⱼ,yⱼ)·ΔAⱼ。
- 建立映射：根据能量守恒，将 M 个光源网格按光通量比例映射到 N 个靶面网格。常用的方法包括：
- 等光通量映射：将光源和靶面都按等光通量划分，然后一一对应。
- 数值优化映射：通过求解最优传输问题（Optimal Transport）获得映射，使光线传输距离最短或变形最小。
- 求解面型：在获得映射关系后，结合斯涅尔定律，通过数值方法（如龙格-库塔法、牛顿迭代法）求解面型坐标。
4.2.3 案例拆解：LED 均匀照明透镜设计
假设我们需要设计一个用于 LED 光源的均匀照明透镜，LED 的发光强度分布近似为朗伯体：I(θ) = I₀·cos(θ)，其中 θ 为与光轴的夹角，最大发光半角为 60°。目标是在距离光源 100mm 的靶面上形成直径 200mm 的均匀圆形光斑。
关键输入数据：
- 光源总光通量：Φ_total = 100 lm
- 靶面半径：R = 100 mm
- 靶面平均照度：E_avg = Φ_total / (π·R²) = 100 / (π·0.01) ≈ 3183 lx
- 光源发光半角：θ_max = 60°
分析判断过程：
首先，将光源的发光空间按等光通量划分为 100 个环带。由于朗伯体的光强分布已知，每个环带的光通量可以通过积分计算。然后，将靶面按等面积划分为 100 个同心圆环，每个圆环的光通量相等（因为目标照度均匀）。最后，建立光源环带与靶面圆环的映射关系：光源的第 i 个环带对应靶面的第 i 个圆环。
实际操作步骤：
- 计算光源第 i 个环带的光通量：ΔΦᵢ = 2π·I₀·∫_{θᵢ₋₁}^{θᵢ} cos(θ)·sin(θ) dθ
- 计算靶面第 i 个圆环的面积：ΔAᵢ = π·(rᵢ² - rᵢ₋₁²)
- 根据能量守恒 ΔΦᵢ = E_avg·ΔAᵢ，求解每个环带对应的靶面半径 rᵢ。
- 根据映射关系 (θᵢ, rᵢ)，结合斯涅尔定律，逐点求解透镜面型。
结果对比：设计完成后，在光学软件中追迹 100 万条光线，得到靶面照度均匀度（最小照度/平均照度）为 92%，满足一般照明要求（>85%）。
风险与代价：由于 LED 光源并非理想点光源，实际光源有一定尺寸，导致边缘光线模糊，照度均匀度略有下降。此外，透镜边缘厚度较薄，注塑成型时可能出现填充不足的问题。
4.3 边缘光线原理：边界条件的“锚点”
4.3.1 核心原理：从边缘到内部的映射策略
边缘光线原理（Edge Ray Principle）指出：光源的边缘光线（即最大发射角的光线）经过光学系统后，应到达靶面的边缘。这一原理在非成像光学中具有重要地位，它保证了光学系统的扩展量（Etendue）守恒，并且为面型求解提供了明确的边界条件。
在自由曲面设计中，边缘光线原理的作用体现在两个方面：
- 确定面型的边界条件：光源边缘光线对应的靶面落点，决定了面型的边界位置和边界处的梯度。
- 保证光线不交叉：如果边缘光线正确映射到靶面边缘，且内部光线按能量守恒单调映射，则光线不会发生交叉，照度分布可控。
4.3.2 判断方法：如何验证边缘光线原理
在光学软件中，可以通过以下方法验证设计是否满足边缘光线原理：
- 在光源处设置一个环形光阑，只允许最大发射角附近的光线通过。
- 追迹这些光线，观察它们在靶面上的落点分布。
- 如果落点集中在靶面边缘附近，且与设计目标一致，则说明边缘光线原理得到满足。
若边缘光线落点偏离靶面边缘，则说明面型边界处的梯度计算有误，需要重新调整边界条件。
4.3.3 实战场景：汽车前照灯近光透镜设计
在汽车前照灯设计中，近光需要形成明暗截止线（Cut-off Line），即靶面上方为暗区、下方为亮区。这一要求可以通过边缘光线原理实现：将光源的上边缘光线映射到靶面的明暗截止线位置，而光源的下边缘光线映射到靶面的最下方。
具体设计约束：
- 光源：LED 发光面尺寸 1mm × 1mm，发光半角 90°
- 靶面：距离透镜 25m 处的垂直平面，明暗截止线高度为 0.7°（约 0.3m）
- 透镜材料：PMMA（n = 1.4935）
分析判断过程：首先，根据边缘光线原理，确定光源上边缘光线（θ = +90°）经过透镜后应到达靶面明暗截止线位置（y = 0.3m）。然后，根据斯涅尔定律，反推透镜上边缘点的法向量。接着，逐步向内求解，直到光源下边缘光线（θ = -90°）映射到靶面最下方（y = -3m）。
实际操作步骤：
- 建立坐标系：光源位于原点，透镜光轴沿 z 轴，靶面位于 z = 25m 处。
- 将光源发光空间按等光通量划分为 200 个角度区间。
- 根据边缘光线原理，设定边界条件：θ = +90° → y = 0.3m；θ = -90° → y = -3m。
- 使用数值求解方法（如四阶龙格-库塔法），从边缘向中心逐点求解透镜面型。
结果对比：设计完成后，在光学软件中模拟，得到明暗截止线的清晰度为：暗区照度 < 0.1 lx，亮区平均照度 > 20 lx，满足法规要求。
复盘经验：边缘光线原理在非成像光学中具有普适性，但需要注意：当光源有一定尺寸时，边缘光线并非严格意义上的“一条光线”，而是一个“边缘光锥”。此时，需要将光源边缘的每个点都视为边缘光线的起点，这增加了设计的复杂度。
4.4 面型如何影响光线方向与照度分布
4.4.1 面型梯度决定光线方向
面型上某一点的梯度 ∇f = (∂f/∂x, ∂f/∂y) 决定了该点的法向量方向，进而通过斯涅尔定律决定了出射光线的方向。具体来说：
- 梯度越大，法向量与光轴的夹角越大，光线偏折越剧烈。
- 梯度方向决定了光线在靶面上的落点方向。例如，若面型在 x 方向有正梯度，则光线向 x 正方向偏折。
因此，面型设计本质上是对梯度场的规划。在数值求解中，我们通常先求解梯度场，再通过积分重建面型。
4.4.2 面型曲率决定照度密度
面型的曲率 κ 决定了光线的汇聚或发散程度。在靶面上，照度 E 与光线密度成正比。如果面型在某区域曲率较大（凸面），则光线汇聚，局部照度升高；如果曲率较小（凹面），则光线发散，局部照度降低。
定量地，对于旋转对称系统，靶面照度与面型曲率的关系可近似表示为：
E(r) ∝ I(θ) / (r · |dr/dθ|)其中 r 为靶面径向坐标，θ 为光源出射角，dr/dθ 反映了面型对光线的“拉伸”或“压缩”程度。当 dr/dθ 较大时，光线在靶面上分布较稀疏，照度较低；反之，照度较高。
4.4.3 面型连续性对照度分布的影响
面型的连续性直接影响照度分布的连续性：
- C⁰ 连续（位置连续）：保证面型没有断裂，但可能出现尖角，导致照度分布出现突变。
- C¹ 连续（切线连续）：保证法向量连续，光线方向连续变化，照度分布平滑。
- C² 连续（曲率连续）：保证曲率连续，照度分布的二阶导数连续，视觉效果最佳。
在实际制造中，通常要求面型至少达到 C¹ 连续，否则在注塑或模压成型时会产生应力集中，且光学性能不稳定。
4.4.4 常见误区与总结复盘
| 误区 | 后果 | 正确做法 |
|---|---|---|
| 只关注斯涅尔定律，忽略能量守恒 | 光线方向正确，但照度分布与目标不符 | 同时求解方向与能量映射 |
| 映射关系不单调 | 光线交叉，照度出现“暗斑”或“亮斑” | 保证映射的雅可比行列式恒正 |
| 面型曲率变化过大 | 制造困难，光学效率下降 | 添加曲率约束或光顺项 |
| 忽略边缘光线原理 | 靶面边缘照度不可控，出现“漏光” | 将边缘光线作为边界条件 |
总结复盘：自由曲面面型设计的本质，是在满足斯涅尔定律（方向约束）和能量守恒（能量约束）的前提下，通过边缘光线原理确定边界条件，最终求解面型梯度场并积分得到面型。面型的梯度决定了光线的方向，面型的曲率决定了照度的密度，面型的连续性决定了照度的平滑度。四者缺一不可，共同构成了目标照度分布反求面型的理论基础。