1. 共面波导色散曲线到底在解决什么问题刚接触CST的人看到“共面波导色散曲线”这几个字第一反应往往是——这玩意儿到底拿来干嘛的我当初也一样翻了好几本微波工程的书公式看了一堆真到软件里点按钮的时候还是懵的。所以这篇东西我不打算从麦克斯韦方程组开始推而是直接告诉你色散曲线本质上就是一条“频率对相位常数”的关系曲线它描述的是电磁波在共面波导结构里传播时不同频率的信号对应的传播特性。为什么这件事重要因为共面波导和普通微带线不一样。微带线的场主要集中在介质基片里而共面波导的场分布在中心导带和两侧地平面之间的槽缝区域介质基片只是其中一部分。这就导致它的色散特性对频率非常敏感尤其在毫米波频段高阶模式、表面波、辐射损耗都会掺和进来。如果你在设计一个宽带共面波导电路比如行波放大器或者宽带混频器不先把色散曲线搞清楚后面匹配网络怎么调都是瞎调。具体来说色散曲线能帮你回答三个核心问题。第一相位常数随频率怎么变这直接决定了你的电路在目标频段内的相位一致性。第二特征阻抗是不是稳定的共面波导的特征阻抗对槽缝宽度和中心导带宽度极其敏感色散曲线能告诉你阻抗随频率的漂移量。第三高阶模式什么时候出现一旦某个频率上出现了非期望的模式你的信号完整性就完蛋了。所以这条曲线不是学术装饰它是设计决策的依据。适合谁来参考这篇内容如果你已经装好了CST打开过界面知道怎么画矩形、怎么设材料但从来没跑过本征模求解器那这篇就是写给你的。如果你已经用过本征模求解器但结果总是不对比如曲线跳变、模式识别混乱那这篇也能帮你排查。至于完全没碰过CST的人建议先花半小时熟悉一下建模界面的基本操作否则后面步骤会卡在很基础的地方。提示共面波导的色散曲线和微带线的色散曲线在低频段趋势相似但在高频段差异巨大。不要拿微带线的经验直接套用尤其是槽缝宽度小于介质厚度的时候。2. 建模前的整体思路与方案选型2.1 为什么用本征模求解器而不是时域求解器CST里画色散曲线最常用的求解器有两个选择时域求解器和本征模求解器。很多人第一次做的时候会想时域求解器不是更直观吗跑一个宽带脉冲做FFT不就能得到频域信息了理论上可以但实际操作中问题很多。时域求解器要得到色散曲线你需要设置两个端口一端激励一端接收然后通过S参数反推相位常数。这个过程涉及相位解缠绕频率高了之后相位跳变非常频繁解缠绕算法稍微出点问题曲线就断了。而且时域求解器跑的是宽带信号共面波导在高频段可能有辐射损耗这些损耗会混在S21里让你分不清到底是传播损耗还是辐射损耗。本征模求解器的逻辑完全不同。它直接求解结构的本征频率和本征场分布每一个本征模式对应一个特定的频率和场型。你只需要设置好周期边界条件让求解器在给定的相位差下找本征频率然后扫描相位差就能得到完整的色散曲线。这个过程不需要端口不需要激励纯粹是求解结构的固有属性。对于共面波导这种开放或半开放结构本征模求解器配合周期边界条件是最干净的做法。我个人的经验是本征模求解器跑色散曲线精度和稳定性都远好于时域求解器反推。尤其是当你要扫描多个模式的时候本征模求解器能直接告诉你每个模式的场分布你一眼就能看出哪个是主模哪个是高阶模。时域求解器反推的话模式识别全靠你自己猜。2.2 周期边界条件的设置逻辑本征模求解器跑色散曲线的核心是周期边界条件。CST里叫“Unit Cell”边界你需要在传播方向上设置相位差扫描。具体来说共面波导的传播方向假设是X轴那么X方向的两个面要设置成周期边界相位差从0度扫到180度。为什么是0到180度因为色散曲线的一个完整布里渊区就是0到180度超过180度就进入下一个布里渊区了曲线是重复的。这里有个细节很多人会搞错相位差的单位是度不是弧度。CST的界面里默认可能是度但如果你从别的软件转过来习惯用弧度一定要检查一下。我见过有人设了3.14度然后纳闷为什么曲线只画了一小段。3.14度对应的相位差几乎为零当然只能得到一个点附近的曲线。另外周期边界条件要求结构在传播方向上必须是严格周期的。也就是说你的共面波导模型在X方向上要能无缝拼接。如果你画了一个有限长的线段然后直接加周期边界求解器会认为这个线段无限重复但线段之间的连接处如果有不连续结果就会出错。所以建模的时候共面波导在传播方向上的长度应该取一个周期单元的长度通常远小于波长但又要足够长以包含结构特征。2.3 端口设置的取舍本征模求解器不需要外部端口但你需要定义端口来提取特征阻抗和模式场。CST里在本征模求解器下可以设置“Port”为“Open (add space)”或者“Waveguide”。对于共面波导我通常建议用“Open (add space)”因为共面波导的场是开放的部分场会辐射到空气中。如果你用“Waveguide”端口相当于强制场在端口处满足波导模式这和实际共面波导的场分布有偏差。端口的尺寸也有讲究。端口的宽度要足够大至少覆盖共面波导的槽缝区域和两侧地平面的一部分。如果端口太小场被截断提取出来的特征阻抗会偏大。我一般会把端口宽度设为中心导带宽度的5到10倍高度设为介质厚度的5倍以上。具体数值要看你的结构但原则是端口边界不能切割场分布的主要区域。注意本征模求解器下端口的数量会影响计算量。如果你只关心主模的色散曲线一个端口就够了。但如果你想同时看多个模式可能需要多个端口来区分。不过CST的本征模求解器本身就能算出多个本征模式端口主要是用来后处理提取参数的。3. 核心参数设置与实操细节3.1 材料参数的设定陷阱共面波导的材料参数看起来简单无非是介质基片的介电常数和损耗角正切加上金属的电导率。但这里有几个坑我踩过不止一次。第一个坑是介电常数的频率相关性。很多资料上给的FR4介电常数是4.4但那是1MHz下的值。到了10GHzFR4的介电常数可能降到4.0甚至更低。如果你做的是宽带色散曲线从1GHz扫到40GHz用常数介电常数算出来的曲线在高频段会有明显偏差。CST里可以设置频率相关的材料你可以导入实测的介电常数曲线或者用Debye模型、Lorentz模型来拟合。对于新手如果只是练手用常数介电常数没问题但要知道这个局限性。第二个坑是金属的厚度和电导率。共面波导的金属层通常很薄几个微米量级。在CST里如果你用“Lossy metal”并设置电导率求解器会考虑趋肤效应。但如果你把金属设成“PEC”理想电导体那就完全没有损耗色散曲线的衰减部分就看不到了。对于色散曲线的相位常数部分PEC和Lossy metal的差异不大但如果你要看衰减常数就必须用Lossy metal。我一般建议先用PEC跑一遍看趋势确认模型没问题再换成Lossy metal跑精确结果。第三个坑是介质基片的厚度。共面波导的色散特性对基片厚度很敏感。基片越薄场越集中在表面色散越弱。基片越厚场越容易进入基片和空气色散越强而且容易激发表面波。CST里设置基片厚度的时候要注意单位。默认可能是毫米但如果你从别的软件导入模型单位可能是米或者微米。我见过有人把500微米设成了500毫米结果跑出来一个巨大的结构求解器直接报内存不足。3.2 网格设置的精度控制本征模求解器的网格设置和时域求解器不太一样。时域求解器有自适应网格加密本征模求解器也有但逻辑不同。本征模求解器是基于四面体网格的网格密度直接影响本征频率的计算精度。对于共面波导槽缝区域的网格必须足够密。因为场主要集中在槽缝附近如果槽缝里只有一两个网格场分布根本算不准。我通常会在槽缝区域手动加密网格设置“Local Mesh”或者“Mesh Refinement”。具体来说槽缝宽度如果是W那么槽缝里的网格尺寸应该小于W/5。如果W是10微米网格尺寸就要小于2微米。这个尺寸在CST里是可以设置的但要注意计算量。网格太密内存和时间都会爆炸。另一个关键是传播方向上的网格。周期边界条件下传播方向的网格应该均匀而且周期长度应该是网格尺寸的整数倍。如果周期长度不是网格尺寸的整数倍求解器会强制调整网格可能导致周期边界处的场不连续。我一般会先确定周期长度然后反推网格尺寸确保整除。还有一个经验本征模求解器的精度设置里“Accuracy”不要设得太高。默认的-30dB或者-40dB已经足够了。设到-60dB计算时间可能翻倍但结果差异在小数点后第三位对于工程应用没有意义。我试过-80dB跑了三个小时结果和-40dB的差异在0.1%以内完全不值得。3.3 相位差扫描的设置方法相位差扫描是本征模求解器跑色散曲线的核心操作。CST里在“Eigenmode”求解器设置里有一个“Phase”选项你可以设置起始相位、终止相位和步长。起始相位设0度终止相位设180度这是标准做法。步长设多少我一般设5度或者10度。5度的话从0到180度有37个点曲线足够平滑。10度的话有19个点对于简单的色散曲线也够用。但如果你要看的模式很多或者曲线有复杂的交叉建议用5度甚至更小。步长太大曲线交叉的地方可能被跳过你会误以为两个模式没有耦合。这里有个隐藏的坑CST的相位扫描是离散的每个相位点独立求解本征模。这意味着求解器在每个相位点都会重新计算本征频率和场分布。如果两个模式的频率非常接近求解器可能在不同的相位点把模式顺序搞混。比如在0度的时候模式1是主模到了90度的时候模式2的频率降到了模式1以下求解器可能把模式2标记为模式1。这会导致色散曲线出现跳变。解决这个问题的办法是在扫描之前先跑一个单点的本征模分析看看在0度和180度分别有哪些模式频率大概是多少。然后根据这个信息在扫描设置里指定你要跟踪的模式数量。CST里可以设置“Number of Modes”你设成3或者4求解器会算出前几个模式。后处理的时候你可以根据场分布来识别哪个是你要的主模。提示如果你发现色散曲线在某个频率附近突然跳变先检查是不是模式识别错了。把那个相位点的场分布调出来看看对比一下相邻相位点的场分布如果场型完全变了那就是模式顺序乱了。4. 完整实操流程与关键步骤4.1 从零开始建立共面波导模型打开CST新建一个“Microwave RF”项目选择“Eigenmode”求解器。单位设成毫米和GHz这是微波频段的常规设置。第一步画介质基片。用“Brick”工具画一个长方体尺寸根据你的设计来。假设基片厚度是0.5毫米长度是传播方向上的一个周期比如2毫米宽度要足够容纳共面波导和两侧地平面比如5毫米。材料选“FR4”或者“Rogers RO4003C”看你的实际设计。第二步画中心导带。在基片上方画一个薄长方体厚度设成0.035毫米这是常见的1盎司铜厚宽度设成0.5毫米长度和基片一样。材料选“Copper”或者“PEC”。第三步画两侧地平面。在中心导带两侧各画一个长方体厚度和中心导带一样宽度比如2毫米长度和基片一样。中心导带和地平面之间的间隙就是槽缝假设槽缝宽度是0.3毫米。那么地平面的内边缘距离中心导带边缘0.3毫米。第四步设置边界条件。在“Boundaries”里X方向设成“Unit Cell”Y方向设成“Open (add space)”Z方向设成“Open (add space)”。X方向的Unit Cell需要设置相位差这个在求解器设置里做。第五步设置端口。在X方向的两个端面上设置“Open (add space)”端口。端口的尺寸要覆盖整个共面波导的横截面包括中心导带、槽缝和部分地平面。4.2 求解器参数的具体配置在“Eigenmode”求解器设置里找到“Phase”选项卡。设置“Phase from”为0度“Phase to”为180度“Phase steps”为37对应5度步长。在“Modes”选项卡里设置“Number of modes”为3。在“Accuracy”选项卡里设置“Accuracy”为-40dB。网格设置里选择“Tet mesh”或者“Hex mesh”。对于共面波导我建议用“Tet mesh”因为槽缝区域的几何形状比较复杂四面体网格适应性更好。在“Mesh Refinement”里对槽缝区域设置局部加密网格尺寸设为槽缝宽度的五分之一。材料设置里如果你要用频率相关的介电常数在“Material”里选择“Frequency dependent”然后导入或者手动输入数据。如果只是练手用常数介电常数就行。4.3 运行求解与后处理提取色散曲线点击“Start”运行求解。求解时间取决于网格数量和模式数量。对于上面这个尺寸的模型大概几分钟到十几分钟。如果跑了半小时还没结束检查一下网格是不是太密了或者模式数量设太多了。求解完成后在“Results”里找到“Eigenmode”结果。CST会列出每个相位点对应的本征频率。你需要把这些数据导出然后手动计算相位常数。相位常数的计算公式是β (相位差 × π) / (180 × 周期长度)其中相位差的单位是度周期长度的单位是米β的单位是弧度每米。如果你想要更直观的单位可以转换成度每毫米或者弧度每毫米。把频率作为横轴相位常数作为纵轴画出来的曲线就是色散曲线。如果你有多个模式每个模式画一条线。主模的曲线应该从原点附近开始随着频率增加单调上升。高阶模的曲线会在某个频率突然出现然后也单调上升。后处理的时候CST可以直接画“Eigenmode”的“Frequency vs Phase”曲线但那个曲线的横轴是相位差不是相位常数。你需要自己转换一下。我一般会把数据导出到Excel或者Python里做转换和画图。Python里用matplotlib画比CST自带的画图工具灵活得多。4.4 特征阻抗的提取方法色散曲线只给了相位常数特征阻抗还需要另外提取。CST的本征模求解器可以计算端口的特征阻抗但需要你在端口设置里指定“Impedance”计算。具体来说在端口的后处理里选择“Line Impedance”或者“Wave Impedance”CST会根据端口的场分布和功率流来计算。对于共面波导特征阻抗的定义有多种功率-电流定义、功率-电压定义、电压-电流定义。不同的定义得到的阻抗值略有差异。CST默认用的是功率-电流定义。如果你要和别的软件对比注意确认定义是否一致。我一般会提取两个阻抗一个是“Characteristic Impedance”一个是“Port Impedance”。前者是传输线的固有阻抗后者是端口看到的阻抗。对于匹配良好的共面波导两者应该接近。如果差异很大说明端口设置有问题可能是端口尺寸不合适或者边界条件不对。5. 常见问题与排查技巧实录5.1 色散曲线不光滑或者跳变这是最常见的问题。曲线跳变的原因通常有三个模式识别错误、网格太粗、相位步长太大。模式识别错误的表现是曲线在某个频率突然从一个模式跳到另一个模式。排查方法是把跳变点附近的场分布调出来看看场型是不是变了。如果变了那就是模式识别问题。解决办法是在求解器设置里增加模式数量然后在后处理里手动选择你要跟踪的模式。网格太粗的表现是曲线整体不光滑有锯齿。排查方法是把网格加密一倍重新跑一遍看看曲线是否变光滑。如果变光滑了那就是网格问题。解决办法是加密网格尤其是槽缝区域。相位步长太大的表现是曲线在某些区域有直线段不够平滑。排查方法是把步长从10度改成5度重新跑。如果曲线变平滑了那就是步长问题。5.2 求解器报错“Out of memory”本征模求解器对内存的需求比时域求解器大得多。如果模型尺寸大、网格密、模式多内存很容易爆。解决办法有几个减少模式数量、降低网格密度、缩小模型尺寸、用64位求解器。我遇到过一次模型尺寸是10毫米×5毫米×0.5毫米网格尺寸是1微米模式数量是5结果内存直接爆了。后来把网格尺寸改成5微米模式数量改成3就顺利跑完了。所以网格尺寸不要盲目追求小够用就行。5.3 特征阻抗提取结果异常特征阻抗提取异常的表现是阻抗值随频率剧烈波动或者阻抗值明显偏离预期。原因可能是端口尺寸不合适、端口模式设置错误、或者边界条件不对。排查方法是先检查端口的场分布看看是不是和预期的共面波导模式一致。如果端口场分布不对调整端口尺寸或者端口类型。如果端口场分布对但阻抗还是不对检查边界条件尤其是Y和Z方向的边界确保没有把场截断。还有一个可能的原因是金属厚度没有考虑。如果你用PEC金属厚度为零特征阻抗会比实际值略低。如果你用Lossy metal金属厚度设对了阻抗会更接近实际。这个差异在低频段不大但在高频段可能达到几个欧姆。5.4 常见问题速查表问题现象可能原因排查方法解决办法曲线跳变模式识别错误检查跳变点场分布增加模式数量手动选择曲线不光滑网格太粗加密网格重跑槽缝区域局部加密曲线有直线段相位步长太大减小步长重跑步长改为5度内存不足网格太密或模式太多检查网格尺寸和模式数降低网格密度减少模式阻抗波动大端口尺寸不合适检查端口场分布调整端口尺寸或类型阻抗偏低金属厚度为零检查材料设置改用Lossy metal求解时间过长精度设太高检查Accuracy设置降到-40dB曲线只画了一小段相位差单位错误检查相位设置确认单位是度注意以上排查方法适用于大多数情况但如果你的结构有特殊之处比如多层介质或者复杂的金属图案可能需要针对具体情况调整。5.5 独家避坑经验第一个经验先跑一个粗模型确认流程再跑精细模型。我刚开始做的时候总想一步到位把网格设得很密结果跑了两个小时发现边界条件设错了全部重来。后来我学乖了先用粗网格跑一遍确认边界条件、端口设置、材料参数都没问题再加密网格跑精确结果。粗模型可能只要几分钟但能帮你省下几个小时。第二个经验保存每个相位点的场分布。CST默认可能只保存最后一个相位点的场分布。如果你要分析模式识别问题需要每个相位点的场分布。在求解器设置里把“Store all modes”或者“Store fields”打开。这样后处理的时候可以逐个查看。第三个经验用Python脚本批量处理数据。CST可以导出Touchstone格式或者CSV格式的数据。如果你要扫描很多相位点手动处理数据很累。写一个Python脚本自动读取CSV计算相位常数画图保存。这个脚本一次写好以后每次跑完直接运行省时省力。第四个经验注意周期长度的选择。周期长度不能太长否则高阶模式会在周期内出现导致色散曲线混叠。周期长度也不能太短否则网格划分困难计算精度下降。一般来说周期长度取中心频率对应波长的十分之一到五分之一比较合适。比如中心频率是10GHz自由空间波长是30毫米介质中波长大概是15毫米周期长度取1.5到3毫米。6. 结果验证与扩展应用6.1 如何验证色散曲线的正确性跑完色散曲线怎么知道结果是对的有几个验证方法。第一个方法是低频极限验证。在频率趋近于零的时候共面波导的相位常数应该趋近于零而且斜率应该等于介质中光速的倒数。你可以把低频段的曲线拟合一条直线看看斜率是不是和理论值一致。如果偏差很大说明模型有问题。第二个方法是与解析公式对比。共面波导有准静态解析公式可以算出低频下的特征阻抗和有效介电常数。你可以用这些公式算一个参考值和CST的结果对比。如果差异在5%以内说明模型基本正确。如果差异很大检查材料参数和几何尺寸。第三个方法是与实测数据对比。如果你有实测的S参数可以通过S参数反推相位常数然后和CST的色散曲线对比。这个方法最可靠但需要实测数据。第四个方法是能量收敛性检查。在CST里本征模求解器会给出每个模式的能量误差。如果能量误差大于1%说明网格不够密结果不可靠。你需要加密网格重新跑。6.2 色散曲线在设计中的实际应用色散曲线跑出来之后怎么用举几个实际例子。如果你在设计一个共面波导定向耦合器色散曲线告诉你两个端口的相位差随频率的变化。你可以根据这个设计耦合器的长度使得在目标频段内相位差接近90度。如果你在设计一个共面波导天线色散曲线告诉你天线的工作频段和带宽。你可以根据色散曲线的斜率判断群速度进而估算天线的辐射效率。如果你在设计一个共面波导滤波器色散曲线告诉你通带和阻带的位置。你可以根据色散曲线的截止频率设计滤波器的尺寸。6.3 从单模到多模的扩展上面讲的都是单模色散曲线。如果你要分析多模共面波导比如耦合共面波导或者多层共面波导色散曲线会变得更复杂。多个模式的曲线会交叉、耦合形成模式耦合区域。分析多模色散曲线的时候关键是模式识别。你需要根据场分布来判断哪个模式是偶模哪个是奇模。CST里可以画每个模式的电场和磁场分布你可以根据对称性来判断。多模色散曲线的另一个关键是模式耦合。当两个模式的频率接近时它们会耦合曲线会出现“避免交叉”现象。这个现象在耦合共面波导中很常见。如果你看到两条曲线在某个频率附近互相排斥那就是模式耦合。6.4 参数扫描与优化色散曲线跑通之后你可以做参数扫描。比如改变槽缝宽度看看色散曲线怎么变。改变介质厚度看看色散曲线怎么变。CST里有“Parameter Sweep”功能可以自动扫描多个参数组合。参数扫描的时候注意每次扫描都要重新计算本征模。如果参数变化很大模式顺序可能改变你需要重新识别模式。我一般会在扫描之前先跑几个单点确认模式顺序稳定再批量扫描。优化的时候你可以把色散曲线的目标值设成优化目标比如“在10GHz到20GHz频段内相位常数偏差小于5%”。CST的优化器会自动调整参数使得色散曲线逼近目标。不过本征模求解器的优化速度比较慢因为每次迭代都要重新求解本征模。如果优化变量很多可能需要很长时间。提示参数扫描和优化的时候建议先用粗网格跑一遍确定大致趋势再用细网格跑精确结果。这样可以节省大量时间。7. 一些个人体会做共面波导色散曲线这件事说难不难说简单也不简单。难的是细节简单的是流程。流程就是建模、设边界、设求解器、跑、后处理五步。但每一步都有坑坑踩多了自然就熟了。我刚开始做的时候最头疼的是模式识别。跑出来的曲线乱七八糟根本分不清哪个是主模。后来我发现先跑一个单点的本征模分析把场分布画出来记住主模的场型然后再跑扫描这个问题就解决了。主模的电场在槽缝里是最强的高阶模的电场分布更复杂有多个瓣。你只要记住这个特征后处理的时候一眼就能认出来。另一个体会是不要迷信自动网格。CST的自动网格对于简单结构还行对于共面波导这种场集中在小区域的结构自动网格往往不够密。我一般会手动加密槽缝区域虽然麻烦一点但结果可靠得多。最后多动手少看书。色散曲线的理论在书上讲得很清楚但软件操作是另一回事。你看十遍书不如自己跑一遍。跑错了排查改设置再跑这个过程才是真正学东西的时候。我跑废了十几个模型才把流程跑通。你如果刚开始不要怕跑错跑错了才有经验。