CN105375482A - Power-supply network flow calculation method for tramcar - Google Patents
Power-supply network flow calculation method for tramcar Download PDFInfo
- Publication number
- CN105375482A CN105375482A CN201510939804.2A CN201510939804A CN105375482A CN 105375482 A CN105375482 A CN 105375482A CN 201510939804 A CN201510939804 A CN 201510939804A CN 105375482 A CN105375482 A CN 105375482A
- Authority
- CN
- China
- Prior art keywords
- voltage
- locomotive
- power
- model
- regenerative braking
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a power-supply network flow calculation method for a tramcar. The power-supply network flow calculation method comprises steps of setting a traction substation initial state as Switch on, wherein a locomotive initial model is a constant power source and initial voltage is system floating voltage, utilizing a current vector iterative method to solve a node voltage equation, performing detection on DC bus voltage and regenerative braking locomotive voltage, if the regenerative braking locomotive voltage is greater than Umax or the DC bus voltage of the traction substation is greater than U0+Beta, automatically modifying the model by a program, replacing the regenerative braking locomotive, the voltage of which crosses the limit, with a power supply modified model, solving the new node voltage equation by the program after the model is modified until the voltage of any node does not cross the limit, comparing the power of each generative braking locomotive with the available braking power, and modifying the locomotive model according to a formulated limitation condition. The invention expresses the overvoltage regenerative braking locomotive through the power supply correction model, brings forward an iteration method for modifying internal resistance and thus guarantees the convergence of the flow calculation.
Description
Technical field:
The present invention relates to net Load flow calculation technical field, relate to a kind of power supply network tidal current computing method about tramcar in particular.
Background technology:
Urban tram DC Traction Network tidal current computing method can be divided into two kinds: the linear system tidal current computing method based on constant current supply model locomotive and the non linear system tidal current computing method based on constant power source model locomotive.Adjacent locomotives can not be taken into account due to locomotive by constant current supply modeling and get influencing each other of stream, based on constant power source model locomotive direct-current traction power supply emulation in closer to truth, therefore, present most of direct-current traction power supply software for calculation adopts constant power source model locomotive.But this model can make power flow equation non-linearization, therefore solve by the method for iteration.When calculating with iterative algorithm the traction power network flow considering regenerative braking, situation will become complicated, and regenerating braking energy superfluous in system causes city rail DC Traction Network Load flow calculation to restrain.Therefore, the problem causing city rail DC Traction Network Load flow calculation to restrain for regenerating braking energy superfluous in prior art is improved.
Summary of the invention:
Object of the present invention is exactly the deficiency for prior art, a kind of power supply network tidal current computing method of tramcar is provided, it represents overvoltage regenerative braking locomotive by power supply correction model, proposes a kind of alternative manner revising internal resistance, thus ensure that the convergence of Load flow calculation.
Technology solution of the present invention is as follows:
A power supply network tidal current computing method for tramcar, comprises the following steps,
A) electric traction institute initial condition is set to Swichon, the locomotive initial model of tramcar is constant power source, and initial voltage is system floating voltage;
B) current phasor iterative method is utilized to obtain the solution of nodal voltage equation;
C) DC bus-bar voltage of electric traction and regenerative braking locomotive voltage are tested, if regenerative braking locomotive voltage is greater than U
maxor traction institute DC bus-bar voltage is greater than U
0+ β, program will revise model automatically, and replaced by the regenerative braking electricity for locomotive source correction model of voltage out-of-limit, after having revised model, program solves new nodal voltage equation again, until arbitrary node voltage is not out-of-limit, U
maxfor voltage threshold limit value, β is the change in voltage allowance allowed.
D) braking power available with it for each regenerative braking locomotive power is compared, according to the limit value condition amendment model locomotive formulated;
E) terminate to calculate.
Preferred as technique scheme, described step c) in power supply correction model be equivalent to the parallel form of current source and internal resistance by Thevenin's theorem.
Preferred as technique scheme, described steps d) in the limit value condition formulated be,
(1) P
brake≤ P
(k), braking power is less than Traction networks demand power, because power supply correction model is to Traction networks injecting power P
(k)time, output voltage can not more than U
max, therefore locomotive is to Traction networks injecting power P
braketime, voltage also can not be out-of-limit, at this moment should change back constant power source by regenerative braking model locomotive, ensures that whole braking energies of this car inject electrical network as regenerated energy;
(2) P
brake>P
(k), braking power is greater than Traction networks demand power, now locomotive brake energy surplus, and unnecessary energy should be braked resistance absorption.The condition that brake resistance drops into is that locomotive voltage reaches U
max, locomotive voltage is tested, if U>U
max-α (α is taken as 1V), so this result meets the requirements, and calculates and terminates; Otherwise, just need to reduce internal resistance, make voltage bring up to U
maxmore than-α.
Beneficial effect of the present invention is: adopt power supply correction model to represent overvoltage regenerative braking locomotive, power supply correction model is in series by ideal voltage source and internal resistance regulating, internal resistance value is adjusted with alternative manner, the locomotive voltage of braking energy surplus can be made to approach locomotive maximum allowable voltage fast, force brake resistance to drop into and absorb superfluous braking energy, make the locomotive voltage of braking energy surplus reach the voltage levvl that brake resistance is dropped into fast, DC Traction Network Load flow calculation can reliable conveyance.
Accompanying drawing illustrates:
Fig. 1 is the flow chart of computational methods of the present invention.
Embodiment:
Embodiment: by particular specific embodiment, embodiments of the present invention are described below, those skilled in the art can understand other advantages of the present invention and effect easily by content disclosed in the present specification.
As shown in Figure 1, a kind of power supply network tidal current computing method of tramcar, comprises the following steps,
A) electric traction institute initial condition is set to Swichon, the locomotive initial model of tramcar is constant power source, and initial voltage is system floating voltage;
B) current phasor iterative method is utilized to obtain the solution of nodal voltage equation;
C) DC bus-bar voltage of electric traction and regenerative braking locomotive voltage are tested, if regenerative braking locomotive voltage is greater than U
maxor traction institute DC bus-bar voltage is greater than U
0+ β, program will revise model automatically, and replaced by the regenerative braking electricity for locomotive source correction model of voltage out-of-limit, after having revised model, program solves new nodal voltage equation again, until arbitrary node voltage is not out-of-limit, U
maxfor voltage threshold limit value, β is the change in voltage allowance allowed.
D) braking power available with it for each regenerative braking locomotive power is compared, according to the limit value condition amendment model locomotive formulated;
E) terminate to calculate.
Concrete computational process is as follows, and DC traction power-supply system can describe with nodal voltage equation below:
[G][U]
(n+1)=[I]
(n)
In formula, n is for calculating iterations, and [G] is the node admittance matrix of m × m, and m is interstitial content, [I]
(n)be n-th iteration node Injection Current column vector, [U]
(n+1)for amount to be asked.
In order to form nodal voltage equation, traction voltage source series resistance model Thevenin's theorem be equivalent to the form of current source shunt resistance.Model locomotive is also converted to ideal current source by constant power source, and its value following formula is asked for.
In formula, i is locomotive unit number, P
ifor i locomotive power, invariable in an iterative process, its value provides by drawing calculating, hereafter this power is called available horsepower.
[I] is the node current column vector corresponding with [G], is made up of traction institute's Injection Current and locomotive Injection Current.Locomotive initial voltage is system floating voltage, and the size of described floating voltage is 1600V, by formula
the locomotive initial current obtained is brought into formula [G] [U]
(n+1)=[I]
(n), obtain new locomotive voltage, then new locomotive voltage brought into formula
upgrade locomotive electric current, until train voltage deviation reaches License Value, iteration stopping.Nodal voltage equation is system of linear equations, adopts QR decomposition method direct solution nodal voltage equation group in the present embodiment.
Represent overvoltage regenerative braking locomotive by power supply correction model in the present embodiment, this model changes power output and voltage by amendment internal resistance, thus restriction excess energy injects Traction networks.Power source model is in series by ideal voltage source and adjustable internal resistance, and voltage source numerical value is set to locomotive maximum permissible voltage U
max, described maximum permissible voltage U
maxsize be 1800V, can ensure that locomotive voltage can not more than U like this
max.Internal resistance initial value is tried to achieve by following formula:
In formula, U
0for traction institute floating voltage, described floating voltage size is 1600V, P
brakefor by drawing the available brake power of this locomotive calculated, wherein, injecting electrical network direction is just.By formula
the Re obtained
(0)the output voltage of power supply can be made to equal U
0time power output be P
brake.
The functional relation of output power of power supply and output voltage is:
By formula
find out, output power of power supply is the binary function of locomotive voltage and internal resistance.When power output P is constant, reduces internal resistance value and can improve output voltage U.
The energy that regenerative braking locomotive injects Traction networks is regenerated energy.Work as Re=Re
(k)time, the regenerative braking locomotive based on power source model injects the power P of Traction networks
(k)=(U
max-U
0) U
0/ Re
(k).Then judge that whether locomotive brake energy is superfluous.
(1) P
brake≤ P
(k), braking power is less than Traction networks demand power, because power supply correction model is to Traction networks injecting power P
(k)time, output voltage can not more than U
max, therefore locomotive is to Traction networks injecting power P
braketime, voltage also can not be out-of-limit, at this moment should change back constant power source by regenerative braking model locomotive, ensures that whole braking energies of this car inject electrical network as regenerated energy.
(2) P
brake>P
(k), braking power is greater than Traction networks demand power, now locomotive brake energy surplus, and unnecessary energy should be braked resistance absorption.The condition that brake resistance drops into is that locomotive voltage reaches U
max, locomotive voltage is tested, if U>U
max-α (α is taken as 1V), so this result meets the requirements, and calculates and terminates; Otherwise, just need to reduce internal resistance, make voltage bring up to U
maxmore than-α.
When output power of power supply is constant, at interval U ∈ (U
max/ 2, U
max) on, output voltage can be expressed as the monotropic function of internal resistance.
Because this function is non-linear, so with alternative manner modified R e.Iterative process is: first give Re initialize, again forms nodal voltage equation and asks for U, then according to U and U
maxdifference modified R e, through successive ignition, finally make U>U
max-α.
Formula
at interval U ∈ (U
max/ 2, U
max) on inverse function be:
Re (U) is elementary function, so at interval U ∈ (U
max/ 2, U
max) on can lead, and then can be micro-in any point in this interval.
So, according to formula
Differential form
Draw following iterative formula:
Wherein, U
(k), P
(k)for Re=Re
(k)time the trend solution that calculates.
In Practical Calculation process, different Re
(k)under the P that calculates
(k)change, but no matter P
(k)how becoming, the variation tendency of U-Re curve is constant, according to formula
The Re calculated
(k+1)always than Re
(k)little, be equivalent to improve power supply capacity, so just can reach the effect improving electric power output voltage.
Traction institute's rectifying device conducting and cut-off state is represented with Switchon and Switchoff in Load flow calculation.When traction institute, busbar voltage is greater than U
0during+β, this traction institute internal resistance is set to a very large number, makes high-impedance state that traction becomes, i.e. Switchoff state; Otherwise change back Switchon by traction institute state.
Described embodiment in order to illustrative the present invention, but not for limiting the present invention.Any those skilled in the art all without prejudice under spirit of the present invention and category, can modify to described embodiment, therefore the scope of the present invention, should listed by claim of the present invention.
Claims (3)
1. a power supply network tidal current computing method for tramcar, is characterized in that: comprise the following steps,
A) electric traction institute initial condition is set to Switchon, the locomotive initial model of tramcar is constant power source, and initial voltage is system floating voltage;
B) current phasor iterative method is utilized to obtain the solution of nodal voltage equation;
C) DC bus-bar voltage of electric traction and regenerative braking locomotive voltage are tested, if regenerative braking locomotive voltage is greater than U
maxor traction institute DC bus-bar voltage is greater than U
0+ β, program will revise model automatically, and replaced by the regenerative braking electricity for locomotive source correction model of voltage out-of-limit, after having revised model, program solves new nodal voltage equation again, until arbitrary node voltage is not out-of-limit, U
maxfor voltage threshold limit value, β is the change in voltage allowance allowed.
D) braking power available with it for each regenerative braking locomotive power is compared, according to the limit value condition amendment model locomotive formulated;
E) terminate to calculate.
2. the power supply network tidal current computing method of a kind of tramcar according to claim 1, is characterized in that: described step c) in power supply correction model be equivalent to the parallel form of current source and internal resistance by Thevenin's theorem.
3. the power supply network tidal current computing method of a kind of tramcar according to claim 1, is characterized in that: described steps d) in formulate limit value condition be,
(1) P
brake≤ P
(k), braking power is less than Traction networks demand power, because power supply correction model is to Traction networks injecting power P
(k)time, output voltage can not more than U
max, therefore locomotive is to Traction networks injecting power P
braketime, voltage also can not be out-of-limit, at this moment should change back constant power source by regenerative braking model locomotive, ensures that whole braking energies of this car inject electrical network as regenerated energy;
(2) P
brake>P
(k), braking power is greater than Traction networks demand power, now locomotive brake energy surplus, and unnecessary energy should be braked resistance absorption.The condition that brake resistance drops into is that locomotive voltage reaches U
max, locomotive voltage is tested, if U>U
max-α, α are taken as 1V, and so this result meets the requirements, and calculate and terminate; Otherwise, just need to reduce internal resistance, make voltage bring up to U
maxmore than-α.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510939804.2A CN105375482A (en) | 2015-12-15 | 2015-12-15 | Power-supply network flow calculation method for tramcar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510939804.2A CN105375482A (en) | 2015-12-15 | 2015-12-15 | Power-supply network flow calculation method for tramcar |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105375482A true CN105375482A (en) | 2016-03-02 |
Family
ID=55377399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510939804.2A Pending CN105375482A (en) | 2015-12-15 | 2015-12-15 | Power-supply network flow calculation method for tramcar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105375482A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055852A (en) * | 2016-07-21 | 2016-10-26 | 成都运达科技股份有限公司 | Subway train operation and traction power supply system coupled simulation modeling method |
CN106655196A (en) * | 2017-02-15 | 2017-05-10 | 河海大学 | Supply network power flow calculation method for rail transit |
CN110034556A (en) * | 2019-03-27 | 2019-07-19 | 沈阳微控新能源技术有限公司 | Regenerating braking energy recycling and control method, system and device |
CN111967708A (en) * | 2020-06-29 | 2020-11-20 | 中铁第一勘察设计院集团有限公司 | Operation stability evaluation method and device based on harmonic resonance of traction power supply system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103390111A (en) * | 2013-07-29 | 2013-11-13 | 北京交通大学 | Calculation method oriented to single-train running energy consumption of urban rail transit |
CN104598700A (en) * | 2015-03-03 | 2015-05-06 | 南车株洲电力机车研究所有限公司 | Method and system for equivalent modeling of train |
-
2015
- 2015-12-15 CN CN201510939804.2A patent/CN105375482A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103390111A (en) * | 2013-07-29 | 2013-11-13 | 北京交通大学 | Calculation method oriented to single-train running energy consumption of urban rail transit |
CN104598700A (en) * | 2015-03-03 | 2015-05-06 | 南车株洲电力机车研究所有限公司 | Method and system for equivalent modeling of train |
Non-Patent Citations (2)
Title |
---|
刑华栋: "《城市轨道交通供电仿真***的研究与开发》", 《中国优秀硕士学位论文全文数据库》 * |
马斌: "《Matlab语言及实践教程》", 31 January 2013 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055852A (en) * | 2016-07-21 | 2016-10-26 | 成都运达科技股份有限公司 | Subway train operation and traction power supply system coupled simulation modeling method |
CN106055852B (en) * | 2016-07-21 | 2019-03-12 | 成都运达科技股份有限公司 | A kind of operation of subway train and tractive power supply system coupling Simulation modeling method |
CN106655196A (en) * | 2017-02-15 | 2017-05-10 | 河海大学 | Supply network power flow calculation method for rail transit |
CN110034556A (en) * | 2019-03-27 | 2019-07-19 | 沈阳微控新能源技术有限公司 | Regenerating braking energy recycling and control method, system and device |
CN111967708A (en) * | 2020-06-29 | 2020-11-20 | 中铁第一勘察设计院集团有限公司 | Operation stability evaluation method and device based on harmonic resonance of traction power supply system |
CN111967708B (en) * | 2020-06-29 | 2024-06-07 | 中铁第一勘察设计院集团有限公司 | Operation stability assessment method and device based on traction power supply system harmonic resonance |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106655196A (en) | Supply network power flow calculation method for rail transit | |
Aguado et al. | Optimal operation of electric railways with renewable energy and electric storage systems | |
Iannuzzi et al. | Stationary ultracapacitors storage device for improving energy saving and voltage profile of light transportation networks | |
CN106055852B (en) | A kind of operation of subway train and tractive power supply system coupling Simulation modeling method | |
Jabr et al. | Solution of DC railway traction power flow systems including limited network receptivity | |
Arboleya et al. | BFS algorithm for voltage-constrained meshed DC traction networks with nonsmooth voltage-dependent loads and generators | |
CN105375482A (en) | Power-supply network flow calculation method for tramcar | |
CN105303021A (en) | Adaptive real-time dynamic mathematic modeling method for direct-current traction power supply system in urban railway traffic | |
CN113300395B (en) | Hybrid energy storage optimal capacity configuration method for flexible traction power supply system | |
CN106549377B (en) | A kind of handover method of alternating current-direct current combined hybrid system electromechanics-electromagnetic transient hybrid simulation | |
Nian et al. | Transient modeling and analysis of VSC based DC microgrid during short circuit fault | |
Tabari et al. | A mathematical model for a stability-enhanced DC distribution system for power system integration of plug-in electric vehicles | |
CN107590343A (en) | A kind of track traffic full line analogue system and emulation mode | |
CN105514987A (en) | 10kV distribution circuit distributed photovoltaic constant volume method and system | |
CN104505866B (en) | Equivalent decoupling method for simulating multi-feed-in direct current fault recovery characteristics | |
CN107884680B (en) | The calculation method of transient under multiterminal flexible direct current system fault condition | |
Jung et al. | A study on peak power reduction using regenerative energy in railway systems through DC subsystem interconnection | |
CN104283208B (en) | The composition decomposition computational methods of the probability available transmission capacity of large-scale power grid | |
CN110048427B (en) | Multi-end flexible direct-current power distribution network power flow protection constant value calculation method | |
CN104155553B (en) | Locomotive electric power feedback equipment test system and test method | |
CN103390890B (en) | Based on the distribution power flow analytical method of current distribution factor | |
Graber et al. | Impact Assessment of Energy Storage Systems Supporting DC Railways on AC Power Grids | |
CN115566685A (en) | Urban rail transit tide tracking method and device | |
CN113852079A (en) | Inversion feedback type regeneration device simulation method suitable for power flow simulation calculation | |
CN103647285B (en) | Based on the cyclization voltage control feasible zone defining method of two-port Thevenin's equivalence |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160302 |
|
RJ01 | Rejection of invention patent application after publication |