CN105631086A - Damping circuit optimization design method for inhibiting numerical oscillation in simulation - Google Patents

Damping circuit optimization design method for inhibiting numerical oscillation in simulation Download PDF

Info

Publication number
CN105631086A
CN105631086A CN201510665187.1A CN201510665187A CN105631086A CN 105631086 A CN105631086 A CN 105631086A CN 201510665187 A CN201510665187 A CN 201510665187A CN 105631086 A CN105631086 A CN 105631086A
Authority
CN
China
Prior art keywords
circuit
simulation
oscillation
damping circuit
numerical
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.)
Pending
Application number
CN201510665187.1A
Other languages
Chinese (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201510665187.1A priority Critical patent/CN105631086A/en
Publication of CN105631086A publication Critical patent/CN105631086A/en
Pending legal-status Critical Current

Links

Landscapes

  • Design And Manufacture Of Integrated Circuits (AREA)

Abstract

本发明设计一种用于仿真中抑制数值振荡的阻尼电路优化设计方法,属电磁暂态领域。针对的是无插值功能的电磁暂态仿真平台中,由开关器件动作使得电感或电容等储能元件的非状态变量产生数值振荡的问题。本发明的核心技术方案是:通过本发明设计的阻尼电路参数优化方法,针对不同的仿真电路,依据具体的电路参数和拓扑结构设计不同的阻尼电路配置方案,基本流程是:1.先确定阻尼电路初值的选取;2.通过自定义模块找寻表征振荡程度的特征信息;3.采用simplex模块设计迭代环,通过多次的非线性迭代不断地更新R、C值,观察特征信息来确定抑制效果最优时的R、C值,达到阻尼电路的参数优化目的。The invention designs an optimal design method of a damping circuit used for suppressing numerical oscillation in simulation, and belongs to the field of electromagnetic transient state. It is aimed at the problem that in the electromagnetic transient simulation platform without interpolation function, the non-state variables of energy storage elements such as inductors or capacitors produce numerical oscillations due to the action of switching devices. The core technical solution of the present invention is: through the damping circuit parameter optimization method designed in the present invention, for different simulation circuits, design different damping circuit configuration schemes according to specific circuit parameters and topological structures. The basic process is: 1. First determine the damping Selection of the initial value of the circuit; 2. Find the characteristic information representing the degree of oscillation through the custom module; 3. Use the simplex module to design the iterative loop, continuously update the R and C values through multiple nonlinear iterations, and observe the characteristic information to determine the suppression The R and C values when the effect is optimal can achieve the purpose of parameter optimization of the damping circuit.

Description

一种用于仿真中抑制数值振荡的阻尼电路优化设计方法An Optimal Design Method for Damping Circuits Used to Suppress Numerical Oscillations in Simulation

技术领域 technical field

本发明属于电磁暂态仿真领域,尤其用于解决仿真中出现的数值振荡问题。 The invention belongs to the field of electromagnetic transient simulation, and is especially used for solving the numerical oscillation problem in simulation.

背景技术 Background technique

目前,电力系统电磁暂态仿真计算软件多种多样,应用的场合也不尽相同。针对不同类型的应用需求,大体上可以分为离线仿真平台和实时仿真器。离线仿真平台包括各种常见的软件包,如PSCAD/EMTDC,EMTP-RV,ATP,MicroTran等。实时仿真器的代表主要有RTDS,RT-LAB等,被广泛应用于工业和学术界。尽管应用场所不同,但不论是离线仿真还是实时仿真所采用的积分方法大都是隐式梯形积分法。隐式梯形积分法具有精度高、稳定性好等优点,但由于其积分特点的影响,在仿真过程中,开关和器件动作,网络结构发生变化,会引起电感电容等储能元件的非状态变量在事件发生后在真解附近不正常地摆动,也即电磁暂态仿真中的数值振荡现象。 At present, there are many kinds of electromagnetic transient simulation calculation software for power system, and the application occasions are also different. According to different types of application requirements, it can be roughly divided into offline simulation platform and real-time simulator. The off-line simulation platform includes various common software packages, such as PSCAD/EMTDC, EMTP-RV, ATP, MicroTran, etc. Representatives of real-time simulators mainly include RTDS, RT-LAB, etc., which are widely used in industry and academia. Although the application places are different, most of the integration methods used in offline simulation or real-time simulation are implicit trapezoidal integration methods. The implicit trapezoidal integration method has the advantages of high precision and good stability. However, due to the influence of its integral characteristics, during the simulation process, the switches and devices operate, and the network structure changes, which will cause non-state variables of energy storage components such as inductors and capacitors. After the event, it swings abnormally around the true solution, that is, the numerical oscillation phenomenon in the electromagnetic transient simulation.

发明内容 Contents of the invention

本发明提出了一种适用于无插值功能的仿真平台中抑制数值振荡的RC阻尼电路设计方法。首先从分析数值振荡的机理入手,确定了阻尼电路的初值配置,并提出了描述数值振荡剧烈程度的量化函数SUM,并将其与simplex算法相结合对RC电路的配置进行优化,大大提高了阻尼电路对数值振荡的抑制效果。 The invention proposes an RC damping circuit design method suitable for suppressing numerical oscillation in a simulation platform without an interpolation function. First, starting from the analysis of the mechanism of numerical oscillation, the initial value configuration of the damping circuit is determined, and a quantization function SUM that describes the severity of numerical oscillation is proposed, and the configuration of the RC circuit is optimized by combining it with the simplex algorithm, which greatly improves the performance of the damping circuit. The damping effect of the damping circuit on numerical oscillations.

本发明的技术方案特征包括以下步骤: Technical solution feature of the present invention comprises the following steps:

步骤1:分析数值振荡产生的机理,针对于不同的电路参数,结合仿真平台所使用的积分方法,考虑电路参数同仿真时间步长的相互配合,确定RC电路的初值配置。 Step 1: Analyze the mechanism of numerical oscillation. For different circuit parameters, combine the integration method used by the simulation platform, and consider the interaction between circuit parameters and simulation time steps to determine the initial value configuration of the RC circuit.

步骤2:根据数值振荡产生时的电气量振荡情况,使用Fortran语言编写自定义模块,生成用于判断振荡程度的特征信号。 Step 2: According to the electrical quantity oscillation when the numerical oscillation is generated, use the Fortran language to write a custom module to generate a characteristic signal for judging the degree of oscillation.

步骤3:将步骤2产生的特征信号作为输入,同simplex模块相结合,设计迭代环,通过非线性迭代得到RC电路的最优配置。 Step 3: Take the characteristic signal generated in step 2 as an input, combine it with the simplex module, design an iterative loop, and obtain the optimal configuration of the RC circuit through nonlinear iteration.

本发明通过上述三个步骤,能够对RC阻尼电路的初值进行确定,并能够通过simplex算法对其参数进一步优化,并找到最优解,有效的解决了仿真中出现的数值振荡问题。 Through the above three steps, the present invention can determine the initial value of the RC damping circuit, further optimize its parameters through the simplex algorithm, and find the optimal solution, effectively solving the numerical oscillation problem in the simulation.

附图说明 Description of drawings

图1为加入RC阻尼电路的RL测试电路。图中E s RL是RL测试电路的具体参数配置,D为二极管,R snb C snb 是阻尼电路的具体参数配置。 Figure 1 is the RL test circuit with RC damping circuit added. In the figure, E s , R , L are the specific parameter configurations of the RL test circuit, D is the diode, R snb , C snb are the specific parameter configurations of the damping circuit.

图2为阻尼电路的戴维南变形电路,R SNB为阻尼电路的等效电阻,V SNB_EQ为阻尼电路的等效历史电压源,其值由上一时刻的电容电压、电容电流、阻尼电阻和仿真步长决定。 Figure 2 is the Thevenin deformation circuit of the damping circuit, RSNB is the equivalent resistance of the damping circuit, V SNB_EQ is the equivalent historical voltage source of the damping circuit, and its value is determined by the capacitor voltage, capacitor current, damping resistance and simulation step at the previous moment long decision.

图3为特征信号SUM函数的产生流程图。图中V L是图1中电感L两端的电压,I L是图1中流过电感L的电流。 Fig. 3 is a flow chart of the generation of the characteristic signal SUM function. In the figure, V L is the voltage across the inductor L in Figure 1, and I L is the current flowing through the inductor L in Figure 1.

图4为RC阻尼电路的优化流程图,其中量化函数SUM即为图2中所述的特征信号,值得说明的是,随着量化函数SUM的不断更新,阻尼电路R和C的取值也会随之更新,从而得到RC电路的最优配置。 Figure 4 is an optimization flow chart of the RC damping circuit, where the quantization function SUM is the characteristic signal described in Figure 2. It is worth noting that, with the continuous update of the quantization function SUM , the values of R and C of the damping circuit will also be Then it is updated to obtain the optimal configuration of the RC circuit.

具体实施方式 detailed description

下面将对本发明涉及的一种基于优化算法simplex的阻尼电路设计方法作详细说明。应该强调的是,下述说明仅仅是示例性的,而不是为了限制本发明的范围及其应用。 A damping circuit design method based on the optimization algorithm simplex involved in the present invention will be described in detail below. It should be emphasized that the following description is only exemplary and not intended to limit the scope of the invention and its application.

本发明所要解决的技术问题是通过优化阻尼电路的参数,对仿真中出现的数值振荡问题进行最大程度的抑制,在不借助附加仿真算法的情况下,使得仿真结果更加准确可信。本发明采用如下技术方案实现: The technical problem to be solved by the present invention is to suppress the numerical oscillation problem in the simulation to the greatest extent by optimizing the parameters of the damping circuit, and to make the simulation result more accurate and credible without resorting to additional simulation algorithms. The present invention adopts following technical scheme to realize:

本发明通过如下三步来实现: The present invention realizes through following three steps:

步骤1:如图1所示,将阻尼电路并联在二极管两端,接下来确定阻尼电路R和C的初值。 Step 1: As shown in Figure 1, connect the damping circuit in parallel to both ends of the diode, and then determine the initial values of R and C of the damping circuit.

R的初值确定:基于二值电阻的理论,二极管导通电阻很小,一般取R on=0.01Ω;二极管关断电阻很大,一般取R off=1MΩ。显然,不论二极管处于通断哪种状态,RC阻尼电路都应尽量减少其对二极管造成的影响,也即,R snb不应过小以免在二极管导通时分流或者是二极管关断时通流,R snb也不应过大,否则将和R off的作用类似以至于削弱暂态阻尼功能的效果。基于以上分析,这里将R snb的初值粗略地取为5kΩ,恰处于0.01Ω和10MΩ的数量级之间。 Determination of the initial value of R: Based on the theory of binary resistance, the on-resistance of the diode is very small, generally R on =0.01Ω; the off-resistance of the diode is large, generally R off =1MΩ. Obviously, regardless of the on-off state of the diode, the RC damping circuit should minimize its impact on the diode, that is, R snb should not be too small to avoid shunting when the diode is on or flowing when the diode is off. R snb should not be too large, otherwise it will be similar to the effect of R off so as to weaken the effect of the transient damping function. Based on the above analysis, the initial value of R snb is roughly taken as 5kΩ, which is just in the order of magnitude between 0.01Ω and 10MΩ.

C的初值确定:由图2可知在关断瞬间,二极管两端的电压也可等效为一个等效电阻R SNB和一个等效历史电压源V SNB_EQ的串联。可以发现,为了使电流i(t+Δt)不真正为零,R SNB应该和R snb取同样数量级,不应太大,又因为这里讨论的仿真情形主要是步长大于10μs的离线仿真,故这里C的初值粗略地选取为0.02μF。可以推断,随着后期优化工作的进行,为了体现RC阻尼电路的效果,C snb值的选取将会向上收敛,因为此时C snb=0.02μF的初值是一个近似于开路的一个保守的配置,而为了抑制电感电压不正常的波动,电容C snb的值就应该足够大来抑制其波动。当然,C snb的取值也不应太大,因为一旦大过一定范围,阻尼电路就相当于对高频电流呈现一个短路的状态,即使后期加入优化算法,也有可能无法收敛到其准确值。 Determination of the initial value of C : It can be seen from Figure 2 that at the moment of turning off, the voltage across the diode can also be equivalent to a series connection of an equivalent resistance R SNB and an equivalent historical voltage source V SNB_EQ . It can be found that, in order to make the current i (t+Δt) not really zero, R SNB should be of the same order of magnitude as R snb , and should not be too large, and because the simulation situation discussed here is mainly an off-line simulation with a step size larger than 10μs, so Here the initial value of C is roughly selected as 0.02μF. It can be inferred that with the optimization work in the later stage, in order to reflect the effect of the RC damping circuit, the selection of the value of C snb will converge upwards, because the initial value of C snb =0.02μF is a conservative configuration close to an open circuit , and in order to suppress the abnormal fluctuation of the inductor voltage, the value of the capacitor C snb should be large enough to suppress its fluctuation. Of course, the value of C snb should not be too large, because once it is larger than a certain range, the damping circuit is equivalent to a short circuit state for high-frequency currents. Even if the optimization algorithm is added later, it may not be able to converge to its exact value.

步骤2:使用Fortran语言编写自定义模块SUM,这里的SUM函数值即表征着数值振荡的剧烈程度。如图3所示,将电感电压和电感电流作为原始信息采集进来,进一步加工,输出的量化函数SUM值即作为特征信息,是下一步迭代更新RC值的输入信息。函数SUM具体算法如式(1)所示。 Step 2: Use the Fortran language to write the custom module SUM , where the SUM function value represents the intensity of numerical oscillation. As shown in Figure 3, the inductor voltage and inductor current are collected as raw information and further processed. The output quantization function SUM value is used as the characteristic information, which is the input information for the next iterative update of the RC value. The specific algorithm of the function SUM is shown in formula (1).

(1) (1)

步骤3:将步骤2得到的量化函数SUM值作为输入量,采用simplex优化模块,按如图4所示的流程设计程序,编写迭代环。通过不停的迭代,根据SUM值的变化找到最合适的RC值,使得抑制振荡的效果达到最佳。 Step 3: Use the quantization function SUM value obtained in step 2 as the input quantity, use the simplex optimization module, design the program according to the flow shown in Figure 4, and write the iterative loop. Through continuous iteration, the most suitable RC value is found according to the change of the SUM value, so that the effect of suppressing the oscillation is optimal.

以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。 The above is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any person skilled in the art within the technical scope disclosed in the present invention can easily think of changes or Replacement should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be determined by the protection scope of the claims.

Claims (2)

1. one kind for suppressing the antihunt circuit Optimization Design of numerical oscillation in emulating, it is characterized in that without in the emulation platform of interpolation function, for the numerical oscillation problem that switching device action causes, on the basis of existing RC antihunt circuit, a kind of brand-new method for optimally designing parameters is proposed, make RC antihunt circuit parameter value possess certain theoretical basis, also make the effect suppressing vibration reach optimum, comprise the following steps:
Step 1: by circuit is made a concrete analysis of, it is considered to the order of magnitude relation between simulation step length, integration method and circuit parameter, completes the initial value configuration of RC circuit;
Step 2: the phenomenon characteristic produced according to numerical oscillation, uses Fortran language to write suitable custom block, gathers current electric parameters and generates a characteristic signal that can characterize degree of oscillation;
Step 3: characteristic signal step 2 generated, as input quantity, optimizes the input quantity of module, and finds R, C value one group optimum by nonlinear iteration repeatedly as simplex so that suppress the effect of vibration to reach optimum.
2. a kind of for suppressing the antihunt circuit Optimization Design of numerical oscillation in emulating based on described in claim 1, it is characterized in that step 1,2 and 3 entirety are as summary of the invention, the design making antihunt circuit has had certain theoretical foundation, can significantly more efficient solving without numerical oscillation problem in interpolation emulation platform, three steps are organic indivisible entirety.
CN201510665187.1A 2015-10-16 2015-10-16 Damping circuit optimization design method for inhibiting numerical oscillation in simulation Pending CN105631086A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510665187.1A CN105631086A (en) 2015-10-16 2015-10-16 Damping circuit optimization design method for inhibiting numerical oscillation in simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510665187.1A CN105631086A (en) 2015-10-16 2015-10-16 Damping circuit optimization design method for inhibiting numerical oscillation in simulation

Publications (1)

Publication Number Publication Date
CN105631086A true CN105631086A (en) 2016-06-01

Family

ID=56046017

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510665187.1A Pending CN105631086A (en) 2015-10-16 2015-10-16 Damping circuit optimization design method for inhibiting numerical oscillation in simulation

Country Status (1)

Country Link
CN (1) CN105631086A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009169733A (en) * 2008-01-17 2009-07-30 Hitachi Ltd Fluid analysis apparatus and method
CN103401416A (en) * 2013-07-31 2013-11-20 西安交通大学 Main circuit structure eliminating right half plane zeros of boost DC-DC (Direct Current - Direct Current) converter and method for determining parameters thereof
JP5500637B2 (en) * 2010-03-09 2014-05-21 国立大学法人東京工業大学 Numerical calculation method, program, and recording medium
CN103986138A (en) * 2014-05-14 2014-08-13 国家电网公司 A Modular Current Limiting Circuit Breaker Power Module

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009169733A (en) * 2008-01-17 2009-07-30 Hitachi Ltd Fluid analysis apparatus and method
JP5500637B2 (en) * 2010-03-09 2014-05-21 国立大学法人東京工業大学 Numerical calculation method, program, and recording medium
CN103401416A (en) * 2013-07-31 2013-11-20 西安交通大学 Main circuit structure eliminating right half plane zeros of boost DC-DC (Direct Current - Direct Current) converter and method for determining parameters thereof
CN103986138A (en) * 2014-05-14 2014-08-13 国家电网公司 A Modular Current Limiting Circuit Breaker Power Module

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
商莹等: "非线性电路暂态仿真中消除数值振荡的改进方法", 《电力系统保护与控制》 *
王毅等: "消除电力系统故障仿真中数值振荡的方法研究", 《华东交通大学学报》 *
许建中等: "超大规模MMC 电磁暂态仿真提速模型", 《中国电机工程学报》 *

Similar Documents

Publication Publication Date Title
Wang et al. A generalized associated discrete circuit model of power converters in real-time simulation
CN105260516B (en) A kind of electromagnetical transient emulation method of the sub-network containing switching characteristic
Chen et al. Model-Free Predictive $ H_ {\infty} $ Control for Grid-Connected Solar Power Generation Systems
CN108988361B (en) A Fast Suppression Method for Chaotic Oscillation in Dual-machine Interconnected Power System
WO2022213479A1 (en) Large-step frequency-shift electro-magnetic transient simulation method and system
CN113158617B (en) General decoupling method and system for electromagnetic transient simulation of voltage source type converter
Ebrahimzadeh et al. Reducing harmonic instability and resonance problems in PMSG-based wind farms
Mahdavi et al. Dynamic transmission network expansion planning considering network losses DG sources and operational costs-part 1: Review and problem formulation
Chalangar et al. A direct mapped method for accurate modeling and real-time simulation of high switching frequency resonant converters
Zheng et al. A semi-implicit parallel leapfrog solver with half-step sampling technique for FPGA-based real-time HIL simulation of power converters
CN108959671A (en) The real-time simulation modeling method of half-bridge and bridge-type modularization multi-level converter
CN103530440A (en) Micro-grid harmonic suppression method based on particle swarm optimization algorithm
Maguire et al. Predicting switch ON/OFF statuses in real time electromagnetic transients simulations with voltage source converters
Khan et al. Singular perturbation‐based model reduction of power electronic circuits
CN103258131A (en) Power circuit component optimization method based on orthogonal learning particle swarm
CN117034854A (en) Transformer electric field time domain simulation method and system
Tavakoli et al. A self‐constructing Lyapunov neural network controller to track global maximum power point in PV systems
Zhang et al. Machine learning based modeling for real-time inferencer-in-the-loop hardware emulation of high-speed rail microgrid
CN106099939A (en) A kind of transformer station reactive apparatus affects the computational methods of sensitivity to busbar voltage
Zhang et al. General Linearized Model of Voltage Source Converter with Fixed Nodal Admittance Matrix
CN105631086A (en) Damping circuit optimization design method for inhibiting numerical oscillation in simulation
CN104467741A (en) Intelligent current tracking control method used for active filter and based on T-S fuzzy modeling
CN112117906B (en) Double-active full-bridge converter optimization method under triple phase shift control
Berger et al. Large-signal modeling of three-phase dual active bridge converters for electromagnetic transient analysis in DC grids
Khan et al. Simulation acceleration of high-fidelity nonlinear power electronic circuits using model order reduction

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20160601