电力系统建模与分析软件：Mathematica二次开发_（3）.电力系统建模方法.docx

PAGE1

PAGE1

电力系统建模方法

1.电力系统的基本模型

在电力系统建模中，基本模型是理解系统运行和分析的基础。这些模型包括发电机模型、负荷模型、线路模型、变压器模型等。本节将详细介绍这些基本模型的原理和在Mathematica中的实现方法。

1.1发电机模型

发电机是电力系统的核心组件之一，其模型包括同步发电机和异步发电机。同步发电机模型通常使用Park变换和Swing方程来描述，而异步发电机模型则涉及转子和定子的动态方程。

1.1.1同步发电机模型

同步发电机模型通常使用Park变换来将三相交流系统转换为两相旋转坐标系下的直流系统。Park变换的原理如下：

ψ

其中，ψd、ψq和ψ0分别是d轴、q轴和0轴的磁链，ψa、ψb和

在Mathematica中，可以使用以下代码实现Park变换：

(*定义Park变换矩阵*)

parkTransformMatrix[θ_]:={

{Cos[θ],Cos[θ-120Degree],Cos[θ+120Degree]},

{-Sin[θ],-Sin[θ-120Degree],-Sin[θ+120Degree]},

{1/2,1/2,1/2}

}

(*定义三相磁链*)

ψa=1;ψb=0.5;ψc=-0.5;

(*计算d轴、q轴和0轴的磁链*)

θ=30Degree;

{ψd,ψq,ψ0}=parkTransformMatrix[θ].{ψa,ψb,ψc};

(*输出结果*)

{ψd,ψq,ψ0}

1.1.2异步发电机模型

异步发电机模型涉及转子和定子的动态方程。定子方程和转子方程分别为：

$$

$$

$$

$$

其中，ψsα、ψsβ是定子磁链，ψrα、ψrβ是转子磁链，Rs、Rr是定子和转子的电阻，

在Mathematica中，可以使用以下代码实现异步发电机模型的动态方程：

(*定义参数*)

Rs=0.1;Rr=0.2;

Lsα=0.05;Lsβ=0.05;

Lrα=0.03;Lrβ=0.03;

ωr=100;

(*定义定子和转子的磁链和电流*)

ψsα=ψsα[t];ψsβ=ψsβ[t];

ψrα=ψrα[t];ψrβ=ψrβ[t];

isα=isα[t];isβ=isβ[t];

irα=irα[t];irβ=irβ[t];

(*定义定子和转子的动态方程*)

eq1=D[ψsα,t]==-Rs*ψsα+Lsα*isα+ωr*ψsβ;

eq2=D[ψsβ,t]==-Rs*ψsβ+Lsβ*isβ-ωr*ψsα;

eq3=D[ψrα,t]==-Rr*ψrα+Lrα*irα+ωr*ψrβ;

eq4=D[ψrβ,t]==-Rr*ψrβ+Lrβ*irβ-ωr*ψrα;

(*初始条件*)

initialConditions={ψsα[0]==1,ψsβ[0]==0,ψrα[0]==0,ψrβ[0]==1};

(*求解动态方程*)

solution=NDSolve[{eq1,eq2,eq3,eq4,initialConditions},{ψsα,ψsβ,ψrα,ψrβ},{t,0,1}];

(*绘制结果*)

Plot[Evaluate[{ψsα[t],ψsβ[t],ψrα[t],ψrβ[t]}/.solution],{t,0,1},

PlotLegends-{ψsα,ψsβ,ψrα,ψrβ},

PlotRange-All]

2.负荷模型

负荷模型用于描述电力系统中负荷的动态行为，常见的负荷模型包括恒定功率负荷、恒定电流负荷和恒定阻抗负荷。

2.1恒定功率负荷

恒定功率负荷模型假设负荷的有功功率和无功功率保持不变，其方程为：

P

Q

其中，P和Q是负荷的有功功率和无功功率，V是负荷点的电压，I是负荷点的电流，?是功率因数角。

在Mathematica中，可以使用以下代码实现恒定功率负荷模型：

(*定义参数*)

V=120;(*电压，单位：V*)

P=1000;(*有功功率，单位：W*)

Q=500;(*无功功率，单位：Var*)

(*计算功率因数角*)

φ=ArcTan[Q/P];

(*计算负荷