1. 共面波导色散仿真为什么值得单独拎出来讲共面波导CPW在射频前端、毫米波封装互连和片上传输线里出现频率极高原因很直接信号线和地平面共面方便探针台在片测试也方便和贴片元件串联并联。但很多人第一次在CST里跑CPW色散特性时都会遇到同一个尴尬——按常规方法建一段传输线两端加波导端口扫频跑完发现高频段相位常数曲线开始抖或者端口模式里混进了不该有的高次模最后提取出来的有效介电常数和理论值对不上。这个问题的根源不在模型画得对不对而在于仿真域的长度和边界条件选择。CPW的色散特性本质上是“无限长均匀传输线”的传播常数随频率的变化你如果只截取有限长度的一段两端端口引入的反射和模式不纯会直接污染相位信息。更麻烦的是当频率升高到毫米波频段CPW的横向尺寸和波长可比高次模和辐射损耗开始出现有限长度的模型根本捕捉不到这些效应。所以正确做法是用周期性边界条件Unit Cell把一段CPW当作无限周期结构的一个单元来处理配合JDMJacobi-Davidson Method求解器做本征模求解直接拿到传播常数和衰减常数。这套流程在CST里不算复杂但有几个关键设置点如果搞错结果会差得离谱。下面我把整个思路、操作细节和踩过的坑完整拆一遍。2. 整体方案设计与工具选型逻辑2.1 为什么选周期性边界而不是波导端口扫频常规的波导端口扫频法Time Domain或Frequency Domain本质上是S参数法你建一段线两端加端口扫频得到S21然后通过相位提取传播常数。这个方法在低频段比如X波段以下够用但到了毫米波频段有几个硬伤。第一端口模式纯度问题。CPW在高频时除了主模偶模还会激励起奇模、表面波模和辐射模。波导端口默认只激励主模但传输线本身的不连续性比如你为了加端口而延伸的馈线会激发高次模这些模在端口处没有被完全吸收反射回来就污染了S21的相位。第二长度选择的两难。线太短相位变化小提取误差大线太长损耗导致S21幅度太小信噪比恶化。而且色散特性本身是每单位长度的相位变化你用有限长度提取本质上是在做差分近似频率越高误差越大。第三辐射损耗无法分离。CPW在高频时辐射损耗不可忽略S参数法只能拿到总损耗没法把导体损耗、介质损耗和辐射损耗分开。周期性边界本征模求解的思路完全不同我不关心波怎么进去怎么出来我直接问“这个无限周期结构支持什么样的传播模式”。CST的本征模求解器会直接给出每个模式的传播常数通过相位因子和品质因数通过衰减物理意义清晰没有端口反射的干扰。2.2 JDM求解器在什么场景下必须用CST的本征模求解器有几种算法可选默认的AKSArnoldi Krylov Subspace在大多数情况下够用但CPW色散仿真有几个特点让它容易翻车。CPW的横截面通常包含高对比度介质比如硅衬底εr11.9和空气εr1以及薄金属层几微米甚至亚微米。这种结构导致本征模问题的矩阵条件数很差AKS求解器需要大量的迭代才能收敛而且容易漏掉某些模式。更关键的是CPW的色散仿真通常需要扫频每个频点都要重新求解本征模问题如果单点求解就慢扫频下来时间不可接受。JDM求解器的优势在于它专门针对大规模稀疏矩阵的极端特征值问题做了优化对于需要精确捕捉多个模式包括高次模的场景收敛速度比AKS快很多。而且JDM支持目标频率附近的模式搜索你可以指定一个频率范围它只找你关心的那几个模式不用把整个谱都算出来。实测下来对于典型的CPW结构衬底厚度100-500μm信号线宽度10-50μm频率到110GHzJDM求解器比AKS快3-5倍而且模式识别更稳定。如果你的模型里还有薄电阻层或者各向异性材料JDM的优势会更明显。2.3 周期单元的长度怎么定这是整个仿真里最容易拍脑袋决定、也最容易出错的地方。周期单元的长度也就是你在传播方向上截取的长度不能随便设它直接决定了你能捕捉到的最高频率和模式数量。基本原则是周期长度必须小于最高分析频率对应的半波长。如果周期长度大于半波长布里渊区折叠会导致模式混叠你看到的传播常数曲线会出现虚假的交叉和分裂。具体计算假设最高分析频率是110GHz衬底是硅εr11.9CPW的有效介电常数大约在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μm10倍这样网格划分更规整数值误差更小。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周期长度。材料设成硅εr11.9损耗角正切0.001。第二步画CPW金属。信号线在衬底上表面宽度20μm厚度用Ohmic Sheet代替零厚度设表面电阻。地平面在信号线两侧间距Gap15μ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步长5GHz21个频点。注意步长不能太大否则模式跟踪会丢模式。5GHz步长对于CPW色散曲线够用如果你要看精细的色散特性比如模式交叉点附近步长要降到1GHz。模式跟踪设置打开“Mode Tracking”参考模式选“Previous Frequency”。这样CST会按频率顺序自动匹配模式避免模式顺序跳变。收敛容差设1e-4最大迭代次数设150。先跑一个频点试试看收敛需要多少次迭代。如果50次以内收敛说明设置没问题如果超过100次检查网格是不是太粗或者材料参数是不是有突变。4.3 结果提取与验证跑完扫频后在“Results”里看“Eigenmode”结果。找到“Phase Factor”或者“Beta”曲线这就是传播常数随频率的变化。提取数据后用上面的公式算εeff。对于典型的CPW信号线20μmGap 15μm硅衬底εeff在低频时应该接近(11.91)/26.45高频时因为场更集中在衬底里εeff会略微上升到110GHz时大约6.8-7.0。如果你的结果偏离这个范围超过10%说明模型或者设置有问题。衰减常数从Q值提取。注意本征模求解器的Q值是无载Q包含了所有损耗。如果你只关心导体损耗把介质损耗角正切设成0再跑一次对比两次结果就能分离出介质损耗。验证方法用同样的结构在另一个仿真工具比如HFSS的本征模求解器里跑一遍对比传播常数曲线。如果两条曲线在低频段重合高频段有偏差说明你的网格或者边界设置在某个频段出了问题。5. 常见问题与排查技巧实录5.1 模式识别混乱怎么办这是最常见的问题扫频跑完传播常数曲线在某些频点突然跳变或者出现不物理的交叉。原因通常是模式跟踪失败CST把不同模式的相位因子混在一起了。排查步骤检查模式数量是不是设少了。如果实际有4个模式你只算3个第4个模式会污染前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 常见问题速查表问题现象可能原因排查方法解决措施传播常数曲线跳变模式跟踪失败检查模式数量和步长增加模式数减小步长高频不收敛网格太粗或场奇异查看迭代次数和场分布加密边缘网格放宽容差εeff偏低地平面太窄或衬底太薄对比解析公式加宽地平面增加衬底厚度εeff虚高周期长度过大计算半波长减小周期长度到λ/10扫频太慢迭代次数多查看收敛曲线用Restart放宽容差Q值异常低辐射损耗大检查截断边界加大截断距离用PML6. 实操心得与进阶技巧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”收敛更稳但每步更慢。实测下来对于硅基CPWKrylov Size15Shift目标频率PreconditionerILU这个组合在速度和稳定性之间平衡最好。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%”这样不同参数之间的网格精度一致结果可比性强。
CST共面波导色散仿真：周期性边界与JDM求解器实战指南
1. 共面波导色散仿真为什么值得单独拎出来讲共面波导CPW在射频前端、毫米波封装互连和片上传输线里出现频率极高原因很直接信号线和地平面共面方便探针台在片测试也方便和贴片元件串联并联。但很多人第一次在CST里跑CPW色散特性时都会遇到同一个尴尬——按常规方法建一段传输线两端加波导端口扫频跑完发现高频段相位常数曲线开始抖或者端口模式里混进了不该有的高次模最后提取出来的有效介电常数和理论值对不上。这个问题的根源不在模型画得对不对而在于仿真域的长度和边界条件选择。CPW的色散特性本质上是“无限长均匀传输线”的传播常数随频率的变化你如果只截取有限长度的一段两端端口引入的反射和模式不纯会直接污染相位信息。更麻烦的是当频率升高到毫米波频段CPW的横向尺寸和波长可比高次模和辐射损耗开始出现有限长度的模型根本捕捉不到这些效应。所以正确做法是用周期性边界条件Unit Cell把一段CPW当作无限周期结构的一个单元来处理配合JDMJacobi-Davidson Method求解器做本征模求解直接拿到传播常数和衰减常数。这套流程在CST里不算复杂但有几个关键设置点如果搞错结果会差得离谱。下面我把整个思路、操作细节和踩过的坑完整拆一遍。2. 整体方案设计与工具选型逻辑2.1 为什么选周期性边界而不是波导端口扫频常规的波导端口扫频法Time Domain或Frequency Domain本质上是S参数法你建一段线两端加端口扫频得到S21然后通过相位提取传播常数。这个方法在低频段比如X波段以下够用但到了毫米波频段有几个硬伤。第一端口模式纯度问题。CPW在高频时除了主模偶模还会激励起奇模、表面波模和辐射模。波导端口默认只激励主模但传输线本身的不连续性比如你为了加端口而延伸的馈线会激发高次模这些模在端口处没有被完全吸收反射回来就污染了S21的相位。第二长度选择的两难。线太短相位变化小提取误差大线太长损耗导致S21幅度太小信噪比恶化。而且色散特性本身是每单位长度的相位变化你用有限长度提取本质上是在做差分近似频率越高误差越大。第三辐射损耗无法分离。CPW在高频时辐射损耗不可忽略S参数法只能拿到总损耗没法把导体损耗、介质损耗和辐射损耗分开。周期性边界本征模求解的思路完全不同我不关心波怎么进去怎么出来我直接问“这个无限周期结构支持什么样的传播模式”。CST的本征模求解器会直接给出每个模式的传播常数通过相位因子和品质因数通过衰减物理意义清晰没有端口反射的干扰。2.2 JDM求解器在什么场景下必须用CST的本征模求解器有几种算法可选默认的AKSArnoldi Krylov Subspace在大多数情况下够用但CPW色散仿真有几个特点让它容易翻车。CPW的横截面通常包含高对比度介质比如硅衬底εr11.9和空气εr1以及薄金属层几微米甚至亚微米。这种结构导致本征模问题的矩阵条件数很差AKS求解器需要大量的迭代才能收敛而且容易漏掉某些模式。更关键的是CPW的色散仿真通常需要扫频每个频点都要重新求解本征模问题如果单点求解就慢扫频下来时间不可接受。JDM求解器的优势在于它专门针对大规模稀疏矩阵的极端特征值问题做了优化对于需要精确捕捉多个模式包括高次模的场景收敛速度比AKS快很多。而且JDM支持目标频率附近的模式搜索你可以指定一个频率范围它只找你关心的那几个模式不用把整个谱都算出来。实测下来对于典型的CPW结构衬底厚度100-500μm信号线宽度10-50μm频率到110GHzJDM求解器比AKS快3-5倍而且模式识别更稳定。如果你的模型里还有薄电阻层或者各向异性材料JDM的优势会更明显。2.3 周期单元的长度怎么定这是整个仿真里最容易拍脑袋决定、也最容易出错的地方。周期单元的长度也就是你在传播方向上截取的长度不能随便设它直接决定了你能捕捉到的最高频率和模式数量。基本原则是周期长度必须小于最高分析频率对应的半波长。如果周期长度大于半波长布里渊区折叠会导致模式混叠你看到的传播常数曲线会出现虚假的交叉和分裂。具体计算假设最高分析频率是110GHz衬底是硅εr11.9CPW的有效介电常数大约在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μm10倍这样网格划分更规整数值误差更小。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周期长度。材料设成硅εr11.9损耗角正切0.001。第二步画CPW金属。信号线在衬底上表面宽度20μm厚度用Ohmic Sheet代替零厚度设表面电阻。地平面在信号线两侧间距Gap15μ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步长5GHz21个频点。注意步长不能太大否则模式跟踪会丢模式。5GHz步长对于CPW色散曲线够用如果你要看精细的色散特性比如模式交叉点附近步长要降到1GHz。模式跟踪设置打开“Mode Tracking”参考模式选“Previous Frequency”。这样CST会按频率顺序自动匹配模式避免模式顺序跳变。收敛容差设1e-4最大迭代次数设150。先跑一个频点试试看收敛需要多少次迭代。如果50次以内收敛说明设置没问题如果超过100次检查网格是不是太粗或者材料参数是不是有突变。4.3 结果提取与验证跑完扫频后在“Results”里看“Eigenmode”结果。找到“Phase Factor”或者“Beta”曲线这就是传播常数随频率的变化。提取数据后用上面的公式算εeff。对于典型的CPW信号线20μmGap 15μm硅衬底εeff在低频时应该接近(11.91)/26.45高频时因为场更集中在衬底里εeff会略微上升到110GHz时大约6.8-7.0。如果你的结果偏离这个范围超过10%说明模型或者设置有问题。衰减常数从Q值提取。注意本征模求解器的Q值是无载Q包含了所有损耗。如果你只关心导体损耗把介质损耗角正切设成0再跑一次对比两次结果就能分离出介质损耗。验证方法用同样的结构在另一个仿真工具比如HFSS的本征模求解器里跑一遍对比传播常数曲线。如果两条曲线在低频段重合高频段有偏差说明你的网格或者边界设置在某个频段出了问题。5. 常见问题与排查技巧实录5.1 模式识别混乱怎么办这是最常见的问题扫频跑完传播常数曲线在某些频点突然跳变或者出现不物理的交叉。原因通常是模式跟踪失败CST把不同模式的相位因子混在一起了。排查步骤检查模式数量是不是设少了。如果实际有4个模式你只算3个第4个模式会污染前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 常见问题速查表问题现象可能原因排查方法解决措施传播常数曲线跳变模式跟踪失败检查模式数量和步长增加模式数减小步长高频不收敛网格太粗或场奇异查看迭代次数和场分布加密边缘网格放宽容差εeff偏低地平面太窄或衬底太薄对比解析公式加宽地平面增加衬底厚度εeff虚高周期长度过大计算半波长减小周期长度到λ/10扫频太慢迭代次数多查看收敛曲线用Restart放宽容差Q值异常低辐射损耗大检查截断边界加大截断距离用PML6. 实操心得与进阶技巧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”收敛更稳但每步更慢。实测下来对于硅基CPWKrylov Size15Shift目标频率PreconditionerILU这个组合在速度和稳定性之间平衡最好。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%”这样不同参数之间的网格精度一致结果可比性强。
