磁滞回线测量方法与Simulink仿真分析研究
第3l卷第3期2014年3月
机电
工
程
V01.3lNo.3Mar.2014
JoumalofMechanical&ElectricalEngineering
DOI:10.3969/j.issn.1001-4551.2014.03.025
磁滞回线测量方法与Simulink仿真分析研究
曹鸿泰,黄汝霖,姚缨英。
(浙江大学电气工程学院,浙江杭州310027)
摘要:针对变压器磁滞回线间接测量存在的电路参数复杂和测量结果准确性难以评价等问题,通过对RC积分器测量方法的电路时域微分方程求解和Simulink仿真结果的分析,提出了一种确定电路参数的方法。借助Simulink饱和变压器的模型,仿真测量了在不同电压下饱和变压器磁滞回线图像,进一步画出了基本磁化曲线;通过对比基本磁化曲线测量值与Simulink中模型值的差异,结果显示利用该方法能仿真得到饱和变压器的磁滞回线,其测量误差可用电路参数估算和控制。根据该方法设计了测量电路,得到了实际变压器的磁滞回线,验证了该方法的正确性。研究结果表明,该方法揭示了RC积分法测量的电路参数,并能预测测量误差,具有准确性高的特点,为非线性饱和变压器磁滞回线测量电路参数设计提供了参考依据。关键词:非线性饱和变压器;磁滞回线;RC积分法;Simulink仿真;电路参数设计中图分类号:TM402
文献标志码:A
文章编号:1001—4551(2014)03—0383—05
AnalysisofhysteresisloopmeasurementandSimulinksimulation
CAOHong—tai,HUANGRu-lin,YAOYing—ying
(College
ofElectrical
Engineering,Zhejiang
University,Hangzhou310027,China)
Abstract:Aimingattheproblemsofcircuitparametersandmeasurementaccuracyofthehysteresisloop,themethod
gratorcircuitparameters
to
determineRCinte—
Wasputforward.Themethod
was
based
on
theresultsofthetimedomaindifferentialequationandSimulinksimula・
tion.Agroupofhysteresisloopsweresimulatedatdifferentvoltagesaccordingtothemagnetization
curve
transformermodel
curve
in
Simulink.Then,thefundamental
error
was
plotted.ComparedwiththeSimulinkmodelcurve.thesimulation
on
wasplotted.The
ofmeasurementWas
hysteresisloop
estimableandcontrollablebythecircuitparameters.Based
was
themethod,anexperimentWasdesigned.Thereal
transformer
measured,andthemethodwasprovedright.The
resultsindicatethattheproposedmethodrevealsthecircuitparametersandpredictsthe
error
ofmeasurement.Themethodischaracterizedofhighaccuracy,andthereferencebasisofcircuitdesignisprovidedforhysteresisloop
measurementofnonlinearsaturabletransfoiTtler.
Keywords:nonlinearsaturabletransformer;hysteresisloop;RCintegrator;Simulinksimulation;circuitparameterdesign
0
引言
成果有:使用特斯拉计的磁滞回线测量Ll。2J,采用反相积分器∞州或锁相放大器¨1等运放电路测量磁滞回线,采用电容积分法PSpice仿真非线性磁芯磁滞回
线M’7J,采用示波器的RC积分法测量磁滞回线‘8引,利
磁滞回线是描绘铁磁性材料磁化性质的重要曲线,也是制造、选择铁磁材料的重要依据,在工程领域
有着重要的应用。测量磁滞回线的方法包括霍尔法等直接测量法和RC积分法等间接测量法。
目前学者对磁滞回线测量已经有一些研究,主要
用Matlab进行电网合闸磁滞回线仿真。1m¨J,也有利用
启发式算法测量和辨识磁滞回线¨2。”J。从中可以看出,直接测量往往需要昂贵的实验设备,操作较为繁
收稿日期:2013—10—22
作者简介:曹鸿泰(1992一)。男,河南南阳人,主要从事电力系统分析与控制方面的研究.E-mail:caohougtaillllll@163.com通信联系人:姚缨英,女,教授,博士生导师.E-mail:yaoyyzju@126.corn
・384・
机电
工
程
第31卷
琐;而间接法通过测量电压电流关系来获得磁滞回线,
其中积分测量法被广泛采用。积分法虽然具有线性度
H・f=N1・il-4-N2・i2
(1)
若变压器负载很小,即:
J7、rl・il》Ⅳ2・i2
(2)
好、准确度高、测量过程简便等优点,但是至今没有明确的参数计算方法。
实物实验往往需要昂贵的成本和较长的周期,而通过计算机软件仿真测量过程则可大大降低实验成本,提高测量效率。Simulink是Matlab中进行动态建
模与仿真的集成软件包,其中有针对电力系统的可视化建模与仿真工具库电力系统模块库(power
system
根据欧姆定律可得R。两端电压为:
旷掣.R12百。1
u1
(3)
Lj,
上式表明:满足条件式(2)时,u。与磁场强度H为
线性关系。
对副边由法拉第电磁感应定律,感应电势为:
blocksets,PSB),是解决电力系统中一些复杂的非线性仿真问题的有效方法。
铲一他・警=一Ⅳ2・警・s
为参考方向,根据KVL,副边回路满足:
(4)
本研究基于RC积分法测量磁滞回线,利用Mat一1ab的Simulink仿真系统中的电力系统模块库(power
system
若取原边电流流进同名端,副边电流流出同名端
blocksets,PSB)里的饱和变压器模型,并对RC
积分测量方法进行深人分析后,通过求解时域微分方
程,提出测量电路参数的确定方法。仿真结果的准确性分析以及简易实验结果表明,本研究提出的参数计
Ⅳ25警=(R2+r2)c:idu2+u2+L2C2可d2u2
当且仅当以下不等式:
(5)
算方法正确而且有效。
1
(R2+r2)c:百du2》“2+£2c:可d2u2
成立时:
N。S
(6)
测量原理
为了测量饱和变压器铁芯的磁滞回线,本研究设
dmB=(尺:+r:)c:百du2
(7)
解常微分方程,得:
u2
计的RC积分器测量电路如图1所示。该电路在变压器原边接入一个电阻(R,),通过电阻上的电压取样Ⅱ。,反映铁芯中磁场强度日的变化;在副边接人RC积
铲杀彘日2蕊而d
㈣
L为,
上式表明:满足条件(7)时u:与B是线性关系。在同时满足条件式(2,7)时,以M。和/2,:分别为X
分器(尺:、c:),通过RC参数的适当选取,利用电容c:
上的电压取样u:,反映铁芯中磁感应强度B的变化。
轴与Y轴作李萨如图形,由于它们分别与H和B呈线性关系,则李萨如图形就表示出饱和变压器的磁滞
回线。
最后,利用虚拟示波器观察“,和“:的李萨如图形即
得到了变压器铁芯的磁滞回线。由于是间接测量,需要分析u。,u:分别与日,日呈线性关系的条件。
圈
2基于Simulink的仿真测试
本研究利用Maflab(R20lOb版本)中的Simulink进行建模,建模步骤如下:
(1)建立一个新的模型窗口,同时打开Simulink的电力电子(SimPowerSystems)模块工具箱,进行仿真
测量模型的搭建。
(2)打开元件(Elements)模块组,复制一个饱和变
压器模型SaturableTransformer到建模窗口中,模型参数不变。
图1
测量饱和变压器磁滞回线电路原理图
2.1饱和变压器模型说明
打开Parameter选项卡,调节参数显示方式为SI
记原边线圈Ⅳ】匝,变压器截面周长为f,记副边线圈Ⅳ2匝,内阻为r:,变压器截面积为S。
对原边由安培环路定律得:
(即国际单位制),以下对饱和变压器的参数作出说明。饱和变压器模型参数如图2所示。
第1行数据表示饱和变压器的设备容量和工作频
第3期
曹鸿泰,等:磁滞回线测量方法与Simulink仿真分析研究
・385・
上的电压。
打开基本仿真(Simulink)模块工具箱,从输出(Sinks)模块组,复制一个XY型通道示波器XYGraph,设置x轴为±6e4,Y轴为±6e2,用于观察原边
取样电阻R.和副边RE积分器的电容C上的电压所形成的李萨如图形,即日一日曲线。
打开电源系统(SimPowerSystems)工具箱,复制一个powergui模块,双击打开,设置Configure
parameters
中的Simulationtype为Discrete,Sampletime设置为
50e6(单位:s)。这是因为仿真类型选择离散型,可以得到更加平滑的磁滞回线,仿真速度也更快。
图2饱和变压器模型参数
在菜单栏Simulation中,设置Configuration
Solver为ode23tb算法。
Para-
meters中的Stoptime至少大于0.04(单位:S)。设置
率,图中参数表示变压器的设备容量是250e6VA,工作频率是50
Hz。
进行适当连接后,即得到利用RC积分器观测磁
第2行数据表示饱和变压器的原边线圈的饱和电压、电阻和电感,图中参数表示变压器原边线圈的饱和
电压是7.35e5V,电阻是4.321
8
滞回线电路的模型。3
n,电感是0.45856
H。
RC积分器的参数设置
积分器测量方法需要满足式(3)和式(8)两个线
第3行数据表示饱和变压器的副边线圈的饱和电压、电阻和电感,图中参数表示变压器副边线圈的饱和电压是3.15e5V,电阻是0.793
8
性关系,则必须有式(2)和式(6)成立。当实验参数不满足线性关系时,虚拟示波器中不能得到正确的磁滞
Q,电感是0.084425
H。
第5行数据表示饱和变压器的模型参数,保持不变即可;第6行数据表示铁芯损耗电阻尺。,保持不变
即可。
回线。下面具体分析测量电路的参数选取。
饱和变压器模型铁芯电路模型中的有功损耗电阻,即为涡流损耗,由式(2)成立,则有i:一0A,这表明应使变压器工作在接近空载状态,此时满足u,与H为
线性关系。
2.2电路模块搭建与参数设定
打开电源(ElectricalSources)模块组,复制交流电压源模块AC
Voltage
Source到建模窗口中。仿真磁滞
在时域上求解微分方程,可找出不等式(6)成立
的条件。电源电压可用U。=A。sin(tot+90)表示,则电
回线,AC参数应当满足超过变压器的原边额定电压,才能使变压器磁路饱和。AC交流电源参数可以分别设定为10e5术sqrt(2)(单位:V),90(单位:deg),50(单位:Hz),Sampletime为0(单位:S)。
打开元件(Elements)模块组,复制3个串联RLC
支路模块到建模窗口中,下面对器件参数选择进行说明。R,为原边取样电阻,不应很大,可取1e3(单位:
容电压可表示为/Z:=K。sin(tot+妒),则不等式(6)可
化简为:
(尺2+r2)C2K。toeos(tot+咖)》K。sin(tot+咖)一L2C2K。t02sin(tot+币)
所示:
(9)
则应使后两项与第一项的比值尽可能小,如下式
Q),L为0(单位:H)、c为iIlf(单位:F),作为原边取样电阻尺。。其余两个串联RLC支路模块参数可以先设置尺:参数为尺为400e3(单位:Q)、L为0(单位:H)、C为inf(单位:F),设置c:参数尺为0(单位:n)、L为0(单位:H)、C为lOe-6(单位:F),作为副边的RC
积分器的电阻R:和电容C:。
打开测量(Measurements)模块组,复制两个电压测量模块VoltageMeasurement到建模窗口中,分别用于测量原边取样电阻尺。和副边RC积分器的电容C:
p(尺:,c:)=黜tan(c£,£+咖)(10)
由于式(10)右边三角函数tan(tot+咖)在一个周期内是发散的,为使比值尽可能小,只能令其系数趋近于0,才能满足式(8),从而电容电压“:与磁感应强度
B为线性关系。
通过Simulink进行饱和变压器磁滞回线的仿真采
用离散仿真模式,因此只需满足正切函数P(R:,C:)在一个周期里有足够长的时间近似为零,就可以得到饱和变压器铁芯的磁滞回线。为了保证仿真准确性,可
・386・
机电工程第3l卷
取上述条件的一个近似等价条件为JPl<0.1在99%
的周期内成立,即:
糍r2toC<—tan(0‰0
(尺2+)2、
暖
R2=1txF,R2=400
中,并将磁滞回线的拐点标出,用折线段拟合,再经过放大得到的磁滞回线图如图4所示,图4中包围面积即为铁芯损耗。由于变压器工作在轻负载状态,包围
面积变化很小。
99扎57川‘3.5耵・.)……~
(11)
H,
Simulink中饱和变压器的磁化曲线数据在hyster-esis.mat文件中,可以绘出模型的基本磁化曲线。根
据仿真结果得到的拐点可以绘出测量结果的基本磁化
通常1"2《R:可忽略,且在该测量中,L:一0.08
∞=314
rad/s,考虑变压器接近似空载运行,可设定R:
为1000倍饱和变压器满容量负载电阻,即:
曲线。由于存在线性关系,可将两幅图的坐标轴比例
(12)
000予。400
kQ
变换,合并为一张图,测量结果与模型值吻合,基本磁
化曲线测量误差对比如图5所示。
864
2
作为积分器电阻值,则解得C:>4.9心。因此,
可用C2=10
klq作为积分器参数。
上述参数仿真结果如图3所示。
基本磁化曲线对比图::饕羹饕誊耄羹
|Irr
IO
65
r
|lll
趔霞翼国龃
O8O6O4O20
;
2
圣一{
一-32
,
0
0
02004006008010012014
励磁电流相对值
:;
-6
图5基本磁化曲线测量误差对比
利用折线插值方法可以得到对应拐点处的相对误差值,其相对误差表明本研究提出的方法的测量结果与饱和变压器模型值非常接近。根据本研究提出的计算RC积分器参数方法得到的测量误差小于1.6%,而
且这个误差由式(15)控制。
因此,只有当R:、c:满足式(11,12)时,才满足测
图3积分器R:=400kn,C2=10心时测得磁滞回线
4仿真结果与分析
为了考察RC积分器测量磁滞饱和变压器磁滞回线的准确性,可以用多组磁滞回线绘出饱和变压器的
量要求中的线性关系,才能满足适用RC积分法的前提条件。
基本磁化曲线,并与Simulink中的磁滞回线的原始数据对比,以说明本研究所提出的实验模型的正确性和
仿真结果的准确度。
5
实验验证
为了进一步验证仿真方法的可操作性与正确性,
笔者分别调整AC电压源电压(单位:V)为9e5木
sqrt(2),9.5e5术sqrt(2),10e5木sqrt(2),10.5e5丰sqrt
根据本研究提出的参数确定方法,笔者设计了RC积分法测量电路,通过实验室测量得到了工频变压器的磁滞回线。
实验所用变压器为上海祥华电器厂生产的BK一
(2),11e5半sqrt(2)和11.5e5丰sqrt(2)进行仿真,可
以得到一组磁滞回线。将这些磁滞回线绘在一张图
50控制变压器,高压侧选择额定电压220V线圈,低
46
确龟
一
,
压侧选择额定电压36V线圈。高压侧经单相调压器
接人市电,串联电阻R.=240Q。低压侧串联RC积分法测量电路,其中R:=200kn,C:=4.7“F。本研究利用Tektronix
TDS
f|
|,
0
2
3054C示波器测量磁滞回线,
其中CHl测量尺。两端电压,CH2测量C:两端电压。
观测李萨如图形,BK一50变压器磁滞回线如图6
U/V
4
所示。
图4饱和变压器磁滞回线仿真结果(局部)
实验中测得磁滞回线与仿真结果相一致,从而验
第3期
!…!…l
120
曹鸿泰,等:磁滞回线测量方法与Simnlink仿真分析研究
…。-
・387・
’。●
80
40>
{0
o~
一40
■≯≥≥毒l
’...一i....1一...1一.一.1一-
o;….;....|....;.
!
i
;
一
一80
∥|l≤§
r’
参考文献(References):
[1]孟祥鹏,张建强,徐静,等.静态法磁滞回线测量实验分
析及微观机理解释[J].实验室研究与探索,2012,31
(9):40—44
∥≥≥≥0
二…:.…:….:….:….
’。
[2]王蕴杰.铁磁材料磁滞回线的Matlab分析方法设计[J].电脑开发与应用,2012,25(7):51-53.
一120
0≥|...∥
图6
[3]
j…。LJ
刘正平,程蔚.软磁材料磁滞回线在线检测系统[J].电测与仪表,2010,47(8):16.19.
o
一I
:一.一.:...・
[4]
吕哗,周洪强.磁滞回线测试系统的硬件电路研究[J].
0f/k/V
电测与仪表,2009,46(7):77—80.[5]SHIRANE
magnetic
T,ITOM.Measurementofhysteresisloopmaterials
on
Off
BK一50变压器磁滞回线
soft
using
lock—in
amplifier[J].IEEE
证了仿真方法的可操作性与RC积分器参数确定方法的正确性。
Transactions
Magnetics,2012,48(4):1437—1440.
[6]黄子平,王文斗.非线性含磁芯线圈的PSpice模拟[J].
强激光与粒子束,2004,16(8):1063-1066.
6结束语
在工程实践过程中,为了得到电力变压器的磁滞回线,需要经过复杂繁琐的操作过程,且往往受到试验条件、测量手段等方面的限制。电力变压器由于
额定电压在千伏以上,难以直接实验测量,而通常采
[7]PEDRA
puter
J,SAINZ
a
L,CORCOLESF,eta1.PSPICEcorn・
trans—,
modelof
nonlinearthree・-phasethree・・leggedTransactions
on
former[J].IEEE19(1):200—207.[8]
PowerDelivery,2004,
汪源浚.用示波器和信号发生器检测磁环B—H磁滞回线的方法[J].教学与科技,2009(3):7一12.
用仿真方法得到其磁滞回线。本研究基于Matlab的Simulink仿真系统,利用电力系统模块库PSB,采用RC积分法实现磁滞回线测量。本研究提出的确定参数方法使得最终测量相对误差在1%左右,具有很高的实用价值。
由于Simulink中饱和变压器模型给出的参数是基
本磁化曲线,原边电压存在较小的取值范围。如何解
[9]李都红,张小青,李敬怡,等.一种测量铁磁材料磁滞回线的方法及仿真[J].变压器,2008,45(4):38—39.
[10]袁兆强,凌燕.基于Matlab的变压器仿真建模及特性分析[J].电力学报,2007,22(2):175—179.
[11]肖伟平,黄绍平.考虑磁滞现象的变压器饱和特性仿真分析[J].高压电器,2005,41(4):254—256.
[12]何汉林,孟爱华,祝甲明,等.基于优化的GRNN和BP
神经网络的磁滞曲线拟合对比分析[J].机电工程,
2013,30(1):116—120.
决因原边电压取值问题引起的仿真错误还有待进一步
的研究。
[13]
MINCHEVS
systems
V.Neuralnetworksformodelingofdynamic
在有关交流电量电路分析中,常常需要结合时域分析与频域分析。本研究借鉴频域分析的思想,并就求解时域微分方程和Simulink仿真方法提出了确定参数的方法,为电路设计提供了参考依据。
with
hysteresis[C]//First
Inter—nationalIEEE
Symposium”Intelligent4247.
Systems”.Varna:[S.n.],2002:
[编辑:李辉]
本文引用格式:
曹鸿泰,黄汝霖,姚缨英.磁滞回线测量方法与Simulink仿真分析研究[J].机电工程,2014,3l(3):383—387.
CAOHong—tai,HUANGRu—lin,YAOYing—ying.AnalysisofhysteresisloopmeasurementandSimulink
simulation[J].Journal
ofMechanical&Electrical
Engineering,2014,31(3):383—387.《机电工程》杂志:http://www.meem.corn.cn
磁滞回线测量方法与Simulink仿真分析研究
作者:作者单位:刊名:英文刊名:年,卷(期):
曹鸿泰, 黄汝霖, 姚缨英, CAO Hong-tai, HUANG Ru-lin, YAO Ying-ying浙江大学电气工程学院,浙江杭州,310027机电工程
Mechanical & Electrical Engineering Magazine2014,31(3)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_jdgc201403025.aspx