CN108920791B - Power electronic switch modeling method and model - Google Patents

Power electronic switch modeling method and model Download PDF

Info

Publication number
CN108920791B
CN108920791B CN201810644176.9A CN201810644176A CN108920791B CN 108920791 B CN108920791 B CN 108920791B CN 201810644176 A CN201810644176 A CN 201810644176A CN 108920791 B CN108920791 B CN 108920791B
Authority
CN
China
Prior art keywords
equivalent
switch
branch
resistance
admittance
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
Application number
CN201810644176.9A
Other languages
Chinese (zh)
Other versions
CN108920791A (en
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.)
CSG Electric Power Research Institute
China Southern Power Grid Co Ltd
Original Assignee
CSG Electric Power Research Institute
China Southern Power Grid 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 CSG Electric Power Research Institute, China Southern Power Grid Co Ltd filed Critical CSG Electric Power Research Institute
Priority to CN201810644176.9A priority Critical patent/CN108920791B/en
Publication of CN108920791A publication Critical patent/CN108920791A/en
Application granted granted Critical
Publication of CN108920791B publication Critical patent/CN108920791B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Abstract

The invention discloses a power electronic switch modeling method which is characterized by comprising the following steps: the conducted switch is equivalent to an inductive resistor series branch circuit; the switched-off switch is equivalent to a capacitor-resistor series branch; a negative resistance device is connected in series in the switching branch. The power electronic switch modeling method disclosed by the invention can enable the admittance matrix in the on and off states after the switch circuit is equivalent to be unchanged, the equivalent inductance parameter and the capacitance parameter are easy to determine, and the power balance is maintained. The embodiment of the invention also discloses a power electronic switch model.

Description

Power electronic switch modeling method and model
Technical Field
The invention relates to the technical field of power electronics, in particular to a power electronic switch modeling method and a model.
Background
In recent years, with the development of power semiconductor devices and automatic control technologies, power electronic systems have been widely used. In the research and development process of the power electronic system, the digital simulation technology plays an important role, and can fully verify the circuit design and the control algorithm of the power electronic system, so that the risk of trial production of products is reduced.
From the perspective of system analysis, circuit simulation involves solving a set of first order differential equations that satisfy kirchhoff's law, i.e., describing the dynamic behavior of a circuit as it changes over time after being subjected to a particular stimulus. Meanwhile, in order to solve the problem by using a digital computer, the system of differential equations needs to be discretized properly so as to obtain the differential equation description of the circuit, so that the iterative calculation of the digital computer is facilitated. By combining the circuit equation with the discretization method, a generalized algorithm, namely an electromagnetic transient Program (EMTP), is generated, and accurate transient simulation calculation can be performed on the circuit containing the distribution parameters and the concentration parameters. In the EMTP method, the differential equation of each circuit element can be regarded as a "norton equivalent circuit" composed of a history current source and an equivalent admittance. And for the whole circuit, obtaining a system admittance matrix according to a node admittance method, thereby forming a differential equation of the whole circuit for iterative computation.
However, due to the complexity of modern power systems, typically only simple circuits (e.g., RLC circuits) have an intuitive differential equation description; complex nonlinear circuits related to power electronic switches and the like generally need to adopt a proper equivalent method to obtain an equivalent model described by simple circuit elements. The power electronic switch is the key of electric energy conversion, so that the research of a simulation model of a power electronic switch device becomes one of important works of electromagnetic transient simulation of a power electronic system.
Generally, there are two power electronic switch equivalent modeling approaches:
(1) the large and small resistance method, i.e. when the switch is on, it is equivalent to a small resistance, and when the switch is off, it is equivalent to a large resistance.
(2) The inductor-capacitor method, i.e. when the switch is on, it is equivalent to a small inductor, and when the switch is off, it is equivalent to a small capacitor. The FPGA (Field-Programmable Gate Array) programming is realized by selecting proper inductance and capacitance parameters.
The inventor finds that the following technical problems exist in the prior art in the process of implementing the invention:
when the large and small resistance method is adopted, the switching action leads to the regeneration of a system admittance matrix and a difference equation set of the circuit, and if necessary, a plurality of iterative operations are required to be carried out at a switching point so as to prevent numerical value oscillation. Along with the increase of the number and frequency of power electronic switches, the time consumption of electromagnetic transient calculation is also increased sharply, and meanwhile, a time-varying system matrix also causes difficulty in Field-Programmable Gate Array (FPGA) programming, so that the calculation efficiency is difficult to be improved by applying the latest parallelization calculation technology; when an inductance-capacitance method is adopted, the selection of inductance and capacitance parameters is very critical and needs to be matched with a simulation step length and circuit port equivalent impedance, but when the method is applied specifically, the matching process is complex due to the lack of a quantification guiding principle, and numerical oscillation is easily generated by the equivalent method of pure inductance and capacitance branches; the absence of compensation elements keeps the circuit power balanced.
Disclosure of Invention
The embodiment of the invention provides a power electronic switch model and a model, which can effectively overcome the defects that in the prior art, an admittance matrix of a system is time-varying, equivalent parameters of inductance and capacitance are difficult to determine, and a compensation element is lacked.
An embodiment of the present invention provides a power electronic switch modeling method, including:
the conducted switch is equivalent to an inductive resistor series branch circuit;
the switched-off switch is equivalent to a capacitor-resistor series branch;
a negative resistance device is connected in series in the switching branch.
As an improvement of the above scheme, after a conducting switch is equivalent to an inductor-resistor series branch, a branch transfer function obtained by laplace transform is as the following formula (1):
Figure BDA0001703120610000021
wherein s is a complex variable, L is an equivalent inductance, and R is an equivalent resistance.
As an improvement of the above scheme, after the turned-on switch is equivalent to an inductor-resistor series branch, the equivalent admittance and the history current of the inductor-resistor series branch obtained according to the electromagnetic transient program are as the following formula (2):
Figure BDA0001703120610000031
wherein, YrlIs equivalent admittance, It-ΔtFor the branch current of the last calculation time step, Ihist.rlFor updated historical current, Δ t is the simulation calculation step length, e is the natural constant, L is the equivalent inductance, and R is the equivalent resistance.
As an improvement of the above scheme, after the turned-off switch is equivalent to a capacitor-resistor series branch, the branch transfer function obtained by laplace transform is as the following formula (3):
Figure BDA0001703120610000032
wherein s is a complex variable, C is an equivalent capacitor, and R is an equivalent resistor.
As an improvement of the above scheme, after the turned-off switch is equivalent to a capacitance-resistance series branch, the equivalent admittance and the history current of the capacitance-resistance series branch obtained according to the electromagnetic transient program are as the following formula (4):
Figure BDA0001703120610000033
wherein, YrcIs equivalent admittance, It-ΔtFor the branch current of the last calculation time step, Ut-ΔtFor the branch voltage of the last calculation time step, Ihist.rcFor updated historical current, Δ t is the simulation calculation step length, e is the natural constant, C is the equivalent capacitance, and R is the equivalent resistance.
As an improvement of the scheme, the admittance matrix in the equivalent circuit is unchanged when the switch is turned on and turned off.
As an improvement of the scheme, the equivalent admittance Y after the equivalent of the conducted switchrlAs in the following formula (5):
Figure BDA0001703120610000034
equivalent admittance Y after switch-off equivalencercIs as the following formula (6);
Figure BDA0001703120610000035
from Yrl=YrcObtaining the formula (7):
Figure BDA0001703120610000041
the equivalent inductance and equivalent capacitance parameters are as in the following equation (8):
Figure BDA0001703120610000042
wherein, YrlIs the equivalent admittance of the switch when it is on, YrcAnd the equivalent admittance when the switch is turned off, delta t is the step length of simulation calculation, e is a natural constant, C is the equivalent capacitance, L is the equivalent inductance, and R is the equivalent resistance.
As an improvement of the above scheme, an absolute value of the resistance value of the negative resistance device is equal to the resistance value of the equivalent resistor.
Another embodiment of the present invention correspondingly provides a power electronic switch model, which includes: capacitance resistance branch circuit module, inductance resistance branch circuit module, first negative resistance device Rcomp1And a second negative resistance device Rcomp2(ii) a The capacitance-resistance branch circuit module comprises a first capacitor C1 and a first resistor R1; the inductance-resistance branch module comprises a first inductance L1 and a second resistance R2;
the first capacitor C1 and the first resistor R1 are connected in series at two ends of a power supply;
the first inductor L1 and the second resistor R2 are connected in series at two ends of a power supply;
the first negative resistance device Rcomp1The capacitor resistor branch circuit module is connected between a power supply and the capacitor resistor branch circuit module in series;
the second negative resistance device Rcomp2The inductor resistor branch circuit module is connected between a power supply and the inductor resistor branch circuit module in series.
The modeling method for the power electronic switch provided by the embodiment of the invention has the following beneficial effects:
the switches in the on and off states are respectively equivalent to an inductance resistance series branch and a capacitance resistance series branch, and the expressions of inductance parameters and capacitance parameters of the equivalent branches are calculated when the admittance matrixes of the two branches are the same, so that the system admittance matrixes before and after the switch states are changed are unchanged, and the method is suitable for Field Programmable Gate Array (FPGA) programming; the expression of the equivalent inductance parameter and the equivalent capacitance parameter is only related to the equivalent resistance and the simulation calculation step length, and the equivalent resistance and the simulation calculation step length are easy to obtain, so the equivalent inductance parameter and the equivalent capacitance parameter are easy to determine; a negative resistance device is connected in series with the switching branch circuit, so that the power balance of the circuit can be kept.
Drawings
Fig. 1 is a schematic flowchart of a power electronic switch modeling method according to an embodiment of the present invention.
Fig. 2 is a schematic structural diagram of a power electronic switch model according to a second embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, a schematic flow chart of a power electronic switch modeling method according to an embodiment of the present invention includes:
s1, enabling the conducted switch to be equivalent to an inductive resistor series branch circuit;
the branch transfer function obtained by laplace transform is shown as the following formula (1):
Figure BDA0001703120610000051
the equivalent admittance and the historical current of the inductance-resistance series branch obtained according to the electromagnetic transient program are as the following formula (2):
Figure BDA0001703120610000052
wherein s is a complex variable, L is an equivalent inductance, R is an equivalent resistance, and Y isrlFor equivalent admittance, It-ΔtFor the branch current of the last calculation time step, Ihist.rlFor updated historical current, Δ t is the simulation calculation step length, and e is a natural constant.
S2, enabling the switched-off switch to be equivalent to a capacitor-resistor series branch;
the branch transfer function obtained by laplace transform is shown as the following formula (3):
Figure BDA0001703120610000053
the equivalent admittance and the historical current of the capacitance-resistance series branch obtained according to the electromagnetic transient program are as the following formula (4):
Figure BDA0001703120610000061
wherein s is a complex variable, C is an equivalent capacitance, R is an equivalent resistance, and Y isrcIs equivalent admittance, It-ΔtFor the branch current of the last calculation time step, Ut-ΔtFor the branch voltage of the last calculation time step, Ihist.rcFor updated historical current, Δ t is the simulation calculation step length, and e is a natural constant.
For the system admittance matrix to be constant when the switch is switched off and on, Y is providedrl=Yrc
Equivalent admittance Y after equivalent of switch onrlAs in the following formula (5):
Figure BDA0001703120610000062
equivalent admittance Y after switch-off equivalencercIs as the following formula (6);
Figure BDA0001703120610000063
from Yrl=YrcObtaining the formula (7):
Figure BDA0001703120610000064
the equivalent inductance and equivalent capacitance parameters are as in the following equation (8):
Figure BDA0001703120610000065
wherein, YrlIs the equivalent admittance of the switch when it is on, YrcAnd the equivalent admittance when the switch is turned off, delta t is the step length of simulation calculation, e is a natural constant, C is the equivalent capacitance, L is the equivalent inductance, and R is the equivalent resistance.
S3, connecting a negative resistance device in series in the switch branch circuit;
wherein, the absolute value of the resistance of the negative resistance device is equal to the resistance of the equivalent resistor.
Referring to fig. 2, a schematic structural diagram of a power electronic switch model provided in the second embodiment of the present invention includes: capacitance resistance branch module, inductance resistance branch module, first negative resistance device Rcomp1And a second negative resistance device Rcomp2(ii) a The capacitance-resistance branch circuit module comprises a first capacitor C1 and a first resistor R1; the inductance-resistance branch module comprises a first inductance L1 and a second resistance R2;
the first capacitor C1 and the first resistor R1 are connected in series at two ends of the power supply;
the first inductor L1 and the second resistor R2 are connected in series at two ends of the power supply;
first negative resistance device Rcomp1The capacitor resistor branch circuit module is connected between the power supply and the capacitor resistor branch circuit module in series;
second negative resistance device Rcomp2And the inductor resistor branch circuit module is connected between the power supply and the inductor resistor branch circuit module in series.
The modeling method for the power electronic switch provided by the embodiment of the invention has the following beneficial effects:
the switches in the on and off states are respectively equivalent to an inductance resistance series branch and a capacitance resistance series branch, and the expressions of inductance parameters and capacitance parameters of the equivalent branches are calculated when the admittance matrixes of the two branches are the same, so that the system admittance matrixes before and after the switch states are changed are unchanged, and the method is suitable for Field Programmable Gate Array (FPGA) programming; the expression of the equivalent inductance parameter and the equivalent capacitance parameter is only related to the equivalent resistance and the simulation calculation step length, and the equivalent resistance and the simulation calculation step length are easy to obtain, so that the equivalent inductance parameter and the equivalent capacitance parameter are easy to determine; a negative resistance device is connected in series with the switching branch circuit, so that the power balance of the circuit can be kept.
It should be noted that the above-described model embodiments are only illustrative, where the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.

Claims (2)

1. A power electronic switch modeling method, comprising:
the conducted switch is equivalent to an inductive resistor series branch circuit;
the switched-off switch is equivalent to a capacitor-resistor series branch;
a negative resistance device is connected in series in the switch branch circuit;
after the switched-on switch is equivalent to an inductive resistor series branch, a branch transfer function obtained through Laplace transformation is as the following formula (1):
Figure FDA0003513830770000011
wherein s is a complex variable, L is an equivalent inductance, and R is an equivalent resistance;
after the conducted switch is equivalent to an inductive resistance series branch, the equivalent admittance and the historical current of the inductive resistance series branch obtained according to the electromagnetic transient program are as the following formula (2):
Figure FDA0003513830770000012
wherein, YrlIs equivalent admittance, It-ΔtFor the branch current of the last calculation time step, Ihist.rlFor the updated historical current, delta t is a simulation calculation step length, e is a natural constant, L is an equivalent inductance, and R is an equivalent resistance;
after the switched-off switch is equivalent to a capacitor-resistor series branch, a branch transfer function obtained through Laplace transformation is as the following formula (3):
Figure FDA0003513830770000013
wherein s is a complex variable, C is an equivalent capacitor, and R is an equivalent resistor;
after the switched-off switch is equivalent to a capacitance-resistance series branch, the equivalent admittance and the historical current of the capacitance-resistance series branch obtained according to the electromagnetic transient program are as the following formula (4):
Figure FDA0003513830770000014
wherein, YrcIs equivalent admittance, It-ΔtFor the branch current of the last calculation time step, Ut-ΔtFor the branch voltage of the last calculation time step, Ihist.rcFor the updated historical current, delta t is a simulation calculation step length, e is a natural constant, C is an equivalent capacitor, and R is an equivalent resistor;
the admittance matrix in the equivalent circuit is unchanged when the switch is switched on and switched off;
equivalent admittance Y after equivalent of switch onrlAs in the following formula (5):
Figure FDA0003513830770000021
equivalent admittance Y after switch-off equivalencercIs as the following formula (6);
Figure FDA0003513830770000022
from Yrl=YrcObtaining the formula (7):
Figure FDA0003513830770000023
the equivalent inductance and equivalent capacitance parameters are as in the following equation (8):
Figure FDA0003513830770000024
wherein Y isrlIs the equivalent admittance of the switch when it is on, YrcThe equivalent admittance when the switch is turned off is shown, delta t is the step length of simulation calculation, e is a natural constant, C is the equivalent capacitance, L is the equivalent inductance, and R is the equivalent resistance;
the absolute value of the resistance value of the negative resistance device is equal to the resistance value of the equivalent resistor.
2. A power electronic switch model, characterized in that it is obtained using the modeling method of claim 1.
CN201810644176.9A 2018-06-21 2018-06-21 Power electronic switch modeling method and model Active CN108920791B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810644176.9A CN108920791B (en) 2018-06-21 2018-06-21 Power electronic switch modeling method and model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810644176.9A CN108920791B (en) 2018-06-21 2018-06-21 Power electronic switch modeling method and model

Publications (2)

Publication Number Publication Date
CN108920791A CN108920791A (en) 2018-11-30
CN108920791B true CN108920791B (en) 2022-05-06

Family

ID=64421421

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810644176.9A Active CN108920791B (en) 2018-06-21 2018-06-21 Power electronic switch modeling method and model

Country Status (1)

Country Link
CN (1) CN108920791B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109614687B (en) * 2018-12-07 2023-04-28 上海交通大学 Constant admittance modeling and real-time simulation method for two-level bridge type converter
CN111725818B (en) * 2019-03-18 2023-05-19 杨利 Grid-connected simulation method and simulation terminal for three-phase weak current network converter group
CN111817333B (en) * 2020-06-03 2022-08-12 中国电力科学研究院有限公司 Alternating current overvoltage calculation method and system based on flexible direct current power transmission system characteristics
CN114818364B (en) * 2022-05-17 2023-03-14 四川大学 Bipolar direct current cable frequency-dependent RLC modeling method based on head wave valley amplitude-frequency characteristic fitting

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793562A (en) * 2014-01-05 2014-05-14 天津大学 Active power distribution network transient state real-time simulation system designing method based on FPGA
CN105742124A (en) * 2016-05-03 2016-07-06 北京邮电大学 Microelectromechanical system switch
WO2017091118A1 (en) * 2015-11-27 2017-06-01 Telefonaktiebolaget Lm Ericsson (Publ) Determining the equivalent series resistance of a power converter
CN106877819A (en) * 2016-11-21 2017-06-20 章策珉 Voltage controlled oscillator based on compound resonator

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103793562A (en) * 2014-01-05 2014-05-14 天津大学 Active power distribution network transient state real-time simulation system designing method based on FPGA
WO2017091118A1 (en) * 2015-11-27 2017-06-01 Telefonaktiebolaget Lm Ericsson (Publ) Determining the equivalent series resistance of a power converter
CN105742124A (en) * 2016-05-03 2016-07-06 北京邮电大学 Microelectromechanical system switch
CN106877819A (en) * 2016-11-21 2017-06-20 章策珉 Voltage controlled oscillator based on compound resonator

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Synchronized switch damping on inductor and negative capacitance;Bilal Mokrani 等;《Journal of Intelligent Material Systems and Structures》;20121231;第18卷(第23期);第2065-2075页 *
一种机载DC/DC浪涌电流抑制电路的设计;杜培德 等;《微电子学》;20180420;第48卷(第2期);第216-221页 *

Also Published As

Publication number Publication date
CN108920791A (en) 2018-11-30

Similar Documents

Publication Publication Date Title
CN108920791B (en) Power electronic switch modeling method and model
Petráš Chaos in the fractional-order Volta’s system: modeling and simulation
Riaza Nondegeneracy conditions for active memristive circuits
CN104410485A (en) A wien-brigde chaotic oscillator based on a firs-order generalized memory resistor
CN101291062B (en) Analysis method of DC power flow in dynamic topology electric network
Petržela Optimal Piecewise-Linear Approximation of the Quadratic Chaotic Dynamics.
He et al. Solutions and memory effect of fractional-order chaotic system: A review
CN110414118A (en) A kind of Boost modeling method and application based on separate type modeling
AbdelAty et al. Optimal charging of fractional-order circuits with Cuckoo search
Mandache et al. Sensitivity analysis of analog circuits based on a modified nodal approach
CN104573164A (en) Hardware realization method for updating equivalent historical current source information of small-step system
Almakhles et al. Single-bit modulator based controller for capacitive power transfer system
CN112597730B (en) T-type three-level converter simulation method and device, electronic equipment and storage medium
WO2022213479A1 (en) Large-step frequency-shift electro-magnetic transient simulation method and system
CN109460575B (en) Hybrid simulation interface system based on Bergeron equivalent circuit and interpolation algorithm
Revankar Topological approach to thyristor-circuit analysis
CN107291981B (en) Simulation method and device of power transmission line
CN112182815A (en) Multi-rate real-time simulation method based on Nonton equivalence
Zhan et al. Taming the stability-constrained performance optimization challenge of distributed on-chip voltage regulation
CN112818537B (en) Photovoltaic grid-connected system stability analysis method and device
CN110858263A (en) Electrical circuit modeling method, simulation test system and simulation terminal
CN108134526A (en) A kind of distributed pulse duration modulation method to be cooperated based on multinuclear
CN104731989B (en) A kind of quick calculation method of the end-to-end equivalent resistance of extensive resistor network
Malyna et al. A comparison of methods for finding steady-state solution in power electronic circuits
Chen et al. Sensitivity analysis for time domain response of transmission lines based on the precise integration method

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