MATLAB Simulink系统建模与仿真 实验报告.docx

《MATLAB Simulink系统建模与仿真 实验报告.docx》由会员分享，可在线阅读，更多相关《MATLAB Simulink系统建模与仿真 实验报告.docx（23页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、MATLAB/Simulink电力系统建模与仿真实验报告姓名：*专业：电气工程及其自动化班级： *学号：*实验一 无穷大功率电源供电系统三相短路仿真1.1 无穷大功率电源供电系统仿真模型构建运行 MATLAB 软件，点击 Simulink 模型构建，根据电路原理图，添加下列模块：（1） 无穷大功率电源模块（Three-phase source）（2） 三相并联 RLC 负荷模块（Three-Phase Parallel RLC Load）（3） 三相串联 RLC 支路模块（Three-Phase Series RLC Branch）（4） 三相双绕组变压器模块（Three-Phase Tran

2、sformer (Two Windings)）（5） 三相电压电流测量模块（Three-Phase V-I Measurement）（6） 三相故障设置模块（Three-Phase Fault）（7） 示波器模块（Scope）（8） 电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下：1.2 无穷大功率电源供电系统仿真参数设置1.2.1 电源模块设置三相电压 110kV，相角0，频率50Hz，接线方式为中性点接地的Y 形接法，电源电阻 0.00529，电源电感 0.000140H，参数设置如下图：1.2.2 变压器模块变压器模块参数采用标幺值设置，功率20MVA，频率

3、 50Hz，一次测采用 Y 型连接，一次测电压 110kV，二次侧采用 Y 型连接，二次侧电压 11kV，经过标幺值折算后的绕组电阻为 0.0033，绕组漏感为 0.052，励磁电阻为 909.09，励磁电感为 106.3，参数设置如下图：1.2.3 输电线路模块根据给定参数计算输电线路参数为：电阻8.5，电感 0.064L，参数设置如下图：1.2.4 三相电压电流测量模块此模块将在变压器低压侧测量得到的电压、电流信号转变成 Simulink 信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签 “Vabc”，电流标签“Iabc”，参

4、数设置如下图：习Vo 1 t a 年 3 年 匹 玉 e-n t I 扯 -as e - t o - 罕 Cun d.F .u.s.e.a.1.a.b.e.1 %.,空 严 a 1 1 a b e 1( u s e a 瓦 立 3 b 拉 吐 t 七 已o11 e吐 5迂n a1)Va 阮厂 Vo :!.t a 巳ein pu习匹 工en t :i 年 匹 玉 en t I声Fus o a 1 ab乞罕 釭 吐 1abe1(us e a 氏 笠3 b 10 Bk t 七 已o11 e吐 s 迂n a1)二己 b 己1.2.5 故障设置模块勾选故障相A、B、C，设置短路电阻0.00001，设置

5、0.02s0.2s 发生短路故障，参数设置如下图：P aT a3 et 己工5F .P.h.a.5.e.F.a.u.1 t.气.,(7 Ph a s e B Faul t(7 Phas e C Fau1 tFau1 t r es 年 t an己 e sRon (o 血 s)IO厂 釭oun d Fau1 t厂 Ex t 釭 na1 con t工o1 f f au 1 t t i 口 i n巳 ：吐 a n 豆 t 还 n s t a t us 1 。 , 1 . . . ) :1 o百 an豆 t 迁 n t i口够 (s) :| O位 0 2 名l己n-Snubb e工5 工 仑5 坛 t

6、anc e 瓦 (o比35)习Snubb er s C忒 炟 豆 t an B e CD lF aT a d) inf妇 s 匹 玉 en t s I on e1.2.6 示波器模块为了得到仿真结果准确数值，可将示波器模块的“Data History”栏设置为下图所示：GeneralI Data h卤 ory ITip: try right clicking on axes厂 Limit data points to last: lsoooF Sto.avedataworkspace; Variable na me 庐习Format:Istructure with time1.3 无穷大功率电

7、源供电系统仿真结果及分析得到以上的电力系统参数后，可以首先计算出在变压器低压母线发生三相短路故障时短路电流周期分量幅值和冲击电流的大小，短路电流周期分量的幅值为Im=10.63kA，时间常数 Ta=0.0211s，则短路冲击电流为 Iim=17.3kA。通过模型窗口菜单中的“Simulation-Configuration Parameters”命令打开设置仿真参数的对话框，选择可变步长的 ode23t 算法，仿真起始时间设置为 0s，终止时间设置为 0.2s， 其他参数采用默认设置。在三相故障模块设置在0.02s 时刻变压器低压母线发生三相短路故障。运行仿真，得到变压器低压侧的三相短路电流波

8、形如下图所示：可见，短路电流周期分量的幅值为10.64kA，冲击电流为 17.39kA，与理论计算相比有差别，这是由于电源模块的内阻设置不同而造成的。实验二 同步发电机突然短路的暂态过程仿真2.1 发电机端突然发生三相短路的 Simulink 仿真模型构建根据给出的发动机参数，添加下列模块：（1）p.u.标准同步发电机模块（Synchronous Machine pu Standard）（2） 常数模块（Constant）（3） 电压测量模块（Voltage Measurement ）（4） 放大器模块（Gain）（5） 信号选择模块（Bus Selector）其他模块选取与前相同，建立 Si

9、mulink 仿真模型如下图所示：2.2 发电机端突然发生三相短路的 Simulink 仿真参数设置2.2.1 同步发电机模块设置同步发电机功率为 200MVA，电压13.8kV，频率50Hz，电抗设置如图，d 轴时间常数选择“Short-circuit”，q 轴时间常数选择“Open-circuit”。时间常数设置如图，定子电阻（p.u.） 0.005，惯性系数 3.2，摩擦系数 0，极对数 32，初始条件栏将由 Powergui 模块自动设定。参数设置如下图：Pra 已习e t 己工5氏 郅 e t = 华 1 I酝1J17 S.ho II det a.i :l. e dp a迂 立 e

10、t e工s 酝 立 type: I苹 1 i en t印 丘JX03沁 IIe 了 L一1 飞o1t . a n d f了已口 Pn ( VA)Vn (V立 3s )f n (Ez ) ：|2OO酕1 3B OO的 农e a Bt a n 它 es Xd Xd XdXC1. XC1. X1 （E心）：| 1 . 0 , 0 生 021;0 60 31二 0 . 1 5 习习d 釭 注 t i !:le 已 on s t a n t 5 : I三t - ,c:ii.工ICU:ii. t n 釭 坛 t i !:le con 尹 n t (s) : IO:De n - ,c:ii.工,c:u :i

11、i. t了i 已1e con s t a n t s 了d.| 1 轧 0 3 生 1 4 了d 了口廿 耘）：:St a t 七工 工 es i s t a n c: e 氐fo. u . ) :仑 e f f . o f i n er t 扫 心 fr 江 t i onf a 吐 凸ran d. 印 1 e p a i r 5 E ( 5)F ( 卯 ) p Q :| 3 2 0 32 -I n it . caon d. 扛 （） t h ( d. 尹：）扫 i b 江(pu)呻 已 扯 b如 己 ( d.尹：）Vf b u) ：54厂 S 立 ul a t e s a t ur a t

12、i on2.2.2 变压器模块设置功率 210MVA，频率50Hz，接线方式为yD1，二次侧中性点接地，绕组参数：一次侧电压 13.8kV，二次侧电压 230kV，电阻均为 0.0027，电感均为 0.08，励磁电阻 500，励磁电感 500。参数设置如下图：PaT已习 e t 己工5 Pn (VA) fn (Hz) 习酝 di n 巳 1 （ 虚 ） Bc nn e 吐 i on :I De l t a (D1)I芒 n di n巨p a丁 立 1e t 红 s V1 Ph - Ph (i 立 3s ) a 氐 环 L) a L1 b u ) 国 a e3 , oo 丐 n di n巳 2

13、( a阮） BEnn e 吐 i on :习L Y巳Wi n di n巳p aT 立 1e t 釭 s i 2 Ph - Ph ( i 立 35) 腔（即 ） L2 (l)l.J.) 厂 :Sa t u:r ab 乓 B廿工仑Y.a 即 e t i 廷 t i on 工 仑 s i s t an 晊 瓦 （ 印 ）55Y.a 即 e t i z a t i on 工廷 Bt a n B e L凹 b u )习踞 a s 匹 竺 e n t s I o n 己2.2.3 Powergui 模块初始化设置双击 Powergui 模块，打开潮流计算和电机初始化窗口，设定同步发电机为平衡节点“Swin

14、g bus”。初始化后，与同步发电机模块输入端口相连的两个常数模块Pm 和 Vf 以及同步发电机模块中的“Init.Cond.”将会自动设置。数据如下图所示：2.3 发电机端突然发生三相短路的 Simulink 仿真结果及分析选择 Ode15 算法，仿真结束时间取 1s。设置故障模块在 0.02025s1s 发生三相短路故障。开始仿真，得到发电机端突然三相短路后的三相定子电流波形图如下图所示：其中，A 相定子电流的冲击电流标幺值为9.1048，和理论计算值存在 0.95%的误差。短路后定子电流的d 轴和q 轴分量 Id、Iq 以及励磁电流 If 的仿真波形如下图所示：St 如 心“ 色rn d

15、 尸 P8 仁 I. i 、I I. I. I. I.I I I I, 牛II I, I I I I. I, I II，,.j 1, I, l, , ， 1111 I| I I, | l l , l ，l , U.,:,:I II,I, I，10001IIW.f,11从 11 ，|。30.40.59.0 70.862.HSt 如 C“ 色n, q 尸 P，旬, r畸 J： 看恤酝旷印因100011 吐 1 片 11,l|。3匮l0.50 70.869/ d pup.2.H扣 出 c 山 色” I过 1 |，I II I I, 1,.f, ，| I I II Ii II 11 11, I1,，:

16、 ,l , .，._，,_. , “,_| | | | Ill, I . o ，. J1111 III 1, I,1 1_,|匾il| 1 . 69.2.H固。 。0 11 吐 ll!l 廿1|。3匮l0.50 70.8现设置在 0.02025s 时发生 BC 两相短路故障。开始仿真，得到发电机端突然两相短路后的三相定子电流仿真波形如下图：实验三 小电流接地系统单相故障3.1 小电流接地系统仿真模型构建3.1.1 中性点不接地系统的仿真模型及计算利用 Simulink 建立一个 10kV 中性点不接地系统仿真模型，添加下列模块：（1） 输电线路模块（Three-Phase PI Section

17、 Line）（2） 信号接收模块（From）（3） 信号输出模块（Demux）（4） 输入加法器模块（Sum）（5） 三相序分量模块（Discrete 3-phase Sequence Analyzer）（6） 万用表模块（MultiMeter）建立中性点不接地系统仿真模型如下图所示：各模块参数设置如下：（1） 三相电源模块：电压 10.5kV，接线方式 Y 形连接，其他参数设置与实验一相同。（2） 输电线路模块 Line1Line4：线路长度分别为130km、175km、1km、150km，其他参数设置相同。下图为 Line1 参数设置。（3） 线路负荷模块：Load1Load3 设置其有功

18、负荷分别为 1MW、0.2MW、2MW，其它参数相同。Load4 设置为纯电阻负荷，有功负荷为1MW。下图分别为 Load1、Load4 参数设置。（4） 三相电压电流测量模块：勾选使用标签，按线路设置标签序号。下图为线路一测量模块的参数设置。（5） 故障模块设置：选择在第三条出线的1km 处（即 Line3 与 Line4 之间）发生A 相金属性单相接地，故障模块的参数设置如下图所示：系统的零序电压 3U0 及每条线路始端的零序电流3I0 采用下图连接方式测量得到：故障点的接地电流 Id 可用下图万用表测量得到：根据以上设置的参数，可以通过计算得到系统在第三条出线的1km 处（即 Line3

19、 与 Line4 之间）发生 A 相金属性单相接地时各线路始端的零序电流有效值为：线路 1：5.75A，线路 2： 13.5A，接地点的电流为 20.18A。3.1.2 中性点经消弧线圈接地系统的仿真模型及计算在上实验基础上，在电源中性点添加一个电感线圈，其他参数不变。仿真模型如下图所示：根据线路参数，如果要使接地点电流近似为0，计算得需要的补偿电感应为 L=0.9566H， 由于完全补偿存在串联谐振过电压问题，因此实际工程常采用过补偿方式，当取过补偿度为10%时，经计算消弧线圈电感为 L=0.8697H。消弧线圈参数设置如下图：3.2 小电流接地系统仿真结果及分析在仿真开始前，选择离散算法，

20、仿真结束时间取0.2s，利用 Powergui 模块设置离散方式，时间为 1x10-5s，系统在 0.04s 时发生A 相金属性单相接地。3.2.1 中性点不接地系统的仿真结果分析设置好参数，运行 10kV 中性点不接地仿真模型，得到系统三相对地电压和线电压的波形，如下图所示。从图中可见，系统在 0.04s 时发生A 相金属性单相接地后，A 相对地电压变为零，BC 相对地电压升高 3 倍，但线电压仍然保持对称故对负荷没有影响。系统的零序电压 3U0 及线路一始端的零序电流 3I0、故障点的接地电流 Id 波形如下图所示：仿真得到的各线路始端零序电流，接地电流Id 的有效值为：线路一：5.83A

21、，线路二：7.99A，线路三：13.86A，Id=20.64A。与理论计算值相比，仿真结果略大，但误差不大于 3%。从上图中可以看出，在中性点不接地方式下，非故障线路的零序电流超前零序电压90；故障线路的零序电流为全系统非故障元件对地电容电流之总和，零序电流滞后零序电压 90；故障线路的零序电流和非故障线路的零序电流相位差为180。故障后的零序分量还可以采用下图所示的“三相序分量模块”方法来得到，下图所示波形为故障线路零序电流的幅值和相位图。由图中可得故障线路零序电流的幅值为I0=6.52A，则3I0 的有效值为 13.83A，与从上图中线路三测量得到的 13.86A 仅相差 0.2%。3.2

22、.2 中性点经消弧线圈接地系统的仿真结果及分析设置好参数，运行 10kv 中性点经消弧线圈接地系统仿真模型，得到系统三相对地电压和线电压的波形与不接地系统仿真图相同。系统的零序电压 3U0 及线路一始端的零序电流 3I0、消弧线圈电流 IL、故障点的接地电流 Id 波形如下图所示：从上图所知，当单相接地故障的暂态过程结束后，故障点的接地电流 Id 的有效值在 2.9A 左右，远小于中性点不接地系统的接地电流，因此补偿的效果十分明显。对于非故障线路来说，其零序电流仍是本身的电容电流，零序电流超前零序电压90， 电容电流的实际方向为由母线流向线路，这与中性点不接地系统是相同的。但是对于故障路线路来

23、说，其零序电流将大于本身的电容电流，并且电容电流的实际方向也是由母线流向线路。因此，在这种情况下无法用电流方向的差别来判断故障线路，也很难用零序电流的大小来找出故障线路。实验四 Simulink 在变压器微机继电保护中的应用举例4.1 变压器仿真模型构建根据双侧电源的双绕组变压器的简单电力系统接线图，利用 Simulink 绘制仿真电路图如下：（1） 电源模块：电源 EM 与电源 EN 电势相位差 10，其他设置相同，下图为电源 EN 参数设置：（2） 变压器模块：勾选“饱和铁心（Saturable core）”，为了简化仿真，变压器两侧的绕组接线方式相同，电压等级也相同，参数设置如下图所示：

24、（3） 三相电压电流测量模块 UM、UN 将在变压器两侧测量到的电压、电流信号转变成Simulink 信号，相当于电压、电流互感器的作用。UM 模块的参数设置如下图所示，UN 模块参数设置与此相仿，只是输出信号分别为“Vabc_N”，“Iabc_N”。（4） 三相断路器模块 QF1 和 QF2 分别来控制变压器投入，故障模块 Fault1 和 Fault2 分别用来仿真变压器保护区内故障和区外故障。4.2 变压器空载合闸时励磁涌流的仿真设置三相断路器模块 QF1 的切换时间为 0s，仿真时间为 0.5s，仿真算法 Ode23t。三相断路器模块 QF2、故障模块 Fault1、Fault2 在仿

25、真中均不动作（设置其切换时间大于仿真时间即可）。为了观察励磁涌流，在仿真中添加下图所示示波器模块，参数设置如下图：运行仿真，得到空载合闸后的三相励磁涌流的波形如下图所示：通过 Powergui 模块中的 FFT Analysis 对励磁涌流波形进行谐波分析，其界面如下图所示：为了比较合闸时励磁涌流与短路电流的大小，设置故障模块 Fault1，使电路在 0.25-0.45s 间发生三相短路，运行仿真，结果如下图所示，在本次仿真中，A 相空载合闸时的励磁涌流峰值比短路电流要稍小，而B、C 相空载合闸时的励磁涌流峰值要比短路电流大。将变压器的二次绕组改为“D11”接线方式时，电源 EM 的 A 相初

26、相位仍设为 0，运行仿真，得到空载合闸后的三相励磁涌流的波形如下图所示：4.3 变压器保护区内、外故障时比率制动的仿真为了仿真比率制动式差动保护在变压器保护区内、外故障时电流的情况，增加运算及示波器模块如下图所示：设置三相断路器模块 QF1、QF2 的切换时间均为 0s，并设置故障模块Fault1，使电路在0.3-0.5s 间发生三相短路，故障模块 Fault2 不动作，运行仿真，得变压器保护区内故障时的电流波形如下图所示，从图中可以明显看出差动电流大于制动电流，保护能够可靠动作。设置故障 Fault2，使电路在 0.3-0.5s 间发生三相短路，故障模块 Fault1 不动作，运行仿真，得变压器保护区外故障的电流波形如下图所示，从图中可以明显看出制动电流大于差动电流，保护制动，不动作。4.4添加模块饱和变压器（Saturable Transformer），勾选“三绕组变压器（Three windings Transformer）”，构造出具有一个初级绕组、两个次级绕组的单相变压器。改造后的仿真模型如下图所示：设置两个次级绕组参数相同，并设置三相断路器模块 QF1、QF2 的切换时间均为 0s，设置故障模块 Fault1 使电路在 0.3-0.5s 间发生 AB 相间短路，故障模块 Fault2 不动作，运行仿真，得到变压器绕组 50%处发生两相短路故障时的电流波形如下图所示：