Ansys Q3D 提取学习指南：从矩阵原理到互连工程工作流

摘要

Ansys Q3D Extractor 是一个面向三维互连结构的准静态场求解器，用于从几何模型提取等效电路参数矩阵，包括电容矩阵 (C)、电感矩阵 (L)、电阻矩阵 (R) 和电导矩阵 (G)。学习 Q3D 的重点在于理解输出矩阵的物理意义、源汇 (source/sink) 的定义对矩阵的影响，以及如何通过矩阵规约 (reduction) 构建最终的电路模型。本指南从矩阵物理定义入手，阐述源/汇建模与边界条件，再介绍矩阵规约与端口定义的原理；随后讨论影响精度的网格与求解器设置；最后围绕 BGA 焊球与 PCB 提取给出典型工作流和工程实践建议，并说明验证方案与常见误区。文中包含必要的公式与步骤说明，以帮助读者系统掌握 Q3D 建模与仿真方法。

1. 软件定位与方法边界

Q3D Extractor 并非通用全波分析器，而是一款寄生参数提取工具，适用于分析短链路、封装互连和局部电路的寄生影响。在频率较低（波长远大于结构尺寸）的场景下，麦克斯韦方程中的位移电流与电导电流可分离处理（电磁场近似为静电场或低频磁场），此时结构的行为可近似视为电容、电感和电阻网络。Q3D 正是在这一假设下工作：使用静电场求解器计算导体之间的电势-电荷关系，生成电容 (C) 矩阵和电导 (G) 矩阵；使用磁准静态求解器计算导体间的自感、互感及电阻，生成电感 (L) 矩阵和电阻 (R) 矩阵。这些矩阵在多端子体系中共同描述了互连的等效线性模型。一般而言，当结构长度小于波长的约 1/4 时，Q3D 提取的寄生参数在电路级模拟中是有效的；对于电磁传播、辐射或谐振主导的长链路，Q3D 就超出了适用范围，需要全波求解器或射线跟踪等方法。

注： Q3D 输出的是在指定源/汇 (source/sink) 布局下获得的完整 RLCG 矩阵，而不是已经构建好的 SPICE 网络。理解矩阵物理含义的关键是：先确定研究对象的电流回路和电压参考，再建立几何和电路端口。只有明确模型的电路拓扑，才能正确解释求解器给出的结果。

2. 矩阵物理定义

2.1 电容矩阵 (C) 与电导矩阵 (G)

静电求解得到的 电容矩阵 C 反映了导体间的静电耦合关系。一般定义为：

其中 Q_i 是第 i 个导体上的总电荷，V_j 是第 j 个导体相对于参考地的电压。电容矩阵 C 本质上是节点电容阵（Maxwell 矩阵）：对角元表示自身电容（自电容），非对角元表示互容（通常为负值）。如果介质有损耗，还会同时求出电导矩阵 G（与泄漏电流关联），可视为并联的泄漏通道，在计算过程中被加到节点导纳矩阵 Y = G + jωC。

需要强调的是，Maxwell 矩阵与电路网络层的对应关系并非一一映射。在等效电路中，通常采用 T 型或星型拓扑将矩阵转为元件（参考第 4 节）。例如，一个简单的两端网络，Maxwell 矩阵的非对角元对应反向的耦合电容（需要取正值），每行之和对应到参考地的电容。因此，在从矩阵到电路转换时，要区分矩阵层和电路层：矩阵层描述了电荷-电压关系，电路层则显示物理元件连接。切记不能将矩阵中的每一项直接等价为一个电容元件。

2.2 电感矩阵 (L) 与电阻矩阵 (R)

磁准静态求解得到的 电感矩阵 L 与 电阻矩阵 R 描述多端导体网络的电流路径阻抗。Q3D 通常对每个导体网络按源/汇配置求解各端口的感应电压响应：即对每对端口施加单位电流激励（各一源一汇），求解结构表面（或体）上的感应电流分布，从而计算阻抗矩阵 Z = R + jωL。具体来说，阻抗矩阵的虚部由磁能存储决定（确定 L 矩阵），实部由导体内的欧姆损耗决定（确定 R 矩阵），都具有对称和半正定的性质。这里的“源/汇”定义很重要：它们指定了电流注入和抽取的位置，从而决定了独立支路的电流方向和参考点。改变源/汇布置实际上改变了定义电流回路的方式，会影响 L 和 R 矩阵的具体数值。在 Q3D 中，求得的感抗矩阵对应于部分电感（partial inductance）意义下的多端网络，需要配合后续的矩阵规约来获得整个闭合回路的感抗值。

简言之，Q3D 的 L 和 R 矩阵来源于磁场仿真，不是简单的几何长度计算；矩阵中的对角元和非对角元分别表示自感和互感效应（或对应的串联感抗），而这些值依赖于所选的源/汇及当前网格解。在使用这些矩阵时，应始终结合源/汇配置来解释结果，而不是将其视作与几何无关的全局常量。

2.3 矩阵与电路的关系

从物理上说，Q3D 给出的 C、L、R、G 矩阵是节点/支路级别的综合结果，需通过 矩阵到电路的映射 转换为 SPICE 等效电路。官方推荐的等效模型常用 T 型网络 表示：左边和右边的支路分别对应引入侧和参考侧的系列阻抗（半系列结构各含 R、L），中间并联支路对应电容矩阵。多个源端时，电容支路会按照端口对称展开。需要注意的是，即便对于单个导体的自感 L_ii，其在电路中往往被分配到左右两个半网络中，与电路拓扑相关；不能简单地将电路中的单个电感元件等同于矩阵对角元。本质上，“矩阵层”与“电路层”有不同的分配方式，不可混淆。

3. 建模与边界条件

3.1 源 (Source) 与汇 (Sink) 定义

在 Q3D 中，每一个多端导体被视为由若干 源端 (source terminal) 与一个 汇端 (sink terminal) 组成。在静/准静态求解中，软件通过源端注入单位电压（或电流）激励，并在汇端抽取回流电流，将整个结构置于接地环境下。源/汇的选择影响每条独立支路的定义：源端之间流过的电流被视为正向，汇端作为共同参考点。通过不同的源/汇组合可以得到不同的阻抗矩阵分量。源/汇的抽象设计使得 Q3D 的矩阵求解可以在不改变物理模型的情况下，仅通过线性代数操作（如“Move Sink”）改变参考，从而得到新的矩阵结果。

对于用户建模而言，务必要明确每个导体的源/汇设置。例如，在 BGA 焊球模型中，一般将焊球顶端作为源端注入电流，焊球底端作为汇端抽取回流（或视焊盘、地平面为汇端）。如果需要将某个导体浮置（无限远参考）或把多个导体共地接地，可使用 “Float at Infinity” 或合并汇端等操作来指定。总之，源/汇决定了矩阵的物理含义：不同的约定会改变 L、R 矩阵的表示方式以及后续电路连接关系。始终把源/汇看作回路的引用点，而不是简单的标记，可以避免在解释结果时出现偏差。

3.2 边界条件与浮置

Q3D 默认采用无限远浮置边界：即所有导体都相对于“无穷远”（视为地）进行解算。这意味着导体矩阵的每一行之和为零（总电荷守恒），行和为零的情况对应于该导体相对于地的电容元件消失。在工程应用中，常需要把某些导体显式“接地”或“浮置”：

• 接地 (Ground)：通常将汇端设为接地端口，或在后续矩阵规约中把导体与地合并。
• 浮置 (Float)：使用“Float at Infinity”命令，可以将选定端口从矩阵中移除对地电容对应项，从而实现自由漂浮的近似。

在封装与 PCB 提取中，经常需要把大面积地平面或参考面作为回流参考；也有场景将整个局部网络浮置（比如只关心信号球之间的耦合）。浮置操作的结果会改变电容矩阵的对角和行和，但只影响等效电路层次的地电容，并不需要重新求解场。理解这一点对于后续的矩阵化简极为重要。

4. 矩阵规约 (Reduction) 与端口定义

4.1 矩阵规约的意义

Q3D 得到的是原始矩阵，对应每个几何导体端口的参数。在实际电路中，我们常需要对矩阵做规约（reduced matrix）来反映端口的合并、串并联连接或接地关系。矩阵规约可以在不重新求解场的前提下，通过线性变换实现等效的电路操作。例如，将同一网络内的两个源串联（Join in Series），对应对阻抗矩阵做等价的串联合并；将两个端口并联（Join in Parallel），对应对电容矩阵进行节点合并（行列和）。Q3D 中常用的规约操作包括：

• Series Join：按串联连接约束电流连续，对 Z = R + jωL 矩阵施加 KCL（基尔霍夫电流定律）条件，得到新的 Z。
• Parallel Join：按并联连接约束电压相等，对电容矩阵施加 KVL（基尔霍夫电压定律）条件，即对矩阵行/列进行合并相加。
• Add/Move Sink：增设或移动汇端，可用于修改回路参考而不改变量源分布。
• Float at Infinity：将端口浮空（行和置零），相当于移除对地电容。

官方说明特别指出，某些规约（如阻抗串联）必须在完整的 Z 矩阵上操作，而不是单独对 L 或 R 处理。这反映了规约的本质：它对应于电路连接关系在矩阵层面的投影。使用矩阵规约时应始终根据预期的电路连接来选择正确的操作顺序和矩阵形式。例如，对电感进行串联合并时，需要先求取阻抗矩阵 Z，再在 Z 上执行加法；而对电容进行并联合并时，可直接在电容矩阵上做行列和。

4.2 端口定义与矩阵输出

在 Q3D 中，每个 Design 会输出若干矩阵对象（典型地命名为 Cmatrix、Ymatrix、Rmatrix、Lmatrix 等）以及经规约后的矩阵（Reduced matrix）。端口的定义实际上决定了矩阵的行列顺序。例如，在包级 BGA 提取中，常将信号球和地球分开定义端口：信号球作为独立源端，所有地球球可并入一个参考端口。随后对信号/地端口做规约，得到针对信号端的 reduced 矩阵。正确的端口划分要结合电路功能：信号端口保留独立节点，地端口合并为统一的回流参考。端口定义和矩阵输出顺序需要与后续电路仿真保持一致，否则会出现“矩阵数值正确但对象不对应”的错误（常见于导出后连接 ADS/Matlab 时）。

典型工作流可概括为：

1. 建立几何模型。
2. 定义导体网络（nets）与端口。
3. 设置求解参数与网格划分。
4. 运行求解，获得原始 R/L/C/G 矩阵。
5. 进行端口重组与矩阵规约（Reduce）。
6. 导出等效电路或子网矩阵。
7. 进入时域仿真或指标分析（NEXT/FEXT）。

5. 网格与求解设置

5.1 网格划分与精度控制

Q3D 使用有限元/有限体积等技术对几何模型离散求解。网格划分的密度和质量直接影响提取结果的准确度和计算成本。一般原则是：重点在于电流分布和场强变化剧烈的区域。例如导体尖端、薄片和拐角处需要较密网格；而大平面区域可适当加粗网格。用户可在 Q3D 的 Adaptive Mesh 选项中设置最大迭代次数和误差公差，软件会逐步细化网格直至收敛。典型地，粗糙网格下的误差可能在 10% 以上，经过数次自适应后误差可降至 1% 左右。下表给出不同网格精度的经验比较：

细节建模对结果敏感时，应进行网格收敛性测试：对比不同细化级别下矩阵变化情况，确认关键参数稳定。如发现误差过大，可增加细化次数、调整局部网格尺寸或手动加细关键区域网格。

5.2 求解选项与收敛

在 Q3D 中，关键的求解设置包括：

• 求解频率：虽然 Q3D 本质上是准静态求解，但软件可输出频率相关的 R/L 矩阵。用户可指定一个或多个提取频率，分别获得对应的阻抗矩阵。通常选择直流和工作频率点进行提取，如关心高频损耗则可设高频点。
• 材料属性：设定导体的电导率和介质常数/损耗角。准确的材料参数对于 R、G 矩阵尤为重要。
• 收敛公差：包括电容收敛阈值和电阻收敛阈值。建议先用默认值，然后如需更高精度可调小阈值或增加迭代次数。
• 并行计算：对于大规模模型，可启用多核并行以加快速度，注意平衡核心数与内存资源。

错误设置可能导致无收敛或不合理结果。常见经验是：如对称结构提取结果应对称，不对称可能因为网格不均；电容矩阵理论上应为对称正定，检查是否因网格过粗出现负值或不对称。有时可尝试修改网格参数或分拆求解器任务以获得稳定解。

6. BGA 与 PCB 提取工作流（工程实践）

6.1 建模注意事项

对 BGA 焊球和多层 PCB 提取，模型通常包含球体、焊盘、过孔、铜柱、地平面等多个导体。关键事项有：

• 主流路径对齐：几何建模不必追求所有机械细节，但需保证电流主路径及耦合路径的尺寸精确。比如焊球半径、间距、连线宽度应准确建模；而倒角、圆角等可省略。
• 材料与层次：对于多层结构，需明确不同介质（如塑料、PCB 介质）的介电常数，并在模型中添加对应的分层结构。地平面可建模为无限薄导体。
• 端口规划：提前规划好哪些导体是信号网络、哪些是地网，合理命名几何对象。例如按行列编号和功能区分命名，有利于矩阵解读。
• 源/汇统一：对整个 BGA 阵列保持一致的源在上、汇在下（或反之）约定。在开始提取前，最好先约定好每个焊球的源/汇方向和顺序，以避免后续混乱。

6.2 提取与规约

推荐的提取步骤：

1. 原始提取：将全部焊球视作独立源端（底部共同汇接）并提取原始矩阵，检查导出的 C、L 矩阵（例如 spicecapMatrix 和 acindMatrix）。此时关注矩阵对称性、量纲和顺序是否符合预期。
2. 功能分类：根据实际功能（信号球 vs 地球）选择要保留的端口编号（keep indices），对矩阵进行规约。例如，将所有地球焊球串联并入一个地汇端，信号焊球端口分开保留；或者把多个信号球并联为一个相同网。通过 Join in Series/Parallel 等操作得到按功能划分后的 reduced 矩阵。
3. 等效网络构建：将规约后的矩阵导出，在电路仿真软件（如 ADS/SPICE）中搭建对应的等效电路网（T 型或提取本身生成的子电路）。进行时域串扰 (NEXT/FEXT) 计算，与预期进行对比验证。

这一流程可以按以下步骤理解：

面向 BGA 的提取与分析流程通常如下：

1. 完成几何模型与端口设置。
2. 在 Q3D 中提取原始矩阵。
3. 检查矩阵的顺序、对称性与量纲。
4. 按信号/地划分执行矩阵规约。
5. 导出规约矩阵或子电路。
6. 在 ADS/MATLAB 中构建等效网络。
7. 进行串扰时域仿真并验证结果。

6.3 验证与误差评估

工程上，应对 Q3D 提取结果进行验证：

• 自洽性检查：检查矩阵基本性质，如对称性、行和规则（浮置时行和为 0）、实对角元大于零等。
• 物理一致性：对比相邻几何参数变化的趋势。例如增大焊球半径应增加互感/互电容；关闭某个焊球（短路）应导致相应行列值变化。确认提取结果与直观物理规律一致。
• 结果收敛：通过网格收敛测试和多次提取验证稳定性。如发现 L 矩阵偏离预期，可尝试更细网格或检查源/汇分配。
• 与实验或理论值对比：对于简单结构可与解析计算（如平行板电容公式）比较；对于复杂结构可与电磁场仿真（全波或其他工具）结果相对比。

通过这些验证，可以估计提取误差并进行必要的模型改进。

7. 常见误区与注意事项

1. 混淆矩阵条目与电路元件：如前所述，Maxwell 矩阵条目不是电路电容电感的直接值。尤其是电容矩阵要通过拓扑转换才能得到实际电容元件（T 型拓扑）。误把矩阵项当成元件值会导致电路混乱。
2. 忽视源/汇重要性：在 L/R 提取中，源/汇确定了电流走向。将其视作可忽略的标签会导致对矩阵的误解。每次更改源/汇都应重新判断矩阵含义。
3. 只会导出矩阵，不会使用或验证：完整的工程能力体现在能够检查矩阵的顺序、量纲和对称性；知道规约会使哪些量消失；理解矩阵变化与物理结构之间的关系，并能将其输入电路仿真进行验证。仅仅运行软件并获得数字结果是不够的。

参考资料

• Ansys Q3D Extractor 官方用户指南 (Ansys Electronics Desktop Help)。
• Ruehli, A.E., “Equivalent circuit models for three-dimensional multiconductor systems” (Partial Element Equivalent Circuit, PEEC)。
• Paul, C.R., Introduction to Electromagnetic Compatibility. Chapters on field-to-circuit coupling.

注：以上来源概括了矩阵物理和耦合方法的原理，并在 Ansys 文档中有对应阐述。