Invention content
It is an object of the present invention to overcome the above deficiencies, it is proposed that based on efficiency and stream deviation it is expected weighted sum square
The parallel operation system optimized control method of battle array.
The technical scheme is that:A kind of based on efficiency and stream deviation it is expected the parallel operation system of weighted sum matrix
Optimal control method, its step are as follows:
(1) the parallel operation system load electric current I of K power module composition is obtainedoutFromIt is divided into according toEquidistantly change toWhen, each power module is in different loads electric currentIn the case of adopt
Collect V output current Datacurr(m') (i) (j), output voltage Datavolt(m') (i) (j) and input power P (m') (i) (j),
And obtain corresponding V power module efficiency eta (m') (i) (j);Wherein:M' is power module serial number;I is load current value pair
The sequence number value answered;J is that output current acquires serial number;M', i, j meet m'={ 1 ... K }, i={ 1 ... U }, j={ 1 ... V }, IN
For the rated current of power module;
(2) output current and stream expectation electric current relative deviation absolute value of power module are obtainedAnd corresponding mathematic expectaionObtain power module
The average value of efficiencyAnd efficiency etam'iWithWeighted suma∈[0,
1];With γm'iMatrix A=(γ is built for elementm'i)K×U;M', i, j meet m'={ 1 ... K }, i={ 1 ... U }, j={ 1 ...
V};
(3) solution matrix A=(γm'i)K×U1- normsCorresponding load current serial number a
Value, i.e.,Corresponding a values are set up, wherein:a∈[1,U],i∈[1,U];The a values are power module
Load current isWhen system efficiency and stream integrated performance index it is optimal;And it obtains efficiency and flows comprehensive performance and refer to
The output current for marking the power module under optimal situation is Iref:
(4) with cycle TsFor the online power module quantity M of interval calculation parallel operation system, and to M online power supplys
The output current of module is acquired, and the output current data of the online power module of m-th of serial number is labeled as Curr (m), m
For the serial number of current online power module;
(5) the output current data array of the online power module of serial number m is obtained:Curr_store (m) (n)=
Curr_store (m) (n+1), Curr_store (m) (T)=Curr (m);Wherein:N=1 ... T-1;M=1,2,3 ... M;T is
Positive integer more than 2, n are the current sample number of the online power module of current sequence number;
(6) the output current average value of the online power module of serial number m is obtained:
Wherein:M=1,2,3 ... M;
(7) load current of the parallel operation system of M online power module compositions is obtainedWith
The equal current load electric current of line power module
(8) judge | Ishare-IrefWhether |≤σ is true;
(9) in step (8) | Ishare-Iref|≤σ is invalid, and it is reference current I to obtain online power module output currentref
When online power module quantity N*, i.e.,N*≤ 1 is arranged N*=2;Conversely, then obtaining parallel operation
System need to adjust online power module amount Δ N*=N*- M, according to Δ N*It is positive and negative, Centralized Controller increases or decreases | Δ N*|
A online power module;
(10) in step (8) | Ishare-Iref|≤σ is set up, then obtains the output current of the online power module of serial number m
Curr_store (m) (n) and stream desired value IshareDeviation θ (m) (n)=Curr_store (m) (n)-Ishare;Wherein:n
=1 ... T;M=1,2,3 ... M;
(11) mathematic expectaion of the online power module deviation θ (m) (n) of serial number m is obtainedIts
In:N=1 ... N;M=1,2,3 ... M;
(12) judgeEstablishment then continues next online power module detection, and invalid then marking serial numbers are m
Online power module current sharing can undesirable, CθIt is expected for deviationMaximum permissible value;
(13) Num label is offline for the undesirable online power module of energy, and from stand-by power supply
Start Num power module work, and continue the operation of step (4), Num is can be undesirable online labeled as current sharing
The quantity of power module.
In step (1)-(2),
(1) t ∈ ((i-1) T, iT], (U >=i >=1), electronic load current isWhen, obtain power supply mould
The stream target reference current of block is:
(2) the power module output current sampled data data of serial number m' are obtained:Datacurr(m') (i) (j), (K >=
M' >=1, U >=i >=1, V >=j >=1), and obtain it and flow relative deviation absolute value δ (m') (i) (j) and be:
(3) power module for obtaining serial number m' existsUnder the conditions of mathematics phases of the δ (m') (i) (j) about j
Hope Em'i:
Em'iIndicate that the power module of serial number m' existsUnder the conditions of relative deviation absolute value average value.
(4) power module for obtaining serial number m' existsCondition efficiency eta (m') (i) (j):
(5) power module for obtaining serial number m' existsUnder the conditions of mathematics phases of the η (m') (i) (j) about j
Hope ηm'i:
ηm'iIndicate that the power module of serial number m' existsUnder the conditions of efficiency average value;
(6) η is obtainedm'iWithWeighted sum γm'iFor:
(7) structure is with γm'iFor matrix A=(γ of elementm'i)K×U:
Therefrom obtain matrix A=(γm'i)K×U1 norm | | A | |1:I=1,2 ... U;
(8) enable | | A | |1Corresponding i values are a, from inequality:
A ∈ [1, U], i=1,2 ... U
The load current I that power module is corresponding when obtaining parallel operation system effectiveness and stream comprehensive performance is optimalref:
The present invention first, obtains the parallel operation system being made of K power module and (covers light in different loads
Carry, semi-load, nominal load and overload conditions) number of each power module efficiency eta and actual current value and the equal flow valuve deviation of expectation
Term hopes Ε;Secondly, seek in the case of different loads power module η andWeighted sumA ∈ [0,1] and with
γ is that element builds K power module in the matrix A of different loads, calculates 1 norm of the matrix | | A | |1And it is corresponding
In the case of efficiency and stream integrated performance index it is optimal when module output current Iref;Again, parallel operation system centralized control
Device seeks the total load current value I of system by obtaining the output current for the power module that works online in real timeoutAnd it divided by obtains in advance
The I obtainedref, obtain parallel operation system on-line operation power module quantity N*;Then, parallel operation system Centralized Controller is real
When control online power module quantity, it is ensured that system work in always current sharing energy optimum point or near;Finally, online fortune is obtained
Capable each power module output current data, calculate the deviation of output current data and the equal flow valuve of target, obtain it
Mathematic expectaion, to differentiate whether the power module of each on-line operation meets the requirements and carry out optimal dispatch control.Due to phase
Its characteristic of the power module of same specification is totally consistent, thus (size of K can determine by user, K of the present invention by measuring K
It fixes tentatively efficiency of the parallel operation system for 10) a power module composition under different loads electric current and flows integrated performance index
The parallel operation system of arbitrary N number of power module composition be can be obtained in the efficiency of different loads and flow comprehensive performance
Index.Then, the mathematic expectaion of the output current and the equal flow valuve deviation of target of each power module of operation is obtained.Ensureing electricity
Control is scheduled to operation power module according to the size of mathematic expectaion under the premise of source module quantity is optimal, it is ensured that Mei Geyun
The performance of capable power module is met the requirements.
The present invention has the advantage that:
(1) present invention covers load current full operating range operating modes, have wide applicability;
(2) present invention can integrate and take into account parallel operation system effectiveness and current sharing energy index, with significant economy with
System reliability;
(3) present invention by efficiency in the case of building full load and stream deviation expectation weighted sum matrix A=
(γmi)K×UAnd calculate A=(γmi)K×U1 norm | | A | |1And its corresponding load current Iref。||A||1Characterize simultaneously alliance
Efficiency and stream integrated performance index is optimal and its corresponding in wire module load current value during electric system flows, for parallel connection
Power supply system efficiency and flow-optimized control provide foundation.
(4) present invention has dynamically adjusts online power module quantity in real time, it is ensured that parallel operation system works in always
Flow near optimal working point.
(5) present invention is in target near optimal working point ensuring parallel operation system effectiveness and flow comprehensive performance
It is that on-line operation power module quantity is optimal down, by the output current and stream desired value I that calculate each power modulesharePartially
Poor mathematic expectaion is unsatisfactory for desired online power module to performance and optimizes scheduling with standby power supply module, realizes in parallel
Power supply system and online power module work near best performance point;
(6) of the present invention to be controlled with the parallel operation system optimization of stream deviation expectation weighted sum matrix based on efficiency
Method has the features such as reliability is high, highly practical;Parallel operation system current sharing energy and efficiency index can be effectively taken into account, is improved
The performance driving economy and reliability of system provide reliable guarantee for the safe efficient operation of parallel operation system.
Specific implementation mode
The embodiment of the present invention is described further below for attached drawing:
The present invention provides the parallel operation system optimization controlling parties that based on efficiency and stream deviation it is expected weighted sum matrix
Method, Fig. 1 show parallel operation system construction drawing, and Fig. 2 show parallel operation system effectiveness and stream all-round property testing system
System structure chart.Fig. 1 includes mainly parallel operation system Centralized Controller, power module and power load.Centralized Controller passes through
Communication bus obtains the IP and its output current in wire module, the quantity and Optimized Operation performance of the online power module of optimal control
Underproof power module;Power module mainly realize the Centralized Controller that powers to the load, receives operation control command and on
Pass output current;Power load includes mainly all kinds of electrical equipments.That whether there is or not communication bus is independently equal for the realization of stream regulatory function
Stream mode and there is communication bus to flow mode, realized by special flow equalizing function module, the present invention does not repeat.Fig. 2 major functions
It is to obtain parallel operation system online module efficiency and stream deviation it is expected weighted sum matrix under different load conditions, and then it is true
The fixed corresponding optimal load electric current of 1 norm of matrix includes mainly host computer (PC machine), program-control electronic load, power module, work(
Rate meter etc..Host computer (PC machine) major function is to obtain online module I P address, input power, module output current, output work
Rate, control program-control electronic load operating current, computational efficiency and the corresponding load of 1 norm of stream deviation expectation weighted sum matrix are electric
Stream;Program-control electronic load is used to adjust the load current of parallel operation system;Power module, which is mainly realized, to be received IP settings, receives
Host computer order data and upload output current, output power are to host computer;Power meter is mainly used for measuring in the defeated of wire module
Enter power.
One, parallel operation system effectiveness and stream comprehensive performance testing system variable declaration are as follows:K tests for parallel operation
The occurrence of system power supply module number, K can be set according to actual conditions.INFor power module rated current;For simultaneously alliance
Electric system output-current rating meetsU is load current point quantity, i.e. parallel operation system load electric current IoutFromIt is divided into according toEquidistantly change to(cover underloading, semi-load, specified load and overload conditions,
U is necessary for the positive integer not less than 20, can be determined according to the maximum load current value that system works by user);
For electronic load at i-th output current, wherein:U≥i≥1;M' is power module serial number, is met:K power module
IP is mapped as m'=1 according to order from small to large, and 2 ... K, i.e. m'=1 are the power module serial number of IP minimums, m=2 IP
Secondary minimal power module serial number ..., and so on m'=K be the maximum power module serial numbers of IP;V is that parallel operation system is in
It need to be to ought individually online power supply power module output current, output voltage and input power data sampling when a certain load current point
Quantity, V can be sized according to actual needs.Datacurr(m') (i) (j), (K >=m' >=1, U >=i >=1, V >=j >=1) are sequence
Number exist for the power module of m'Under the conditions of j-th of current sampling data;Datavolt(m') (i) (j), (K >=m'
>=1, U >=i >=1, V >=j >=1) exist for the power module of serial number m'Under the conditions of the sampling of j-th output voltage
Data;P (m') (i) (j), (K >=m' >=1, U >=i >=1, V >=j >=1) are that the power module of serial number m' exists
Under the conditions of j-th of input power sampled data;η (m') (i) (j), (K >=m' >=1, U >=i >=1, V >=j >=1) are serial number m'
Power module existUnder the conditions of j-th of efficiency data calculating, meet:ηm'iPower module for serial number m' existsUnder the conditions of V η (m') (i) (j) mathematic expectaion, meet:Iref(i) it is module
?Under the conditions of flow target reference, meet:Wherein:U≥i≥1;δ(m')
(i) power module that (j) is serial number m' existsUnder the conditions of j-th sample rate current with flow reference target electric currentRelative deviation absolute value, meet:Em'iFor serial number m'
Power module existUnder the conditions of V δ (m') (i) (j) mathematic expectaion, meet:
α is weighting coefficient, is met:α∈[0,1];γm'iPower module for serial number m' existsUnder the conditions of efficiency and
It flows deviation and it is expected weighted sum, meet:A=(γm'i)K×UFor efficiency and stream deviation it is expected weighting
And matrix.
Define the last moment that t=0 is parallel operation system no-load running;T is two neighboring load current interval time;
Then t ∈ ((i-1) T, iT], (U >=i >=1) be parallel operation system load electric currentRun time.Due toIt needs to acquire 3V sample data to each power module in operational process, thus, host computer need to acquire 3 altogether
× K × V data.Assuming that the time that host computer acquires a data is T1, then system work inState needs
Ttotal=3 × K × V × T1Time, thus must satisfy T >=Ttotal.Again due to current sharing energy data reliability and sampling number
With sampling time T1Correlation, thus T and T need to be considered according to actual demand1Size, it is ensured that the reliability of current sharing energy index.
First, by control engineering knowledge it is found that the overshoot that the performance of evaluation system can be responded by system step, adjustment
Time and steady-state deviation index are weighed.Thus, parallel operation system electronic load byStep isWhen, we again may be by measuring the electric current output of power module and between stream target reference
Dynamic response evaluates the current sharing energy of power module.By mathematical statistics knowledge it is found that flowing the mathematics phase of relative deviation ordered series of numbers
Hope characterization is the global consistency between actual value and desired value, embodies the accuracy during its step response, can reflect
Power module current sharing energy index;Secondly, parallel operation system is while meeting current sharing energy index, it should take into account system fortune
Capable economic benefit.System operation economic benefit is determined by operational efficiency, embodies the energy consumption index of system in the process of running.
To make parallel operation system consider efficiency and flow index, efficiency need to be built and flow deviation mathematic expectaion weighted sum and its
Corresponding matrix;Finally, by seeking efficiency and stream deviation mathematic expectaion weighted sum matrix A=(γm'i)K×U1- norms ‖
A‖1Determine parallel operation overall system performance optimal value and its corresponding load current.This is to be based on A=(γm'i)K×U1- models
Number ‖ A ‖1Physical significance show which kind of load current is parallel operation system be under efficiency and stream comprehensive performance is best.
T ∈ ((i-1) T, iT], (U >=i >=1), electronic load current isThe then stream of power module
Target reference current is:
Obtain the power module output current sampled data data of serial number m':Datacurr(m') (i) (j), (K >=m' >=
1, U >=i >=1, V >=j >=1), thus, flowing relative deviation absolute value δ (m') (i) (j) is:
The power module for seeking serial number m' existsUnder the conditions of mathematic expectaions of the δ (m') (i) (j) about j
Em'iFor:
Em'iPhysical significance be:The power module of serial number m' existsUnder the conditions of relative deviation absolute value
Average value, Em'iIt is smaller show power moduleUnder the conditions of current sharing can be better.
The power module for seeking serial number m' existsCondition efficiency eta (m') (i) (j) is:
The power module for seeking serial number m' existsUnder the conditions of mathematic expectaions of the η (m') (i) (j) about j
ηm'iFor:
ηm'iPhysical significance be:The power module of serial number m' existsUnder the conditions of efficiency average value,
ηm'iShow more greatly power moduleUnder the conditions of economic performance it is better, it is more energy saving;
Seek ηm'iWithWeighted sum γm'iFor:
γm'iPhysical significance be:System effectiveness performance and current sharing can be unified to examine to a comprehensive performance indicator
Consider, is conducive to the Design of Performance of system;
Structure is with γmiFor matrix A=(γ of elementm'i)K×U:
Solution matrix A=(γm'i)K×U1 norm ‖ A ‖1:
||A‖1Physical significance be:The parallel operation system effectiveness that be made of K power module and stream comprehensive performance are most
It is excellent.
Enable ‖ A ‖1Corresponding i values are a, that is, meet inequality:
The then parallel operation system effectiveness load current I corresponding with power module when stream comprehensive performance is optimalrefFor:
Two, the figure variable declaration of parallel operation system optimization control structure is as follows:
TsOnline power module quantity is calculated for Centralized Controller and acquires the period of power module output current data;M
For online power module quantity;IoutFor the load current of parallel operation system;Curr (m) is the online power module of serial number m
Output current sampled value, m=1,2, ┄, M;M is the quantity of current online power module, IrefFor parallel operation system effectiveness
Load current corresponding to online power module when stream comprehensive performance is optimal;IshareFor parallel operation system operation when it is electric online
Source module output current flows desired value;Δ I is IshareWith IrefAbsolute value of the difference;σ is IshareWith IrefDifference absolute value
Maximum permissible value.Curr_store (m) (n) be serial number m online power module output current store array, m=1,2,
┄, M;N=1,2, ┄, T;N is the current acquisition number of the online power module of current sequence number, and T is the positive integer more than 2,The average value of array Curr_store (m) (n) is stored for the output current of the online power module of serial number m;θ(m)
(n) Curr_store (m) (n) for the online power module for being serial number m and stream desired value IshareDeviation:For θ (m)
(n) mathematic expectaion;CθForMaximum permissible value;
In t=KTs, the Centralized Controller of K=0,1,2,3 ... moment, parallel operation system starts to adopt by communication bus
Collect the output current Curr (m), m=1,2 of M online power modules, ┄, M;
Update the output current data of the online power module of serial number m:
Curr_store (m) (n)=Curr_store (m) (n+1), (11)
Curr_store (m) (T)=Curr (m), (12)
Wherein:M=1,2, ┄, M, n=1,2, ┄, T-1;
Calculate the online power module output current average value of serial number m:
Wherein:M=1,2,3 ... M;
Calculate parallel operation system load electric current Iout, meet:
Calculate online power module output current desired value Ishare, meet:
It calculates the online power module output current storage data Curr_store (m) (n) of serial number m and flows desired value
IshareDeviation:
θ (m) (n)=Curr_store (m) (n)-Ishare, (16)
Wherein:N=1 ... T;M=1,2,3 ... M;
Calculate the mathematic expectaion of the online power module deviation θ (m) (n) of serial number m;
Wherein:N=1 ... T;M=1,2,3 ... M;
Calculate online power module output current desired value IshareWith IrefAbsolute value of the difference Δ I, meet:
Δ I=| Ishare-Iref|, (18)
Judge whether Δ I meets inequality:
Δ I < σ, (19)
In the case where inequality (19) meet, judges the output current of the online power module of serial number m and flow mesh
Whether the mathematic expectaion of scale value deviation meets inequality:
If the online power module of serial number m meets inequality (20), illustrate that online power module performance is qualified;It is no
Then, online power module performance is unqualified, needs to cut qualified power module work from spare module.
Ungratified in inequality (19), calculating parallel operation system load electric current is IoutWhen, it efficiency and flows
Online power module quantity N under comprehensive performance optimal conditions*, meet:
It is ungratified in inequality (19), calculate the online power module regulated quantity Δ N of parallel operation system*, full
Foot:
ΔN*=N*- M, (22)
Centralized Controller increases (reduction) | Δ N*| a online power module, it is ensured that parallel operation system effectiveness and flow comprehensive
Close best performance.
The present invention provides the parallel operation system optimization controlling parties that based on efficiency and stream deviation it is expected weighted sum matrix
Method includes the following steps:
(1) the parallel operation system load electric current I of K power module composition is obtained in advanceoutFromAccording to interval
ForEquidistantly change toWhen (cover underloading, semi-load, specified load and overload conditions, U to meet and be necessary for
Positive integer not less than 20;INFor the rated current of power module), each power module is in different loads electric currentIn the case of acquire V output current Datacurr(m') (i) (j), output voltage Datavolt(m') (i) (j) and
Input power P (m') (i) (j), and calculate corresponding V power module efficiency eta (m') (i) (j) (V can be by user according to reality
Border determines size).Wherein:M' is power module serial number;I is the corresponding sequence number value of load current value;J is that output current acquires sequence
Number;M', i, j meet m'={ 1 ... K }, i={ 1 ... U }, j={ 1 ... V };
(2) it calculates power module output current and flows expectation electric current relative deviation absolute valueAnd corresponding mathematic expectaion(Em'iIt is smaller to show
The current sharing of power module can be better);Calculate the average value of power module efficiencyAnd efficiency etam'iWithWeighted suma∈[0,1](γm'iShow that the comprehensive performance of power module is better more greatly);With
γm'iMatrix A=(γ is built for elementm'i)K×U;M', i, j meet m'={ 1 ... K }, i={ 1 ... U }, j={ 1 ... V };
(3) solution matrix A=(γm'i)K×U1- normsCorresponding load current serial number a
Value, i.e.,Corresponding a values are set up, wherein:a∈[1,U],i∈[1,U];The a values are that power module is negative
Carrying electric current isWhen system efficiency and stream integrated performance index it is optimal;It enables efficiency and flows the optimal feelings of integrated performance index
Power module output current is I under conditionrefMeet
(4) with cycle TsFor the online power module quantity M of interval calculation parallel operation system, and to M online power supply moulds
The output current of block is acquired, and the output current data of the online power module of first serial number is labeled as Curr (1), when
Preceding online power module serial number m, enables m=1;
(5) the output current data array of the online power module of update serial number m, i.e.,:Curr_store (m) (n)=
Curr_store (m) (n+1), Curr_store (m) (T)=Curr (m);Wherein:N=1 ... T-1;M=1,2,3 ... M;T is
Positive integer more than 2;
(6) the output current average value of the online power module of serial number m is calculated:
Wherein:M=1,2,3 ... M;
(7) load current of the parallel operation system of M online power module compositions is calculatedWith
The equal current load electric current of line power module
(8) judge | Ishare-IrefWhether |≤σ is true?If it is, entering step (16);Conversely, then entering step
(9);
(9) it is reference current I to calculate online power module output currentrefWhen online power module quantity N*, i.e.,
(10) judge N*Whether≤1 is true?If it is, entering step (11);Conversely, entering step (12);
(11) N is set*=2;This is because N*<It is single supply module for power supply when 2, does not have flow equalizing function.
(12) online power module amount Δ N need to be adjusted by calculating parallel operation system*=N*-M;
(13) judge Δ N*Whether > 0 is true, if it is, entering step (14);Conversely, entering step (15);
(14) Centralized Controller increases Δ N*A online power module, subsequently into step (4);
(15) Centralized Controller reduces Δ N*A online power module, subsequently into step (4);
(16) it calculates the output current Curr_store (m) (n) of the online power module of serial number m and flows desired value
IshareDeviation θ (m) (n)=Curr_store (m) (n)-Ishare;Wherein:N=1 ... T;M=1,2,3 ... M;
(17) mathematic expectaion of the online power module deviation θ (m) (n) of serial number m is calculatedIts
In:N=1 ... N;M=1,2,3 ... M;
(18) m=1 is initialized;
(19) unqualified online power module quantity Num=0 is initialized;
(20) judge≤Cθ?(CθIt is expected for deviationMaximum permissible value) if so, entering step (23);Conversely, into
Enter step (21);
(21) the online power module current sharing that marking serial numbers are m can be undesirable;
(22) more new variables Num=Num+1;
(23) m=m+1 is updated;
(24) judge m<=M?If so, entering step (20);Otherwise, it enters step (25);
(25) the Num undesirable online power module of current sharing energy is offline, and start Num from stand-by power supply
A power module work;Subsequently into step (4).