做仿真最怕什么最怕换一个尺寸整个模型推倒重来。我之前用Lumerical FDTD Solutions做光学仿真一开始也是老老实实GUI里画结构鼠标点得飞起。直到有一次要扫20组不同宽度的波导每组还得改两次位置我在界面上整整折腾了一下午改到后面眼都花了才痛下决心把脚本捡起来好好学。学完之后只有一个感觉如果早半年学会之前那些无效加班的时间足够我多刷两部剧了。这篇东西就是针对“用FDTD脚本创建结构”这个主题把建模阶段最核心的脚本操作做个完整的梳理。内容不绕弯子全部是实操向的从最基础的语法讲起到怎么建复杂结构、怎么做参数化扫描再到我踩过的那些坑一整套都会覆盖。适合刚接触FDTD脚本、或者GUI建结构建到想骂人的朋友看完基本能上手干活。1. FDTD脚本建模到底解决什么问题1.1 GUI建模的三大痛点你肯定遇到过先说个真实对比。用GUI建一个简单的矩形波导确实不难选中结构按钮一拖就行。但一旦碰上下面这几种情况GUI就开始让人抓狂了。第一重复性操作。比如你有一个周期光栅周期是500纳米总共要200个周期你用GUI画一个矩形再复制199次试试不是不能做但是每次复制都要手动设置偏移量、检查有没有重叠、确认位置对不对稍不注意就错一个。第二结构尺寸需要参数化。很多时候你不是只仿一个尺寸而是要扫描一组尺寸看看某个物理量随尺寸的变化趋势。GUI里画的结构虽然也能设变量但变量的定义、传递、循环扫描逻辑最终还是得靠脚本完成。你总不能一组一组手动改完再手动跑仿真。第三可复现性。GUI建模有个致命伤——你今天拖出来的结构明天想恢复一模一样的很难还原到像素级一致。但脚本不一样代码摆在那里随时跑随时得到一模一样的结果。这一条对学术研究尤其重要审稿人万一要求你复现数据脚本直接甩过去比截图有说服力得多。1.2 表面是省时间本质是结构化思维刚开始用脚本的时候我只是觉得它省事。但后来我意识到它的价值远不止“省事”这么简单。GUI建模的思维方式是“点哪里、拖哪里”跟着感觉走。而脚本建模逼着你先想清楚我的结构有哪几部分每部分的关键尺寸是什么这些尺寸之间有没有关联结构的几何关系是怎样的这个思考过程本身就是一次模型逻辑的梳理。用脚本建出来的模型不只是几何堆砌而是一套可调节、可扩展的参数系统。后续要做优化、要做容差分析、要跑批量仿真都是水到渠成的事。所以哪怕你只是做一个最简单的结构我也建议你从脚本入手哪怕只是用脚本建一个矩形也要让自己习惯这种建模方式。2. 脚本创建结构的基础语法与核心命令2.1 Lumerical脚本的几个底层逻辑Lumerical FDTD的脚本语言语法上跟MATLAB有很多相似之处核心是面向对象的操作方式——每个对象结构、监视器、光源都是一个对象通过set和get这两个函数来设置和获取属性。我第一次写脚本的时候最大的困惑就是我该从哪里知道一个结构有哪些属性答案其实很简单在GUI里选中一个结构属性编辑器里显示的每一项就是它对应的属性名。脚本里的set(属性名, 值)本质上就是在程序化地操作属性编辑器里的每个字段。另外要注意Lumerical的脚本里单位系统默认为米。你设置一个宽度为500纳米的波导要写500e-9而不是写500。这个转换初学者经常忘记导致结构尺寸差了十万八千里。2.2 高频结构创建命令速查先澄清一个概念Lumerical建模中你创建的结构可以分成两大类。一类是几何实体比如矩形、圆、多边形它们是有物理形状的另一类是材料对象它们是把某个材料属性赋给某个几何区域。实际使用中这两类常常是同时出现的。先看最基础的结构创建命令族。常用的有addrect添加矩形addcircle添加圆addpoly添加多边形addring添加圆环addcylinder添加圆柱以addrect为例一个最基础的矩形可以这样写addrect; set(name, waveguide); set(x, 0e-6); set(y, 0e-6); set(z, 0e-6); set(x span, 2e-6); set(y span, 0.5e-6); set(z span, 0.22e-6); set(material, Si (Silicon) - Palik);这段代码做了几件事创建一个矩形命名为waveguide位置在原点尺寸是长2微米、宽0.5微米、高0.22微米材料是硅Palik数据。这里有一个很关键的点设置完几何尺寸后一定要设置材料。如果不设置材料默认是一个无材料的结构仿真时不会产生任何物理效果。另外addrect之后新添加的结构会自动成为“当前选中对象”后续的set操作都是作用在这个结构上。如果你想切换到别的结构进行操作需要用select(name)来指定。这段代码有没有可以简化的地方Lumerical脚本还支持另一种写法把set操作合并到addrect的括号参数里。但对于新手来说分行写更清晰排查问题的时候也更容易定位。3. 从简单到复杂手把手搭一个硅波导耦合结构3.1 先搭底层基底、波导、包层理论说了半天不如动手来一个完整案例。我挑一个硅光领域非常常见的结构绝缘体上硅SOI条波导上面覆盖一层二氧化硅包层。这个结构虽然简单但覆盖了多结构创建、材料设置、坐标设定等核心操作。第一步先创建基底。基底通常是厚层的二氧化硅尺寸要比波导大不少避免仿真区域边界对结果产生影响。addrect; set(name, substrate); set(x, 0e-6); set(y, 0e-6); set(z, -1e-6); set(x span, 10e-6); set(y span, 10e-6); set(z span, 2e-6); set(material, SiO2 (Glass) - Palik);这里把z设在-1e-6基底厚度设为2微米这样基底上表面正好在z 0的位置方便后续在上面叠加波导结构。基底搞定后第二步创建硅波导。此时波导的底面要和基底上表面接触所以波导的z坐标要设为0e-6加上波导自身高度的一半。addrect; set(name, si_waveguide); set(x, 0e-6); set(y, 0e-6); set(z, 0.11e-6); set(x span, 5e-6); set(y span, 0.5e-6); set(z span, 0.22e-6); set(material, Si (Silicon) - Palik);注意这里z设为0.11e-6是因为波导中心高度是0.22e-6的一半即0.11微米。这样波导底面恰好位于z 0e-6也就是基底上表面。这个“中心坐标加半高”的计算是初学者最容易搞混的地方。第三步创建包层。包层覆盖整个波导结构厚度不需要太厚够用就行但是要完全覆盖波导。addrect; set(name, cladding); set(x, 0e-6); set(y, 0e-6); set(z, 0.6e-6); set(x span, 10e-6); set(y span, 10e-6); set(z span, 1e-6); set(material, SiO2 (Glass) - Palik);三层结构搭完硅波导就被包裹在二氧化硅中间了。跑一下仿真能看到基本的模式分布。3.2 循环批量生成光栅阵列也能一键搞定单结构创建会了接下来解决“重复性操作”的痛点。假设我们要建一个均匀光栅周期为700纳米占空比50%总共20个周期。如果用GUI一个一个复制没半小时搞不定。用脚本几行代码的事。思路是在一个for循环里计算每个齿的位置和尺寸然后创建矩形。# 清除之前的结构 deleteall; # 参数定义 period 700e-9; duty_cycle 0.5; num_periods 20; tooth_width period * duty_cycle; etch_depth 0.07e-6; x_start -num_periods * period / 2; for (i 0; i num_periods; i) { x_center x_start i * period tooth_width / 2; addrect; set(name, grating_tooth_ num2str(i)); set(x, x_center); set(y, 0e-6); set(z, 0.11e-6); set(x span, tooth_width); set(y span, 0.5e-6); set(z span, etch_depth); set(material, Si (Silicon) - Palik); }这里有一个关键细节周期是700纳米占空比50%所以每个齿的宽度是350纳米。齿与齿之间的间距也是350纳米。循环变量i从0到19每个齿的中心坐标通过x_start i * period tooth_width / 2计算出来。这段逻辑很直白但有几个细节我特别提醒一下。第一变量名不要和Lumerical的内置函数重名比如不要用length、size这类名字否则很容易出诡异问题。第二循环里每次addrect之后都要重新set(name, ...)否则循环结束所有结构都叫同一个名字后面想select某个具体的齿就麻烦了。如果你想做非均匀光栅比如周期从700纳米渐变到800纳米也只需要在循环里把周期设为一个变量让它随着i递增即可。脚本建模的灵活性在这个场景下体现得淋漓尽致。4. 参数化建模与自动化扫描让脚本真正替你干活4.1 变量传递与表达式计算脚本建模的杀手锏在于参数化。你不需要把每个尺寸都写成死数字而是可以定义变量、用表达式计算尺寸然后用这些变量去驱动多个结构。举一个最常见的场景波导宽度变化对有效折射率的影响。你想扫一组宽度从300纳米到500纳米步长20纳米。如果手动建模每组宽度都要重建模型、重新设置监视器、重新跑仿真工作量惊人。但用脚本你只需要在最外层套一个参数扫描循环。在动手写扫描之前要先明白Lumerical脚本的一个特性FDTD仿真区域、光源、监视器这些对象和结构对象一样都可以用set命令来配置。所以参数化不只是针对结构尺寸光源波长范围、仿真区域大小、监视器位置……所有的一切都可以用变量来驱动。以波导宽度扫描为例核心思路是在循环里修改波导结构的y span属性然后重新运行仿真再提取结果。这样每一次循环就是一次完整的“建模-仿真-提取数据”流程。4.2 完整的自动扫描脚本框架下面这种写法是我自己的习惯经过多次验证稳定可靠# 扫描参数定义 widths linspace(300e-9, 500e-9, 11); results matrix(length(widths), 2); for (i 1:length(widths)) { # 选择波导并修改宽度 select(si_waveguide); set(y span, widths(i)); # 运行仿真 run; # 提取有效折射率 select(index_monitor); neff getdata(index_monitor, neff); results(i, 1) widths(i); results(i, 2) neff; } # 保存结果 write(neff_vs_width.mat, results);这段代码有三个地方特别重要。第一select(si_waveguide)之后再用set修改属性这确保了每次循环操作的都是同一条波导。如果你在循环里用addrect新创建一个结构老的还在就会同时存在多条波导结果自然不对。第二run命令是Lumerical脚本里的核心命令之一它执行当前FDTD仿真区域内的所有仿真任务。这里要注意run会跑完整个仿真如果仿真时间太长扫描会非常耗时。建议先用粗网格、小仿真区域做验证确保流程跑通了再上细网格。第三getdata函数用于从监视器提取数据。具体能提取什么数据取决于监视器类型。比如模式监视器可以提取neff功率监视器可以提取T透过率。提取前最好先查一下对应监视器支持的数据类型否则会报错。4.3 扫描结果的可视化与保存跑完扫描之后结果如果只是存在终端里等于白跑。Lumerical脚本支持把结果保存成多种格式最常用的是.mat格式方便后续用MATLAB或者Python做后处理。上面的代码里write(neff_vs_width.mat, results)就是把整个二维数组保存为.mat文件。这个文件在Lumerical的安装目录或者你的工作目录里可以找到。如果你想在Lumerical里直接画图也可以用plot命令plot(results(:, 1), results(:, 2), Width (m), Effective Index, neff vs width);这样会弹出一个绘图窗口方便快速检查趋势是否正确。注意plot接受的数据是二维数组第一个参数是x轴第二个参数是y轴后面的字符串是坐标轴标签和标题。5. 常见脚本报错与排查实录5.1 坐标体系不统一结构长在奇怪的地方最典型的报错场景你建了一个结构结果渲染出来看到它跑到了老远的地方或者干脆看不到。这种问题十有八九出在坐标设定上。Lumerical的默认单位是米你写500表示500米。如果用户从GUI导入一个尺寸为500nm的结构脚本里重建时候写成set(x span, 500)那这个结构就是500米宽——直接脱离仿真区域了。解决办法有两种一是自己做好单位换算纳米就写500e-9二是用Lumerical提供的set(x span, 500e-9)明确带上单位。我个人建议全程使用SI单位虽然写起来长一点但排查问题的时候一眼就能看出来哪里不对劲。另外一个小细节结构中心坐标x和跨距x span是两回事。很多人会把这两个概念搞混导致结构永远偏在一边。x决定结构的中心位置x span决定结构在x方向上的总长度。5.2 变量被全局覆盖循环结果全错有时候你会遇到这种情况扫描完一堆宽度最后发现所有仿真结果一模一样。这种问题通常是因为循环里没有正确修改波导的属性或者说改完又被后面覆盖了。我之前踩过的坑是这样的为了省时间我在仿真的模型里既建立了si_waveguide又建立了si_waveguide_copy结果循环里select(si_waveguide)和select(si_waveguide_copy)混着用数据改来改去最后测到的全是同一根波导的模场分布。排查这类问题有个技巧在循环里插入一行?si_waveguide y span num2str(getnamed(si_waveguide, y span));打印一下当前波导的实际宽度。如果打印出来的宽度根本没变那就是select或者set出了问题缩小排查范围会非常快。5.3 快速调试脚本的三个实用小技巧调试脚本最原始的方法是打印输出。Lumerical脚本里?就是打印命令的简写比如?x_center; ?getnamed(si_waveguide, y span);这样可以在命令行看到变量当前的值快速定位问题。第二个技巧是用deleteall清理场景。脚本跑到一半如果结构乱了手动删除非常麻烦。可以在脚本开头加一个deleteall;每次运行前自动清空所有结构。这样做的好处是脚本可以被反复执行而不会产生垃圾结构。第三个技巧是分段执行。Lumerical的脚本编辑器支持选中部分代码执行你不需要每次都从头跑到尾。比如你只想调试结构创建这部分就选中这几行右键执行。这个习惯能让你快速迭代不用每次等完整的仿真跑完。5.4 性能调优结构复杂了脚本卡成PPT最后一个问题是性能层面的。当你的结构非常复杂比如上千个光栅齿每个齿都是一个独立的矩形实体FDTD在建模阶段就要处理大量的几何对象可能导致仿真区域网格化变得异常缓慢。解决方案之一是使用Lumerical的“结构组”功能把一个周期单元定义为一个组然后用阵列方式复制。脚本里对应的是addstructuregroup命令。这种方式不仅能减少几何对象的数量还能让结构更加清晰有序。另外一个方案是调整网格精度。FDTD仿真的本质是时间域差分结构的几何精度和网格精度直接相关。如果你发现网格化时间过长可以考虑使用“保形网格”功能或者适当放宽网格精度。但要注意网格变粗会导致数值误差变大需要在精度和速度之间做一个平衡。我的经验是先粗跑探路确定结果大致没问题再上细网格出正式数据。最后分享一点个人心得从GUI建模切换到脚本建模最初阶段确实会有点不习惯明明点几下就能出的结构非要用代码写半天。但一旦你完成第一个参数化扫描这个观念就会瞬间扭转过来——脚本建模的价值不在“建”这一个动作上而在“改”的灵活性和“批量跑”的效率上。每次写脚本的时候我都会给自己一条准则如果某个结构要改两次以上就值得写成参数化如果要改五次以上就值得写成可扫面的循环。这十年的仿真生涯里这条准则帮我省下的时间不计其数。最后分享一个小技巧我习惯给每个结构的名字加上前缀比如s_代表结构structurem_代表监视器monitors_代表光源sourcef_代表仿真区域FDTD。虽然这是个人习惯但在脚本越来越长的时候这个简单的命名规范能让你在几十行代码里快速定位想要操作的对象。你可以试试说不定会喜欢这种组织方式。