[工学]电力系统稳态分析课程设计.doc
《[工学]电力系统稳态分析课程设计.doc》由会员分享,可在线阅读,更多相关《[工学]电力系统稳态分析课程设计.doc(40页珍藏版)》请在三一文库上搜索。
1、内蒙古科技大学电力系统稳态分析课程设计说明书 题 目:基于Matlab的两机五节点 网络潮流仿真计算牛拉法 姓 名: 学 号: 学 院:信息工程学院 专 业:电气工程及其自动化 班 级:10 电气 目录 目录II摘 要III内蒙古科技大学课程设计任务书IV第一章 电力系统潮流计算及MATLAB概述- 1 -第一节 潮流计算简介- 1 -第二节 MATLAB概述- 1 -第三节 矩阵的简单基本运算- 2 -第二章 潮流计算的数学模型- 3 -第一节 导纳矩阵的原理及计算方法- 3 -第二节 电力系统节点分类- 8 -第三节 潮流计算的约束条件- 9 -第三章 牛顿拉夫逊法概述- 10 -第一节
2、牛顿-拉夫逊法- 10 -第二节 用牛顿-拉夫逊法求解潮流计算的过程- 12 -第四章 牛顿拉夫逊法潮流具体计算- 17 -第一节 按照设计要求先列出手算过程- 17 -第五章 牛顿拉夫逊直角坐标潮流计算Matlab程序及运行结果- 23 -第一节 程序框图- 23 -第一节 Matlab程序- 24 -第二节 Matlab程序运行结果- 30 -第四节 本程序的符号说明- 34 -总结及感想- 35 -参考文献及资料- 36 -摘 要潮流计算的目的在于:确定电力系统的运行方式;检查系统中各元件是否过电压或过载;为电力系统继电保护的整定提供依据;为电力系统的稳定计算提供初值,为电力系统规划和经
3、济运行提供分析基础。牛顿迭代法(Newtons method)又称为牛顿-拉夫逊方法(Newton-Raphson method),牛顿-拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。MATLAB是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界,主要用于矩阵运算采用迭代法,通过建立矩阵的修正方程来依次迭代,逐步逼近真值来计算出电力网的电压,功率分布。本文采用牛顿-拉夫逊法解算电力稳态潮流,用手算和计算机算法对其进行设计。使用MATLAB软件进行编程,在很大程度
4、上节省了内存,减少了计算量。通过对本题计算我们了解了一些工程计算和解决工程问题的方法。 关键词: 电力系统潮流计算;牛顿拉夫逊法潮流计算;潮流计算的数学模型 MATLAB矩阵运算;程序; 内蒙古科技大学课程设计任务书课程名称电力系统稳态分析课程设计设计题目基于Matlab的两机五节点网络潮流仿真计算牛拉法指导教师刘景霞时间1周一、教学要求 电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点,提高学生综合能力为目标,尽可能结合实际工程进行。设计内容的安排要充分考虑学校现有的设备,设计时间及工程实际需要,并使学生初步学会运用所学知识解决工程中的实际问题。二、设计资料及参数(一)设计原始资
5、料1、待设计电气设备系统图2、电力系统网络各元件参数3、电力系统电气元件的使用规范4、电力工程电气设计手册(二)设计参考资料1、电力系统稳态分析,陈珩,中国电力出版社,2007,第三版2、电力系统分析,韩祯祥,浙江大学出版社,2005,第三版3、电力系统分析课程实际设计与综合实验,祝书萍,中国电力出版社,2007,第一版三、设计要求及成果1.根据给定的参数或工程具体要求,收集和查阅资料;学习相关软件(软件自选)。2.在给定的电力网络上画出等值电路图。3.运用计算机进行潮流计算。4.编写设计说明书。基本要求:1.编写潮流计算程序;2.在计算机上调试通过(?);3.运行程序并计算出正确结果(?);
6、4.写出课程设计报告(包括以下内容)(1份)(1)程序框图;(2)源程序;(3)符号说明表;(4)算例及计算结果5.编写计算说明书(1份)。四、进度安排根据给定的参数或工程具体要求,收集和查阅资料(半天)学习软件(MATLAB或C语言等)(一天半)编程计算复杂系统潮流计算(三天)编写计算设计书(一天)五、评分标准课程设计成绩评定依据包括以下几点:1) 工作态度(占10%);2) 基本技能的掌握程度(占20%);3) 程序编写是否合理是否有运行结果(40%);4) 课程设计说明书编写水平(占30%)。5) 分为优、良、中、合格、不合格五个等级。考核方式:设计期间教师现场检查;评阅设计报告系统接线
7、图其中节点1为平衡节点,节点2、3、4、5为PQ节点。VI 第一章 电力系统潮流计算及MATLAB概述第一节 潮流计算简介潮流计算是电力系统分析中的一种最基本的计算,它的任务是在给定的接线方式和运行条件下,确定系统的运行状态,如各母线上的电压(幅值和相角)、网络中的功率分布及功率损耗等,是电力系统的稳态计算。潮流计算是对电力系统正常运行状况的分析和计算,即电力系统中的电压、电流、功率的计算,即潮流计算;潮流计算方法很多:高斯塞德尔法、牛顿拉夫逊法、P-Q分解法、直流潮流法,以及由高斯塞德尔法、牛顿拉夫逊法演变的各种潮流计算方法。第二节 MATLAB概述目前电子计算机已广泛应用于电力系统的分析计
8、算,潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求:(1)计算速度快(2)内存需要少(3)计算结果有良好的可靠性和可信性(4)适应性好,亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强(5)简单实用。 MATLAB是一种交互式、面向对象的程序设计语言,广泛应用于工业界与学术界,主要用于矩阵运算,同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。 MATLAB程序设计语言结构完整,且具有优良的移植性,它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题,特别是关于矩阵和矢量的计算。MATLAB与C
9、语言和FORTRAN语言相比更容易被掌握。通过M语言,可以用类似数学公式的方式来编写算法,大大降低了程序所需的难度并节省了时间,从而可把主要的精力集中在算法的构思而不是编程上。 另外,MATLAB提供了一种特殊的工具:工具箱(TOOLBOXES).这些工具箱主要包括:信号处理(SIGNAL PROCESSING)、控制系统(CONTROL SYSTEMS)、神经网络(NEURAL NETWORKS)、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟(SIMULATION)等等。不同领域、不同层次的用户通过相应工具的学习和应用,可以方便地进行计算、分析及设计工作。 MATLAB
10、设计中,原始数据的填写格式是很关键的一个环节,它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计,主要应从使用的角度出发,原则是简单明了,便于修改。第三节 矩阵的简单基本运算矩阵是MATLAB数据存储的基本单元,而矩阵的运算是MATLAB语言的核心,在MATLAB语言系统中几乎一切运算均是以对矩阵的操作为基础的。矩阵的基本数学运算包括矩阵的四则运算、与常数的运算、逆运算、行列式运算、秩运算、特征值运算等基本函数运算,这里进行简单介绍。矩阵的加、减、乘运算符分别为“+,*” ,用法与数字运算几乎相同,但计算时要满足其数学要求 在MATLAB中矩阵的除法有两种形式:左除“”和右除“
11、/”。矩阵的转置,矩阵的求逆在MATLAB应用非常广泛。如:位阵的生成 eye(n) 、eye(m,n) ;零矩阵的生成 zeros(n) zeros(m,n) ;全1矩阵的生成 ones(n) ones(m,n) ;随机元素矩阵函数rand(n,m) 、rand(n) ;对角矩阵diag(V) % V=1 2 3 4; 伴随矩阵compan(P) % p=1,a1,a2,.,an ;上三角矩阵 下三角矩阵triu(B) tril(B) 矩阵的行列式det(A) ;矩阵求逆inv(A) ;矩阵的秩rank(A) ;矩阵三角分解L,U = lu(A) 矩阵奇异值分解cond(A) ;矩阵的范数N
12、 = norm(A,选项) ;矩阵的特征多项值与特征向量V,D=eig(A) ;矩阵的特征多项式、特征方程和特征根P = poly(A) 、V = roots(P) 第二章 潮流计算的数学模型第一节 导纳矩阵的原理及计算方法一 自导纳和互导纳的确定方法电力网络的节点电压方程: (2-1)为节点注入电流列向量,注入电流有正有负,注入网络的电流为正,流出网络的电流为负。根据这一规定,电源节点的注入电流为正,负荷节点为负。既无电源又无负荷的联络节点为零,带有地方负荷的电源节点为二者代数之和。为节点电压列向量,由于节点电压是对称于参考节点而言的,因而需先选定参考节点。在电力系统中一般以地为参考节点。如
13、整个网络无接地支路,则需要选定某一节点为参考。设网络中节点数为(不含参考节点),则,均为n*n列向量。为n*n阶节点导纳矩阵。节电导纳矩阵的节点电压方程: 展开为: : (2-2)是一个n*n阶节点导纳矩阵,其阶数就等于网络中除参考节点外的节点数。 节点导纳矩阵的对角元素 (i=1,2,n)成为自导纳。自导纳数值上就等于在i节点施加单位电压,其他节点全部接地时,经节点i注入网络的电流,因此,它可以定义为: (2-3)节点i的自导纳数值上就等于与节点直接连接的所有支路导纳的总和。节点导纳矩阵的非对角元素 (j=1,2,n;i=1,2,。,n;j=i)称互导纳,由此可得互导纳数值上就等于在节点i施
14、加单位电压,其他节点全部接地时,经节点j注入网络的电流,因此可定义为: (2-4) 节点j,i之间的互导纳数值上就等于连接节点j,i支路到导纳的负值。显然,恒等于。互导纳的这些性质决定了节点导纳矩阵是一个对称稀疏矩阵。而且,由于每个节点所连接的支路数总有一个限度,随着网络中节点数的增加非零元素相对愈来愈少,节点导纳矩阵的稀疏度,即零元素数与总元素的比值就愈来愈高。二非标准变比变压器等值电路变压器型等值电路更便于计算机反复计算,更适宜于复杂网络的潮流计算.双绕组变压器可用阻抗与一个理想变压器串联的电路表示.理想变压器只是一个参数,那就是变比。现在变压器阻抗按实际变比归算到低压侧为例,推导出变压器
15、型等值电路. a 双绕组变压器原理图b变压器阻抗归算到低压侧等值模型流入和流出理想变压器的功率相等 (2-6)式中, 是理想变压器的变比,和 分别为变压器高,低绕组的实际电压.从图b直接可得: (2-7) 从而可得: (2-8)式中,又因节点电流方程应具有如下形式: (2-9)将式(1-8)与(1-9)比较,得: 因此可得各支路导纳为: (2-10)由此可得用导纳表示的变压器型等值电路:图 c 潮流计算的基本方程(a) 潮流计算用的电网结构图 在潮流问题中,任何复杂的电力系统都可以归纳为以下元件(参数)组成。(1)发电机(注入电流或功率)(2)负荷(注入负的电流或功率)(3)输电线支路(电阻,
16、电抗)(4)变压器支路(电阻,电抗,变比)(5)母线上的对地支路(阻抗和导纳)(6)线路上的对地支路(一般为线路充电点容导纳)集中了以上各类型的元件的简单网络如图(a).采用导纳矩阵时,节点注入电流和节点电压构成以下线性方程组 (2-11)其中 可展开如下形式 (2-12)由于实际电网中测量的节点注入量一般不是电流而是功率,因此必须将式中的注入电流用节点注入功率来表示。节点功率与节点电流之间的关系为 (2-13)式中,因此用导纳矩阵时,PQ节点可以表示为把这个关系代入式中 ,得(2-14)式(3-4 )就是电力系统潮流计算的数学模型-潮流方程。它具有如下特点:(1)它是一组代数方程,因而表征的
17、是电力系统的稳定运行特性。(2)它是一组非线性方程,因而只能用迭代方法求其数值解。(3)由于方程中的电压和导纳既可以表为直角坐标,又可表为极坐标,因而潮流方程有多种表达形式-极坐标形式,直角坐标形式和混合坐标形式。 a。取 ,得到潮流方程的极坐标形式: (2-15)b。 取 , ,得到潮流方程的直角坐标形式: (2-16)c。取, ,得到潮流方程的混合坐标形式: (2-17)(4)它是一组n个复数方程,因而实数方程数为2n个但方程中共含4n个变量:P,Q,U和,i=1,2,n,故必须先指定2n个变量才能求解。第二节 电力系统节点分类用一般的电路理论求解网络方程,目的是给出电压源(或电流源)研究
18、网络内的电流(或电压)分布,作为基础的方程式,一般用线性代数方程式表示。然而在电力系统中,给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的,一般是给出发电机母线上发电机的有功功率(P)和母线电压的幅值(U),给出负荷母线上负荷消耗的有功功率(P)和无功功率(Q)。主要目的是由这些已知量去求电力系统内的各种电气量。所以,根据电力系统中各节点性质的不同,很自然地把节点分成三类: PQ节点对这一类点,事先给定的是节点功率(P,Q),待求的未知量是节点电压向量(U,),所以叫PQ节点。通常变电所母线都是PQ节点,当某些发电机的输出功率P。Q给定时,也作为PQ节点。PQ节点上的发电机称之为
19、PQ机(或PQ给定型发电机)。在潮流计算中,系统大部分节点属于PQ节点。 PV节点这类节点给出的参数是该节点的有功功率P及电压幅值U,待求量为该节点的无功功率Q及电压向量的相角。这类节点在运行中往往要有一定可调节的无功电源。用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线做PU节点处理。PU节点上的发电机称为PU机(或PV给定型发电机) 平衡节点在潮流计算中,这类节点一般只设一个。对该节点,给定其电压值,并在计算中取该节点电压向量的方向作为参考轴,相当于给定该点电压向量的角度为零。也就是说,对平衡节点给定的运行参数是U和,因此有城为U节点,而待求量是该
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 工学 电力系统 稳态 分析 课程设计
链接地址:https://www.31doc.com/p-1977763.html