CN112199914B - 一种电力电子开关恒导纳模型建立方法及系统 - Google Patents
一种电力电子开关恒导纳模型建立方法及系统 Download PDFInfo
- Publication number
- CN112199914B CN112199914B CN202011036925.3A CN202011036925A CN112199914B CN 112199914 B CN112199914 B CN 112199914B CN 202011036925 A CN202011036925 A CN 202011036925A CN 112199914 B CN112199914 B CN 112199914B
- Authority
- CN
- China
- Prior art keywords
- equivalent
- circuit
- value
- power electronic
- inductance
- 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 93
- 230000010354 integration Effects 0.000 claims abstract description 31
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims description 61
- 230000010355 oscillation Effects 0.000 abstract description 20
- 238000004422 calculation algorithm Methods 0.000 abstract description 11
- 230000009471 action Effects 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 8
- 238000013016 damping Methods 0.000 abstract description 6
- 238000004088 simulation Methods 0.000 description 23
- 230000001052 transient effect Effects 0.000 description 14
- 230000008569 process Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 7
- 230000004044 response Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 239000003990 capacitor Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000011960 computer-aided design Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/33—Design verification, e.g. functional simulation or model checking
- G06F30/3308—Design verification, e.g. functional simulation or model checking using simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Computing Systems (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开的一种电力电子开关恒导纳模型建立方法及系统,所述建立方法,首先在导通和关断等效支路中分别引入电阻,形成具有阻尼作用的一阶动态电路形式,进而通过终值定理及对连续时间系统的特征根分析确定等效电路的最优参数,提高精度,最后采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型,因指数积分法具备的L‑稳定性,可有效抑制开关动作时的数值振荡,因指数积分法具备的高阶精度,又能实现开关动作前后仿真算法的一致性,避免了状态切换带来的复杂性,本发明提供了一种精度高、数值振荡小且仿真简单的电力电子开关恒导纳开关模型建立方法。
Description
技术领域
本发明涉及电力电子开关性能分析技术领域,特别涉及一种电力电子开关恒导纳模型建立方法及系统。
背景技术
近年来,大规模新能源发电、高压直流输电技术,以及柔性交流输电(FACTS)技术快速发展,导致电网中电力电子装置的比例快速提高,当前电网的电力电子化趋势也越来越明显。包含大量电力电子设备的换流系统的准确高效以及实时仿真对电网的研究非常重要。作为换流器的核心器件,开关模型的暂态误差以及数值积分算法的选取关系着电磁暂态仿真的效率,因此研究换流器开关模型有着重要的意义。
在电磁暂态仿真算法中,后退欧拉法虽然不会带来非原生数值振荡,但是计算精度较低。常用的隐式梯形积分法一般仅具有A-稳定,开关动作时存在因非状态量突变而引起的非原生数值振荡问题。因此,现有的仿真过程中,一般在没有开关动作时采用隐式梯形积分法,而在动作时切换为后退欧拉法,但这又增加了仿真的复杂性。
一种基于响应匹配的换流器恒导纳建模方法,该方法基于连续和离散系统零、极点响应匹配原理,采用含参数化历史电流源的恒导纳开关模型,并将其视为线性离散时间系统进行分析。但是该方法中,历史电流由电压和电流两项构成,在开关导通和关断时会存在因非状态量突变而引起的非原生数值振荡(类似隐式梯形积分存在的问题)。
一种基于参数化历史电流源的广义小步长模型,通过参数设置消除了开关模型初始暂态误差,较好地解决了虚拟功率损耗的问题,但是该模型没有解决非原生数值振荡问题。
基于传统小步长开关模型,提出采用合理参数设置的方法来消除虚拟损耗的问题。但是由于该模型的支路差分方程存在非状态量,无法抑制数值振荡。
降采样方法阻尼电流振荡,但其中的恒导纳模型存在严重的虚拟功率损耗问题,其应用范围大大受限。
响应匹配方法,实现步骤如下:
(1)提出开关模型支路组合形式;
(2)得到离散时间系统下的差分方程;
(3)利用终值定理匹配开关模型稳态特性;
(4)通过谱半径分析匹配开关模型暂态特性;
响应匹配的方法基于连续和离散系统零、极点响应匹配原理,采用含参数化历史电流源的恒导纳开关模型,并将其视为线性离散时间系统进行分析。但是该方法中,历史电流由电压和电流两项构成,在开关导通和关断时会存在因非状态量突变而引起的非原生数值振荡(类似隐式梯形积分存在的问题)。
如何提供一种精度高、数值振荡小且仿真简单的电力电子开关恒导纳开关模型建立方法成为一个亟待解决的技术问题。
发明内容
本发明的目的是提供一种电力电子开关恒导纳模型建立方法及系统,以提供一种精度高、数值振荡小且仿真简单的电力电子开关恒导纳开关模型建立方法。
为实现上述目的,本发明提供了如下方案:
一种电力电子开关恒导纳模型建立方法,所述建立方法包括如下步骤:
将电力电子开关电路中的闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路;
基于KCL定律,建立所述等效电路的微分方程;
确定所述微分方程的状态矩阵;
根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值;
根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型。
可选的,所述基于KCL定律,建立所述等效电路的微分方程,具体包括:
基于KCL定律,建立所述等效电路的微分方程为:
其中,i1(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的上桥臂的电流,u2(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的下桥臂的电压,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值,Ud1表示上桥臂输入电压,Ud2表示下桥臂输入电压,io表示负载电流。
可选的,所述确定所述微分方程的状态矩阵,具体包括:
确定所述微分方程的状态矩阵A为:
其中,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值。
可选的,等效电路中与等效电感元件串联的等效电阻元件的电阻值为1mΩ,等效电路中与等效电容元件并联的等效电阻元件的电阻值为1MΩ。
可选的,所述根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值,具体包括:
以等效电容元件的电容值和等效电感元件的电感值为变量,建立状态矩阵的特征根分布等高线;
根据所述特征根分布等高线,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值。
一种电力电子开关恒导纳模型建立系统,所述建立系统包括:
等效电路建立模块,用于将电力电子开关电路中的闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路;
微分方程建立模块,用于基于KCL定律,建立所述等效电路的微分方程;
状态矩阵建立模块,用于确定所述微分方程的状态矩阵;
数值求解模块,用于根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值;
导纳模型建立模块,用于根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型。
可选的,所述微分方程建立模块,具体包括:
微分方程建立子模块,用于基于KCL定律,建立所述等效电路的微分方程为:
其中,i1(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的上桥臂的电流,u2(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的下桥臂的电压,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值,Ud1表示上桥臂输入电压,Ud2表示下桥臂输入电压,io表示负载电流。
可选的,所述状态矩阵建立模块,具体包括:
状态矩阵建立子模块,用于确定所述微分方程的状态矩阵为:
其中,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值。
可选的,等效电路中与等效电感元件串联的等效电阻元件的电阻值为1mΩ,等效电路中与等效电容元件并联的等效电阻元件的电阻值为1MΩ。
可选的,所述数值求解模块,具体包括:
特征根分布等高线建立子模块,用于以等效电容元件的电容值和等效电感元件的电感值为变量,建立状态矩阵的特征根分布等高线;
数值求解子模块,用于根据所述特征根分布等高线,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明公开了一种电力电子开关恒导纳模型建立方法及系统,所述建立方法包括如下步骤:将电力电子开关电路中的闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路;基于KCL定律,建立所述等效电路的微分方程;确定所述微分方程的状态矩阵;根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值;根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型。本发明导通和关断等效支路中分别引入电阻,形成具有阻尼作用的一阶动态电路形式,进而通过终值定理及对连续时间系统的特征根分析确定等效电路的最优参数,提高精度,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型,因指数积分法具备的L-稳定性,可有效抑制开关动作时的数值振荡,因指数积分法具备的高阶精度,又能实现开关动作前后仿真算法的一致性,避免了状态切换带来的复杂性,本发明提供了一种精度高、数值振荡小且仿真简单的电力电子开关恒导纳开关模型建立方法。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种电力电子开关恒导纳模型建立方法的流程图;
图2为本发明提供的差分等效电路图;图2(a)为导通状态时的差分等效电路图,图2(b)为关断状态时的差分等效电路图;
图3为本发明提供的三相桥式电路图;
图4为本发明提供的三相桥式电路中一个桥臂的等效电路图;
图5为本发明提供的一种电力电子开关恒导纳模型建立方法的原理图;
图6为本发明提供的不同算法的仿真结果对比图。
具体实施方式
本发明的目的是提供一种电力电子开关恒导纳模型建立方法及系统,以提供一种精度高、数值振荡小且仿真简单的电力电子开关恒导纳开关模型建立方法。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对发明作进一步详细的说明。
如图1和5所示,本发明提供一种电力电子开关恒导纳模型建立方法,所述建立方法包括如下步骤:
步骤101,将电力电子开关电路中的闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路。
小L/C恒导纳模型因存在着开关器件的暂态过程而更接近电力电子器件的实际工作过程。同时,我们也知道,实际的开关导通或关断后,并不是理想的短路和断路,而是存在着一定的导通和关断电阻。此外,为了阻尼数值振荡,二值电阻模型中有时也会引入缓冲电路(subber电路)。因此,为了更准确的描述开关器件的导通/关断过程及其工作特性,本发明提出一种能够综合二值电阻模型和小L/C恒导纳模型特点的开关器件新型等效支路拓扑,即:通过在小L/C恒导纳模型导通和关断等效支路中分别引入电阻Ron和Roff,形成具有阻尼作用的一阶动态电路形式。
该等效支路中的闭合等效支路中的Ron与L、关断等效支路中的Roff与C分别存在串、并联两种组合形式。而组合形式不同,数值特性也不相同。电磁暂态仿真中,为了抑制非原生数值振荡,通常在开关动作时将仿真算法切换为后退欧拉法。表1为采用后退欧拉法时不同等效支路拓扑的差分方程。
表1采用后退欧拉法时不同等效支路拓扑的差分方程
从表1中可以看出:RL并联和RC串联两种拓扑的历史电流源中,存在着非状态变量,这将会给仿真带来非原生数值振荡,而RL串联、RC并联则没有这一问题。因此,本发明采用RL串联(R+L)、RC并联(R//C)分别作为开关器件导通和关断状态的等效支路拓扑。此外,通过上面的分析也可以看出,snubber缓冲电路电路由于采用的是R与C的串联形式,存在着非原生数值振荡问题。
图2为采用RL串联(R+L)、RC并联(R//C)分别作为开关器件导通和关断状态的等效支路拓扑时开关器件在导通和关断状态时的差分等效电路。其中,图2(a)为导通状态时的差分等效电路图,图2(b)为关断状态时的差分等效电路图。图2中:Yon和Yoff分别为开关导通和关断时的等效差分导纳。
其中,以三相桥式电路中的一个桥臂(图3中虚线框内的桥臂)为例,其预设工作状态为上桥臂导通、下桥臂断开,三相桥式电路如图3所示,预设工作状态的桥臂的等效电路如图4所示。
步骤102,基于KCL定律,建立所述等效电路的微分方程。
假定上、下桥臂的中点电压为Uo。当前电路状态对应的时间系统以电压Ud1、Ud2和负载电流io为输入量,以桥臂电流i1、i2和负载电压Uo为输出量。开关的状态切换过程即为输出响应的暂态收敛过程。
将上桥臂开关闭合、下桥臂开关关断的时域系统定义为系统1,以RL串联支路的电流,RC并联支路的电压为状态量,对图4所示的等效电路应用KCL定律,推导得到所述等效电路的微分方程为:
其中,i1(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的上桥臂的电流,u2(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的下桥臂的电压,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值,Ud1表示上桥臂输入电压,Ud2表示下桥臂输入电压,io表示负载电流。
步骤103,确定所述微分方程的状态矩阵。
对于系统1而言,推导得到状态矩阵A为:
将上桥臂开关关断、下桥臂开关闭合的时域系统定义为系统2。同理,推导可得状态矩阵为:
步骤104,根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值。
时域系统的暂态响应特性取决于系统状态矩阵A的特征根,对于图4所示的等效电路而言,矩阵A的特征根越大,对应系统的暂态过程越短,稳定性越好,开关切换时的虚拟损耗也就越低。
以L、C为自变量,以步长1μs为例,做特征根分布等高线图状态矩阵的特征根越大对应系统暂态过程越快,对应系统暂态过程收敛速度最快,开关切换过程中的虚拟功率损耗最低。
求解EICAM模型的电阻Ron和Roff参数值。
采用R+L形式作为开关导通状态时的等效支路,其时域状态矩阵为:
上式经拉式变换得到其传递函数为:
令输入信号u(t)为一单位阶跃信号,经拉式变换得:
根据终值定理可得:
采用R//C形式作为开关关断状态时的等效支路,其时域状态矩阵为:
上式经拉式变换得到其传递函数为:
令输入信号i(t)为一单位阶跃信号,经拉式变换得:
根据终值定理可得:
由式(7)知,在直流稳态情况下,电感相当于短路,R+L串联支路等效为一纯电阻支路Ron。类似的,由式(11)知,在直流稳态情况下,电容相当于开路,R//C并联支路也可等效为一纯电阻支路Roff。因此,开关器件指数恒导纳模型中导通和关断状态下的Ron和Roff阻值可分别取为二值电阻模型中的阻值,即:Ron=1mΩ,Roff=1MΩ。
步骤105,根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型。
对电阻R+电感L串联支路差分方程的推导:
支路方程为:
其解析解为:
对式(13)离散差分化,其中第二项采用后退欧拉法,得到支路差分方程:
与图2(a)相对应,其中导通态等效导纳Yon以及历史电流为:
电阻R+电感L串联支路差分方程为下式:
推导电阻R//电容C并联支路差分方程:
支路方程为:
其解析解为:
对式(18)离散差分化,其中第二项采用后退欧拉法,得到支路差分方程:
与图2(b)相对应,其中导通态等效导纳Yoff以及历史电流为:
为了提高仿真效率,需使得开关模型状态切换时导通态等效导纳Yon等于关断态等效导纳Yoff,即:
Ysw=Yon=Yoff (22)
电阻R//电容C并联支路差分方程为下式:
理论上,对于一阶动态电路来说,如果其时间常数很小,则其端口电压或支路电流在一个仿真步长内可看作为一个不变的直流量。经过这样处理,对于开关器件来说,一方面,进入稳态后,根据终值定理可以证明其等效支路的阻抗特性与二值电阻模型一致;另一方面,为了减小虚拟损耗而对暂态过程持续时间长短进行的调节,通过仿真步长与等效支路中的R、L、C之间的优化协调配置即可简单实现,与外电路无关。
本发明以图3所示的全控型三相桥式电路作为仿真对象,并在pscad(PSCAD/EMTDC,电磁暂态仿真软件)中搭建相同的算例电路,其中,PSCAD的全拼为Power SystemsComputerAided Design,EMTDC的全拼为Electromagnetic Transients including DC,通过对本发明方法、传统的基于L/C开关模型的仿真结果和在PSCAD中的仿真结果进行对比来验证本发明提出的模型。其中,无特殊说明时仿真步长均采用1μs,电源电压Ud1、Ud2均为50V,交流侧电阻为1Ω、电感为1*10-4H,三角载波频率f=10kHz,调制波频率为50Hz,利用三角载波和正弦波进行比较得出控制上、下桥臂IGBT开断的脉冲信号。
为验证基于指数积分恒导纳(EICAM)模型的有效性,本发明算例采用PSCAD/EMTDC软件中开关采用二值电阻模型以及基于后退欧拉法差分化的传统L/C模型来与EICAM模型进行对比。基于EICAM模型以及传统电容/电感模型的换流器开关S1的电压、电流和功率损耗仿真波形与PSCAD模型的对比结果及其放大图如图6所示。
通过图6可以看出,基于后退欧拉法的传统电感/电容模型精度较低,开关状态切换时的暂态过程较长,粗略估计为十几微秒,开关模型的虚拟功率损耗较高;而本发明提出的EICAM模型电压、电流误差小,暂态过程较短,收敛速度更快,虚拟功率损耗也更低。
下面从理论上分析本发明提供的方法的优点:
1、抑制数值振荡
指数积分恒导纳模型(EICAM)对RL串联、RC并联形成的差分方程历史项中均不含非状态变量,这说明能够有效抑制仿真中的非原型数值振荡。
2、较好地解决了虚拟功率损耗的问题
该模型综合二值电阻模型和传统L/C恒导纳模型的特点,通过在传统L/C恒导纳模型导通和关断等效支路中分别引入电阻Ron、Roff,形成具有阻尼作用的一阶动态电路形式。采用指数积分法相较于传统数值积分算法仿真精度更高,收敛速度更快,通过终值定理及对连续时间系统的特征根分析确定开关模型的最优参数。
3、提高计算精度
以R+L导通等效支路为例,来说明指数积分恒导纳模型的特点。传统积分算法和指数积分法的差分方程均可写成以下形式:
i(t+Δt)=u(t+Δt)/Req+ihis(t+Δt) (23)
式(24)为R+L串联支路采用指数积分法得到的差分方程。
将式(24)中具有指数形式的两项系数进行泰勒展开可得:
式(27)为R+L串联支路采用后退欧拉法得到的差分方程。
分析式(25)(26)可知,泰勒展开式中的一阶项分别与式(27)中的系数一致,所以两者之间存在一阶泰勒近似的关系。
式(28)为R+L串联支路采用隐式梯形法得到的差分方程。
将式(24)中具有指数形式的两项系数采用其他方式进行泰勒展开可得:
分析式(29)(30)可知,第一项系数的泰勒展开式中的一阶项与式(28)中第一项的系数一致,所以两者之间存在一阶泰勒近似的关系,式(28)中第二项的系数与存在一阶泰勒近似的关系,系数出现不完全一致的原因为指数积分法采用欧拉法的格式。
通过上述过程可知,R+L串联支路采用欧拉法、隐式梯形法得到的差分方程各项的系数与指数积分法中具有指数形式的两项系数存在一阶泰勒近似的关系。所以,指数积分法的仿真计算精度相比于前两种方法会有大大地提高,同时仿真计算过程中也会有更快的收敛速度。
基于上述分析,本发明的方法具有如下优点:
1、本发明提出一种基于指数积分的电力电子开关恒导纳模型即EICA模型,该模型综合PSCAD二值电阻模型和传统L/C恒导纳模型的特点,通过在传统L/C恒导纳模型导通和关断等效支路中分别引入电阻Ron、Roff,形成具有阻尼作用的一阶动态电路形式。
2、指数积分恒导纳模型(Exponential Integral Constant Admittance Model,EICAM),因其具备的L-稳定性,可有效抑制开关动作时的数值振荡,同时,因其具备的高阶精度,又能实现开关动作前后仿真算法的一致性,避免了算法切换带来的复杂性。
3、通过终值定理及对连续时间系统的特征根分析确定开关模型的最优参数。
本发明还提供一种电力电子开关恒导纳模型建立系统,所述建立系统包括:
等效电路建立模块,用于将电力电子开关电路中的闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路;
微分方程建立模块,用于基于KCL定律,建立所述等效电路的微分方程。
所述微分方程建立模块,具体包括:
微分方程建立子模块,用于基于KCL定律,建立所述等效电路的微分方程为:
其中,i1(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的上桥臂的电流,u2(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的下桥臂的电压,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值,Ud1表示上桥臂输入电压,Ud2表示下桥臂输入电压,io表示负载电流。
状态矩阵建立模块,用于确定所述微分方程的状态矩阵。
所述状态矩阵建立模块,具体包括:
状态矩阵建立子模块,用于确定所述微分方程的状态矩阵为:
/>
其中,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值。
数值求解模块,用于根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值。
所述数值求解模块,具体包括:特征根分布等高线建立子模块,用于以等效电容元件的电容值和等效电感元件的电感值为变量,建立状态矩阵的特征根分布等高线;数值求解子模块,用于根据所述特征根分布等高线,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值。
其中,等效电路中与等效电感元件串联的等效电阻元件的电阻值为1mΩ,等效电路中与等效电容元件并联的等效电阻元件的电阻值为1MΩ。
导纳模型建立模块,用于根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明公开了一种电力电子开关恒导纳模型建立方法,所述建立方法包括如下步骤:将电力电子开关电路中的闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路;基于KCL定律,建立所述等效电路的微分方程;确定所述微分方程的状态矩阵;根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值;根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型。本发明导通和关断等效支路中分别引入电阻,形成具有阻尼作用的一阶动态电路形式,进而通过终值定理及对连续时间系统的特征根分析确定等效电路的最优参数,提高精度,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型,因指数积分法具备的L-稳定性,可有效抑制开关动作时的数值振荡,因指数积分法具备的高阶精度,又能实现开关动作前后仿真算法的一致性,避免了状态切换带来的复杂性,本发明提供了一种精度高、数值振荡小且仿真简单的电力电子开关恒导纳开关模型建立方法。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
Claims (10)
1.一种电力电子开关恒导纳模型建立方法,其特征在于,所述建立方法包括如下步骤:
将电力电子开关电路中闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路;
基于KCL定律,建立所述等效电路的微分方程;
确定所述微分方程的状态矩阵;
根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值;
根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型,具体包括:
构建RL串联电路的支路方程;
对RL串联电路的支路方程进行解析,获得RL串联电路的解析解;
对RL串联电路的解析解进行离散化,获得RL串联电路的支路差分方程;
求解RL串联电路的支路差分方程,获得预设工作状态下电力电子开关电路中闭合的开关的导通态等效导纳;
构建RC并联电路的支路方程;
对RC并联电路的支路方程进行解析,获得RC并联电路的解析解;
对RC并联电路的解析解进行离散化,获得RC并联电路的支路差分方程;
求解RC并联电路的支路差分方程,获得预设工作状态下电力电子开关电路中断开的开关的导通态等效导纳。
2.根据权利要求1所述的电力电子开关恒导纳模型建立方法,其特征在于,所述基于KCL定律,建立所述等效电路的微分方程,具体包括:
基于KCL定律,建立所述等效电路的微分方程为:
其中,i1(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的上桥臂的电流,u2(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的下桥臂的电压,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值,Ud1表示上桥臂输入电压,Ud2表示下桥臂输入电压,io表示负载电流。
3.根据权利要求1所述的电力电子开关恒导纳模型建立方法,其特征在于,所述确定所述微分方程的状态矩阵,具体包括:
确定所述微分方程的状态矩阵A为:
其中,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值。
4.根据权利要求2或3所述的电力电子开关恒导纳模型建立方法,其特征在于,等效电路中与等效电感元件串联的等效电阻元件的电阻值为1mΩ,等效电路中与等效电容元件并联的等效电阻元件的电阻值为1MΩ。
5.根据权利要求1所述的电力电子开关恒导纳模型建立方法,其特征在于,所述根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值,具体包括:
以等效电容元件的电容值和等效电感元件的电感值为变量,建立状态矩阵的特征根分布等高线;
根据所述特征根分布等高线,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值。
6.一种电力电子开关恒导纳模型建立系统,其特征在于,所述建立系统包括:
等效电路建立模块,用于将电力电子开关电路中闭合的开关等效为RL串联电路,断开的开关等效为RC并联电路,建立预设工作状态的电力电子开关电路的等效电路;
微分方程建立模块,用于基于KCL定律,建立所述等效电路的微分方程;
状态矩阵建立模块,用于确定所述微分方程的状态矩阵;
数值求解模块,用于根据所述状态矩阵,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值;
导纳模型建立模块,用于根据等效电路中等效电容元件的电容值和等效电感元件的电感值,采用指数积分法建立预设工作状态的电力电子开关电路的恒导纳模型,具体包括:
构建RL串联电路的支路方程;
对RL串联电路的支路方程进行解析,获得RL串联电路的解析解;
对RL串联电路的解析解进行离散化,获得RL串联电路的支路差分方程;
求解RL串联电路的支路差分方程,获得预设工作状态下电力电子开关电路中闭合的开关的导通态等效导纳;
构建RC并联电路的支路方程;
对RC并联电路的支路方程进行解析,获得RC并联电路的解析解;
对RC并联电路的解析解进行离散化,获得RC并联电路的支路差分方程;
求解RC并联电路的支路差分方程,获得预设工作状态下电力电子开关电路中断开的开关的导通态等效导纳。
7.根据权利要求6所述的电力电子开关恒导纳模型建立系统,其特征在于,所述微分方程建立模块,具体包括:
微分方程建立子模块,用于基于KCL定律,建立所述等效电路的微分方程为:
其中,i1(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的上桥臂的电流,u2(t)表示三相桥式电路中的一个桥臂的预设工作状态为上桥臂导通下桥臂关断时的下桥臂的电压,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值,Ud1表示上桥臂输入电压,Ud2表示下桥臂输入电压,io表示负载电流。
8.根据权利要求6所述的电力电子开关恒导纳模型建立系统,其特征在于,所述状态矩阵建立模块,具体包括:
状态矩阵建立子模块,用于确定所述微分方程的状态矩阵A为:
其中,L表示等效电路中等效电感元件的电感值、C表示等效电路中等效电容元件的电容值,Ron表示等效电路中与等效电感元件串联的等效电阻元件的电阻值,Roff表示等效电路中与等效电容元件并联的等效电阻元件的电阻值。
9.根据权利要求7或8所述的电力电子开关恒导纳模型建立系统,其特征在于,等效电路中与等效电感元件串联的等效电阻元件的电阻值为1mΩ,等效电路中与等效电容元件并联的等效电阻元件的电阻值为1MΩ。
10.根据权利要求6所述的电力电子开关恒导纳模型建立系统,其特征在于,所述数值求解模块,具体包括:
特征根分布等高线建立子模块,用于以等效电容元件的电容值和等效电感元件的电感值为变量,建立状态矩阵的特征根分布等高线;
数值求解子模块,用于根据所述特征根分布等高线,确定状态矩阵的特征根最大时的等效电路中等效电容元件的电容值和等效电感元件的电感值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011036925.3A CN112199914B (zh) | 2020-09-28 | 2020-09-28 | 一种电力电子开关恒导纳模型建立方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011036925.3A CN112199914B (zh) | 2020-09-28 | 2020-09-28 | 一种电力电子开关恒导纳模型建立方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112199914A CN112199914A (zh) | 2021-01-08 |
CN112199914B true CN112199914B (zh) | 2024-06-04 |
Family
ID=74007502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011036925.3A Active CN112199914B (zh) | 2020-09-28 | 2020-09-28 | 一种电力电子开关恒导纳模型建立方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112199914B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375876A (zh) * | 2014-10-11 | 2015-02-25 | 清华大学 | 一种输入量突变情况下的0+误差免疫电磁暂态仿真算法 |
CN104462661A (zh) * | 2014-11-14 | 2015-03-25 | 上海交通大学 | 电力电子开关插值实时仿真方法 |
CN109543339A (zh) * | 2018-12-07 | 2019-03-29 | 上海交通大学 | 三电平换流器的定导纳建模与实时仿真方法 |
CN109614687A (zh) * | 2018-12-07 | 2019-04-12 | 上海交通大学 | 两电平桥式换流器的定导纳建模与实时仿真方法 |
CN110334476A (zh) * | 2019-07-19 | 2019-10-15 | 华北电力大学 | 一种电磁暂态仿真方法及系统 |
WO2020113874A1 (zh) * | 2018-12-07 | 2020-06-11 | 上海交通大学 | 电力电子换流器的定导纳建模与实时仿真方法 |
CN111709209A (zh) * | 2020-06-16 | 2020-09-25 | 华北电力大学 | 基于支路指数积分形式的电磁暂态仿真方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9665672B2 (en) * | 2015-04-27 | 2017-05-30 | Christian Dufour | Method and system for reducing power losses and state-overshoots in simulators for switched power electronic circuit |
-
2020
- 2020-09-28 CN CN202011036925.3A patent/CN112199914B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375876A (zh) * | 2014-10-11 | 2015-02-25 | 清华大学 | 一种输入量突变情况下的0+误差免疫电磁暂态仿真算法 |
CN104462661A (zh) * | 2014-11-14 | 2015-03-25 | 上海交通大学 | 电力电子开关插值实时仿真方法 |
CN109543339A (zh) * | 2018-12-07 | 2019-03-29 | 上海交通大学 | 三电平换流器的定导纳建模与实时仿真方法 |
CN109614687A (zh) * | 2018-12-07 | 2019-04-12 | 上海交通大学 | 两电平桥式换流器的定导纳建模与实时仿真方法 |
WO2020113874A1 (zh) * | 2018-12-07 | 2020-06-11 | 上海交通大学 | 电力电子换流器的定导纳建模与实时仿真方法 |
CN110334476A (zh) * | 2019-07-19 | 2019-10-15 | 华北电力大学 | 一种电磁暂态仿真方法及系统 |
CN111709209A (zh) * | 2020-06-16 | 2020-09-25 | 华北电力大学 | 基于支路指数积分形式的电磁暂态仿真方法及系统 |
Non-Patent Citations (5)
Title |
---|
一种基于指数积分的电力电子开关恒导纳模型;许明旺等;电网技术;20211018(第04期);1519-1529 * |
一种基于离散相似的电磁暂态仿真方法研究;宋文达等;中国电机工程学报;20200715;第40卷(第21期);6885-6894 * |
基于响应匹配的电力电子换流器恒导纳建模;徐晋等;中国电机工程学报;20190606;第39卷(第13期);3879-3889 * |
电力电子变压器的电磁暂态高效建模仿真方法研究;许明旺;中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑);20230315;第46卷(第03期);C042-508 * |
面向大规模新能源接入的电力系统暂态多时间尺度指数积分方法;富晓鹏;中国博士学位论文全文数据库(工程科技Ⅱ辑);20171215(第12期);C042-138 * |
Also Published As
Publication number | Publication date |
---|---|
CN112199914A (zh) | 2021-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pejovic et al. | A new algorithm for simulation of power electronic systems using piecewise-linear device models | |
CN106649927B (zh) | 一种基于fpga的电力电子元件实时仿真组合建模方法 | |
CN108229021B (zh) | 基于实时数字仿真器的模块化多电平换流器闭锁建模方法 | |
CN110518801B (zh) | 一种双重移相调制下双有源全桥变换器小信号建模方法 | |
CN110676841B (zh) | 基于直接法的电力电子化电力系统暂态稳定性分析方法 | |
CN109245532B (zh) | 一种升降压变换器的分数阶滑模控制方法 | |
CN108429286A (zh) | 一种基于自抗扰控制的并网电流调节器 | |
CN103746583A (zh) | 一种mmc的桥臂等效仿真方法 | |
CN103296905A (zh) | 三相电压型功率因数校正变换器的自适应控制方法 | |
Gao et al. | Fast simulation model of voltage source converters with arbitrary topology using switch-state prediction | |
Li et al. | Influence of non-ideal factors on the boundary control of buck converters with curved switching surfaces | |
CN112199914B (zh) | 一种电力电子开关恒导纳模型建立方法及系统 | |
CN106208791B (zh) | 一种九开关管三桥臂变流器模型搭建方法及其装置 | |
CN107546998A (zh) | 一种基于双环预测控制的切换型控制方法 | |
CN111898282B (zh) | 一种改进的模块化多电平换流器戴维南等效建模方法 | |
CN111709209B (zh) | 基于支路指数积分形式的电磁暂态仿真方法及系统 | |
CN112001142B (zh) | 一种半桥型模块化多电平换流器的实时仿真方法 | |
CN111628495B (zh) | 含vsc-mmc交直流电网系统的电磁暂态仿真方法 | |
CN112928942A (zh) | 一种三相逆变电路及其控制方法 | |
Pejovic et al. | An algorithm for solving piecewise-linear networks that include elements with discontinuous characteristics | |
CN112434411A (zh) | 采用变阶变步长3s-dirk算法的电磁暂态仿真方法 | |
Liu et al. | Simulation of switched-mode power conversion circuits with extended impedance method | |
Rajagopalan et al. | User-friendly dedicated power electronic converter simulator | |
Tse et al. | Stepwise quadratic state-space modeling technique for simulation of power electronics circuits | |
Wang et al. | Reduced-order model of VSC-HVDC applied to passive network based on optimal reduction algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |