FDTD结构组脚本进阶：从复制粘贴到理解，自定义任意旋转体（含锥体/圆台）

FDTD结构组脚本进阶：从复制粘贴到理解，自定义任意旋转体（含锥体/圆台）

在FDTD仿真中，几何建模的灵活性往往决定了设计能力的上限。当标准库中的基础几何体无法满足需求时，结构组脚本中的addcustom功能便成为突破限制的利器。本文将带您从几何原理出发，彻底掌握旋转体参数化建模的核心逻辑。

1. 旋转体建模的数学基础

旋转体（Surface of Revolution）的生成遵循严格的数学规则：给定一条平面曲线（母线），绕固定轴旋转一周形成的曲面。在FDTD中，这个原理通过三个关键参数实现：

• 旋转轴（first axis）：决定母线围绕哪根坐标轴旋转
• 母线方程（equation 1）：定义在旋转轴与另一坐标轴构成的平面内的曲线
• 生成方式（create 3D object by）：选择"revolution"即启用旋转体模式

以圆锥生成为例，其母线是直线方程y=kx+b。当这条直线绕y轴旋转时，斜率k决定锥角大小，截距b影响锥顶位置。官网脚本中这段代码正是该原理的实现：

?theta = atan((r_bottom - r_top)/z_span) -- 计算锥体半角 ?ht = r_top/tan(theta) -- 虚拟锥顶到截面的距离 ?eqn = num2str(r_top/ht)+"*(x+"+num2str((z_span/2+ht)*1e6)+")" -- 构建直线方程

2. 参数化建模的四大核心要素

2.1 坐标系与旋转轴设定

first axis参数决定了旋转体的空间取向。常见配置包括：

注意：设置旋转轴后，方程变量必须对应平面内的坐标轴。例如选择y为旋转轴时，方程中的变量应为y和z。

2.2 母线方程构建技巧

母线不限于直线，任何可解析表达的曲线都能生成复杂旋转体。以下是常见曲线类型及对应方程：

1. 多项式曲线：y = ax³ + bx² + cx + d
   • 可创建渐变过渡的瓶状结构
2. 三角函数曲线：y = a*sin(b*x) + c
   • 适合生成波纹状表面
3. 分段函数：组合多个曲线段
   • 实现阶梯状或复杂轮廓

-- 示例：创建正弦波纹旋转面 amplitude = 0.2e-6 -- 振幅200nm period = 1e-6 -- 周期1μm ?eqn = num2str(amplitude)+"*sin(2*pi*"+num2str(1/period)+"*x)"

2.3 Span范围的几何意义

Span参数定义了建模的有效区域，相当于数学上的定义域。设置时需要特别注意：

• 必须完全包含母线曲线
• 超出部分会被硬性截断
• 对于旋转体，非旋转方向的span应至少两倍于母线最大半径

常见错误场景：

• 圆锥母线斜率为0.5，x方向延伸1μm，则y最大值应为0.5μm
• 若设置y span=0.8μm（即±0.4μm），顶部会被截平

2.4 材料与空间定位

通过脚本可以动态设置材料属性，实现更灵活的建模流程：

material = "Si (Silicon) - Palik" index = 3.45 -- 自定义折射率 set("material", material) if(get("material")=="<Object defined dielectric>") { set("index", index) -- 自定义介电材料 }

3. 复杂旋转体实战案例

3.1 平顶高斯透镜建模

结合高斯函数与线性函数，可以创建具有平顶特征的透镜结构：

-- 参数定义 waist_radius = 1e-6 -- 1μm束腰半径 flat_radius = 0.7e-6 -- 平顶区域半径 height = 0.5e-6 -- 最大高度500nm -- 分段函数构建 ?eqn = "("+num2str(height)+"*exp(-(x-"+ num2str(flat_radius)+")^2/"+ num2str(2*waist_radius^2)+")) * (x>"+ num2str(flat_radius)+") + "+ num2str(height)+"*(x<="+ num2str(flat_radius)+")"

3.2 多级微结构阵列

通过脚本批量生成参数化结构组，创建复杂功能表面：

-- 多锥体阵列参数 Nx = 5 -- x方向数量 Ny = 5 -- y方向数量 pitch = 2e-6 -- 间距2μm base_radius = 0.5e-6 -- 基底半径500nm aspect_ratio = 1.5 -- 高宽比 deleteall; -- 清除现有结构 for i = 1,Nx do for j = 1,Ny do addcustom; -- 位置设置 x_pos = (i-1-(Nx-1)/2)*pitch y_pos = (j-1-(Ny-1)/2)*pitch set("x",x_pos); set("y",y_pos); -- 动态尺寸 r_top = base_radius*(0.8+0.4*rand()) height = aspect_ratio*2*r_top*(0.9+0.2*rand()) -- 旋转体参数 set("first axis","z"); set("x span",height); set("y span",4*r_top); set("z span",4*r_top); ?theta = atan(r_top/height); ?eqn = num2str(tan(theta))+"*x"; set("equation 1",eqn); set("create 3D object by","revolution"); end end

4. 高级技巧与调试方法

4.1 曲面平滑度控制

旋转体表面质量由两个参数决定：

1. 网格划分密度：通过mesh order控制
2. 旋转分段数：在custom属性中设置

set("mesh order",3); -- 更高阶的网格划分 set("number of segments",64); -- 旋转方向细分段数

注意：增加分段数会显著提升计算资源消耗，需权衡精度与效率。

4.2 复杂交叠结构处理

当多个旋转体相交时，布尔运算规则决定最终形态：

4.3 常见错误排查

1. 方程不生效：

   • 检查变量名是否匹配旋转轴平面坐标
   • 确认方程语法符合Lumerical表达式规范

2. 结构显示不完整：

   • 验证span范围是否足够包含整个曲线
   • 检查旋转轴设置是否正确

3. 表面出现锯齿：

   • 增加number of segments参数值
   • 提升网格划分密度

-- 调试技巧：可视化母线曲线 ?x_values = linspace(-span/2,span/2,100); ?y_values = execute(eqn); plot(x_values,y_values,"母线验证");

掌握这些原理后，您可以将标准脚本转化为真正的创作工具。某次在设计特殊光场调控器件时，正是通过参数化方程构建的非对称旋转体，实现了传统建模工具难以达到的精确场分布控制。