基于TDR分析的PCB信号完整性评估

简单来说，TDR就是“电路板的雷达”。它向PCB走线发送高速阶跃信号，并在时间轴上观测反射信号从哪里、以多大强度返回。观察反射波形，就能看到过孔的阻抗不连续、连接器的寄生电容、线宽变化、参考平面（地平面）的缺失——所有这些都会以“台阶”的形式出现在波形中。

问得好。VNA的S参数是频域数据，看的是“整个传输通道”的特性。而TDR是时域的，所以能直观地看出“问题出在哪个位置”。例如，在连接器后方30mm处的过孔阻抗跳变，在TDR波形上一目了然。实际工作中两者会区分使用，但要定位空间上哪里有问题，TDR绝对更方便。

TDR设备产生阶跃函数（上升时间约20~35ps）。这个阶跃信号沿着传输线前进，在阻抗变化的地方，一部分会反射回来。反射系数 $\rho$ 由下式表示：

这里 $Z_0$ 是系统的基准阻抗（通常50 $\Omega$），$Z_L$ 是不连续点处的阻抗。如果 $\rho > 0$，则阻抗升高（感性不连续）；如果 $\rho < 0$，则阻抗降低（容性不连续）。

没错。从TDR波形读取时刻 $t$ 的反射系数 $\rho(t)$，则该位置处的阻抗为：

将其沿时间轴绘制的图就是阻抗剖面。SI工程师通过观察这个剖面来判断“这里偏离了50 $\Omega$”、“过孔处跳变到55 $\Omega$”、“连接器处下降到45 $\Omega$”。对于DDR5内存走线设计，通常的标准是控制在 $\pm 5\%$（47.5~52.5 $\Omega$）以内。

时间和位置的转换很简单。因为信号在传输线上往返，所以时间差 $\Delta t$ 对应的空间距离为：

这里 $v_p$ 是传输线中的信号传播速度。对于FR-4基板（$\varepsilon_r \approx 4.0$），$v_p \approx 1.5 \times 10^8$ m/s（约为光速的一半）。TDR设备的上升时间 $t_r$ 决定了空间分辨率：

例如，$t_r = 20$ ps 时，$\Delta x_{\min} \approx 1.5$ mm。也就是说，1.5mm以下的短不连续无法分辨。对于DDR5或PCIe Gen5/6级别的高速设计，就必须使用上升时间在20ps以下的TDR设备。

现场常见的有这4种：

正是如此。$\rho > 0$（波形向上跳变）意味着阻抗升高=感性不连续。$\rho < 0$（向下凹陷）意味着阻抗降低=容性不连续。把这个记在脑子里，解读TDR波形的速度就会大幅提升。

当然可以。将实测TDR与仿真TDR重叠，验证FEM模型的合理性——这是SI工程师的基本工作流程。获取仿真TDR的方法主要有两种：

• 直接法（时域）：使用FDTD或FIT输入阶跃信号，直接观测反射波。CST Studio Suite擅长此道
• 间接法（频域→时域转换）：用FEM通过频率扫描计算S参数（$S_{11}$），再用逆FFT（IFFT）转换为TDR波形。HFSS/SIwave的标准方法

实际工作中间接法更常用。原因是频域分析的网格要求较低，且可以利用现有的S参数数据。

如果有 $S_{11}(f)$ 数据，按以下步骤转换：

1. 将 $S_{11}(f)$ 乘以阶跃函数的频谱 $V_{\text{step}}(f) = 1/(j2\pi f)$
2. 应用逆FFT（IFFT）得到时域的反射波形 $v_{\text{refl}}(t)$
3. 将反射波加到入射阶跃上构成TDR电压波形：$V_{\text{TDR}}(t) = V_0 + v_{\text{refl}}(t)$
4. 转换为阻抗：$Z(t) = Z_0 \cdot V_{\text{TDR}}(t) / (2V_0 - V_{\text{TDR}}(t))$

频率数据的上限决定了带宽，下限决定了TDR波形的时间长度。例如，$f_{\max} = 20$ GHz 时，有效上升时间约为 $0.35/f_{\max} \approx 17.5$ ps。

现场工作中区分使用很重要。传输线建模有三个层次：

实际工作中，主流是混合方法：“走线部分用2D截面分析提取RLGC→仅对过孔和连接器进行3D全波分析→在电路仿真器中连接整体并合成TDR波形”。全部用3D分析往往会浪费计算资源。

全波分析（FEM/FDTD/MoM）直接求解麦克斯韦方程组，理论上最精确。但结果质量依赖于网格质量、端口设置、材料参数精度。另一方面，电路模型提取是将全波结果简化为RLC等效电路。优点是：

• 可以在SPICE类电路仿真器（HSpice, ADS）中快速进行瞬态分析
• 可以与驱动器/接收器的IBIS模型结合，进行系统级的TDR/眼图评估
• 参数化研究（改变走线长度、调整终端电阻等）速度极快

需要注意的是，电路模型“仅在原始全波分析的频率范围内有效”。将基于20 GHz分析创建的等效电路用于56 Gbaud NRZ设计是危险的。

TDR测量需要以下物品：

理想情况下，校准到DUT的入口，即SMA连接器的尖端。这称为“连接器尖端校准”。如果校准不充分，连接器本身的阻抗波动就会与DUT的特性混在一起。量产评估中在基板上设置SMA焊盘的原因正在于此——用手持探针接触的方法，地线电感会在波形中产生额外的台阶，导致无法判断是原本的阻抗波动还是探头的影响。

“重叠”只是第一步。相关性验证（模型与实测的关联验证）步骤如下：

1. 获取实测TDR波形——使用校准后的测量系统记录阻抗剖面
2. 获取仿真TDR波形——从同一DUT的3D模型计算S参数，并通过IFFT转换为TDR
3. 对齐时间轴——对齐测量面（参考平面）的位置。调整延迟使电长度一致
4. 重叠波形并评估差异——确认阻抗绝对值差是否在 $\pm 2.5\,\Omega$（50 $\Omega$ 系统下 $\pm 5\%$）以内
5. 调整不一致处的参数——主要的调整参数有以下3个：

是的，通过相关性验证的模型可以作为“已验证模型”对待。对于下一代产品的设计变更——例如增加过孔数量或改变走线层——就可以仅通过仿真进行评估。无需每次都制作原型进行实测，从而大幅缩短设计周期。这就是V&V（验证与确认）在SI中的实践形式。

时间门控是一种“在TDR波形上只截取特定时间窗口”的技术。例如，在测量连接器→走线→过孔→走线→连接器这样的路径时，有时希望去除连接器的反射，只看走线和过孔的特性。设置时间窗口（门）并屏蔽连接器部分，就可以只提取DUT（被测器件）本体的阻抗。

类似但严格来说不同。

直接法（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)$ のオーダーで誤差が減少（滑らかな解の場合）

TDR解析（SI応用）の実務的な解析フローと注意点を解説する。

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水準以上で確認したか
• 力の釣り合い（反力合計）を検証したか
• 結果が物理的に妥当な範囲か確認したか
• 既知の理論解またはベンチマーク問題と比較したか

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

TDR解析（SI応用）に対応する主要な商用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. 要素値、積分点値）に差異

TDR解析（SI応用）のツール選定においては以下を考慮:

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

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

TDR解析（SI応用）における最新の研究動向と先進的手法を見ていこう。

• 等幾何解析 (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: 収束失敗 → ステップ設定見直し