CN104079228B - 双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 - Google Patents
双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 Download PDFInfo
- Publication number
- CN104079228B CN104079228B CN201410323683.4A CN201410323683A CN104079228B CN 104079228 B CN104079228 B CN 104079228B CN 201410323683 A CN201410323683 A CN 201410323683A CN 104079228 B CN104079228 B CN 104079228B
- Authority
- CN
- China
- Prior art keywords
- axle
- delta
- alpha
- rotor
- omega
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000013016 damping Methods 0.000 title claims abstract description 40
- 230000001052 transient effect Effects 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 238000002347 injection Methods 0.000 claims abstract description 27
- 239000007924 injection Substances 0.000 claims abstract description 27
- 230000009466 transformation Effects 0.000 claims abstract description 25
- 230000004907 flux Effects 0.000 claims abstract description 9
- 238000004804 winding Methods 0.000 claims description 94
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000014509 gene expression Effects 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 10
- 230000010354 integration Effects 0.000 description 6
- 230000005284 excitation Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000010355 oscillation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000000889 atomisation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Landscapes
- Control Of Eletrric Generators (AREA)
Abstract
本发明提供了一种双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法。(1).从三相电路方程和磁链方程出发,经派克变换得到电路方程,使用带阻尼的隐式梯形积分方法对电路方程进行差分化,得到双馈异步电机定子侧和转子侧的戴维南等值电路;(2).提出双馈异步电机与外部电磁暂态网络的接口方法,以定子侧和转子侧三相注入电流源和伴随导纳矩阵的形式与外部网络相联;(3).使用带阻尼的隐式梯形积分方法对双馈异步电机的转子运动方程进行差分化,得到多质量块的双馈异步电机转动方程,通过迭代计算实现了双馈异步电机的带阻尼的隐式梯形积分建模和计算。本发明方法提高了双馈异步电机模型的电磁暂态仿真精度,增强了模型的稳定性。
Description
技术领域
本发明属于电力领域,具体涉及一种双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法。
背景技术
双馈异步电机也称交流励磁电机,定子接50Hz工频电网,转子接交流励磁电源,随着交流励磁自动控制系统对转子励磁电流的频率、幅值大小和相位的调节,双馈异步电机能够调节转子转速,在不同的工况下运行,而定子输出电压和频率可以维持不变,运行方式灵活。双馈即是指有两个能量流动通道,双馈异步电机的定子、转子都可以与电网交换能量,实现能量从定子和转子到电网的两个通道流动。
随着风能在全球范围内的大规模开发,双馈型风力发电机组由于具有能够变速运行、风能转换效率高、具有最大风能捕获能力并减小了风电机组机械部件所受应力、实现有功、无功功率解耦控制、改善风电场功率因数及电压稳定等优点,已经成为世界风电市场上的主流机型之一。作为双馈型风力发电机组关键设备之一的双馈异步电机,其定子侧直接通过变压器与电网相连,转子侧则通过背靠背换流器与电网相连接,接入所需低频励磁电流。
由于双馈型风力发电机组含有高频开断的电力电子换流器,在通过数字仿真研究双馈型风力发电机组的控制策略、动态特性、故障响应及低电压穿越能力等方面时通常需要使用微秒级的仿真步长,在电磁暂态环境下进行建模和仿真计算。在现有的电力系统电磁暂态仿真研究中,通常使用隐式梯形积分法或欧拉法对双馈异步电机进行建模。这种建模方法在仿真应用中对一些现象会导致电机的数值振荡,存在着精度和稳定性较差的问题,仿真结果不理想。
发明内容
本发明的目的在于提供一种双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,该方法通过在积分方法中引入阻尼系数,对电机的电压方程、磁链方程和转子运动方程进行差分化处理后联入外部网络,综合了隐式梯形积分法和后退欧拉法的优点,能够有效抑制开关动作引发的双馈异步电机及系统的数值振荡,增强了仿真方法的稳定性。
为实现上述目的,本发明的技术方案是:一种双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,包括如下步骤,
步骤S1:从三相电路方程和磁链方程出发,经派克变换得到dq0坐标轴下的电路方程,使用带阻尼的隐式梯形积分方法对标幺化后的dq0轴电路方程进行差分化,得到双馈异步电机定子侧和转子侧的dq0轴戴维南等值电路;
步骤S2:提出双馈异步电机与外部电磁暂态网络的接口方法,以定子侧和转子侧三相注入电流源和伴随导纳矩阵的形式与外部网络相联;
步骤S3:使用带阻尼的隐式梯形积分方法对双馈异步电机的转子运动方程进行差分化,得到多质量块的双馈异步电机转动方程,计算转子角和转速,通过迭代计算实现了双馈异步电机的带阻尼的隐式梯形积分建模和计算。
在本发明实施例中,在所述步骤S1中,将双馈异步电机的磁链方程代入三相电路方程,使用派克变换和标幺化处理得到d轴、q轴和0轴的标幺化向量方程;以d轴为例,其方程如下式所示:
式中,电压,电流,旋转电势,电阻,
电抗,ω0=2πf0,f0为基波频率;、分别是定子绕组d轴、转子绕组d轴的电压;、分别是定子绕组d轴、转子绕组d轴的电流;、分别是定子绕组d轴、转子绕组d轴的旋转电势;、分别是定子绕组和转子绕组的电阻,、和分别为定子绕组d轴电抗、转子绕组d轴电抗和d轴互感抗;
记为带阻尼的隐式梯形积分方法的阻尼系数,用带阻尼隐式梯形积分方法对上式进行差分可得:
式中,为当前计算时刻,为计算步长,为上一计算时刻,,为历史项,该项中的电流、电压和旋转电势各量均为上一时刻的值,对当前计算时刻而言是已知量。
记,其中,,;将差分方程展开整理后可得定子d轴和转子d轴的戴维南等值电路:
式中,、分别为定子绕组的d轴等值电阻、d轴等值电压源,、分别为转子绕组的d轴等值电阻、d轴等值电压源;记定子绕组和转子绕组的d轴等值历史项的表达式分别为和,戴维南等值电路中各量的详细表达式如下:
同理可得q轴和0轴的戴维南等值电路;将dq0轴的等值电路方程按照定子侧等值电路和转子侧等值电路分别排列如下,
定子侧:
转子侧:
。
式中,、分别是定子绕组q轴、转子绕组q轴的电压,、分别是定子绕组0轴、转子绕组0轴的电压;、分别是定子绕组q轴、转子绕组q轴的电流,、分别是定子绕组0轴、转子绕组0轴的电流;、分别是定子绕组q轴、转子绕组q轴的等值电压源,、分别是定子绕组0轴、转子绕组0轴的等值电压源;、分别是定子绕组和转子绕组0轴的等值电阻。
在本发明实施例中,在所述步骤S2中,首先对定子侧的戴维南等值电路进行诺顿等值,得到dq0轴注入电流源,如下所示:
通过反派克变换将定子侧dq0量转换成abc三相量,产生三相注入电流源和伴随导纳矩阵的并联模型,其电流源为:
式中为反派克变换中使用的同步轴和定子轴之间的夹角,、和为定子abc三相的注入电流,
dq0坐标轴下的定子绕组等值电阻矩阵为对角矩阵,三个对角元依次为、和;通过派克变换得到abc坐标轴下的三相等值电阻矩阵,记其逆为三相并联导纳矩阵,有:
三相等值电阻矩阵为对称矩阵,对角线元素记为,非对角元素记为,,;
然后对转子侧的戴维南等值电路进行诺顿等值,得到dq0轴注入电流源,如下所示:
通过反派克变换将转子侧dq0量转换成abc三相量,产生三相注入电流源和伴随导纳矩阵的并联模型,其电流源为:
式中为反派克变换中使用的同步轴和转子轴之间的夹角,、和为转子abc三相的注入电流;
dq0坐标轴下的转子绕组等值电阻矩阵为对角矩阵,三个对角元依次为、和;通过派克变换得到abc坐标轴下的三相等值电阻矩阵,记其逆为三相并联导纳矩阵,有:
上式中的对角线元素记为,非对角元素记为,,。
在本发明实施例中,在所述步骤S3中,双馈异步电机的转子运动方程使用标幺化的多质量块系统描述,以满足风力发电场合的应用需求,如下式所示:
式中,为双馈异步电机各质量块的转子角,为各质量块的转速,为惯性时间常数的对角线矩阵,为阻尼系数的三对角线矩阵,为弹性系数的三对角线矩阵,和为各质量块上的电磁转矩和机械转矩,ω0=2πf0,f0为基波频率;为带阻尼的隐式梯形积分方法的阻尼系数,为当前计算时刻,为计算步长,为上一计算时刻;
用带阻尼的隐式梯形积分法对转子运动方程进行差分,可得:
式中:
。
相较于现有技术,本发明具有以下有益效果:本发明方法克服双馈异步电机的电磁暂态建模方法在一些现象中存在电机仿真计算数值振荡的问题,在对双馈异步电机标幺化的电路方程和转子运行方程差分化过程中在积分环节引入了阻尼系数,积分方法能够随变化而变化:当时,积分方法为纯隐式梯形积分方法;当时,积分方法为后退欧拉法;当在0和1之间时,积分方法介于梯形积分法和后退欧拉法之间;通过对阻尼系数数值的灵活选择,能够综合隐式梯形积分法和后退欧拉法的优点,在保证仿真精度的同时有效抑制双馈异步电机的数值振荡;得到计及阻尼效果的定子侧和转子侧戴维南等值电路方程、转动方程后,通过诺顿等值和反派克变换,最终得到定子侧和转子侧与系统的接口方法,即以三相注入电流源和伴随导纳矩阵并联的形式并入系统,计算转动方程求解转子角和转速,并通过迭代计算保证仿真结果的精确性。
附图说明
图1(a)、1(b)、1(c)分别为双馈异步电机的定子侧的dq0轴的戴维南等值电路图。
图1(d)、1(e)、1(f)分别为双馈异步电机的转子侧的dq0轴的戴维南等值电路图。
图2(a)为双馈异步电机的定子侧与系统的接口方法示意图。
图2(b)为双馈异步电机的转子侧与系统的接口方法示意图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
本发明一种双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,包括如下步骤,
步骤S1:从三相电路方程和磁链方程出发,经派克变换得到dq0坐标轴下的电路方程,使用带阻尼的隐式梯形积分方法对标幺化后的dq0轴电路方程进行差分化,得到双馈异步电机定子侧和转子侧的dq0轴戴维南等值电路;
步骤S2:提出双馈异步电机与外部电磁暂态网络的接口方法,以定子侧和转子侧三相注入电流源和伴随导纳矩阵的形式与外部网络相联;
步骤S3:使用带阻尼的隐式梯形积分方法对双馈异步电机的转子运动方程进行差分化,得到多质量块的双馈异步电机转动方程,计算转子角和转速,通过迭代计算实现了双馈异步电机的带阻尼的隐式梯形积分建模和计算。
以下为本发明的具体实施例。
本发明的双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,包括如下步骤:
(1).从三相电路方程和磁链方程出发,经派克变换得到dq0坐标轴下的电路方程,使用带阻尼的隐式梯形积分方法对标幺化后的dq0轴电路方程进行差分化,得到双馈异步电机定子侧和转子侧的dq0轴戴维南等值电路;
(2).提出双馈异步电机与外部电磁暂态网络的接口方法,以定子侧和转子侧三相注入电流源和伴随导纳矩阵的形式与外部网络相联;
(3).使用带阻尼的隐式梯形积分方法对双馈异步电机的转子运动方程进行差分化,得到多质量块的双馈异步电机转动方程,计算转子角和转速,通过迭代计算实现了双馈异步电机的带阻尼的隐式梯形积分建模和计算。
在所述步骤(1)中,将双馈异步电机的磁链方程代入三相电路方程,使用派克变换和标幺化处理得到d轴、q轴和0轴的标幺化向量方程。以d轴为例,其方程如下式所示:
式中,电压,电流,旋转电势,电阻,
电抗,ω0=2πf0,f0为基波频率。、分别是定子绕组d轴、转子绕组d轴的电压;、分别是定子绕组d轴、转子绕组d轴的电流;、分别是定子绕组d轴、转子绕组d轴的旋转电势;、分别是定子绕组和转子绕组的电阻,、和分别为定子绕组d轴电抗、转子绕组d轴电抗和d轴互感抗。
记为带阻尼的隐式梯形积分方法的阻尼系数,用带阻尼隐式梯形积分方法对上式进行差分可得:
式中,为当前计算时刻,为计算步长,为上一计算时刻,,为历史项,该项中的电流、电压和旋转电势各量均为上一时刻的值,对当前计算时刻而言是已知量。
记,其中,,;将差分方程展开整理后可得定子d轴和转子d轴的戴维南等值电路:
式中,、分别为定子绕组的d轴等值电阻、d轴等值电压源,、分别为转子绕组的d轴等值电阻、d轴等值电压源。记定子绕组和转子绕组的d轴等值历史项的表达式分别为和,戴维南等值电路中各量的详细表达式如下:
同理可以得到q轴和0轴的戴维南等值电路。
将dq0轴的等值电路方程按照定子侧等值电路和转子侧等值电路分别排列后,如图1(a)至图1(f)所示:
定子侧:
转子侧:
。
式中,、分别是定子绕组q轴、转子绕组q轴的电压,、分别是定子绕组0轴、转子绕组0轴的电压;、分别是定子绕组q轴、转子绕组q轴的电流,、分别是定子绕组0轴、转子绕组0轴的电流;、分别是定子绕组q轴、转子绕组q轴的等值电压源,、分别是定子绕组0轴、转子绕组0轴的等值电压源;、分别是定子绕组和转子绕组0轴的等值电阻。
在所述步骤(2)中,首先对定子侧的戴维南等值电路进行诺顿等值,得到dq0轴注入电流源,如下所示:
使用反派克变换将定子侧dq0量转换成abc三相量,产生三相注入电流源和伴随导纳的并联模型,如图2(a)所示。其电流源为:
式中为反派克变换中使用的同步轴和定子轴之间的夹角,、和为定子abc三相的注入电流;
dq0坐标轴下的定子绕组等值电阻矩阵为对角矩阵,三个对角元依次为、和;通过派克变换得到abc坐标轴下的三相等值电阻矩阵,记其逆为三相并联导纳矩阵,有:
三相等值电阻矩阵为对称矩阵,对角线元素记为,非对角元素记为,,;
然后对转子侧的戴维南等值电路进行诺顿等值,得到dq0轴注入电流源,如下所示:
使用反派克变换将转子侧dq0量转换成abc三相量,产生三相注入电流源和伴随导纳矩阵的并联模型,如图2(b)所示。其电流源为:
式中为反派克变换中使用的同步轴和转子轴之间的夹角,、和为转子abc三相的注入电流;
dq0坐标轴下的转子绕组等值电阻矩阵为对角矩阵,三个对角元依次为、和;通过派克变换得到abc坐标轴下的三相等值电阻矩阵,记其逆为三相并联导纳矩阵,有:
上式中的对角线元素记为,非对角元素记为,,。
在所述步骤(3)中,双馈异步电机的转子运动方程使用标幺化的多质量块系统描述,以满足风力发电等场合的应用需求,如下式所示:
式中,为双馈异步电机各质量块的转子角度(弧度),为各质量块的转速,为惯性时间常数的对角线矩阵,为阻尼系数的三对角线矩阵,为弹性系数的三对角线矩阵,和为各质量块上的电磁转矩和机械转矩,ω0=2πf0,f0为基波频率。
用带阻尼的隐式梯形积分法对转子运动方程进行差分,可得:
式中:
。
通过以下实施例对双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法做进一步描述。
如图2(a)至图2(b)所示,前述步骤3将双馈异步电机处理为三相注入电流源和伴随导纳的并联模型,分别接入定子侧网络和转子侧网络,参与全网计算。三相注入电流源为时变量,需要每步进行计算;相空间的伴随导纳矩阵和为非时变量,仅与电机参数和阻尼系数有关,因此可以在计算的预处理阶段求出双馈异步电机的伴随导纳矩阵和,将其分别并入网络方程导纳矩阵的对应位置,实现网络方程的联立,以节省计算时间。
在电磁暂态仿真阶段,双馈异步电机的求解方法为:
1)预测旋转电势、和、,预测定子侧和转子侧电压、和、,预测转子角和转速,根据前述步骤(1)的方法,计算双馈异步电机定子绕组和转子绕组的dq0轴戴维南等值电路的等值电压源,如图1(a)至图1(f)所示;
2)根据前述步骤(2)的方法,计算双馈异步电机定子绕组和转子绕组的abc三相等值注入电流源,如图2(a)至图2(b)所示。
3)将等值电流源并入网络方程右端项的对应位置,求解网络方程,得到双馈异步电机定子侧和转子侧的三相电压;
4)将双馈异步电机的定子侧和转子侧三相电压转换为dq0坐标系统下的变量:、、和、、;根据图1(a)至图1(f)所示的dq0轴戴维南等值电路,计算定子电流、、和转子电流、、;
5)根据双馈异步电机的dq0轴电流,计算电磁功率;
6)根据前述步骤(3)的方法,计算转动方程的转子角和转速;
7)设定收敛判据为,根据转速的预测值和计算值之间的误差是否小于阀值进行判断。如果满足则本时步计算完成;如果不满足,根据计算结果更新预测量,返回4),继续进行求解直至收敛。
需要声明的是,本发明内容及具体实施方式意在证明本发明所提供技术方案的实际应用,不应解释为对本发明保护范围的限定。本领域技术人员在本发明的精神和原理启发下,可作各种修改、等同替换、或改进。但这些变更或修改均在申请待批的保护范围内。
Claims (4)
1.一种双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,其特征在于:包括如下步骤,
步骤S1:从三相电路方程和磁链方程出发,经派克变换得到dq0坐标轴下的电路方程,使用带阻尼的隐式梯形积分方法对标幺化后的dq0轴电路方程进行差分化,得到双馈异步电机定子侧和转子侧的dq0轴戴维南等值电路;
步骤S2:提出双馈异步电机与外部电磁暂态网络的接口方法,以定子侧和转子侧三相注入电流源和伴随导纳矩阵的形式与外部网络相联;
步骤S3:使用带阻尼的隐式梯形积分方法对双馈异步电机的转子运动方程进行差分化,得到多质量块的双馈异步电机转动方程,计算转子角和转速,通过迭代计算实现了双馈异步电机的带阻尼的隐式梯形积分建模和计算。
2.根据权利要求1所述的双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,其特征在于:在所述步骤S1中,将双馈异步电机的磁链方程代入三相电路方程,使用派克变换和标幺化处理得到d轴、q轴和0轴的标幺化向量方程;以d轴为例,其方程如下式所示:
式中,电压v=[vdvD]T,电流i=[idiD]T,旋转电势u=[uduD]T,电阻 电抗 ω0=2πf0,f0为基波频率;vd、vD分别是定子绕组d轴、转子绕组d轴的电压;id、iD分别是定子绕组d轴、转子绕组d轴的电流;ud、uD分别是定子绕组d轴、转子绕组d轴的旋转电势;Rs、Rr分别是定子绕组和转子绕组的电阻,Xd、XD和Xmd分别为定子绕组d轴电抗、转子绕组d轴电抗和d轴互感抗;
记α为带阻尼的隐式梯形积分方法的阻尼系数,用带阻尼隐式梯形积分方法对上式进行差分可得:
式中,t为当前计算时刻,△t为计算步长,t-△t为上一计算时刻,h(t-△t)=[hd(t-△t)hD(t-△t)]T,为历史项,该项中的电流、电压和旋转电势各量均为上一时刻的值,对当前计算时刻而言是已知量;
记 其中 将差分方程展开整理后可得定子d轴和转子d轴的戴维南等值电路:
vd(t)=ed-Rdid(t)
vD(t)=eD-RDiD(t)
式中,Rd、ed分别为定子绕组的d轴等值电阻、d轴等值电压源,RD、eD分别为转子绕组的d轴等值电阻、d轴等值电压源;记定子绕组和转子绕组的d轴等值历史项的表达式分别为Histd(t-△t)和HistD(t-△t),戴维南等值电路中各量的详细表达式如下:
Rd=Rdd-RdD*RdD/RDD
ed=ud(t)-RdD/RDD*uD(t)+Histd(t-△t)+RdD/RDD*vD(t)
Histd(t-△t)=hd(t-△t)-RdD/RDD*hD(t-△t)
RD=RDD-RdD*RdD/Rdd
eD=uD(t)-RdD/Rdd*ud(t)+HistD(t-△t)+RdD/Rdd*vd(t)
HistD(t-△t)=hD(t-△t)-RdD/Rdd*hd(t-△t)
同理可得q轴和0轴的戴维南等值电路;将dq0轴的等值电路方程按照定子侧等值电路和转子侧等值电路分别排列如下,
定子侧:
vd(t)=ed-Rdid(t)
vq(t)=eq-Rdiq(t)
v0s(t)=e0s-R0si0s(t)
转子侧:
vD(t)=eD-RDiD(t)
vQ(t)=eQ-RDiQ(t)
v0r(t)=e0r-R0ri0r(t)
式中,vq、vQ分别是定子绕组q轴、转子绕组q轴的电压,v0s、v0r分别是定子绕组0轴、转子绕组0轴的电压;iq、iQ分别是定子绕组q轴、转子绕组q轴的电流,i0s、i0r分别是定子绕组0轴、转子绕组0轴的电流;eq、eQ分别是定子绕组q轴、转子绕组q轴的等值电压源,e0s、e0r分别是定子绕组0轴、转子绕组0轴的等值电压源;R0s、R0r分别是定子绕组和转子绕组0轴的等值电阻。
3.根据权利要求2所述的双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,其特征在于:在所述步骤S2中,首先对定子侧的戴维南等值电路进行诺顿等值,得到dq0轴注入电流源,如下所示:
通过反派克变换将定子侧dq0量转换成abc三相量,产生三相注入电流源和伴随导纳矩阵的并联模型,其电流源为:
式中θ1为反派克变换中使用的同步轴和定子轴之间的夹角,isa,source、isb,source和isc,source为定子abc三相的注入电流,
dq0坐标轴下的定子绕组等值电阻矩阵为对角矩阵,三个对角元依次为Rd、Rd和R0s;通过派克变换得到abc坐标轴下的三相等值电阻矩阵Requivs,记其逆为三相并联导纳矩阵Gequivs,有:
三相等值电阻矩阵为对称矩阵,对角线元素记为R1,非对角元素记为R2,
然后对转子侧的戴维南等值电路进行诺顿等值,得到dq0轴注入电流源,如下所示:
通过反派克变换将转子侧dq0量转换成abc三相量,产生三相注入电流源和伴随导纳矩阵的并联模型,其电流源为:
式中θ2为反派克变换中使用的同步轴和转子轴之间的夹角,ira,source、irb,source和irc,source为转子abc三相的注入电流;
dq0坐标轴下的转子绕组等值电阻矩阵为对角矩阵,三个对角元依次为RD、RD和R0r;通过派克变换得到abc坐标轴下的三相等值电阻矩阵Requivr,记其逆为三相并联导纳矩阵Gequivr,有:
上式中的对角线元素记为R3,非对角元素记为R4,
4.根据权利要求1所述的双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法,其特征在于:在所述步骤S3中,双馈异步电机的转子运动方程使用标幺化的多质量块系统描述,以满足风力发电场合的应用需求,如下式所示:
式中,δ为双馈异步电机各质量块的转子角,ω为各质量块的转速,TJ为惯性时间常数的对角线矩阵,D为阻尼系数的三对角线矩阵,K为弹性系数的三对角线矩阵,Te和Tm为各质量块上的电磁转矩和机械转矩,ω0=2πf0,f0为基波频率,α为带阻尼的隐式梯形积分方法的阻尼系数,t为当前计算时刻,△t为计算步长,t-△t为上一计算时刻;
用带阻尼的隐式梯形积分法对转子运动方程进行差分,可得:
δ(t)=MM·δ(t-△t)+NN·ω(t-△t)+PP·△T(t-△t)+QQ·△T(t)+RR
式中:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410323683.4A CN104079228B (zh) | 2014-07-09 | 2014-07-09 | 双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410323683.4A CN104079228B (zh) | 2014-07-09 | 2014-07-09 | 双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104079228A CN104079228A (zh) | 2014-10-01 |
CN104079228B true CN104079228B (zh) | 2016-06-29 |
Family
ID=51600309
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410323683.4A Active CN104079228B (zh) | 2014-07-09 | 2014-07-09 | 双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104079228B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104361159A (zh) * | 2014-10-31 | 2015-02-18 | 国家电网公司 | 一种大规模电力系统暂态稳定时间空间并行仿真方法 |
CN105470950B (zh) * | 2015-12-03 | 2020-06-30 | 中国电力科学研究院 | 故障分析用永磁直驱风电场电磁暂态等值模型建立方法 |
CN106354973A (zh) * | 2016-09-14 | 2017-01-25 | 南方电网科学研究院有限责任公司 | 异步感应电机实时仿真系统 |
CN108573094B (zh) * | 2018-03-21 | 2019-03-15 | 华北电力大学 | 同步发电机的vbr电磁暂态仿真模型的建立方法及系统 |
CN110807244A (zh) * | 2019-09-25 | 2020-02-18 | 中国电力科学研究院有限公司 | 一种电力电子开关的电磁暂态仿真方法及装置 |
CN110912186B (zh) * | 2019-11-27 | 2021-04-06 | 国网冀北电力有限公司电力科学研究院 | 双馈型风力发电虚拟同步发电机建模方法 |
CN111697889B (zh) * | 2020-05-06 | 2021-11-05 | 南方电网科学研究院有限责任公司 | 一种基于时域变换的异步电动机仿真建模方法及装置 |
CN112214899A (zh) * | 2020-10-16 | 2021-01-12 | 哈尔滨理工大学 | 双轴励磁同步发电机的2s-dirk电磁暂态建模方法 |
CN112383252B (zh) * | 2020-10-30 | 2022-05-06 | 华北电力科学研究院有限责任公司 | 双馈发电机组励磁控制系统标幺方法及装置 |
CN115021638A (zh) * | 2022-06-15 | 2022-09-06 | 南方电网科学研究院有限责任公司 | 一种高效同步机的电磁暂态建模方法、系统及设备 |
CN114977165A (zh) * | 2022-06-15 | 2022-08-30 | 南方电网科学研究院有限责任公司 | 基于饱和效应的同步机电磁暂态建模方法、系统及设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609575A (zh) * | 2012-01-19 | 2012-07-25 | 浙江大学 | 一种基于隐式数值积分的电力系统暂态稳定仿真方法 |
CN102436534B (zh) * | 2012-01-13 | 2013-07-24 | 浙江省电力公司电力科学研究院 | 一种电力系统仿真方法 |
EP2553802B1 (de) * | 2010-03-30 | 2014-01-22 | Sew-Eurodrive GmbH & Co. KG | Verfahren zum ansteuern eines optimalen betriebspunktes bei einer synchronmaschine und eine umrichtergespeiste synchronmaschine |
-
2014
- 2014-07-09 CN CN201410323683.4A patent/CN104079228B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2553802B1 (de) * | 2010-03-30 | 2014-01-22 | Sew-Eurodrive GmbH & Co. KG | Verfahren zum ansteuern eines optimalen betriebspunktes bei einer synchronmaschine und eine umrichtergespeiste synchronmaschine |
CN102436534B (zh) * | 2012-01-13 | 2013-07-24 | 浙江省电力公司电力科学研究院 | 一种电力系统仿真方法 |
CN102609575A (zh) * | 2012-01-19 | 2012-07-25 | 浙江大学 | 一种基于隐式数值积分的电力系统暂态稳定仿真方法 |
Non-Patent Citations (2)
Title |
---|
基于机电暂态-电磁暂态混合仿真的电网合环分析计算系统;朱雨晨等;《电力系统保护与控制》;20121201;第40卷(第23期);第73-79页 * |
直流电网模型和仿真的发展与挑战;李亚楼等;《电力系统自动化》;20140225;第38卷(第4期);第127-135页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104079228A (zh) | 2014-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104079228B (zh) | 双馈异步电机的带阻尼的隐式梯形积分电磁暂态建模方法 | |
Chang et al. | Fault current analysis of type-3 WTs considering sequential switching of internal control and protection circuits in multi time scales during LVRT | |
Wang et al. | Dynamic stability enhancement and power flow control of a hybrid wind and marine-current farm using SMES | |
Fan et al. | Modal analysis of a DFIG-based wind farm interfaced with a series compensated network | |
Bongiorno et al. | On control of static synchronous series compensator for SSR mitigation | |
Mohammadpour et al. | SSR analysis of a DFIG-based wind farm interfaced with a gate-controlled series capacitor | |
Shao et al. | Equivalent modeling and comprehensive evaluation of inertia emulation control strategy for DFIG wind turbine generator | |
Zeng et al. | Co-ordinated control strategy for hybrid wind farms with PMSG and FSIG under unbalanced grid voltage condition | |
CN106356884A (zh) | 一种基于虚拟同步机的光伏并网控制方法、装置及系统 | |
CN109495031A (zh) | 一种基于ESO-CPC的激光追踪控制系统电机Simulink仿真方法 | |
Wu et al. | Small signal stability of synchronous motor-generator pair for power system with high penetration of renewable energy | |
Tan et al. | Multi-time scale model reduction strategy of variable-speed pumped storage unit grid-connected system for small-signal oscillation stability analysis | |
CN106294959A (zh) | 模型参考自适应控制与双馈风机降阶模型相结合的建模仿真方法 | |
Li et al. | DFIG-based wind farm electromagnetic dynamic model and impact on protection relay of transmission network | |
CN106055817B (zh) | 用于电磁暂态分析的双馈异步发电机双工况模型建模方法 | |
Song et al. | High voltage ride-through control method for DFIG-based wind turbines based on Resonant Controller | |
CN114792055A (zh) | 一种基于暂态电抗后电势的异步电动机等效惯量评估方法 | |
Yan et al. | Impact of DFIG-based WTs on subsynchronous oscillation damping of SG based on the net electrical damping method | |
Huang et al. | Improving rotational inertia of power system with variable speed synchronous condenser | |
CN110518629A (zh) | 一种风电机组的虚拟同步控制方法、装置及控制器 | |
Wang et al. | Available inertia estimation of induction machine and its effect on frequency response of power systems with wind penetration | |
Moussavi et al. | Small signal stability assessment of power systems with large-scale wind farms | |
Sai Sindhura et al. | Control And Modeling Of Doubly Fed Induction Machine For Wind Turbines | |
Arvani et al. | Power oscillation damping controller for the power system with high wind power penetration level | |
Jing et al. | Control strategy of current balance based on VSG of DFIG under unbalanced grid voltage conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |