一、从一张仿真图说起
这是 COMSOL 算出来的动态BL(x)。看图的人通常会卡在：一条 BL(x) 为什么变成两条？这个和我们平时说的 BL(x) 是什么关系？
往深一步问，做大信号建模的人迟早会碰上：
- Le(x) 仿真值和 Klippel LSI 实测值为什么总对不上？
- Le 到底取哪个频率的值？1 kHz 还是 10 kHz？
- BL(x) 和 Le(x) 之间到底有没有耦合？
这其实是同一个问题：你把磁通匝链 ψ(x, i) 拆得太碎了。
二、第一性原理：从 ψ(x, i) 出发
音圈匝链的总磁通：
$$\psi(x, i) = \psi_{pm}(x) + L(x, i) \cdot i$$
第一项是永磁体在位置 x 处的匝链，与电流无关；第二项是电流自身产生的匝链,受位置（及饱和时的电流）调制。
用能量法（co-energy）对 x 求偏导，得到完整电磁力：
$$F_{em} = i \cdot \frac{\partial \psi_{pm}}{\partial x} + \frac{1}{2} i^2 \cdot \frac{\partial L}{\partial x}$$
定义：
$$BL_{static}(x) \equiv \frac{\partial \psi_{pm}(x)}{\partial x}$$
这是传统洛伦兹力意义下的 BL(x)。第二项 $\tfrac{1}{2}i^2 \cdot \partial L/\partial x$ 是磁阻力（reluctance force），只取决于 Le(x) 的形状，与永磁体无关。
下面所有问题都从这条公式里长出来。
三、动态 BL 为什么会"分叉"
工程上很多软件喜欢把总电磁力强行写成 BL·i 形式：
$$BL_{dyn}(x, i) = BL_{static}(x) + \frac{1}{2} i \cdot \frac{\partial L(x)}{\partial x}$$
回到开篇那张图：仿真扫描时电流取了正负两个值，BL_static 部分两次相同，但 磁阻力 项 $\tfrac{1}{2}i \cdot \partial L/\partial x$ 随电流反号——于是两条曲线在大位移处分得越来越开。
动态 BL 分叉的程度就是 $\partial L/\partial x$ 的代理指标——不用单独看 Le(x) 曲线就能定性判断磁阻力调制的严重程度。
橙色曲线在 x 方向的斜率即传统 BL(x)；绿色曲线在 i 方向的斜率即电感 L。曲面在 i 方向的轻微"扭曲"就是 L(x) 随位置变化的几何表征——也是 动态 BL 分叉的源头。
四、Le(x) 非线性的物理来源
Le(x) 随位置变化，本质是音圈周围磁阻路径的几何变化：音圈在磁隙中心时磁路里有铁极片（高 μ_r），外移时更多被空气包围（μ_r ≈ 1），自感下降——这就是 Le(x) 钟形曲线的来源。大电流下极片局部 B 超过饱和点，μ_r 从数千掉到接近 1，Le 骤降，得到 Le(x, i)。短路环、铜帽本质是被短路的二次绕组，反向磁通把 Le(x) 压低并平坦化。
五、Le(x, f)：不能回避的频率依赖
“Le 仿真和实测对不上"困扰过无数人，绝大多数案例的根源就一句话：大家都在偷偷做单频近似，但用的不是同一个单频。
铁极片是导电的。音圈交流磁通在极片里激发涡流，涡流反向屏蔽内部磁场，结果就是经典的趋肤效应：等效 Le 随频率下降，等效串联损耗 R 随频率上升。
不同电路模型对应不同精度：
| 模型 | 元件结构 | 适用范围 |
|---|---|---|
| R + L | Re + Le | 低频近似、工程粗算 |
| LR-2 (Thorborg) | Re + Le + (L₂ ∥ R₂) | 中频精度提升 |
| Wright | Re + K·(jω)^n | 分数阶，涡流物理更贴 |
| 多阶有理 | 高阶网络 | 大信号高保真 |
两个事实：
- 实测 Le——你拿到的 “Le = 0.5 mH” 是某个模型在某个频段拟合的等效值。换模型、换频带，数值就变。LSI 给的 Le(x) 是在特定小信号扰动频率下、按某个固定电路结构拟合出来的曲线。
- 仿真 Le——FEM 磁准静态求解给出的是给定频率下的复电感 $\hat L(f)$。要得到工程意义上的 Le，必须扫频后按相同的电路模型拟合。
严格写法：
$$L_e(x, f), \quad R_e(x, f)$$
仿真和实测都该展到这个维度。直接拿"COMSOL 算的 1 kHz 单点 L 值"对比"LSI 在工作频段拟合出的 Le(x)"，谈不上"差异”，只能叫"不可比"。
六、BL 与 Le 耦合的工程后果
把 BL(x) 和 Le(x) 当成两条独立曲线分别优化是工程上的偷懒。它们是同一个 ψ(x, i) 曲面的两个偏导。几个常见现象就自然了：
二次谐波失真：磁阻力项含 i²，必然产生 2 倍频分量。即使 BL_static(x) 关于 0 完美对称，只要 $\partial L/\partial x \neq 0$，二次谐波就跑不掉。
直流位移（DC offset）：大信号下 Le(x) 非对称调制产生平均力，导致音圈中性位置漂移。这是为什么平坦化 Le(x)（铜环、铜帽）能显著减少跳圈和压缩感。
电感反电动势的"位置项"：
$$u_L = \frac{d(L \cdot i)}{dt} = L \cdot \frac{di}{dt} + i \cdot \frac{\partial L}{\partial x} \cdot v$$
第二项是位置、速度、电流三者的耦合，常被简化丢掉，但在大冲程下是真实的能量交换通道。
七、给工程师的可执行结论
- 看到 动态BL(x) 分叉不要慌——那是磁阻力被强行折叠进 BL 的几何痕迹，分叉幅度就是 $\partial L/\partial x$ 的代理指标。
- 谈 Le 必须问清楚：哪个频率测的？什么电路模型拟合的？仿真和实测在比较之前先对齐这两个维度。
- 严格的 Le 仿真流程：每个 x 位置 FEM 扫频得 $\hat L(x, f)$，再按目标电路模型拟合到等效 Le(x)——和实测 LSI 的处理保持一致才可比。
- 短路环、铜帽、低 Bi 合金极芯这些设计，本质上都是在平坦化 $\partial L/\partial x$，从源头降低磁阻力带来的是真。