聊到“波导端口模式的色散特性曲线”很多刚把CST Studio Suite用得比较熟的朋友都会在某个瞬间卡住。端口模式大家肯定都见过——建模时加一个Waveguide PortCST会自动列出一串模式、给出截止频率甚至能看场图但真想把某个模式的传播常数β随频率变化的完整曲线画出来界面上却找不到一个写着“Dispersion Curve”的现成按钮。这篇文章就是来填这个坑的。我会把两种常用做法完整走一遍一种是用频域求解器配合端口模式直接导出传播常数适合工程快速评估另一种是用本征模求解器加相位扫描手动构建色散图适合研究、论文或周期性结构需要完整ω-β关系的场景。两种方法各有适用边界我会把原理、步骤、坑都讲透。1. 色散曲线到底是什么你真正需要哪一种1.1 端口模式、传播常数与 ω-β 图三者的关系很多新手会把“端口模式”和“色散曲线”当成两个独立的东西其实它们是同一件事的两面。波导端口模式是CST在端口截面上求解二维本征值问题得到的场分布。每个模式对应一组横向场结构比如矩形波导的TE10、TE20、TM11等。这个本征值问题得到的是每一个模式的截止波数 kc它只和截面尺寸、填充介质有关和纵向传播无关。而色散特性曲线说的是这个模式在波导里“能不能传、传得多快”的频率关系。具体而言无耗波导中沿纵向传播常数 β 满足β sqrt(k0² - kc²)其中 k0 ω√(με) 是自由空间波数。当 k0 kc 时 β 是实数模式可以传播当 k0 kc 时 β 变成纯虚数模式截止只能以衰减形式存在。把每个模式在频率 f 下的 β 画出来就是 ω-β 图也就是工程上常说的色散曲线。所以端口模式计算结果里包含的截止频率本质上就是色散曲线与频率轴的交点。你要的完整曲线其实就是把端口模式这个“截面本征解”沿纵向“展开”到各个频率点的连续结果。理解了这一层后面所有操作逻辑都顺了。1.2 三种典型需求对应三条不同路线我在实际工程里被问过很多次“怎么画色散曲线”但细问之下大家的需求其实分三种第一种只想知道某个波导端口支持哪些模式、截止频率多少、场分布长什么样。这种最简单加好Waveguide Port在端口属性里就能看到模式列表跑一次求解器后场图也有了。第二种想知道某个工作频带内模式传播常数 β 随频率如何变化用于判断带宽、相移长度、群速度等。这种推荐用频域求解器加波导端口后处理中直接提取传播常数或者用一小段波导的S参数反推。第三种面向周期性结构比如慢波结构、光子晶体、频率选择表面需要画出第一布里渊区内的完整色散图。这种要走本征模求解器加相位扫描的路线手动构建 ω-β 曲线。很多人拿着第三种需求去问第一种问题的答案自然找不到按钮。先确认自己的场景再选路线能省大量时间。下面两张表把对应关系列清楚你的真实需求推荐做法适用场景模式列表、截止频率、场图波导端口 时域/频域求解器常规微波器件、馈电结构设计某模式 β 随频率变化曲线频域求解器 端口模式后处理或 S21 反推滤波器、移相器、匹配段设计周期结构第一布里渊区色散图本征模求解器 相位扫描慢波结构、超材料、FSS、论文分析1.3 为什么CST不给你一个“一键出曲线”的按钮CST是一款通用高频电磁仿真软件不是专用模式分析工具。它把“求解端口模式”“求解本征频率”“计算S参数”这些底层能力拆成一个个模块方便在不同场景下组合使用。但这也意味着像“波导模式的完整色散曲线”这种跨了截面本征解和纵向传播关系的中间产物需要用户自己把几个模块串起来。再加上不同人对“色散曲线”的定义不同——有人要的是 β 曲线有人要的是频率-波数关系图有人只要截止频率。软件不可能用一个按钮覆盖所有口径。明白了这一点你就不会在界面上盲目地找“Dispersion”按钮了而是按下面的路线去组合求解器。2. 路线一用端口模式结果直接导出传播常数2.1 频域求解器下的模式计算与结果位置先推荐最直接的一条路用频域求解器Frequency Domain Solver算一段两端带波导端口的均匀波导。频域求解器在每个频点都会重新计算端口模式所以传播常数 γ α jβ 这个量在求解过程中天然存在只是藏得有点深。操作上这样走先建一个波导模型比如标准的矩形波导中间填真空或介质。在两端各加一个Waveguide Port。添加端口时注意端口平面必须垂直于传播方向端口尺寸要比波导截面稍微大一点避免截断带来的误差。在端口属性里可以设置计算多少个模式我一般会根据需要设 3 到 5 个宁可多算两个也不用后补。然后选择频域求解器设置好扫频范围直接运行。求解完成后在导航树的 2D/3D Results 下展开 Port Modes能看到每个模式在不同频率下的电场/磁场分布。而传播常数相关的结果在部分版本中会以 1D 曲线的形式出现在结果树里名字通常带 Propagation Constant、Gamma 或者 Complex Propagation 之类的关键词。如果版本里找不到直接输出还有一个变通办法查看端口的模式阻抗 Z(ω)。对无耗波导模式阻抗和 β 之间满足明确的解析关系比如矩形波导TE模式的波阻抗 Z_TE ωμ/β反解一下就能得到 β。这个方法我从老版本开始就用过算是兼容性最稳的路径。2.2 找不到传播常数输出时的替代做法从S21反推CST版本差异很大有时候结果树里确实翻不到现成的传播常数曲线。这时候我建议别纠结直接用均匀波导段的S参数反推方法可靠且物理意义清晰。取一段长度 L 的均匀波导两端加匹配良好的波导端口。只考虑单模传输频段内的情况这一段可以等效为理想传输线S21 与传播常数之间满足γ -ln(S21) / L其中 S21 是复数。得到 γ 后实部就是衰减常数 α虚部就是相位常数 β。在无耗模型里 S21 幅度接近 1α 几乎为 0主要关心虚部。具体做法是在CST里跑完频域求解后把 S21 的幅度和相位导出成 ASCII 文件然后在外部脚本里做计算。要注意的是导出的相位通常是 -180° 到 180° 的包裹值必须先做 unwrap 处理才能得到随频率连续增加的相位曲线。这一点是新手最容易踩的坑——不 unwrap 的话β 曲线会在每个 2π 处跳变画出来全是锯齿。用一段 WR-90 尺寸的矩形波导举例宽边 22.86mm、窄边 10.16mm取 L 50mm空气填充。在 8 到 12 GHz 频段内只有 TE10 单模传播S21 的相位会近似线性增长。用 Python 处理时大致是这样import numpy as np import matplotlib.pyplot as plt # 假设从CST导出了三列数据freq(GHz), |S21|, angle(S21)(deg) data np.loadtxt(s21_export.txt, delimiter,) f_ghz data[:, 0] phase_deg data[:, 2] L 0.05 # 波导段长度单位米 phase_unwrap np.unwrap(np.deg2rad(phase_deg)) beta -phase_unwrap / L # 注意S21相位随长度衰减加负号 plt.plot(f_ghz, beta) plt.xlabel(Frequency (GHz)) plt.ylabel(beta (rad/m)) plt.title(TE10 Propagation Constant from S21) plt.grid(True) plt.show()这个方法的限制在于必须确保端口匹配良好且处于单模工作区。如果工作频带内有多个模式传播S21 就是多模叠加的结果不能简单用一个 γ 描述。这时候还是回去用端口模式的后处理结果更可靠。2.3 这种方法最适合什么局限在哪路线一最大的好处是快而且完全基于实际模型损耗、色散介质、非标准截面都能自动考虑进去。比如介质填充波导、脊波导这类没有解析公式的结构用这个方法能直接拿到真实的结构色散关系这是它作为工程方法的核心价值。局限也同样明显。首先如果模型里有一段不连续结构比如弯头、阶梯、滤波器腔体那 S21 反推法就不适用了因为此时测到的相位变化不只是均匀段传播造成的。其次多模区域时端口模式的混合会让传播常数提取变得混乱不同模式之间会耦合。最后S21 反推方法只能给出可传播模式β 为实数的频段截止区以下的衰减特性不容易准确提取。所以我的使用习惯是设计均匀传输线或移相段时用这条路线快速出曲线做复杂器件或要研究截止区特性时切换到本征模相位扫描法。3. 路线二本征模求解器 相位扫描构建色散图3.1 原理把“传播”换成“边界条件”很多做模式分析的人一听到“色散图”首先想到的就是光子晶体那套周期结构仿真方法。其实均匀波导也可以借用这个思路只是概念上要转一个弯。设想一段沿 z 方向均匀的波导任意模式沿 z 的场变化是 e^{-jβz}。如果我截取一个长度为 L 的元胞并要求元胞两端场满足周期边界条件即一端到另一端相位变化为 φ那么自然有φ βL反过来如果我在CST里把 z 方向设置成周期性边界条件并指定两个周期面之间的相位差 φ然后求解本征频率那我得到的就是“在给定 β φ/L 时该模式能够存在的频率”。把φ从 0 扫到 360°就得到一组(f, β)数据点连起来就是色散曲线。这个方法的关键在于本征模求解器本质上是在求解“满足边界条件的谐振频率”而周期性边界条件把纵向传播常数转化成了边界上的相位约束。横截面上的本征解加上纵向的相位条件正好完整描述了波导模式。这个方法对均匀波导、周期加载波导、慢波结构都适用是通用性最强的一条路。3.2 建模与边界设置的关键细节具体操作从建模开始。以矩形波导为例画一个长方体尺寸为宽 a、高 b、长度 L材料设为真空。这里 L 的选择有讲究L 太小会让网格数量激增因为纵向也要剖分L 太大则会让 β 轴范围变窄即第一布里渊区边界 π/L 变小曲线容易被折叠。我一般取工作波长的 1/20 到 1/30比如 10 GHz 时取 L 0.5 到 1 mm效果就不错。边界条件设置是这条路线最核心的环节。横向 x 和 y 方向对于金属波导直接设为 Electric 边界切向电场为零即可等效于理想金属壁。如果结构有对称性比如矩形波导可以沿对称面设置 Magnetic 边界来减半求解区域但要注意有些模式的对称性与边界条件不匹配会漏掉模式。新手我不建议用对称面先跑通无对称版本再说。纵向 z 方向设置为 Periodic 边界CST 会要求填写 Phase Shift。这个值不要写死而是填一个参数比如 phi_deg。这样后面就能通过参数扫描批量改变相位差实现自动扫点。相位差的范围一般从 0 到 360°步长 5° 到 10°。对于均匀波导如果把 L 取得够小0 到 360° 对应的 β 范围足够覆盖工程需要的波数区域。3.3 参数扫描与批量求解建模和边界设好后选择 Eigenmode Solver。在求解器设置里需要指定频率范围这个范围应当覆盖你关心的所有模式的工作频段比如上一节那个例子可以设 5 到 15 GHz。模式数量也要设够我一般设为 3 到 6 个根据你想看几条分支来定。然后就是在参数扫描里做文章。把 phi_deg 设为扫描参数从 0 扫到 360步长 5。每次求解器会在该相位差下算出若干个本征频率。注意Eigenmode Solver 的频率范围会滤掉低于下限和高于上限的模式所以范围别设得太窄否则高阶模式会被切掉。整个扫描跑完后你会得到大量数据每组相位差对应一组模式频率。把这些数据按模式编号整理就能得到多条随 β 变化的频率分支。这里有一个比较麻烦的问题模式的编号顺序在不同相位点之间可能会变化即某个模式的频率曲线在跨越与其他模式简并或交叉点时编号会发生跳变。如果直接用编号画图会发现曲线不连续、出现“断枝”。解决办法是用频率连续性做分组即从上一相位点的模式频率出发找下一相位点中频率最接近的值作为同一模式的延续。这个分组逻辑我后面会展开讲。3.4 数据提取从本征频率还原 β扫描求解完成后数据可以从结果树中导出。在 Navigation Tree 里找到 Eigenmode 相关结果右键选择 Export通常可以导出每个模式在不同参数下的本征频率。导出的格式是文本文件包含参数名、参数值和频率值。拿到数据后β 的计算很简单β (π/180) * phi_deg / L注意这里 phi_deg 要先从度转弧度L 的单位是米得到 β 的单位是 rad/m。频率单位如果是 Hz画图时往往除以 1e9 转成 GHz 更直观。我习惯把数据整理成两列β 和 f然后直接丢进绘图脚本。示例代码如下import numpy as np import matplotlib.pyplot as plt # 假设导出的文件每行是: phi_deg, mode1_freq_Hz, mode2_freq_Hz, ... data np.loadtxt(eigenmode_result.txt) phi_deg data[:, 0] L 0.001 # 元胞长度单位米这里按1mm示例 beta np.deg2rad(phi_deg) / L modes data[:, 1:] # 每列对应一个模式 plt.figure(figsize(6, 4.5)) for i in range(modes.shape[1]): plt.plot(beta, modes[:, i] / 1e9, o-, ms3, labelfMode {i1}) plt.xlabel(beta (rad/m)) plt.ylabel(Frequency (GHz)) plt.legend() plt.grid(True) plt.show()从这个图上可以直接读出几个关键信息曲线和频率轴的交点就是截止频率曲线在截止附近的斜率趋近于 0对应群速度趋近于 0高频端曲线逐渐逼近自由空间色散直线。只要数据质量没问题这些特征都会清晰可见。4. 验证与实用化让曲线可信、有用4.1 用解析公式做交叉验证无论你走哪条路线拿到色散曲线后我强烈建议做一步交叉验证——用理论公式检查曲线长对不对。这一步骤能帮你快速发现边界条件设置错误、模式编号错乱、网格太粗等问题比盲目信任仿真结果可靠得多。矩形波导是最容易验证的结构。对空气填充的矩形波导TE_mn 模式的截止频率为f_c (c/2) * sqrt((m/a)² (n/b)²)其中 c 是光速a 是宽边b 是窄边。传播区内的 β 为β sqrt((2πf/c)² - (mπ/a)² - (nπ/b)²)以 WR-90 矩形波导为例a 22.86mmb 10.16mm前几个模式的截止频率如下模式截止频率GHzTE106.56TE2013.12TE0114.76TE11 / TM1116.48如果你用本征模相位扫描法画的曲线里TE10 分支的截止点明显偏离 6.56 GHz那一定哪里出了错最可能是横向边界条件没设对或者 L 取值太大导致布里渊区折叠干扰判断。对于圆波导TE11 的截止频率是 f_c 1.8412 * c / (2πr)。比如半径 r 10mm 时TE11 截止约 8.79 GHz。这种解析值就像照妖镜一比对就能发现问题。4.2 从色散曲线提取相速度、群速度和截止频率色散曲线不只是论文里的配图它在工程中是实打实可以用的。最常见的是提取相速度和群速度判断带宽和信号失真特性。相速度定义为v_p ω / β群速度定义为v_g dω / dβ在无耗波导中TE 模式的相速度大于光速群速度小于光速这是正常的不违反相对论因为相速度不传递能量。实际工作中我经常用群速度曲线来判断一个传输段的色散严重程度。如果在一定带宽内 v_g 变化剧烈信号经过后会明显展宽这在脉冲信号传输、高精度移相器设计中都是要尽量避免的。从数值结果中提取 v_g 时不需要解析微分直接用相邻数据点做差分即可v_g ≈ 2π * (f[i1] - f[i]) / (β[i1] - β[i])在截止频率附近β 快速变化而 f 变化缓慢差分值趋近于 0对应群速度趋近于 0。这一特征在慢波结构中会被刻意利用——设计周期结构时让工作点落在色散曲线变平的区域就能实现低群速度甚至“慢光”效应。4.3 数据导出与画图脚本CST 本身有后处理模板可以把 β、f 等量直接画成图但说实话真正做设计时我还是喜欢导出数据用外部脚本画。原因是在外部脚本里可以灵活地对模式分组、做差分、拟合、对比理论曲线CST 的后处理界面适合查看不适合反复调试。导出数据时注意单位问题。CST 内部频率按 Hz 保存传播常数按 1/m 保存。导出文本后我通常会加一个简单的文件头备注列名避免过两周再回来时自己都忘了哪列是哪列。这种事我在项目里发生过不止一次提醒各位养成好习惯。如果要在论文里用画图时建议用散点加连线的方式既能体现数据点又容易看趋势。线宽、字体大小按期刊要求调整即可。5. 常见问题与避坑指南5.1 模式窜扰导致曲线“断枝”本征模相位扫描法跑完后最常见的异常就是曲线存在断点或者某条分支在中途突然跳到另一条分支上。发生这种情况通常不是求解器算错了而是模式编号在不同相位点之间发生了变化。当一个模式的色散曲线接近另一个模式时求解器输出的排序可能交换。处理办法有两个层面。第一层在后处理分组时用“延续性”逻辑即下一相位点的第 j 个模式选择与上一相位点某个模式频率差值最小的那个进行匹配。这个逻辑写几十行 Python 就能实现。第二层干脆在画图前人工检查场图在可疑频点查看模式的电场/磁场分布确认分支对应的是同一个物理模式。两种方式结合可以解决绝大多数断枝问题。5.2 简并模式的处理与场图判认方形截面的波导、圆波导、对称性较高的结构里会出现简并模式比如正方形波导的 TE10 与 TE01 完全简并圆波导的 TE21 与 TM01 也可能在某段频率范围内靠得很近。简并模式的色散曲线完全重合画出来就是一条线叠着一条线。这时候只看频率曲线无法区分模式必须借助场图。CST 的端口模式结果里能看每个模式的电场和磁场矢量分布。矩形波导中TE10 的电场只有一个横向分量且沿宽边呈半正弦分布TE01 的电场沿窄边变化。看一眼场图就能确认身份。工程上如果简并模式同时被激励往往会导致极化混乱或能量泄漏所以提前识别并抑制是很重要的。5.3 损耗介质的 α 曲线与单位换算如果用有损介质填充波导比如常见的高频板材、陶瓷填充色散曲线不再只是 β 实部一条线衰减常数 α 也会随频率变化。此时传播常数 γ α jβ 的实部 α 表示能量衰减。CST 的频域求解器端口模式结果中如果给出的是 γ 的复数形式实部单位一般是 Np/m奈培每米。工程上更常用的是 dB/m换算关系是乘以 8.686。另外α 曲线也可以由材料的损耗角正切粗略估算但精确值还是以全波仿真为准。看 α 曲线时要特别注意截止频率附近的特性——很多模式在接近截止时 α 会急剧上升这就是通常说的“截止衰减”。5.4 网格、频率范围与计算量的平衡相位扫描法要跑很多个频率点计算量不小网格一密就是十几甚至几十个求解任务排队。我的策略是分两步走第一步用默认网格或稍粗的网格跑通全流程确认曲线形态、模式分支数量都合理第二步再对感兴趣的频段和相位区间局部加密网格细化曲线。不要一上来就追求最高精度否则一次扫描可能跑几个小时才发现边界条件设错了。频率范围也要精打细算。范围设得太大Eigenmode Solver 会在每个相位点求解很多不需要的高阶模式白白增加计算时间范围设得太小又会漏掉需要关注的模式分支。建议先根据解析公式估算目标模式的截止频率范围然后把频率范围设为从略低于最低截止频率到比最高关心频率高 20% 左右。网格方面本征模求解器对网格质量比 S 参数求解更敏感尤其是场变化剧烈的区域。矩形波导这种简单结构用默认六面体网格就够了波纹波导、脊波导这类结构要关注金属边缘的网格加密。跑完一组扫描后用同一相位点对比两套网格密度的结果如果频率变化小于 0.1%基本可以认为收敛了。我个人的体会是色散曲线这件事本质上是在“模式的本征值”和“纵向传播条件”之间搭桥。CST 里没有现成按钮是因为解决方案本来就该由用户组合模块来完成。多跑几次之后你会慢慢形成自己的模板和脚本再做类似任务时效率会非常高。最后分享一个小技巧把建模、边界条件、参数扫描设置、导出脚本这些步骤固化成一个可复用的工程模板下次遇到不同波导结构时只要改尺寸和频率范围就能直接出一张完整的色散图能省下大量重复操作的时间。 
长河君 · 2026/10/19 1:47:05 · 阅读 1.2k CST波导端口模式色散特性曲线绘制：两种实用方法详解
聊到“波导端口模式的色散特性曲线”很多刚把CST Studio Suite用得比较熟的朋友都会在某个瞬间卡住。端口模式大家肯定都见过——建模时加一个Waveguide PortCST会自动列出一串模式、给出截止频率甚至能看场图但真想把某个模式的传播常数β随频率变化的完整曲线画出来界面上却找不到一个写着“Dispersion Curve”的现成按钮。这篇文章就是来填这个坑的。我会把两种常用做法完整走一遍一种是用频域求解器配合端口模式直接导出传播常数适合工程快速评估另一种是用本征模求解器加相位扫描手动构建色散图适合研究、论文或周期性结构需要完整ω-β关系的场景。两种方法各有适用边界我会把原理、步骤、坑都讲透。1. 色散曲线到底是什么你真正需要哪一种1.1 端口模式、传播常数与 ω-β 图三者的关系很多新手会把“端口模式”和“色散曲线”当成两个独立的东西其实它们是同一件事的两面。波导端口模式是CST在端口截面上求解二维本征值问题得到的场分布。每个模式对应一组横向场结构比如矩形波导的TE10、TE20、TM11等。这个本征值问题得到的是每一个模式的截止波数 kc它只和截面尺寸、填充介质有关和纵向传播无关。而色散特性曲线说的是这个模式在波导里“能不能传、传得多快”的频率关系。具体而言无耗波导中沿纵向传播常数 β 满足β sqrt(k0² - kc²)其中 k0 ω√(με) 是自由空间波数。当 k0 kc 时 β 是实数模式可以传播当 k0 kc 时 β 变成纯虚数模式截止只能以衰减形式存在。把每个模式在频率 f 下的 β 画出来就是 ω-β 图也就是工程上常说的色散曲线。所以端口模式计算结果里包含的截止频率本质上就是色散曲线与频率轴的交点。你要的完整曲线其实就是把端口模式这个“截面本征解”沿纵向“展开”到各个频率点的连续结果。理解了这一层后面所有操作逻辑都顺了。1.2 三种典型需求对应三条不同路线我在实际工程里被问过很多次“怎么画色散曲线”但细问之下大家的需求其实分三种第一种只想知道某个波导端口支持哪些模式、截止频率多少、场分布长什么样。这种最简单加好Waveguide Port在端口属性里就能看到模式列表跑一次求解器后场图也有了。第二种想知道某个工作频带内模式传播常数 β 随频率如何变化用于判断带宽、相移长度、群速度等。这种推荐用频域求解器加波导端口后处理中直接提取传播常数或者用一小段波导的S参数反推。第三种面向周期性结构比如慢波结构、光子晶体、频率选择表面需要画出第一布里渊区内的完整色散图。这种要走本征模求解器加相位扫描的路线手动构建 ω-β 曲线。很多人拿着第三种需求去问第一种问题的答案自然找不到按钮。先确认自己的场景再选路线能省大量时间。下面两张表把对应关系列清楚你的真实需求推荐做法适用场景模式列表、截止频率、场图波导端口 时域/频域求解器常规微波器件、馈电结构设计某模式 β 随频率变化曲线频域求解器 端口模式后处理或 S21 反推滤波器、移相器、匹配段设计周期结构第一布里渊区色散图本征模求解器 相位扫描慢波结构、超材料、FSS、论文分析1.3 为什么CST不给你一个“一键出曲线”的按钮CST是一款通用高频电磁仿真软件不是专用模式分析工具。它把“求解端口模式”“求解本征频率”“计算S参数”这些底层能力拆成一个个模块方便在不同场景下组合使用。但这也意味着像“波导模式的完整色散曲线”这种跨了截面本征解和纵向传播关系的中间产物需要用户自己把几个模块串起来。再加上不同人对“色散曲线”的定义不同——有人要的是 β 曲线有人要的是频率-波数关系图有人只要截止频率。软件不可能用一个按钮覆盖所有口径。明白了这一点你就不会在界面上盲目地找“Dispersion”按钮了而是按下面的路线去组合求解器。2. 路线一用端口模式结果直接导出传播常数2.1 频域求解器下的模式计算与结果位置先推荐最直接的一条路用频域求解器Frequency Domain Solver算一段两端带波导端口的均匀波导。频域求解器在每个频点都会重新计算端口模式所以传播常数 γ α jβ 这个量在求解过程中天然存在只是藏得有点深。操作上这样走先建一个波导模型比如标准的矩形波导中间填真空或介质。在两端各加一个Waveguide Port。添加端口时注意端口平面必须垂直于传播方向端口尺寸要比波导截面稍微大一点避免截断带来的误差。在端口属性里可以设置计算多少个模式我一般会根据需要设 3 到 5 个宁可多算两个也不用后补。然后选择频域求解器设置好扫频范围直接运行。求解完成后在导航树的 2D/3D Results 下展开 Port Modes能看到每个模式在不同频率下的电场/磁场分布。而传播常数相关的结果在部分版本中会以 1D 曲线的形式出现在结果树里名字通常带 Propagation Constant、Gamma 或者 Complex Propagation 之类的关键词。如果版本里找不到直接输出还有一个变通办法查看端口的模式阻抗 Z(ω)。对无耗波导模式阻抗和 β 之间满足明确的解析关系比如矩形波导TE模式的波阻抗 Z_TE ωμ/β反解一下就能得到 β。这个方法我从老版本开始就用过算是兼容性最稳的路径。2.2 找不到传播常数输出时的替代做法从S21反推CST版本差异很大有时候结果树里确实翻不到现成的传播常数曲线。这时候我建议别纠结直接用均匀波导段的S参数反推方法可靠且物理意义清晰。取一段长度 L 的均匀波导两端加匹配良好的波导端口。只考虑单模传输频段内的情况这一段可以等效为理想传输线S21 与传播常数之间满足γ -ln(S21) / L其中 S21 是复数。得到 γ 后实部就是衰减常数 α虚部就是相位常数 β。在无耗模型里 S21 幅度接近 1α 几乎为 0主要关心虚部。具体做法是在CST里跑完频域求解后把 S21 的幅度和相位导出成 ASCII 文件然后在外部脚本里做计算。要注意的是导出的相位通常是 -180° 到 180° 的包裹值必须先做 unwrap 处理才能得到随频率连续增加的相位曲线。这一点是新手最容易踩的坑——不 unwrap 的话β 曲线会在每个 2π 处跳变画出来全是锯齿。用一段 WR-90 尺寸的矩形波导举例宽边 22.86mm、窄边 10.16mm取 L 50mm空气填充。在 8 到 12 GHz 频段内只有 TE10 单模传播S21 的相位会近似线性增长。用 Python 处理时大致是这样import numpy as np import matplotlib.pyplot as plt # 假设从CST导出了三列数据freq(GHz), |S21|, angle(S21)(deg) data np.loadtxt(s21_export.txt, delimiter,) f_ghz data[:, 0] phase_deg data[:, 2] L 0.05 # 波导段长度单位米 phase_unwrap np.unwrap(np.deg2rad(phase_deg)) beta -phase_unwrap / L # 注意S21相位随长度衰减加负号 plt.plot(f_ghz, beta) plt.xlabel(Frequency (GHz)) plt.ylabel(beta (rad/m)) plt.title(TE10 Propagation Constant from S21) plt.grid(True) plt.show()这个方法的限制在于必须确保端口匹配良好且处于单模工作区。如果工作频带内有多个模式传播S21 就是多模叠加的结果不能简单用一个 γ 描述。这时候还是回去用端口模式的后处理结果更可靠。2.3 这种方法最适合什么局限在哪路线一最大的好处是快而且完全基于实际模型损耗、色散介质、非标准截面都能自动考虑进去。比如介质填充波导、脊波导这类没有解析公式的结构用这个方法能直接拿到真实的结构色散关系这是它作为工程方法的核心价值。局限也同样明显。首先如果模型里有一段不连续结构比如弯头、阶梯、滤波器腔体那 S21 反推法就不适用了因为此时测到的相位变化不只是均匀段传播造成的。其次多模区域时端口模式的混合会让传播常数提取变得混乱不同模式之间会耦合。最后S21 反推方法只能给出可传播模式β 为实数的频段截止区以下的衰减特性不容易准确提取。所以我的使用习惯是设计均匀传输线或移相段时用这条路线快速出曲线做复杂器件或要研究截止区特性时切换到本征模相位扫描法。3. 路线二本征模求解器 相位扫描构建色散图3.1 原理把“传播”换成“边界条件”很多做模式分析的人一听到“色散图”首先想到的就是光子晶体那套周期结构仿真方法。其实均匀波导也可以借用这个思路只是概念上要转一个弯。设想一段沿 z 方向均匀的波导任意模式沿 z 的场变化是 e^{-jβz}。如果我截取一个长度为 L 的元胞并要求元胞两端场满足周期边界条件即一端到另一端相位变化为 φ那么自然有φ βL反过来如果我在CST里把 z 方向设置成周期性边界条件并指定两个周期面之间的相位差 φ然后求解本征频率那我得到的就是“在给定 β φ/L 时该模式能够存在的频率”。把φ从 0 扫到 360°就得到一组(f, β)数据点连起来就是色散曲线。这个方法的关键在于本征模求解器本质上是在求解“满足边界条件的谐振频率”而周期性边界条件把纵向传播常数转化成了边界上的相位约束。横截面上的本征解加上纵向的相位条件正好完整描述了波导模式。这个方法对均匀波导、周期加载波导、慢波结构都适用是通用性最强的一条路。3.2 建模与边界设置的关键细节具体操作从建模开始。以矩形波导为例画一个长方体尺寸为宽 a、高 b、长度 L材料设为真空。这里 L 的选择有讲究L 太小会让网格数量激增因为纵向也要剖分L 太大则会让 β 轴范围变窄即第一布里渊区边界 π/L 变小曲线容易被折叠。我一般取工作波长的 1/20 到 1/30比如 10 GHz 时取 L 0.5 到 1 mm效果就不错。边界条件设置是这条路线最核心的环节。横向 x 和 y 方向对于金属波导直接设为 Electric 边界切向电场为零即可等效于理想金属壁。如果结构有对称性比如矩形波导可以沿对称面设置 Magnetic 边界来减半求解区域但要注意有些模式的对称性与边界条件不匹配会漏掉模式。新手我不建议用对称面先跑通无对称版本再说。纵向 z 方向设置为 Periodic 边界CST 会要求填写 Phase Shift。这个值不要写死而是填一个参数比如 phi_deg。这样后面就能通过参数扫描批量改变相位差实现自动扫点。相位差的范围一般从 0 到 360°步长 5° 到 10°。对于均匀波导如果把 L 取得够小0 到 360° 对应的 β 范围足够覆盖工程需要的波数区域。3.3 参数扫描与批量求解建模和边界设好后选择 Eigenmode Solver。在求解器设置里需要指定频率范围这个范围应当覆盖你关心的所有模式的工作频段比如上一节那个例子可以设 5 到 15 GHz。模式数量也要设够我一般设为 3 到 6 个根据你想看几条分支来定。然后就是在参数扫描里做文章。把 phi_deg 设为扫描参数从 0 扫到 360步长 5。每次求解器会在该相位差下算出若干个本征频率。注意Eigenmode Solver 的频率范围会滤掉低于下限和高于上限的模式所以范围别设得太窄否则高阶模式会被切掉。整个扫描跑完后你会得到大量数据每组相位差对应一组模式频率。把这些数据按模式编号整理就能得到多条随 β 变化的频率分支。这里有一个比较麻烦的问题模式的编号顺序在不同相位点之间可能会变化即某个模式的频率曲线在跨越与其他模式简并或交叉点时编号会发生跳变。如果直接用编号画图会发现曲线不连续、出现“断枝”。解决办法是用频率连续性做分组即从上一相位点的模式频率出发找下一相位点中频率最接近的值作为同一模式的延续。这个分组逻辑我后面会展开讲。3.4 数据提取从本征频率还原 β扫描求解完成后数据可以从结果树中导出。在 Navigation Tree 里找到 Eigenmode 相关结果右键选择 Export通常可以导出每个模式在不同参数下的本征频率。导出的格式是文本文件包含参数名、参数值和频率值。拿到数据后β 的计算很简单β (π/180) * phi_deg / L注意这里 phi_deg 要先从度转弧度L 的单位是米得到 β 的单位是 rad/m。频率单位如果是 Hz画图时往往除以 1e9 转成 GHz 更直观。我习惯把数据整理成两列β 和 f然后直接丢进绘图脚本。示例代码如下import numpy as np import matplotlib.pyplot as plt # 假设导出的文件每行是: phi_deg, mode1_freq_Hz, mode2_freq_Hz, ... data np.loadtxt(eigenmode_result.txt) phi_deg data[:, 0] L 0.001 # 元胞长度单位米这里按1mm示例 beta np.deg2rad(phi_deg) / L modes data[:, 1:] # 每列对应一个模式 plt.figure(figsize(6, 4.5)) for i in range(modes.shape[1]): plt.plot(beta, modes[:, i] / 1e9, o-, ms3, labelfMode {i1}) plt.xlabel(beta (rad/m)) plt.ylabel(Frequency (GHz)) plt.legend() plt.grid(True) plt.show()从这个图上可以直接读出几个关键信息曲线和频率轴的交点就是截止频率曲线在截止附近的斜率趋近于 0对应群速度趋近于 0高频端曲线逐渐逼近自由空间色散直线。只要数据质量没问题这些特征都会清晰可见。4. 验证与实用化让曲线可信、有用4.1 用解析公式做交叉验证无论你走哪条路线拿到色散曲线后我强烈建议做一步交叉验证——用理论公式检查曲线长对不对。这一步骤能帮你快速发现边界条件设置错误、模式编号错乱、网格太粗等问题比盲目信任仿真结果可靠得多。矩形波导是最容易验证的结构。对空气填充的矩形波导TE_mn 模式的截止频率为f_c (c/2) * sqrt((m/a)² (n/b)²)其中 c 是光速a 是宽边b 是窄边。传播区内的 β 为β sqrt((2πf/c)² - (mπ/a)² - (nπ/b)²)以 WR-90 矩形波导为例a 22.86mmb 10.16mm前几个模式的截止频率如下模式截止频率GHzTE106.56TE2013.12TE0114.76TE11 / TM1116.48如果你用本征模相位扫描法画的曲线里TE10 分支的截止点明显偏离 6.56 GHz那一定哪里出了错最可能是横向边界条件没设对或者 L 取值太大导致布里渊区折叠干扰判断。对于圆波导TE11 的截止频率是 f_c 1.8412 * c / (2πr)。比如半径 r 10mm 时TE11 截止约 8.79 GHz。这种解析值就像照妖镜一比对就能发现问题。4.2 从色散曲线提取相速度、群速度和截止频率色散曲线不只是论文里的配图它在工程中是实打实可以用的。最常见的是提取相速度和群速度判断带宽和信号失真特性。相速度定义为v_p ω / β群速度定义为v_g dω / dβ在无耗波导中TE 模式的相速度大于光速群速度小于光速这是正常的不违反相对论因为相速度不传递能量。实际工作中我经常用群速度曲线来判断一个传输段的色散严重程度。如果在一定带宽内 v_g 变化剧烈信号经过后会明显展宽这在脉冲信号传输、高精度移相器设计中都是要尽量避免的。从数值结果中提取 v_g 时不需要解析微分直接用相邻数据点做差分即可v_g ≈ 2π * (f[i1] - f[i]) / (β[i1] - β[i])在截止频率附近β 快速变化而 f 变化缓慢差分值趋近于 0对应群速度趋近于 0。这一特征在慢波结构中会被刻意利用——设计周期结构时让工作点落在色散曲线变平的区域就能实现低群速度甚至“慢光”效应。4.3 数据导出与画图脚本CST 本身有后处理模板可以把 β、f 等量直接画成图但说实话真正做设计时我还是喜欢导出数据用外部脚本画。原因是在外部脚本里可以灵活地对模式分组、做差分、拟合、对比理论曲线CST 的后处理界面适合查看不适合反复调试。导出数据时注意单位问题。CST 内部频率按 Hz 保存传播常数按 1/m 保存。导出文本后我通常会加一个简单的文件头备注列名避免过两周再回来时自己都忘了哪列是哪列。这种事我在项目里发生过不止一次提醒各位养成好习惯。如果要在论文里用画图时建议用散点加连线的方式既能体现数据点又容易看趋势。线宽、字体大小按期刊要求调整即可。5. 常见问题与避坑指南5.1 模式窜扰导致曲线“断枝”本征模相位扫描法跑完后最常见的异常就是曲线存在断点或者某条分支在中途突然跳到另一条分支上。发生这种情况通常不是求解器算错了而是模式编号在不同相位点之间发生了变化。当一个模式的色散曲线接近另一个模式时求解器输出的排序可能交换。处理办法有两个层面。第一层在后处理分组时用“延续性”逻辑即下一相位点的第 j 个模式选择与上一相位点某个模式频率差值最小的那个进行匹配。这个逻辑写几十行 Python 就能实现。第二层干脆在画图前人工检查场图在可疑频点查看模式的电场/磁场分布确认分支对应的是同一个物理模式。两种方式结合可以解决绝大多数断枝问题。5.2 简并模式的处理与场图判认方形截面的波导、圆波导、对称性较高的结构里会出现简并模式比如正方形波导的 TE10 与 TE01 完全简并圆波导的 TE21 与 TM01 也可能在某段频率范围内靠得很近。简并模式的色散曲线完全重合画出来就是一条线叠着一条线。这时候只看频率曲线无法区分模式必须借助场图。CST 的端口模式结果里能看每个模式的电场和磁场矢量分布。矩形波导中TE10 的电场只有一个横向分量且沿宽边呈半正弦分布TE01 的电场沿窄边变化。看一眼场图就能确认身份。工程上如果简并模式同时被激励往往会导致极化混乱或能量泄漏所以提前识别并抑制是很重要的。5.3 损耗介质的 α 曲线与单位换算如果用有损介质填充波导比如常见的高频板材、陶瓷填充色散曲线不再只是 β 实部一条线衰减常数 α 也会随频率变化。此时传播常数 γ α jβ 的实部 α 表示能量衰减。CST 的频域求解器端口模式结果中如果给出的是 γ 的复数形式实部单位一般是 Np/m奈培每米。工程上更常用的是 dB/m换算关系是乘以 8.686。另外α 曲线也可以由材料的损耗角正切粗略估算但精确值还是以全波仿真为准。看 α 曲线时要特别注意截止频率附近的特性——很多模式在接近截止时 α 会急剧上升这就是通常说的“截止衰减”。5.4 网格、频率范围与计算量的平衡相位扫描法要跑很多个频率点计算量不小网格一密就是十几甚至几十个求解任务排队。我的策略是分两步走第一步用默认网格或稍粗的网格跑通全流程确认曲线形态、模式分支数量都合理第二步再对感兴趣的频段和相位区间局部加密网格细化曲线。不要一上来就追求最高精度否则一次扫描可能跑几个小时才发现边界条件设错了。频率范围也要精打细算。范围设得太大Eigenmode Solver 会在每个相位点求解很多不需要的高阶模式白白增加计算时间范围设得太小又会漏掉需要关注的模式分支。建议先根据解析公式估算目标模式的截止频率范围然后把频率范围设为从略低于最低截止频率到比最高关心频率高 20% 左右。网格方面本征模求解器对网格质量比 S 参数求解更敏感尤其是场变化剧烈的区域。矩形波导这种简单结构用默认六面体网格就够了波纹波导、脊波导这类结构要关注金属边缘的网格加密。跑完一组扫描后用同一相位点对比两套网格密度的结果如果频率变化小于 0.1%基本可以认为收敛了。我个人的体会是色散曲线这件事本质上是在“模式的本征值”和“纵向传播条件”之间搭桥。CST 里没有现成按钮是因为解决方案本来就该由用户组合模块来完成。多跑几次之后你会慢慢形成自己的模板和脚本再做类似任务时效率会非常高。最后分享一个小技巧把建模、边界条件、参数扫描设置、导出脚本这些步骤固化成一个可复用的工程模板下次遇到不同波导结构时只要改尺寸和频率范围就能直接出一张完整的色散图能省下大量重复操作的时间。