Construction method of composite electrochemical polarization model of lithium ion battery pack
Technical Field
The invention relates to the technical field of charge state estimation of lithium ion battery packs, in particular to a construction method of a composite electrochemical polarization model of a lithium ion battery pack.
Background
The power battery system is one of the most important components, and has great influence on the performance and safety of the whole machine, and the lithium ion battery becomes the most widely used power battery for the machine at present by virtue of excellent performance; for a lithium ion Battery, not only a Battery core material with excellent performance needs to be continuously searched, but also a Battery Management System (BMS) needs to manage the Battery core material so as to prevent dangerous conditions such as spontaneous combustion and spontaneous explosion of the lithium ion Battery under certain extreme working conditions; the SOC (State of Charge) estimation of the battery is the basis of the BMS, the estimation precision is an important basis for judging the quality of a battery management system, the accurate SOC estimation is favorable for the full and reasonable use of the battery, the service life of the battery is prolonged, and the safety of the battery is improved; the accurate lithium ion battery model is the premise of improving the SOC estimation precision, so the accurate establishment of the battery model is the key of a battery management system; the battery is a nonlinear system, a common battery model can be divided into an electrochemical model and an equivalent circuit model, the electrochemical model is modeled according to an electrochemical reaction mechanism, the internal action of the battery is described through a mathematical relationship, and the internal action of the battery is described through pure mathematics to fully describe the characteristics of the battery, but the traditional pure electrochemical model is difficult to simulate the dynamic performance of the battery; the equivalent circuit model is established by forming a circuit network by using a resistor, a capacitor and a voltage source according to the dynamic characteristics and the working principle of the battery; however, for most equivalent circuit models, the Open Circuit Voltage (OCV) is directly measured or several estimated specified SOC values are measured in experiments, so that in the equivalent circuit model, simply using OCV as a function of SOC is not only time-consuming and labor-consuming, but is even prone to error; worse yet, the discrete representation of the OCV over multiple SOCs is not intuitive for kalman filtering or kalman filtering-based SOC estimation techniques.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: aiming at the problem that the SOC value cannot be accurately estimated due to insufficient accuracy of the existing lithium ion battery model, the construction method of the composite electrochemical polarization model of the lithium ion battery pack is provided.
The technical scheme adopted by the invention is as follows:
a construction method of a composite electrochemical polarization model of a lithium ion battery pack comprises the following steps:
s1, establishing an equivalent circuit model of the lithium ion battery pack;
s2, constructing a state space equation of the equivalent circuit model based on the kirchhoff circuit law;
s3, calculating the closed circuit voltage U of the lithium ion battery pack at the moment k according to the state space equation of the equivalent circuit modelL(k) The discrete expression of (1);
s4, solving the closed circuit voltage U of the lithium ion battery pack at the moment k by adopting a least square methodL(k) The parameters to be identified in the discrete expression are solved to calculate the parameters in the state space equation of the equivalent circuit model, so as to obtain the composite electrochemical polarization model of the lithium ion battery pack.
Further, the equivalent circuit model is composed of: polarization resistance RPAnd a polarization capacitor CPForming a first-order RC circuit; internal resistance to charging RcAnd internal discharge resistance RdRespectively reverse-serially connected with diodes DcAnd a diode DdThen connected in parallel, and then sequentially connected in series with a first-order RC circuit, an ohmic internal resistance Ro and an open-circuit voltage U of the lithium ion battery packOC(ii) a The closed circuit voltage of the equivalent circuit model is UL。
Further, in S2, based on kirchhoff' S circuit law, a state space equation for constructing the equivalent circuit model is:
in the formula, R0Ohmic internal resistance of the lithium ion battery pack; rPIs the polarization resistance of the lithium ion battery pack; cPIs the polarization capacitance of the lithium ion battery pack; rcdIs the internal resistance of the lithium ion battery pack during charging and discharging, Rcd=RcAt discharge time Rcd=Rd(ii) a k is the time of SOC estimation of the lithium ion battery pack; u shapeL(k) The closed circuit voltage of the lithium ion battery pack at the moment k is obtained; i isL(k) The output current of the lithium ion battery pack at the moment k is obtained; u shapeP(k) Is a polarization resistance RPThe voltage at the moment k at both ends; open circuit voltage U of lithium ion battery packOCEquivalence by Nernst model, K1、K2、K3For the three constants selected by the Nernst model, SOC (k) is the SOC estimate for the lithium ion battery pack at time k;
further, in step S3, the closed-circuit voltage U of the lithium ion battery pack at time k is calculated according to the state space equation of the equivalent circuit modelL(k) The discrete formula of (2) is as follows:
s3-1, polarization resistance R in the state space equationPVoltage U across the terminals at time kP(k) By first order backward difference instead, to obtain UP(k) The discretized equivalent differential equation of (1):
in the formula of UP(k) And UP(k-1) are respectively a polarization resistance RPVoltages of two ends at the time k and the time k-1, wherein T is a parameter detection period of the lithium ion battery pack, namely a signal sampling time interval;
s3-2, combined with equivalent UP(k) Is differential expression and UP(k) Obtaining the polarization resistance R by the discretization equivalent differential formulaPVoltage U across the terminals at time kP(k) The discrete equation of (c):
in the formula IL(k) And IL(k-1) respectively representing the output current of the lithium ion battery pack at the k moment and the k-1 moment;
s3-3, polarization resistance RPVoltage U across the terminals at time kP(k) Substituting the discrete expression into the state space equation to obtain the closed circuit voltage U of the lithium ion battery pack at the moment kL(k) The discrete equation of (c):
UL(k)=a1+a2UL(k-1)+a3ln(SOC(k))+a4ln(1-SOC(k))+a5IL(k)+a6IL(k-1)
in the formula, the parameter a to be identified1、a2、a3、a4、a5、a6The development of (a) is as follows:
further, in step S4, the closed-circuit voltage U of the lithium ion battery pack at time k is solved by using a least square methodL(k) The parameters to be identified in the discrete expression are solved to obtain the parameters to be identified, and the parameters in the state space equation of the equivalent circuit model are calculated by using the parameters to be identified, so that the process of obtaining the composite electrochemical polarization model of the lithium ion battery pack is as follows:
s4-1, the closed circuit voltage U of the lithium ion battery pack at the moment k is measuredL(k) The discrete expression of (a) is regarded as a six-element set consisting of independent variables and dependent variables, each element contains N units, and N is 1,2,3, …, k:
{UL(k)UL(k-1)SOC(k)SOC(k-1)IL(k)IL(k-1)}
wherein the matrix of dependent variables is: y ═ UL(1)UL(2)…UL(N)]T(ii) a The matrix of arguments is: x ═ β (1) β (2) … β (N)]T;
S4-2, according to UL(k) The discrete expression of (a) can be obtained:
β(k)=[1UL(k-1)ln(SOC(k))ln(1-SOC(k))IL(k)IL(k-1)]
then there is a parameter matrix to be identified: a ═ a1a2a3a4a5a6]T;
S4-3, estimating by using a least square method, and solving a parameter matrix A (X) needing to be identifiedTX)-1XTY=X-1Y, obtaining the parameter a to be identified1、a2、a3、a4、a5、a6;
S4-4, identifying the parameter a according to the requirement1、a2、a3、a4、a5、a6The expansion form of the equivalent circuit model is further deformed to obtain a parameter calculation formula in a state space equation of the equivalent circuit model:
solving the obtained parameters a needing to be identified1、a2、a3、a4、a5、a6Substituting the calculation result into the state space equation of the equivalent circuit model to obtain the composite electrochemical polarization model of the lithium ion battery pack.
In summary, due to the adoption of the technical scheme, the invention has the beneficial effects that:
according to the method, the representation accuracy and the calculation complexity of the lithium ion battery pack are comprehensively considered, the advantages of different equivalent circuit models are combined, and a lithium ion battery pack composite electrochemical polarization model is provided and constructed in a mode of combining electrochemistry and equivalent circuits; the provided composite electrochemical polarization model of the lithium ion battery pack realizes accurate mathematical expression of working conditions and working processes of the lithium ion battery pack by simulating different effects in the lithium ion battery pack in group cascade connection; a lithium ion battery parameter identification method based on a least square method of a lithium ion battery pack composite electrochemical polarization model is provided; provides a useful reference and reference for the combined application of the electrochemical model and the equivalent circuit model in the battery management.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained according to the drawings without inventive efforts.
Fig. 1 is a flow chart of a method for constructing a composite electrochemical polarization model of a lithium ion battery pack according to the present invention.
Fig. 2 is a schematic diagram of an equivalent circuit model of the lithium ion battery pack according to the present invention.
Detailed Description
The features and properties of the present invention are described in further detail below with reference to examples.
The method for constructing a lithium ion battery pack composite Electrochemical Polarization Model (S-EPM Model) provided in this embodiment, as shown in fig. 1, includes:
s1, establishing an equivalent circuit model of the lithium ion battery pack; as shown in fig. 2, the equivalent circuit model is composed of: polarization resistance RPAnd a polarization capacitor CPForming a first-order RC circuit; internal resistance to charging RcAnd internal discharge resistance RdRespectively reverse-serially connected with diodes DcAnd a diode DdThen connected in parallel, and then sequentially connected in series with a first-order RC circuit, an ohmic internal resistance Ro and an open-circuit voltage U of the lithium ion battery packOC(ii) a The closed circuit voltage of the equivalent circuit model isUL。
S2, constructing a state space equation of the equivalent circuit model based on the Kirchoff circuit law:
in the formula, R0The ohmic internal resistance of the lithium ion battery pack can represent the voltage drop of the two ends of the positive electrode and the negative electrode caused by the ohmic effect in the charging and discharging processes of the lithium ion battery pack; rPIs the polarization resistance, C, of a lithium ion battery packPIs the polarization capacitance, R, of a lithium ion battery packPAnd CPThe formed first-order RC circuit represents the relaxation effect in the charge and discharge process of the lithium ion battery pack, and further realizes the expression of the transient response of the lithium ion battery pack; rcdIs the internal resistance of the lithium ion battery pack during charging and discharging, Rcd=RcAt discharge time Rcd=Rd(ii) a k is the time of SOC estimation of the lithium ion battery pack; u shapeL(k) The closed circuit voltage of the lithium ion battery pack at the moment k is obtained; i isL(k) The output current of the lithium ion battery pack at the moment k is obtained; u shapeP(k) Is a polarization resistance RPThe voltage at the moment k at both ends; open circuit voltage U of lithium ion battery packOCThe Nernst model is used for equivalence, the Nernst model can provide a good fitting effect in the whole charging and discharging process of the lithium ion battery pack, and K is1、K2、K3The three constants selected by the Nernst model can be determined by different combination tests according to the special working conditions of the grouped work of the lithium ion battery pack; SOC (k) is the SOC estimated value of the lithium ion battery pack at the time k.
S3, calculating the closed circuit voltage U of the lithium ion battery pack at the moment k according to the state space equation of the equivalent circuit modelL(k) The discrete expression of (1);
s3-1, polarization resistance R in the state space equationPVoltage U across the terminals at time kP(k) By first order backward difference instead, to obtain UP(k) The discretized equivalent differential equation of (1):
in the formula of UP(k) And UP(k-1) are respectively a polarization resistance RPVoltages of two ends at the time k and the time k-1, wherein T is a parameter detection period of the lithium ion battery pack, namely a signal sampling time interval;
s3-2, combined with equivalent UP(k) Is differential expression and UP(k) Obtaining the polarization resistance R by the discretization equivalent differential formulaPVoltage U across the terminals at time kP(k) The discrete equation of (c):
in the formula IL(k) And IL(k-1) respectively representing the output current of the lithium ion battery pack at the k moment and the k-1 moment;
s3-3, polarization resistance RPVoltage U across the terminals at time kP(k) Substituting the discrete expression into the state space equation to obtain the closed circuit voltage U of the lithium ion battery pack at the moment kL(k) The discrete equation of (c):
UL(k)=a1+a2UL(k-1)+a3ln(SOC(k))+a4ln(1-SOC(k))+a5IL(k)+a6IL(k-1)
in the formula, the parameter a to be identified1、a2、a3、a4、a5、a6The development of (a) is as follows:
s4, solving the closed circuit voltage U of the lithium ion battery pack at the moment k by adopting a least square methodL(k) The parameters to be identified in the discrete expression are calculated by utilizing the parameters to be identified obtained by solvingObtaining parameters in a state space equation of the effective circuit model to obtain a lithium ion battery pack composite electrochemical polarization model:
s4-1, the closed circuit voltage U of the lithium ion battery pack at the moment k is measuredL(k) The discrete expression of (a) is regarded as a six-element set consisting of independent variables and dependent variables, each element contains N units, and N is 1,2,3, …, k:
{UL(k)UL(k-1)SOC(k)SOC(k-1)IL(k)IL(k-1)}
wherein the matrix of dependent variables is: y ═ UL(1)UL(2)…UL(N)]T(ii) a The matrix of arguments is: x ═ β (1) β (2) … β (N)]T;
S4-2, according to UL(k) The discrete expression of (a) can be obtained:
β(k)=[1UL(k-1)ln(SOC(k))ln(1-SOC(k))IL(k)IL(k-1)]
then there is a parameter matrix to be identified: a ═ a1a2a3a4a5a6]T;
S4-3, estimating by using a least square method, and solving a parameter matrix A (X) needing to be identifiedTX)-1XTY=X-1Y, obtaining the parameter a to be identified1、a2、a3、a4、a5、a6;
S4-4, identifying the parameter a according to the requirement1、a2、a3、a4、a5、a6The expansion form of the equivalent circuit model is further deformed to obtain a parameter calculation formula in a state space equation of the equivalent circuit model:
solving the obtained parameters a needing to be identified1、a2、a3、a4、a5、a6Parameter calculation substituted into the state space equation of the equivalent circuit modelAnd substituting the calculation result into a state space equation of the equivalent circuit model to obtain the composite electrochemical polarization model of the lithium ion battery pack.
In conclusion, the invention comprehensively considers the representation accuracy and the calculation complexity of the lithium ion battery pack, combines the advantages of different equivalent circuit models, and provides and constructs the composite electrochemical polarization model of the lithium ion battery pack by combining electrochemistry and equivalent circuits; the provided composite electrochemical polarization model of the lithium ion battery pack realizes accurate mathematical expression of working conditions and working processes of the lithium ion battery pack by simulating different effects in the lithium ion battery pack in group cascade connection, and experiments show that the error can be controlled within 0.5 percent, and the precision is higher; a lithium ion battery parameter identification method based on a least square method of a lithium ion battery pack composite electrochemical polarization model is provided; provides a useful reference and reference for the combined application of the electrochemical model and the equivalent circuit model in the battery management.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.