三相变压器建模与仿真研究设计

三相变压器的建模与仿真（Matlab）

摘 要

研究三相变压器地暂态过程，建立一个完善的变压器仿真模型，对变压器保护方案的设计具有非常重要地意义。

本文在Matlab的编程环境下，分析了当前的变压器仿真的方法。在单相情况下，分析了在饱和和不饱和的励磁涌流现象，和单相励磁涌流的特征。在三相情况下，在用分段拟和加曲线压缩法的基础上，分别用两条修正的反正切函数，和两条修正的反正切函数加上两段模拟饱和情况的直线两种方法建立了Yd11、Ynd11、Yny0和Yy0四种最常用接线方式下三相变压器的数学仿真模型，并在Matlab下仿真实现。通过对三相励磁涌流和磁滞回环波形分析，三相励磁涌流的特征分析，总结出影响三相变压器励磁涌流地主要因素。最后，分析了两种方法的优劣，建立比较完善的变压器仿真模型。

建立比较完善的变压器仿真模型，利用MATLAB对变压器的励磁涌流、内部故障和外部故障进行仿真。

1．2 本论文的主要工作

针对前面所讨论的三相变压器建模问题，本论文进行了重点、深入的研究，进行了理论分析和仿真计算，并提出了相对较好地变压器仿真模型。本论文主要包括一下几个方面：

（1）简要分析了当前的变压器仿真的方法，比较了相互之间的优缺点。

（2）在单相情况下，分析了在饱和和不饱和的励磁涌流现象，分析了单相励磁涌流的特征。单相程序在Matlab上仿真实现。

（3）在三相情况下，在用分段拟和加曲线压缩法的基础上，分别用两条修正的反正切函数，和两条修正的反正切函数加上两段模拟饱和情况的直线两种方法建立了Yd11、Ynd11、Yny0和Yy0四种最常用接线方式下三相变压器的数学仿真模型，并在Matlab下仿真实现。通过对三相励磁涌流和磁滞回环波形分析，三相励磁涌流的特征分析，总结出影响三相变压器励磁涌流地主要因素。最后，分析了两种方法的优劣，建立比较完善的变压器仿真模型。

三相变压器建模及仿真研究 第二章 变压器的基本原理

2.1变压器的工作原理

变压器是一种静止的电器，用于将一种形式的交流电能改变成另一种形式的交流电能，其形式的改变是多种多样的。既可以改变电压、电流；也可以改变等效阻抗或电源相数、频率等。

i1i2 u1 e1?1 ?2 e1u2

图1－1 变压器示意图

以单相为例，研究变压器的基本工作原理。图1－1所示为一台变压器的示意图。它由铁芯和线圈组成。接电源的原边线圈成为初级线圈；接负载的副边线圈称次级线圈。设原、副边线圈匝数分别为

?1、

?2。根据电磁感应现象，电能可从原边输送到副边，

但原、副边具有不同的电压和电流。

变压器内部的磁场分布的情况是非常复杂的，但是我们总可以把它们折算为等效的两部分磁通。其中一部分磁通?沿铁芯闭合，同时与原、副绕组相链，是变压器能量变换和传递的主要因素，称为主磁通或互感磁通；另一部分磁通?1?主要是通过非磁性介质（空气或油），它仅与原绕组全部相链（只与原绕组部分匝数相链的露刺痛已折算为全部原绕组相链而数值减少的等效磁通），故称它为原绕组的漏磁通。

根据电磁感应定律，当磁通?和?1?随时间变化时，分别在它们所交链的绕组内感应电动势：

1

三相变压器建模及仿真研究 d??e1???1?dt?d?? ?e2???2 （2.1）

dt?d?1??e???1?1?dt?式中e1、e2是主磁通在原、副绕组所感应的电动势瞬时值；e1?是原绕组漏磁通在原边感应的电动势瞬时值。

所以，

e1e2??1?2，设变压器的变比为k??1?2，则U1?kU2,I1?I2k。所以利用变

压器可以在传输电能的同时改变其电压和电流。

2.2 三相变压器的等效电路及连接组问题

现在电力系统都采用三相制，所以实际上使用得最广泛的是三相变压。从运行原理来看，三相变压器在对称负载下运行时，各相的电压、电流大小相等，相位彼此互差120，故可任取一相分析，即三相问题可简化为单相问题。

根据变压器原、副绕组电动势的相位关系，把变压器绕组的连接分成各种不同组号称为绕组的连接组。在不同的连接组下，三相变压器的等效电路略有不同。现以Yd11连接组为例，做三相等效电路等效电路图如图1－5所示。

uaLsrs?ubLsrsucLsrs

图1－5 Yd11连接组三相等效电路

在三相变压器中，用大写字母A、B、C表示高压绕组的手段，用X、Y、Z表示高压

2

三相变压器建模及仿真研究 绕组的末端；低压绕组首、末端则应用对应的小写字母a、b、c和x、y、z表示。星形连接的中点用字母O表示。

不论原绕组或副绕组，我国主要采用星形和三角形两种连接方式。

为了形象地表示原、副边电动势相位地关系，采用所谓的时钟表示法：即把高压绕组的电动势向量作为时钟的长针并指向12，低压绕组的电动势相量作为时钟的短针，其所指数字作为单相变压器连接组的组好。

在我国生产的变压器中，以Yd11、Ynd11、Yny0、Yy0（n表示中性点接地）四种连接组为主。

第三章 变压器仿真的方法

从20世纪60年代开始，人们就花费大量的精力去解决变压器的计算机模型问题。由于变压器的非线性特性，这被证明是困难的课题。不像线性系统一样，没有一般的解决方案可以解决非线性方程。即便是数字式的解决方案，也只能很困难的解决某一类的非线性方程，在稳定的状态下，存在好的变压器模型。然而，在瞬变的状态下，还没有完全令人满意的变压器模型。

变压器的性能主要取决于其铁心的磁化特性，即铁心的磁滞回环，因此对铁心磁滞回环的拟合是最基础、最重要的工作。在变压器特性的数值仿真计算中，对磁滞回环的拟合提出了以下几个要求：①具有较高的精度；②在大范围内不分段，具有光滑性，否则会引起变压器特性仿真计算过程的不稳定；②具有稠密性，因为通过实验只能得到有限条磁化曲线，而变压器仿真中需要知道B—H平面中的任意一条曲线。由于铁心的饱和特性、磁滞现象等非线性因素的影响，很难用数学模型精确地描述铁心的动态磁化过程。

变压器通过铁心磁场作用建立一次侧和二次侧的电磁联系。因此变压器暂态建模的关键是对铁心动态磁化过程的数学描述。根据对磁化特性曲线描述的不同，现有研究用的变压器模型大致有下列4种：

（1）模型A——基于基本励磁曲线的静态模型； （2）模型B――基于暂态励磁特性曲线的动态模型；

（3）模型C――基于暂态励磁特性曲线的非线性时域等效电路模型； （4）模型D――基于ANN的变斜率BP算法创建的模型。

3

三相变压器建模及仿真研究 3.1 基于基本励磁曲线的静态模型

基于基本励磁曲线的变压器模型只考虑饱和引起的非线性，即采用如图2－1所示的基本磁化曲线作为变压器暂态工作特性曲线进行二次侧电流的计算。等效电路图如图2－2所示。

图2－1 基本励磁曲线

图2－2 静态模型

由磁通守恒和KCL定律可以得到以下基本方程组：

?N1(i1?i)0?Ni22?did???R2i2?L22 （2.1） ?N2dtdt???f(i0)??式中i1为一次侧电流；i0为励磁电流；i2为二次侧电流；?为主磁通；N1、N2为一、二次侧匝数；R2、L2为二次侧负载。

由方程组（2－1）中的第一和第三个方程得到i2?N1N2(i1?i0)，将

d?dt?d?di0?di0dt代入

4

三相变压器建模及仿真研究 方程组第二方程，整理可得：

di0?di???2d???N1R2i1?N1L21?N1R2i0?/?N2?N1L2? （2.2） dtdtdt0????HlN1di0dt因?＝Bs和i0?

，故有

d?di0?sN1dBldHdi1dt，代入式（2.2）可得：

2?(N1R2i1?N1L2?N1R2i0)/(N2sN1dH?N1L2)

（2.3）

用四阶龙格－库塔法或隐式梯形公式就可以求解一阶常微分方程式（2－2）或式（2－3），从而建立了变压器仿真数学模型。

3.2基于暂态磁化特性曲线的动态模型

这类变压器模型建立在对动态磁化特性曲线的数学描述之上。

暂态磁化特性曲线??f(i0)的描述，最常用的是采用极限回环压缩法。即假定铁心磁化曲线的主磁滞回环和次磁滞回环具有相似性，由主磁滞回环压缩生成次磁滞回环。例如用反正切函数拟合主磁滞回环，其表达式为：

hni?C(??)i ?(i0)??arcta （2.4） 00式中?、h、?和C为常数。

在上升轨迹???d???0??dt?和下降轨迹?????d??(0)(0)?0?的转折点(i0,?)?dt?将主磁滞回环按压缩

??)/?压缩生成系数K??b?(0)?(?i0(0)?)?/?arctanh(i0(0)?C)??向直线??(?i0?22??2??次级回环的下降支或上升支。

图2－3所示为动态磁化特性曲线，其中f1(i0)，f2(i0)为极限磁滞回环，(i0(0)，?(0))为转折点(假设?从上升变成下降)，则f3(i0)为经过该点的次级回环下降支。在f1(i0)，

??f2(i0)形成的回环内的f3(i0)部分为经过该转折点的暂态磁化轨迹。

5

三相变压器建模及仿真研究

图2－3 局部磁滞回环轨迹

3.3非线性时域等效电路模型

该模型用几个电路元件分别模拟造成变压器非线性的因素。

因为引起变压器非线性的主要因素有饱和、涡流和磁滞，所以用三个电路元件模拟这些因素，并将各元件流过的电流线性叠加，得到励磁电流。其表达式为：

i0?im?ih?ie (2.5)

式中 im 为磁化电流；ih为磁滞电流；ie为涡流电流。

因为剔除了其它影响因素而单独进行考虑，故im可以用无磁滞曲线(基本磁化曲线)来表示，这是一个仅仅与磁链有关的表达式。其表达式可以表示为

(??) (2.6) im?f(?)?g? 磁滞是由交变电流产生，其大小和电压以及频率有关。但实验表明，在50Hz到400Hz内，磁滞随频率的变化而改变得很小，故频率的影响一般用一个常数表示。磁滞电流部分的表达式为：

??1 ih?khV??d??kh??dt????1 (2.7)

6

三相变压器建模及仿真研究 其中?为斯坦梅茨(Steinmetz)系数，由铁磁材料的特性决定。的一个常数。

kh设定为在50Hz下

涡流电流和磁通、磁通变化率以及频率有关。但是在电流频率不超过400Hz的情况下，涡流电流不会因频率改变而显著变化。因此可以不考虑频率变化对涡流的影响，表达式如下：

ie?re?1V?re?1?d? （2.8） dt 综合以上各式可得到考虑了饱和、磁滞和涡流影响的变压器励磁电流暂态数学模型，其表达式为：

i0?g(?)???(Kh(d???2d? )?re?1)dtdt令 Lm?g(?),Rm?1?Gm?kh(d???2)?re?1 dtd? 则有 i0?Lm?vdt?Rm?1? （2.9）

dt其等效电路如图2－4所示。

图2－4 非线性时域等效电路模型

7

三相变压器建模及仿真研究

第四章 三相变压器的仿真

电力系统中的变压器通常是三相的，而三相变压器的磁路结构型式、绕组接线方式(Y结、D结)、中性点接地与否等多种因素对励磁涌流的大小和波形有着较大影响，故本文仅对电力系统中最常见的Yd11、Ynd11、Yny0、Yy0（n表示中性点接地）接线的三相三柱心式变压器进行仿真研究。为简化分析，在研究变压器空载合闸哲态过程时忽略铁心的损耗，认为励磁支路为纯电感支路。

4. 1 三相变压器仿真的数学模型

首先对各种不同连接组情况下，根据电路原理的基础知识，建立三相变压器的数学模型。

4.1.1三相变压器Yd11连接组模式

图11为Yd11接线的变压器的三相接线图和单相等效电路。

uaLsrsubLsrsucLsrsLsrsr1L1RDLDuc

图4－1 Yd11 接线得变压器空载合闸时三项接线图和单相等效电路

8

三相变压器建模及仿真研究 当Y侧空载合闸后其暂态方程如下：

dia?u?(r?r)i?(L?L)?uN?eas1as1?adt?di? ?ub?(rs?r1)ib?(Ls?L1)b?uN?eb （4.1）

dt?dic??uc?(rs?r1)ic?(Ls?L1)dt?uN?ec? 式中，uN为Y侧中性点电压，其它符号见图3－1。

考虑到一次为Y接线，二次为D接线，所以有：

ia?ib?ic?0 (4.2)

diD?rDiD) (4.3) dt ea?eb?ec?3(LD而ua＋ub＋uc＝0，将式（4.1）三式相加并计及式（4.2）、式（4.3），化简得：

uN??(LD又由单相等效电路可知：

?ia?iD?ima;? ?ib?iD?imb; （4.5）

?i?i?iDmc?cdiDdt?rDiD) （4.4）

式（4.5）三式相加得到：

iD??(ima?imb?imc) （4.6）

31而

ej?d?dt?NjSdBjljdH2?dimjdt?Kj?dimjdt?Kj?imj(j?a,b,c) （4.7）

?式中 imj?―――电流imj的导数。

将式（4.4）、（4.5）、（4.6）、（4.7）代入方程式（4.1），计及一、二次绕组漏抗近似相等（r1＝rD，L1=LD）,经化简得：

9

三相变压器建模及仿真研究 ?3ua??3r1?2rs???3ub???rs?????3uc????rs?rs3r1?2rs?rs??ima???rs??imb?3r1?2rs???imc?rs??????

??3L1?3Ka?2Ls?Ls?Ls??i???Ls3L1?3Kb?2Ls?L??ma?i??s??mb????Ls?Ls3L?2L???1?3Kcs????imc??

若忽略系统阻抗，即rs＝0，Ls＝0，Ls0＝0，则上式可化简为：

?ua??r00??i ?u??1??ma????L1?Ka00??i?ma??b???0r10??imb???0L1?Kb0??i????mb???u?c???00r1????imc????00L???1?Kc????imc??2式中：Kjj?NjSdBl(j?a,b,c)????动态感应系数

jdH Ls，Ls0――――电源内部等值正序电感与零序电感 S，la，lb，lc――――变压器铁心截面积与各相磁路长度 L1，Na，Nb，Nc――――一次绕组漏抗和各相匝数 i??ma，imb，i?mc――――电流ima，imb，imc的导数

4.1.2三相变压器Ynd11连接组模式

Ynd11接线的三相变压器Yn侧空载合闸时，其暂态方程为：

??ua?(rs?r1)ia?(Ls?L1)dia?(Ls0?Ls)diD?ea?dtdt??ub?(rs?r1)ib?(LdibdiDs?L1)?(Ls?L)?eb (4.9) ?dt0sdt?di?u?(rcdi?cs?r1)ic?(Ls?L1)dt?(Ls0?Ls)Ddt?ec考虑到一次为Yn接线，二次为D接线，所以：

ia?ib?ic?3iD (4.10)

10

4.8） （

三相变压器建模及仿真研究 ea?eb?ec?3(LDdiDdt?rDiD) (4.11)

又ua?ub?uc?0，则式（4.9）三式相加得：

0?(Ls0?L1)di0dt?(rs?r1)iD?LDdiDdt?rDiD (4.12)

同样将式（4.5）三式相加得：

i0?13(ima?imb?imc)?iD (4.13)

将式（4.5）、（4.7）、（4.13）代入方程式（4.9）并联立式（4.12），若不计系统阻抗且认为变压器一、二次绕组漏抗相等，则可得其空载合闸状态方程：

?ua??r1????ub??0??u??0c????r??0???1?L1??????0r10r1000r1r1r1??ima????r1??imb?r1??imc??6r1????iD00L1?KcL1?????

?Ka00L1L1?Kb0L1*L1??ima????*L1??imb???L1??i*??mc??6L1???iD? (4.14)

4.1.3三相变压器Yny0连接组模式

Yny0接线的三相变压器Yn侧空载合闸时，其暂态方程与Ynd11接线一样，如式（4.9）所示。因为一次为Yn接线，二次为y接线，所以

?ia?ib?ic?3i0; (4.15) ??iD?0因而由单相等效电路可得：

ia?ima, ib?imb, ic?imc (4.16) 则 ima?imb?imc?3i0 (4.17)

同样，根据类似的推导过程并计及前述各假设，可得该接线三相变压器空载合闸状态方程：

11

三相变压器建模及仿真研究 ?ua??r1??? ?ub???0???uc???000??ima??L1?Kar10??imb??????0r1???imc??????0L1?Kb000??ima?????0??imb? (4.18)

????L1?Kc?i??mc??0?4.1.4三相变压器Yy0连接组模式

Yy0接线的三相变压器空载合闸时，其暂态方程与Ynd11接线一样，同样如式（4.1）所示。考虑到一次为Y接线，二次为y接线，则：

?ia?ib?ic?0; (4.19) ??iD?0从而可得：

ia?ima, ib?imb, ic?imc

又ua?ub?uc?0，将式（4.1）三式相加并计及式（4.19），化简得：

3uN??(ea?eb?ec) (4.20) 同样，根据类似的推导过程，可得Yy0接线得三相变压器空载合闸状态方程为：

?3ua??3r1???3ub??0?????3uc???003r10??3r1???imc?0??ima????0??imb????ima???*??Kc??imb??*?3L1?2Kc???imc?

?3L1?2Ka???Ka??Ka??Kb3L1?2Kb?Kb?Kc* (4.21)

至此，式（4.8）、（4.14）、（4.18）、（4.21）和各相动态磁化曲线Bj?f(Hj)及

Hj?(ImjNj)/lj(j?a,b,c)构成了Yd11、Ynd11、Yny0、Yy0接线三相变压一次侧空

载合闸的基本方程。

4.2电源电压的描述

根据前述假设，电源电压u（相电压）可用式（4.22）描述。

12

三相变压器建模及仿真研究 ?ua?(Um/1.73)?sin(?t??)?? ?ub?(Um/1.73)?sin(?t???120) (4.22)

???uc?(Um/1.73)?sin(?t???120)式中，Um为电源线电压峰值，取1.1倍额定电压。α为A相空载合闸初相角。在用Matlab仿真得过程中，α的设定并非是一个可以输入的变量，如果需要改变初相角，可在程序内部直接改变相电压u。

4.3铁心动态磁化过程简述

根据试验所得到变压器铁心磁化曲线数据分段拟和其极限磁滞回环是我们的基本原理。由试验所得到的数据可以帮助我们界定程序中一些参数，而如何选择界定函数将很大程度上影响试验仿真得结果。在这次的试验计划中，我们将选择两种方式（即选择不同的函数逼近）进行仿真，然后分别讨论两种方案的优劣，得出最佳的方案。

第一种是比较简单的模式，基本上不考虑曲线进入饱和区的情况（尽管饱和区是不可回避的问题，但这样做亦不失其合理性，这一点将在后面被讨论到。），采用两条修正的反正切函数做为极限磁滞回环。然后，对于主区间内的动态磁滞回环，根据不同的转折点和运行趋势对极限磁滞回环向极限磁滞回环拟合。极限磁滞回环的数学描述由于和第二种情况相近，只是将第二种方式的饱和区考虑在外，所以具体方法将不再赘述，可以参考4.1.3.1部分。对于暂态局部磁滞回环的描述，具体方法可以参考4.1.3.2部分。对于剩磁的处理的处理，具体方法可以参考4.1.3.3。

第二种是比较复杂得模式，需要在第一种的情况下考虑饱和的问题。这种方法不但描述了铁心的饱和特性，而且能够反映铁心的磁滞特性。基本原理是，首先，格局试验所得的变压器铁心磁化曲线数据分段拟合其极限磁滞回环：

（1）对于未饱和时主区间内的两条极限磁滞回环，采用修正的反正切函数加以拟合；

（2）对于饱和后主区间外的磁化曲线，认为其已进入线性可逆区（直线段），采用两条平行的直线段加以描述。然后，对于主区间内的动态磁滞回环，我们根据其不同的转折点和运行趋势对极限磁滞回环向饱和后的两条平行直线进行压缩，就可得变压器铁心实际运行的动态磁化轨迹。

13

三相变压器建模及仿真研究 下面，就以第二种方法为例，详细的解释一下极限磁滞回环的描述、暂态局部磁滞回环的描述及剩磁的处理等等问题。

4.3.1极限磁滞回环的数学描述

A．主区间内【－HZ,HZ】极限磁滞回环可用下式所示的修正反正切函数表示。

????tan?1?(H?C)??H (4.23)

????tan?1?(H?C)??H (4.24) 式中，参数?、?、?、C可根据实测磁滞回环数据由非线性的曲线拟合程序求得。 因此：

?1?? ????2dH??1????H?C?????dB?1?? ????2dH??1????H?C?????dB??? (4.25) ????? (4.26) ??

B．饱和后的磁化曲线（|H|>HZ）拟合为两条平行的直线段。 当H>HZ时，

B(H)?Bz?当H

B(H)??Bz?因而，

dBdH?Bs?BzHs?HzBs?BzHs?Hz?(H?Hz) (4.28)

Bs?BzHs?Hz?(H?Hz) (4.27)

(4.29)

14

三相变压器建模及仿真研究 4.3.2暂态局部磁滞回环的描述

由于铁心材料电磁性能的复杂性，对动态局部磁滞回环的精确仿真是比较困难的。但因极限 磁滞回环已描述了磁滞的基本轮廓，故根据不同转折点对其进行压缩就可近似模拟动态磁化过程中的某一段上升轨迹和下降轨迹。它分两种情况模拟。

A．dB/dH<0，运行点下降轨迹

由于极限磁滞回环左侧描述了减磁过程，将极限磁滞回环左侧回线在纵轴方向按比例地朝直线B?(?H??2)/?压缩，可得一簇下降曲线。

对通过某一转折点（H(0),B(0)）的运行点下降轨迹可由左极限磁滞回环按压缩系数KX向直线B?(?H??2)/?压缩而得（如图4－2所示）。此处KX＝(B－B2)/(B1－B2)则通

(0)

过该点的下降轨迹为：

?B?0??(?H(0)B???tan???1?H??0???1?????1??tan??H?C??????H??(4.30) ???2???2??C?2??2??)?因此：

dBdH?B?0??(?H(0)??2) ???1??tan?H???0??C????2????1?????H?C???2??? (4.31)

15

三相变压器建模及仿真研究

图4－2 局部磁滞回环的模拟

B．dB/dH>0，运行点上升轨迹

同样，将极限磁滞回环右侧向直线B?(?H?（H，B）的运行点上升轨迹

?B?0?(0)

(0)

?2)/?方向压缩，可得通过转折点

?(?H(0)B???tan???1?H??0???1?????1??tan??H?C??????H??(4.32) ???2???2??C?2??2??)?

dBdH?B?0??(?H(0)??2)???tan???1?H??0??C????2????1?????H?C???2??? (4.33)

4.3.3剩磁的处理

铁心中原始剩磁的存在，只体现为铁心的动态磁化轨迹（B－H曲线）在纵轴上起始

16

三相变压器建模及仿真研究 工作点的不同，运用如上所述的方法，可很方便地考虑剩磁的影响。

至此，变压器空载合闸暂态特征方程中的dB/dH（即Kj）已可求解，电源电压u又为已知，选用定步长四阶龙格－库塔法不难求解上列各微分方程。这样就建立了上述接线三相变压器空载合闸励磁涌流仿真数学模型。

第五章 三相变压器仿真的程序流程及功能介绍

由于在这次仿真模型设计中，先后采用了两种不同的思路。则分别按不同的思路列出程序流程图。

5.1分段拟和加曲线压缩法（两段修正的反正切函数）

5.1.1方法一

励磁涌流的计算可用图4－3所示的流程图来描述，图中主要部分的功能如下： 第一部分：主要是处理原始剩磁，确定初始化磁化轨迹。

T为时间，当时间刚大于零时，即程序开始运行时，沿第一部分，主要功能就是处理剩磁，将变量赋以初值，确定初始化磁化轨迹。

第二部分：主要是描述极限磁滞回环内的工作轨迹，确定动态磁滞回环。 其主要原则是：在迭代求解过程中，如某一时刻的磁感应强度变化方向（

dBdt）发

生改变，将前于该时刻的（H，B）作为下一局部磁滞回环的转折点，并按磁感应强度变化方向进行曲线压缩以确定下一时刻的局部磁化曲线。当然，如果磁链变化方向不变，则所运行的局部磁滞回环也不变。当磁链方向发生了变化，则改变局部磁滞回环。

17

三相变压器建模及仿真研究

(0)开始T>0 NH(0)?HBFirst(0)?BFirstYdBdBdtprewnow*dBdtnow?0NdtYY?0NHY(0)?HnowB?BnowdBdtnow?0YdBdtnow?0YNN按式（4-32）（4-33）解按式（4-30）（4-31）解用龙格库塔解NT?TYmax回复图4-3 方法一的仿真程序图

18

三相变压器建模及仿真研究 5.1.2方法二

励磁涌流的计算可用图4－4所示的流程图来描述，图中主要部分的功能如下： 第一部分：主要是处理原始剩磁，确定初始化磁化轨迹。

T为时间，当时间刚大于零时，即程序开始运行时，沿第一部分，主要功能就是处理剩磁，将变量赋以初值，确定初始化磁化轨迹。

第二部分：主要是描述极限磁滞回环内的工作轨迹，确定动态磁滞回环。 其主要原则是：在迭代求解过程中，如某一时刻的磁感应强度变化方向（

dBdt）发

生改变，将前于该时刻的（H，B）作为下一局部磁滞回环的转折点，并按磁感应强度变化方向进行曲线压缩以确定下一时刻的局部磁化曲线。当然，如果磁链变化方向不变，则所运行的局部磁滞回环也不变。当磁链方向发生了变化，则改变局部磁滞回环。

第三部分：主要是描述铁心深度饱和后以直线表示的工作轨迹，以便和第二部分结合。

19

三相变压器建模及仿真研究 T>0 开始NH(0)?HFirstB(0)?BFirstYYN zH?HNdBdtnowH?HY?0NzYdBdBdtprewprewdB*Ynowdt?0NH(0)dt*dBdtnow?0NNdBdtprewdB*Ynowdt?0Y(0)?HnowB?BnowdBdt(0)now?0YH(0)?HnowB?BnowdBdtnow?0YNN按式（4-28）（4-29）解按式（4-27）（4-29）解按式（4-32）（4-33）解按式（4-30）（4-31）解用龙格库塔解 N

T?Tmax Y回复图4－4 方法二的仿真程序流程图

20

三相变压器建模及仿真研究 5.2 三相变压器仿真的计算实例及结果分析

由于在这次仿真模型设计中，先后采用了两种不同的思路，则按不同的思路分别仿真并进行结果分析。

5.2.1方法一：用两段修正的反正切函数拟和压缩

某一额定容量为5000kVA、额定电压为35/6.3kV三相心式变压器，其一、二次绕组匝数744/232匝，一次绕组阻抗0.7925＋j0.07336?，铁心截面积721.14cm2，铁心磁路长度la?lc?291.5cm，lb?163.5cm。当改变压器高压侧空载合闸时，运用上述仿真模型就可对其励磁涌流进行仿真计算。仿真时铁心各相剩磁按Bra??Brb??Brc处理，原始剩磁Br分别取0.9Bm、0.7Bm、0.35Bm，A相合闸初相角?取0?、

?6仿真所得到的典

型涌流曲线如图4－5所示。因为就不同的连接组而言，波形是大致相同的，所以以Yd11为例，其余的波形将不再赘述。

图4－5 Yd11连接组模式的三相变压器典型涌流波形

21

三相变压器建模及仿真研究 仿真结果说明，励磁涌流波形为尖顶波，偏于时间轴一侧，含有大量的非周期分量和高次谐波，具有明显的间断角特性；而且在三相涌流重，往往有一相近似为周期性涌流。应用傅立叶变换（DFT）对其进行谐波分析可知，励磁涌流中二次谐波含量最大。

图4－6显示是在Yd11连接组模式的三相变压器磁滞回环的波形图。

图4－6 Yd11连接组模式的三相变压器磁滞回环波形

5.2.2方法二：用两段修正的反正切函数加两段直线拟和压缩

依然使用前面所述的变压器做试验仿真。假定某一额定容量为5000kVA、额定电压为35/6.3kV三相心式变压器，其一、二次绕组匝数744/232匝，一次绕组阻抗0.7925＋j0.07336?，铁心截面积721.14cm2，铁心磁路长度la?lc?291.5cm，lb?163.5cm。当改变压器高压侧空载合闸时，运用上述仿真模型就可对其励磁涌流进行仿真计算。仿真时铁心各相剩磁按Bra??Brb??Brc处理，原始剩磁Br分别取0.9Bm、0.7Bm、0.35Bm，A相合闸初相角?取0?、

?6仿真所得到的典型涌流曲线如图4－7所示。同样，因为就

不同的连接组而言，波形是大致相同的，所以以Yd11为例，其余的波形将不再赘述。

22

三相变压器建模及仿真研究

图4－7 Yd11连接组模式的三相变压器典型涌流波形

图4－8显示是在Yd11连接组模式的三相变压器磁滞回环波形图。

图4－8 Yd11连接组模式的三相变压器磁滞回环波形

5.1.3两种方法的比较分析

本文是通过两种方法对变压器进行仿真设计的，就原理而论，四条曲线拟和的方法相对复杂一些，但考虑问题也相对周密一些，可是在同一套参数下，得到的图形却和我们理想中的图形有一定的差距，主要是电流收敛的速度太快了。

23

三相变压器建模及仿真研究 用两条曲线拟和用四条曲线拟和

图4－9 两种方法比较下的励磁涌流波形图

我们可以看到，用了四条曲线拟和的方法中，励磁涌流衰减的速度要高于用了两条曲线拟和的方法。综合对波形的分析，用两条修正的反正切函数拟和的方法更相符于实际情况，相比之下，是较好的方法。

5.2影响变压器励磁涌流的主要因素及结果分析

虽然本文采用了两种方法去进行仿真，但实际上也是从不同的角度对变压器的励磁涌流进行分析，所以从两个方法出发进行分析影响变压器励磁涌流的主要因素，得到的结果都是近似的。所以下面仅以两条曲线进行模拟的方法为基础，在Yd11连接组模式下，进行影响变压器励磁涌流的主要因素分析。

变压器的励磁涌流不仅与其铁心材料、磁路结构形式、绕组的接线方式、中性点接地与否有关，而且与其铁心原始剩磁、合闸初相角以及电源电压幅值、系统等值内阻抗的大小等多种因素有关。下面对影响变压器励磁涌流的主要因素――剩磁、合闸初相角做仿真分析研究。

24

三相变压器建模及仿真研究 5.2.1剩磁对变压器励磁涌流的影响

以上述Yd11接线的变压器为例，在其它条件保持不变情况下(合闸初相角?＝30?)，铁心剩磁Br取0.9Bm、0.7Bm、0.35Bm进行仿真。所得典型涌流波形(A相)如图35所示，由图35可知，铁心剩磁越大，有的相励磁涌流越呈尖顶波，涌流幅值越大，其间断角、二次谐波含量越小；但有的相则相反。

0.9Bm 0.7Bm 0.35Bm

图4－10 Yd11不同剩磁下变压器典型涌流波形图

5.2.2合闸初相角对变压器励磁涌流的影响

同样以Yd11接线的变压器为例，在其它条件保持不变(剩磁Br＝0.7Bm)，合闸初相角?分别取0?、30?等进行仿真，所得典型涌流波形(A相)如图36所示。由图36可知，不同的?不仅对涌流幅值、间断角及二次谐波含量有显著影响，而且影响励磁涌流幅值

25

三相变压器建模及仿真研究 出现的时间。在其它条件相同时，?＝0?时合闸，Yny0、Yy0联结组变压器的A相涌流可达最大；而?＝30?时合闸，Ydll、Yndll接线的变压器的A相涌流可达最大。

??0 ???30?

图4－11 Yd11不同合闸初相角时变压器典型涌流波形

当然，变压器励磺涌流还与其它多种因素有关。根据仿真计算可得，空载合闸电源电压幅值越大、系统等值阻抗越小，则励磁涌流越大；铁心饱和磁通越小(即饱和特性越显著)、铁心磁路越长、铁心截面越小，则励磁涌流越大。

由上述仿真研究可知，变压器的励磁涌流具有明显的间断角特性，且二次谐波含量最大。

5.3三相变压器励磁涌流的特征

三相变压器空载合闸时，三相绕组都会产生励磁涌流。对于Yd11接线方式的三相变压器，引入每相差动保护的电流是两个变压器绕组电流之差，其励磁涌流也应该是两个

26

三相变压器建模及仿真研究 绕组励磁涌流的差值，即ie?A?R?ie?A?ie?B、ie?B?R?ie?B?ie?C、

ie?C?R?ie?C?ie?A。两个励磁

涌流相减后，涌流的时域特征和频域特征都有所变化。下面结合一个算例来说明它们的特点。计算条件：?m?1.1，?sat?1.15；三相的剩磁?r?A?0.7，?r?B??r?C??0.7；A相的合闸角?A?0。由于三相电压是对称的，故?B?4?/3，?C?2?/3。ie?A、ie?B和ie?C的波形如图6 (a)所示，ie?A?R、ie?B?R、ie?C?R的波形分别如图6 (b)、(c)和(d)所示。

ie?A23??43?2??tie?Bie?C(a)ie?A?R2?(b)ie?B?R2??t?t(c)2??tie?C?R(d)

图4－12 三相变压器励磁涌流波形

在图6 (a)中，要注意ie?A、ie?B和ie?C最大值出现的时刻：ie?A是正向涌流，在?t??时达到最大值；ie?B是反向涌流，故在?t?2?/3(即?t??B?2?)时达到最大值；ie?C也是反向涌流，最大值发生在?t?4?/3处。ie?A、ie?B、ie?C的间断角和二次谐波分别为；ie?A?R、

ie?B?R、ie?C?R的间断角和二次谐波则分别为78.6?、49.6?、78.6?和14.8％、37.6％、14.8％。

结合上面的算例，对于一般情况，三相变压器励磁涌流有以下特点：

1．由于三相电压之间有120o的相位差，因而三相励磁涌流不会相同，任何情况下空

27

三相变压器建模及仿真研究 载投入变压器，至少在两相中要出现不同程度的励磁涌流。

2．某相励磁涌流(ie?B?R)可能不再偏离时间轴的一侧，变成了对称性涌流。其它两相仍为偏离时间轴一侧的非对称性涌流。对称性涌流的数值比较小。非对称性涌流仍含有大量的非周期分量，但对称性涌流中无非周期分量。

3. 与单相变压器励磁涌流相比，其中一相或两相励磁涌流的二次谐波显著减小，但至少有一相励磁涌流仍含有大量的二次谐波。

4. 励磁涌流的波形仍然是间断的，但间断角显著减小，其中又以对称性涌流的间断角最小。但对称性涌流有另外一个特点：励磁涌流的正向最大值与反向最大值之间的相位相差120?。这个相位差称为‘波宽’，显然稳态故障电流的波宽为180?。

附录 Matlab程序

附录A.在Yd11接线方式下两段反正切函数拟和极限磁滞回环的程序

% 文件名：transa.m

% 目的：仿真电流互感器的暂态过程

%%%%%%%%%% 初始化%%%%%%%%%%%% clear; %清除变量 %磁滞回线起始点 B0a=0; I0a=0; B0b=0; I0b=0; B0c=0; I0c=0;

%变量：磁通 Ba=0; Bb=0; Bc=0;

%变量：励磁电感

K=[1242.496;1242.496;1242.496]; %RK法参数设定，此处为默认值 options=[];

deltaT=0.0001; %采样时间间隔

28

三相变压器建模及仿真研究 Tmax=0.4; %仿真时间

tspan = [0,deltaT]; %RK法时间宽度 yzero = [0;0;0]; %RK法初始值

%励磁电流变化率，象征磁滞回线的改换方向的点 DI0a=0; DI0b=0; DI0c=0;

%%%%%%%%%%%%%%%%仿真开始，是个循环过程%%%%%%%%%%5 for i=1:Tmax/deltaT

% 四阶RK解微分方程,注意此为定步长RK

[t,num_y] =rk4fixed(@transtatea,tspan,yzero,2,K); % 更新参数，为下次RK做准备

tspan=[t(length(t),1),t(length(t),1)+deltaT]; tt(i)=[t(length(t),1)]; yzero=num_y(2,:); yzero=yzero';

Im(i,:)=num_y(2,:); %取得励磁电流变化率

DI1a=diff(num_y(:,1))./diff(t); DI1b=diff(num_y(:,2))./diff(t); DI1c=diff(num_y(:,3))./diff(t); %取得当前状态的励磁电感值

[Ka,B0a,I0a,Ba]=getLm(Ba,Im(i,1),B0a,I0a,DI1a,DI0a,i); [Kb,B0b,I0b,Bb]=getLm(Bb,Im(i,2),B0b,I0b,DI1b,DI0b,i); [Kc,B0c,I0c,Bc]=getLm(Bc,Im(i,3),B0c,I0c,DI1c,DI0c,i); %为计算励磁电感更新参数 DI0a=DI1a; DI0b=DI1b; DI0c=DI1c;

%输出中间变量，电感和磁通 K=[Ka;Kb;Kc]; kk(i)=Ka; bb(i)=Ba; end;

%文件名：getLm.m

%功能:按磁滞回线计算动态电感 %(I0,B0):回环开始点

function [L,B0,I0,B]=getLm(B,I,B0,I0,DI1,DI0,begin) alfa=0.00377;

29

三相变压器建模及仿真研究 beta=20;

gama=2*0.016867; C=0.2227;

% 处理初始剩磁 if (begin==1) Hzero=I0; Bzero=B0; if(DI1<0)

sgn=1;

B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1/alfa

*(gama*I-pi/2);

L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(beta*(I+

C)))+gama/alfa;

else

sgn=-1;

B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(I-C))-pi/2)+1/alfa*(

gama*I+pi/2);

L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(beta*(I-C))*(beta*(I-C )))+gama/alfa; end;

% 结束处理起始磁通 else

% 处理回环时的状态 if(DI0*DI1<0) I0=I; B0=B;

if(DI1>0)

B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(I-C))-pi/2)+1

/alfa*(gama*I+pi/2);

L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(beta*(I-C))*(bet

a*(I-C)))+gama/alfa;

else

B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1/alf

a*(gama*I-pi/2);

L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(beta*(I

+C)))+gama/alfa;

30

三相变压器建模及仿真研究 end else

if(DI1<0)

B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1

/alfa*(gama*I-pi/2);

L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(bet

a*(I+C)))+gama/alfa;

else

B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(I-C))-pi/2)+1/a

lfa*(gama*I+pi/2);

L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(beta*(I-C))*(beta

*(I-C)))+gama/alfa;

end; end;

%结束处理回环时状态 end;

%文件名：rk4fixed.m

%功能：四阶定步长龙格－库塔法。

function [T,X]=rk4fixed(Fcn,Tspan,X0,N,K)

% Matlab implementation of a fixed-step RK4 algorithm. h = (Tspan(2)-Tspan(1))/N; halfh = 0.5*h; neqs=size(X0);

X=zeros(neqs(1),N); T=zeros(1,N); X(:,1)=X0; T(1)=Tspan(1); Td = Tspan(1); Xd = X0; for i=2:N,

RK1 = feval(Fcn,Td,Xd,K); Thalf = Td + halfh;

Xtemp = Xd + halfh*RK1;

RK2 = feval(Fcn,Thalf,Xtemp,K); Xtemp = Xd + halfh*RK2;

RK3 = feval(Fcn,Thalf,Xtemp,K);

31

三相变压器建模及仿真研究 Tfull = Td + h;

Xtemp = Xd + h*RK3;

RK4 = feval(Fcn,Tfull,Xtemp,K);

X(:,i) = Xd + h*(RK1+2.0*(RK2+RK3)+RK4)/6; T(i) = Tfull; Xd = X(:,i); Td = T(i); end

X=X';T=T'; End

%文件：trastate.m

% 目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程

function dIm = transtatea(t,Im,K) R1=0.7925; L1=0.07336; Um=100000;

dIm=(inv([(L1+K(1)),0,0;0,(L1+K(2)),0;0,0,(L1+K(3))]))*(Um*[sin(314*t);sin(314*t+pi*2/3);sin(314*t+pi*4/3)]-[R1,0,0;0,R1,0;0,0,R1]*Im);

附录B.在Ynd11接线方式下两段反正切函数拟和极限磁滞回环的程序

% 文件名：transa.m

% 目的：仿真电流互感器的暂态过程

%%%%%%%%%% 初始化%%%%%%%%%%%% clear; %%清除变量 %磁滞回线起始点 B0a=0; I0a=0; B0b=0; I0b=0; B0c=0; I0c=0; B0d=0; I0d=0;

%变量：磁通 Ba=0;

32

三相变压器建模及仿真研究 Bb=0; Bc=0; Bd=0;

K=[1242.496;1242.496;1242.496;0]; %变量：励磁电感 options=[];%RK法参数设定，此处为默认值 deltaT=0.0001; %采样时间间隔 Tmax=0.4; %仿真时间

tspan = [0,deltaT]; %RK法时间宽度 yzero = [0;0;0;0]; %RK法初始值

%励磁电流变化率，象征磁滞回线的改换方向的点 DI0a=0; DI0b=0; DI0c=0; DI0d=0;

%%%%%%%%%%%%%%%%仿真开始，是个循环过程%%%%%%%%%%5 for i=1:Tmax/deltaT

% 四阶RK解微分方程,注意此为定步长RK [t,num_y] =rk4fixed(@transtatea,tspan,yzero,2,K); % 更新参数，为下次RK做准备

tspan=[t(length(t),1),t(length(t),1)+deltaT]; tt(i)=[t(length(t),1)]; yzero=num_y(2,:); yzero=yzero';

Im(i,:)=num_y(2,:); %取得励磁电流变化率

DI1a=diff(num_y(:,1))./diff(t); DI1b=diff(num_y(:,2))./diff(t); DI1c=diff(num_y(:,3))./diff(t); DI1d=diff(num_y(:,4))./diff(t); %取得当前状态的励磁电感值

[Ka,B0a,I0a,Ba]=getLm(Ba,Im(i,1),B0a,I0a,DI1a,DI0a,i); [Kb,B0b,I0b,Bb]=getLm(Bb,Im(i,2),B0b,I0b,DI1b,DI0b,i); [Kc,B0c,I0c,Bc]=getLm(Bc,Im(i,3),B0c,I0c,DI1c,DI0c,i); [Kd,B0d,I0d,Bd]=getLm(Bd,Im(i,4),B0d,I0d,DI1d,DI0d,i); %为计算励磁电感更新参数 DI0a=DI1a; DI0b=DI1b; DI0c=DI1c; DI0d=DI1d;

%输出中间变量，电感和磁通 K=[Ka;Kb;Kc;0]; end;

33

三相变压器建模及仿真研究

%文件名：getLm.m

功能:按磁滞回线计算动态电感 可以使用§1中的getLm.m程序。

%文件名：rk4fixed.m

功能：四阶定步长龙格－库塔法。 可以使用§1中的rk4fixed.m程序。

%文件：trastate.m

% 目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程

function dIm = transtatea(t,Im,K) R1=0.7925; L1=0.07336; Um=100000;

dIm=(inv([(L1+K(1)),0,0,L1;0,(L1+K(2)),0,L1;0,0,(L1+K(3)),L1;L1,L1,L1,6*L1]))*(Um*[sin(314*t);sin(314*t+pi*2/3);sin(314*t+pi*4/3);0]-[R1,0,0,R1;0,R1,0,R1;0,0,R1,R1;R1,R1,R1,6*R1]*Im);

附录C.在Yny0接线方式下两段反正切函数拟和极限磁滞回环的程序

%文件名：transa.m

目的：仿真电流互感器的暂态过程 可以使用§1中的transa.m程序。

%文件名：getLm.m

功能:按磁滞回线计算动态电感 可以使用§1中的getLm.m程序。

34

三相变压器建模及仿真研究 %文件名：rk4fixed.m

功能：四阶定步长龙格－库塔法。 可以使用§1中的rk4fixed.m程序。

%文件：trastate.m

% 目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程

function dIm = transtatea(t,Im,K)

R1=0.7925; L1=0.07336; Um=100000;

dIm=(inv([(L1+K(1)),0,0;0,(L1+K(2)),0;0,0,(L1+K(3))]))*(Um*[sin(314*t);sin(314*t+pi

*2/3);sin(314*t+pi*4/3)]-[R1,0,0;0,R1,0;0,0,R1]*Im);

附录D.在Yy0接线方式下两段反正切函数拟和极限磁滞回环的程序

% 文件名：transa.m

目的：仿真电流互感器的暂态过程 可以使用§1中的transa.m程序。

%文件名：getLm.m

功能:按磁滞回线计算动态电感 可以使用§1中的getLm.m程序。

%文件名：rk4fixed.m

功能：四阶定步长龙格－库塔法。 可以使用§1中的rk4fixed.m程序。

35

三相变压器建模及仿真研究 %文件：trastate.m

% 目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程

function dIm = transtatea(t,Im,K)

R1=0.7925; L1=0.07336; Um=100000;

dIm=(inv([(3*L1+2*K(1)),-K(2),-K(3);-K(1),(3*L1+2*K(2)),-K(3);-K(1),-K(2),(3*L1

+2*K(3))]))*(3*Um*[sin(314*t);sin(314*t+pi*2/3);sin(314*t+pi*4/3)]-[3*R1,0,0;0,3*R1,0;0,0,3*R1]*Im);

附录E.在Yd11接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序

% 文件名：transa.m

目的：仿真电流互感器的暂态过程 可以使用§1中的transa.m程序。

%文件名：getLm.m

%功能:按磁滞回线计算动态电感 %(I0,B0):回环开始点

function [L,B0,I0,B]=getLm(B,I,B0,I0,DI1,DI0,begin)

alfa=0.00377; beta=20;

gama=2*0.016867; C=0.2227; slope=9; Hz=1.2; Bz=380;

% 处理初始剩磁 if (begin==1) Hzero=I0; Bzero=B0;

36

三相变压器建模及仿真研究 if(DI1<0)

B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(I+C))+pi/2)+1

/alfa*(gama*I-pi/2);

L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(beta*(I+C))*(bet

a*(I+C)))+gama/alfa;

else

B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(I-C))-pi/2)

+1/alfa*(gama*I+pi/2);

L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(beta*(I-C))*(

beta*(I-C)))+gama/alfa;

end;

% 结束处理起始磁通 else

% 处理回环时的状态 if (I>Hz)

if(DI0*DI1<0) I0=I; B0=B; end

B=Bz+slope*(I-Hz); L=slope; elseif(I<-Hz)

if(DI0*DI1<0) I0=I; B0=B; end

B=-Bz+slope*(I+Hz); L=slope; else

if(DI0*DI1<0) I0=I; B0=B; if(DI1>0)

B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(

I-C))-pi/2)+1/alfa*(gama*I+pi/2);

L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(bet

a*(I-C))*(beta*(I-C)))+gama/alfa;

else

B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(

I+C))+pi/2)+1/alfa*(gama*I-pi/2);

L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(bet

37

三相变压器建模及仿真研究 a*(I+C))*(beta*(I+C)))+gama/alfa;

end else

if(DI1<0)

B=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*(atan(beta*(

I+C))+pi/2)+1/alfa*(gama*I-pi/2);

L=(alfa*B0-(gama*I0-pi/2))/(alfa*(atan(beta*(I0+C))+pi/2))*beta/(1+(bet

a*(I+C))*(beta*(I+C)))+gama/alfa;

else

B=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*(atan(beta*(

I-C))-pi/2)+1/alfa*(gama*I+pi/2);

L=(alfa*B0-(gama*I0+pi/2))/(alfa*(atan(beta*(I0-C))-pi/2))*beta/(1+(bet

a*(I-C))*(beta*(I-C)))+gama/alfa;

end; end; end;

%结束处理回环时状态 end;

%文件名：rk4fixed.m

功能：四阶定步长龙格－库塔法。 可以使用§1中的rk4fixed.m程序。

%文件：trastate.m

目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§1中的trastate.m程序。

附录F.在Ynd11接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序

% 文件名：transa.m

目的：仿真电流互感器的暂态过程 可以使用§2中的transa.m程序。

38

三相变压器建模及仿真研究

%文件名：getLm.m

功能:按磁滞回线计算动态电感 可以使用§5中的getLm.m程序。

%文件名：rk4fixed.m

功能：四阶定步长龙格－库塔法。 可以使用§2中的rk4fixed.m程序。

%文件：trastate.m

目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§2中的trastate.m程序。

附录G.在Yny0接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序

% 文件名：transa.m

目的：仿真电流互感器的暂态过程 可以使用§3中的transa.m程序。

%文件名：getLm.m

功能:按磁滞回线计算动态电感 可以使用§5中的getLm.m程序。

%文件名：rk4fixed.m

功能：四阶定步长龙格－库塔法。

39

三相变压器建模及仿真研究 可以使用§3中的rk4fixed.m程序。

%文件名：trastate.m

目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§3中的trastate.m程序。

附录H.在Yy0接线方式下两段反正切函数加两段直线拟和极限磁滞回环的程序

% 文件名：transa.m

目的：仿真电流互感器的暂态过程 可以使用§4中的transa.m程序。

%文件名：getLm.m

功能:按磁滞回线计算动态电感 可以使用§5中的getLm.m程序。

%文件名：rk4fixed.m

功能：四阶定步长龙格－库塔法。 可以使用§4中的rk4fixed.m程序。

%文件：trastate.m

目的：得到系统微分方程 dIm=inv(L)*(U-r*IM),此文件适合于仿真涌流过程 可以使用§4中的trastate.m程序。

40

本文来源：https://www.bwwdw.com/article/c9kg.html