1. 项目概述：为什么要在SPICE里搞定磁滞？

在电路仿真领域，SPICE（Simulation Program with Integrated Circuit Emphasis）是工程师的“数字实验室”，从一颗简单的电阻到复杂的射频芯片，几乎都能在它搭建的虚拟世界里先行验证。然而，这个强大的工具在面对一个看似基础却异常顽固的现象——磁滞（Hysteresis）时，往往会显得力不从心。磁滞现象广泛存在于磁性元件（如电感、变压器铁芯）和某些特殊材料（如铁电存储器、形状记忆合金）中，其核心特征是系统的当前状态不仅取决于当前输入，还严重依赖于历史状态。简单说，它“记得”自己走过的路。当你尝试用SPICE里标准的线性或非线性元件模型去拟合一个真实的磁芯电感时，仿真的波形和实测结果可能会大相径庭，损耗算不准，饱和点预测错误，导致整个电源或电机驱动设计在样机阶段翻车。

因此，“Modeling Hysteresis in SPICE”这个课题，远不止是建立一个数学曲线那么简单。它关乎仿真结果的保真度，关乎设计一次成功的可能性，是连接理想模型与物理现实的关键桥梁。无论是处理开关电源中磁性元件的铁损，还是模拟比较器、施密特触发器的回差特性，亦或是研究新兴存储器件的物理行为，一个精准的磁滞模型都至关重要。近年来，随着系统级仿真和硬件在环测试的普及，对模型精度要求越来越高，像level 1 spice这类基础模型库已无法满足需求，而wrenai mdl (modeling definition language)等新一代建模语言的兴起，也反映出业界对复杂物理现象（包括磁滞）进行标准化、高精度描述的迫切需求。

本文将从一个资深模拟电路设计者的角度，深入拆解在SPICE中构建磁滞模型的完整路径。我不会只停留在理论公式，而是聚焦于如何选择模型、如何用现有SPICE元件搭建、会遇到哪些坑，以及如何验证模型的正确性。目标是为读者提供一套从理论到实践、可直接复现的解决方案。

2. 磁滞现象核心原理与SPICE建模挑战

2.1 磁滞的本质：为什么它如此特殊？

要建模，必须先理解对象。磁滞，本质上是一种记忆效应和能量耗散现象。以最经典的铁磁材料磁化曲线为例（B-H曲线）：当你施加一个从零开始增大的磁场强度H时，磁感应强度B会沿着一条初始磁化曲线上升。当H减小到零时，B并不会回到零点，而是保留一个值，称为剩磁Br。需要施加一个反向的矫顽力Hc才能使B降为零。继续反向增大H，B会沿反向饱和，如此循环，形成一个闭合的回线。

这个回线揭示了两个关键点：

1. 多值性与路径依赖：对于同一个H值，B可能对应两个不同的值（上升支和下降支），具体是哪一个，取决于你是从哪个状态“走”过来的。这是与普通电阻、电容等无记忆元件最根本的区别。
2. 能量损耗：磁滞回线所包围的面积，代表了一个磁化周期内以热能形式耗散的能量（磁滞损耗）。回线越“胖”，损耗越大。

在电路中，这直接翻译为电感器的电压-电流关系或磁通-磁动势关系是非单值、有记忆的。SPICE的标准电感模型L，其特性是V = L * dI/dt，描述的是一个线性、无记忆的理想关系，完全无法刻画磁滞。

2.2 SPICE建模的经典思路与取舍

在SPICE框架内实现磁滞，通常有几条路径，各有优劣：

思路一：基于受控源的数学函数拟合这是最灵活也是最常见的方法。核心思想是使用SPICE中的行为建模元件（如伯克利SPICE3中的B源，或商业仿真器中的BV,BI源），通过电压或电流控制源，直接描述磁滞的输入-输出关系。例如，构建一个以电流I（对应H）为输入，电压V（对应B或磁通Φ）为输出的受控源，其控制函数是一个包含记忆状态的复杂数学表达式，如Preisach模型、Jiles-Atherton模型等的离散化实现。

• 优点：精度高，理论上可以拟合任意形状的磁滞回线。
• 挑战：实现复杂，需要在模型中内置状态变量（如前一个时间步的输入输出值）来实现“记忆”，这对SPICE求解器的收敛性是个考验。同时，计算量较大。

思路二：利用现成的非线性元件组合有些工程师会巧妙地利用SPICE库中已有的、具备滞回特性的元件来“拼凑”模型。最典型的例子是使用具有滞回特性的电压比较器（如施密特触发器）模型为核心，配合其他电路构建一个磁滞系统。或者，使用多个二极管、电阻、电容和受控源搭建一个能模拟磁滞回线形状的等效电路。

• 优点：有时能利用仿真器对标准模型更好的优化和收敛处理。
• 挑战：物理意义不直接，参数调整困难，且能模拟的回线形状和动态特性非常有限，通常只适用于特定场景（如数字电路中的滞回比较）。

思路三：调用外部模型或C代码高级的SPICE仿真器（如PSpice、LTspice、SIMetrix）支持通过外部子电路、MODEL语句定义更复杂的特性，甚至允许嵌入C语言或其它编程语言编写的模型（这就是wrenai mdl这类语言想标准化的领域）。你可以用C语言实现一个完整的Jiles-Atherton磁滞模型，然后编译成.DLL或.SO文件供SPICE调用。

• 优点：功能强大，运行效率可能更高，模型可复用性好。
• 挑战：门槛极高，需要跨领域的编程和编译知识，且模型移植性差，在不同仿真器间可能无法直接使用。

对于大多数电路设计者，思路一（行为建模）在灵活性、实现难度和精度之间取得了最佳平衡，是我们接下来重点探讨的方法。

注意：选择建模方法前，必须明确你的需求。你是需要精确计算铁芯损耗（要求回线面积准确）？还是只需要模拟数字电路的滞回开关特性（关心阈值点）？前者必须采用物理意义明确的磁滞模型（如Jiles-Atherton），后者可能一个简单的施密特触发器模型就足够了。盲目追求复杂模型只会增加不必要的仿真时间和收敛风险。

3. 基于Jiles-Atherton模型的行为级实现详解

在众多磁滞模型中，Jiles-Atherton (J-A) 模型因其物理意义相对清晰、参数较少（通常5个），且能较好地描述软磁材料的磁滞特性，在工程中应用广泛。我们将以此为例，展示如何在SPICE中用行为源实现它。

3.1 Jiles-Atherton模型快速回顾

J-A模型通过一组微分方程描述磁化过程。其核心思想是将总磁化强度M分解为可逆部分和不可逆部分。模型的关键方程如下：

1. 无磁滞磁化强度：M_an = M_s * (coth(He/a) - a/He)，其中He = H + alpha * M，这是朗之万函数，描述了没有磁滞时的理想磁化曲线。M_s是饱和磁化强度，a是形状参数，alpha是分子场系数。
2. 微分磁化率：dM/dH = (M_an - M) / (delta * k - alpha*(M_an - M)) + c * dM_an/dH。这个方程是核心，它决定了磁化强度M随磁场H变化的速率。其中，k是不可逆损耗系数，c是可逆系数，delta是方向符号（dH/dt > 0时为+1，否则为-1）。

我们的目标就是在SPICE中，构建一个以磁场H（或与之成正比的电流I）为输入，以磁化强度M或磁感应强度B（B = μ0*(H+M)）为输出的受控源，其行为由上述方程控制。

3.2 在LTspice中的分步实现

我们选择LTspice进行演示，因为它免费、强大且广泛使用。其行为电压源B源或行为电流源Bi源支持丰富的数学函数和if语句，适合实现J-A模型。

步骤1：定义模型参数和状态变量首先，我们需要在子电路或模型定义中，声明J-A模型的5个核心参数：Ms,a,alpha,k,c。同时，我们需要至少一个状态变量来存储“记忆”——即上一个时间点的磁化强度M。在SPICE中，我们可以用一个电容两端的电压来模拟状态变量，因为电容电压具有连续性（记忆性）。

.subckt JA_Hysteresis Pin Nin Pout Nout params: Ms=1.6e6 a=50 alpha=1e-3 k=10 c=0.1 * Pin, Nin: 输入端口，电压差代表H（或比例于H） * Pout, Nout: 输出端口，电压代表B（或M） * 定义内部节点用于状态存储 C_state M_state 0 1 IC=0 ; 电容C_state用于存储状态M，初始值为0

这里，输入电压V(Pin, Nin)被设定为与磁场强度H成正比。电容C_state上的电压V(M_state)代表磁化强度M。

步骤2：实现无磁滞磁化曲线M_an使用B源计算朗之万函数。注意，LTspice没有内置的coth函数，但可以用cosh(x)/sinh(x)或(exp(x)+exp(-x))/(exp(x)-exp(-x))实现。

B_Man Man 0 V=V(Pin,Nin)/1000 + {alpha}*V(M_state) ; 计算He，这里假设输入电压1V对应1A/m，需按实际比例调整 B_Man1 Man1 0 V={Ms} * ( (exp(V(Man)/{a})+exp(-V(Man)/{a}))/(exp(V(Man)/{a})-exp(-V(Man)/{a})) - {a}/V(Man) ) * 增加一个平滑处理，防止V(Man)接近零时的除零错误 B_Man_safe Man_safe 0 V=if(abs(V(Man))>1e-6, V(Man1), 0)

步骤3：实现微分方程并更新状态这是最关键的步骤。我们需要构建一个电路，使其动态满足J-A的微分方程dM/dH = ...。由于SPICE求解的是代数微分方程组，我们可以将微分方程转化为积分形式。 我们可以利用B源和电容的天然积分特性。让一个受控电流源Bi向状态电容C_state充电，而Bi的大小正是由dM/dH * dH/dt决定。由于dH/dt近似为输入电压的变化率，我们可以这样构建：

* 计算方向符号 delta B_delta delta 0 V=sgn(ddt(V(Pin,Nin))) ; ddt是LTspice中对时间求导的函数 * 计算dM_an/dH (近似) B_dMandan dMandan 0 V=ddt(V(Man_safe))/ddt(V(Pin,Nin)) * 防止分母为零，做一个保护 B_denom denom 0 V={delta}*{k} - {alpha}*(V(Man_safe)-V(M_state)) B_denom_safe denom_safe 0 V=if(abs(V(denom))>1e-3, V(denom), 1e-3*{delta}) * 计算dM/dH B_dMdH dMdH 0 V=(V(Man_safe)-V(M_state)) / V(denom_safe) + {c}*V(dMandan) * 将dM/dH 转化为流入电容的电流：I = C * dV/dt = dM/dt = (dM/dH) * (dH/dt) * dH/dt 用输入电压的导数近似 Bi_update M_state 0 I={dMdH} * ddt(V(Pin,Nin))

步骤4：生成输出最后，将状态变量M（即V(M_state)）与输入H结合，输出磁感应强度B。

B_Out Pout Nout V=4*3.1416e-7*( V(Pin,Nin)/1000 + V(M_state) ) ; B = μ0*(H + M)，单位换算需注意 .ends JA_Hysteresis

3.3 参数获取与模型验证

模型搭建好了，但参数Ms, a, alpha, k, c从哪来？通常有两个途径：

1. 材料手册：部分磁性材料供应商会提供J-A模型参数。
2. 曲线拟合：如果你有材料的实测B-H回线数据，可以使用参数优化工具（如LTspice的.step param结合.measure指令进行粗略扫描，或使用MATLAB、Python的优化算法）来反推这五个参数，使模型仿真曲线最接近实测曲线。

验证方法：

1. 静态（DC）扫描：在输入端施加一个从负最大值到正最大值缓慢变化的电压（模拟缓慢变化的H），观察输出B的波形。它应该画出一个漂亮的、光滑的磁滞回线。检查剩磁Br、矫顽力Hc是否与预期相符。
2. 动态（AC）测试：施加一个正弦波激励，观察不同频率和幅度下的B-H回线。回线面积应随频率增加而略微增大（如果模型包含了涡流损耗的简化等效，否则J-A模型本身是准静态的）。
3. 损耗计算：对一个周期的回线进行积分，计算面积，估算磁滞损耗，与材料手册给出的损耗参数进行对比。

实操心得：在实现微分方程时，最大的“坑”是除零保护和导数计算。ddt()函数在仿真中可能产生噪声，导致求解器不稳定。我的经验是，在分母表达式外加一个if(abs(denom)>small_value, denom, small_value*sign)的判断，这个small_value需要根据模型参数尺度仔细调整，太小不起作用，太大会引入误差。另外，初始条件IC的设置也很关键，错误的初始磁化状态可能导致仿真一开始就发散。

4. 针对不同应用场景的模型简化与变体

J-A模型虽然经典，但并非万能。在实际工程中，我们需要根据具体应用场景对模型进行简化或调整。

4.1 模拟施密特触发器：最简单的磁滞

对于数字电路或比较器中的滞回，我们不需要完整的磁滞回线，只需要两个不同的阈值电压。这可以用一个简单的正反馈电路来实现，SPICE中甚至有内置的施密特触发器逻辑门模型。但用行为源构建一个更通用的电压比较磁滞模型也很有教学意义：

.subckt Simple_Hyst in out params: Vhigh=2.5 Vlow=1.5 Vout_high=5 Vout_low=0 * 内部使用一个电压控制开关（或比较器）与正反馈 B1 fb 0 V=V(out) ; 反馈电压 B2 comp 0 V=if(V(in) > Vhigh, Vout_high, if(V(in) < Vlow, Vout_low, V(fb))) Rout out 0 1k ; 一个负载电阻 Bout out 0 V=V(comp) .ends Simple_Hyst

这个模型没有物理意义上的磁滞能量损耗，只有状态切换的延迟，适用于逻辑行为仿真。

4.2 功率磁芯损耗建模：关注能量损耗

在开关电源仿真中，我们有时更关心磁芯的总损耗（磁滞损耗+涡流损耗），而非精确的B-H回线形状。这时可以采用损耗电阻与理想电感并联的经典模型。但这个电阻是非线性的，其阻值随频率和磁通密度变化。我们可以用一个行为电阻来模拟：

.subckt Lossy_Inductor L=100uH Rdc=0.1 params: Kh=1e-5 Ke=1e-9 L1 1 2 {L} Rdc 1 2 {Rdc} * 行为电阻模拟损耗，其阻值是频率和电流的函数 B_Rac 2 3 V=V(1,2) ; 感知电感电压 B_Freq freq 0 V=abs(ddt(V(1,2))/(6.28*V(1,2)+1e-6)) ; 粗略估算频率 B_Rloss Rloss 0 V={Kh}*V(freq)*I(L1) + {Ke}*V(freq)*V(freq)*I(L1) ; 损耗电阻值，公式为示例 R_ac 3 4 R={V(Rloss)} .ends Lossy_Inductor

这里Kh和Ke是分别对应磁滞和涡流损耗的经验系数，需要通过实测损耗数据拟合得到。这种方法计算速度快，在系统级仿真中非常实用。

4.3 利用Level 1 SPICE模型与宏模型

对于一些非常早期的仿真需求，或者受限于仿真器功能，可能需要利用最基础的level 1 spice二极管、晶体管模型来搭建磁滞功能。例如，可以用两个背靠背的齐纳二极管配合电阻电容，构建一个具有滞回特性的压控模块。这种方法极其粗糙，精度很低，且参数与物理量之间没有直接关系，通常只用于概念验证或教学演示，不推荐用于严肃的工程设计。

5. 仿真调试、收敛性问题与实战技巧

即使模型方程正确，在SPICE中也可能无法成功仿真。以下是几个常见的“拦路虎”及解决之道。

5.1 收敛性失败与解决方案

磁滞模型因其强烈的非线性和记忆性，对SPICE求解器（特别是瞬态分析.tran）是巨大的挑战。

• 问题现象：仿真在开始或某个时间点突然停止，报错“Time step too small”、“Singular matrix”或“No convergence”。
• 排查与解决：
  1. 设置合理的初始条件：给状态电容（如我们例子中的C_state）一个合理的初始电压（IC参数），这相当于给磁芯一个初始磁化状态。不要让它从零开始一个剧烈变化的激励。
  2. 放宽仿真精度：在.tran指令中，尝试添加选项.tran 0 1m 0 1u startup uic。startup选项会使仿真器在开始时用非常小的时间步长缓慢“启动”，帮助克服初始非线性。uic表示使用用户定义的初始条件。
  3. 调整求解器选项：在仿真器设置中（如LTspice的Control Panel -> SPICE），可以尝试将求解器从默认的Modified trap改为Gear方法。Gear法对刚性方程（stiff equation，即变化速率差异巨大的方程）有时更稳定。也可以适当增大Trtol（相对误差容限，如从7改为8或9）或Abstol（绝对电流误差容限）、Vntol（绝对电压误差容限）。
  4. 添加串联小电阻或并联大电容：在模型的关键节点（如受控源输出端、状态节点）对地并联一个很大的电阻（如1GΩ）和一个很小的电容（如1pF）。这可以为求解器提供直流路径和高频衰减路径，有助于稳定数值计算。在电感或电流源支路串联一个极小电阻（如1μΩ）也能防止节点浮空。
  5. 平滑处理非线性函数：如前所述，对coth(x)、1/x这类在零点附近剧变的函数，一定要用if语句进行平滑处理，避免数学上的奇点。

5.2 仿真速度优化

行为模型，尤其是包含ddt()（求导）和复杂函数运算的模型，会显著拖慢仿真速度。

• 技巧：
  1. 简化模型：在满足精度要求的前提下，使用最简模型。对于系统级仿真，4.2节的损耗电阻模型可能比全功能的J-A模型快几个数量级。
  2. 减少ddt()的使用：求导运算开销大。考虑是否能用其他方式近似。例如，对于低频激励，dH/dt可以用(H_now - H_previous)/timestep在行为源内自己实现，但要注意引入延迟。
  3. 使用查表法：如果磁滞回线形状固定，可以先离线计算好一系列（H, B）数据点，然后在SPICE中使用B源的table或poly功能进行插值。这能极大提升速度，但失去了动态变化的灵活性。
  4. 合理设置最大时间步长：在.tran指令中设置一个合适的最大步长max timestep，防止求解器在平滑区域浪费计算资源。

5.3 模型验证清单

在将磁滞模型用于关键设计前，请完成以下检查：

1. 静态DC扫描：输入一个非常缓慢的三角波，输出是否形成封闭、光滑、对称的回线？
2. 对称性：正负半周的回线是否镜像对称？不对称可能意味着模型实现有误或初始条件不对。
3. 幅度缩放：改变输入幅度，回线是否按比例缩放？饱和区是否出现？
4. 频率响应：在合理频率范围内（如1Hz到100kHz），回线形状是否基本稳定？（纯J-A模型应基本稳定，若包含涡流等效，面积应随频率增加）。
5. 能量守恒检查：计算一个周期内输入的能量（∫V*I dt）和模型内部损耗（如果有的话），是否合理？
6. 与实测数据对比：这是黄金标准。将仿真得到的B-H回线与器件数据手册或实测曲线重叠，比较关键点（饱和点、剩磁、矫顽力）。

最后，记住一点，所有模型都是对现实的近似。在SPICE中建模磁滞，是在计算复杂度、仿真速度、模型精度和工程需求之间寻找最佳平衡点。从简单的滞回比较器到复杂的物理磁滞模型，选择最适合你当前任务的那一个，远比追求最复杂的模型更重要。我的经验是，对于大多数电源和电机驱动设计，一个经过合理参数拟合的J-A行为模型或一个等效损耗模型，已经能极大地提升仿真置信度，帮助你在设计阶段预见并解决很多潜在的磁性元件问题。