南昌大学电力系统分析实验报告.doc

《南昌大学电力系统分析实验报告.doc》由会员分享，可在线阅读，更多相关《南昌大学电力系统分析实验报告.doc（17页珍藏版）》请在三一文库上搜索。

1、南昌大学实验报告实验类型： 验证 综合 设计 创新 实验日期： 实验成绩： 一、实验项目名称 电力网数学模型模拟实验二、实验目的与要求：本实验通过对电力网数学模型形成的计算机程序的编制与调试，获得形成电力网数学模型：节点导纳矩阵的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩阵。通过实验教学加深学生对电力网数学模型概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问

2、题结合起来。三、主要仪器设备及耗材计算机、软件（已安装，包括各类编程软件C语言、C+、VB、VC等、应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等）四、实验步骤1、 将事先编制好的形成电力网数学模型的计算程序原代码由自备移动存储设备导入计算机。（1） 编程思想 I 无变压器支路的节点导纳矩阵计算方法。以下语句用于输入题目已知节点之间的导纳和阻抗值：z=input(请输入由节点号对应的阻抗形成的矩阵：z=);y=input(请输入由节点号对应导纳形成的矩阵：y=);%其中，yii即为i节点对地导纳yii0，两节点之间无直接相连的通路则yij输入为0，zij输入为inf（即无穷大）

3、，zii也输入为inf。以下语句用于计算无变压器支路的节点导纳矩阵：for(i=1:n)for(j=1:n)Y(i,i)=sum(y(i,:),2)+sum(1./z(i,:),2);if j=iY(i,j)=Y(i,i);elseY(i,j)=-1/z(i,j); endendend%其中，对角线元素Yii=yii0+jyij=yii0+j1zij非对角线元素Yij=-yik II变压器支路的等值电路1:k图1 变压器支路的型等值电路kzk*kz/(1-k)kz/(k-1)jjizi当节点a，b间接有变压器支路时（见图1），当然可以用型等值电路，然后按照上述原则形成导纳矩阵。但在实际应用程序

4、中，往往直接计算变压器支路对导纳矩阵的影响。根据图1-5-1可以写出节点a，b的自导纳和节点间的互导纳增量分别如下：节点a的自导纳改变量式（I-1）： Yii=1kz+k-1kz=1z （I-1） 节点j的自导纳改变量式（I-2）： Yjj=1kz+1-kk2z=1k2z （I-2） 增加节点i，j间的互导纳式（I-3）： Yij=Yji=-1kz （I-3） 注意：在输入电路已知参数时，将k以矩阵形式输入，若i、j两节点之间无变压器，则kij=1,（kii=1），若i、j之间有变压器则将靠近变压器的一端记为节点j，靠近阻抗的一端记为节点i，输入人员将kij输入为变压器变比，而kji以1输入。

5、程序段中for(i=1:n)for(j=1:n)if k(i,j)=1 y(j,i)=(1-k(i,j)/(k(i,j)*k(i,j)*z(i,j); y(i,j)=(k(i,j)-1)/(k(i,j)*z(i,j); z(i,j)=k(i,j)*z(i,j); z(j,i)=z(i,j);endendend%判断两个节点之间是否有变压器，如果k1则对节点之间的导纳和阻抗矩阵进行修改，得到新的yij和zij（2）源代码如下：clcclear n=input(请输入节点数：n=); z=input(请输入由节点号对应的阻抗形成的矩阵：z=);y=input(请输入由节点号对应导纳形成的矩阵：y=

6、);k=input(请输入变比：k=);Y=zeros(n);for(i=1:n)for(j=1:n)if k(i,j)=1 y(j,i)=(1-k(i,j)/(k(i,j)*k(i,j)*z(i,j); y(i,j)=(k(i,j)-1)/(k(i,j)*z(i,j); z(i,j)=k(i,j)*z(i,j); z(j,i)=z(i,j);endendendfor(i=1:n)for(j=1:n)Y(i,i)=sum(y(i,:),2)+sum(1./z(i,:),2);if j=iY(i,j)=Y(i,i);elseY(i,j)=-1/z(i,j); endendend2、 在相应的编程

7、环境下对程序进行组织调试。3、 应用计算例题验证程序的计算效果。将例题中的已知条件如下依次输入:n=5z=inf 0.105i inf inf inf;0.105i inf 0.024+0.065i 0.03+0.08i inf;inf 0.024+0.065i inf 0.018+0.05i inf;inf 0.03+0.08i 0.018+0.05i inf 0.184i;inf inf inf 0.184i infy=0 0 0 0 0 ;0 0 0+0.0160i 0+0.0200i 0;0 0+0.0160i 0 0+0.0130i 0;0 0+0.0200i 0+0.0130i 0

8、 0;0 0 0 0 0k=1.0000 1.0500 1.0000 1.0000 1.0000;1.0000 1.0000 1.0000 1.0000 1.0000;1.0000 1.0000 1.0000 1.0000 1.0000;1.0000 1.0000 1.0000 1.0000 1.0000;1.0000 1.0000 1.0000 0.9600 1.00004、 对调试正确的计算程序进行存储、打印。5、 完成本次实验的实验报告。五、实验数据及处理结果 运行自行设计的程序，把结果与手工计算结果相比较，验证所采用方法及所编制程序运行的正确性。1手算节点导纳矩阵：详细计算过程见课本P

9、73例4-1其结果为：0.0000-j9.52380.0000+j9.07030000.0000+j9.07039.1085-i33.1002j13.5388-4.9989j10.9589-4.106900j13.5388-4.998911.3728-j31.2151j17.7053-6.373900j10.9589-4.1069j17.7053-6.373910.4835-j34.52830.0000+j5.66120000.0000+j5.66120.0000-j5.43482用matlab程序计算节点导纳矩阵实验结果： 对比可知通过比较，发现得到的节点导纳矩阵中大部分的数据是相同的，只是

10、少部分的结果有微小误差。这也就说明了手算能得到同样的结果，但我们发现用一般手算计算节点导纳矩阵比用计算机编程计算节点导纳矩阵繁琐了很多，而一个系统的节点数不可能是几个，而是由许多的节点构成的，若在此时用一般手算计算节点导纳矩阵，这就显得非常复杂，在实际操作过程当中也没有那么多的时间去让我们去用手算计算节点导纳矩阵。所以在电力系统中我们往往都是采用计算机编程计算节点导纳矩阵，这样大大提高了工作效率。六、思考讨论题或体会或对改进实验的建议1什么是输入阻抗？什么是转移阻抗？网络化简的方法有哪些？答：输入阻抗：在节点i单独注入电流而所有其他节点的注入电流都等于0时，在节点i产生的电压同注入电流之比，即

11、等于节点i的自阻抗，即输入阻抗。转移阻抗：节点i单独施加Ei时，该电势与其在节点j产生的短路电流Ij之比即等于节点i、j之间的转移阻抗。网络化简方法：星网变换法，高斯消去法。2.简述节点导纳矩阵的形成的过程，节点导纳矩阵的阶数与系统的节点数有什么关系？节点导纳矩阵的互导纳yij在数值上等于什么？举例说明当网络结构发生变化时，如何修改？答：节点导纳矩阵阶数等于系统的节点数； 节点导纳矩阵的互导纳yij在数值上等于联接节点i、j的支路导纳的负值；节点导纳矩阵的修改（1） 从网络的原有节点i引出一条导纳为yik 的支路，同时增加一个节点k，由于节点数加1，导纳矩阵将增加一行一列。新增的对角线元素Yk

12、k=yik。新增的非对角线元素中，只有Yik =Yki = -yik，其余的元素都为矩阵的原有部分，只有节点i的自导纳应增加Yii=yik。（2） 在网络的原有节点i、j之间增加一条导纳为yij的支路。由于只增加支路，不增加节点，故导纳矩阵的阶次不变。因而只要对与节点i、j有关的元素分别增添以下的修改增量即可：Yii=Yjj=yij ， Yij=Yji=-yij其余元素都不用修改。（3） 在网络的原有节点i、j之间切除一条导纳为yij的支路，这种情况可以当做是在i、j节点间增加一条导纳为-yij的支路来处理，因此，导纳矩阵中有关元素的修正增量为Yii=Yjj=-yij ， Yij=Yji=yi

13、j举例：如将例4-1中节点4、5间变压器变比由=0.96改为=0.98后作节点导纳矩阵的修改：变比修改后，相当于先切除变比为的变压器，再接入变比为的变压器。与之相关的4、5节点的导纳值要发生变化，而网络的其他元素不发生变化。利用例4-1中结果，与节点4、5有关的导纳矩阵元素的修正增量应为3.节点电压方程的求解方法有哪些？答：直接求解法和高斯消去法4.简述节点阻抗矩阵元素的物理意义及其形成。答：节点导纳矩阵中，Yik =IiVk|Vj=0,jk当k=i时，上述公式说明，当网络中除节点i以外，所有节点都接地时，从节点i注入网络的电流同施加于节点i的电压之比，即等于节点i的自导纳Yii。换句话说，自

14、导纳Yii是节点i以外所有节点都接地时节点i对地的总导纳。当ki时，公式说明，当网络中除节点k以外所有接地都接地时，从节点i流入网络的电流同施加于节点k的电压之比，即等于节点k、i之间的互导纳Yik 。在这种情况下，节点i的电流实际上是自网络流出并流入地中的电流，所以Yik 应等于节点k、i之间的支路导纳的负值。七、参考资料1. 电力系统分析何仰赞 华中科技大学出版社2. 电力系统稳态分析陈珩 中国电力出版社3. 电力系统暂态分析李光琦 中国电力出版社4. 电力系统计算 水利电力出版社南昌大学实验报告学生姓名： 学 号： 专业班级： 实验类型： 验证 综合 设计 创新 实验日期： 11.30

15、实验成绩： 一、实验项目名称 电力网数学模型模拟实验二、实验目的与要求：本实验通过对电力网数学模型形成的计算机程序的编制与调试，获得形成电力网数学模型：节点导纳矩阵的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩阵。通过实验教学加深学生对电力网数学模型概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。三、主要仪器设备及耗材计算机、软件（已安装，包括各类

16、编程软件C语言、C+、VB、VC等、应用软件MATLAB等）、移动存储设备（学生自备，软盘、U盘等）四、实验步骤6、 将事先编制好的形成电力网数学模型的计算程序原代码由自备移动存储设备导入计算机。7、 在相应的编程环境下对程序进行组织调试。8、 应用计算例题验证程序的计算效果。9、 对调试正确的计算程序进行存储、打印。10、 完成本次实验的实验报告。五、实验数据及处理结果 运行自行设计的程序，把结果与手工计算结果相比较，验证所采用方法及所编制程序运行的正确性。实验程序：clear clc;z=inf,0.10+0.40i,0.3i,0.12+0.50i;0.10+0.40i,inf,inf,0

17、.08+0.40j;0.3i,inf,inf,inf;0.12+0.50i,0.08+0.40i,inf,inf; （各支路的阻抗）y=0,0.01528i,0,0.01920j;0.01528i,0,0,0.01413j;0,0,0,0;0.01920i,0.01413i,0,0; （各支路的导纳）Y=zeros(4,4);for(i=1:4), for(j=1:4), if i=j Y(i,j)=Y(i,j) else Y(i,j)=-1.0/z(i,j) end endend for (i=1:4), for(j=1:4), Y(i,i)=Y(i,i)+y(i,j)+1.0/z(i,j)

18、 end endY(1,1)=Y(1,1)-1.0/z(3,1)+1.1*1.1/z(3,1) （对变压器支路的两个节点进行修正）Y(1,3)=Y(1,3)*1.1Y(3,1)=Y(3,1)*1.1运行结果，实验结果导纳截图：手工计算结果：Y11=y140+y120+1/z12+1/z14+1.1*1.1/z13=j0.01920+j0.01528+1/(0.10+j0.40)+ 1/(0.12+j0.50)+ 1.1*1.1/j0.3=1.0241-j8.2429Y12=Y21=-1/z12=1/(0.10+j0.40)=-0.5882+j2.3529Y13=Y31=-1.1/z13=-1.

19、1/j0.3=j3.667Y14=Y41=-1/z14=1/(0.12+j0.50)=-0.4539+j1.8911Y22=y210+y240+1/z21+1/z24=j0.01528+j0.01413+1/(0.10+j0.40)+1/(0.08+j0.40)=1.069-j4.7274Y24=Y42=-1/z24=-1/(0.08+j0.40)=-0.4808+j2.4038Y33=1/z13=1/(j0.3)=-j3.333Y44=y410+y420+1/z14+1/z42=j0.01920+j0.01413+1/(0.12+j0.50)+1/(0.08+j0.40)=0.9346-j4

20、.2616六、思考讨论题或体会或对改进实验的建议1什么是输入阻抗？什么是转移阻抗？网络化简的方法有哪些？答: 输入阻抗就是指节点的自阻抗即节点i单独注入电流时，在节点i产生电压与电流之比；转移阻抗就是指节点i单独注入电势，与在短路点j产生的电流之比；网络化简方法有高斯消去法：也就是带有节点电流移置的星网变换； 2.简述节点导纳矩阵的形成的过程，节点导纳矩阵的阶数与系统的节点数有什么关系？节点导纳矩阵的互导纳yij在数值上等于什么？举例说明当网络结构发生变化时，如何修改？答：（1）节点导纳矩阵的形成过程：列写节点电压方程：以零电位点为计算节点电压参考点，根据基尔霍夫电流定律，写出各独立节点电流平

21、衡方程，经整理后形成矩阵方程：。矩阵Y就是所求节点导纳矩阵。（2）节点导纳矩阵阶数等于系统的节点数（注：不包括接地点）。（3）节点导纳矩阵的互导纳数值上等于节点i、j间支路导纳的负值，即。（4）对于网络结构发生的变化主要有这几种形式：从网络的原有节点i引出一条支路；在网络的原有节点ij之间增加一条支路；在网络的原有节点之间切除一条支路。下面如将例中节点1、3间变压器变比由K13=1.1改为K=1后作节点导纳矩阵的修改：变比修改后，相当于先切除变比为1.1的变压器，再接入变比为1的变压器。与之相关的1、3节点的导纳值Y33、Y11、Y31、Y13要发生变化，而网络的其他元素不发生变化。利用上述的

22、计算结果Y11=Y11-1.1*1.1/+1/z13=1.0241-j8.2429-1.1*1.1/j0.3+1/j0.3=1.0241-j7.5429Y33= Y33= =1/(j0.3)=-j3.333Y13=Y31= -1/z13=-j3.3333.节点电压方程的求解方法有哪些？答;节点电压方程的求解方法有：直接求解法和高斯消去法4.简述节点阻抗矩阵元素的物理意义及其形成。答：节点阻抗矩阵元素的物理意义：在节点k单独注入电流，在节点k产生的电压同注入电流之比，即等于节点k的自阻抗;在节点i产生的电压同节点k注入电流之比，即等于节点k和节点i之间的互阻抗。节点阻抗矩阵的形成：将节点电压方程

23、：改写成,式中就是阻抗矩阵七、实验体会通过本次实验，对节点导纳矩阵的求解有了更深的认识，特别是怎么求解一个网络的支路节点导纳和有变压器的支路的导纳有了更深的了解。为以后在求解类似的问题中有了一定的基础。通过实验呢更是加深了对求解的过程，应先怎么办，再怎么做有了了解。还有在本次实验中对于MATLAB的使用有了一定的基础，虽然时间很短，但在实验之后对于用MATLAB软件来求解电力系统中问题总是有一定的方便性认识更深。八、参考资料5. 电力系统分析何仰赞 华中科技大学出版社6. 电力系统稳态分析陈珩 中国电力出版社7. 电力系统暂态分析李光琦 中国电力出版社8. 电力系统计算 水利电力出版社实验素材

24、：南昌大学实验报告学生姓名： 学号： 专业班级： 实验类型： 验证 综合 设计 创新 实验日期： 12.7 实验成绩： 一、实验项目名称 电力系统短路计算实验二、实验目的与要求：目的：通过实验教学加深学生的基本概念，掌握电力系统的特点，使学生通过系统进行物理模拟和数学模拟，对系统进行电力系统计算和仿真实验，以达到理论联系实际的效果。通过电子计算机对电力系统短路等计算的数学模拟，分析电力系统的故障计算方法、实现工程计算的功能。提高处理电力系统工程计算问题的实际能力，以及实现对电力系统仿真的过程分析。要求：l、 使学生掌握对电力系统进行计算、仿真试验的方法，了解实验对电力系统分析研究的必要性和意义

25、。2、使学生掌握使用实验设备计算机和相关计算软件、编程语言。3、应用电子计算机完成电力系统的短路计算。4、应用电子计算机及相关软件对电力系统进行仿真。 三、主要仪器设备及耗材1.每组计算机1台、相关计算软件1套四、实验步骤1. 将事先编制好的形成电力网数学模型的计算程序原代码由自备移动存储设备导入计算机。2. 在相应的编程环境下对程序进行组织调试。3. 应用计算例题验证程序的计算效果。4. 对调试正确的计算程序进行存储、打印。5. 完成本次实验的实验报告。五、实验数据及处理结果运行自行设计的程序，把结果与例题的计算结果相比较，验证所采用的短路电流计算方法及程序运行的正确性。如果采用的是近似计算

26、方法，还需分析由于近似所产生的误差是否在运行范围内。实验程序：clear clc;z=0.2i,inf,0.51i,inf;inf,4i,0.59i,inf;0.51i,0.59i,inf,1.43i;inf,inf,1.43i,inf;y=0,0,0,0;0,0,0,0;0,0,0,0;0,0,0,0;f=4;Y=zeros(4,4);for(i=1:4), for(j=1:4), if i=j Y(i,j)=Y(i,j) else Y(i,j)=-1.0/z(i,j) end endend for (i=1:4), for(j=1:4), Y(i,i)=Y(i,i)+y(i,j)+1.0/

27、z(i,j) end endZ=inv(Y);If=1/Z(f,f);实验结果：If=0-0.48902i 实验例题所给结果短路电流：If = - j0.4895，与程序运行结果在误差允许范围之内，故验证了该程序的正确性。六、思考讨论题或体会或对改进实验的建议1. 理解课本上讲述的同步电机突然三相短路的物理分析。答：同步电机稳态对称运行（包括稳态对称短路）时，电枢磁势大小不随时间变化，而在空间以同步速旋转，同转子没有相对运动，故不会在转子绕组中感应电流。突然短路时，定子电流发生急剧的变化，电枢反应磁通也随着变化，定转子间电流会相互影响，这是同步电机突然短路暂态过程区别于稳态短路的显著特点。我们

28、在进行磁感分析时，对于每个绕组都遵守磁链守恒原则。对于无阻尼绕组同步电机突然三相短路时，短路后定子侧将出现：基频电流，由三相对称绕组的基频电流产生的交变磁链，用以抵消转子主磁场对定子各相绕组产生的交变磁链；直流：对各绕组产生的不变、，来维持定子绕组的磁链初值不变；倍频电流：定子各相直流产生的恒定磁势，当转子旋转时，因转子纵横轴磁阻不同，转子每转过，磁阻经历一个变化周期，为适应磁阻的变化，产生倍频电流与直流共同作用，才能维持定子侧磁链初值不变。转子侧产生：附加直流分量：为抵消定子电流对转子产生的强烈电枢反应影响，维持磁链不变，该附加直流与原直流同向，加强了励磁绕组的磁场，而且产生磁通的一部分也要

29、穿过定子绕组，故激起定子基频电流大大超过稳态短路电流转子基频交流：为抵消定子的直流和倍频电流产生的电枢反应，该基频电流在转子产生脉振磁场，分解为正反两方向磁场，又来影响定子侧的以及在定子侧产生两倍频的交变磁链，定子侧的就是为了抵消该磁链而产生的。2. 简述无阻尼绕组同步电机突然三相短路时，短路电流所含各种分量以及各自由电流衰减时间常数的确定，有阻尼绕组同步电机突然三相短路时有什么不同之处？答：（1）无阻尼绕组同步电机突然三相短路时，定子侧短路电流有：基频电流：抵消转子主磁场对定子各相绕组产生的交变磁链；直流：维持定子绕组的磁链初值不变。倍频电流：为适应磁阻的变化，倍频电流与直流共同作用，才能维

30、持定子侧磁链初值不变。转子侧短路电流有：附加直流分量：为抵消定子电流对转子产生的强烈电枢反应影响，维持磁链不变；转子基频交流：为抵消定子的直流和倍频电流产生的电枢反应。（2）在实用计算中，各自由电流衰减常数确定，常采用以下的简化原则：为维持磁链守恒而在短路瞬间出现的自由电流，若它产生的磁通对于本绕组相对静止，那该电流将按本绕组时间常数衰减，一切同改自由电流发生依存关系的均按同一时间常数衰减；某绕组时间常数是该绕组电感同电阻之比,而忽略其他绕组电阻的影响。（3）有阻尼绕组同步电机突然三相短路时的不同：有阻尼绕组电机，在转子纵轴向有励磁绕组和阻尼绕组，横轴向也有阻尼绕组。突然短路时，定子基频电流突然增大，电枢反应磁通也突然增加，励磁绕组、阻尼绕组为维持磁链守恒，都会感应出自由直流，由此来抵消电枢反应磁通的增加。转子各绕组自由直流产生磁通一部分又进入定子侧，由此定子侧会有基频电流自由分量。注意，转子纵轴向的励磁绕组和阻尼绕组间存在互感关系，突然短路瞬间它们当中任一绕组的磁链守恒都是靠两绕组的自有电流共同维持的。