CN108509762A - A kind of the physicochemical change performance parameter analogy method and device of battery - Google Patents
A kind of the physicochemical change performance parameter analogy method and device of battery Download PDFInfo
- Publication number
- CN108509762A CN108509762A CN201810273896.9A CN201810273896A CN108509762A CN 108509762 A CN108509762 A CN 108509762A CN 201810273896 A CN201810273896 A CN 201810273896A CN 108509762 A CN108509762 A CN 108509762A
- Authority
- CN
- China
- Prior art keywords
- lithium ion
- battery
- coordinate system
- ion battery
- parameter
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/36—Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
- G01R31/3644—Constructional arrangements
- G01R31/3648—Constructional arrangements comprising digital calculation means, e.g. for performing an algorithm
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
Abstract
The present invention relates to the physicochemical change performance parameter analogy methods and device of a kind of battery.This method includes：Structure state equation corresponding with physicochemical change during battery use, wherein state equation include performance parameter corresponding with physicochemical change during battery use；Obtain actual measurement parameter corresponding with lithium ion battery；According to actual measurement parameter corresponding with lithium ion battery and state equation, performance parameter curve corresponding with physicochemical change during battery use is simulated.The performance parameter curve corresponding with physicochemical change during battery use obtained by the above method, can be very good the variation of the corresponding performance parameter in reaction cell discharge process.Moreover, the result is only to be combined part Experiment parameter with state equation and obtain, without being obtained by a large amount of experiment test.It is easy to operate, it is at low cost.
Description
Technical field
The present invention relates to battery technology process fields more particularly to a kind of physicochemical change performance parameter of battery to simulate
Method and device.
Background technology
As energy shortage has become society, economic development restraining factors, the application of lithium ion is increasingly extensive.Therefore,
How quickly and accurately to judge that the behaviour in service of battery, assessment, detection safety become the important of lithium ion battery functionization
Research topic.Wherein, battery mathematical model of structure on the basis of data accumulation and abstract operation principle refers to conductivity cell research and development
With the effective means of application.
Research and the most widely used model are parameter fitting models at present.Its main feature is that data through a large number of experiments
Accumulation and certain computer disposal, build the mapping of battery important performance characteristic (voltage, state-of-charge etc.) and pumping signal
Functional relation.Its major defect is to need to carry out a large amount of experiment test to obtain pumping signal and performance parameter relationship, behaviour
The result made acquired in cumbersome and experimental data is limited, once exceeding scope of experiment, then acquired result would be possible to not
Accurately.
Invention content
In order to solve the above technical problems, the present invention provides a kind of physicochemical change performance parameter analogy methods of battery
And device.
In a first aspect, the present invention provides a kind of physicochemical change performance parameter analogy method of battery, this method packet
It includes：
Corresponding with physicochemical change during the battery use state equation of structure, wherein state equation include and battery
The corresponding performance parameter of physicochemical change during use；
Obtain actual measurement parameter corresponding with lithium ion battery；
According to actual measurement parameter corresponding with lithium ion battery and state equation, simulation and physics during battery use
The curve of the corresponding performance parameter of chemical change.
The beneficial effects of the invention are as follows：It is obtained by the above method corresponding with physicochemical change during battery use
Performance parameter curve, can be very good the variation of the corresponding performance parameter in reaction cell discharge process.Moreover, should
As a result only it is that part Experiment parameter is combined with state equation and is obtained, without being obtained by a large amount of experiment test.
It is easy to operate, it is at low cost.
Second aspect, the present invention provides a kind of physicochemical change performance parameter simulator of battery, the device packets
It includes：
State equation construction unit, for building state equation corresponding with physicochemical change during battery use,
Wherein state equation includes performance parameter corresponding with physicochemical change during battery use；
Acquiring unit, for obtaining actual measurement parameter corresponding with lithium ion battery；
Processing unit, for according to actual measurement parameter corresponding with lithium ion battery and state equation, simulation to make with battery
With the curve of the corresponding performance parameter of physicochemical change in the process.
The beneficial effects of the invention are as follows：It is obtained by above-mentioned apparatus corresponding with physicochemical change during battery use
Performance parameter curve, can be very good the variation of the corresponding performance parameter in reaction cell discharge process.Moreover, should
As a result only it is that part Experiment parameter is combined with state equation and is obtained, without being obtained by a large amount of experiment test.
It is easy to operate, it is at low cost.
Description of the drawings
Fig. 1 is a kind of physicochemical change performance parameter analogy method flow signal of battery provided in an embodiment of the present invention
Figure；
Fig. 2 is cylindrical lithium ion battery Structural abstraction schematic diagram provided by the invention；
The quasi- two dimensional model of Fig. 3 lithium ion batteries provided by the invention builds schematic diagram；
Fig. 4 is certain the 14450 model cylindrical type ferrous phosphate lithium ion battery of company realized using the method for the present invention
Discharge curve simulates schematic diagram；
Fig. 5 is that certain company's 14450 model cylindrical type ferrous phosphate lithium ion battery electric discharge temperature rise curve simulates schematic diagram；
Fig. 6 is certain the 14450 model cylindrical type ferrous phosphate lithium ion battery of company realized using the method for the present invention
Lithium concentration simulates schematic diagram in discharge process electrolyte；
Fig. 7 is certain the 14450 model cylindrical type ferrous phosphate lithium ion battery of company realized using the method for the present invention
Lithium concentration simulates schematic diagram in discharge process anode；
Fig. 8 is a kind of physicochemical change performance parameter simulator structural representation of battery provided in an embodiment of the present invention
Figure.
Specific implementation mode
In being described below, for illustration and not for limitation, it is proposed that such as specific system structure, interface, technology it
The detail of class, to understand thoroughly the present invention.However, it will be clear to one skilled in the art that there is no these specific
The present invention can also be realized in the other embodiments of details.In other situations, omit to well-known device, circuit and
The detailed description of method, in case unnecessary details interferes description of the invention.
Fig. 1 is a kind of physicochemical change performance parameter analogy method flow signal of battery provided in an embodiment of the present invention
Figure.As shown in Figure 1, this method includes：
Step 110, state equation corresponding with physicochemical change during battery use is built.
Step 120, the actual measurement parameter corresponding with lithium ion battery of acquisition.
Specifically, actual measurement parameter may include experiment test parameter, fabrication process parameters and battery core characterisitic parameter etc..It is real
It can be the internal resistance numerical value under battery different temperatures, SOC states to test test parameter, and GITT methods may be used in these internal resistance numerical value
It is tested.Fabrication process parameters include mainly cell height, diameter, quality, membrane thicknesses, coating thickness etc..Core material
Characterisitic parameter includes mainly electrochemical reaction constant, lithium ion apparent diffusion coefficient etc..
Step 130, according to actual measurement parameter corresponding with lithium ion battery and state equation, simulation was used with battery
The curve of the corresponding performance parameter of physicochemical change in journey.
Optionally, in a specific example, such as physicochemical change includes at least during battery use：Solid phase
Mass transfer, mass transfer in liquid phase, heat transfer reaction, electrode polarization reaction and electrochemical reaction.
So, mathematical equation is then selection equation corresponding with these electrochemical reactions.For example, in structure and solid phase
It is uniform-spherical it is assumed that diffusion of the lithium ion in solid phase mass transfer based on active material when the corresponding equation of state of mass transfer
The corresponding mathematical formulae of Fick's second law, performance corresponding with physicochemical change during battery use can directly be used
Parameter and the state equation that structure solid phase mass transfer is combined with corresponding actual measurement parameter of lithium ion etc. obtained in advance.
And other state equations, it is built again by similar mode.It will specifically introduce in greater detail below.
Optionally, during structure is with battery use when the corresponding state equation of physicochemical change, from principle angle
For, corresponding performance parameter is as unknown number during physicochemical change during battery uses, with time, the change in space
Change is continuous and reversible, namely corresponding certain mapping relations.The essence of this mapping relations is each performance parameter in lithium
The characteristic being mutually coupled in ion battery discharge process carries out this mapping relations quantitative mathematics according to the principle of calculus
Description obtains function expression, the i.e. structure of completion status equation.
According to the real change occurred in lithium ion battery discharge process：The positive insertion and cathode deintercalation of lithium ion exist
Completed in solid positive and negative pole material, corresponding structure description lithium ion in positive and negative pole material concentration at any time, change in location
Partial differential equation；Lithium ion is from cathode to the moving process of anode, and concentration is at any time in the electrolytic solution for corresponding structure description lithium ion
Between, the partial differential equation of change in location；During the insertion of lithium ion and deintercalation, simultaneous and lithium ion equivalent electronics
Insertion and deintercalation, the process kinetics of this variation are described using Budm-Volmer equations；Companion while chemical change carries out
With the generation and transmission of heat, thus temperature at any time, the functional relation of change in location describes with varied heat quantity conductive equation.
Aforementioned each process constantly carries out, and potential difference is generated between the positive and negative anodes of lithium ion battery, can be carried in a certain range to outside
For stable voltage output, the i.e. functional equation of voltage output.
For the step angle (math equation establishes angle) specifically executed, it can specifically include：
Step 1101, two-dimensional coordinate system is radially established according to lithium ion battery 3 d structure model, determine lithium from
The parameter of each structural unit in the sub- battery 3 d structure model coordinate value in two-dimensional coordinate system respectively.
Specifically, lithium ion battery 3 d structure model it is practical be exactly lithium ion battery entity structure.Lithium ion battery
Overall structure is：Positive and negative electrode slurry is respectively coated on collector, and battery core is filled with electricity after being closely wound by anode, diaphragm, cathode
Solution liquid encapsulates.
It can be followed successively according to the structural unit of collector, cathode, diaphragm, anode, collector in the radial direction in battery
Circle permutation, as shown in Fig. 2.For a cylindrical battery, as shown in Figure 2 a, the constitutional repeating unit of Fig. 2 b can be abstracted as
It is wound.Its interior main matter composition is as shown in Figure 2 c, as minimum constitutional repeating unit.Each minimum repetitive structure
Unit is to contain collector, cathode, diaphragm, anode and collector.
And two-dimensional coordinate system is radially established according to lithium ion battery 3 d structure model, reality is exactly will be in battery
Portion's constitutional repeating unit is stretched by winding method, then using preset as coordinate origin, builds a two-dimensional coordinate system.
It is further alternative, after lithium ion battery 3 d structure model is radially unfolded, it is in turn divided into：Positive polar region
Domain, diaphragm and negative regions.As shown in figure 3, positive and negative active material in some constitutional repeating unit has actually been depicted in Fig. 3
The structural schematic diagram being separated from each other.Also it is both the structural schematic diagram of lithium-ion anode and cathode separation, as shown in Figure 3a.It will be in figure
Anode region is as positive pole zone Lp (left side), and cathode zone is as negative regions Ln (right side), and intermediate is then diaphragm Ls, such as
Shown in Fig. 3 b.
Meanwhile as above, it is contemplated that active material is assumed to be uniform-spherical, as shown in Figure 3c, so also needing to introducing one
A radial dimension r.As shown in Figure 3d.
So, it may include two kinds radially to establish two-dimensional coordinate system then according to lithium ion battery 3 d structure model
Mode：
The first, preset is as origin, using lithium ion battery radial direction as abscissa, time shaft using in positive pole zone
For ordinate, the first two-dimensional coordinate system is established；
Second, using lithium ion radial direction as abscissa, in lithium ion cell positive region or in negative regions
The diametric(al) perpendicular to lithium ion radial direction of active material is ordinate, establishes the second two-dimensional coordinate system.
It should be noted that positive pole zone and negative regions have active material, and all it is uniform orbicule.So
The radius of active material is equal, then can be made with positive pole zone or the diametric(al) of the active material of negative regions
For ordinate direction, the second two-dimensional coordinate system is established.
Determine coordinate value of the parameter in the first two-dimensional coordinate system in lithium ion battery 3 d structure model；Or determine lithium
Coordinate value of the parameter in the second two-dimensional coordinate system in ion battery 3 d structure model.Specifically, as shown in Figure 2 b, lithium from
Sub- battery 3 d structure model includes the structural unit of multiple repetitions, and the parameter included by each structural unit is shown in Fig. 2 c
's：Collector, cathode, diaphragm, anode and collector.When two-dimensional coordinate system is the first two-dimensional coordinate system, lithium-ion electric is determined
Coordinate value of the parameter in two-dimensional coordinate system in the 3 d structure model of pond, specifically includes：
Respectively determine each performance parameter in 1 constitutional repeating unit different location (plus plate current-collecting body, cathode, diaphragm,
Anode and negative current collector) coordinate value in the first two-dimensional coordinate system.
Alternatively, when two-dimensional coordinate system is the second two-dimensional coordinate system, the ginseng in lithium ion battery 3 d structure model is determined
Coordinate value of the number in two-dimensional coordinate system, specifically includes：It determines and is parallel to the second two-dimensional coordinate system cross in each active material
Coordinate diametrically any point and two-end-point (point number depend on Mathematical during finite difference set step-length)
Coordinate value in second two-dimensional coordinate system, it is practical to be also to determine the position of each active material and the radius of active material.
Step 1102, it according to coordinate value and math equation corresponding with physicochemical change during battery use, obtains
Take state equation corresponding with physicochemical change during battery use.
Radius r appeared in following state equation is the radius of active material described above, when t is corresponding
Between parameter t.
Specifically, in the corresponding state equation of structure solid phase mass transfer, concrete principle includes：
According to Fick's first law, the radial flow at any point is on the spherical surface at centre of sphere r,
Wherein, D is diffusion coefficient, and C is lithium concentration, and r is active material ball-type radius.
In for 2 spherical shells of dr, reaction particle concentration pace of change is
By 2 Taylor series expansions and high-order term is omitted, is obtained
Wherein, t is time, CS, kFor the extension concentration of lithium ion in solids, s is solid (solid) english abbreviation, k generations
Table p (positive- positive pole zones), n (negative regions negative-), s (separator- diaphragm areas).
During building mass transfer in liquid phase, relate generally to be diffusion, electromigration and electrification of the lithium ion in electrolyte
Learn reaction (insertion of electronics and deintercalation process).
It is the prior art that mass transfer in liquid phase, which builds state equation, and it is as follows only to enumerate final equation here：
Wherein, εE, kFor correction coefficient, CE, kFor the diffusion concentration of lithium ion in the electrolytic solution, jkIt is lithium ion in active material
Expect the insertion/deintercalation flow on surface；For the effective diffusion cofficient of lithium ion；F is Ferrari constant, and x is that the first two dimension is sat
Abscissa (the i.e. extension of plus plate current-collecting body-positive electrode-diaphragm-electrolyte-diaphragm-negative material-negative current collector in mark system
Direction), t+For the transport coefficient of lithium ion in the electrolytic solution.
The concrete principle of state equation when building electrochemical reaction process is：
The electrochemical reaction process of inside lithium ion cell is according to the classical electrochemical reaction equation side Butler-Volmer
Journey is described, and specific formula is as follows：
Wherein, JkFor the flow of lithium ion insertion/deintercalation；i0, kFor insertion/deintercalation current density, α is the ratio of electrode material
Surface area, R are english abbreviation (positive, negative), the η that ideal gas constant, n and p are respectively positive and negative anodesK, ECFor electricity
Chemically react overpotential, CeFor the concentration of lithium ion in the electrolytic solution, CS, surfIt is lithium ion in positive and negative pole material active material table
Bright concentration, CS, max0.5 power that the maximum concentration for being lithium ion in positive and negative anodes active material, 0.5 represent parameters.
Output voltage equation is：
VCell=U- ηEC-ηe-ηOhm(formula 7)
U=Up(θp)-Un(θn) (formula 8)
θk=Cs,k,surf/Cs,k,max(formula 9)
Wherein, VCellIt is battery equilibrium electromotive force, η for cell output voltage, UECFor electrochemical reaction overpotential, ηeFor liquid
Phase overpotential, ηOhmOhmic polarization overpotential, Up(θp) it is positive electrode electromotive force (for the letter of positive electrode material state-of-charge
Number), Un(θn) it is negative material electromotive force (for the function of negative electrode material state-of-charge), θkFor electrode material state-of-charge
(english abbreviation that k can be p or n is positive and negative anodes)
Structure liquid phase overpotential state equation concrete principle be：
Li+ can cause Gibbs free energy by the regions high concentration β phase low concentration α regional diffusions under the conditions of constant temperature, level pressure
Variation.For 1-1 type electrolyte LiFP6,t++t-=1, the activity factor takes 1, obtains
ΔGDFBe temperature for gibbs free energy change amount, T, z be this process be related to charge variation amount, Δ EDFFor
Potential change caused by diffusion, while this special relationship by the energy of reversible change, can obtain
The formula describes, and in liquid phase section, the variation of potential is caused by 2 parts, spreads (Δ EDF) and electrochemistry it is anti-
Answer (Δ ENernst).Entire electrode zone is further divided into n concentration thin layer, then Concentration overpotential can be expressed as n thin layer
Between liquid phase overpotential (ηe,Ce) it is cumulative,
It is respectively positive and negative afflux body interface at concentration layer starting and ending, therefore entire cell area liquid concentration is excessively electric
Gesture can approximate expression be
Wherein,WithRespectively represent the lithium ion liquid concentration at plus plate current-collecting body and negative current collector.
Heat transfer conditions equation concrete principle is：
Battery heat production is divided into hot 3 parts of electrochemical heat, ohm heat and Entropy Changes to be described respectively, further considers battery
Heat exchange with environment obtains.According to varied heat quantity conductive equation
Wherein, ρ vCpDensity, volume and the specific heat capacity of battery are indicated respectively.Consider that temperature sensing can not be deep into battery
Portion, analog result lacks apparent judgment criteria, therefore ignores the conduction of heat on inside battery and surface in model, that is, omits
▽(λk▽ T), only consider that temperature changes over time.Q is battery heat rate, and battery heat production is divided into electrochemical heat, ohm heat
It is described respectively with hot 3 parts of Entropy Changes, further considers that the heat exchange of battery and environment obtains,
Wherein U=Up-UnFor battery open circuit voltage, dUk/ dT is the temperature coefficient of electrode material equilibrium potential, can be by Entropy Changes
Experiment obtains.And the two is fitted to the empirical function of SOC, can directly be read use.Wherein h batteries are with environment heat exchange
Number, A battery surfaces product, TambFor environment temperature, V be cell output voltage, I is electric current.
By above-mentioned formula as can be seen that the simple mathematical description structure model of state equation that the present invention is established, phase
Compared with the critical data description structure mould of physicochemical change in certain courses of work by lithium ion battery in the prior art
For type (such as principle model in the prior art, its shortcoming is that passing through physical chemistry in the lithium ion battery course of work
The strict mathematical description structure model of variation, evens penetrate into quantum rank, therefore professional too strong, is that it imports practical application
Maximum bottleneck), simpler, professional low, calculation amount is small, more extensive in practical applications.
And when being solved and simulated especially by above-mentioned state equation, it can be realized using MATLAB processors with this.Mould
The mathematical form of type equations and saturation equation is analyzed, partial differential equation of second order involved in calculating (solid, liquid phase concentration governing equation), often
The differential equation (energy-balance equation), common algebra equation (electrochemical reaction equation and output voltage equation).For not Tongfang
The mathematical characteristics of journey use different computational methods, for partial differential equation using the discrete rear programming iterative solutions of FDM, ordinary differential
Equation carries out Analytical Solution using classical quadravalence R-K methods iterative solution, algebraic equation by the symbolic operation function of MATLAB.
It is integrated finally by MATLAB processors.Specifically workflow is realized using MATLAB, it can be understood as the prior art.This
In only do simple introduction, it is specific as follows without doing excessive description：
Read the parameters such as multiplying power (discharge current), parameter to be solved and input temp.It is calculated according to electrochemical reaction equation
Lithium ion flow (hole wall flow), the variation of solid phase and liquid concentration and overpotential caused by then further calculating thus, this
The variation of a little parameters represents the progress of electrochemical reaction and the generation of fuel factor in battery, therefore can export new temperature.
One calculates cycle and terminates, if reaching battery system blanking voltage, the end that exits the program calculates.If not up to ending electricity
Pressure, the then calculating parameter recycled upper one are inputted as the calculating of subsequent cycle.When reaching blanking voltage, calculation process knot
Beam.It is final to obtain performance corresponding with physicochemical change during battery use, and simulate and physics during battery use
The corresponding performance parameter curve of chemical change.Here performance parameter include temperature of the battery in discharge process, voltage and
The variation etc. of lithium concentration.
Fig. 4 is certain the 14450 model cylindrical type ferrous phosphate lithium ion battery of company realized using the method for the present invention
Discharge curve simulates schematic diagram.Under the conditions of being 30 DEG C for environment temperature, the discharge curve of different multiplying, dotted line be analog result,
Curve is measured curve.It can be seen that having preferable accordance.Initial rapid decrease, plateau and corner feature obtain
It presents.
Fig. 5 is certain the 14450 model cylindrical type ferrous phosphate lithium ion battery of company realized using the method for the present invention
The temperature rise curve that discharges simulates schematic diagram.As shown in figure 5, under the conditions of environment temperature is 30 DEG C, the discharge battery surface of different multiplying
Temperature variations.Dotted line is analog result, and curve is measured curve.Preferably meet it can be seen that simulation has with actual measurement
Property presents the characteristic that temperature in measured data is begun to ramp up, middle stable, later stage rise rapidly well.
Fig. 6 is certain the 14450 model cylindrical type ferrous phosphate lithium ion battery of company realized using the method for the present invention
Lithium concentration simulates schematic diagram in discharge process electrolyte.The variation of lithium concentration can not be measured directly in electrolyte,
Simulation result is presented in this.Output at 4 interfaces (positive and negative electrode afflux body interface and positive and negative electrode diaphragm interface) liquid phase Li+ at any time
Between variation (Fig. 6 left).It therefrom it can be seen that, is stepped up from plus plate current-collecting body to negative current collector liquid concentration, and big
After about 500s, the processes such as electrochemistry, electromigration and diffusion reach stable state, and concentration is kept constant.Therefore, it is outputed within 500s
The distribution situation of liquid phase Li+ concentration in entire electrode zone is partly diaphragm region as shown in Fig. 6 (Fig. 6 is right), wherein between dotted line
Domain, due to not having electrochemical reaction, concentration is by spreading control thus being linearly distributed.It can be seen that the result that simulation obtains meets theory
On changing rule.
Fig. 7 is certain the 14450 model cylindrical type ferrous phosphate lithium ion battery of company realized using the method for the present invention
Lithium concentration simulates schematic diagram in discharge process anode.Equally, lithium concentration also can not be measurement originally in solid phase, emulation knot
Fruit is following as shown in 7 left sides, is different moments, the situation of change of 3 representative locations solid concentrations of positive active particles, electric discharge is just
Phase, due to electrochemical reaction so that particle surface Li+ concentration increases rapidly, and due to needing by solid-state diffusion mistake at the centre of sphere
Journey, so rising slow.In discharge process, Li+ is embedded in positive electrode, and solid-state diffusion process ecto-entad carries out, therefore
Grain surface concentration always above the centre of sphere at, until electric discharge terminate concentration just reach unanimity.Fig. 7 (right figure) then illustrates an activity
The distribution of Li+ inside particle is the centre of sphere at interior joint 20, therefrom it can also be seen that with electric discharge progress, inside and outside particle
Concentration difference gradually reduce, surface concentration change first quick and back slow.
The physicochemical change performance parameter analogy method of a kind of battery provided in an embodiment of the present invention, by acquisition and electricity
The corresponding performance parameter curve of physicochemical change, can be very good in reaction cell discharge process therewith during the use of pond
The variation of corresponding performance parameter.Moreover, the result is only to be combined part Experiment parameter with state equation and obtain, without
It needs to obtain by a large amount of experiment test.It is easy to operate, it is at low cost.
Correspondingly, the embodiment of the present invention additionally provides a kind of physicochemical change performance parameter simulator structure of battery
Schematic diagram.Specifically as shown in figure 8, the device includes：State equation construction unit 801, acquiring unit 802 and processing unit
803。
Wherein, state equation construction unit 801, for building shape corresponding with physicochemical change during battery use
State equation, wherein state equation include performance parameter corresponding with physicochemical change during battery use.
Acquiring unit 802, for obtaining actual measurement parameter corresponding with lithium ion battery.
Processing unit 803, for according to actual measurement parameter corresponding with lithium ion battery and state equation, simulating and electricity
The curve of the corresponding performance parameter of physicochemical change during the use of pond.
Optionally, state equation construction unit 801 is specifically used for：
Two-dimensional coordinate system is radially established according to lithium ion battery 3 d structure model, determines lithium ion battery three-dimensional
Coordinate value of the parameter in two-dimensional coordinate system in structural model；
According to coordinate value and math equation corresponding with physicochemical change during battery use, acquisition makes with battery
With the corresponding state equation of physicochemical change in the process.
It is further alternative, after lithium ion battery 3 d structure model is radially unfolded, it is in turn divided into：Positive polar region
Domain, diaphragm and negative regions, positive pole zone and negative regions include the active material of unifonn spherical；
State equation construction unit 801 is specifically used for：
Preset is as origin using in positive pole zone, and using lithium ion battery radial direction as abscissa, time shaft is vertical sits
Mark, establishes the first two-dimensional coordinate system；
And using lithium ion radial direction as abscissa, with the work in lithium ion cell positive region or in negative regions
Property material the diametric(al) perpendicular to lithium ion radial direction be ordinate, establish the second two-dimensional coordinate system.
Further alternative, physicochemical change includes at least during battery use：Solid phase mass transfer, mass transfer in liquid phase, biography
Thermal response, electrode polarization reaction and electrochemical reaction.
Further alternative, lithium ion battery 3 d structure model includes the structural unit of multiple repetitions, each structure list
Parameter included by member is：Collector, cathode, diaphragm, anode and collector；When two-dimensional coordinate system is the first two-dimensional coordinate system
When, state equation construction unit 801 is specifically used for：Determine the parameter in the lithium ion battery 3 d structure model described
Coordinate value in one two-dimensional coordinate system；Or determine that the parameter in the lithium ion battery 3 d structure model is two-dimentional described second
Coordinate value in coordinate system.Specifically, including determining collector, cathode, diaphragm, anode and collector in each structural unit
Coordinate value in the first two-dimensional coordinate system respectively.
Further alternative, when two-dimensional coordinate system is the second two-dimensional coordinate system, state equation construction unit 801 is specifically used
In：
Determine the diameter two-end-point for being parallel to the second two-dimensional coordinate system abscissa in each active material respectively second
Coordinate value in two-dimensional coordinate system.
The function performed by each component in the device physics during a kind of battery use of above-described embodiment
It is had been described in detail in chemical change performance parameter analogy method, which is not described herein again.
The physicochemical change performance parameter simulator of a kind of battery provided in an embodiment of the present invention, by acquisition and electricity
The corresponding performance parameter curve of physicochemical change, can be very good in reaction cell discharge process therewith during the use of pond
The variation of corresponding performance parameter.Moreover, the result is only to be combined part Experiment parameter with state equation and obtain, without
It needs to obtain by a large amount of experiment test.It is easy to operate, it is at low cost.
Reader should be understood that in the description of this specification reference term " one embodiment ", " is shown " some embodiments "
The description of example ", " specific example " or " some examples " etc. mean specific features described in conjunction with this embodiment or example, structure,
Material or feature are included at least one embodiment or example of the invention.In the present specification, above-mentioned term is shown
The statement of meaning property need not be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described
It may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other, this
The technical staff in field can be by the spy of different embodiments or examples described in this specification and different embodiments or examples
Sign is combined.
Although the embodiments of the present invention has been shown and described above, it is to be understood that above-described embodiment is example
Property, it is not considered as limiting the invention, those skilled in the art within the scope of the invention can be to above-mentioned
Embodiment is changed, changes, replacing and modification.
Claims (10)
1. a kind of physicochemical change performance parameter analogy method of battery, which is characterized in that the method includes：
Structure state equation corresponding with physicochemical change during lithium ion battery use, wherein the state equation includes
Performance parameter corresponding with physicochemical change during the battery use；
Obtain actual measurement parameter corresponding with lithium ion battery；
According to the actual measurement parameter corresponding with lithium ion battery and the state equation, simulation was used with the battery
The curve of the corresponding performance parameter of physicochemical change in journey.
2. according to the method described in claim 1, it is characterized in that, the structure and physicochemical change during battery use
Corresponding state equation, specifically includes：
Two-dimensional coordinate system is radially established according to the lithium ion battery 3 d structure model, and determines the lithium-ion electric
Coordinate value of the parameter in the two-dimensional coordinate system in the 3 d structure model of pond；
According to the coordinate value and math equation corresponding with physicochemical change during the battery use, described in acquisition
State equation.
3. according to the method described in claim 2, it is characterized in that, the lithium ion battery 3 d structure model radially
After expansion, it is in turn divided into：Positive pole zone, diaphragm and negative regions, the positive pole zone and the negative regions include equal
Even spherical active material；
It is described that two-dimensional coordinate system is radially established according to the lithium ion battery 3 d structure model, and determine the lithium from
Coordinate value of the parameter in the two-dimensional coordinate system in sub- battery 3 d structure model, specifically includes：
Using preset in the positive pole zone as origin, using the lithium ion battery radial direction as abscissa, time shaft is
Ordinate establishes the first two-dimensional coordinate system；
And using the lithium ion radial direction as abscissa, in the lithium ion cell positive region or cathode area
The diametric(al) perpendicular to the lithium ion radial direction of active material in domain is ordinate, establishes the second two-dimensional coordinate
System；
Determine coordinate value of the parameter in first two-dimensional coordinate system in the lithium ion battery 3 d structure model；Or really
Coordinate value of the parameter in second two-dimensional coordinate system in the fixed lithium ion battery 3 d structure model.
4. according to the method described in claim 1-3, which is characterized in that physicochemical change is at least wrapped during battery use
It includes：Solid phase mass transfer, mass transfer in liquid phase, heat transfer reaction, electrode polarization reaction and electrochemical reaction.
5. according to the method described in claim 3, it is characterized in that, the lithium ion battery 3 d structure model includes multiple heavy
Multiple structural unit, the parameter included by each structural unit are：Collector, cathode, diaphragm, anode and collector；It is described true
Coordinate value of the parameter in first two-dimensional coordinate system in the fixed lithium ion battery 3 d structure model, specifically includes：
Determine that collector, cathode, diaphragm, anode and collector in each structural unit are sat in first two dimension respectively
Coordinate value in mark system.
6. according to the method described in claim 3, it is characterized in that, in the determination lithium ion battery 3 d structure model
Coordinate value of the parameter in second two-dimensional coordinate system, specifically include：
Determine the diameter two-end-point for being parallel to the second two-dimensional coordinate system abscissa in each active material respectively described
Coordinate value in second two-dimensional coordinate system.
7. a kind of physicochemical change performance parameter simulator of battery, which is characterized in that described device includes：
State equation construction unit, for building state equation corresponding with physicochemical change during battery use, wherein
The state equation includes performance parameter corresponding with physicochemical change during the battery use；
Acquiring unit, for obtaining actual measurement parameter corresponding with lithium ion battery；
Processing unit, for actual measurement parameter corresponding with lithium ion battery and the state equation according to, simulation and institute
State the curve of the corresponding performance parameter of physicochemical change during battery use.
8. device according to claim 7, which is characterized in that the state equation construction unit is specifically used for：
Two-dimensional coordinate system is radially established according to the lithium ion battery 3 d structure model, and determines the lithium-ion electric
Coordinate value of the parameter in the two-dimensional coordinate system in the 3 d structure model of pond；
According to the coordinate value and math equation corresponding with physicochemical change during the battery use, described in acquisition
State equation.
9. device according to claim 8, which is characterized in that the lithium ion battery 3 d structure model is radially
After expansion, it is in turn divided into：Positive pole zone, diaphragm and negative regions, the positive pole zone and the negative regions include equal
Even spherical active material；
State equation construction unit is specifically used for：
Using preset in the positive pole zone as origin, using the lithium ion battery radial direction as abscissa, time shaft is
Ordinate establishes the first two-dimensional coordinate system；
And using the lithium ion radial direction as abscissa, in the lithium ion cell positive region or cathode area
The diametric(al) perpendicular to the lithium ion radial direction of active material in domain is ordinate, establishes the second two-dimensional coordinate
System；
Determine coordinate value of the parameter in first two-dimensional coordinate system in the lithium ion battery 3 d structure model；Or really
Coordinate value of the parameter in second two-dimensional coordinate system in the fixed lithium ion battery 3 d structure model.
10. according to claim 7-9 any one of them devices, which is characterized in that physicochemical change during battery use
It includes at least：Solid phase mass transfer, mass transfer in liquid phase, heat transfer reaction, electrode polarization reaction and electrochemical reaction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810273896.9A CN108509762A (en) | 2018-03-29 | 2018-03-29 | A kind of the physicochemical change performance parameter analogy method and device of battery |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810273896.9A CN108509762A (en) | 2018-03-29 | 2018-03-29 | A kind of the physicochemical change performance parameter analogy method and device of battery |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108509762A true CN108509762A (en) | 2018-09-07 |
Family
ID=63379497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810273896.9A Pending CN108509762A (en) | 2018-03-29 | 2018-03-29 | A kind of the physicochemical change performance parameter analogy method and device of battery |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108509762A (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109388074A (en) * | 2018-09-20 | 2019-02-26 | 深圳大学 | A kind of battery pack emulation mode and system, BMS calibration system |
| CN109738806A (en) * | 2018-12-29 | 2019-05-10 | 蜂巢能源科技有限公司 | The method, apparatus of simulated battery heat production rate, medium |
| CN110967618A (en) * | 2019-03-21 | 2020-04-07 | 宁德时代新能源科技股份有限公司 | Circuit parameter detection method and detection device |
| CN111044922A (en) * | 2019-12-31 | 2020-04-21 | 国联汽车动力电池研究院有限责任公司 | Method for testing heat dissipation uniformity of lithium ion battery in environmental chamber |
| CN111766518A (en) * | 2020-06-18 | 2020-10-13 | 合肥国轩高科动力能源有限公司 | Quantitative determination method for reversible lithium separation of lithium ion battery |
| CN112526357A (en) * | 2020-11-25 | 2021-03-19 | 上海空间电源研究所 | Lithium ion battery power matching performance evaluation method |
| CN112964999A (en) * | 2021-03-18 | 2021-06-15 | 潍柴动力股份有限公司 | Battery state of charge acquisition method, device, equipment, medium and program product |
| CN113253131A (en) * | 2021-05-17 | 2021-08-13 | 蜂巢能源科技有限公司 | Method and device for determining charging and discharging performance of battery core, storage medium and electronic equipment |
| CN113609689A (en) * | 2021-08-09 | 2021-11-05 | 中国工程物理研究院激光聚变研究中心 | Simulation method, quality prediction method and system for fusion target pellet in complex cavity |
| CN116027199A (en) * | 2022-12-08 | 2023-04-28 | 帕诺(常熟)新能源科技有限公司 | Method for detecting short circuit in whole service life of battery cell based on electrochemical model parameter identification |
| CN116029143A (en) * | 2023-02-07 | 2023-04-28 | 上海玫克生储能科技有限公司 | Early warning diagnosis method, system, electronic equipment and medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102661983A (en) * | 2012-04-27 | 2012-09-12 | 苏州市汉微环保科技有限公司 | Value iterative fitting method of electrochemical impedance spectrums of porous electrodes |
| CN103020445A (en) * | 2012-12-10 | 2013-04-03 | 西南交通大学 | SOC (State of Charge) and SOH (State of Health) prediction method of electric vehicle-mounted lithium iron phosphate battery |
| CN104899439A (en) * | 2015-06-02 | 2015-09-09 | 吉林大学 | Mechanism modeling method for lithium ion battery |
| CN104991980A (en) * | 2014-10-31 | 2015-10-21 | 吉林大学 | Electrochemical mechanism modeling method for lithium-ion battery |
-
2018
- 2018-03-29 CN CN201810273896.9A patent/CN108509762A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102661983A (en) * | 2012-04-27 | 2012-09-12 | 苏州市汉微环保科技有限公司 | Value iterative fitting method of electrochemical impedance spectrums of porous electrodes |
| CN103020445A (en) * | 2012-12-10 | 2013-04-03 | 西南交通大学 | SOC (State of Charge) and SOH (State of Health) prediction method of electric vehicle-mounted lithium iron phosphate battery |
| CN104991980A (en) * | 2014-10-31 | 2015-10-21 | 吉林大学 | Electrochemical mechanism modeling method for lithium-ion battery |
| CN104899439A (en) * | 2015-06-02 | 2015-09-09 | 吉林大学 | Mechanism modeling method for lithium ion battery |
Non-Patent Citations (4)
| Title |
|---|
| 庞辉: "基于电化学模型的锂离子电池多尺度建模及其简化方法", 《物理学报》 * |
| 李冠宇: "小功率 LiFePO4锂离子电池简化电化学-热耦合模型研究", 《中国优秀硕士学位论文全文数据库（电子期刊） 工程科技II辑》 * |
| 李冠宇: "磷酸亚铁锂-石墨体系恒流放电的理论模拟", 《中国固态离子学暨电池材料青年学术论坛-论文摘要集》 * |
| 王铭 等: "锂离子电池模型研究进展", 《电源技术》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109388074A (en) * | 2018-09-20 | 2019-02-26 | 深圳大学 | A kind of battery pack emulation mode and system, BMS calibration system |
| CN109738806A (en) * | 2018-12-29 | 2019-05-10 | 蜂巢能源科技有限公司 | The method, apparatus of simulated battery heat production rate, medium |
| CN110967618A (en) * | 2019-03-21 | 2020-04-07 | 宁德时代新能源科技股份有限公司 | Circuit parameter detection method and detection device |
| CN111044922A (en) * | 2019-12-31 | 2020-04-21 | 国联汽车动力电池研究院有限责任公司 | Method for testing heat dissipation uniformity of lithium ion battery in environmental chamber |
| CN111766518A (en) * | 2020-06-18 | 2020-10-13 | 合肥国轩高科动力能源有限公司 | Quantitative determination method for reversible lithium separation of lithium ion battery |
| CN112526357B (en) * | 2020-11-25 | 2023-04-18 | 上海空间电源研究所 | Lithium ion battery power matching performance evaluation method |
| CN112526357A (en) * | 2020-11-25 | 2021-03-19 | 上海空间电源研究所 | Lithium ion battery power matching performance evaluation method |
| CN112964999A (en) * | 2021-03-18 | 2021-06-15 | 潍柴动力股份有限公司 | Battery state of charge acquisition method, device, equipment, medium and program product |
| CN112964999B (en) * | 2021-03-18 | 2022-10-25 | 潍柴动力股份有限公司 | Battery state of charge acquisition method, device, equipment, medium and program product |
| CN113253131A (en) * | 2021-05-17 | 2021-08-13 | 蜂巢能源科技有限公司 | Method and device for determining charging and discharging performance of battery core, storage medium and electronic equipment |
| CN113253131B (en) * | 2021-05-17 | 2022-08-30 | 蜂巢能源科技股份有限公司 | Method and device for determining charging and discharging performance of battery core, storage medium and electronic equipment |
| CN113609689A (en) * | 2021-08-09 | 2021-11-05 | 中国工程物理研究院激光聚变研究中心 | Simulation method, quality prediction method and system for fusion target pellet in complex cavity |
| CN113609689B (en) * | 2021-08-09 | 2023-07-04 | 中国工程物理研究院激光聚变研究中心 | Simulation method, quality prediction method and system of fusion target pill in complex cavity |
| CN116027199A (en) * | 2022-12-08 | 2023-04-28 | 帕诺(常熟)新能源科技有限公司 | Method for detecting short circuit in whole service life of battery cell based on electrochemical model parameter identification |
| CN116027199B (en) * | 2022-12-08 | 2023-09-29 | 帕诺(常熟)新能源科技有限公司 | Method for detecting short circuit in whole service life of battery cell based on electrochemical model parameter identification |
| CN116029143A (en) * | 2023-02-07 | 2023-04-28 | 上海玫克生储能科技有限公司 | Early warning diagnosis method, system, electronic equipment and medium |
| CN116029143B (en) * | 2023-02-07 | 2023-09-19 | 上海玫克生储能科技有限公司 | Early warning diagnosis method, system, electronic equipment and medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108509762A (en) | A kind of the physicochemical change performance parameter analogy method and device of battery | |
| CN107066722B (en) | Electrochemical model-based combined estimation method for state of charge and state of health of power battery system | |
| Kong et al. | Pseudo-two-dimensional model and impedance diagnosis of micro internal short circuit in lithium-ion cells | |
| Fu et al. | Modeling, validation and analysis of mechanical stress generation and dimension changes of a pouch type high power Li-ion battery | |
| Jiang et al. | A microscopic investigation of ion and electron transport in lithium-ion battery porous electrodes using the lattice Boltzmann method | |
| Tröltzsch et al. | Characterizing aging effects of lithium ion batteries by impedance spectroscopy | |
| CN110457742A (en) | A kind of lithium ion battery battery chemo-thermo-mechanical coupling model and modeling method | |
| US10539621B2 (en) | Method and apparatus for identifying a battery model | |
| Rajabloo et al. | An inverse method for estimating the electrochemical parameters of lithium-ion batteries | |
| Ender | An extended homogenized porous electrode model for lithium-ion cell electrodes | |
| CN110165314A (en) | Battery battery core performance parameter acquisition methods and acquisition device | |
| CN109738806B (en) | Method, device and medium for simulating heat generation rate of battery | |
| CN109946622B (en) | Lithium deposition prediction method and device for lithium ion battery | |
| CN109581064A (en) | Measure method, apparatus, the medium of battery low temperature direct resistance | |
| CN113253131B (en) | Method and device for determining charging and discharging performance of battery core, storage medium and electronic equipment | |
| CN113702845A (en) | Method and equipment for evaluating core parameters of retired lithium battery | |
| Chen et al. | Electrochemical model parameter identification of lithium-ion battery with temperature and current dependence | |
| CN111999665B (en) | Lithium ion battery aging test method based on micro-mechanism automobile driving condition | |
| Chayambuka et al. | Determination of state-of-charge dependent diffusion coefficients and kinetic rate constants of phase changing electrode materials using physics-based models | |
| Yang et al. | Effects of temperature differences among cells on the discharging characteristics of lithium‐ion battery packs with series/parallel configurations during constant power discharge | |
| CN114280493B (en) | Battery internal health state diagnosis method and system based on simplified P2D model | |
| CN109033619A (en) | A kind of transient temperature model modelling approach of 18650 type lithium battery discharge cycles | |
| Esfahanian et al. | Design parameter study on the performance of lead-acid batteries | |
| CN114996932A (en) | Modeling method of lithium ion battery overcharge thermal runaway model | |
| Wang et al. | Construction of electrochemical model for high C-rate conditions in lithium-ion battery based on experimental analogy method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180907 |
|
| RJ01 | Rejection of invention patent application after publication |