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 PDFInfo
- 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
- damping circuit
- oscillation
- simulation
- numerical oscillation
- 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
Links
Landscapes
- Design And Manufacture Of Integrated Circuits (AREA)
Abstract
The invention designs a damping circuit optimization design method for inhibiting numerical oscillation in simulation, and belongs to the field of electromagnetic transient. The invention aims to solve the problem of numerical oscillation generated by non-state variables of energy storage elements such as an inductor, a capacitor and the like through an action of a switch device in an electromagnetic transient simulation platform without an interpolation function. The core technical scheme provided by the invention is that through the damping circuit parameter optimization method designed by the invention, for different simulation circuits, different damping circuit configuration schemes are designed according to specific circuit parameters and topological structures. A basic process comprises the steps of 1, determining the selection of an initial value of a damping circuit; 2, searching for feature information representing an oscillation degree through a custom module; and 3, designing an iterative loop by adopting a simplex module, continuously updating R and C values through multi-time nonlinear iterations, and observing characteristic information to determine R and C values when the inhibition effect is optimal, thereby achieving the purpose of damping circuit parameter optimization.
Description
Technical field
The invention belongs to electromagnetic transient simulation field, be particularly useful for the numerical oscillation problem solving occur in emulation.
Background technology
At present, electromagnetic transient in power system simulation calculation software is varied, and the occasion of application is also not quite similar. For different types of application demand, off-line simulation platform and real-time simulator generally can be divided into. Off-line simulation platform includes various common software kit, such as PSCAD/EMTDC, EMTP-RV, ATP, MicroTran etc. The representative of real-time simulator mainly has RTDS, RT-LAB etc., is widely used in industry and academia. Although application places is different, but the integration method that whether off-line simulation real-time simulation adopts is mostly hiding-trapezium integral method. Hiding-trapezium integral method has the advantages such as precision height, good stability, but the impact due to its integration feature, in simulation process, switch and device activity, network structure changes, the non-state variable that can cause the energy-storage travelling wave tubes such as inductance capacitance swings abnormally after an event occurs near true solution, namely the numerical oscillation phenomenon in electromagnetic transient simulation.
Summary of the invention
The present invention proposes a kind of suitable in without the RC antihunt circuit method for designing suppressing numerical oscillation the emulation platform of interpolation function. First start with from the mechanism analyzing numerical oscillation, it is determined that the initial value configuration of antihunt circuit, and propose the quantization function describing numerical oscillation severe degreeSUM, and the configuration to RC circuit that it combined with simplex algorithm is optimized, and substantially increases the inhibition of antihunt circuit logarithm value vibration.
Technical scheme feature comprises the following steps:
Step 1: analyze the mechanism that numerical oscillation produces, be directed to different circuit parameters, in conjunction with the integration method that emulation platform uses, it is considered to circuit parameter cooperates with simulation time step-length, it is determined that the initial value configuration of RC circuit.
Step 2: electric parameters oscillatory condition when producing according to numerical oscillation, uses Fortran language to write custom block, generates the characteristic signal for judging degree of oscillation.
Step 3: characteristic signal step 2 produced, as input, combines with simplex module, design iteration ring obtains the optimal allocation of RC circuit by nonlinear iteration.
The present invention passes through above three step, it is possible to the initial value of RC antihunt circuit is determined, and by simplex algorithm, its parameter can be optimized further, and finds optimal solution, effectively solves the numerical oscillation problem occurred in emulation.
Accompanying drawing explanation
Fig. 1 is that the RL adding RC antihunt circuit tests circuit. In figureE s ��R��LIt is the RL design parameter configuration testing circuit,DFor diode,R snb ��C snb It it is the design parameter configuration of antihunt circuit.
Fig. 2 is the Dai Weinan flexible circuit of antihunt circuit,R SNBFor the equivalent resistance of antihunt circuit,V SNB_EQFor the equivalent history voltage source of antihunt circuit, its value was determined by the capacitance voltage in a upper moment, capacitance current, damping resistance and simulation step length.
Fig. 3 is characteristic signalSUMThe generation flow chart of function. In figureV LIt it is inductance in Fig. 1LThe voltage at two ends,I LIt is that Fig. 1 flows through inductanceLElectric current.
Fig. 4 is the Optimizing Flow figure of RC antihunt circuit, wherein quantization functionSUMIt is the characteristic signal described in Fig. 2, what deserves to be explained is, along with quantization functionSUMContinuous renewal, the value of antihunt circuit R and C also can update therewith, thus obtaining the optimal allocation of RC circuit.
Detailed description of the invention
A kind of antihunt circuit method for designing based on optimized algorithm simplex that the present invention relates to will be elaborated below. It is emphasized that the description below is merely exemplary, rather than in order to limit the scope of the present invention and application thereof.
The technical problem to be solved is the parameter by optimal damping circuit, the numerical oscillation problem occurred farthest is suppressed, in emulation when not by additional simulation algorithm so that simulation result is more accurately credible. The present invention adopts the following technical scheme that realization:
The present invention is realized by following three steps:
Step 1: as it is shown in figure 1, antihunt circuit to be connected in parallel on diode two ends, next determine the initial value of antihunt circuit R and C.
The initial value of R is determined: based on the theory of two-value resistance, diode current flow resistance is only small, generally takesR on=0.01 ��; Diode closes resistance break very greatly, generally takesR off=1M ��. Obviously, no matter diode is in which kind of state of break-make, RC antihunt circuit all should reduce its impact that diode is caused as far as possible, that is,R snbShould too small in case when diode current flow shunting or diode turn off time through-flow,R snbShould be not excessive yet, otherwise will be withR offEffect similar to such an extent as to weaken transient state damping function effect. Based on above analysis, here willR snbInitial value be taken as 5k �� roughly, be just between the order of magnitude of 0.01 �� and 10M ��.
CInitial value determine: as shown in Figure 2 at shutdown moment, the voltage at diode two ends also can be equivalent to an equivalent resistanceR SNBWith an equivalent history voltage sourceV SNB_EQSeries connection. It is found that in order to make electric currenti(t+ �� t) is not really zero,R SNBShould be withR snbTake same number level, should be too big, again because emulation scenarios discussed here is mainly the step-length off-line simulation more than 10 �� s, therefore hereCInitial value be chosen for 0.02 �� F roughly. It is inferred that along with the carrying out of later stage Optimization Work, in order to embody the effect of RC antihunt circuit,C snbChoosing of value will upwards restrain, because nowC snbThe initial value of=0.02 �� F is a conservative configuration being similar to open circuit, and in order to suppress the abnormal fluctuation of inductive drop, electric capacityC snbValue just it sufficiently large should be suppressed to fluctuate. Certainly,C snbValue also should be too big because once greater than certain limit, antihunt circuit is equivalent to high frequency electric presents a short-circuit state, even if the later stage adds optimized algorithm, it is also possible to its exact value cannot be converged to.
Step 2: use Fortran language to write custom blockSUM, hereSUMNamely functional value characterizes the severe degree of numerical oscillation. As it is shown on figure 3, inductive drop and inductive current are come in as raw information collection, it is processed further, the quantization function of outputSUMValue, namely as characteristic information, is the input information of next step iteration renewal RC value. FunctionSUMSpecific algorithm is as the formula (1).
(1)
Step 3: the quantization function that step 2 is obtainedSUMIt is worth as input quantity, adopts simplex to optimize module, by flow scheme design program as shown in Figure 4, write iteration loop. By ceaselessly iteration, according toSUMMost suitable RC value is found in the change of value so that suppress the effect of vibration to reach the best.
The above; being only the present invention preferably detailed description of the invention, but protection scope of the present invention is not limited thereto, any those familiar with the art is in the technical scope that the invention discloses; the change that can readily occur in or replacement, all should be encompassed within protection scope of the present invention. Therefore, protection scope of the present invention should be as the criterion with 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.
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)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009169733A (en) * | 2008-01-17 | 2009-07-30 | Hitachi Ltd | Fluid analyzing device 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 | 国家电网公司 | Modularized current-limiting breaker power module |
-
2015
- 2015-10-16 CN CN201510665187.1A patent/CN105631086A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009169733A (en) * | 2008-01-17 | 2009-07-30 | Hitachi Ltd | Fluid analyzing device 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 | 国家电网公司 | Modularized current-limiting breaker power module |
Non-Patent Citations (3)
| Title |
|---|
| 商莹等: "非线性电路暂态仿真中消除数值振荡的改进方法", 《电力系统保护与控制》 * |
| 王毅等: "消除电力系统故障仿真中数值振荡的方法研究", 《华东交通大学学报》 * |
| 许建中等: "超大规模MMC 电磁暂态仿真提速模型", 《中国电机工程学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106849057B (en) | Distributed wind power supply optimization method based on modern interior point method and Sensitivity Analysis Method | |
| CN106548265B (en) | Power transmission network reliability assessment method based on cascading failure accident chain search | |
| CN107230999B (en) | Regional distributed photovoltaic maximum capacity access evaluation method | |
| CN114285090B (en) | New energy source limit absorption capacity assessment method based on single station-partition-whole network | |
| CN107092984A (en) | A kind of network function end node propagation prediction method based on cascading failure | |
| CN104102956A (en) | Distribution network expansion planning method based on strategy adaption differential evolution | |
| CN114336608B (en) | Unit blocking optimization method and system considering dynamic capacity increase and reconstruction | |
| CN108229021A (en) | Modularization multi-level converter locking modeling method based on Real Time Digital Simulator | |
| CN104166880A (en) | Power system cascading failure mode prediction method | |
| Wang et al. | A review of power system transient stability analysis and assessment | |
| CN104868468A (en) | UPFC optimization configuration method based on life cycle cost | |
| CN105048473A (en) | Low-voltage reactive equipment voltage sensitivity statistics method in automatic voltage control of transformer station | |
| Zeng et al. | Resilience assessment for power systems under sequential attacks using double DQN with improved prioritized experience replay | |
| CN110797889B (en) | Energy storage power station arrangement method for solving tidal current congestion problem | |
| CN111651878A (en) | Large power grid static voltage stability optimization decision method and system considering situation assessment | |
| CN105631086A (en) | Damping circuit optimization design method for inhibiting numerical oscillation in simulation | |
| CN105373862A (en) | Terminal station planning scheme reliability evaluation optimal selection method based on annual load curve | |
| CN104868481A (en) | Capacitance compensating method for minimizing whole-day active power loss of distribution network | |
| CN103390249A (en) | Power distribution scheduling aid decision making method based on multiple dimensions | |
| CN105552906B (en) | A kind of area power grid load nargin analysis method based on prim al- dual interior point m ethod | |
| Li et al. | Optimal Attack Strategy of Power Grid based on Double Q-learning Algorithm | |
| CN110797906B (en) | Method and device for determining maximum installed capacity of wind power base and storage medium | |
| Li et al. | Deep reinforcement learning-based fast prediction of strategies for security control | |
| Zhou et al. | Distribution network reconfiguration considering the random character of wind power generation | |
| CN102841971B (en) | GIS (gas insulated switchgear)-based distribution network tuning parameter processing method and device |
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 |
Application publication date: 20160601 |
|
| WD01 | Invention patent application deemed withdrawn after publication |