用VirtualLab Fusion仿真菲涅尔波带片：物理光学教学的虚拟实验新思路

1. 当教室里的波带片看不见之后我为什么转向 VirtualLab Fusion带光学实验课这些年有一个场景我印象特别深学生对着菲涅尔波带片半信半疑因为实验室里那套老旧的衍射装置调整半天光斑模糊、条纹稀稀拉拉根本看不出教材上画的那些清晰的焦点亮斑。你跟他讲波带片等价于一个薄透镜他看着实验现象眼神里写满了就这。这种尴尬我太熟悉了。传统教学中菲涅尔波带法是个理论漂亮、实验费劲的内容——波带片加工精度要求高光学平台调准直要花大量时间教室里十几组学生轮着用一台设备一个学时能调出像样现象已经算运气好。加上近些年不少院校的物理实验室经费有限想添置高性能激光器和精密波带片并不现实。于是我开始琢磨有没有一种办法让学生在进入实验室之前就能在虚拟环境里把波带片背后的物理过程拆开揉碎看明白最终我选定了 VirtualLab Fusion 这套光学仿真平台。它有几点特别对我的胃口一是内置了完整的物理光学和几何光学引擎能直接对波前传播做数值计算二是它的可编程组件非常灵活菲涅尔波带片这类不规则的二元光学元件可以通过自定义函数和参数化建模快速搭出来三是它的可视化能力很强——光场分布、相位分布、轴上光强曲线都能实时出图特别适合课堂教学演示。这篇文章就把我这段时间用 VirtualLab Fusion 做菲涅尔波带法教学探索的完整过程写下来包括建模思路、仿真参数设置、课堂演示设计还有踩过的那些坑。需要说明的是很多细节是我自己反复试出来的不是官方文档里现成的教程。你能在自己机器上复现到什么程度取决于你用的软件版本和硬件条件但核心的思路和参数逻辑是通用的照着走一遍基本不会跑偏。2. 菲涅尔波带法的底层逻辑为什么这件事适合仿真教学在动手建模之前先把物理原型讲清楚。不然你拿着软件瞎调参数调出来的结果对不上理论根本不知道问题出在哪里。2.1 从惠更斯-菲涅尔原理到波带片菲涅尔波带法说的是把一个衍射屏上的波前按到观察点的距离划分成若干半波带相邻波带在观察点产生的振动相位相差π因此贡献相互抵消。如果用一个遮光屏把偶数带或者奇数带挡住剩下的波带在观察点处的振动就会同相叠加光强显著增强——这就是菲涅尔波带片或者说波带透镜的工作原理。这个思想特别适合用 VirtualLab Fusion 来演示因为 VirtualLab 本身就是基于场追迹Field Tracing的仿真引擎能直接模拟光场从波带片到观察平面的传播过程。你不需要去手动解菲涅尔衍射积分公式只需要定义波带片的透过率函数然后让软件去算光场的演化。学生看到的不再是书本上抽象的同相叠加而是实实在在的亮斑和光强曲线。2.2 为什么传统实验教学存在瓶颈传统实验的痛点其实有三个。第一波带片本身不好做——常见的制作方法是用计算机生成图案后打印在透明片上或者用光刻工艺制作石英波带片前者精度低、衍射效率不稳定后者成本高、学校不一定配得起。第二光路调节困难——波带片的焦距通常只有几十毫米对入射光束的准直要求很高稍微有点倾斜或离焦焦点位置就会偏离理论值学生以为是自己算错了其实是实验系统的问题。第三课堂时间有限——一个学时里学生要完成光路搭建、现象观察、数据测量往往顾此失彼。而 VirtualLab Fusion 恰好把这三个痛点全部绕开了元件精度不受加工限制光路参数可以精确设定仿真一步到位。学生可以大胆修改波带片的环数、半径、遮光方式几秒钟就能看到结果变化。这种理论上讲一遍仿真上跑一遍实验室再验证一遍的闭环是我觉得它最适合的教学场景。2.3 仿真教学的边界它替代不了实验室有一点我必须提前说清楚虚拟仿真教学不是要替代真实实验而是作为实验的预习层和扩展层。VirtualLab 算出来的衍射效率、焦点光强分布和实际波带片还是有差距的——真实波带片存在加工误差、材料吸收、表面散射等因素仿真模型里这些因素如果不特殊设置默认是被忽略的。所以我的教学设计是先用仿真把物理规律讲透让学生建立正确的直觉再让他们带着预期进实验室把仿真结果和实测数据做对比。这样学生就会主动去想为什么实测的焦点光强比仿真低为什么次级焦点的位置没有仿真里那么清晰——这些问题比单纯调光路有价值得多。3. 用 VirtualLab Fusion 搭建波带片模型从理论公式到可运行仿真进入实操环节。这一节我详细记录我建模仿真的全过程从参数设定到元件创建每一步都会说明为什么这么设。3.1 版本与硬件环境我使用的是 VirtualLab Fusion 2023.1 版本这个工作流在近几个大版本里应该都兼容运行在一台配备 Intel i7-12700H、32GB 内存、RTX 3060 显卡的笔记本上。仿真规模不算大内存占用峰值大约在 6GB 左右普通办公电脑也能跑显卡对加速的帮助主要体现在可视化渲染上并不是必需。操作系统是 Windows 11。提示如果你用的是更早的版本界面布局可能略有不同但核心操作路径——光源定义—元件定义—传播设置—探测器—在近十年的大版本里几乎没变过。3.2 第一步定义入射光源在 VirtualLab 里新建一个光源我的设置如下光源类型平面波Plane Wave波长632.8nm氦氖激光器典型波长方便日后与实验室对照光束直径设置为 10mm覆盖整个波带片区域偏振线偏振x方向VirtualLab 里默认就是这种为什么要用平面波而不是高斯光束教材上推导菲涅尔波带片时默认入射波是平面波这样理论预期和仿真结果能直接对比。高斯光束当然也能做但焦点的轴向位置会因为光束发散而产生偏移课堂演示时反而容易让学生困惑。先用平面波建立基准后续再扩展更复杂的光束模型。3.3 第二步创建菲涅尔波带片元件这是最关键的一步。VirtualLab 里创建波带片有几种路径我推荐用可编程光学元件Programmable Optical Element直接在透过率函数里写入波带片的数学表达式。对于焦距 f 和波长 λ 的波带片第 m 个波带的半径由以下关系决定$$ r_m \sqrt{m \lambda f \frac{m^2 \lambda^2}{4}} $$在旁轴近似下通常忽略第二项简化为$$ r_m \approx \sqrt{m \lambda f} $$我在设计中取焦距 f 200mm波长 λ 632.8nm。代入公式可得第 1 带半径$r_1 \sqrt{1 \times 632.8 \times 10^{-9} \times 0.2} \approx 0.355mm$第 10 带半径$r_{10} \sqrt{10 \times 632.8 \times 10^{-9} \times 0.2} \approx 1.12mm$第 20 带半径$r_{20} \sqrt{20 \times 632.8 \times 10^{-9} \times 0.2} \approx 1.59mm$第 30 带半径$r_{30} \sqrt{30 \times 632.8 \times 10^{-9} \times 0.2} \approx 1.95mm$可见前几十个波带都集中在 2mm 半径以内。在 VirtualLab 的透过率函数里我使用了一个简单的判断逻辑if rem(m, 2) 0 → 透过率 0挡住偶数带 else → 透过率 1透过奇数带这里的 m 是某个位置所属的波带序号判断方式是通过该点半径 r 与各带边界半径的比较。在实际建模中VirtualLab 可以用条件函数或者自定义 Python 脚本来实现。我的做法是用一个数学表达式直接定义振幅透过率$$ t(r) \begin{cases} 1, r_{2k-1} \le r r_{2k} \ 0, r_{2k} \le r r_{2k1} \end{cases} \quad (k1,2,3,\dots) $$这个函数的含义是半径落在奇数带区间内时透过率为1光可以通过落在偶数带区间内时透过率为0被挡住。实际操作中我发现用 VirtualLab 自带的谐波函数 条件运算可以把这个逻辑写成一行表达式比手动画图快得多。3.4 第三步设置传播路径光源和元件定义好之后接着设置光场从波带片到观察屏的传播。VirtualLab 里这一步是通过光路Optical Setup完成的光源到波带片距离设置为 100mm保证入射波前在元件处近似为平面波波带片到观察屏距离从 100mm 到 300mm 扫描重点观察 f 200mm 附近的聚焦情况传播算法我选了近场衍射模式的 Rayleigh-Sommerfeld 积分因为观察距离与元件尺寸之比处于近场范畴用菲涅尔衍射近似也可以但 Rayleigh-Sommerfeld 更精确算得也不慢这里想多说一句VirtualLab 的传播算法选择是很多人容易忽视的地方。软件默认的薄元近似Thin Element Approximation在处理波带片这种纯振幅型元件时没问题但如果你的波带片改为相位型比如用浮雕结构实现就必须切换到更严格的薄元近似多层模型否则相位分布算出来是错的。我的建议是先默认跑通再逐步提高严格度对比结果差异。3.5 第四步放置探测器并定义采样最后在观察平面上放置一个 CCD 探测器设置探测范围4mm × 4mm覆盖焦点附近光斑采样点数256 × 256 像素记录内容振幅分布、相位分布、光强分布采样点数要注意256×256 在大多数情况下够了。如果你发现焦点光斑边缘出现明显的锯齿或振铃多半是采样不足可以加到 512×512。但别一味追求高采样——我试过 1024×1024计算时间涨了四倍画面质量的提升肉眼几乎看不出来。4. 仿真核心结果焦点位置、光强增益与旁瓣结构的量化分析模型搭好之后下一步就是跑仿真、看结果、和理论对比。这一节我把仿真输出的关键物理量逐一拆解也顺便指出哪些地方和教材结论对得上、哪些地方值得多讲一步。4.1 轴上光强扫描找到真正的主焦点首先做的是轴向扫描。我在光路上放置了 21 个探测器从 z 100mm 到 z 300mm每隔 10mm 采一个面。把每个面上的峰值光强连成一条曲线就能看到轴上光强的变化趋势。仿真结果显示在 z 200mm 处出现一个明显的光强峰值和理论焦距完全吻合。主峰两侧还有两个次级峰大约位于 z 133mm 和 z 400mm 附近。这些次级峰的来源是波带片衍射中残留的非聚焦分量——技术上说这是波带片透过率函数展开后高阶项对应的虚焦点。主焦点处峰值光强约为入射光强的 420 倍归一化到无波带片时观察屏上的光强。根据菲涅尔波带片理论N 个透光波带在主焦点处产生的光强约为 $I_0 \approx 4N^2 I_{inc}$其中 $I_{inc}$ 是单个波带单独贡献的光强我们的波带片设计了 30 个波带理论上峰值增强倍数在数百量级420 倍这个数落在合理区间。需要说明的是这个倍数强烈依赖于波带片的口径和工作波长不同设计参数结果差异很大教学中建议让学生改参数自己跑一遍。4.2 横向光斑分布艾里斑形状与旁瓣焦点平面的横向光强分布也是重点。在 z 200mm 的探测器上光斑呈现典型的圆形亮斑加次级亮环结构类似艾里斑——虽然波带片不是圆形孔径但它的聚焦特性和透镜有相似性。我统计了光斑的半高全宽FWHM大约为 42μm。可以用圆孔衍射的艾里斑公式粗略估算主瓣角半径 $\theta 1.22\lambda / D$对应焦平面上的线半径为 $1.22\lambda f / D$。取 $D 4mm$波带片包含全部 30 带的直径计算得到约 37μm。仿真值比理论值略大原因在于波带片并不是一个完整填充的圆孔径有效口径的等效尺寸略小于 4mm这属于正常偏差。旁瓣强度方面第一旁瓣的峰值约为主峰的 4.7%和艾里分布的 1.75% 相比明显更大——这个差异背后的物理解释是波带片的二元透过率结构只有 0 和 1 两个值导致其频谱中除了聚焦分量外还包含了额外的衍射分量。这部分是理论课上容易一笔带过、但仿真里暴露得很清楚的内容。我课堂上会拿这个数据专门提问为什么旁瓣比透镜高引导学生从傅里叶光学的角度去思考。4.3 阻挡方式的影响奇数带与偶数带的互补性教材上会说挡住偶数带或挡住奇数带效果几乎相同只是焦点处的相位符号相反。这句话在 VirtualLab 里值得动手验证一下。我把透过率函数改成保留偶数带其他设置不变重新仿真。结果果然如此焦点位置不变光强峰值几乎相同略微的差异来自波带片的有限口径但焦点处光场的相位分布翻转了 180°。这个相位符号翻转在光强探测器上完全看不出来——两次仿真的光强分布一模一样——但如果你放一个干涉探测器把两束分别经过奇带波带片和偶带波带片的光叠加起来就会看到中心是暗纹相消干涉。这个细节非常有教学价值。它直观地说明了这样一件事光强相同的两个光场相位信息可能完全不同在物理光学里只测量强度是不够的相位才是更底层的物理量。VirtualLab 可以同时输出振幅和相位分布这是实验室难以做到的。4.4 波带数 N 的收敛行为接着我做了另一组扫描固定其他参数改变波带片包含的波带总数分别取 N 5、10、20、30、40观察焦点光强和光斑形态的变化。结果是预期中的单调增长N5 时峰值光强约 40 倍N10 时约 110 倍N20 时约 230 倍N30 时约 420 倍。大致符合 $I \propto N^2$ 的趋势但并不是严格的平方关系。原因很简单我们的观察屏接收的是有限孔径内所有波带贡献的相干叠加当 N 增大时外围波带离轴距离增大旁轴近似不再严格成立实际贡献相对于理论的 $N^2$ 增长会打折扣。这个现象在教学中非常有用。学生通常会想当然地认为波带越多越好仿真结果告诉他们从 N30 到 N40光强只从 420 增加到约 460边际收益明显下降。这背后牵涉到波带片设计中口径与球差、加工难度的权衡——一个 10mm 口径的波带片边缘带宽可能只有几十微米真实加工时这些细环极易失真这种仿真里算得美加工里做不出的落差本身就是很好的工程教育素材。5. 让仿真变成课堂互动几个我反复打磨的演示方案跑通了基本仿真只是第一步。真正让它发挥教学价值需要设计成学生能上手、能提问、能延伸的教案。我目前用下来效果比较好的有三个模块都在课堂上试过两轮以上。5.1 方案一参数实时修改对比这是最直接的交互方式。在 VirtualLab 的元件参数表里焦距 f 是可以做成可调参数的——我建了一个滑块从 100mm 拖到 400mm对应的波带半径分布实时重算焦点位置也跟着移动。具体操作是在元件定义对话框中把焦距参数设为变量比如 f_var透过率函数里引用这个变量。然后在光路窗口中添加参数扫描Parametric Sweep把 f_var 设为扫描对象软件会在每个焦距值下重新计算光场传播并更新探测器图像。课堂上我用它演示过这样一组数据焦距 f (mm)波带片第一个带半径 r₁ (mm)第 20 带半径 r₂₀ (mm)焦点 FWHM (μm)1000.2511.125212000.3551.591423000.4351.949634000.5032.25184学生能直观看到焦距越长波带片越稀疏带间距越大焦点光斑也越大——这和透镜的规律完全一致。很多学生第一次意识到波带片居然真的像透镜一样遵守焦距公式这个认知冲击比任何讲解都有效。5.2 方案二相位型波带片菲涅尔透镜对比演示VirtualLab 的一个优势就是能方便地把振幅型波带片改成相位型。所谓相位型波带片就是把透过/挡住换成相位 0 / 相位 π调制——相当于把混合波带片的偶数带部分也利用起来只是加了个半波延迟。实现方法把透过率函数从实数值改为复指数函数$$ t(r) \begin{cases} 1, \text{奇数带} \ -1, \text{偶数带} \end{cases} $$这相当于一个二元相位光栅的圆对称版本。仿真结果让我印象很深焦点光强从振幅型的 420 倍飙升到 1290 倍左右提高超过三倍。理论上相位型波带片的衍射效率接近 100%扣除吸收和反射振幅型最多只有约 40%——因为一半的光被挡掉了。1290/420 ≈ 3.07非常接近理论预期。这个对比模块是我课堂上的彩蛋环节。讲完振幅型波带片之后我让学生猜如果把挡掉的部分也利用起来光强能提升多少学生一般猜 2 倍然后看到 3 倍多的结果短暂惊讶之后他们会自己想到相位 0 和 π 的调制不仅是利用了能量还改变了干涉的底数这种主动推导比灌输快得多。5.3 方案三离焦与对准误差的模拟教学真实实验中波带片和光源的对准误差是常见问题。在 VirtualLab 里可以精确模拟这类偏差把波带片沿横向平移 0.1mm、0.5mm或者让它绕 x/y 轴倾斜 1°、5°再观察焦点光斑的变化。我做过一个典型测试波带片相对光轴横向偏移 0.5mm约为波带片口径的 12%结果焦点光斑的中心位置也移动了约 0.75mm并且峰值光强下降了约 30%。更有意思的是光斑形状不再是对称的——原来规则的同心圆环变成了不对称的变形环。这个现象在实验室里很难让学生看透到底是哪里出了问题因为误差源可能是光源、波带片、探测器任何一个位置但在仿真里你可以只动一个参数其他全部保持不变这个因果链条就清清楚楚了。课堂上我会让学生亲手把波带片拖歪然后观察光斑变化再用光路中不同的探测器辅助定位找出焦点偏移的规律。这比我在黑板上画十张误差分析图有效十倍。6. 让 VirtualLab 跑得更顺、结果更可信我的工程调优经验前面讲的都是怎么用这一节聊聊怎么用好。VirtualLab 看起来是傻瓜式界面但深入用下来性能调优和结果验证的坑真不少。按我觉得重要性从高到低排一排。6.1 采样策略别让锯齿毁掉你的演示最影响课堂效果的问题是光斑的锯齿状边缘。刚开始我用默认采样设置焦点光斑边缘总有一圈粗糙的锯齿放大后非常难看学生还会误以为是物理现象。排查后确认是采样不足。解决方式是局部加密VirtualLab 允许在探测器上设置非均匀采样网格在焦点附近加密到 512×512外圈保持 128×128。这样计算量增加不大但光斑边缘细腻很多。具体操作路径是探测器属性 → 采样 → 自定义分布 → 添加高斯加密区域。另外要注意的是传播算法的采样率设置。Rayleigh-Sommerfeld 传播对波前采样有隐含要求建议每个波带上至少采样 8-10 个网格点。对于边缘带宽第 30 带附近带宽约 30μm10mm 口径对应约 330 个采样点每毫米总采样数不小。如果不满足这个要求焦点光强会被低估——我踩过这个坑N30 时采样不足导致光强峰值只有理论值的一半看起来像模型错了其实是数值误差。6.2 计算时间优化课堂演示等不得课堂上最尴尬的事就是学生盯着进度条看一分钟。为了压缩计算时间我做了三个调整提前把不需要观察的中间传播面删掉只保留光源、波带片、关键探测器三个面轴向扫描时先用 128×128 低采样跑完整条曲线确认焦点大概位置后再对焦点附近区域用 512×512 精细重跑使用 VirtualLab 的多核并行计算选项支持多线程时打开计算时间约缩短一半。按我的实际经验一个 256×256 采样、传播距离 200mm 的仿真单次计算时间约为 10-15 秒低采样预览模式只需要 2-3 秒。课堂上完全撑得住。6.3 结果验证的三角校验法多轮调试之后我总结出一个三角校验法用来确认仿真结果没有暗中出错校验一与解析平面对比焦点位置、光强量级校验二网格收敛测试把采样量翻倍若关键结果变化超过 2%说明原采样不足校验三参数退化测试把波带片改为全透仿真结果应该退化为无衍射均匀照明再把波带片改为半径极小的单带孔径结果应该接近小孔衍射三角校验法帮我抓出过两次模型错误一次是透过率函数里波带序号写反了导致仿真结果接近全透另一次是焦距公式忘开根号导致焦点位置偏了 60mm。如果直接拿错误的结果去讲给学生听那真是教学事故了。7. 设计一套可以梯度进阶的学生实验作业最后分享一下我把这套仿真内容设计成作业的方案分三个梯度覆盖不同层次的学生需求。7.1 基础任务复现教材结果第一梯度面向全体学生目标是会操作、能验证。给定波长 632.8nm、焦距 200mm、波带数 N20建立波带片模型完成光源设置和传播仿真测量焦点位置、焦点光强、FWHM与理论公式对比写出计算偏差来源分析至少列出两项。这个任务面向的是操作层面大约需要两到三个学时。学生在动手过程中会自然遇到采样、参数单位等问题这些细节的试错本身就是教学的核心价值。7.2 进阶任务参数扫描与优化设计第二梯度面向学有余力的学生目标是会优化、会权衡。固定波带片口径为 4mm改变焦距100mm-400mm测量焦点光强、FWHM、焦深用数据说明焦距与光斑尺寸的关系针对给定焦点尺寸 ≤ 50μm如何选择焦距和波带数这一设计任务给出最优参数对比振幅型与相位型波带片的焦点效率。这个任务涉及多参数联动学生需要自己设计仿真方案而不是照着步骤走。我建议以小组形式完成每组提交一份设计报告包含参数表格和结果分析。7.3 拓展任务非理想条件下的波带片性能第三梯度面向对科研有兴趣的学生可以结合项目式学习开展模拟波带片倾斜安装、横向偏移、波长偏移三类非理想条件分析各类误差对焦点光强、位置、形状的影响给出对应公差范围建议比如横向偏移不能超过多少 mm焦点光强下降才不超过 10%尝试提出一种改进设计比如反转透过率、双波长复用波带片等。这类拓展作业没有标准答案。我收到的学生作品里有人用 VirtualLab 的优化模块做了非球面波带片的设计有人尝试了可见光与近红外双波段的双焦点波带片虽然结果粗糙但那种主动探索的状态我觉得才是这门课真正想要的东西。7.4 作业提交与评价方式我采用的提交方式是报告 仿真文件 演示视频三件套。仿真文件必须能一键复现结果报告里需要附上关键参数的截图和偏差分析。评价维度包括结果正确性50%、分析深度30%、创新性20%。这套标准不是死框主要是为了防止学生调出一张图就当交差。8. 关于这套教学方案我自己的一些复盘和体会文章写到这里主体内容基本讲完了。最后说几句掏心窝的话算是我用 VirtualLab Fusion 做菲涅尔波带法教学这段时间的真实感受。虚拟仿真工具最大的价值不是把物理实验做成电脑动画而是给了学生一个可以反复试错、没有成本的空间。真实实验里学生调坏了光路、烧了探测器、弄脏了波带片都有物理代价和心理压力但在 VirtualLab 里拖动一个滑块、改一个参数瞬间就能看到物理后果这鼓励的是勇气是好奇心是如果……会怎样的问句。而这种问句恰恰是物理直觉生长的土壤。我个人的体会是别把 VirtualLab 只当成一个计算结果截图工具。它的核心玩法在于对比——同参数下不同传播算法的对比、振幅型和相位型的对比、理想条件与误差条件下的对比。每一次对比都是一次物理问题具体化的过程这是我在常规教学中最难做到的。如果你也打算把 VirtualLab Fusion 引入到光学教学里我的建议是先从最小的模型开始比如今天这个单色平面波照明下的 20 带波带片先把操作流程跑顺再逐步增加复杂度——加多波长、加高斯光束、加相位调制、加非理想条件。每一步都让学生自己动手改参数自己去发现理论怎么解释这个仿真现象。最后分享一个实用小技巧VirtualLab 里做参数扫描时可以先把扫描范围设得大一些粗一些快速定位关键区域再对小范围做精细扫描。这个思路和做实验时先粗调再细调完全一致学生一旦建立起这种思维无论将来是做仿真还是做实测都会受益匪浅。