CN105117543B - 一种考虑多种闭锁模态的基于全桥子模块mmc的等效仿真方法 - Google Patents

一种考虑多种闭锁模态的基于全桥子模块mmc的等效仿真方法 Download PDF

Info

Publication number
CN105117543B
CN105117543B CN201510516122.0A CN201510516122A CN105117543B CN 105117543 B CN105117543 B CN 105117543B CN 201510516122 A CN201510516122 A CN 201510516122A CN 105117543 B CN105117543 B CN 105117543B
Authority
CN
China
Prior art keywords
mrow
msub
equivalent
mfrac
submodule
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.)
Expired - Fee Related
Application number
CN201510516122.0A
Other languages
English (en)
Other versions
CN105117543A (zh
Inventor
徐政
刘高任
李海英
肖晃庆
唐庚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
NR Electric Co Ltd
Original Assignee
Zhejiang University ZJU
State Grid Corp of China SGCC
NR Electric Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, State Grid Corp of China SGCC, NR Electric Co Ltd filed Critical Zhejiang University ZJU
Priority to CN201510516122.0A priority Critical patent/CN105117543B/zh
Publication of CN105117543A publication Critical patent/CN105117543A/zh
Application granted granted Critical
Publication of CN105117543B publication Critical patent/CN105117543B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Inverter Devices (AREA)

Abstract

本发明公开了一种考虑多种闭锁模态的基于全桥子模块MMC的等效仿真方法,包括:获取对应桥臂的桥臂电流、桥臂各全桥子模块IGBT的开关信号;根据运行参数构建全桥子模块的等效电路;根据全桥子模块等效电路建立相应的仿真系统,并进行仿真。本发明MMC等效仿真方法利用数值计算以及电磁仿真的内部机理,将子模块等效为戴维南等效支路,避免了大量电力电子器件的频繁开断,大幅减少可仿真计算时间;同时引入真实二极管器件,考虑了全桥子模块的各种闭锁方式,能够在保证仿真速度的前提下对全桥子模块各闭锁状态进行精确模拟。

Description

一种考虑多种闭锁模态的基于全桥子模块MMC的等效仿真 方法
技术领域
本发明属于电力系统仿真技术领域,具体涉及一种考虑多种闭锁模态的基于全桥子模块MMC的等效仿真方法。
背景技术
2001年,德国学者R.Marquardt首先提出了模块化多电平换流器(modularmultilevel converter,MMC)拓扑,该型拓扑的桥臂采用基本运行单元级联的形式,避免大量开关器件直接串联,不存在一致触发等问题,因此在近年来得到了学界与工业界的广泛关注,并被迅速应用到工程实际当中。
如图1所示,MMC拓扑的桥臂由多个子模块串联而成,子模块的结构根据使用场合的需要可分为不同的类型,目前比较常用的有半桥子模块、全桥子模块以及箝位双子模块三类。其中全桥子模块因为可靠性强、处理直流故障能力优秀而更具应用前景。
系统仿真建模是对MMC拓扑研究的基础。MMC拓扑包含大量的电力电子器件,在正常运行时,这些电力电子器件将频繁开断,这会对系统的仿真计算产生很大负担。特别是在实际工程中,随着MMC电压等级以及容量的增大,单个桥臂所需要串联的子模块数量将随之增大。如果每个子模块都采用传统方法用真实模型搭建,一方面搭建难度大,另一方面后续的仿真计算也将花费很长时间。因此,有学者提出将子模块中的电容器用时域戴维南等效支路替代,进而将子模块等效为一个戴维南等效支路,最终将子模块戴维南等效支路级联构成整个桥臂的戴维南等效支路。这种算法大大减少了MMC导纳矩阵的维数,加快了仿真速度,但无法对二极管的插值进行精确地处理,导致子模块闭锁时系统仿真特性与真实特性出现较大偏差。此外,现有研究一般仅均针对半桥子模块,对于全桥子模块的研究较少。全桥子模块不同于半桥子模块,其闭锁特性与闭锁方式也有很大不同。因此,有必要提出一种考虑全桥子模块闭锁特性的快速仿真建模方法。
发明内容
针对现有技术所存在的上述技术问题,本发明提供了一种考虑多种闭锁模态的基于全桥子模块MMC的等效仿真方法,该方法将子模块等效为戴维南等效支路,同时引入真实二极管器件,能够在保证仿真速度的前提下对全桥子模块各闭锁状态进行精确模拟。
一种考虑多种闭锁模态的基于全桥子模块MMC的等效仿真方法,所述的MMC为三相六桥臂结构,每个桥臂均由若干个全桥子模块级联构成;所述的全桥子模块由四个带反并联二极管的IGBT管T1~T4、一个均压电阻和一个电容组成;其中,IGBT管T1的集电极与电容的正极、均压电阻的一端以及IGBT管T3的集电极相连,IGBT管T1的发射极与IGBT管T2的集电极相连并构成全桥子模块的正极,IGBT管T2的发射极与电容的负极、均压电阻的另一端以及IGBT管T4的发射极相连,IGBT管T3的发射极与IGBT管T4的集电极相连并构成全桥子模块的负极;四个IGBT管T1~T4的基极均接收外部设备提供的开关信号;所述的等效仿真方法包括如下步骤:
(1)获取MMC的桥臂运行参数,所述的桥臂运行参数包括桥臂电流和桥臂各子模块内IGBT管的开关状态;
(2)根据所述的桥臂运行参数构建MMC的桥臂等效电路,并确定桥臂等效电路中各器件的参数;
所述的桥臂等效电路由四个等效电阻R1~R4、四个等效电压源U1~U4以及六个等效二极管D1~D6构成;其中,等效电压源U1的正极对应为桥臂等效电路的正极,等效电压源U1的负极与等效电阻R1的一端相连,等效电阻R1的另一端与等效二极管D1的阳极和等效二极管D2的阴极相连,等效二极管D1的阴极与等效电压源U2的正极和等效二极管D3的阴极相连,等效电压源U2的负极与等效电阻R2的一端相连,等效电阻R2的另一端与等效二极管D2的阳极和等效二极管D4的阳极相连,等效二极管D3的阳极与等效二极管D4的阴极、等效电压源U3的正极以及等效电压源U4的负极相连,等效电压源U3的负极与等效电阻R3的一端相连,等效电压源U4的正极与等效电阻R4的一端相连,等效电阻R3的另一端与等效二极管D5的阳极相连,等效电阻R4的另一端与等效二极管D6的阴极相连,等效二极管D5的阴极与等效二极管D6的阳极相连并构成桥臂等效电路的负极;
(3)根据所述的桥臂等效电路建立MMC的仿真系统,并对该系统进行仿真。
进一步地,所述等效二极管D1~D6的通态电阻均为1×10-6Ω,关断电阻均为1×109Ω。所述的等效二极管为理想二极管,其通态电阻很小,可忽略不计,其关断电阻为一个很大的数值,保证了器件关断后通过极小的电流。
进一步地,所述等效电阻R1~R4的电阻值以及等效电压源U1~U4的电压值的表达式如下:
其中:Reqi(t)为t时刻等效电阻Ri的电阻值,Rsmji(t)为t时刻桥臂中第j个子模块对应等效电阻Ri的等效电阻值,ueqi(t)为t时刻等效电压源Ui的电压值,usmji(t)为t时刻桥臂中第j个子模块对应等效电压源Ui的等效电压值,i为自然数且1≤i≤4,N为桥臂的子模块级联个数,t表示时间。
进一步地,对于桥臂中的任一子模块,当该子模块处于正常运行状态下,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm2(t)=0
usm2(t)=0
Rsm3(t)=0
usm3(t)=0
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
当该子模块处于全闭锁运行状态下,即该子模块中的四个IGBT管T1~T4均关断,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm3(t)=0
usm3(t)=0
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
当该子模块处于左半部分闭锁运行状态下,即该子模块中的IGBT管T1和T2关断;若该子模块在闭锁前的输出电压为正电平或零电平2,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usmj2(t)=0
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
若该子模块在闭锁前的输出电压为负电平或零电平1,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usm2(t)=0
Rsm3(t)=0
usm3(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
当该子模块处于右半部分闭锁运行状态下,即该子模块中的IGBT管T3和T4关断;若该子模块在闭锁前的输出电压为正电平或零电平1,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usm2(t)=0
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
若该子模块在闭锁前的输出电压为负电平或零电平2,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usm2(t)=0
Rsm3(t)=0
usm3(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
其中:Rc为子模块电容的等效阻值,Rb为子模块均压电阻的阻值,ΔT为仿真步长,C0为子模块电容的容值,ic(t)为t时刻流经子模块电容的电流值,ic(t-ΔT)为t-ΔT时刻流经子模块电容的电流值,uc(t)为t时刻子模块电容的电压值,uc(t-ΔT)为t-ΔT时刻子模块电容的电压值,uceq(t-ΔT)为t-ΔT时刻子模块电容的等效历史电压源幅值,R1(t)为t时刻子模块中IGBT管T1的等效电阻值,R2(t)为t时刻子模块中IGBT管T2的等效电阻值,R3(t)为t时刻子模块中IGBT管T3的等效电阻值,R4(t)为t时刻子模块中IGBT管T4的等效电阻值,RA(t)、RE(t)、RM(t)和RN(t)均为中间变量且均由公式给定,iarm(t)为t时刻的桥臂电流,iblk2(t)为t时刻流经等效电阻R2的电流,iblk3(t)为t时刻流经等效电阻R3的电流,iblk4(t)为t时刻流经等效电阻R4的电流。
当t时刻子模块中的IGBT管T1导通,则R1(t)为0.01Ω,IGBT管T1关断,则R1(t)为1MΩ;当t时刻子模块中的IGBT管T2导通,则R2(t)为0.01Ω,IGBT管T2关断,则R2(t)为1MΩ;当t时刻子模块中的IGBT管T3导通,则R3(t)为0.01Ω,IGBT管T3关断,则R3(t)为1MΩ;当t时刻子模块中的IGBT管T4导通,则R4(t)为0.01Ω,IGBT管T4关断,则R4(t)为1MΩ。
若子模块的输出电压为正电平,则子模块中IGBT管T1和T4导通,IGBT管T2和T3关断;若子模块的输出电压为负电平,则子模块中IGBT管T1和T4关断,IGBT管T2和T3导通;若子模块的输出电压为零电平1,则子模块中IGBT管T2和T4关断,IGBT管T1和T3导通;若子模块的输出电压为零电平2,则子模块中IGBT管T2和T4导通,IGBT管T1和T3关断。
本发明MMC等效仿真方法利用数值计算以及电磁仿真的内部机理,能够有效地将全桥型子模块以及由该类型子模块构成的的桥臂等效为由受控电压源、可调电阻以及二极管构成的简单支路,大大减少了系统的节点数及相应的仿真运算量;本发明不仅能够大幅度地提升MMC的仿真速度,而且考虑了全桥子模块的各种闭锁方式,本方法搭建的MMC仿真模型在正常运行以及闭锁时均具备有很高的仿真精度。
附图说明
图1为MMC的拓扑结构示意图。
图2为全桥型子模块的结构示意图。
图3为本发明桥臂等效电路的结构示意图。
图4为本发明仿真系统的结构示意图。
图5为本发明等效模型与真实模型关于阀侧电压的对比示意图。
图6为本发明等效模型与真实模型关于阀侧电流的对比示意图。
图7为本发明等效模型与真实模型关于有功功率的对比示意图。
图8为本发明等效模型与真实模型关于子模块电容电压的对比示意图。
具体实施方式
为了更为具体地描述本发明,下面结合附图及具体实施方式对本发明的技术方案及其相关原理进行详细说明。
本实施方式中,MMC为三相全桥模块化多电平换流器。如图1所示,MMC包含有三个相,每个相有两个桥臂,每个桥臂由若干个全桥子模块构成。子模块为全桥型子模块,如图2所示,由四个IGBT(T1、T2、T3、T4)、四个二极管(D1、D2、D3、D4)、一个均压电阻(Rb)和一个电容(C)组成。其中,均压电阻(Rb)与子模块电容(C)并联;第一IGBT(T1)的发射极与第一二极管(D1)的阳极,第二IGBT(T2)的集电极、第二二极管(D2)的阴极相连并构成子模块的一端;第三IGBT(T3)的发射极与第三二极管的阳极(D3),第四IGBT(T4)的集电极、第四二极管(D4)的阴极相连并构成子模块的另一端。IGBT的基极接收外部设备提供的开关信号。
全桥子模块可以输出正电平、负电平以及零电平。当第一、第四IGBT导通,第二、第三IGBT关断时,子模块端口输出电压为正,记为正电平;当第一、第四IGBT关断,第二、第三IGBT导通时,子模块端口输出电压为负,记为负电平;当第一、第三IGBT导通,第二、第四IGBT关断时,子模块端口输出电压为零,记为零电平1;当第一、第三IGBT关断,第二、第四IGBT导通时,子模块端口输出电压为零,记为零电平2。
本发明考虑多种闭锁模态的基于全桥子模块MMC的等效仿真方法,包括如下步骤:
(1)获取对应桥臂的桥臂电流、桥臂各全桥子模块IGBT的开关信号。
获取桥臂的运行参数包括桥臂输入电流Iarm、桥臂所有子模块第一开关信号Φ1、第二开关信号Φ2、第三开关信号Φ3、第四开关信号Φ4;第一开关信号Φ1、第二开关信号Φ2、第三开关信号Φ3、第四开关信号Φ4分别为子模块中T1、T2、T3、T4接收的开关信号。
(2)根据运行参数构建桥臂的等效电路。
如图3所示,由全桥子模块组成的桥臂等效电路由四个等效电阻Req1、Req2、Req3、Req4,四个等效电压源ueq1、ueq2、ueq3、ueq4以及6个等效二极管D、D、D、D、D、D构成。其中第一等效电压源的正极对应等效桥臂的一端,其负极与第一等效电阻相连。第一等效电阻的另一端与第一等效二极管的阳极和第二等效二极管的阴极相连,第一等效二极管的阴极与第二等效电压源的正极和第三等效二极管的阴极相连,第二等效电压源的负极与第二等效电阻的一端相连,第二等效电阻的另一端与第二等效二极管的阳极和第四等效二极管的阳极相连。第三等效二极管的阳极与第四等效二极管的阴极、第五等效二极管的阳极和第六等效二极管的阴极相连。第五等效二极管的阴极与第三等效电压源的正极相连,第三等效电压源的负极与第三等效电阻的一端相连;第六等效二极管的阳极与第四等效电压源的负极相连,第四等效电压源的正极与第四等效电阻的一端相连;第三等效电阻的另一端与第四等效电阻的另一端相连并对应桥臂的另一端。桥臂等效电路可以真实地模拟稳态以及闭锁状态下的真实桥臂,并且所使用的原件简单,便于后续的求解。
等效二极管的通态电阻为1e-6ohm,关断电阻为1e9ohm。等效二极管为理想二极管,其通态电阻很小,可忽略不计,其关断电阻为一个很大的数值,保证了器件关断后通过极小的电流。
等效桥臂中各子模块电容可由一个等效的历史电压源及一个与电容大小、仿真步长有关的等效电阻来替代,计算公式如下所示:
ucj(t)=Rcicj(t)+ucjeq(t)
其中:RC为子模块电容的等效阻值,ΔT为仿真步长,C0为子模块电容的容值;ucjeq(t)为t时刻第j个子模块电容的历史电压源,icj(t)为t时刻第j个子模块的电容电流,ucj(t)为t时刻第j个子模块电容的电压值。
等效桥臂中的等效电阻与等效电压源的计算公式如下所示:
其中:Req1、Req2、Req3、Req4分别为第一、二、三、四等效电阻的等效阻值,ueq1(t)、ueq2(t)、ueq3(t)、ueq4(t)分别为第一、二、三、四等效电压源的等效电压值,Rsmj1、Rsmj2、Rsmj3、Rsmj4分别为第j个子模块中对应第一、二、三、四等效电阻的等效阻值,usmj1(t)、usmj2(t)、usmj3(t)、usmj4(t)分别为第j个子模块中对应第一、二、三、四等效电压源的等效电压值。
当MMC中第j个子模块处于正常运行状态下,等效桥臂计算公式如下所示:
Rsmj2=0
usmj2(t)=0
Rsmj3=0
usmj3(t)=0
Rsmj4=0
usmj4(t)=0
另有:
其中:iarm(t)为桥臂的输入电流;R1j为第j个子模块的第一IGBT以及第一二极管的等效电阻,R2j为第j个子模块的第二IGBT以及第二二极管的等效电阻,R3j为第j个子模块的第三IGBT以及第三二极管的等效电阻,R4j为第j个子模块的第四IGBT以及第四二极管的等效电阻,Rbj为第j个子模块的均压电阻的阻值。电容的等效电阻的阻值和等效电压源的电压值的计算公式可以将子模块中的电容等效为一个戴维南等效支路,为后续的仿真电路的等效提供了基础。
当第j个子模块中的第一IGBT导通时,R1j的取值为0.01Ω;关断时,R1j的取值为1MΩ。当第j个子模块中的第二IGBT导通时,R2j的取值为0.01Ω;关断时,R2j的取值为1MΩ。当第j个子模块中的第三IGBT导通时,R3j的取值为0.01Ω;关断时,R3j的取值为1MΩ。当第j个子模块中的第四IGBT导通时,R4j的取值为0.01Ω;关断时,R4j的取值为1MΩ。
当MMC第j个子模块处于全闭锁运行状态下,此时该子模块第一、第二、第三、第四IGBT全部处于闭锁关断状态。等效桥臂计算公式如下所示:
Rsmj1=0.02
usmj1(t)=0
Rsmj3=0
usmj3(t)=0
Rsmj4=0
usmj4(t)=0
另有:
其中,iblk2为流经第二等效电阻的电流。
当MMC第j个子模块处于左半部分闭锁运行状态下,此时该子模块第一、第二IGBT全部处于闭锁关断状态。
若子模块初始触发状态为正电平或零电平2,则:
Rsmj1=0.02
usmj1(t)=0
Rsmj2=0
usmj2(t)=0
Rsmj4=0
usmj4(t)=0
另有:
其中,iblk3为流经第三等效电阻的电流。
若子模块初始触发状态为负电平或零电平1,则:
Rsmj1=0.02
usmj1(t)=0
Rsmj2=0
usmj2(t)=0
Rsmj3=0
usmj3(t)=0
另有:
其中,iblk4为流经第四等效电阻的电流。
当MMC第j个子模块处于右半部分闭锁运行状态下,此时该子模块第三、第四IGBT全部处于闭锁关断状态。
若子模块初始触发状态为正电平或零电平1,则:
Rsmj1=0.02
usmj1(t)=0
Rsmj2=0
usmj2(t)=0
Rsmj4=0
usmj4(t)=0
另有:
其中,iblk3为流经第三等效电阻的电流。
若子模块初始触发状态为负电平或零电平2,则:
Rsmj1=0.02
usmj1(t)=0
Rsmj2=0
usmj2(t)=0
Rsmj3=0
usmj3(t)=0
另有:
其中,iblk4为流经第四等效电阻的电流。
(3)根据桥臂等效电路建立相应的仿真系统,并进行仿真。
根据步骤(1)至(2),遍历MMC中的所有桥臂;从而得到MMC的快速仿真电路。以下利用电力系统电磁暂态仿真软件对MMC的仿真系统进行仿真,该仿真系统如图4所示,其中PPC-1和PPC-2分别为两侧交流系统的公共负荷连接点,该仿真系统的仿真参数如表1所示:
表1
仿真场景:在1.5s时,直流线路中间发生双极短路故障,1.502s时两端换流站所有子模块同时闭锁。分别使用真实桥臂模型以及本发明提出的等效建模方法搭建了模型进行仿真,图5为阀侧电压的对比图,图6为直流电压的对比图,图7为桥臂电流的对比图,图8为子模块电容电压的对比图。从仿真波形可见,本发明桥臂等效模型和真实桥臂模型关于主要电气运行参数(如阀侧电压、直流电压、桥臂电流及子模块电容电压)几乎一致,故也验证了本发明等效建模方法的有效性。

Claims (6)

1.一种考虑多种闭锁模态的基于全桥子模块MMC的等效仿真方法,所述的MMC为三相六桥臂结构,每个桥臂均由若干个全桥子模块级联构成;所述的全桥子模块由四个带反并联二极管的IGBT管T1~T4、一个均压电阻和一个电容组成;其中,IGBT管T1的集电极与电容的正极、均压电阻的一端以及IGBT管T3的集电极相连,IGBT管T1的发射极与IGBT管T2的集电极相连并构成全桥子模块的正极,IGBT管T2的发射极与电容的负极、均压电阻的另一端以及IGBT管T4的发射极相连,IGBT管T3的发射极与IGBT管T4的集电极相连并构成全桥子模块的负极;四个IGBT管T1~T4的基极均接收外部设备提供的开关信号;所述的等效仿真方法包括如下步骤:
(1)获取MMC的桥臂运行参数,所述的桥臂运行参数包括桥臂电流和桥臂各子模块内IGBT管的开关状态;
(2)根据所述的桥臂运行参数构建MMC的桥臂等效电路,并确定桥臂等效电路中各器件的参数;
所述的桥臂等效电路由四个等效电阻R1~R4、四个等效电压源U1~U4以及六个等效二极管D1~D6构成;其中,等效电压源U1的正极对应为桥臂等效电路的正极,等效电压源U1的负极与等效电阻R1的一端相连,等效电阻R1的另一端与等效二极管D1的阳极和等效二极管D2的阴极相连,等效二极管D1的阴极与等效电压源U2的正极和等效二极管D3的阴极相连,等效电压源U2的负极与等效电阻R2的一端相连,等效电阻R2的另一端与等效二极管D2的阳极和等效二极管D4的阳极相连,等效二极管D3的阳极与等效二极管D4的阴极、等效电压源U3的正极以及等效电压源U4的负极相连,等效电压源U3的负极与等效电阻R3的一端相连,等效电压源U4的正极与等效电阻R4的一端相连,等效电阻R3的另一端与等效二极管D5的阳极相连,等效电阻R4的另一端与等效二极管D6的阴极相连,等效二极管D5的阴极与等效二极管D6的阳极相连并构成桥臂等效电路的负极;
(3)根据所述的桥臂等效电路建立MMC的仿真系统,并对该系统进行仿真。
2.根据权利要求1所述的等效仿真方法,其特征在于:所述等效二极管D1~D6的通态电阻均为1×10-6Ω,关断电阻均为1×109Ω。
3.根据权利要求1所述的等效仿真方法,其特征在于:所述等效电阻R1~R4的电阻值以及等效电压源U1~U4的电压值的表达式如下:
<mfenced open = '' close = ''> <mtable> <mtr> <mtd> <mrow> <msub> <mi>R</mi> <mrow> <mi>e</mi> <mi>q</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>m</mi> <mi>j</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <msub> <mi>u</mi> <mrow> <mi>e</mi> <mi>q</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mi>m</mi> <mi>j</mi> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中:Reqi(t)为t时刻等效电阻Ri的电阻值,Rsmji(t)为t时刻桥臂中第j个子模块对应等效电阻Ri的等效电阻值,ueqi(t)为t时刻等效电压源Ui的电压值,usmji(t)为t时刻桥臂中第j个子模块对应等效电压源Ui的等效电压值,i为自然数且1≤i≤4,N为桥臂的子模块级联个数,t表示时间。
4.根据权利要求3所述的等效仿真方法,其特征在于:对于桥臂中的任一子模块,当该子模块处于正常运行状态下,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
<mrow> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>E</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> 1
<mrow> <msub> <mi>R</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mo>(</mo> <mrow> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> <mo>+</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> <mo>+</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>R</mi> <mi>E</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>R</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow>
<mrow> <msub> <mi>R</mi> <mi>M</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>N</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>R</mi> <mi>N</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>+</mo> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>&amp;times;</mo> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mo>(</mo> <mi>t</mi> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow>
Rsm2(t)=0
usm2(t)=0
Rsm3(t)=0
usm3(t)=0
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
<mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mi>A</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mi>E</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>i</mi> <mrow> <mi>a</mi> <mi>r</mi> <mi>m</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>R</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> </mfrac> <mo>&amp;CenterDot;</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>N</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
当该子模块处于全闭锁运行状态下,即该子模块中的四个IGBT管T1~T4均关断,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
<mrow> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
Rsm3(t)=0
usm3(t)=0
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
<mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>i</mi> <mrow> <mi>b</mi> <mi>l</mi> <mi>k</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
当该子模块处于左半部分闭锁运行状态下,即该子模块中的IGBT管T1和T2关断;若该子模块在闭锁前的输出电压为正电平或零电平2,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usmj2(t)=0
<mrow> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
<mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>i</mi> <mrow> <mi>b</mi> <mi>l</mi> <mi>k</mi> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
若该子模块在闭锁前的输出电压为负电平或零电平1,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usm2(t)=0
Rsm3(t)=0
usm3(t)=0
<mrow> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>4</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>4</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
<mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>i</mi> <mrow> <mi>b</mi> <mi>l</mi> <mi>k</mi> <mn>4</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
当该子模块处于右半部分闭锁运行状态下,即该子模块中的IGBT管T3和T4关断;若该子模块在闭锁前的输出电压为正电平或零电平1,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usm2(t)=0
<mrow> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
Rsm4(t)=0
usm4(t)=0
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
<mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>i</mi> <mrow> <mi>b</mi> <mi>l</mi> <mi>k</mi> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
若该子模块在闭锁前的输出电压为负电平或零电平2,则t时刻其对应等效电阻Ri的等效电阻值Rsmi(t)以及对应等效电压源Ui的等效电压值usmi(t)的计算公式如下:
Rsm1(t)=0.02
usm1(t)=0
Rsm2(t)=0
usm2(t)=0
Rsm3(t)=0
usm3(t)=0
<mrow> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>4</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>u</mi> <mrow> <mi>s</mi> <mi>m</mi> <mn>4</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
uceq(t-ΔT)=Rcic(t-ΔT)+uc(t-ΔT)
uc(t)=Rcic(t)+uceq(t-ΔT)
<mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>T</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>C</mi> <mn>0</mn> </msub> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>i</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>i</mi> <mrow> <mi>b</mi> <mi>l</mi> <mi>k</mi> <mn>4</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>c</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;Delta;</mi> <mi>T</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>R</mi> <mi>c</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow>
其中:Rc为子模块电容的等效阻值,Rb为子模块均压电阻的阻值,ΔT为仿真步长,C0为子模块电容的容值,ic(t)为t时刻流经子模块电容的电流值,ic(t-ΔT)为t-ΔT时刻流经子模块电容的电流值,uc(t)为t时刻子模块电容的电压值,uc(t-ΔT)为t-ΔT时刻子模块电容的电压值,uceq(t-ΔT)为t-ΔT时刻子模块电容的等效历史电压源幅值,R1(t)为t时刻子模块中IGBT管T1的等效电阻值,R2(t)为t时刻子模块中IGBT管T2的等效电阻值,R3(t)为t时刻子模块中IGBT管T3的等效电阻值,R4(t)为t时刻子模块中IGBT管T4的等效电阻值,RA(t)、RE(t)、RM(t)和RN(t)均为中间变量且均由公式给定,iarm(t)为t时刻的桥臂电流,iblk2(t)为t时刻流经等效电阻R2的电流,iblk3(t)为t时刻流经等效电阻R3的电流,iblk4(t)为t时刻流经等效电阻R4的电流。
5.根据权利要求4所述的等效仿真方法,其特征在于:当t时刻子模块中的IGBT管T1导通,则R1(t)为0.01Ω,IGBT管T1关断,则R1(t)为1MΩ;当t时刻子模块中的IGBT管T2导通,则R2(t)为0.01Ω,IGBT管T2关断,则R2(t)为1MΩ;当t时刻子模块中的IGBT管T3导通,则R3(t)为0.01Ω,IGBT管T3关断,则R3(t)为1MΩ;当t时刻子模块中的IGBT管T4导通,则R4(t)为0.01Ω,IGBT管T4关断,则R4(t)为1MΩ。
6.根据权利要求4所述的等效仿真方法,其特征在于:若子模块的输出电压为正电平,则子模块中IGBT管T1和T4导通,IGBT管T2和T3关断;若子模块的输出电压为负电平,则子模块中IGBT管T1和T4关断,IGBT管T2和T3导通;若子模块的输出电压为零电平1,则子模块中IGBT管T2和T4关断,IGBT管T1和T3导通;若子模块的输出电压为零电平2,则子模块中IGBT管T2和T4导通,IGBT管T1和T3关断。
CN201510516122.0A 2015-08-21 2015-08-21 一种考虑多种闭锁模态的基于全桥子模块mmc的等效仿真方法 Expired - Fee Related CN105117543B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510516122.0A CN105117543B (zh) 2015-08-21 2015-08-21 一种考虑多种闭锁模态的基于全桥子模块mmc的等效仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510516122.0A CN105117543B (zh) 2015-08-21 2015-08-21 一种考虑多种闭锁模态的基于全桥子模块mmc的等效仿真方法

Publications (2)

Publication Number Publication Date
CN105117543A CN105117543A (zh) 2015-12-02
CN105117543B true CN105117543B (zh) 2017-12-08

Family

ID=54665531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510516122.0A Expired - Fee Related CN105117543B (zh) 2015-08-21 2015-08-21 一种考虑多种闭锁模态的基于全桥子模块mmc的等效仿真方法

Country Status (1)

Country Link
CN (1) CN105117543B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106909702A (zh) * 2015-12-23 2017-06-30 加拿大欧泊实时技术有限公司北京代表处 一种模块化多电平整流器及变电站、建模仿真平台、方法
CN106202664A (zh) * 2016-06-30 2016-12-07 全球能源互联网研究院 一种半桥型功率子模块的仿真方法
CN106407616B (zh) * 2016-11-11 2019-08-06 南方电网科学研究院有限责任公司 一种全桥子模块等效仿真方法
CN106682352B (zh) * 2017-01-11 2020-02-21 南方电网科学研究院有限责任公司 一种箝位双子模块等效仿真方法
CN107748313B (zh) * 2017-10-16 2019-12-03 华北电力大学 基于与或逻辑的识别hbsm-mmc内部短路故障的方法
CN108133095A (zh) * 2017-12-14 2018-06-08 广东电网有限责任公司电力科学研究院 一种双半桥子模块mmc建模仿真方法及装置
CN108376992B (zh) * 2018-03-16 2021-04-30 云南电网有限责任公司 一种混合子模块mmc等效仿真方法
CN109787265B (zh) * 2019-01-31 2022-11-01 中国能源建设集团天津电力设计院有限公司 一种平抑新能源出力波动的柔性直流输电系统协调控制方法
CN110674611A (zh) * 2019-09-20 2020-01-10 全球能源互联网集团有限公司 一种锂电池储能型mmc的等效仿真方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995933A (zh) * 2014-06-18 2014-08-20 上海傲亚微电子有限公司 一种新型晶体管小信号等效电路模型
CN104298837A (zh) * 2014-11-12 2015-01-21 东南大学 器件等效电路模型参数提取方法及焊盘寄生参数提取方法
CN104753523A (zh) * 2013-12-25 2015-07-01 上海华虹宏力半导体制造有限公司 带寄生效应的圆形高压场效应管等效电路及仿真方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104753523A (zh) * 2013-12-25 2015-07-01 上海华虹宏力半导体制造有限公司 带寄生效应的圆形高压场效应管等效电路及仿真方法
CN103995933A (zh) * 2014-06-18 2014-08-20 上海傲亚微电子有限公司 一种新型晶体管小信号等效电路模型
CN104298837A (zh) * 2014-11-12 2015-01-21 东南大学 器件等效电路模型参数提取方法及焊盘寄生参数提取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DC-Side Harmonic Currents Calculation and DC-Loop Resonance Analysis for an LCC–MMC Hybrid HVDC Transmission System;Zheren Zhang 等;《IEEE TRANSACTIONS ON POWER DELIVERY》;20150430;第30卷(第2期);第642-651页 *
Impacts of Three MMC-HVDC Configurations on AC System Stability Under DC Line Faults;Geng Tang 等;《IEEE TRANSACTIONS ON POWER SYSTEMS》;20141130;第29卷(第6期);第3030-3040页 *

Also Published As

Publication number Publication date
CN105117543A (zh) 2015-12-02

Similar Documents

Publication Publication Date Title
CN105117543B (zh) 一种考虑多种闭锁模态的基于全桥子模块mmc的等效仿真方法
CN105117549B (zh) 一种考虑多种闭锁模态的基于箝位双子模块mmc的等效仿真方法
CN102663174B (zh) 一种mmc的仿真方法及应用
CN103746583B (zh) 一种mmc的桥臂等效仿真方法
CN108229021B (zh) 基于实时数字仿真器的模块化多电平换流器闭锁建模方法
CN102969888B (zh) 基于rtds的mmc自定义子模块的设计方法
CN103914599B (zh) 一种模块化多电平换流器的戴维南等效整体建模方法
CN103605850B (zh) 一种带子模块闭锁功能的mmc等效建模方法
CN103296885B (zh) Pwm功率变换器模拟仿真电路及仿真方法
CN105045952B (zh) 模块化多电平换流器多维度建模方法与仿真方法
CN103593520A (zh) 一种模块化多电平换流器的等值仿真计算建模方法
CN103336866A (zh) 一种电磁暂态仿真中含负电阻支路的处理方法
CN107942729A (zh) 基于实时数字仿真器的模块化多电平换流器桥臂建模方法
CN106909702A (zh) 一种模块化多电平整流器及变电站、建模仿真平台、方法
CN108133095A (zh) 一种双半桥子模块mmc建模仿真方法及装置
CN104950694A (zh) 一种联合rtds和rt-lab的mmc仿真系统
CN103729502B (zh) 一种提高电力系统电磁暂态仿真速度的方法
CN103956741B9 (zh) 一字链及支链式的三相对称多电源非环网电力系统直算法
CN107257205B (zh) 一种mmc功率模块非线性特征仿真模型
CN108376992A (zh) 一种混合子模块mmc等效仿真方法
CN108460181A (zh) 一种适用于半、全桥子模块混合型mmc戴维南等效模型的线性均压排序算法
CN111898282B (zh) 一种改进的模块化多电平换流器戴维南等效建模方法
CN106407616B (zh) 一种全桥子模块等效仿真方法
CN103279626B (zh) 一种仿真模型建立方法及系统
CN113158617B (zh) 电压源型换流器电磁暂态仿真通用解耦方法与系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Xu Zheng

Inventor after: Liu Gaoren

Inventor after: Li Haiying

Inventor after: Xiao Huangqing

Inventor after: Tang Geng

Inventor before: Xu Zheng

Inventor before: Liu Gaoren

Inventor before: Xiao Huangqing

Inventor before: Tang Geng

COR Change of bibliographic data
TA01 Transfer of patent application right

Effective date of registration: 20151231

Address after: 310027 Hangzhou, Zhejiang Province, Xihu District, Zhejiang Road, No. 38, No.

Applicant after: Zhejiang University

Applicant after: State Grid Corporation of China

Applicant after: Nanrui Relay Protection Electricity Co., Ltd., Nanjing

Address before: 310027 Hangzhou, Zhejiang Province, Xihu District, Zhejiang Road, No. 38, No.

Applicant before: Zhejiang University

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20171212

Address after: 310027 Hangzhou, Zhejiang Province, Xihu District, Zhejiang Road, No. 38, No.

Co-patentee after: State Grid Corporation of China

Patentee after: Zhejiang University

Co-patentee after: Nanrui Relay Protection Electricity Co., Ltd., Nanjing

Co-patentee after: State Grid Zhejiang Electric Power Company

Address before: 310027 Hangzhou, Zhejiang Province, Xihu District, Zhejiang Road, No. 38, No.

Co-patentee before: State Grid Corporation of China

Patentee before: Zhejiang University

Co-patentee before: Nanrui Relay Protection Electricity Co., Ltd., Nanjing

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171208

Termination date: 20190821