CN105375482A - Power-supply network flow calculation method for tramcar - Google Patents

Power-supply network flow calculation method for tramcar Download PDF

Info

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
Application number
CN201510939804.2A
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201510939804.2A priority Critical patent/CN105375482A/en
Publication of CN105375482A publication Critical patent/CN105375482A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject 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

A kind of power supply network tidal current computing method of tramcar
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.
I i ( n ) = P i / U i ( n )
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:
Re ( 0 ) = P B r a k e ( U m a x - U 0 ) U 0
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:
P ( U ) = - 1 Re ( U - U m a x 2 ) 2 + U m a x 2 4 Re
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.
U ( Re ) = U m a x 2 4 - Re · P + U m a x 2
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 ) = - 1 P U 2 + U m a x P U
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 Re ( U ) = - 1 P U 2 + U m a x P U Differential form
d Re = ( - 2 P U + U m a x P ) d U
Draw following iterative formula:
Re ( k + 1 ) = Re ( k ) + ( - 2 P ( k ) U ( k ) + U m a x P ( k ) ) ( U m a x - U ( k ) )
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 Re ( k + 1 ) = Re ( k ) + ( - 2 P ( k ) U ( k ) + U m a x P ( k ) ) ( U m a x - U ( k ) ) 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-α.
CN201510939804.2A 2015-12-15 2015-12-15 Power-supply network flow calculation method for tramcar Pending CN105375482A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
刑华栋: "《城市轨道交通供电仿真***的研究与开发》", 《中国优秀硕士学位论文全文数据库》 *
马斌: "《Matlab语言及实践教程》", 31 January 2013 *

Cited By (6)

* Cited by examiner, † Cited by third party
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