MATLAB变压器仿真讲解

1、扬州大学专业软件应用综合设计报告水能 学院 13 级 电气专业题目变压器综合仿真设计二学生 某某某 学号 131504207指导教师 张建华2015年 12月 30 日目录 TOC o 1-5 h z 一、设计题目2二、正文2 HYPERLINK l bookmark6 o Current Document 1、引言22、设计依据及框图3设计平台3 HYPERLINK l bookmark11 o Current Document 设计思想4 HYPERLINK l bookmark13 o Current Document 设计结构框图或流程图 6 HYPERLINK l bookmark1

2、5 o Current Document 4各模块功能简介63、软件调试分析104、结语235、参考文献256、致谢25变压器综合仿真设计二摘要：随着变压器技术的进步， 传统仿真已经受到了很大的限制。并且当下要推动变压器技术的发展，已经不能再依靠传统仿真。因此，对于变压器的计算机仿真技术势在必行。本为通过MATLA瞅件，对变压器的运行特性进行了仿真。主要仿真的内容包括：变压 器磁路电流畸变以及变压器负载运行特性曲线研究。仿真所用到的方法为数值计算方法，通过插值的方法实现了对曲线的拟合。仿真时，结合实际情况可输入不同参数便于研究。文中给出了各种运行特性的仿真结果图，并且结合理论对其做了简单的分析

3、，验证了仿真方法的准确性和可行性。关键字：变压器；MATLA昉真分析；曲线拟合1引言设随着科学技术进步，电工电子新技术的不断发展，新型电气备不断涌现，人们使用 电的频率越来越高，人与电的关系也日益紧密，对于电性能和电气产品的了解，已成为人 们必需的生活常识。变压器是一种静止的电气设备，它是利用电磁感应原理把一种电压的交流电能转变成 同频率的另一种电压的交流电能，以满足不同负载的需要。在电力系统中，变压器是一个 重要的电气设备，它对电能的经济传输，灵活分配和安全使用具有重要的作用，此外，也 使人们能够方便地解决输电和用电这一矛盾。由于计算机仿真技术的出现，传统的物理仿 真系统逐渐的被计算机仿真系

4、统代替。计算机仿真系统所具有的效率高、精度高、重复性 和通用性好、容易改变仿真参数等优点，还可以实现物理仿真无法实现的有危险性的或者 是成本昂贵的仿真。在我国电力行业发展迅速的今天，变压器的仿真技术不能够再依托于 传统的物理仿真系统，而是需要能够采用能够促进变压器技术发展的仿真技术。对变压器特性的仿真涉及到很多方面，比如变压器空载励磁电流在饱和和磁滞影响时的特性、变压器磁滞回环在不同电压等级下的数据仿真、变压器空载合闸时的过电流现象、变压器在突发短路时的过电流现象，还有基本的比如效率特性、外特性、短路试验、空载试验等。在学习完本课程后，运用 MATLAB1关仿真技术对变压器进行仿真研究，本文的

5、仿真主 要以变压器磁路电流畸变以及变压器负载运行特性曲线为主要研究对象，通过结合实际进行曲线拟合、波形分析，得出相应结论。2设计依据及框图设计平台计算机技术的发展使得大量的数据计算变得方便快捷，一些因为需要不断的迭代而数据量庞大的数学算法也可以在实际中得到应用。不仅将工作者从繁忙的数据计算中解脱出来， 而且还可以做到不同精度的计算。MATLAB软件在数值计算方面独占鳌头，由于其提供了数据视图，文字处理的同一环境而受到欢迎。MATLAB勺中文意思为矩阵实验室，起源于 20世纪80年代，由其开创者 Cleve Moler 开发。经过后期的不断完善，MATLAB1终走向正轨，并且由MATHWORKS

6、司以商品形式发布。从MATLAB勺开创到现今，随着其版本的更替，功能也变得愈加强大。其核心编写所采 用的语言最终也从 FORTRA语言变为了 C语言。MATLAB! matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。 它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功 能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortr

7、an ）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数 值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图 像处理、信号检测、金融建模设计与分析等领域。MATLAB勺基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似， 故用MATLA味解算问题要比用 C, FORTRA除语言完成相同的事情简捷得多，并且 MATLAB 也吸收了像Maple等软件的优点，使 MATLAB

8、为一个强大的数学软件。在新的版本中也加 入了对C, FORTRANC+, JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序 导入到MATLAE0数库中方便自己以后调用，此外许多的MATLABS好者都编写了一些经典的程序，用户可以直接进行下载就可以用。其主要功能有：数值分析、数值和符号计算、工程与科学绘图、控制系统的设计与仿真、 数字图像处理、数字信号处理、通讯系统设计与仿真、财务与金融工程等MATLA由产生之日起就具有方便的数据可视化功能，以将向量和矩阵用图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用于科学计算和工程

9、绘图。新版本的MATLAB寸整个图形处理功能作了很大的改进和完善，使它不仅在一般数据可视化软件都具有的功能（例如二维曲线和三维曲面的绘制和处理等）方面更加完善，而且对于一些其他软件所没有的功能（例如图形的光照处理、色度处理以及四维数据的表现等），MATLAB同样表现了出色的处理能力。同时对一些特殊的可视化要求，例如图形对话等，MATLAB也有相应的功能函数，保证了用户不同层次的要求。另外新版本的 MATLA比着重在图形用户界面（GUI）的制作上作了很大的改善，对这方面有特殊要求的用户也可以得到满足。MATLAB对许多专门的领域都开发了功能强大的模块集和工具箱。一般来说，它们都是 由特定领域的专

10、家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。领域，诸如数据采集、数据库接口、概率统计、样条拟合、优化算法、偏微 分方程求解、神经网络、小波分析、信号处理、图像处理、系统辨识、控制系统设计、LMI控制、鲁棒控制、模型预测、模糊逻辑、金融分析、地图工具、非线性控制设计、实时快速 原型及半物理仿真、嵌入式系统开发、定点仿真、DSP与通讯、电力系统仿真等，都在工具箱（Toolbox ）家族中有了自己的一席之地。MATLAB包括拥有数百个内部函数的主包和三十几种工具包。工具包又可以分为功能性 工具包和学科工具包。功能工具包用来扩充MATLAB勺符号计算，可视化建模仿真

11、，文字处理及实时控制等功能。 学科工具包是专业性比较强的工具包，控制工具包，信号处理工具包，通信工具包等都属于此类。开放性使MATLAET受用户欢迎。除内部函数外，所有MATLA在包文件和各种工具包都 是可读可修改的文件，用户通过对源程序的修改或加入自己编写程序构造新的专用工具 包。现在的MATLAB已经发生了质的飞跃。完善的数值计算系统和简单的程序编写环境， 使得MATLA瞅件不需要太过专业的程序编写技能基础就可以轻松的使用。MATLAB矩阵作为运算的基本单位，使得矩B府！算变得方便快捷。此外，MATLA所提供的丰富的函数可以很容易的实现各种数值算法。MATLAB最为突出的特点就是编程语言的

12、简洁、直观。而且其对语法要求不是特别严格，不像其他编程语言，比如C语言等。MATLA的言程序文件为文本文件， 后缀为.m, 一般称为M文件。M文件具有保存和容易 修改命令的优点。 MATLA邈供了专门的 M文件编辑器。通过 M文件还可以自行的定义具有 具体功能的函数，使的程序的编写得到简化。MATLAB提供的数据可视功能为实现仿真的结果分析提供了方便，通过图像可以很容易 的对大量数据的变化实现分析。设计思想对变压器特性的仿真涉及到很多方面，比如变压器空载励磁电流在饱和和磁滞影响时 的特性、变压器磁滞回环在不同电压等级下的数据仿真、变压器空载合闸时的过电流现象、变压器在突发短路时的过电流现象，还

13、有基本的比如效率特性、外特性、短路试验、空载 试验等。变压器动态特性的分析主要是为计算和分析变压器在空载合闸和突然短路或者是其他4 故障时出现的暂态过电流和过电压。在电力系统中，变压器从发电厂到输配电网中都充当 着重要的角色。变压器的运行特性直接影响到电力系统的正常运行，所以对变压器的运行 特性进行研究是非常必要的。电力变压器的保护受到变压器励磁涌流的影响，单相变压器 空载合闸产生励磁涌流的大小与变压器合闸角有关，通过对励磁电流特点的研究可以制定 出防止励磁涌路引起误动的方法。随着人们对变压器技术的不断探索，变压器技术已从刚开始的基本感应定律发展到能够 对变压器的结构进行合理的设计、在理论上对

14、变压器运行时的各种现象进行详细分析、使变压器的理论模型达到更高的精确程度。在20世纪80年代末期，国外学者引入了磁化曲线和各频率下的等效激磁并联电导， 建立了变压器的频域模型， 但是在计算中存在了一定的误差， 并且数据参数的获取非常复杂。另一种方法为引入磁滞回线，建立时域模型，这种方法较为精确，但是在拟合磁化曲线建立较为准确的曲线模型存在着困难。有的学者通过传统的T型等效电路的励磁支路来实现变压器非线性的特性，但是在仿真中要求对铁芯损耗进行估 计，而且采用的计算方法相当的复杂，难以实现。另外有学者对变压器的非线性特性进行了详细的分析，但是并没有给出具体的模型和仿真结果。很多的文献中提到了分段磁

15、化曲线的简单模拟方法，但是通过这种方法的仿真结果不够精确，模型的准确程度不够。后来出现的通过拟合曲线的方法来准确的反应铁芯磁化的非线性。变压器的Jiles.Atherton 模型详细的说明了铁芯的磁化过程，并且通过磁学理论明了铁芯磁滞现象的原理。在这个基础上学者通过MATLA瞅件进行了仿真，仿真所涉到的五个模型参数包括饱和磁化强度、表明滞后磁 化曲线形状的参数、可逆磁化系数、磁畴对运动阻碍作用的参数、磁畴间相互作用的参数， 通过实验的方法可以获得。为了能够使计算速度更快， 没有采用经典的四阶龙格库塔算法来求解，而是采用了较为简单的欧拉算法。后来较为常见的描述变压器磁滞回线的方法是将磁滞回线分为

16、了主磁滞环和饱和区，这种描述可已较为精确的考虑到磁滞特性。通过拟合曲线的方法，将实验数据拟合成为接近实际的磁化曲线。主磁滞环采用修正的反正切函数加以拟合，饱和区特性将其视为线性可逆， 即磁化曲线用直线表示。在拟合出最大的磁滞回环后，通过对曲线按照比例的压缩，就可以得到一系列的反应不同程度的磁滞回线。通过龙格库塔算法求解非常系数的微分方程，考虑到了励磁阻抗随着饱和程度而变化地情况。随着仿真软件所提供的功能越来越强大，变压器的仿真可以省去编写程序的负担，而是直接用仿真软件中的功能块来实现，通过对模块的合理的联结和参数设置就可以完成仿真。在MATLA邸境下，基于单相变压器数学模型的基础，建立单独的功

17、能模块， 通过有机的整合功能模块搭建单相变压器的仿真模型。变压器以磁链为状态变量的数学模型与以电流作为状态变量的数学模型比较，前者在物理概念上更加的清晰。而且在采用数值积分的方法来计算时，此模型要更加稳定。上述所提到的利用数值计算来实现仿真的方法虽然都较为的成熟了，但是数值计算方法在仿真时对于初值的依赖较强。 如果初值的设定不合理，仿真程序很容易出现错误。甚至会得到发散的仿真结果，使得程序无止境的循环下去。本文将通过MATA既件对变压器的稳态和动态特性进行仿真。稳态特性仿真主要是考虑变压器磁路电流畸变以，动态特性仿真为变压器负载运行特性曲线研究。方法依托于MATLAB仿真软件的数值计算能力。

18、在对变压器磁路电流畸变数据处理时将用到三次样条插值拟合曲线的方法。动态特性仿真同样是利用了MATLAB勺数值计算功能。通过对仿真结果的分析，验证数值算法仿真的正确性和可行性。主要针对单相变压器负载运行时的参数分析，关注负载运行效率问题，研究负载运行效率与不同性质负载及不同负载系数之间的关系，通过仿真得出相关波形再进行讨论研究。设计结构框图或流程图变压器综合仿真不同负载时的外特性负载特性与效率曲线变压器电流畸变研究图1变压器综合仿真框图图2变压器综合仿真流程图各模块功能简介SIMULINK是一个对动态系统（包括连续系统、离散系统和混合系统）进行建模、仿真和综合分析的集成软件包，是MATLAB勺一

19、个附加组件，其特点是模块化操作、易学易用,而且能够使用 MATLAB1供的丰富的仿真资源。在SIMULINK环境中，用户不仅可以观察现实世界中非线性因素和各种随机因素对系统行为的影响，而且也可以在仿真进程中改变感兴趣的参数，实时地观察系统行为的变化，因此已在许多领域，如通信、信号处 理、DSR电力、金融、生物系统等获得重要应用。对电类专业的学生来说，无论是学 习专业课程或相关课程设计，还是在今后的工作中，SIMULINK都是一个重要的仿真建模工具。SimPowerSystems是在Simulink环境下进行电力电子系统建模和仿真先进工具。SimPowerSystems是Simulink下面的一

20、个专用模块库，包含电气网络中常见的元器件和 设备，以直观易用的图形方式对电气系统进行模型描述。模型可与其它Simulink模块的相连接，进行一体化的系统级动态分析。一、SimPowerSystems专用模块库的特点：.使用标准电气符号进行电力系统的拓扑图形建模和仿真；.标准的AC和DC电机模型模块；变压器；传输线；信号和脉冲发生器；HVDC控制； IGBT模块和大量设备模型，有断路器，二极管，IGBT , GTO, MOSFET和晶闸管；3.使用Simulink强有力的变步长积分器和零点穿越检测功能，给出高度精确的电力系统仿真计算结果.为快速仿真和实时仿真提供了模型离散化方法；.提供多种分析方

21、法，可以计算电路的状态空间表达、计算电流和电压的稳态解、设定 或恢复初始电流/电压状态、电力机械的潮流计算；.提供了扩展的电气系统网络设备模块，如电力机械，功率电子元件，控制测量模块和3相元器件；.提供36个功能演示模型，可直接运行仿真；二、SimPowerSystems专用模块库的强大功能：SimPowerSystems 中的模块SimPowerSystems中模块的数学模型基于成熟的电磁和机电方程，用标准的电气符号表 示。它们可以同标准的 Simulink模块一起使用建立包含电气系统和控制回路的模型。连接通过与SimPowerSystems提供的测量模块实现。SimPowerSystems

22、拥有近100个模块，分别位于 7个子模块库中。这些库模块涵盖了以 下应用范围：.电气网络(Electrical Sources & Elements )RLC支路和负载，兀型传输线，线性和饱和变压器，浪涌保护，电路分离器，互感，分7布参数传输线，3相变压器（2个和3个绕组），AC和DC电压源，受控电压源和受控 电流源。.电力机械（Machines）完整或是简化形式的异步电动机，同步电动机，永磁同步电动机，直流电动机，激磁系 统和水轮机涡轮机/调速系统模型。.电力电子（Power Electronics ）.二极管，简化/复杂晶闸管，GTO,开关，MOSFET, IGBT和通用型桥接管模型。.控

23、制和测量模块（Measurements）电压、电流和电抗测量，RMS测量，有功和无功功率计算，计时器，万用表，傅立叶分析，HVDC控制，总谐波失真，abc到dq0和dq0到abc轴系变换，3相V-I测量，3 相脉冲和信号发生，3相序列分析，3相PLL和连续/离散同步6-, 12-脉冲发生器。.三相网络元器件（Electrical Sources & Elements ）3相RLC负载和支路，3相断路器，3相，3相电抗，兀型传输线，AC电压源，6-脉冲 二极管和晶闸桥管，整流二极管，Y-A/Y-A/Y-A/Y-Y-A可配置3相变压器。三、SimPowerSystems常用模块库简介. Conti

24、nuous:（连接器元件库有 10种模块）.Electrical Sources:（电源元件库有 7种电源功能模块）DC Voltage Source:直流电压源AC V oltage Source:交流电压源AC Current Source:交流电流源Controlled V oltage Source:受控电压源Controlled Current Source:受控电流源3-Phase Source:三相电源3-Phase Programmable Voltage Source:三相可编程电压源Elements:（线路元件库有 24种模块）Machines:（电机元件库有 16种模块）

25、Measurements:（电路测量模块元件库有5种模块）而本文所采用的仿真方法大多是是利用m文件进行仿真研究。下面就这种方式及相关知识进行介绍。Matlab输入命令的常用方式有两种：一种是直接在Matlab的命令窗门中逐条输入Matlab命令；二是m文件工作方式。当命令行很简单时，使用逐条输入方式还是比较方8便的。但当命令行很多时（比如说几十行乃至全成百上千行命令），显然再使用这种方式输入MATLAB命令，就会显得杂乱无章，不易于把握程序的具体走向，并且给程序的修改和维 护带来了很大的麻烦。这时，建议采用Matlab命令的第二种输入形式 m文件工作方式。m文件工作方式，指的是将要执行的命令全

26、部写在一个文本文件中，这样既能使程序显得简洁明了，又便于对程序的修改与维护。m文件直接采用 Matlab命令编写，就像在 Matlab的命令窗口直接输入命令一样，因此调试起来也十分方便，并且增强了程序的交互性。m文件与其他文本文件一样，可以在任何文本编辑器中进打编辑、存储、修改和读取。利用m文件还可以根据白己的需要编写一些函数，这些函数也可以橡Matlab提供的函数一样进行调用。从某种意义上说，这也是对MATLAB的二次开发。m文件有两种形式：一种是命令方式或称脚本方式；另一种就是函数文件形式。两种形式的文件扩展名均是.mo1、M文件当遇到输入命令较多以及要重复输入命令的情况时，利用命令文件就

27、显得很方便了。将所有要执行的命令按顺序放到一个扩展名为.m的文本文件中，每次运行时只需在命令窗口输入m文件的文件名就可以了。需要注意的是，m文件最好直接放在 Matlab的默认搜索路径下（一般是Matlab安装目录的子目录 work中），这样就不用设置 m文件的路径了，否则 应当用路径操作指令 path重新设置路径。另外， m文件名不应该与 Matlab的内置函数名以 及工具箱中的函数重名，以免发生执行错误命令的现象。Matlab对命令文件的执行等价于从命令窗口中顺序执行文件中的所有指令。命令文件可以访问Matlab工作空间里的任何变量及数据。命令文件运行过程中产生的所有变量都等价于从Matl

28、ab工作空间中创建这些变量。因此，任何其他命令文件和函数都可以自由地访问这些变量。这些变量一旦产生就一直保存在内存中，只有对它们重新赋值，它们的原有值才会变化。关机后，这里变量也就全部 消失了。另外，在命令窗口中运行clear命令，也可以把这些变量从工作空间中删去。当然，在Matlab的工作空间窗口中也可以用鼠标选择想要删除的变量，从而将这些变量从工作空 间中删除。2、M函数m函数文件是一个特殊的 m文件，其常见格式如下：funcdon返问变量列表=函数名（输入变量列表）注释说明语句段函数体语句。需要说明的是，这里输入变量的个数以及输出变量的个数是由MATLAB本身提供的两个保留变量 narg

29、in 和 nargout 来给出的，它们分另1J是Number of function input arguments 和Number of function output arguments的缩写形式。输入变量要用逗号隔开，输出变量多于个时，要用方括号括起来。用户可以借助于help命令显示其中的注释说明语句段。通过这样的方法就可以建立函数文件或者称m函数，其调用方法与一放的Matlab函数的调用方法相同。函数文件相当于对 Matlab进行了二次开发。其作用与其他高级语言子函数的作用基本相同，都是为了实现特定目的而由用户自己编写的子函数。函数文件与命令文件有着鲜明的区别：函数文件的第一行必须包

30、含function字符；命令文件无此要求。函数文件的第一行必须指定函数名、输入参数及输出参数，命令文件无此要求。一个函数文件可以合 0个、1个或多个输入参数和返回值。函数文件要在文件的开头定义函数名，如function y1,y2=func(x,a,b,c),则该函数文件名必须存为func.m ,而命令文件无此要求。命令文件的变量在文件执行结束以后仍然保存在内存中而不会丢失，而函数文件的变量仅在函数运行期间有效(除非用global把变量说明成全局变量，否则函数文件中的变量均为局部变量)，当函数运行完毕后，这些变量也就消失了。需要说明的是，调用函数时所用的输入输出变量名并不要求与编写函数文件时所

31、用的输入输 出变量名相同。3软件调试分析励磁电流的研究图3磁通、磁化曲线、励磁电流波形图铁芯在反复磁化时，磁化特性曲线就会偏离基本磁化曲线，励磁电流会滞后于磁通的变化。磁滞损耗使得励磁电流的波形超前于磁通波形。实际上铁芯中涡流产生的损耗也会使得10励磁电流超前于磁通波形。空载励磁电流波形除了受到铁芯磁化的非线性影响外，还受到铁芯磁滞现象的影响。如 果铁芯在剩磁为零时开始磁化，那么铁芯磁化特性将沿着基本磁化曲线进行，但是当磁通开始发生转折后，磁化特性将不再按照基本磁化曲线进行，而是磁通变化滞后于励磁电流的变化，这种现象为铁芯的磁滞现象。由于铁芯磁滞想象的存在，铁芯的磁化呈现为基本磁化曲线两侧的回

32、线。在这种情况下，同一个数值的磁通会对应两个励磁电流。下图以为MATLAB程序绘制的磁磁滞回线。有一实际例子磁滞回线磁通和电流数据为分别为：。1=卜1.0 ,-0.90 , -0.79 , -0.50 , 0.30 , 0.70 , 0.85 , 0.91 , 0.95 , 0.97 , 1.0 , （）2=-1.0 , -0.97 , -0.91 , -0.85 , -0.70 , -0.30 , 0.50 , 0.79 , 0.90 , 1.0 , i = -0.05 : 0.01 : 0.05 ,单位为.p.u 。 其中（P与i对应磁滞回线中的下降曲线，。2与i对应磁滞回线中的上升曲线。

33、当输入磁通信号时，依旧用磁通信号的数据对磁滞回线进行插值计算。将插值计算的结果可视化输出，滞现象的影响下，耗，称之为磁滞损耗。形。_融发回海jIIIF=:F S言:t :：；J ：“一；i|I-I- .1 J . j3 - r1 5*=i&*j-.1-n DBJans - D4DDEI2 QO4 DCE OCB I加时电液即U图4磁滞回线图即得到受到铁芯磁滞影响的励磁电流波形图。励磁电流波形超前于磁通波形。当改变磁通信号幅值时，从图中知道励磁电流在铁芯磁这说明铁芯的磁滞现象在铁芯中造成了损 就会得到受不同程度磁滞影响的励磁电流波,唐里滞悬的加码屯赧漉带图5考虑磁滞后的电流波形图变压器磁路电流畸

34、变程序设计变压器铁芯磁化数据如表所示，设磁路中磁通按正弦规律变化，分析励磁电流的变化规律。铁磁材料磁化数据1234567891011121314151617I1.381.541.711.912.162.502.933.494.225.216.528.3610.914.822.437.864.00.400.480.560.640.720.800.880.961.041.121.201.281.361.441.521.601.68(1)研究磁路中按正弦规律变化的磁通，研究磁通与时间的关系，即虫=。M文件具体如下:idata=1.38,1.54,1.71,1.91,2.16,2.50,2.93,3.

35、49,4.22,5.21,6.52,8.36,10.91,14.8,22.4,37.8,64.0;phidata=0.40:0.08:1.68;t=0:0.001:pi;phi=2*sin(t);11i=spline(phidata,idata,phi);plot(t,phi)xlabel(t/s);ylabel(/wb)title( =f(t)可得仿真波形，从仿真波形可以看出磁路中磁通按正弦规律变化。=f(t)21.81.61.41.210.80.60.40.20.52.53.51.5 t/s图6磁路磁通图(2)研究磁通 与电流的关系，即二杷)。M文件具体如下：idata=1.38,1.54

36、,1.71,1.91,2.16,2.50,2.93,3.49,4.22,5.21,6.52,8.36,10.91,14.8,22.4,37.8,64.0;phidata=0.40:0.08:1.68;t=0:0.001:pi;phi=2*sin(t);i=spline(phidata,idata,phi);plot(i,phi)xlabel(I/A)ylabel(/Wb)title( =f(I)axis(0 400 0 3)12=f(I)2.50 LL1LJ050100150200250300350400I/A图7磁通电流关系图(3)研究电流的畸变规律，即 I=f(t)idata=1.38,1

37、.54,1.71,1.91,2.16,2.50,2.93,3.49,4.22,5.21,6.52,8.36,10.91,14.8,22.4,37.8,64.0;phidata=0.40:0.08:1.68;t=0:0.001:pi;phi=2*sin(t);i=spline(phidata,idata,phi);plot(t,i)xlabel(t/s)ylabel(I/A)title(I=f(t)axis(0,3.5,0,400)13400I=f(t)350 300 250 A 200 150 100 50 000.511.522.533.5t/s图8畸变电流图 小结：由于变压器的铁芯存在着饱

38、和现象，即铁芯磁化的非线性现象，使得磁通与。与励磁电流之间的关系。=f(i 0),呈现非线性。磁通和电压按照同样的规律变化，即正弦变化，只是。要滞后电压0.5*pi的电角度。根据磁化曲线可以求得对应每一时刻的励磁电流，从 而绘制励磁电流波形图。上面的波形图为正弦变化地磁通曲线、基本磁化曲线和励磁电流 曲线。图中励磁电流波形偏离了正弦波，呈现为尖顶波。3、变压器负载运行变压器的负载运行是指原绕组接入电源电压，副绕组接负载时的工作状况。这时，变压器的副边也有电流流通， 原边的接入电路与空载相比相应增大，副边端电压将受到负载的影响而发生变化，这是负载运行与空载运行的主要区别。在变压器负载运行时，由于

39、变压器的一次侧漏阻抗很小，即使负载时的一次侧电流变为了空载电流的很多倍， 依旧存在了 l1nZ1U的关系。所以忽略一次侧电压在漏阻抗上的电压降落，认为U = E1=4.44fN1m,所以负载运行时的磁通与空载运行时的磁通相差无几。二次侧电流产生的磁动势势必会影响到励磁磁动势，要使得缠身磁通的励磁磁动势不发生变化，14一次侧必须有中和二次侧电流产生的磁动势的电流。所以负载运行时，一次侧电流可以将其 视为励磁分量和负载分量之和。这里直接给出变压器的基本方程，如下：U1 =E1 II乙5=乐一|左2虫一 k E2I2 T =k-由ZmU2 = iI2zl将二次侧量折合到一次侧，基本方程变为下式:U1

40、 二一由 I1Z1U 2 = E2 - I2z2Ei =住2，1 I2 = II0I =20 ZmU2 =以(3-1)(3-2)根据式（2-10），可以得到变压器的 T形等效电路，如图9。负载时的一次侧电流要比励 磁电流大很多，如果忽略励磁电流， 便可简化T形等效电路为变压器的简化等效电路。如图10。另Zk =乙+Z；,电路图变为如图11。 TOC o 1-5 h z hR1jX1R2jX2OIIU1丫 山=也2U2 ZL HYPERLINK l bookmark46 o Current Document jXm幸 * oX._图9 T形等效电路15变压器在空载时的二次侧电压为其额定电压。而当

41、变压器二次侧接负载时，变压器的二次侧电压会发生变化，用 AU表示：AU =(U2N -U2 )/U2N *100%(3-3)根据给出的等效电路图和简化向量图，给出另一个AU的计算公式：U = P * (R *cos 平2+Xk *sin 平2 )(3-4)其中P为负载因数，R和为短路电阻和短路电抗，邛2为短路阻抗的阻抗角。变压器带感性负载时电压降低，带容性负载时电压升高。变压器效率受到铁损和铜损的影响。由于变压器空载和负载时铁芯中主磁通基本不变，所以我们认为变压器铁损基本不变。而变压器铜损随着变压器负载的变化会发生变化。设额定电流时，变压器的空载损耗为pkN，那么我们可以将铜损表达如下：Pcu

42、=B2*PkN(3-5)这里我们给出变压器效率的计算公式：=1 -初 + P2* PkN) /(P* Sn*cos % + P0 + P2* PkN)(3-6 )本次课程设计中，针对上述理论结合实际利用m文件编程，所编写的程序应用性强，可适用于不同规格型号的变压器，只需提供输入相关参数即可。下面是各程序的主体：(1)求解变压器相关参数数值的程序% 俞入基本参数SN=input(请输入额定容量 SN=);U1N=input(请输入一次侧额定电压U1N=);U2N=input(请输入一次侧额定电压U2N=);r1=input(请输入一次绕组漏阻值r1=);16r2=input（,请输入二次绕组漏阻

43、值r2=）;x1=input（,请输入一次绕组漏抗值x1=）;x2=input（,请输入二次绕组漏抗值x2=）;rm=input（请输入励磁电阻 rm=）;xm=input（请输入励磁电抗 xm=）;ZL=input（,请输入负载阻抗 ZL=）;%十算额定电流和变比I1N=SN/U1N;I2N=SN/U2N;k=U1N/U2N;%十算T型等效电路中的未知参数，其中Z1为原边漏阻抗，rr2为副边漏阻抗折算值，xx2*rr2为副边漏电抗折算值%ZZ劝副边漏阻抗折算值，ZZL为负载阻抗折算值，Zm为励磁阻抗，Zd为输入阻抗Z1=r1+j*x1;rr2=kA2*r2;xx2=kA2*x2;ZZ2=rr

44、2+j*xx2;ZZL=kA2*ZL;Zm=rm+j*xm;Zd=Z1 + 1/（1/Zm+1/（ZZ2+ZZL）;%U1I为原边电压，I1I原边电流，E1I励磁电动势，I22I折算后副边电流，I2I副边电流，U22I折算后副边电压%U2I副边电压U1I=U1N;I1I=U1I/Zd;E1I=-（U1I-I1I*Z1）;I22I=E1I/（ZZ2+ZZL）;I2I=k*I22I;U22I=I22I*ZZL;U2I=U22I/k;%cospsil输入侧功率因数，cospsi2负载侧功率因数，P1输入有功功率，P2输出有功功率17cospsi1=cos(angle(Zd);cospsi2=cos(

45、angle(ZL);P1=abs(U1I)*abs(I1I)*cospsi1;P2=abs(U2I)*abs(I2I)*cospsi2;eta=P2/P1;%ImI励磁电流，pFe铁损耗，pcu1原边铜损耗，pcu2副边铜损耗ImI=E1I/Zm;pFe=abs(ImI)A2*rm;pcu1=abs(I1I)A2*r1;pcu2=abs(I2I)A2*i2;%俞出数据disp(原边电流=),disp(abs(I1I)disp(副边电流=),disp(abs(I2I)disp(副边电压=),disp(abs(U2I)disp(原边功率因数=),disp(cospsi1)disp(原边功率=),d

46、isp(P1)disp(副边功率因数=),disp(cospsi2)disp(副边功率=),disp(P2)disp(效率=),disp(eta)disp(励磁电流=),disp(abs(ImI)disp(铁损耗=),disp(pFe)disp(原边铜损耗=),disp(pcu1)disp(副边铜损耗=),disp(pcu2)下面通过一个具体事例进行分析，有一单相变压器参数：SN=10kVA, U1n/ Un=380/220V,门=0.14Q , r2=0.035 Q , x1=0.22 Q , x2=0.055 Q , rm=30Q , xm=310Q。通过输入相关参数可以得到仿真结果为：请

47、输入额定容量 SN=10e3请输入一次侧额定电压U1N=380请输入二次侧额定电压 U2N=22018请输入一次绕组漏阻值r1=0.14请输入二次绕组漏阻值r2=0.035请输入一次绕组漏抗值 x1=0.22请输入二次绕组漏抗值 x2=0.055请输入励磁电阻rm=30请输入励磁电抗xm=310请输入负载阻抗ZL=4+j*3原边电流=25.5752副边电流=42.7447副边电压=213.7237原边功率因数=0.7725原边功率=7.5072e+003副边功率因数=0.8000副边功率=7.3085e+003效率=0.9735励磁电流=1.1998铁损耗=43.1836原边铜损耗=91.57

48、2519副边铜损耗=63.9489可以看出，上面的程序对于求解变压器参数很便捷，只要输入一些已知的基本参数便可求 出另外的一些重要参数，另外程序在编写的过程中也体现出“自动”的概念，及实现了输入参数与程序相结合的方式，而不仅仅简单的把数据直接编写在程序中，实现了一定的“智能化”。（2）变压器不同负载性质运行时的效率的程序M文件程序主体：c0=input（ 请输入阻性功率因数c0=）;s0=input（请输入阻性因数s0=）;c1=input（ 请输入感性功率因数c1=）;s1=input（请输入感性因数s1=）;c2=input（,请输入容性功率因数c2=）;s2=input（请输入容性因数s

49、2=）;xk=input（请输入短路电抗xk=）;rk=input（请输入短路阻抗rk=）;x=0:0.1:1.5;du2=x*（rk*c2+xk*s2）;du1=x*（rk*c1+xk*s1）;du0=x*（rk*c0+xk*s0）;u0=1-du0;u1=1-du1;u2=1-du2;plot（x,u2,x,x,u1,o,x,u0,x,1,-.）;xlabel（负载系数3 ）;ylabel(u2/.p.u);20title（,变压器的外特性,）;legend（,感性负载,容性负载,纯电阻负载,3）下面以一个例子加以说明，若有一变压器不同负载参数如下 c2=0.8,s2=0.6,xk=0.0

50、504,rk=0.024;c1=0.8,s1=-0.6;c0=1,s0=0;对其进行仿真可以得至 U下面的波形。u -DT2U119 8 7 口,9,9,91 o o O载 载载负 负负阻 性性电 感容纯,5 O数系载负21图12外特性图（3）变压器运行效率与负载系数的程序M文件程序如下: f=inline(x*sn*c/(x*sn*c+p0+x.A2*pkn),sn,c,p0,pkn,x);c=input(请输入功率因数值c：);sn=input(,请输入额定容量值sn：);p0=input（,请输入空载损耗值p0:）;pkn=input（,请输入额定损耗值pkn :）;n=1;for x=0:0,01:1.5h(n)=f(sn,c,p0,pkn,x);n=n+1;endx=0:0.01:1.5plot(x,h)axis(0 1.5 0 1)ymax=max(h)xm=find(diff(sign(diff(h)=-2)+1;xmax=xm(1)*0.01plot(x,h,xm(1)*0.01,max(h),r*)xlabel(负载系数3 )ylabel(变压器效率刀,)title( 效率特性曲线,)text(xm(1)