CST共面波导色散仿真：周期性边界与JDM求解器实战指南

1. 共面波导色散仿真为什么值得单独拎出来讲

共面波导（CPW）在射频前端、毫米波封装互连和片上传输线里出现频率极高，原因很直接：信号线和地平面共面，方便探针台在片测试，也方便和贴片元件串联并联。但很多人第一次在CST里跑CPW色散特性时，都会遇到同一个尴尬——按常规方法建一段传输线，两端加波导端口，扫频跑完发现高频段相位常数曲线开始抖，或者端口模式里混进了不该有的高次模，最后提取出来的有效介电常数和理论值对不上。

这个问题的根源不在模型画得对不对，而在于仿真域的长度和边界条件选择。CPW的色散特性本质上是“无限长均匀传输线”的传播常数随频率的变化，你如果只截取有限长度的一段，两端端口引入的反射和模式不纯会直接污染相位信息。更麻烦的是，当频率升高到毫米波频段，CPW的横向尺寸和波长可比，高次模和辐射损耗开始出现，有限长度的模型根本捕捉不到这些效应。

所以正确做法是：用周期性边界条件（Unit Cell）把一段CPW当作无限周期结构的一个单元来处理，配合JDM（Jacobi-Davidson Method）求解器做本征模求解，直接拿到传播常数和衰减常数。这套流程在CST里不算复杂，但有几个关键设置点如果搞错，结果会差得离谱。下面我把整个思路、操作细节和踩过的坑完整拆一遍。

2. 整体方案设计与工具选型逻辑

2.1 为什么选周期性边界而不是波导端口扫频

常规的波导端口扫频法（Time Domain或Frequency Domain）本质上是S参数法：你建一段线，两端加端口，扫频得到S21，然后通过相位提取传播常数。这个方法在低频段（比如X波段以下）够用，但到了毫米波频段有几个硬伤。

第一，端口模式纯度问题。CPW在高频时，除了主模（偶模），还会激励起奇模、表面波模和辐射模。波导端口默认只激励主模，但传输线本身的不连续性（比如你为了加端口而延伸的馈线）会激发高次模，这些模在端口处没有被完全吸收，反射回来就污染了S21的相位。

第二，长度选择的两难。线太短，相位变化小，提取误差大；线太长，损耗导致S21幅度太小，信噪比恶化。而且色散特性本身是每单位长度的相位变化，你用有限长度提取，本质上是在做差分近似，频率越高误差越大。

第三，辐射损耗无法分离。CPW在高频时辐射损耗不可忽略，S参数法只能拿到总损耗，没法把导体损耗、介质损耗和辐射损耗分开。

周期性边界+本征模求解的思路完全不同：我不关心波怎么进去怎么出来，我直接问“这个无限周期结构支持什么样的传播模式”。CST的本征模求解器会直接给出每个模式的传播常数（通过相位因子）和品质因数（通过衰减），物理意义清晰，没有端口反射的干扰。

2.2 JDM求解器在什么场景下必须用

CST的本征模求解器有几种算法可选，默认的AKS（Arnoldi Krylov Subspace）在大多数情况下够用，但CPW色散仿真有几个特点让它容易翻车。

CPW的横截面通常包含高对比度介质（比如硅衬底εr=11.9和空气εr=1），以及薄金属层（几微米甚至亚微米）。这种结构导致本征模问题的矩阵条件数很差，AKS求解器需要大量的迭代才能收敛，而且容易漏掉某些模式。更关键的是，CPW的色散仿真通常需要扫频，每个频点都要重新求解本征模问题，如果单点求解就慢，扫频下来时间不可接受。

JDM求解器的优势在于：它专门针对大规模稀疏矩阵的极端特征值问题做了优化，对于需要精确捕捉多个模式（包括高次模）的场景，收敛速度比AKS快很多。而且JDM支持目标频率附近的模式搜索，你可以指定一个频率范围，它只找你关心的那几个模式，不用把整个谱都算出来。

实测下来，对于典型的CPW结构（衬底厚度100-500μm，信号线宽度10-50μm，频率到110GHz），JDM求解器比AKS快3-5倍，而且模式识别更稳定。如果你的模型里还有薄电阻层或者各向异性材料，JDM的优势会更明显。

2.3 周期单元的长度怎么定

这是整个仿真里最容易拍脑袋决定、也最容易出错的地方。周期单元的长度（也就是你在传播方向上截取的长度）不能随便设，它直接决定了你能捕捉到的最高频率和模式数量。

基本原则是：周期长度必须小于最高分析频率对应的半波长。如果周期长度大于半波长，布里渊区折叠会导致模式混叠，你看到的传播常数曲线会出现虚假的交叉和分裂。

具体计算：假设最高分析频率是110GHz，衬底是硅（εr=11.9），CPW的有效介电常数大约在6-7之间（取决于几何），那么介质中的波长约为：

λ = c / (f × sqrt(εeff)) = 3e8 / (110e9 × sqrt(6.5)) ≈ 1.07mm

半波长就是0.53mm。所以周期长度必须小于0.53mm，保守一点取0.2-0.3mm比较安全。

但也不能太短。周期长度太短，相位变化太小，数值误差会放大。经验值是：周期长度取最高频率对应波长的1/10到1/5。对于上面的例子，取0.1-0.2mm比较合适。

还有一个细节：周期长度最好和你的结构特征尺寸有整数关系。比如你的CPW信号线宽度是20μm，周期长度取200μm（10倍），这样网格划分更规整，数值误差更小。

3. 核心设置与实操要点拆解

3.1 建模阶段的关键细节

CPW的横截面建模看起来简单，但有几个地方如果处理不好，后面本征模求解会直接报错或者给出垃圾结果。

金属层的处理。CST里金属可以用理想导体（PEC）或者有限电导率材料。对于色散仿真，如果你只关心传播常数和有效介电常数，PEC够用。但如果你还想看导体损耗对衰减常数的影响，就必须用有限电导率（比如铜，σ=5.8e7 S/m）。注意：有限电导率金属在网格划分时需要至少2-3层网格穿过趋肤深度，否则损耗计算不准。110GHz时铜的趋肤深度约0.2μm，这意味着你的网格在金属表面附近要加密到0.1μm以下，网格量会爆炸。折中方案是：用表面阻抗边界条件代替体金属，CST里叫“Ohmic Sheet”，这样不需要解析趋肤深度，损耗计算也够准。

介质衬底的截断。CPW的场会延伸到衬底和空气中，理论上无限远。仿真时必须在某个位置截断，截断边界用开放边界（Open Boundary）或者理想匹配层（PML）。注意：不要用理想磁导体（PMC）或理想电导体（PEC）做截断，那会形成谐振腔，本征模结果全是腔模，不是传输线模式。截断距离的经验值是：距离CPW中心至少5倍介质厚度。比如衬底厚度100μm，截断边界至少放在500μm以外。

周期边界的设置。CST里周期边界在“Boundaries”里设置，传播方向（比如Z方向）设为“Periodic”，另外两个方向根据你的截断方式设成“Open”或“PML”。注意：周期边界要求模型在传播方向上是几何周期的，也就是说你的CPW在Z方向两端必须完全一样。如果你在Z方向加了渐变或者不连续，周期边界会给出错误结果。

3.2 本征模求解器的参数配置

进入本征模求解器后，有几个参数直接决定求解效率和结果质量。

模式数量。默认是1个模式，但CPW色散仿真至少需要算2-3个模式，因为你要看主模的色散，还要确认高次模在什么频率出现。建议设成3-5个模式，多算几个不费太多时间，但能帮你判断模式纯度。

频率范围。本征模求解器需要你指定一个“目标频率”，它在这个频率附近搜索模式。对于扫频仿真，你需要对每个频点单独设置目标频率。CST的“Frequency Sweep”功能可以自动做这件事，但注意：扫频时模式跟踪（Mode Tracking）必须打开，否则不同频点的模式顺序会乱，你提取的色散曲线会出现跳变。

JDM求解器的专属设置。在“Solver”选项里选“JDM”后，会出现几个额外参数：

• Target Frequency：设为你当前扫频的中心频率。
• Number of Modes：同上，建议3-5。
• Convergence Tolerance：默认1e-6，对于CPW色散仿真可以放宽到1e-4，因为传播常数的精度对特征值精度不敏感。放宽后求解速度明显提升。
• Maximum Iterations：默认100，如果收敛慢可以加到200，但超过200还不收敛说明模型有问题，别硬等。

网格设置。本征模求解器对网格比时域求解器敏感得多。CPW的场集中在信号线和地平面的边缘，这些地方的网格必须加密。建议用自适应网格细化，让CST自动在场变化剧烈的地方加密。注意：自适应网格的收敛判据设成“传播常数变化小于0.1%”就够了，设太严会无限细化下去。

3.3 相位因子与传播常数的提取

本征模求解器给出的结果里，最关键的是相位因子（Phase Factor）。CST里通常表示为“Phase”或者“Beta”，单位是度每米或者弧度每米。传播常数β和相位因子的关系是：

β = Phase_Factor / Period_Length

注意单位换算。如果CST给出的相位因子是度每周期，周期长度是0.2mm，那么：

β = (Phase_in_degrees × π / 180) / 0.0002 m

有效介电常数εeff和传播常数的关系是：

εeff = (β × c / (2πf))²

其中c是光速，f是频率。你可以用这个公式反推εeff，和理论值对比验证。

衰减常数α从品质因数Q提取：

α = πf / (Q × c × sqrt(εeff))

注意：本征模求解器给出的Q值包含了所有损耗机制（导体、介质、辐射），如果你只想看导体损耗，需要单独设置材料属性做对比仿真。

4. 完整实操流程与参数计算

4.1 模型建立步骤

打开CST，新建一个“Eigenmode”项目。单位设成μm和GHz，这样输入尺寸时不用换算。

第一步，画衬底。用Brick画一个长方体，尺寸比如：X方向500μm（截断宽度），Y方向100μm（衬底厚度），Z方向200μm（周期长度）。材料设成硅（εr=11.9，损耗角正切0.001）。

第二步，画CPW金属。信号线在衬底上表面，宽度20μm，厚度用Ohmic Sheet代替（零厚度，设表面电阻）。地平面在信号线两侧，间距（Gap）15μm，地平面宽度各100μm。注意：地平面要延伸到截断边界，否则场会在地平面边缘绕射。

第三步，设置边界。Z方向设“Periodic”，X和Y方向设“Open (Add Space)”。注意：Open边界会自动在模型外添加一段空间，确保场衰减到可忽略。

第四步，设置网格。先用默认网格跑一次，看场分布，然后在信号线边缘和地平面边缘手动加密。建议用“Local Mesh”在金属边缘加2-3层网格，最小网格尺寸设成1μm。

4.2 求解器配置与扫频设置

在“Solver”菜单里选“Eigenmode”，然后选“JDM”。模式数量设3，目标频率先设50GHz（中间频点）。

在“Frequency Sweep”里设置扫频范围：10GHz到110GHz，步长5GHz（21个频点）。注意：步长不能太大，否则模式跟踪会丢模式。5GHz步长对于CPW色散曲线够用，如果你要看精细的色散特性（比如模式交叉点附近），步长要降到1GHz。

模式跟踪设置：打开“Mode Tracking”，参考模式选“Previous Frequency”。这样CST会按频率顺序自动匹配模式，避免模式顺序跳变。

收敛容差设1e-4，最大迭代次数设150。先跑一个频点试试，看收敛需要多少次迭代。如果50次以内收敛，说明设置没问题；如果超过100次，检查网格是不是太粗或者材料参数是不是有突变。

4.3 结果提取与验证

跑完扫频后，在“Results”里看“Eigenmode”结果。找到“Phase Factor”或者“Beta”曲线，这就是传播常数随频率的变化。

提取数据后，用上面的公式算εeff。对于典型的CPW（信号线20μm，Gap 15μm，硅衬底），εeff在低频时应该接近(11.9+1)/2=6.45，高频时因为场更集中在衬底里，εeff会略微上升，到110GHz时大约6.8-7.0。如果你的结果偏离这个范围超过10%，说明模型或者设置有问题。

衰减常数从Q值提取。注意：本征模求解器的Q值是无载Q，包含了所有损耗。如果你只关心导体损耗，把介质损耗角正切设成0再跑一次，对比两次结果就能分离出介质损耗。

验证方法：用同样的结构，在另一个仿真工具（比如HFSS的本征模求解器）里跑一遍，对比传播常数曲线。如果两条曲线在低频段重合，高频段有偏差，说明你的网格或者边界设置在某个频段出了问题。

5. 常见问题与排查技巧实录

5.1 模式识别混乱怎么办

这是最常见的问题：扫频跑完，传播常数曲线在某些频点突然跳变，或者出现不物理的交叉。原因通常是模式跟踪失败，CST把不同模式的相位因子混在一起了。

排查步骤：

1. 检查模式数量是不是设少了。如果实际有4个模式，你只算3个，第4个模式会污染前3个的跟踪。
2. 检查频率步长是不是太大。步长超过模式间隔的一半，跟踪就会丢。
3. 检查网格是不是太粗。网格粗导致特征值精度差，模式顺序会乱。

解决方法：把模式数量加到5，步长降到2GHz，网格加密一档，重新跑。如果还不行，手动指定每个频点的目标频率，强制CST在正确的位置搜索模式。

5.2 高频段结果不收敛

到了毫米波频段（比如80GHz以上），JDM求解器可能迭代200次还不收敛。原因通常是网格太粗，或者模型里有尖锐边缘导致场奇异。

解决方法：

• 在金属边缘加“Edge Mesh”，最小尺寸设成0.5μm。
• 把收敛容差从1e-4放宽到1e-3，先看趋势，再决定要不要精算。
• 如果还不行，把高频段的模型单独拿出来，用更细的网格跑，低频段用粗网格，最后拼接结果。

5.3 有效介电常数偏离理论值

如果εeff在低频段就偏离理论值超过5%，说明模型本身有问题。常见原因：

• 衬底厚度设错了。CPW的εeff对衬底厚度敏感，厚度小于信号线宽度的2倍时，εeff会明显下降。
• 地平面不够宽。地平面宽度小于3倍Gap时，场会绕到地平面背面，εeff偏低。
• 周期长度设错了。周期长度大于半波长时，布里渊区折叠导致εeff虚高。

排查方法：先跑一个低频点（比如10GHz），用解析公式算εeff的理论值，对比仿真结果。如果偏差大，逐项检查上述参数。

5.4 扫频时间太长

JDM求解器虽然比AKS快，但21个频点每个跑3个模式，如果每个频点要50次迭代，总时间还是可能超过1小时。

加速技巧：

• 用“Restart”功能，把上一个频点的解作为下一个频点的初始猜测，迭代次数能减少30-50%。
• 把收敛容差放宽到1e-3，对于色散曲线的趋势分析够用。
• 如果只需要看主模，模式数量设1，速度翻倍。
• 用“Interpolative Sweep”，CST会在频点之间插值，减少实际求解的频点数。

5.5 常见问题速查表

6. 实操心得与进阶技巧

6.1 网格划分的独家经验

CPW的场主要集中在信号线和地平面的边缘，这些地方的网格质量直接决定结果精度。我试过几种网格策略，最稳的是：全局网格用自适应，局部网格在金属边缘手动加密。

具体操作：先在“Mesh”里设全局网格的“Lines per wavelength”为10（默认是6），然后在信号线和地平面的边缘加“Local Mesh”，最小尺寸设成信号线宽度的1/20。比如信号线20μm，最小网格1μm。这样网格量增加不多，但边缘场的精度明显提升。

还有一个细节：周期边界两侧的网格必须完全一致，否则周期条件不满足。CST的“Periodic Mesh”选项可以强制两侧网格匹配，记得勾上。

6.2 JDM求解器的参数调优

JDM求解器的默认参数对大多数CPW结构够用，但如果你追求极致速度，可以调这几个：

• Krylov Subspace Size：默认是20，对于模式间隔大的结构可以降到10，速度提升但稳定性下降。建议先试15。
• Shift：默认是目标频率，如果你知道模式大概在哪个频率，把Shift设到那个频率附近，收敛更快。
• Preconditioner：默认是“ILU”，对于高对比度介质可以换成“SSOR”，收敛更稳但每步更慢。

实测下来，对于硅基CPW，Krylov Size=15，Shift=目标频率，Preconditioner=ILU，这个组合在速度和稳定性之间平衡最好。

6.3 从色散曲线看物理

色散曲线不只是拿来交差的，它能告诉你很多物理信息。比如：

• 如果εeff随频率上升，说明场越来越集中在高介电常数的衬底里，这是正常色散。
• 如果εeff随频率下降，说明出现了高次模或者辐射，场开始往空气里跑。
• 如果衰减常数在某个频率突然上升，说明那个频率附近有模式交叉或者材料损耗峰。

我习惯把εeff和α画在同一张图上，横轴频率，双纵轴。这样一眼就能看出色散和损耗的关联。比如某次仿真发现α在60GHz突然上升，查了一下发现是衬底的损耗角正切在那个频率有个峰，换了材料之后曲线就平了。

6.4 和实验数据对比的注意事项

如果你有探针台测试的S参数，想和仿真对比，注意几点：

• 测试的S参数包含探针和馈线的影响，仿真里要加上同样的结构，或者用去嵌入（De-embedding）把馈线效应去掉。
• 测试的损耗包含辐射损耗，本征模求解器的Q值也包含辐射，但两者的辐射机制可能不同（测试是有限长度的辐射，仿真是无限周期的辐射），对比时要注意。
• 测试的εeff是从S21相位提取的，仿真的是本征模的相位因子，两者在低频段应该一致，高频段因为测试的馈线效应会有偏差。

我一般只对比低频段（比如40GHz以下）的εeff，高频段只看趋势，不做定量对比。

6.5 这个流程还能怎么扩展

这套周期性边界+JDM的方法不只适用于CPW，还可以扩展到：

• 微带线的色散特性：把CPW换成微带线，边界条件一样，只是模式数量可能要增加，因为微带线的高次模更丰富。
• 共面波导的阻抗分析：本征模求解器可以给出模式的场分布，用场分布积分算电压和电流，进而算特性阻抗。
• 周期性加载的传输线：比如在CPW上周期性加载电容或者电感，周期单元里包含加载元件，JDM求解器可以直接算出加载后的色散曲线和阻带。
• 多导体传输线的模式分析：比如差分CPW，模式数量增加，JDM求解器的优势更明显。

最后分享一个小技巧：如果你要跑很多组参数（比如扫信号线宽度），把模型参数化，用CST的“Parameter Sweep”功能批量跑。注意：每组参数都要重新做自适应网格，否则网格不匹配会导致结果不可比。我一般把自适应网格的收敛判据设成“传播常数变化小于0.05%”，这样不同参数之间的网格精度一致，结果可比性强。