预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10
亲,该文档总共48页,到这已经超出免费预览范围,如果喜欢就直接下载吧~
如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
编号:时间:2021年x月x日书山有路勤为径,学海无涯苦作舟页码:第PAGE48页共NUMPAGES48页第PAGE\*MERGEFORMAT48页共NUMPAGES\*MERGEFORMAT48页课程设计报告书题目:电力系统稳定分析和计算学院专业学生姓名学生学号指导教师课程编号课程学分起始日期封面纸推荐用210g/m2的绿色色书编辑完后需将全文绿色说明文字删除,格式不变C2课程设计题目:电力系统稳定分析和计算姓名:指导教师:一个220kV分网结构和参数如下:500kV站(#1)的220kV母线视为无穷大母线,电压恒定在230kV。图中,各变电站负荷均等效接于220kV母线。变电站参数和负荷如下表:编号类型220kV最大负荷,MVA#1500kV站平衡节点#2220kV站200+j110#3220kV站380+j140#4220kV站240+j130#5220kV站300+j150#6220kV站190+j80各线路长度如图所示。所有线路型号均为LGJ-2*300,基本电气参数为:正序参数:r=0.054Ω/km,x=0.308Ω/km,C=0.0116µF/km;零序参数:r0=0.204Ω/km,x0=0.968Ω/km,C0=0.0078µF/km;40ºC长期运行允许的最大电流:1190A。燃煤发电厂G有三台机组,发电机与升压变之间采用单元接线。电厂220kV侧采用单母分段接线,正常运行时分段开关闭合。发电机组主要参数如下表:机组台数单台容量(MW)额定电压(EV)功率因数升压变容量MVA升压变Vs%XdXd’XqTd0’TJ=2H330010.50.8535010.52.00.231.7415当发电机采用三阶模型时,励磁环节(含励磁机和励磁调节器)模型如下(不考虑PSS):上图中参数如下:TR=0,KA=25,TA=0.1,Te=0.15,KE=1,KF=0.05,TF=0.7发电厂升压变均采用Y/Y0接线,变比10.5kV/242kV。不计内阻和空载损耗。发电厂按PV方式运行,高压母线电压定值为1.05VN。发电厂厂用电均按出力的7%考虑。稳定仿真中不考虑发电厂的调速器和原动机模型。负荷采用恒阻抗模型,负序阻抗标幺制取0.2。设计的主要内容:1、手工进行参数计算和标幺化,形成潮流计算的网架参数;2、用Matlab编制潮流计算程序,要求采用P-Q分解潮流计算方法。3.考虑该电厂开机三台,均为额定容量输出的方式,用编制的程序计算潮流。3、用PowerWorld软件进行潮流计算并与自己编制的软件计算结果进行校核和分析;4、设#2和#4母线之间双回线路中一回的中点分别发生以下2种故障:1)1s时发生三相短路,1.1s同时切除故障线路三相;2)1s时单相接地短路,1.1s时同时切除故障线路三相。1.9s时三相重合闸。因重合于永久性故障,2.1s时再次切除故障线路。试手工计算序网以及用网络变换法求解转移阻抗;5、针对问题4,用Matlab编制稳定计算程序(三台机可并联等值成一台机),发电机采用二阶经典模型(注:用ode45函数既可求解),用软件计算出摇摆曲线,要求输出发电机功角,角速度。6、用PowerWorld软件的分析问题4,并与编程计算结果进行比较校核。7、发电机采用三阶模型,用PowerWorld作为分析工具,对问题4的两种故障方案进行稳定计算,给出摇摆曲线,并计算故障的极限切除时间。8、比较两种模型的仿真结果,分析发电机模型选择对于稳定计算结果的影响。9、编制课程设计报告。设计要求和设计成果:1、2位同学为一组,自行分工,但任务不能重复;2、每位同学对自己的设计任务编写课程设计说明书一份;3、一组同学共同完成一份完整的设计报告;2、设计说明和报告应包含:以上设计任务每一部分的计算过程和结果分析;所编制的潮流和稳定计算源程序(主要语句应加注释);潮流计算结果(潮流图)稳定计算的功角曲线等;网络变换法求解转移阻抗的变换过程图。成绩评定成绩按五级制评分,将根据课程设计答辩情况给分。评定内容包括:设计过程中的态度(20%)、课程设计报告质量(40%)和答辩过程回答问题情况(40%)。附注:ODE函数说明Matlab提供了一阶常微分方程组求解的系列函数:ode**。包括:ode45,ode23,ode113等,还有针对刚性系统的ode15s,ode23s等。这里可采用ode45编程(大家也可选择和对比其它函数,不同编号采用的数值积分算法不同)。函数形式:[t,y]=ode45(odefun,[t1,tf],x0,options);说明: