潮流上机实验报告.docx
《潮流上机实验报告.docx》由会员分享,可在线阅读,更多相关《潮流上机实验报告.docx(32页珍藏版)》请在三一办公上搜索。
1、潮流上机实验报告电力系统潮流上机课程设计报告 院系 班级: 学号: 学生姓名: 指导教师: 设计周数 成绩: 日期:年月日课程课程设计报告 一、课程设计的目的与要求 培养学生的电力系统潮流计算机编程能力,掌握计算机潮流计算的相关知识 二、设计正文 1.手算: 要求应用牛顿-拉夫逊法或P-Q分解法手算求解,要求精度为0.001MW。 &=1.00,节点2为PQ节点,负荷功率S=0.8+j0.6,节节点1为平衡节点,电压U12点3是PV节点,P3=0.4,U3=1.1,两条支路分别为Z13=0.01+j0.04,Z12=0.05+j0.2,对地支路y30=j0.33。 1 课程课程设计报告 2.计
2、算机计算:编写潮流计算程序,要求如下: 2.1据给定的潮流计算任务书整理潮流计算的基础数据:节点的分类,线路模型,等值变压器模型,电压等级的归算,标幺值的计算; 2.2基础数据的计算机存储:节点数据,支路数据; 2.3用牛顿-拉夫逊法计算; 2.4根据所选潮流计算方法画流程图,划分出功能模块,有数据输入模块,导纳阵形成模块,解线性方程组模块,计算不平衡功率模块,形成雅可比矩阵模块,解修正方程模块,计算线路潮流,网损,PV节点无功功率和平衡节点功率,数据输出模块; 2.5据上述模块编制程序并上机调试程序,得出潮流计算结果; 2.6源程序及其程序中的符号说明集、程序流图 简单系统如下图所示,支路数
3、据如下: z12=0.1+j0.41,z13=j0.3,z14=0.12+j0.5,z24=0.08+j0.40 ,y10,4=y40,1=j0.0192,y20,4=y40,2=j0.01413 y10,2=y20,1=j0.01528k=1.1 节点数据如下: &=-0.30-j0.18,S&=-0.55-j0.13, S21&=0.5,U3=1.10,U4=1.050o S31)节点导纳阵 #include #include #include #include LF.h /form node conductance matrix intMakeY( intnB, intnL, Line*
4、 sL, double* YG, double* YB ) inti,j,l; double r,x,d1,g,b,t; for(i=0;inB;i+) for(j=0;jnB;j+) YGij=0.0; YBij=0.0; for(i=0;inL;i+) r=sLi.R; x=sLi.X; g=r/(r*r+x*x); b=-x/(r*r+x*x); switch(sLi.Type) case 1:/Line 2 课程课程设计报告 break; case 2:/Transformer b*=1/sLi.K; break; g*=1/sLi.K; YGsLi.NumIsLi.NumI+=g;
5、YGsLi.NumJsLi.NumJ+=g; YGsLi.NumIsLi.NumJ-=g; YGsLi.NumJsLi.NumI-=g; YBsLi.NumIsLi.NumI+=b+sLi.B; YBsLi.NumIsLi.NumJ-=b; YBsLi.NumJsLi.NumI-=b; YBsLi.NumJsLi.NumJ+=b+sLi.B; printf(实部:n); for(i=0;inB;i+) for(j=0;jnB;j+) printf(%lft,YGij); printf(n); for(i=0;inB;i+) for(j=0;jnB;j+) printf(%lft,YBij);
6、printf(n); printf(虚部:n); /* Check the Y matrix */ ofstreamfout(out.txt); fout -Y Matrix- endl; for(i=0;inB;i+) for(j=0;jnB;j+) fout YGij +j YBij t; foutendl; fout.close; return 0; 2)计算功率不平衡量 #include #include #include #include LF.h /form delta p and delta q intCalDeltaPQ( intnpv, intnpq, Bus* bus, d
7、ouble* YG, double* YB, int* p_Jtobus, double* deltaf ) intk,i,j; for(k=0;knpv+npq*2;k+) i=p_Jtobusk; if(knpv) 3 课程课程设计报告 deltafk=busi.GenP-busi.LoadP; for(j=0;jnpv+npq+1;j+) deltafk-=busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); printf(PV节点%d的有功功率是%lfn,i,delta
8、fk); if(k=npv) deltafk=busi.GenP-busi.LoadP; for(j=0;jnpv+npq+1;j+) deltafk-=busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); printf(PQ节点%d的有功功率是%lfn,i,deltafk); if(k=npv+npq) deltafk=busi.GenQ-busi.LoadQ; for(j=0;jnpv+npq+1;j+) deltafk-=busi.Volt*busj.Volt*(YGij
9、*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); 3)雅各比矩阵的计算 /* Purpose: for undergraduate course Task: Load Flow Copyright NCEPU, Liu Chongru */ #include #include #include #include LF.h /form Jacobian matrix intFormJacobian( intnpv, intnpq, Bus* bus, double* YG, double* YB, int* p_Jtobus
10、, double* Jac ) intnp = npv+npq,j,k,i,m; /TODO for(k=0;knpv+npq*2;k+) i=p_Jtobusk; qi=0; 4 printf(PQ节点%d的无功功率是 %lfn,i,deltafk); return 0; double a14,q14; ai=0; 课程课程设计报告 if(knp)/H N if(j!=i) for(j=0;jnp+1;j+) ai+=busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); qi+=busj.Vol
11、t*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); j=p_Jtobusm; if(j!=i) for(m=0;mnpv+npq*2;m+) if(mnp) Jackm=busi.Volt*busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase);/Form H else Jackm=busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Ph
12、ase-busj.Phase);/Form N else for(j=0;jnp+1;j+) if(mnp) Jackm=-busi.Volt*ai;/Form H else if(j=i) else Jackm=busi.Volt*qi+2*busi.Volt*busi.Volt*YGij;/Form N if(j!=i) ai+=busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); qi+=busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(
13、busi.Phase-busj.Phase); j=p_Jtobusm; if(j!=i) for(m=0;mnpv+npq*2;m+) if(mnp) Jackm=-busi.Volt*busj.Volt*(YGij*cos(busi.Phase-busj.Phase)+YBij*sin(busi.Phase-busj.Phase); /Form J else Jackm=busi.Volt*busj.Volt*(YGij*sin(busi.Phase-busj.Phase)-YBij*cos(busi.Phase-busj.Phase); /Form L else if(j=i) if(m
14、np) 5 课程课程设计报告 Jackm=busi.Volt*qi; Jackm=busi.Volt*ai-2*busi.Volt*busi.Volt*YBij; else for(i=0;inp+npq;i+) for(int j=0;jnp+npq;j+) printf(%d %d %f ,i,j,Jacij); printf(n); /Output the matrix to check the Jacobian matrix ofstreamfout(out.txt,ios:app); fout -Jacobian Matrix- endl; for(i=0; inp+npq;i+ )
15、 for(j=0; jnp+npq; j+ ) foutJacij t; foutendl; fout.close; return 0; 4)线路损耗 /8.calculate the power flow double* p_Pij, *p_Qij, *p_Pji, *p_Qji; p_Pij = new doublenL; p_Qij = new doublenL; p_Pji = new doublenL; p_Qji = new doublenL; int x1,x2; for( i=0; inL; i+ ) x1=linei.NumI; x2=linei.NumJ; if(linei
16、.Type=1) p_Piji=busx1.Volt*busx1.Volt*(-YGx1x2)-busx1.Volt*busx2.Volt*(-YGx1x2)*cos(busx1.Phase-busx2.Phase)+(-YBx1x2)*sin(busx1.Phase-busx2.Phase); p_Qiji=-busx1.Volt*busx1.Volt*(linei.B+(-YBx1x2)-busx1.Volt*busx2.Volt*(-YGx1x2)*sin(busx1.Phase-busx2.Phase)-(-YBx1x2)*cos(busx1.Phase-busx2.Phase); p
17、_Pjii=busx2.Volt*busx2.Volt*(-YGx2x1)-busx2.Volt*busx1.Volt*(-YGx2x1)*cos(busx2.Phase-busx1.Phase)+(-YBx2x1)*sin(busx2.Phase-busx1.Phase); p_Qjii=-busx2.Volt*busx2.Volt*(linei.B+(-YBx2x1)-busx2.Volt*busx1.Volt*(-YGx2x1)*sin(busx2.Phase-busx1.Phase)-(-YBx2x1)*cos(busx2.Phase-busx1.Phase); else p_Piji
18、=busx1.Volt*busx1.Volt*(-YGx1x2)/linei.K-busx1.Volt*busx2.Volt*(-YGx1x2)*cos(busx1.Phase-busx2.Phase)+(-YBx1x2)*sin(busx1.Phase-busx2.Phase); p_Qiji=-busx1.Volt*busx1.Volt*(-YBx1x2)/linei.K+linei.B)-busx1.Volt*bu6 课程课程设计报告 sx2.Volt*(-YGx1x2)*sin(busx1.Phase-busx2.Phase)-(-YBx1x2)*cos(busx1.Phase-bus
19、x2.Phase); p_Pjii=busx2.Volt*busx2.Volt*(-YGx2x1*linei.K)-busx2.Volt*busx1.Volt*(-YGx2x1)*cos(busx2.Phase-busx1.Phase)+(-YBx2x1)*sin(busx2.Phase-busx1.Phase); p_Qjii=-busx2.Volt*busx2.Volt*(-YBx2x1)*linei.K+linei.B)-busx2.Volt*busx1.Volt*(-YGx2x1)*sin(busx2.Phase-busx1.Phase)-(-YBx2x1)*cos(busx2.Pha
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 潮流 上机 实验 报告
链接地址:https://www.31ppt.com/p-3639163.html