SI对应PCB设计 — 叠层阻抗控制与布线拓扑的电磁场仿真

粗略来说，首先通过叠层设计控制阻抗。接着优化布线拓扑——禁止T型分支，基本采用菊花链。最小化过孔数以减少不连续点。最后作为串扰对策，确定信号间的间距规则。对于PCIe Gen5级别，推荐3W（线宽的3倍）以上的间距。

当然。所以才会用2D/3D电磁场仿真进行事前验证。例如在汽车ECU电路板中，CAN-FD、Ethernet、摄像头用的FPD-Link混在一起，串扰的影响用手算根本无法估算。不做仿真就试制的话，一次试制就会浪费数百万日元。

一个参考标准是“布线长度超过信号波长的1/10时，就要当作传输线来处理”。从上升时间 $t_r$ 求必要带宽的话是这样的：

例如PCIe Gen5的上升时间大约是15ps，所以带宽大约是23GHz。在FR-4电路板上的波长大约是8mm，布线长度超过0.8mm时，传输线效应就不能忽略了。在现代高速数字设计中，几乎所有的信号都是SI的对象。

PCB布线的SI分析有两个层面。首先是传输线理论（电报方程）：

这里 $R$, $L$, $G$, $C$ 是单位长度的电阻、电感、电导、电容（RLGC矩阵）。精确求出这些值就是2D截面场求解器的工作。

没错。特性阻抗是：

无损耗时可简化为 $Z_0 = \sqrt{L/C}$。信号的传播延迟是：

$l$ 是布线长度，$\varepsilon_{\text{eff}}$ 是有效介电常数，$c$ 是光速。FR-4的 $\varepsilon_r \approx 4.0$，所以 $\varepsilon_{\text{eff}}$ 在微带线中约为3.0~3.3，在带状线中约为4.0。

像过孔、连接器、封装这样的3D结构无法用电报方程描述。这时就需要求解全波的麦克斯韦方程组：

在实际工作中，通常构建混合模型：布线部分用RLGC模型，过孔和连接器部分用3D全波分析求得的S参数。

完全不是。叠层设计是同时决定SI、PI（电源完整性）和EMC三要素的最重要决策事项。具体来说：

• 信号层与参考平面的间距决定了特性阻抗
• 电源/GND平面对的间距决定了电源阻抗
• 叠层的对称性控制了电路板的翘曲

例如8层板中典型的SI优先叠层是：

根据目标阻抗反推。例如想要50$\Omega$的微带线，线宽设为5mil，在FR-4（$\varepsilon_r=4.0$）的情况下，用场求解器计算介质厚度约为3.5mil（约90$\mu$m）。不过电路板厂家的标准半固化片厚度是有限的，所以实际工作中是从可用的材料来调整线宽。

这里有个重要点——叠层的公差直接关系到阻抗。介质厚度有 $\pm$10%的偏差，阻抗就会有 $\pm$5%左右的变动。必须与电路板厂家进行规格确认。

这里是SI设计中最反直觉的地方。高频信号的返回电流，并不是回到最近的GND引脚，而是在信号布线下方的接地平面上，与信号平行地流动。这是因为电流会选择电感最小的路径。

这正是问题所在。参考平面上如果有开槽或过孔禁布区，返回电流就会绕行。这样会导致环路面积变大：

• 环路电感增大 → 阻抗不连续 → 产生反射
• 环路作为天线工作 → 成为EMI辐射源
• 与相邻信号的磁场耦合增大 → 串扰恶化

实际工作中常见的是，信号横穿电源平面的分割区域。这绝对要避免，不得已时，必须在分割处附近打缝合过孔来确保回流路径。

首先从2D截面场求解器开始。输入布线的截面形状（包括梯形轮廓）和介质的叠层结构，求解拉普拉斯方程或麦克斯韦方程组的2D版本：

这样可以高精度地提取RLGC参数。实际制造中，由于铜箔蚀刻，布线截面会变成梯形（上边窄），这与理想矩形截面的近似公式有很大偏差。2D求解器可以精确地建模这种梯形效应。

代入电报方程，生成传输线的SPICE模型（W-element、RLGC矩阵）。然后将其与IBIS-AMI模型结合，进行电路级的通道仿真。这是计算效率最高、布线长度在数十厘米以上时的标准方法。

像过孔、连接器、BGA封装、差分对的弯曲部分这样的3D结构，无法用电报方程描述。这时就要用3D全波求解器。主要方法有：

• FEM（有限元法）：以Ansys HFSS为代表。精确建模任意形状。自适应网格细化保证精度。
• FDTD（时域有限差分法）：CST MWS、Cadence Clarity 3D。一次计算获取宽频带特性。
• MoM（矩量法）：Cadence Sigrity、Keysight ADS Momentum。针对平面结构（PCB多层）优化。

实际工作中，求解一个过孔的S参数需要数分钟到数小时。所以对整个布线进行3D求解是不现实的，“3D部件用S参数提取 → 布线部分用2D模型 → 级联连接”这种混合方法是常规做法。

对于$N$端口网络，用频率函数描述入射波与反射波、透射波之间的关系。对于2端口（1输入1输出）的情况：

• $S_{11}$（回波损耗）：在输入端反射回来的比例。建议在$-20$dB以下。
• $S_{21}$（插入损耗）：传输到输出端的比例。表示通道的总损耗。
• $S_{12}$（反向传输）：对于被动的PCB布线，与$S_{21}$相等。
• $S_{22}$（输出端回波损耗）：输出端的阻抗匹配度。

对于差分信号，则是4端口，要看$S_{dd11}$（差分反射）、$S_{dd21}$（差分插入损耗）、$S_{cd21}$（模式转换=EMI指标）。

可以。将通道的全部S参数（Tx→封装→布线→过孔→连接器→Rx的全段级联）与IBIS-AMI模型的收发均衡器结合，就可以估计统计眼图和BER（误码率）。

在PCIe Gen5/6和DDR5中，通道的插入损耗预算由标准规定，例如PCIe Gen5的上限是-28dB@16GHz。超过这个值，均衡器就无法完全补偿了。

布线拓扑是“一个网络如何分叉连接到多个接收端”的构成。这是对SI设计影响最大的部分。

准确说是桩线谐振。信号在T型分支点到末端（=桩线长度$l_s$）之间往返。当桩线长度达到波长的1/4时会发生谐振，该频率处的插入损耗会急剧增加：

例如桩线长3mm时，大约在14GHz会出现陷波。PCIe Gen4的基频是8GHz，所以三次谐波会直接撞上。因此原则上禁止T型分支，要使用菊花链或Fly-by。

直接法（LU分解、Cholesky分解）または反復法（CG法、GMRES法）により連立方程式を解く。大規模問題では前処理付き反復法が効果的なんだ。

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

現在の所属: Ansys Inc.

Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。

現在の所属: Dassault Systèmes SIMULIA

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

現在の所属: COMSOL AB

異なるソルバー間でモデルを変換する際は、要素タイプの対応関係、材料モデルの互換性、荷重・境界条件の表現差異に注意が必要になるんだ。特に高次要素や特殊要素（コヒーシブ要素、ユーザー定義要素等）はソルバー間で直接変換できない場合が多い。

メッシュ収束性の確認、境界条件の妥当性検証、材料パラメータの感度分析がすごく大事なんだ。

• メッシュ依存性の検証: 少なくとも3水準のメッシュ密度で収束性を確認
• 境界条件の妥当性: 物理的に意味のある拘束条件の設定
• 結果の検証: 理論解、実験データ、既知ベンチマーク問題との比較

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

形状関数 $N_i$ を用いて未知量を近似:

これを数式で表すとこうなるよ。

これを数式で表すとこうなるよ。

連続体の支配方程式を離散化すると、以下の代数方程式系が得られる:

ここで $[K]$ は全体剛性マトリクス（または同等のシステムマトリクス）、$\{u\}$ は未知節点変数ベクトル、$\{F\}$ は外力ベクトルなんだ。

• 完全積分: 全ての項を正確に積分。剛性過大評価の傾向（ロッキング）
• 低減積分: 積分点数を削減。計算効率向上だが、アワーグラスモード発生のリスク
• 選択的低減積分 (B-bar法): 体積項と偏差項を分離して積分。ロッキング回避

• h-refinement: メッシュを細分化（要素サイズ h を小さく）して精度向上
• p-refinement: 要素の多項式次数を上げて精度向上
• hp-refinement: h と p を同時に最適化

収束速度: 二次要素で $O(h^2)$ のオーダーで誤差が減少（滑らかな解の場合）

PCB信号品質設計の実務的な解析フローと注意点を解説する。

1. 前処理 (Pre-processing)

• CADデータのインポートと形状簡略化
• 材料特性の定義
• メッシュ生成（要素タイプ・サイズの決定）
• 境界条件と荷重条件の設定

2. 求解 (Solving)

• ソルバー設定（解法、収束基準、出力制御）
• ジョブ投入と計算実行
• 収束モニタリング

3. 後処理 (Post-processing)

• 結果の可視化（変位、応力、その他の物理量）
• 結果の検証と妥当性確認
• レポート作成

• 応力集中部: 最低3層以上の要素を配置
• 応力勾配の大きい領域: 要素サイズを周囲の1/3〜1/5に
• 荷重印加点近傍: 局所細分化
• 遠方領域: 粗いメッシュで計算効率を確保

• 過拘束に注意: 剛体移動の拘束は6自由度のみ
• 対称条件の活用: 計算規模の削減
• 荷重の等価分配: 集中荷重 vs. 分布荷重の選択

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

現在の所属: Ansys Inc.

Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。

現在の所属: Dassault Systèmes SIMULIA

• メッシュ収束性を3水準以上で確認したか
• 力の釣り合い（反力合計）を検証したか
• 結果が物理的に妥当な範囲か確認したか
• 既知の理論解またはベンチマーク問題と比較したか

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

PCB信号品質設計に対応する主要な商用CAEツールの機能比較と、各製品の歴史的背景を詳述する。

Ansoft Corporationが開発した3D高周波電磁界シミュレータ。2008年にAnsysがAnsoftを買収。

現在の所属: Ansys Inc.

Computer Simulation Technology (ドイツ) が開発。2016年にDassault Systèmesが買収しSIMULIAに統合。

現在の所属: Dassault Systèmes SIMULIA

1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。

現在の所属: COMSOL AB

• 要素タイプの非互換: ソルバー固有要素は中立フォーマットで表現不可
• 材料モデルの差異: 同名でも内部実装が異なる場合がある
• 境界条件の再定義: 多くの場合、手動での再設定が必要
• 結果データの比較: 出力変数の定義（節点値 vs. 要素値、積分点値）に差異

PCB信号品質設計のツール選定においては以下を考慮:

• 解析規模: 数万〜数億DOFへのスケーラビリティ
• 物理モデル: 必要な構成則・要素タイプの対応状況
• ワークフロー: CADとの連携、自動化の容易さ
• コスト: 初期投資 + 年間保守 + 教育コスト
• サポート: 技術サポートの質とレスポンス

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。

PCB信号品質設計における最新の研究動向と先進的手法を見ていこう。

• 等幾何解析 (IGA): NURBS基底関数を直接使用し、CAD-CAE間のシームレスな連携を実現
• 粒子法 (SPH, MPM): メッシュフリー手法による大変形・破壊の追跡
• 位相場法 (Phase-Field): 界面の暗示的表現による複雑な界面追跡
• 機械学習支援: サロゲートモデル、物理インフォームドニューラルネットワーク (PINN)

症状: ソルバーが指定反復回数内に収束せず異常終了

考えられる原因:

• メッシュ品質の不足（過度に歪んだ要素）
• 材料パラメータの不適切な設定
• 不適切な初期条件
• 非線形性が強すぎる（荷重ステップの不足）

対策:

• メッシュ品質チェックを実施（アスペクト比、ヤコビアン）
• 材料パラメータの単位系を確認
• 荷重を複数ステップに分割（サブステップ数の増加）
• 収束判定基準の緩和（ただし精度に注意）

症状: 応力/変位/温度等が物理的に非現実的な値

考えられる原因:

• 境界条件の誤設定
• 単位系の混在（SI単位と工学単位の混同）
• 不適切な要素タイプの選択
• 応力特異点の存在

対策:

• 反力の合計を確認（力の釣り合い）
• 単位系の一貫性を確認
• 要素タイプの適切性を再検討
• 特異点除去またはサブモデリング

症状: 計算が想定時間の何倍もかかる

対策:

• メッシュの粗密分布の最適化
• 対称性の活用（1/2, 1/4モデル）
• ソルバー設定の最適化（反復法、前処理の選択）
• 並列計算の活用

症状: Out of Memory エラー

対策:

• アウトオブコア解法の使用
• メッシュ規模の削減
• 64bit版ソルバーの使用確認
• メモリ割り当ての増加

• FATAL 2012: 特異剛性マトリクス → 拘束条件の見直し
• USER WARNING 5291: 要素品質不良 → メッシュ修正
• SYSTEM FATAL 3008: メモリ不足 → MEM設定の調整

• Excessive distortion: 要素の過大変形 → NLGEOM確認、メッシュ改善
• Zero pivot: 拘束不足 → 境界条件追加
• Time increment too small: 収束失敗 → ステップ設定見直し