《电力系统的MATLAB-SIMULINK仿真与应用.ppt》由会员分享，可在线阅读，更多相关《电力系统的MATLAB-SIMULINK仿真与应用.ppt（161页珍藏版）》请在三一办公上搜索。
1、第6章 电力系统稳态与暂态仿真,6.1 Powergui模块 6.2 电力系统稳态仿真 6.3 电力系统电磁暂态仿真 6.4 电力系统机电暂态仿真 习题,6.1 Powergui模块Powergui模块为电力系统稳态与暂态仿真提供了有用的图形用户分析界面。通过Powergui模块，可以对系统进行可变步长连续系统仿真、定步长离散系统仿真和相量法仿真，并实现以下功能：(1)显示测量电压、测量电流和所有状态变量的稳态值；(2)改变仿真初始状态；(3)进行潮流计算并对包含三相电机的电路进行初始化设置；(4)显示阻抗的依频特性图；,(5)显示FFT分析结果；(6)生成状态空间模型并打开“线性时不变系统”
2、(LTI)时域和频域的视窗界面；(7)生成报表，该报表中包含测量模块、电源、非线性模块和电路状态变量的稳态值，并以后缀名.rep保存；(8)设计饱和变压器模块的磁滞特性。6.1.1 主窗口功能简介MATLAB提供的Powergui模块在SimPowerSystems库中，图标如图6-1所示。,图6-1 Powergui模块图标,双击Powergui模块图标将弹出该模块的主窗口，如图6-2所示。该主窗口包含“仿真类型”(Simulation Type)和“分析工具”(Analysis Tools)两块内容，简介如下。1.仿真类型1)“相量法仿真”(Phasor simulation)单选框点击该
3、单选框后，在该单选框下方的“频率”(Frequency)文本框中输入指定的频率，进行相量法分析。若未选中该单选框，“频率”文本框显示为灰色。,图6-2 Powergui模块主窗口,2)“离散系统仿真”(Discretize electrical model)单选框点击该单选框后，在“采样时间”(Sample time)文本框中输入指定的采样时间(Ts0)，按指定的步长对离散化系统进行分析。若采样时间等于0，表示不对数据进行离散化处理，采用连续算法分析系统。若未选中该单选框，“采样时间”文本框显示为灰色。3)“连续系统仿真”(Continuous)单选框点击该单选框后，采用连续算法分析系统。4)
4、“显示分析信息”(Show message during analysis)复选框选中该复选框后，命令窗口中将显示系统仿真过程中的相关信息。,2.分析工具1)“稳态电压电流分析”(Steady-State Voltages and Currents)按键打开稳态电压电流分析窗口，显示模型文件的稳态电压和电流。2)“初始状态设置”(Initial States Setting)按键打开初始状态设置窗口，显示初始状态，并允许对模型的初始电压和电流进行更改。3)“潮流计算和电机初始化”(Load Flow and Machine Initialization)按键打开潮流计算和电机初始化窗口。,4)
5、“LTI视窗”(Use LTI Viewer)按键打开窗口，使用“控制系统工具箱”(Control System Toolbox)的LTI视窗。5)“阻抗依频特性测量”(Impedance vs.Frequency Measurement)按键打开窗口，如果模型文件中含阻抗测量模块，该窗口中将显示阻抗依频特性图。6)“FFT分析”(FFT Analysis)按键打开FFT分析窗口。7)“报表生成”(Generate Report)按键打开窗口，产生稳态计算的报表。,8)“磁滞特性设计工具”(Hysteresis Design Tool)按键打开窗口，对饱和变压器模块和三相变压器模块的铁芯进行磁
6、滞特性设计。9)“计算RLC线路参数”(Compute RLC Line Parameters)按键打开窗口，通过导线型号和杆塔结构计算架空输电线的RLC参数。6.1.2 稳态电压电流分析窗口打开“稳态电压电流分析”窗口如图6-3所示。该窗口中含有以下内容：,图6-3“稳态电压电流分析”窗口,(1)“稳态值”(Steady state value)列表框：显示模型文件中指定的电压、电流稳态值。(2)“单位”(Units)下拉框：选择将显示的电压、电流值是“峰值”(Peak)还是“有效值”(RMS)。(3)“频率”(Frequency)下拉框：选择将显示的电压、电流相量的频率。该下拉框中列出模型
7、文件中电源的所有频率。(4)“状态”(States)复选框：显示稳态下电容电压和电感电流的相量值。默认状态为不选。(5)“测量”(Measurements)复选框：显示稳态下测量模块测量到的电压、电流相量值。默认状态为选中。,(6)“电源”(Sources)复选框：显示稳态下电源的电压、电流相量值。默认状态为不选。(7)“非线性元件”(Nonlinear elements)复选框：显示稳态下非线性元件的电压、电流相量值。默认状态为不选。(8)“格式”(Format)下拉框：在下拉列表框中选择要观测的电压和电流的格式。“浮点格式”(floating point)以科学计数法显示5位有效数字；“最
8、优格式”(best of)显示4位有效数字并且在数值大于9999时以科学计数法表示；最后一个格式直接显示数值大小，小数点后保留2位数字。默认格式为“浮点格式”。(9)“更新稳态值”(Update Steady State Values)按键：重新计算并显示稳态电压、电流值。,6.1.3 初始状态设置窗口仿真时，常常希望仿真开始时系统处于稳态，或者仿真开始时系统处于某种初始状态，这时，就可以使用“初始状态设置”按键。打开“初始状态设置”窗口如图6-4所示。该窗口中含有以下内容：(1)“初始状态”(Initial state values for simulation)列表框：显示模型文件中状态变
9、量的名称和初始值。(2)“设置到指定状态”(Set selected state)文本框：对“初始状态”列表框中选中的状态变量进行初始值设置。(3)“设置所有状态量”(Reset all States)：选择从“稳态”(To Steady State)或者“零初始状态”(To Zero)开始仿真。,(4)“加载状态”(Reload States)：选择从“指定的文件”(From File)中加载初始状态或直接以“当前值”(From Diagram)作为初始状态开始仿真。,图6-4“初始状态设置”窗口,(5)“应用”(Apply)按键：用设置好的参数进行仿真。(6)“返回”(Revert)按键：
10、返回到“初始状态设置”窗口打开时的原始状态。(7)“保存初始状态”(Save Initial States)按键：将初始状态保存到指定的文件中。(8)“格式”(Format)下拉框：选择观测的电压和电流的格式。格式类型见节。默认格式为“浮点格式”。(9)“分类”(Sort values by)下拉框：选择初始状态值的显示顺序。“默认顺序”(Default order)是按模块在电路中的顺序显示初始值；“状态序号”(State number)是按状态空间模型中状态变量的序号来显示初始值；“类型”(Type)是按电容和电感来分类显示初始值。默认格式为“默认顺序”。,6.1.4 潮流计算和电机初始化
11、窗口打开“潮流计算和电机初始化”窗口如图6-5所示。该窗口中含有以下内容：(1)“电机潮流分布”(Machines load flow)列表框：显示“电机”(Machines)列表框中选中电机的潮流分布。(2)“电机”(Machines)列表框：显示简化同步电机、同步电机、非同步电机和三相动态负荷模块的名称。选中该列表框中的电机或负荷后，才能进行参数设置。,图6-5 潮流计算和电机初始化窗口,(3)“节点类型”(Bus type)下拉框：选择节点类型。对于“PV节点”(P&V Generator)，可以设置电机的端口电压和有功功率；对于“PQ节点”(P&Q Generator)，可以设置电机的
12、有功和无功功率；对于“平衡节点”(Swing Bus)，可以设置终端电压UAN的有效值和相角，同时需要对有功功率进行预估。如果选择了非同步电机模块，则仅需要输入电机的机械功率；如果选择了三相动态负荷模块，则需要设置该负荷消耗的有功和无功功率。(4)“终端电压UAB”(Terminal voltage UAB)文本框：对选中电机的输出线电压进行设置(单位：V)。,(5)“有功功率”(Active power)文本框：设置选中的电机或负荷的有功功率(单位：W)。(6)“预估的有功功率”(Active power guess)文本框：如果电机的节点类型为平衡节点的话，设置迭代起始时刻电机的有功功率。
13、(7)“无功功率”(Reactive power)文本框：设置选中的电机或负荷的无功功率(单位：var)。(8)“电压UAN的相角”(Phase of UAN voltage)文本框：当电机的节点类型为平衡节点时，该文本框被激活。指定选中电机a相相电压的相角。,(9)“负荷频率”(Load flow frequency)下拉框：对潮流计算的频率进行设置，通常为60 Hz或者50 Hz。(10)“负荷潮流初始状态”(Load flow initial condition)下拉框：常常选择默认设置“自动”(Auto)，使得迭代前系统自动调节负荷潮流初始状态。如果选择“从前一个结果开始”(Start
14、 from previous solution)，则负荷潮流的初始值为上次仿真结果。如果改变电路参数、电机的功率分布和电压后负荷潮流不收敛，就可以选择这个选项。,(11)“更新电路和测量结果”(Update Circuit&Measurements)按键：更新电机列表，更新电压相量和电流相量，更新“电机潮流分布”列表框中的功率分布。其中的电机电流是最近一次潮流计算的结果。该电流值储存在电机模块的“初始状态参数”(Initial conditions)文本框中。(12)“更新潮流分布”(Update Load Flow)按键：根据给定的参数进行潮流计算。,6.1.5 LTI视窗打开“LTI视窗”
15、窗口如图6-6所示。该窗口中含有以下内容：(1)“系统输入”(System inputs)列表框：列出电路状态空间模型中的输入变量，选择需要用到LTI视窗的输入变量。(2)“系统输出”(System outputs)列表框：列出电路状态空间模型中的输出变量，选择需要用到LTI视窗的输出变量。(3)“打开新的LTI视窗”(Open New LTI Viewer)按键：产生状态空间模型并打开选中的输入和输出变量的LTI视窗。(4)“打开当前LTI视窗”(Open in current LTI Viewer)按键：产生状态空间模型并将选中的输入和输出变量叠加到当前LTI视窗。,图6-6“LTI视窗”
16、窗口,图6-7“阻抗依频特性测量”窗口,6.1.6 阻抗依频特性测量窗口打开“阻抗依频特性测量”窗口如图6-7所示。,该窗口中含有以下内容：(1)图表：窗口左上侧的坐标系表示阻抗频率特性，左下侧的坐标系表示相角频率特性。(2)“测量模块”(Measurement)列表框：列出模型文件中的阻抗测量模块，选择需要显示依频特性的阻抗测量模块。使用“Ctrl”键可选择多个阻抗显示在同一个坐标中。(3)“范围”(Range)文本框：指定频率范围(单位：Hz)。该文本框中可以输入任意有效的MATLAB表达式。(4)“对数阻抗”(Logarithmic Impedance)单选框：坐标系纵坐标的阻抗以对数值
17、形式表示。,(5)“线性阻抗”(Linear Impedance)单选框：坐标系纵坐标的阻抗以线性形式表示。(6)“对数频率”(Logarithmic Frequency)单选框：坐标系横坐标的频率以对数值形式表示。(7)“线性频率”(Linear Frequency)单选框：坐标系横坐标的频率以线性形式表示。(8)“网格”(Grid)复选框：选中该复选框，阻抗频率特性图和相角频率特性图上将出现网格。默认设置为无网格。,(9)“更新后保存数据”(Save data when updated)复选框：选中该复选框后，该复选框下面的“工作间变量名”(Workspace variable name)
18、文本框被激活，数据以该文本框中显示的变量名被保存在工作间中。复数阻抗和对应的频率保存在一起。其中频率保存在第1列，阻抗保存在第2列。默认设置为不保存。(10)“显示/保存”(Display/Save)按键：开始阻抗依频特性测量并显示结果，如果选择了“更新后保存数据”(Save data when updated)复选框，数据将保存到指定位置。6.1.7 FFT分析窗口打开“FFT分析”窗口如图6-8所示。,图6-8“FFT分析”窗口,该窗口中含有以下内容：(1)图表：窗口左上侧的图形表示被分析信号的波形，窗口左下侧的图形表示该信号的FFT分析结果。(2)“结构”(Structure)下拉框：列
19、出工作间中带时间的结构变量的名称。使用下拉菜单选择要分析的结构变量。这些结构变量名可以由“示波器”(Scope)模块产生。打开示波器模块参数对话框，选中“数据历史”(Data history)标签页，如图6-9所示，在“变量名”(Variable name)文本框中输入该结构变量的名称，在“存储格式”(Format)下拉框中选择“带时间的结构变量”(Structure with time)。,图6-9 由“示波器”模块产生结构变量,这些结构变量名也可以由“到工作间”(To Workspace)模块产生，如图6-10所示。打开该模块的对话窗口，并在“变量名”(Variable name)文本框中
20、输入该结构变量的名称，“存储格式”(Save format)为“结构变量”(Structure)。,图6-10 由“到工作间”模块产生结构变量(a)“到工作间”模块图标；(b)“到工作间”模块参数对话框,(3)“输入变量”(Input)下拉框：列出被选中的结构变量中包含的输入变量名称，选择需要分析的输入变量。(4)“信号路数”(Signal number)下拉框：列出被选中的输入变量中包含的各路信号的名称。例如，若要把a、b、c三相电压绘制在同一个坐标中，可以通过把这三个电压信号同时送入示波器的一个通道来实现，这个通道就对应一个输入变量，该变量含有3路信号，分别为a相、b相和c相电压。(5)“
21、开始时间”(Start time)文本框：指定FFT分析的起始时间。(6)“周期个数”(Number of cycles)文本框：指定需要进行FFT分析的波形的周期数。,(7)“显示FFT窗/显示完整信号”(Display FFT window/Display entire signal)下拉框：选择“显示完整信号”(Display entire signal)，将在左上侧插图中显示完整的波形；选择“显示FFT窗”(Display FFT window)将在左上侧插图中显示指定时间段内的波形。(8)“基频”(Fundamental frequency)文本框：指定FFT分析的基频(单位：Hz)
22、。(9)“最大频率”(Max Frequency)文本框：指定FFT分析的最大频率(单位：Hz)。(10)“频率轴”(Frequency axis)下拉框：在下拉框中选择“赫兹”(Hertz)使频谱的频率轴单位为Hz，选择“谐波次数”(Harmonic order)使频谱的频率轴单位为基频的整数次倍数。,(11)“显示类型”(Display style)下拉框：频谱的显示类型可以是“以基频或直流分量为基准的柱状图”(Bar(relative to Fund.or DC)、“以基频或直流分量为基准的列表”(list(relative to Fund.or DC)、“指定基准值下的柱状图”(Bar
23、(relative to specified base)、“指定基准值下的列表”(List(relative to specified base)四种类型。(12)“基准值”(Base value)文本框：当“显示类型”下拉框中选择“指定基准值下的柱状图”或“指定基准值下的列表”时，该文本框被激活，输入谐波分析的基准值。(13)“显示”(Display)按键：显示FFT分析结果。,6.1.8 报表生成窗口打开“报表生成”窗口如图6-11所示。该窗口中含有以下内容：(1)“报表中包含的内容”(Items to include in the report)：包括“稳态”(Steady state)
24、复选框、“初始状态”(Initial states)复选框和“电机负荷潮流”(Machine load flow)复选框，这三个复选框可以任意组合。(2)“报表中的频率”(Frequency to include in the report)下拉框：选择报表中包含的频率。可以是60 Hz或者全部，默认为60 Hz。,(3)“单位”(Units)下拉框：选择以“峰值”(Peak)还是“有效值”(Units)显示数据。(4)“格式”(Format)下拉框：与节相关内容相同。(5)“报表生成”(Create Report)按键：生成报表并保存。,图6-11“报表生成”窗口,6.1.9 磁滞特性设计工
25、具窗口打开“磁滞特性设计工具”窗口如图6-12所示。该窗口中含有以下内容：(1)“磁滞曲线”(Hysteresis curve for file)图表：显示设计的磁滞曲线。(2)“分段”(Segments)下拉框：将磁滞曲线做分段线性化处理，并设置磁滞回路第1象限和第4象限内曲线的分段数目。左侧曲线和右侧曲线关于原点对称。(3)“剩余磁通”(Remanent flux Fr)文本框：设置零电流对应的剩磁。(4)“饱和磁通”(Saturation flux Fs)文本框：设置饱和磁通。,图6-12“磁滞特性设计工具”窗口,(5)“饱和电流”(Saturation current Is)文本框：设
26、置饱和磁通对应的电流。(6)“矫顽电流”(Coercive current Ic)文本框：设置零磁通对应的电流。(7)“矫顽电流处的斜率”(dF/dI at coercive current)文本框：指定矫顽电流点的斜率。(8)“饱和区域电流”(Saturation region currents)文本框：设置磁饱和后磁化曲线上各点所对应的电流值，仅需设置第1象限值。注意该电流向量的长度必须和“饱和区域磁通”的向量长度相同。,(9)“饱和区域磁通”(Saturation region fluxes)文本框：设置磁饱和后磁化曲线上各点所对应的磁通值，仅需要设置第1象限值。注意该向量的长度必须和“
27、饱和区域电流”的向量长度相同。(10)“变压器额定参数”(Transfo Nominal Parameters)文本框：指定额定功率(单位：VA)、一次绕组的额定电压值(单位：V)和额定频率(单位：Hz)。(11)“参数单位”(Parameter units)下拉框：将磁滞特性曲线中电流和磁通的单位由国际单位制(SI)转换到标幺制(p.u.)或者由标幺制转换到国际单位制。,(12)“放大磁滞区域”(Zoom around hysteresis)复选框：选中该复选框，可以对磁滞曲线进行放大显示。默认设置为“可放大显示”。6.1.10 计算RLC线路参数窗口打开“计算RLC线路参数”窗口如图6-1
28、3所示。该窗口可分为三个子窗口，左上窗口输入常用参数(单位、频率、大地电阻和文件注释)，右上窗口输入线路的几何结构，下方窗口输入导线的特性。,图6-13“计算RLC线路参数”窗口,1.常用参数子窗口(1)“单位”(Units)下拉框：在下拉菜单中，选择以“米制”(metric)为单位时，以厘米作为导线直径、几何平均半径GMR和分裂导线直径的单位，以米作为导线间距离的单位；选择以“英制”(english)为单位时，以英寸作为导线直径、几何平均半径GMR和分裂导线直径的单位，以英尺作为导线间距离的单位。(2)“频率”(Frequency)文本框：指定RLC参数所用的频率(单位：Hz)。(3)“大地
29、电阻”(Ground resistivity)文本框：指定大地电阻(单位：m)。输入0表示大地为理想导体。,(4)“注释”(Comments)多行文本框：输入关于电压等级、导线类型和特性等的注释。该注释将与线路参数一同被保存。2.线路几何结构子窗口(1)“导线相数”(Number of phase conductors(bundle)文本框：设置线路的相数。(2)“地线数目”(Number of ground wires(bundle)文本框：设置大地导线的数目。(3)导线结构参数表：输入导线的“相序”(Phase number)、“水平挡距”(X)、“垂直挡距”(Y tower)、“挡距中央
30、的高度”(Y min)、“导线的类型”(Conductor(bundle)type)共五个参数。,3.导线特性子窗口(1)“导线类型的个数”(Number of conductor types or bundle types)文本框：设置需要用到导线类型(单导线或分裂导线)的数量。假如需要用到架空导线和接地导线，该文本框中就要填“2”。(2)“导线内电感计算方法”(Conductor internal inductance evaluated from)下拉框：选择用“直径/厚度”(T/D ratio)、“几何平均半径”(Geometric Mean Radius(GMR)或者“1英尺(米)间
31、距的电抗”(Reactance Xa at 1-foot spacing)进行内电感计算。(3)“考虑导线集肤效应”(Include conductor skin effect)复选框：选中该复选框后，在计算导线交流电阻和电感时将考虑集肤效应的影响。若未选中，电阻和电感均为常数。,(4)导线特性参数表：输入导线“外径”(Conductor Outside diameter)、“T/D”(Conductor T/D ratio)、“GMR”(Conductor GMR)、“直流电阻”(Conductor DC resistance)、“相对磁导率”(Conductor relative perm
32、eability)、“分裂导线中的子导线数目”(Number of conductors per bundle)、“分裂导线的直径”(Bundle diameter)、“分裂导线中1号子导线与水平面的夹角”(Angle of conductor 1)共八个参数。(5)“计算RLC参数”(compute RLC parameters)按键：点击该按键后，将弹出RLC参数的计算结果窗口。,(6)“保存”(Save)按键：点击该按键后，线路参数以及相关的GUI信息将以后缀名.mat被保存。(7)“加载”(Load)按键：点击该按键后，将弹出窗口，选择“典型线路参数”(Typical line dat
33、a)或“用户定义的线路参数”(User defined line data)将线路参数信息加载到当前窗口。,6.2 电力系统稳态仿真6.2.1 连续系统仿真以例4.4为例，说明Powergui模块在电力系统稳态仿真分析中的应用。【例6.1】计算例4.4的潮流分布，并利用Powergui模块实现连续系统的稳态分析。解：例4.4看似一个极简单的电路，可是它的潮流分析却不简单。已知电源侧的电压幅值、相角和额定电压下的负荷大小，这意味着在求出负荷侧电压前，负荷未知，因此按潮流计算的理论进行分析，这种电路必须经过一些假设，并反复迭代计算，直到逼近真实解。但是，利用Powergui模块却很简单。,图6-1
34、4 例6.1的系统仿真图,(1)重新布置系统仿真图，如图6-14所示。相对例4.4，本例添加了新的模块，新增模块的名称及提取路径如表6-2所示。,表6-2 例6.1新增模块的名称及提取路径,(2)Powergui仿真。打开Powergui模块窗口，选中“连续系统仿真”(Continuous)单选框后，点击“稳态电压电流分析”按键，出现稳态电压电流分析窗口，如图6-15所示。,图6-15 例6.1的稳态电压电流分析窗口,图中，状态变量用电流和电压的符号加上电感或电容的模块名表示，例如“I1_110Mvar”表示110Mvar负荷上的电流大小，“Uc_input：Pi Line”表示PI形线路左侧
35、并联电容器上的电压大小，“Uc_output：Pi Line”表示PI形线路右侧并联电容器上的电压大小，“I1_section _1：Pi Line”表示第一段PI形线路串联电感上的电流大小。电压源上电压的名称与系统电压源名称一致，如“Vs”表示电压源Vs上的电压大小。测量模块测得的电压值用测量模块的名称表示，如“V1”表示电压表模块V1测得的电压大小(电源侧电压)，“V2”表示电压表模块V2测得的电压大小，“I2”表示电流表模块I2测得的电流大小。,由图6-15可见，PI形电路左侧的电压相量为244.88 0.19 kV，PI形电路右侧的电压相量为166.413.66kV，PI形电路上的电流
36、为529.21.0886.12 A，负荷侧电流为610.1786.15 A。因此，负荷大小为,(6-1),图6-16所示为直接通过测量模块得到的PI形电路两侧电压和实际负荷大小。图中波形从上到下依次为PI形电路左侧电压、PI形电路右侧电压、负荷侧有功功率、负荷侧无功功率。该结果与Powergui所得结论一致。利用测量模块还可以得到电源侧电流的大小，这样就可以很容易求到线路各处的功率分布、功率损耗和电压损耗了。读者可以动手试一试，并和电力系统分析课程中潮流计算的相关内容进行比较。仿真时，特别是在对含电机的电路或者节点较多的电力系统，往往希望仿真开始时系统处于某种指定的初始状态。下例将说明对系统初
37、始状态的设置方法。,图6-16 例6.1的仿真波形图,【例6.2】改变例6.1的初始状态，使仿真开始时，电路处于零初始状态。解：打开Powergui模块窗口，点击“初始状态设置”按键，出现初始状态设置窗口，如图6-17所示。点击按键“to zero”将全部状态变量的初始值设为0，回到模型文件主窗口，重新开始仿真，得到仿真波形如图6-18所示。,图6-17 例6.2的初始状态设置窗口,图6-18 例6.2的仿真波形图,由仿真波形图可见，当仿真从零初始状态开始时，线路两侧的电压在仿真初期不是理想正弦波。简单地说，就是合闸引起了一个暂态过程。从“波”的角度来讲，就是合闸后在线路上引起了“波”过程，“
38、波”在阻抗不连续点(线路两端)不断进行折射和反射，形成了图6-18所示的波形。如果系统中存在电机模块，电机模块的初始状态需要通过“潮流计算和电机初始化”按键进行设置。,【例6.3】利用Powergui模块设置例4.2同步发电机的初始状态。解：例4.2中直接在同步发电机参数对话框中输入了初始状态，这些初始状态实际上是由Powergui模块自动产生的。首先在同步发电机参数对话窗口(见图6-19)的初始状态文本框中输入初始状态0 0 0 0 0 0 0 0 1，确认后关闭该对话框。,图6-19 例6.3的同步发电机参数对话窗口,点击“潮流计算和电机初始化”按键，进入潮流计算和电机初始化设置窗口，如图
39、6-20所示。该窗口中只显示了一个电机模块，名为SM_p.u.，选中该模块，并点击“更新潮流分布”按键，更新负荷潮流，窗口左侧将显示更新后的负荷潮流分布情况。现在重新打开同步发电机参数对话窗口(见图6-21)，可以看到，初始状态已经自动更新了。同步发电机的初始状态设置好后，剩下的仿真工作就可交给读者了。,图6-20 例6.3的潮流计算和电机初始化设置窗口,图6-21 更新后的同步发电机参数对话窗口,6.2.2 离散系统仿真连续系统仿真通常采用变步长积分算法。对小系统而言，变步长算法通常比定步长算法快，但是对含大量状态变量或非线性模块(如电力电子开关)的系统而言，采用定步长离散算法的优越性更为明
40、显。对系统进行离散化时，仿真的步长决定了仿真的精确度。步长太大可能导致仿真精度不足，步长太小又可能大大增加仿真运行时间。判断步长是否合适的唯一方法就是用不同的步长试探并找到最大时间步长。对于50 Hz或60 Hz的系统，或者带有整流电力电子设备的系统，通常2050 s的时间步长都能得到较好的仿真结果。,对于含强迫换流电力电子开关器件的系统，由于这些器件通常都运行在高频下，因此需要适当地减小时间步长。例如，对运行在8 kHz左右的脉宽调制(PWM)逆变器的仿真，需要的时间步长为1 s。【例6.4】将例4.4中的PI形电路的段数改为10，对系统进行离散化仿真并比较离散系统和连续系统的仿真结果。解：
41、(1)重新布置系统仿真图，如图6-22所示。,图6-22 例6-4的系统仿真图,(2)参数设置。双击例4-4模型文件中PI形电路模块，打开参数对话框，将分段数改为10，如图6-23所示。打开Powergui模块，选择“离散系统仿真”单选框，设置采样时间为25e-6 s，如图6-24所示。仿真时该系统将以25 s的采样率进行离散化。,图6-23 例6-4的PI形电路参数对话窗,图6-24 例6-4的Powergui模块参数对话窗,图6-25 例6-4的仿真参数设置对话窗,由于系统离散化了，因此在该系统中无连续的状态变量，所以不需要采用变步长的积分算法进行仿真。打开菜单SimulationConf
42、iguration parameters对话框，按图6-25设置仿真参数，选择“定步长”(Fixed-step)和“离散”(discrete(no continuous states)选项并设置步长为25 s。,图6-26 例6-4的仿真运行时间,(3)仿真运行时间比较。为了得到仿真运行时间，在MATLAB命令窗口输入如下命令：tic;sim(gcs);toc仿真结束后，仿真所用的时间将以秒为单位显示在MATLAB命令窗口中，如图6-26所示。,可见，离散化系统后，仿真运行时间为0.188 s。将离散系统的采样时间设为0并回到连续系统的仿真状态，仿真算法改为连续积分算法ode23tb，可以得到
43、连续系统仿真需要的运行时间为0.219 s。因此，离散积分算法比连续积分算法更快。(4)仿真精度比较。为了比较两种方法的精确度，执行以下三种仿真：连续系统仿真，Ts=0 s；离散系统仿真，Ts=25 s；离散系统仿真，Ts=50 s。,如图6-27所示，双击并打开V2示波器模块，选择“参数”(Parameters)项，在打开的窗口中选择“数据历史”(Data history)，去掉“仅保留最新的数据点”(Limit data points to last)复选框，这样可以观察到整个仿真过程中的波形变化。选中“将数据保存到工作空间”(Save data points to workspace)复
44、选框，将变量名指定为V2，格式为“列”(Array)。,图6-27 例6-4示波器V2的参数设置(a)波形；(b)参数标签页,开始连续系统仿真，仿真结束时间选为0.02 s。仿真结束后，在MATLAB命令窗口中输入命令：V2C=V2;这样，电压V2被保存在变量V2C中。重新开始仿真，将系统离散化，设置仿真步长Ts=25 s，注意仿真参数中的步长设置也要改为25 s，仿真结束时间为0.02 s。仿真结束后，将电压V2保存在变量V2d25中。再次仿真，设置仿真步长为Ts=50 s。仿真结束后，将电压V2保存在变量V2d50中。,在MATLAB命令窗口中输入如下语句，可画出三种情况下的电压波形，如图
45、6-28所示。plot(V2C(:,1),V2C(:,2),V2d25(:,1),V2d25(:,2),V2d50(:,1),V2d50(:,2),图6-28 三种仿真方法波形比较,使用图形窗口中的放大功能，将目标集中到0.0045 s附近观察三种仿真的差别。如图6-29所示，25 s下的仿真结果与50 s的仿真结果一致，连续系统的仿真结果除了步长不同，结果也相同。可见，本例中，选择50 s的步长不但可以提高计算速度而且不影响仿真的精确度。,图6-29 放大后三种仿真波形比较,6.2.3 相量法仿真相量是代表特定频率下的正弦电压和电流的复数，可以用直角坐标或者极坐标表示。相量法是电力系统正弦稳
46、态分析的主要手段。它只关心系统中电压电流的相角和幅值，不需要求解电力系统状态方程，不需要特殊的算法，因此计算速度快得多。必须清楚的是，相量法给出的解是在特定频率下的解。【例6.5】用相量法分析例6.4。解：(1)参数设置。打开Powergui模块，选择“相量法分析”单选框，并在“频率”对话框中将频率改为50 Hz。关闭Powergui模块，模型文件主窗口中的Powergui模块图标显示为“相量法”(Phasors)分析，如图6-30所示。,图6-30 例6.5的Powergui模块相量法分析图标,图6-31 例6.5的电压测量模块V1,打开电压测量模块V1，选择“幅值相角”(Magnitude
47、-Angle)模式，如图6-31所示。电压测量模块V2也选择幅值相角模式。,(2)仿真。开始仿真，得到输电线路送端V1和受端V2的电压幅值和相角，如图6-32所示。可见，V1侧电压幅值为1 p.u.，相角为0.19；V2侧电压幅值为0.67 p.u.，相角为3.66。这和图6-16稳态分析的结论一致。,图6-32 例6.5的仿真结果(a)V1；(b)V2,6.3 电力系统电磁暂态仿真SIMULINK的电力系统暂态仿真过程通过机械开关设备，如“断路器”(circuit breakers)模块或者电力电子设备的开断实现。6.3.1 断路器模块SimPowerSystems库提供的断路器模块可以对开
48、关的投切进行仿真。断路器合闸后等效于电阻值为Ron的电阻元件。Ron是很小的值，相对外电路可以忽略。断路器断开时等效于无穷大电阻，熄弧过程通过电流过零时断开断路器完成。开关的投切操作可以受外部或内部信号的控制。,外部控制方式时，断路器模块上出现一个输入端口，输入的控制信号必须为0或者1，其中0表示切断，1表示投合；内部控制方式时，切断时间由模块对话框中的参数指定。如果断路器初始设置为1(投合)，SimPowerSystems库自动将线性电路中的所有状态变量和断路器模块的电流进行初始化设置，这样仿真开始时电路处于稳定状态。断路器模块包含Rs-Cs缓冲电路。如果断路器模块和纯电感电路、电流源和空载
49、电路串联，则必须使用缓冲电路。,图6-33 单相断路器模块图标,带有断路器模块的系统进行仿真时需要采用刚性积分算法，如ode23tb、odel5s，这样可以加快仿真速度。1.单相断路器模块外部控制方式、带缓冲电路和不带缓冲电路的单相断路器模块图标如图6-33所示。,双击断路器模块，弹出该模块的参数对话框如图6-34。该对话框中含有如下参数：(1)“断路器电阻”(Breaker resistance Ron)文本框：断路器投合时的内部电阻(单位：)。断路器电阻不能为0。(2)“初始状态”(Initial state)文本框：断路器初始状态。断路器为合闸状态，输入1，对应的图标显示投合状态；输入0
50、，表示断路器为断开状态。(3)“缓冲电阻”(Snubber resistance Rs)文本框：并联缓冲电路中的电阻值(单位：)。缓冲电阻值设为inf时，将取消缓冲电阻。,图6-34 单相断路器模块参数对话框,(4)“缓冲电容”(Snubber capacitance Cs)文本框：并联缓冲电路中的电容值(单位：F)。缓冲电容值设为0时，将取消缓冲电容；缓冲电容值设为inf时，缓冲电路为纯电阻性电路。(5)“开关动作时间”(Switching times)文本框：采用内部控制方式时，输入一个时间向量以控制开关动作时间。从开关初始状态开始，断路器在每个时间点动作一次。例如，初始状态为0，在时间向
