发电机励磁系统的数学模型.docx
- 文档编号:27406214
- 上传时间:2023-06-30
- 格式:DOCX
- 页数:25
- 大小:117.04KB
发电机励磁系统的数学模型.docx
《发电机励磁系统的数学模型.docx》由会员分享,可在线阅读,更多相关《发电机励磁系统的数学模型.docx(25页珍藏版)》请在冰豆网上搜索。
发电机励磁系统的数学模型
课程设计报告
课程名称电力系统自动装置原理
设计题目发电机励磁系统数学建模
与PID控制仿真
设计时间2016-2017学年第一学期
专业年级电气133班
姓名晓
学号2012012154
提交时间2016年12月30日
成绩
指导教师帝伊谭亲跃
水利与建筑工程学院
发电机励磁系统数学建模与PID控制仿真
摘要:
本文主要进行了发电机励磁系统的数学建模和PID控制仿真。
励磁系统在电力系统的规划与控制领域都有非常重要的作用,精确的模型结构与参数是选择有效控制手段和整个电力系统仿真准确性的基础。
文过对励磁系统建模与仿真的研究,在整理系统稳定性判断理论发展的基础上,运用MATLAB软件仿真,论证了PID励磁调节可有效地改进励磁控制品质,仿真试验是调整励磁系统参数的有效措施。
关键字:
电力系统、励磁系统、根轨迹、PID、仿真
第一章绪论
1.1本课题研究意义
供给同步发电机励磁电流的电源与其附属设备统称为励磁系统。
它一般由励磁功率单元和励磁调节器两个主要部分组成。
励磁功率单元向同步发电机转子提供励磁电流;而励磁调节器则根据输入信号和给定的调节准则控制励磁功率单元的输出。
励磁系统的自动励磁调节器对提高电力系统并联机组的稳定性具有相当大的作用。
尤其是现代电力系统的发展导致机组稳定极限降低的趋势,也促使励磁技术不断发展。
在电力系统正常运行或事故运行中,同步发电机的励磁控制系统起着重要的作用。
优良的励磁控制系统不仅可以保证发电机可靠运行,提供合格的电能,而且还可有效地提高系统的技术指标。
励磁控制系统承担着如下重要任务:
(1)维持发电机端电压在给定值,当发电机负荷发生变化时,通过调节磁场的强弱来恒定机端电压。
(2)合理分配并列运行机组之间的无功分配。
(3)提高电力系统的稳定性,包括静态稳定性和暂态稳定性与动态稳定性。
(4)改善电力系统的运行条件。
(5)水轮发电机组的强行减磁[1]。
同步发电机的励磁控制系统是一个自动控制系统。
一般说来,对于自动控制系统的基本要:
首先,系统必须是稳定的;其次是系统的暂态性能应满足生产工艺所要求的暂态性能指标;其三是系统的稳态误差要满足生产的工艺要求[2]。
其中,稳定性是控制系统的首要条件,一个不稳定的系统是无法完成预期控制任务的。
因此,如何判别一个系统是否稳定以与怎样改善其稳定性乃是系统分析与设计的一个首要问题。
在经典控制理论中,对于单输入单输出线形定常系统,应用劳斯判据和胡维茨判据等代数方法间接判定系统的稳定性,而用根轨迹法与频域中的奈奎斯特判据和波德图则是更为有效的方法,它不仅用于判定系统是否稳定,还能指明改善系统稳定性的方向。
但这些方法在绘图和计算时需要花费大量的时间和精力。
MATLAB是1980年推出的用于工程计算和数值分析的交互式语言。
经过多年的完善,它已成为当前最受流行的软件,集数值分析、矩阵运算、信号处理和图形显示于一体[3]。
MATLAB有很强的绘图功能,只要写两三句代码就能得到所需要的图形。
1.2本文主要容
本课题通过对电力系统的基本知识的学习,和以往电力系统励磁控制方法的学习、总结、研究,提出了基于matlab的同步发电机励磁控制系统的仿真,文章的主要容是:
第一章,具体说明了同步发电机励磁控制系统的重要作用与其稳定性研究的意义。
第二章,通过查阅资料和之前所学过自动控制原理的基本知识,构建本文所采用的同步发电机励磁控系统数学模型,即建立了分析发电机励磁系统稳定性的传递函数。
第三章,利用控制理论中的根轨迹法研究励磁系统的稳定性。
第四章,采用文中建立的模型,进行仿真研究数值分析,是本文的重点。
第五章,对本文所做工作进行了总结并指出了本文存在的一些不足之处和下一步需要继续工作的方向。
第二章发电机励磁系统的数学模型
2.1励磁系统数学模型的发展
励磁控制对电力系统输送功率的能力和电力系统稳定性有着重要的影响。
人们对励磁控制系统的认识,是随着电力系统的不断扩大、计算机技术、控制理论等的不断进步而逐渐加深的。
在20世纪30年代,北美创建了第一个大型水电站,长距离的输电线以与慢速的动作的继电器和线路开关,使得稳定问题突出。
当时,前联、美国和加拿大等国都对此进行了深入的研究,但都没有考虑控制系统的影响,发电机是利用暂态电抗表示的,这一方法一直沿用了许多年。
从40年代到50年代,前联和北美的学者在研究励磁调节器对系统的影响时,用一阶惯性环节的比例放大器来模拟实际的励磁调节系统。
到了50年代末期,数字计算机的出现,使计算速度得到了很大的提高,虽然当时仍用暂态电抗后的恒定电势来表示发电机,但却可以计算更多的发电机和更大规模的电网。
进入60年代后,随着计算程序和计算机技术的发展,人们用更符合实际、更精确的发电机模型代替了以前的恒定电势,但同时却出现了许多不同标准的数学模型。
在此基础上,美国电气电子工程师学会(IEEE)电力生产委员会励磁系统分委会进行了大量的工作,提出了用于模拟当时存在的各种不同励磁控制系统的计算机模型和通用的专业术语,并于1968年在IEEE的学术刊物上发表。
60年代以后对励磁系统模型的研究有了更大的发展。
因为随着时代的进步出现了许多新型的励磁调节器,也采用了新的控制策略,已有的模型已不能满足要求,新的模型就不断被开发出来。
为此,IEEE于1981年又推出了新一版的励磁系统数学模型。
它比1968年版的数学模型更加详细、准确,同时也推出了新的交流励磁机模型。
1992年,IEEE的标准委员会再次对模型进行了更新,提出了附加控制特性的模型,并用标准推荐的准则将这些模型规格化。
我国在80年代前一直采用E恒定的模型,没有励磁系统模型。
N80年代初,中国电力科学研究院在电力系统分析综合程序里,开发了两种励磁模型,不但能模拟一般的直流励磁机励磁系统,也能模拟自并励和它励可控硅励磁系统。
由于IEEE的模型并不完全适合我国的情况,因此,中国电机工程学会大电机专委会励磁分委会,于1989年成立了励磁系统数学模型专家小组,对国的大型发电机励磁系统的数学模型进行了深入、广泛的研究,在1991年发表了适合于我国电力系统稳定计算的励磁系统数学模型。
又在1997年颁布了《同步电机励磁系统电力系统研究用模型》,此后,电科院又结合实际,提出了一组更为通用的新型励磁系统模型,EFM-FV共10种励磁模型。
新模型吸收了IEEE模型的精华,并溶入了新的东西,形成了自己独特的风格。
这些模型已被编进了中国版的BPA暂态稳定程序和PSASP电力系统分析综合程序里,成为电力分析计算的基础[4]。
2.2发电机励磁系统原理与分类
根据我国国家标准,同步发电机励磁系统是指“向同步发电机提供励磁的所有部件的总和”。
励磁系统分为直流励磁机励磁系统、交流励磁机励磁系统和静止励磁机励磁系统三类。
静止励磁机系统即是(晶闸管)励磁系统、与电力系统稳定计算有关的部件有励磁机、功率整流器(可控和不可控)、自动电压调节器(AVR)、电力系统稳定器(PSS)与各种限制和保护。
例如过磁通(伏赫)限制、低励磁限制和保护,过励磁限制和保护,高起始励磁系统的励磁机磁场电流瞬时过流限制等。
励磁系统数学模型由各个部件的模型组合而成。
发电机励磁系统的调节原理框图如图2-1所示。
按不同的分类标准,励磁系统有不同的种类。
按励磁系统电源供给方式的不同,励磁系统可分为三类:
(1)直流励磁机励磁系统:
包括他励和自励励磁系统。
(2)交流励磁机励磁系统:
又可分为他励静止整流器方式和他励旋转整流器方式。
这一类励磁系统采用与主机同轴的交流电机作为交流励磁电源,经可控硅整流后供给励磁绕组励磁电流,由于励磁电源来自主机以外的独立电源,故又称他励励磁系统。
(3)静止励磁系统:
又可分为交流侧自并励方式和交流侧串并联自复励方式。
2.3发电机励磁系统的数学模型
通常将励磁功率单元(励磁电源)和励磁调节器叫做励磁系统,而同步发电机和励磁系统组成同步发电机励磁控制系统。
建立同步发电机励磁控制系统的数学模型,是为了分析它本身的稳定性和动、静态性能,以与励磁系统的整定调试;特别是为了分析计算励磁控制系统对电力系统稳定性的影响,附加励磁控制(PSS)的设计与其参数的整定调试。
为简单起见,我们建立直流励磁机励磁系统的传递函数。
2.3.1励磁机的传递函数
以他励直流励磁机为例,假设其转速恒定。
励磁调节器的输出加于励磁绕组输入端、输出为励磁机电压,如图2-2所示。
励磁机绕组两端的电压方程为
(2-1)
式中—励磁机励磁绕组的磁链
—励磁机励磁绕组的电阻
—励磁机励磁绕组的电流
—励磁机励磁绕组的输入电
用磁通代换磁链,并且假定磁通与N匝键链,则可得
(2-2)
对应不同的运行点,采用饱和系数SE来表达iEE与uEE之间的非线性关系。
通常用图2-3所示的励磁机的饱和特性曲线来计与其饱和影响。
定义饱和函数为
(2-3)
SE随运行点而变,时非线性的,在整个运行围可用某一线性函数来近似的表示。
如果气隙特性的斜率是1/G,则可写出励磁机电压与励磁电流间的关系式,即
(2-4)
在恒定转速下,电压与气隙磁通成正比,即
(2-5)
又有
(2-6)
故可得
(2-7)
表示为典型的传递函数为
(2-8)
所以他励直流机的传递函数框图如图2-4所示。
2.3.2励磁调节器各单元的传递函数
励磁调节器主要由、综合放大与功率放大等单元组成。
这里以电子模拟式励磁调节器为例。
(一)电压测量比较单元的传递函数
电压测量比较单元由测量变压器、整流滤波电路与测量比较电路组成。
其中电压测量的整流滤波电路略有延时,可用一阶惯性环节来近似描述。
比较电路一般可以忽略它们的延时。
因此,测量比较电路的传递函数可表示为
(2-9)
式中—电压比例系数;
—电压测量回路的时间常数。
(二)综合放大单元的传递函数
综合放大单元在在电子型调节器中是由运算放大器组成,在电磁型调节器中则采用磁放大器。
它们的传递函数通常都可视为放大系数为的一阶惯性环节,其传递函数为
(2-10)
式中—电压放大系数;
—放大器的时间常数。
对于运算放大器,由于其响应快,可近似地认为。
此外,放大器具有一定的工作围,输出电压
综合放大单元的框图和工作特性如图2-5所示。
(三)功率放大单元的传递函数
电子型励磁调节器的功率放大单元是晶闸管整流器。
包括触发器在的晶闸管整流器的传递函数为
(2-11)
可展开为泰勒级数,略去高次项得到简化后的传递函数
(2-12)
2.3.3同步发电机的传递函数
同步发电机是电力系统中物理过程最复杂的的元件,既有机械运动过程又有电磁暂态过程,并且包含变量众多。
因此只能是根据某种目的,按照某种要求来建立相应的数学模型,这里要建立的是分析发电机励磁控制系统所用的传递函数,故发电机的近似传递函数为:
(2-13)
表示发电机的放大倍数,表示其时间常数,忽略饱和现象。
2.3.4励磁稳定器
为了提高励磁控制系统的稳定性,改善其调节品质,通常设有串、并联校正单元。
串联校正单元又叫做PID调节器。
模拟式PID传递函数为:
(2-14)
为积分环节(亦称滞后环节),它可以提高稳态增益,保证发电机的电压精度。
为微分环节(亦称超前环节),可以提高励磁电压初始上升速度,低频震荡区增益较低,可以提高励磁控制系统的稳定性。
数字式PID传递函数为:
(2-15)
并联校正单元又称为励磁系统稳定器(ESS),其模型如图2-6所示。
其输入信号可以是发电机的励磁电压(仅用于有刷励磁系统)或交流励磁机的励磁电流(有刷或无刷系统均有使用)。
输出信号的嵌入点可因调节器的不同而不同。
并联校正单元模型参数有两个和。
都应通过测量或辨识取得。
2.4励磁控制系统的传递函数
求得励磁控制系统各单元的传递函数后,可组成励磁控制系统的传递函数框图,如图2-7所示。
图2-7励磁控制系统的传递函数框图
忽略励磁机的饱和特性和放大器的饱和限制,则由图2-7可得
(2-16)
上式即为空载时同步发电机励磁控制系统的传递函数。
第三章励磁控制系统的稳定性
1948年,W.R.Evans提出了一种求特征根的简单方法,并且在控制系统的分析与设计中得到广泛的应用。
这一方法不直接求解特征方程,用作图的方法表示特征方程的根与系统某一参数的全部数值关系,当这一参数取特定值时,对应的特征根可在上述关系图中找到。
这种方法叫根轨迹法。
根轨迹法具有直观的特点,利用系统的根轨迹可以分析结构和参数已知的闭环系统的稳定性和瞬态响应特性,还可分析参数变化对系统性能的影响。
在设计线性控制系统时,可以根据对系统性能指标的要求确定可调整参数以与系统开环零极点的位置,即根轨迹法可以用于系统的分析与综合。
3.1传统方法绘制根轨迹
设某励磁控系统的参数如下:
=0s,
=8.38s,
=0.69s,
=0.04s,
=1,
=1
由图2-7得系统的开环传递函数为
其中
开环极点为:
s=-0.12,s=-1.45,s=-25
为了确定根轨迹的形状,根据根轨迹绘制原则,进行以下计算步骤:
1、根轨迹的连续性
闭环系统特征方程的系数是增益Kg的函数。
当Kg从0到无穷变化时,这些系数是连续变化的。
故特征方程的根是连续变化的,即根轨迹曲线是连续曲线。
2、根轨迹的对称性:
系统特征方程的系数是实数,其根必为实根或共轭复根。
即根轨迹位于复平面的实轴上或对称于实轴。
3、根轨迹的分支数:
3阶特征方程有3个根。
当Kg从0到无穷大变化时,3个根在复平面连续变化组成3支根轨迹。
即根轨迹的分支数等于闭环特征根的数目。
4、根轨迹的起点和终点:
⑴根轨迹起始于开环极点,终止于开环零点。
在该传递函数中开环零点数目为0小于开环极点数目3,所以有3条根轨迹终止于无限远处的零点(无限零点)。
⑵由于根轨迹是当Kg从0变到∞时闭环极点的轨迹,所以根轨迹的起点是对应于系统参数Kg=0时特征根在S平面上的位置;而根轨迹的终点则是对应于Kg=∞时特征根在S平面上的分布位置。
5、实轴上的根轨迹:
实轴上具有根轨迹的区间是:
其右侧开环实极点数和实零点数的总和为奇数。
所以该传递函数的根轨迹区间为:
(-∞,-25),(-1.45,-0.12)
6、根轨迹的会合点和分离点:
⑴若实轴上两相邻开环极点之间有根轨迹,则这两相邻极点之间必有分离点,该传递函数在(-1.45,-0.12)区间上必然存在分离点。
⑵在分离点处根轨迹的切线方向与实轴正方向的夹角,称为分离角。
设有l条根轨迹分支进入分离点又离开,则分离角为:
闭环特征方程为:
用给定值代入,得
由=0,与K>0
解得s=-0.775,这就是根轨迹在实轴上的分离点。
7.根轨迹的渐近线
系统有3个开环极点,0个开环零点,当k→∞时,有3条根轨迹分支沿着它们的渐近线趋于无限零点,而且渐近线的条数为3条,渐近线就是决定这3条根轨迹趋向无穷远处的方位。
渐近线包括两个容:
渐近线的倾角;
令k=0,±1,则=,=,=,
渐近线与实轴的交点
=
=-8.86
8、根轨迹的出射角和入射角:
根轨迹的出射角是指起始于开环极点的根轨迹在起点处的切线与水平正方向的夹角。
根轨迹的入射角是指终止于开环零点的根轨迹在终点处的切线与水平正方向的夹角。
9、根轨迹和虚轴的交点
根轨迹和虚轴相交时,系统处于临界稳定状态。
则闭环特征方程至少有一对纯虚根。
这时的增益Kc称为临界根轨迹增益。
由项的辅助多项式可计算根轨迹与虚轴的交叉点,解得
3.2用MATLAB绘制根轨迹
除用常规的方法计算出根轨迹外,还可以使用MATLAB提供的根轨迹函数,可方便、准确地绘制控制系统的根轨迹图,并可利用根轨迹图对控制系统进行分析。
源程序为:
num=1;
den=conv([1,0.12],conv([1,1.45],[1,25]));
rlocus(num,den)
axisequal
所得根轨迹图如图3-1所示
图3-1根轨迹图
第四章PID在发电机励磁系统中的应用
4.1同步发电机的励磁系统的动态指标
同步发电机的励磁系统的动态指标通常是采用机组额定转速下零起升压的参数来衡量的,同步发电机励磁控制系统的动态特性是指在外界干扰信号作用下,该系统从一个稳定工作状态变化到另一个稳定工作状态的时间响应特性。
我国同步发电机励磁系统国家标准中对同步发电机励磁自动控制系统动态特性的超调量、调节时间和摆动次数有明确规定,在我国大中型同步发电机励磁自动控制系统技术要求(GB/T7409-1997)对同步发电机动态响应的技术规定为:
1.同步发电机在空载额定电压下,当电压给定阶跃响应为10%时发电机电压超调量应不大于阶跃量的50%,摆动次数不超过3次,调节时间不超过10s。
2.当同步发电机突然零启动升压时,自动电压调节器应保证发电机端电压超调量不得超过额定值的15%,调节时间不应大于10s,电压摆动次数不应大于3次。
下面我们将通过对比加入PID控制前后发电机励磁系统的动态响应指标来观察PID调节对发电机励磁系统性能改善的作用。
4.2无PID调节的励磁系统
在simulink中建立如图4-1所示的仿真图。
图4-1无PID的励磁系统仿真原理图
4.2.1源程序
我们也可以直接用M语言编写程序来生成如图4-1所示的模型,并进行指标计算。
%没有PID的励磁系统
%0s开始零起升压
%传递函数计算
%绘制根轨迹
%指标计算
TA=0;
Td0=8.38;
TE=0.69;
TR=0.04;
KE=1;
KG=1;
KA=25;
KR=1;%KA*KR<241系统才会稳定
sysa=tf([KA],[TA,1]);
syse=tf([1],[TE,KE]);
sysg=tf([KG],[Td0,1]);
sysr=tf([KR],[TR,1]);
sysopen=series(series(sysa,series(syse,sysg)),sysr);%开环传递函数
sysclose=feedback(series(sysa,series(syse,sysg)),sysr);%闭环传递函数
figure
rlocus(sysopen);%绘制根轨迹
figure
[y,t]=step(sysclose,60);%阶跃响应,仿真时间60s
plot(t,y)
%计算0s时100%阶跃响应的指标,上升时间,超调量,调整时间,摆动次数。
tck=0;%参考时间值,即阶跃发生的时间/s
yck=1;%参考电压值
yfn=y(end);%稳态电压值
ess=(yfn-yck)/yck%稳态误差
%计算上升时间
forn=1:
length(t)
ify(n)>=yfn
break
end
end
tr=t(n)-tck
%计算超调量
overshoot=(max(y)-yfn)/yfn
%计算调整时间
forn=length(t):
-1:
1
ifabs(y(n)-yfn)>=0.05*yfn
break
end
end
ts=t(n)-tck
%计算摆动次数
bd=0;
z=y-1.05*yfn;
forn=1:
length(t)-1
ifz(n)==0&&z(n+1)==0
break
elseifz(n)==0
ifz(n-1)*z(n+1)<0
bd=bd+1;
end
elseifz(n)*z(n+1)<0
bd=bd+1;
end
end
4.2.2数值计算结果
令KA取不同的值,观察阶跃响应的各项指标
当同步发电机突然零启动升压时,自动电压调节器应保证发电机端电压超调量不得超过额定值的15%,调节时间不应大于10s,电压摆动次数不应大于3次。
表4-1指标计算汇总表
KA
稳态误差
上升时间/s
超调量
调节时间/s
摆动次数/次
5
-0.1667
3.6301
0.0276
3.0137
0
25
-0.0385
0.9836
0.3355
3.7074
2
45
-0.0217
0.6742
0.4871
4.7191
6
65
-0.0152
0.5634
0.5904
4.9765
8
ab
cd
图4-2KA取不同值的阶跃响应曲线
a.KA=5;b.KA=25;c.KA=45;d.KA=65
由以上计算结果和图4-2都可以看到,KA取值大了会导致不稳定,KA取值小了会导致误差增大,这两者矛盾。
我们采用的解决办法是加入PID控制。
4.3有PID调节的励磁系统
在目前控制系统设计中,大多采用微机控制技术,此时使用的是数字PID控制器,它是将模拟PID控制算法离散化,通过程序实现,不需要像模拟控制系统那样用硬件电路来实现,因此使系统设计更灵活、方便,由于PID控制算法具有直观的物理解释,并且能满足大多数系统的要求,因此至今PID控制仍然是常规控制系统设计应用最普遍的控制算法。
PID是以它的三种纠正算法而命名的。
这三种算法都是用加法调整被控制的数值,其输入为误差值(设定值减去测量值后的结果)或是由误差值衍生的信号。
这三种算法是:
1.比例,控制当前,误差值和一个正值的常数P(表示比例)相乘。
P只是在控制器的输出和系统的误差成比例的时候成立。
比如说,一个电热器的控制器的比例尺围是10°C,它的预定值是20°C。
那么它在10°C的时候会输出100%,在15°C的时候会输出50%,在19°C的时候输出10%,注意在误差是0的时候,控制器的输出也是0。
2.积分,控制过去,将误差值过去一段时间和(误差和)乘以一个正值的常数I。
I从过去的平均误差值来找到系统的输出结果和预定值的平均误差。
一个简单的比例系统会震荡,会在预定值的附近来回变化,因为系统无法消除多余的纠正。
通过加上一个负的平均误差比例值,平均的系统误差值就会总是减少。
所以,最终这个PID回路系统会在预定值稳定下来。
3.微分,控制将来,计算误差的一阶导,并和一个正值的常数D相乘。
这个导数的控制会对系统的改变做出反应。
导数的结果越大,那么控制系统就对输出结果做出更快速的反应。
这个D参数也是PID被称为可预测的控制器的原因。
D参数对减少控制器短期的改变很有帮助。
一些实际中的速度缓慢的系统可以不需要D参数。
在系统中加入PID模块,,在Simulink中建立如图4-3所示的模型。
图4-3有PID的励磁系统仿真原理图
4.3.1源程序
%simulink仿真之后处理ScopeData的函数
%计算0s时100%阶跃响应的指标,上升时间,超调量,调整时间,摆动次数。
tck=0;%参考时间值,即阶跃发生的时间/s
yck=1;%参考电压值
y=ScopeData(:
2);%提取y轴数据
t=ScopeData(:
1);%提取x轴数据
%根据绝对时刻寻找阶跃发生后的数据
forn=1:
length(t)
ift(n)>=tck
break
end
end
y=y(n:
end,1);%提取阶跃发生后的y轴数据
t=t(n:
end,1);%提取阶跃发生后的x轴数据
yfn=y(end);%稳态电压值
ess=(yfn-yck)/yck
plot(t,y)
%计算上升时间
forn=1:
length(t)
ify(n)>=yfn
break
e
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 发电机 系统 数学模型