comsol磁可调双带吸收器。
搞电磁超材料的兄弟应该都懂，玩双带吸收器最头疼的就是怎么动态调谐。传统结构一旦加工成型，吸收峰就焊死在固定频段了。最近在COMSOL里折腾磁可调方案发现个骚操作——在铁氧体基底里埋钇铁石榴石（YIG）阵列，外挂静磁场就能实现双峰平移，实测效果贼带劲。
先看核心结构：基底用环氧树脂掺铁氧体粉末，中间嵌YIG圆柱阵列。注意这里要搞两套参数化几何，用COMSOL的Java API写个自动生成阵列的脚本：
for (int i=0; i<numRows; i++) { for (int j=0; j<numCols; j++) { model.geom("geom1").feature().create("cyl"+i+j, "Cylinder"); model.geom("geom1").feature("cyl"+i+j).set("r", 0.8e-3); model.geom("geom1").feature("cyl"+i+j).set("pos", new double[]{i*2e-3, j*2e-3, 0}); model.geom("geom1").feature("cyl"+i+j).set("h", 1.5e-3); } }这段循环生成YIG柱的关键在于设置间距参数时留出1.2mm的间隙，防止磁耦合过强导致频带粘连。柱高1.5mm刚好让TE波在垂直方向形成驻波共振。
comsol磁可调双带吸收器。
材料参数设置更讲究，YIG的磁导率张量必须绑定外部磁场变量：
mu_r = [mu_xx mu_xy 0; -mu_xy mu_yy 0; 0 0 mu_zz]; mu_xx = @(H0) 1 + (deltaM*gamma*H0)/(gamma^2*H0^2 - omega^2); mu_xy = @(H0) (deltaM*omega)/(gamma^2*H0^2 - omega^2);这里的gamma是旋磁比，deltaM对应饱和磁化强度。重点是通过H0参数绑定全局变量，后续扫参时直接在COMSOL里拖滑块调磁场强度，不用重新编译材料属性。
仿真时发现个坑：当磁场超过800A/m时低频吸收峰会突然劈裂成两个子峰。后来发现是YIG柱的尺寸接近磁畴壁移动临界值，得在频域求解器里开启非线性磁化选项：
study.step("frequency").set("plist", "linspace(2e9,20e9,180)"); study.step("frequency").set("punit", "Hz"); study.step("frequency").set("preconditioner", "auto"); study.step("frequency").set("nonlinear", "on"); # 关键开关！打开非线性选项后计算时间翻倍，但能准确捕捉磁饱和效应。实测双峰调谐范围从6-8GHz扩展到5.5-12GHz，Q值还能保持在35以上。
最后给个实用技巧：用参数化扫描同时扫磁场和频率，然后用结果表的矩阵操作提取最佳参数组合：
ResultTable rt = model.result().table("data"); double[][] S11 = rt.getRealPrimitive(); double maxAbsorp = 0; for (int i=0; i<S11.length; i++) { for (int j=0; j<S11[0].length; j++) { if (S11[i][j] < -20 && (i*0.1e9 + j*0.05e9) > maxAbsorp) { maxAbsorp = ... // 找出双带都满足-20dB的参数点 } } }这种磁电协同设计比传统机械调谐快三个量级，下次做智能隐身衣可以考虑这方案。不过注意别把磁场加太大，YIG柱崩了得重画几何，别问我是怎么知道的...
