時間領域反射率測定(TDR)
理论与物理
TDR概述与原理
TDR是测量什么的装置?和电磁分析有什么关系?
简单来说,这是一种将阶跃波——即急剧上升的脉冲——入射到传输线中,观察反射波时间波形的测量方法。在有阻抗不连续的地方信号会发生反射,因此通过分析反射波的振幅和到达时间,就能知道“哪里”有“什么样”的不连续。
原来如此,和雷达原理相似呢。但这和CAE有什么关系呢?
这在高速基板的过孔和连接器设计验证中是必不可少的。例如PCIe Gen5/6和DDR5的布线,一个过孔的阻抗不连续就可能毁掉信号质量。通过3D FEM分析预测TDR波形,在试制前发现问题,这是最大的优点。
能在试制前发现问题确实很重要!具体在什么场景下使用呢?
列举一些典型用途:
- PCB布线特性阻抗验证:检查是否控制在50 Ω或100 Ω差动设计值的±5%以内
- 过孔不连续性评估:优化通孔过孔的残桩长度或反焊盘直径
- 连接器阻抗剖面:定位焊点处的容性不连续
- 线缆/线束断线位置定位:汽车或飞机大规模线束诊断
反射系数与传输线理论
能用公式解释一下反射发生的机制吗?
当负载阻抗为 $Z_L$,传输线特性阻抗为 $Z_0$ 时,反射系数 $\rho$(希腊字母“rho”)定义如下:
我们来整理一下 $\rho$ 值的含义:
- $\rho = 0$:完全匹配($Z_L = Z_0$)。无反射。理想状态
- $\rho = +1$:开路端($Z_L = \infty$)。与入射波同相全反射
- $\rho = -1$:短路端($Z_L = 0$)。与入射波反相全反射
- $0 < \rho < 1$:感性不连续。TDR波形呈阶跃式上升
- $-1 < \rho < 0$:容性不连续。TDR波形呈阶跃式下降
TDR波形向上跳变表示阻抗向高方向不连续,向下跳变表示向低方向,对吧。很直观!
阻抗剖面重构
能从反射系数求出阻抗值本身吗?
当然可以。从反射波形读取时刻 $t$ 的反射系数 $\rho(t)$,则该位置的阻抗可通过下式重构:
这里 $\rho(t)$ 是反射波 $V_r$ 与入射波 $V_i$ 之比,可通过实测求得:
例如,用 $Z_0 = 50\,\Omega$ 的同轴电缆进行TDR测量时,某位置 $\rho = 0.2$ 则:
$$Z = 50 \times \frac{1 + 0.2}{1 - 0.2} = 50 \times \frac{1.2}{0.8} = 75\,\Omega$$
即知道存在75 Ω的不连续。这是实际工作中汽车天线电缆常见的模式。50 Ω系和75 Ω系的连接错误一目了然。
时间轴对应空间位置呢。因为是往返时间,距离要减半吗?
没错。到不连续点的距离 $d$ 由传播速度 $v_p$ 和往返时间 $\Delta t$ 决定:
$c$ 是光速,$\varepsilon_{\text{eff}}$ 是有效介电常数。FR4基板的话 $\varepsilon_{\text{eff}} \approx 3.8$ 左右,传播速度约为光速的51%。1 ns的往返延迟对应约76 mm的距离。
电报方程与TDR的关系
我想再深入了解一下传输线的物理,支配方程是什么?
传输线由电报方程(Telegrapher's equations)描述。使用单位长度的电感 $L$、电容 $C$、电阻 $R$、电导 $G$:
无损耗时($R = G = 0$),该联立方程变为以速度 $v_p = 1/\sqrt{LC}$ 传播的波动方程。特性阻抗为:
TDR利用的就是这种波动在不连续点反射的现象。CAE仿真中,主流方法是直接用3D FEM/FDTD求解麦克斯韦方程得到S参数,然后从中计算TDR波形。电报方程是1D近似,但像过孔、连接器这样的3D结构无法用1D模型表达,因此需要全3D分析。
上升时间与空间分辨率
TDR的“分辨率”是怎么决定的?脉冲越尖锐似乎看得越细?
很敏锐的直觉。TDR的空间分辨率 $\Delta d$ 由阶跃脉冲的上升时间 $t_r$(10%〜90%)决定:
我们看个具体例子:
| 上升时间 $t_r$ | 空间分辨率(FR4) | 对应频率带宽 |
|---|---|---|
| 100 ps | 约 7.7 mm | 〜3.5 GHz |
| 35 ps | 约 2.7 mm | 〜10 GHz |
| 15 ps | 约 1.2 mm | 〜23 GHz |
PCIe Gen5(32 GT/s)的布线验证需要 $t_r \leq 20$ ps,要分辨过孔的残桩长度(0.5〜2 mm左右)就需要这个级别的分辨率。仿真中可以将上升时间作为设计参数自由设定,这是相对于实测的一大优点。
TDR核心方程总结
- 反射系数 $\rho = (Z_L - Z_0)/(Z_L + Z_0)$:表示阻抗不连续的大小和极性。$|\rho|$ 越大反射越强,正号表示感性(Z增加方向),负号表示容性(Z减小方向)的不连续。
- 阻抗重构 $Z(t) = Z_0(1+\rho)/(1-\rho)$:从TDR波形的反射振幅反算该位置的阻抗值。从而可以将传输线沿线的阻抗空间分布“可视化”。
- 距离转换 $d = v_p \Delta t / 2$:将往返延迟时间转换为单程距离。传播速度 $v_p$ 取决于介质的有效介电常数,因此基板材料的准确介电常数数据必不可少。
- 特性阻抗 $Z_0 = \sqrt{L/C}$:由传输线单位长度的电感与电容之比决定。通过截面形状(线宽、介质厚度、接地距离)来控制。
假设条件与适用范围
- TEM(横电磁波)近似:对于微带线等准TEM线,频率升高时精度下降,频率色散不可忽略
- 单次反射假设:仅在多重反射可忽略时,上述简单的阻抗重构才成立。存在多个邻近不连续时,多重反射的影响会导致精度下降
- 无损耗近似:对于有损传输线,脉冲在传播过程中会钝化,表观上升时间增大,分辨率下降
- 线性性:包含非线性器件(如ESD保护元件)时,叠加原理不成立,通常的TDR解释不适用
TDR的起源——雷达技术转用于传输线诊断的历史
TDR的原理与雷达完全相同:“脉冲发射 → 反射接收 → 根据到达时间计算距离”。1960年代,惠普(现是德科技)发布了首款商用TDR“HP 1415A”,首次能够以数值方式定位同轴电缆的断线/短路位置。最初主要用途是通信运营商的电缆维护,但随着1990年代GHz频段数字设计的普及,逐渐演变为PCB的SI(信号完整性)验证工具。如今,实测TDR设备的带宽已达到70 GHz以上,可以分辨100 μm以下的结构。
数值解法与实现
FDTD法进行TDR仿真
如何在仿真中再现TDR?因为是时域分析,所以FDTD很自然吗?
解读得很好。FDTD(有限差分时域法)是TDR仿真最直观的方法。直接在时域和空域对麦克斯韦方程进行差分近似:
使用FDTD进行TDR的要点:
- 激励源:在输入端口施加高斯脉冲或阶跃函数。上升时间要与实测条件匹配
- 电压/电流监视器:在输入端口分离并记录入射波和反射波
- 吸收边界条件:使用PML(完全匹配层)对计算区域边界进行无反射处理
- 时间步长:必须满足CFL条件 $\Delta t \leq \Delta x / (c\sqrt{3})$
CST Studio Suite是这种方式代表,其Transient Solver正是如此。
FDTD能直接输出时间波形,看来很容易进行TDR式的观察。
FEM与从S参数到TDR的转换
像Ansys HFSS这样的FEM求解器呢?是频域的吧?
FEM是在频域求得S参数后,通过逆傅里叶变换(IFFT)转换为TDR阻抗剖面。步骤如下:
- 在宽频带(DC〜数十GHz)求得S参数 $S_{11}(f)$
- 应用窗函数(如Kaiser-Bessel)抑制吉布斯振荡
- 通过IFFT得到时域反射响应 $s_{11}(t)$
- 由 $\rho(t) = s_{11}(t)$ 根据 $Z(t) = Z_0(1+\rho)/(1-\rho)$ 计算阻抗剖面
FEM方式的优点是,由于每个频点可以独立求解,因此可以使用自适应网格细化。以HFSS为例,它会在每个频率自动细化网格以保证S参数的收敛。实际工作中这种方式更主流。
是在频域求解再转换到时域呢。S参数的频率范围越宽,TDR的分辨率就越高吗?
没错。通过IFFT得到的TDR上升时间大致与带宽 $BW$ 成反比:$t_r \approx 0.35/BW$。所以如果求解DC〜20 GHz,就能得到约 $t_r \approx 17.5$ ps 的分辨率。不过,高频段S参数的精度直接关系到整个TDR波形的可靠性,因此必须特别注意网格和求解的收敛性。
网格策略
TDR分析有特别的网格注意事项吗?
TDR分析,特别是高速基板的过孔结构,有以下要点:
| 区域 | 网格尺寸参考 | 理由 |
|---|---|---|
| 过孔筒壁面 | 最高频率的 $\lambda/20$ 以下 | 精确捕捉圆柱面的电流分布 |
| 反焊盘间隙 | 间隙宽度的1/3以下 | 准确评估容性耦合 |
| 焊盘与焊盘连接部 | 导体厚度的2〜3倍网格层 | 再现趋肤效应和电流集中 |
| 介质内部(远处) | $\lambda/10$ 左右即可 | 控制计算成本 |
对于20 GHz分析,FR4($\varepsilon_r \approx 4.3$)时 $\lambda_{\min} \approx 7.2$ mm,因此过孔周边需要 $\leq 0.36$ mm 的网格。
端口设置与去嵌
听说端口设置错误会导致TDR波形异常?
这是实际工作中非常常见的错误。端口设置的注意事项:
- 波端口:能准确激励波导模式,但如果端口面不平坦会出错。微带线时端口面要包含足够的空气层
- 集总端口:简便,但高频时模式纯度下降。10 GHz以下实用
- 去嵌:扣除从端口到DUT(被测器件)之间的传输线段。不做这个处理,TDR波形会叠加多余的延迟和反射
例如,在HFSS中通过 Port → De-embed Distance 指定馈线长度。CST的时域求解器会自动进行端口位置校准,但必须确认参考平面偏移。
FDTD vs FEM——TDR分析中的选择
FDTD好比“用摄像机拍视频”——可以实时看到波的传播随时间如何推进。FEM则是“拍多张照片做成动画”的感觉——先求每个频率下的稳态,再用IFFT合成为动画。实际工作中,如果需要宽频带的TDR波形就选FDTD(CST),如果需要S参数精确收敛就选FEM(HFSS)。两者最终得到的TDR剖面是相同的。
関連トピック(Cross-topics)
- マイクロストリップ線路 — TDRで測定する主要な伝送路構造
- Sパラメータ — TDR変換の元データ。周波数領域の反射・透過特性
- シグナルインテグリティ(SI) — TDRが使われる設計検証の全体像
- 表皮効果 — 高周波での導体損失。TDR波形の鈍化に影響