Simulink变压器饱和模型与励磁涌流仿真搭建指南

先说我为什么折腾这个模型。之前做变压器差动保护算法验证，我用Simulink搭了一个看似很标准的变压器模型，空载合闸跑出来，涌流压根看不到，电流就是一个小小的尖峰，跟教科书里画的完全不同。排查了半天，问题出在我用的励磁支路是固定电感，铁心永远不饱和，磁通再大电流也涨不上去，自然不会有励磁涌流。把这个支路换成带饱和特性的模型后，经典的偏置尖顶波一下子就出来了。

这个经历让我彻底明白一件事：研究励磁涌流、剩磁、磁滞这些问题，饱和模型不是可选项，是必选项。这篇内容我就把这套模型从物理原理到Simulink搭建方法，从参数换到波形调通，完整梳理一遍，最后附上我调试时踩过的坑和排查思路。

1. 为什么变压器饱和模型和励磁涌流模型难建

1.1 励磁涌流的本质：铁心饱和瞬间的电流爆发

变压器空载合闸为什么会产生涌流？核心机理就一句话：电压-磁通关系决定了磁通不能突变，而铁磁材料的饱和迫使电流在饱和段急剧上升。

展开说，变压器接上电压源后，铁心里的磁通必须满足法拉第电磁感应定律，也就是电压等于磁通对时间的导数，u(t)=dΦ/dt。如果合闸瞬间电压正好过零，那么稳态磁通的正弦分量是从负的最大值开始增长的，但从物理上讲，合闸前铁心里可能已经有剩磁，磁通不能从那个状态跳变到新的稳态，于是实际磁通等于稳态磁通再加上一个直流偏置分量。

这个偏置非常关键。假设稳态磁通幅值是Φm，剩磁是Φr，那么合闸后半个周波，磁通峰值可能达到2Φm+Φr，这已经远超变压器设计的饱和磁通Φsat了。一旦磁通越过饱和拐点，铁心的励磁电感会急剧下降，可能只有线性区电感的几十分之一甚至更低。这时候励磁电流不再按正常空载电流的百分数缓慢变化，而是被很小的电感限制不住，瞬间飙升到额定电流的好几倍甚至十几倍，波形呈尖顶状，偏在时间轴一侧，这就是励磁涌流。

很多初学者不理解，为什么线性模型看不到涌流。我把道理说直白一点：线性模型里的电感是常数，电流始终和磁通成正比。磁通再大，电流也只是线性增长，最多等于正常励磁电流乘一个系数，远远到不了“涌流”的量级。饱和模型的本质是让励磁电感在磁通超过拐点后“塌陷”，电感一小，电流才会爆发。所以做涌流仿真，第一条原则就是——励磁支路必须是非线性的。

1.2 磁滞、剩磁、饱和三者的耦合关系

铁心的磁化曲线不是一条单值曲线，而是一条带回线的闭合曲线。磁滞回线包围的面积，代表一个交变周期内铁心损耗的大小。表面上看这事跟涌流没什么关系，但其实影响很大。

磁滞的第一个作用是决定剩磁。当变压器正常运行时，铁心工作磁通在某个幅值附近交变。如果这时拉开断路器，电压降为零，磁通沿磁滞回线回退，但不会回到零，而是停在一个剩余磁通密度Br上，这个Br就是剩磁。剩磁的方向和大小，取决于断开瞬间电压的相位以及铁心材料特性。在实际工程里，通常剩磁标幺值在0.5到0.9之间，保守起见仿真中取0.6到0.8是很常见的。

磁滞的第二个作用是影响动态磁导率。磁滞回线的斜率决定了铁心在不同磁通水平下的动态电感，这会影响涌流的峰值和衰减速度。不过在Simulink的官方变压器模型里，饱和特性通常用单值磁化曲线近似，磁滞被忽略，剩磁用一个初始磁通标幺值来表示。这种近似对大多数继电保护研究完全够用，因为大家主要关心涌流峰值、二次谐波含量、衰减时间常数这些宏观指标，磁滞回线的具体形状影响不大。真正需要精细磁滞建模的场景，是研究铁心局部过饱和、损耗细致分布、或者某些特殊保护原理时才会用到Jiles-Atherton模型。

1.3 为什么值得把模型做扎实

励磁涌流不只是变压器自身的电磁暂态现象，它直接决定了差动保护的整定策略。差动保护需要区分内部故障电流和励磁涌流，常用的制动原理包括二次谐波制动、间断角闭锁、波形对称识别，这些算法都要在涌流环境下验证。如果模型里涌流波形不对，那么验证出来的保护逻辑就没有说服力。

此外，变压器饱和还会向系统注入谐波，在特定系统参数下可能引发谐波谐振或过电压。这些研究都需要一个能准确反映饱和非线性的变压器模型。所以一个可靠、可复现、参数方便的饱和模型，是整个仿真研究的地基，值得花时间做扎实。

2. Simulink建模方案选型：从磁路等效到可视化

2.1 三条可行的建模路线

在Simulink里建变压器饱和模型，我实际用下来有三条路线，适合不同需求。

路线A：使用Simscape Electrical（原SimPowerSystems）自带的变压器模型，比如Three-Phase Transformer或单相Transformer，勾选饱和选项后填入磁化曲线。这是最简单快捷的方法，适合大多数电力系统仿真、继电保护研究。官方模型在内部已经做好了磁路-电路耦合，你只需要把参数填对。

路线B：用Simscape Foundation library里的磁路域元件自己搭磁路模型。这种做法的自由度最高，能从几何尺寸、磁路截面积、气隙等物理层面建立模型，但这属于电机设计或电磁装置研发的范畴，做电力系统暂态研究一般不需要这么底层。

路线C：用S-Function或MATLAB Function自己实现J-A磁滞模型，把变压器当成状态方程嵌进Simulink。这条路适合研究磁滞细节，比如铁心损耗建模、剩磁精确计算、特殊励磁条件下的动态特性。

我的建议是，先走A，凡事先跑通再往细里做。如果课题方向明确要研究磁滞现象，再从A切到C。B通常不是电力系统研究方向的首选。

2.2 官方模型的适用范围和隐藏问题

官方变压器模型在勾选饱和选项后，需要填一条电流-磁通的曲线，默认按标幺值输入。模型内部把这条曲线转换为非线性电感，查表插值求电流。

但官方模型有一个隐藏问题：它的非线性电感不考虑磁滞损耗，所以在仿真中涌流衰减速度可能偏慢，因为缺少了磁滞损耗这个天然阻尼。如果只是定性研究涌流特征，问题不大；如果是定量计算涌流衰减到某个值的时间，就需要额外处理。常见做法是在励磁支路并联一个代表铁损的等效电阻，阻值按空载损耗数据折算。

另外要注意：官方模型的处理是基于查表插值，曲线拐点处的点如果分布不合理，会导致数值抖动。我一般会在拐点附近多加几个点，让曲线过渡平滑一些；但点太多又会明显拖慢仿真速度，这个平衡需要自己拿捏。

2.3 参数从哪里来：铭牌参数到磁化曲线的换算

实际拿到的变压器参数往往只有铭牌数据：额定容量、额定电压、空载损耗、空载电流百分数、短路损耗、阻抗电压。这些数据足以支撑磁化曲线的构建。

思路是这样的：空载电流百分数给出了额定电压下的励磁电流，这就是磁化曲线工作点的电流坐标；空载损耗可以用来折算铁损电阻；阻抗电压可以用来约束饱和区的斜率上限，防止曲线斜率过于极端导致数值发散。

从这两三个信息出发，构建一条“工作点+拐点+饱和区点”的磁化曲线，是完全够用的。我后面在3.4节给一个完整的计算示例，直接把每一步算给你看。

3. 核心实操：搭建变压器饱和与励磁涌流仿真模型

3.1 搭建基础仿真系统

打开Simulink，新建一个模型。这里我以单相变压器为例，三相变压器的搭法完全类似，只是把元件换成三相版本。

需要的元件有：

• 单相交流电压源（Simscape > Electrical > Specialized Power Systems > Fundamental Blocks > Electrical Sources > AC Voltage Source）
• 单相断路器（Breaker），通过外部信号控制合闸/跳闸
• 单相双绕组变压器（Transformer，Two Windings），在参数对话框中勾选Simulate saturation
• 电压测量和电流测量模块
• Scope示波器

搭建顺序没什么特别的，电源接断路器，断路器接变压器一次侧，二次侧可以开路（这正是空载合闸的条件）或者接一个很小的负载。注意变压器二次侧开路时，相当于空载运行，这是励磁涌流研究的标准工况。

合闸角的控制是涌流大小的第一变量。用Step信号控制断路器时，Step的发生时刻直接决定合闸瞬间的电压相位角。这里有个细节：为了让合闸角可预测，建议把交流电压源的初始相位设成0，然后通过Step时间精确计算合闸相位。50Hz系统一个周期是20ms，每毫秒对应18°电角度。要合在电压过零瞬间（0°），就在整周期时刻合闸；要合在电压峰值瞬间（90°），就在过零时刻的基础上加5ms。

我在操作中习惯把Step信号的跳变时间留成变量，比如用常量模块加一个Math运算，这样方便批量跑不同合闸角的工况。

3.2 剩磁初始化：两个办法

剩磁怎么设置，是很多人卡住的地方。我总结下来有两条路。

方法一：直接在变压器饱和参数设置里填初始磁通。官方变压器模型中有一个Initial flux参数，填标幺值就行，比如0.7代表70%额定磁通的剩磁。这是最简单直接的做法，但有一个坑：这个初始磁通是在仿真开始的时刻生效的。如果你合闸前变压器已经接在电网上，那么这个初始磁通和当时电压积分得到的磁通不匹配，仿真一启动就会产生虚假振荡。

方法二：先让变压器正常运行一段时间，在某个时刻跳闸，铁心里自然留下剩磁，然后再在另一个时刻合闸。这个方法最真实，也是我实际推荐的。实现方式就是两个断路器控制信号，第一个信号控制“正常运行-跳闸”的时序，第二个信号控制“跳闸-合闸”的时序。

我个人的实测体会是：研究剩磁影响时，优先用先跳闸后合闸的时序仿真，而不是手动填Initial flux。手动填剩磁最大的问题是你不知道剩磁方向和合闸电压相位的组合关系，经常糊里糊涂跑出来一个大涌流或小涌流，分析时还得回头核对参数设置。而用断路器时序，这个物理过程是自然发生的，剩磁方向和大小由跳闸瞬间决定，结果更有说服力。

3.3 磁滞特性建模：J-A模型的Simulink实现思路

如果要精细模拟磁滞，绕不开Jiles-Atherton模型。简单说，这个模型把铁磁材料的磁化强度M分解为可逆分量和不可逆分量，用5个参数描述：饱和磁化强度Ms、形状参数a、磁畴耦合系数α、钉扎系数k、可逆系数c。

核心方程是磁化强度对磁场强度的微分关系，H_e = H + αM，M_an = Ms·(coth(H_e/a) - a/H_e)，然后dM/dH按J-A方程计算。在Simulink里实现这个方程，通常用MATLAB Function或S-Function写成微分方程，输出B-H曲线，再通过B = μ0(H+M)和电磁感应关系连接到电路侧。

实际操作中，有几个地方要特别注意。第一，求解器要用变步长刚性求解器，推荐ode15s或ode23tb，因为J-A方程本身是刚性的，用普通RK算法很容易发散。第二，计算步长要足够小，我遇到过参数不收敛的情况，排查下来是k值和步长不匹配。第三，建议先用MATLAB脚本离线跑一遍J-A方程，把B-H曲线画出来看一眼，确认磁性回线形状正常，再放进Simulink联调，别一上来就做全系统联调，否则出了问题根本不知道是磁滞模型的参数问题还是电路拓扑的问题。

3.4 参数转换计算示例

我举一个具体的例子，把参数转换的每一步算清楚。

假设有一台变压器，额定容量1000kVA，电压10kV/0.4kV，空载电流百分数I0=1.5%，空载损耗P0=2kW，阻抗电压Uk=6%。

先算高压侧额定电流：Ie = 1000 / (√3 × 10) = 57.7A（这是三相变压器的算法，单相变压器直接用容量除以电压）。空载电流有效值就是0.015 × 57.7 = 0.866A，标幺值为0.015p.u.。

接下来构建磁化曲线的三个点（标幺值）：

• 第1点：工作点，磁通1.0p.u.，电流0.015p.u.，对应额定电压下的正常运行状态。
• 第2点：拐点附近，磁通1.1p.u.，电压高10%，励磁电流明显增大，保守可取0.04p.u.。这里不要照搬空载电流百分数，因为电压超过额定值10%后，变压器已经接近或进入饱和区，电流随电压非线性地快速上升。
• 第3点：饱和区，磁通1.2p.u.，电流可以取0.3到0.5p.u.。这个点电流给这么大，是因为铁心已经深度饱和，励磁电感剧烈下降，电流由饱和深度主导，不再是正常工作状态下的额定励磁电流概念。

这三个点加原点，就构成了一条可用于Simulink的电流-磁通曲线。填入官方变压器模型时，直接填标幺值坐标，Simulink内部会按照基准值转换回实际物理单位。

标幺化基准值怎么取？电压基准取相电压有效值，Ue = 10000/√3 = 5773V，角频率314rad/s，那么磁通基准Φb = Ue/(2πf) = 5773/314 ≈ 18.4Wb。电流基准就是额定电流57.7A。所以实际第1个饱和点：电流0.866A，磁通18.4Wb。仿真时如果模型里直接填的就是标幺值，就不会有单位换算问题。

这里还要提醒一句：变压器饱和磁通密度本身有一个物理范围，硅钢片一般在1.7T到2.0T，如果模型里的参数需要填特斯拉而不是标幺值，记得通过铁心有效截面积换算：Φ = B × A，A是铁心有效截面积，不要直接拿B值当磁通用。

4. 仿真结果解读：涌流波形怎么看

4.1 典型励磁涌流波形的特征

仿真跑通之后，Scope上会看到非常典型的励磁涌流波形：电流偏置在时间轴一侧，正半周有一个很高的尖顶波，峰值可能达到几十安培甚至几百安培，负半周几乎为零。这个形态背后对应的是铁心只在磁通超过饱和点的半个周期进入饱和，另一个半周处于线性区，励磁电流非常小。

波形里还有一个重要特征叫间断角，就是电流过零点附近有一段几乎为零的时间。这个间断角的大小与饱和深度直接相关，饱和越深，间断角越小。差动保护里的间断角闭锁原理，检测的就是这个特征。

建议把电压波形和磁通波形一起放到Scope里观察。磁通波形能够直观看到偏置加饱和的效果：稳态磁通在正常位置正弦变化，实际磁通被抬升一个直流偏移，进入饱和区的那段磁通被“削平”，这个削平现象就是饱和的标志。把磁通和电流对比看，很快能判断模型是否正常工作。

我在验证模型时还有一个习惯，就是先看磁通波形是否连续无跳变。磁通是电压的积分，在断路器动作瞬间可能会有一个小阶跃，但如果出现大幅跳变，那就是初始状态设置有问题，需要回头查Initial flux和断路器时序。

4.2 合闸角、剩磁与涌流幅值的关系

合闸角和剩磁是决定涌流幅值的两个关键变量。我做一个四组对照，结论特别清晰。

这个表格不是拍脑袋写的，是通过四组对照仿真得到的典型结果。物理上很好理解：合闸角为0°时，稳态磁通从负峰值开始变化，加上正向剩磁后，半个周波内磁通达到最高点，饱和最严重；反过来，如果剩磁方向与偏置方向相反，反而抵消了一部分磁通偏移，饱和程度降低，涌流就小。

在仿真里做这个对照实验时，可以用我前面说的先跳闸后合闸的方法，也可以通过Initial flux设置不同的剩磁初值。不管用哪种方法，这四组数据都会给你一个清晰的涌流幅值变化趋势，对理解涌流机理很有帮助。

4.3 涌流衰减与多机互联场景

励磁涌流的衰减不是一瞬间完成的。它由回路电阻和电感决定，时间常数τ=L/R。由于饱和后电感很小，涌流初期衰减较快；但磁通的直流偏置衰减慢，所以涌流又会在一个相对较长的时间尺度上逐渐减小。完整衰减可能需要几十个周波，尤其是大型变压器，电阻小，时间常数大，涌流持续时间更长。

在多台变压器并联的场景下，还会出现和应涌流现象。简单说，一台变压器合闸产生的涌流流过系统阻抗，会在另一台本已稳态运行的变压器中感应出涌流。这个现象只有多变压器互联的模型才能复现，我在做这类扩展仿真时，会把状态变量输出到工作区，方便后续在不同拓扑下复用模型。

5. 常见问题与排查技巧实录

5.1 仿真一开始就震荡：初始化问题

这是最常见的坑。表现是仿真启动后，电压电流波形剧烈振荡，有时候直接报错。原因基本可以锁定在初始磁通设置与合闸前变压器带电状态不匹配。

如果你手动设置了Initial flux，但合闸前电压源已经接在变压器上，仿真第一步就会因为磁通状态与电压积分不一致而产生虚假振荡。解决思路很简单：Initial flux只适用于“合闸前变压器不带电”的工况。如果合闸前变压器要带电运行，就改用先跳闸再合闸的时序，不要用Initial flux。

5.2 数值发散和代数环

官方非线性电感基于查表插值，如果磁化曲线不够平滑，特别在拐点附近，仿真会频繁出现步长缩小到接近零的情况，表现为仿真速度极慢或者直接提示数值发散。

对策有三条：

• 磁化曲线拐点区域多给几个点，让曲线过渡平滑。
• 求解器改用ode15s或ode23tb，这两个是刚性求解器，对这类非线性电感问题稳定很多。
• 如果提示存在代数环，检查一下电压测量信号通路。通常在测量电压后加一个一阶惯性环节或者受控电压源缓冲，能有效切断电流-电压-电流的代数依赖环。

5.3 参数理解错位导致结果异常

我见过不少人把“空载电流百分数1.5%”当成“容量百分数”，这完全不对。空载电流百分数是额定电流的百分数，而且严格来说是励磁电流。大型变压器空载电流可能只有0.5%到1%，小型变压器可能到4%到5%，不同变压器的数值差很多。填磁化曲线时，一定要按铭牌上的空载电流百分数来折算。

另外，饱和磁通密度和磁通的换算是另一个易错点。如果模型需要填特斯拉单位，记得通过铁心有效截面积把磁通密度换算成磁通量，不要直接拿B值当磁通用。

5.4 调试技巧：小步跑，逐步加复杂度

我踩过最大的坑就是一步到位，线性模型、饱和、剩磁、磁滞全放一起。一旦结果不对，根本看不出是哪一层的参数出了问题。

现在我调试这套模型的流程是：

1. 先用线性模型跑通电压、电流、磁通的基本关系，确认电路连接正确。
2. 加饱和特性，看涌流有没有出现，波形是否是半周尖顶波，有没有间断角。
3. 加剩磁设置，对比不同合闸角下的涌流大小，验证机理。
4. 最后再加磁滞模型或铁损电阻，修正涌流衰减速度。

每一步都有明确的验证点，出了问题能快速定位到刚加进去的模块，调试时间至少省一半。Scope里建议同时显示标幺值和实际值，分析波形时心里更有数。我自己会在模型里加一组输出下标幺值的计算模块，这样无论参数怎么改，波形对比都不会被实际电气量的量级差异干扰。

最后再说一个扩展方向。这套饱和模型不只用来观察涌流波形，更实用的玩法是把差动保护的逻辑模块接到二次侧电流互感器后面，直接做保护算法的闭环验证。很多继电保护方向的课题和实际工程问题，就是从这样一个变压器饱和模型开始的，祝你也能调出漂亮的涌流波形。