CN107145629A - A kind of method for optimizing battery electrode thickness - Google Patents
A kind of method for optimizing battery electrode thickness Download PDFInfo
- Publication number
- CN107145629A CN107145629A CN201710209193.5A CN201710209193A CN107145629A CN 107145629 A CN107145629 A CN 107145629A CN 201710209193 A CN201710209193 A CN 201710209193A CN 107145629 A CN107145629 A CN 107145629A
- Authority
- CN
- China
- Prior art keywords
- model
- electrode
- battery
- thickness
- electrochemical
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Battery Electrode And Active Subsutance (AREA)
- Secondary Cells (AREA)
Abstract
A kind of method for optimizing battery electrode thickness, comprises the following steps：1) maximized using energy density E or power density P is maximized and is used as the optimization aim for this method；2) electrode structural designs parameter, the kinetic parameter and thermal physical property parameter of electrode material of base batteries specification are obtained；3) battery electrochemical thermal coupling model is set up, the electrochemical heat coupling model is the coupling of accurate a two-dimentional electrochemical model and a three-dimensional thermal model；4) validity of model is verified；5) the battery electrode thickness after being optimized.The present invention can effectively shorten the construction cycle of new material or new product, reduce development cost, have certain directive significance for the exploitation of new material or new product.
Description
Technical field
The present invention relates to field of lithium, more particularly to a kind of method for optimizing battery electrode thickness.
Background technology
At present, continuing to develop with lithium ion battery technology, lithium ion battery applications market is constantly widened, correspondingly objective
Family requires that emphasis is also more and more diversified to the specifications and characteristics of lithium ion battery.
Battery new product development is carried out using traditional means of experiment and its optimization design consumes a large amount of manpower and materials, and is imitated
Rate is relatively low.The optimization design for carrying out battery using computer simulation means can greatly shorten the battery design cycle.Moreover, at present
Analogue simulation has become the important means that research inside battery characteristic instructs battery design.Electrochemical heat coupling model is also increasingly
Maturation, such as Jie Li (J Power Source, 2014, DOI:10.1016/j.jpowsour.2013.01.007) set up ginseng
The electrochemical heat coupling model of number dynamic responses, can battery discharge behavior preferably under simulation different multiplying.
Calculating optimal value obtains optimized algorithm and emerged in an endless stream at present, but application characteristic and computational efficiency are each has something to recommend him, wherein
BOBYQA methods (quasi- two-dimensional linear border optimization) are a kind of New Type of Numerical optimized algorithm that Powell was proposed in 2009, algorithm
Advantage be that need not solve the derivative of object function, can conveniently processing black box problem.
The B of patent CN 102170022 disclose a kind of design method of lithium ion battery, the reality for calculating lithium battery
Performance parameter, that is, set up it is constructed under the conditions of first shape lithium battery and the second shape lithium battery between performance data mould
Whether type, meet according to the actual performance parameter of another shaped cells of battery actual performance parameter prediction of one of which shape
It is required that.The patent is mainly used in after cell shapes change, the prediction of battery performance, but is not suitable for the different rule of same shape
The battery optimization design that lattice or new material are introduced.
The content of the invention
The technical problem to be solved in the present invention overcomes the deficiencies in the prior art, and there is provided a kind of optimization battery electrode thickness
Method, to shorten lithium ion battery design cycle, reduction development cost.
In order to solve the above technical problems, technical scheme proposed by the present invention is：A kind of method for optimizing battery electrode thickness,
Comprise the following steps：1) maximized using energy density E or power density P is maximized as the optimization aim for this method, its
The middle maximized object function of energy density isThe maximized object functions of power density P are：
Wherein, V, M represent discharge voltage and electrode quality respectively, and i is discharge current, and t is discharge time；According to client's
Demand, selects battery specifications close therewith as the base batteries specification of this method, battery electrode is carried out on this basis thick
The optimization of degree；
2) electrode structural designs parameter, the kinetic parameter of electrode material and the hot physical property ginseng of base batteries specification are obtained
Number；
3) battery electrochemical thermal coupling model is set up, the electrochemical heat coupling model includes an accurate two-dimentional electrochemistry mould
Type and a Three Dimensional Thermal coupling model；
1. accurate two-dimentional electrochemical model：The different parts of electrode are represented with line segment, one-dimensional geometrical model are set up, one-dimensional
The differential or partial differential equation of loading description discharge process, net is carried out by finite element theory to geometrical model on geometrical model
After formatting, the differential or partial differential equation of description discharge process are calculated, obtains different in discharge process under a thickness of electrode
Discharge time corresponding discharge voltage；
The one-dimensional geometrical model includes herein below：
A, lithium ion is described in solid phase internal transmission on negative pole and positive pole with Fick's first law：
Boundary condition equation is：
B, describes liquid phase lithium ion course of dissolution on positive pole, negative pole and barrier film：
T in formula+、F、jlo4,iLithium ion carry-over factor, Faraday constant and local current densities are represented respectively；
C：Electrochemical reaction is described with Butler-Volmer equations on positive pole and negative pole
αa,i、αc,i, R and T represent that anode electromigration number, cathodic electromigration number, universal gas constant and battery are actual respectively
Temperature.Wherein αa,i、αc,iValue with R is respectively 0.5,0.5 and 8.314 (Jmol-1·K-3)。
Wherein current density j:,i, accounting equation is as follows：
Subscript surf represents particle surface,Represent just
Pole solid phase particles surface lithium concentration or negative pole solid phase particles surface lithium concentration.
Overpotential ηiEqual to solid phase potentialSubtract liquid phase potentialEquilibrium potential U is subtracted againiI.e.：
Wherein equilibrium potential UiCalculating consider the influence of temperature change, i.e.,：
TrefReference temperature is represented, value is 298K, Uref,iRepresent positive pole open-circuit voltage U under reference temperatureref,pOr ginseng
Negative pole open-circuit voltage U at a temperature of examiningref,n；Represent positive pole equilibrium potential temperature coefficientOr negative pole equilibrium potential temperature
Spend coefficient
D：Electron charge conservation is calculated with Ohm's law on negative current collector, negative pole, positive pole and plus plate current-collecting body, i.e.,：
Absorbing boundary equation is：
k1,iRepresent positive pole reaction rate k1,pOr positive pole reaction rate k1,n,For solid phase potential, UappFor battery discharge
Voltage.
L represents thickness in formula, and L subscript n cc, n, sep, p represent negative current collector, negative pole, barrier film, positive pole respectively；
E：With concentrated solution theoretical description ionic charge conservation, i.e.,：
Absorbing boundary equation is：
Liquid phase potential, f electrolyte activation coefficients, value is 1, t+For lithium ion carry-over factor；
F：Battery heat in charge and discharge process is described on negative pole and positive pole, the battery heat includes electrochemical reaction
Heat, ohm heat and polarization heat；Ohm heat is described on negative current collector, the description polarization heat on barrier film, just
Ohm heat is described on the collector of pole；
Electrochemical reaction heat：
Ohm heat：
Polarize heat：
k2Represent ionic conductivity,
C, D, r, ε, t, R in A, B, C, D, E, F formulaiRepresent lithium concentration respectively, it is electrolyte diffusion coefficient, anti-
Answer interface radius, volume fraction, discharge time and positive and negative pole material grain diameter；Subscript 12 represents solid phase, liquid phase, subscript respectively
I represents equation institute loading area i.e. positive pole or negative pole, and T represents temperature；
2. Three Dimensional Thermal coupling model：Three-dimensional batteries geometrical model is set up according to the three-dimensional dimension of base batteries specification, with electricity
The average heat of electrode obtained by chemical model is calculated is thermal source Q, and the thermal field of heat convection is loaded on 3-D geometric model, is obtained
To three-dimensional thermal model, carried out by finite element theory after gridding, calculate thermal field equation, and then obtain the change of temperature, by this
Temperature change is fed back in electrochemical model, realizes the bidirectional couple of electrochemical model and thermal model；
The three-dimensional batteries geometrical model includes herein below：
Thermal source Q accounting equation is：
The equation for describing the thermal field of heat convection is as follows：
ρ, C in formulap, K represent the density, specific heat capacity, thermal conductivity factor of electrode material respectively, and the relevant parameter of electrode material
As shown in table 5.Boundary condition is represented with Newtonian Cooling formula, i.e.,：
Wherein h is natural heat-exchange coefficient, value 7.17W/ (Km2), TambFor environment temperature, value is 298K.
4) verification step 3) in electrochemical heat coupling model validity, discharged by modeling battery with different multiplying
The maximum temperature of voltage and discharge process with minimum temperature result of calculation compared with experimental measurements under the same terms, phase
To error ＜ 1%；
5) optimization aim is set according to customer demand, using thickness of electrode as optimized variable, in the reasonable design limiting bar of electrode
Under part, establish the optional scope of thickness of electrode, by BOBYQA algorithms (quasi- two-dimensional linear border optimization) to step 3) in mould
Type is calculated, the battery electrode thickness after being optimized；The BOBYQA algorithms are with the thickness of electrode of base batteries specification
For initial value, substitute into electrochemical heat coupling model, the value of calculating target function feeds back in BOBYQA algorithms, is based on
The next thickness of electrode design of BOBYQA algorithms selections substitutes into the value of calculating target function in electrochemical heat coupling, then feeds back to
In BOBYQA algorithms, next iteration is judged whether to according to whether object function maximizes, until object function is maximized,
And provide the thickness of electrode design of last time iteration, that is, the thickness of electrode design after optimizing.
The method of above-mentioned optimization battery electrode thickness, it is preferred that step 5) in the reasonable design limiting condition of electrode include
Negative pole theoretical capacity is 1.1~1.2 times of positive pole theoretical capacity, 50 DEG C of maximum temperature ＜ in discharge process；I.e. 50 DEG C of maximum temperature T ＜ in three-dimensional batteries geometrical model.
Compared with prior art, the advantage of the invention is that：The present invention is by setting up electrochemical heat coupling model, in model
Middle increase optimization module, only need to prepare less battery sample and model is verified, model can provide setting optimization aim and limit
Under fixed condition, thickness of electrode optimal design parameters, and need not be during the designing and developing of new product to all new models
Lithium battery all carries out complicated, cumbersome electrochemical property test, therefore, it is possible to effectively shorten the exploitation week of new material or new product
Phase, development cost is reduced, there is certain directive significance for the exploitation of new material or new product.
Brief description of the drawings
Fig. 1 optimizes the schematic flow sheet of the method for battery electrode thickness for the present invention.
Fig. 2 is the structural representation of accurate two-dimentional electrochemical model in the present invention.
Fig. 3 is the three-dimensional thermal model and its grid schematic diagram of embodiment 1 in the present invention.
Fig. 4 is the simulation electric discharge in the embodiment of the present invention 1 when battery different multiplying under the conditions of room temperature natural cooling is discharged
Curve and experiment discharge curve.
Fig. 5 is 1C electric discharge finish time battery infrared thermal imaging figures in the embodiment of the present invention 1.
Embodiment
For the ease of understanding the present invention, present invention work more comprehensively, is meticulously described below in conjunction with preferred embodiment,
But protection scope of the present invention is not limited to embodiment in detail below.
Unless otherwise defined, the implication that all technical terms used hereinafter are generally understood that with those skilled in the art
It is identical.Technical term used herein is intended merely to describe the purpose of specific embodiment, is not intended to the limitation present invention
Protection domain.
Embodiment
A kind of method for optimizing battery electrode thickness, comprises the following steps：1) maximized with energy density E or power is close
Spend P to maximize as the optimization aim for this method, the wherein maximized object function of energy density is
The maximized object functions of power density P are：：
Wherein, V, M represent discharge voltage and electrode quality respectively, and i is discharge current, and t is discharge time；According to client's
Demand, selects battery specifications close therewith as the base batteries specification of this method, battery electrode is carried out on this basis thick
The optimization of degree；
2) electrode structural designs parameter, the kinetic parameter of electrode material and the hot physical property ginseng of base batteries specification are obtained
Number；
3) battery electrochemical thermal coupling model is set up, the electrochemical heat coupling model includes an accurate two-dimentional electrochemistry mould
Type and a Three Dimensional Thermal coupling model；
1. accurate two-dimentional electrochemical model：As shown in Fig. 2 electrochemical model is that the different parts of electrode are represented with line segment,
One-dimensional geometrical model is set up, the differential or partial differential equation of loading description discharge process on one-dimensional geometrical model, by having
The first thought of limit is carried out after gridding to geometrical model, is calculated the differential or partial differential equation of description discharge process, is obtained at one
Cell discharge voltage curve under thickness of electrode；
The one-dimensional geometrical model includes herein below：
A, lithium ion is described in solid phase internal transmission on negative pole and positive pole with Fick's first law：
Boundary condition equation is：
B, describes liquid phase lithium ion course of dissolution on positive pole, negative pole and barrier film：
T in formula+、F、jlo4,iLithium ion carry-over factor, Faraday constant and local current densities are represented respectively；
C：Electrochemical reaction is described with Butler-Volmer equations on positive pole and negative pole
αa,i、αc,i, R and T represent that anode electromigration number, cathodic electromigration number, universal gas constant and battery are actual respectively
Temperature.Wherein αa,i、αc,iValue with R is respectively 0.5,0.5 and 8.314 (Jmol-1·K-3)。
Wherein current density j:,i, accounting equation is as follows：
Overpotential ηiEqual to solid phase potentialSubtract liquid phase potentialEquilibrium potential U is subtracted againiI.e.：
Wherein equilibrium potential UiCalculating consider the influence of temperature change, i.e.,：
TrefReference temperature is represented, value is 298K, Uref,iRepresent positive pole open-circuit voltage U under reference temperatureref,pOr ginseng
Negative pole open-circuit voltage U at a temperature of examiningref,n；Represent positive pole equilibrium potential temperature coefficientOr negative pole equilibrium potential temperature
Spend coefficient
D：Electron charge conservation is calculated with Ohm's law on negative current collector, negative pole, positive pole and plus plate current-collecting body, i.e.,：
Absorbing boundary equation is：
k1,iRepresent positive pole reaction rate k1,pOr positive pole reaction rate k1,n,Solid phase potential, UappFor battery discharge electricity
Pressure.
L represents thickness in formula, and L subscript n cc, n, sep, p represent negative current collector, negative pole, barrier film, positive pole respectively；
E：With concentrated solution theoretical description ionic charge conservation, i.e.,：
Absorbing boundary equation is：
For liquid phase potential, f electrolyte activation coefficients, value is 1, t+For lithium ion carry-over factor；
F：Battery heat in charge and discharge process is described on negative pole and positive pole, the battery heat includes electrochemical reaction
Heat, ohm heat and polarization heat；Ohm heat is described on negative current collector, the description polarization heat on barrier film, just
Ohm heat is described on the collector of pole；
Electrochemical reaction heat：
Ohm heat：
Polarize heat：
k2Represent ionic conductivity,
C, D, r, ε, t, R in A, B, C, D, E, F formulaiRepresent lithium concentration respectively, it is electrolyte diffusion coefficient, anti-
Answer interface radius, volume fraction, discharge time and positive and negative pole material grain diameter；Subscript 12 represents solid phase, liquid phase, subscript respectively
I represents equation institute loading area i.e. positive pole or negative pole, and T represents temperature；
2. Three Dimensional Thermal coupling model：Three-dimensional batteries geometrical model is set up according to the three-dimensional dimension of base batteries specification, with electricity
The average heat of electrode obtained by chemical model is calculated is thermal source Q, and the thermal field of heat convection is loaded on 3-D geometric model, is obtained
To three-dimensional thermal model, carried out by finite element theory after gridding, three-dimensional thermal model and its grid schematic diagram as shown in Figure 3；
Thermal field equation is calculated, and then obtains the change of temperature, this temperature change is fed back in electrochemical model, electrochemical model is realized
With the bidirectional couple of thermal model；
The three-dimensional batteries geometrical model includes herein below：
Thermal source Q accounting equation is：
The equation for describing the thermal field of heat convection is as follows：
ρ, C in formulap, K represent the density, specific heat capacity, thermal conductivity factor of electrode material respectively, and the relevant parameter of electrode material
As shown in table 5.Boundary condition is represented with Newtonian Cooling formula, i.e.,：
Wherein h is natural heat-exchange coefficient, value 7.17W/ (Km2), TambEnvironment temperature, value is 298K.
4) verification step 3) in electrochemical heat coupling model validity, discharged by modeling battery with different multiplying
The maximum temperature of voltage and discharge process with minimum temperature result of calculation compared with experimental measurements under the same terms, phase
To error ＜ 1%；
5) optimization aim is set according to customer demand, using thickness of electrode as optimized variable, in the reasonable design limiting bar of electrode
Under part, establish the optional scope of thickness of electrode, by BOBYQA algorithms (quasi- two-dimensional linear border optimization) to step 3) in mould
Type is calculated, the battery electrode thickness after being optimized；The BOBYQA algorithms are with the thickness of electrode of base batteries specification
For initial value, substitute into electrochemical heat coupling model, the value of calculating target function feeds back in BOBYQA algorithms, is based on
The next thickness of electrode design of BOBYQA algorithms selections substitutes into the value of calculating target function in electrochemical heat coupling, then feeds back to
In BOBYQA algorithms, next iteration is judged whether to according to whether object function maximizes, until object function is maximized,
And provide the thickness of electrode design of last time iteration, that is, the thickness of electrode design after optimizing.
The method of above-mentioned optimization battery electrode thickness, it is preferred that step 5) in the reasonable design limiting condition of electrode include
Negative pole theoretical capacity is 1.1~1.2 times of positive pole theoretical capacity, 50 DEG C of maximum temperature ＜ in discharge process；I.e. 50 DEG C of maximum temperature T ＜ in three-dimensional batteries geometrical model.
Embodiment 1
By taking ferric phosphate lithium cell system as an example, the battery size requirement of client is as shown in table 1, and battery applications are usually that 1C is put
Electricity is, it is necessary to which the battery energy density is maximum, i.e., under 1C discharging conditions, battery core energy density is maximized.
The battery three-dimensional dimension parameter of table 1
| Width | Thickness | Highly |
| 100mm | 12mm | 115mm |
According to customer demand, selection requires the battery electrode design ginseng selected in close battery specifications, this example therewith
Number, i.e. battery structure parameter, as shown in table 2,1C electric discharges are that battery core energy density is 159.40Wh/kg under the design.
The cell electrode structure design parameter of table 2
| Parameter name | Negative pole | Positive pole | Barrier film | Negative current collector | Plus plate current-collecting body |
| Solid volume fraction | 0.55 | 0.43 | \ | \ | \ |
| Liquid phase volume fraction | 0.33 | 0.332 | 0.54 | \ | \ |
| Thickness of electrode (μm) | 40 | 70 | 25 | 12 | 20 |
| Mean particle radius (μm) | 6.0 | 0.08 | \ | \ | \ |
| Maximum lithium concentration (mol m-3) | 31370 | 22806 | \ | \ | \ |
Remarks：Slash represents to be not present or does not consider the parameter
Battery electrochemical thermal coupling simulation model is set up based on this design, and verifies model validation.It is electric in this example
The foundation of pond electrochemical heat Coupling Simulation Model is to be based on COMSOL platforms, and model manipulation step is as follows：
1) by test measurement or literature survey method obtain electrode structural designs parameter needed for model, electrode material it is dynamic
Mechanics parameter and thermal physical property parameter.
As shown in table 2, the kinetic parameter of battery electrode material is as shown in table 3, battery for cell electrode structure design parameter
The thermal physical property parameter of electrode material is as shown in table 5.
The kinetic parameter of the battery electrode material of table 3
2) as shown in Fig. 2 electrochemical model is that the different parts of electrode are represented with line segment, one-dimensional geometrical model is set up,
The differential or partial differential equation of loading description discharge process on geometrical model, and geometrical model is entered by finite element theory
After row grid is cutd open, then the equation of description discharge process is solved, when can obtain different electric discharges in the discharge process under thickness of electrode design
Between corresponding discharge voltage.Shown in dependent equation following A, B, C, D, E, F, area is calculated needed for physical field loading area, i.e. equation
Domain distribution is as shown in table 4.
A：Lithium ion is described with Fick's first law in solid phase internal transmission, equation is as follows：
Boundary condition equation is：
C, D, r, ε, t, R in formulaiRespectively represent lithium concentration, diffusion coefficient, reaction interface radius, volume fraction, when
Between and positive and negative pole material grain diameter；Subscript 1,2, i represent solid phase, liquid phase and different battery compositions respectively.I is added by equation
The region of load, it is specific as shown in table 4.
B：Liquid phase lithium ion course of dissolution：
T in formula+、F、jLoc, iTransference number of ions, Faraday constant and local current densities are represented respectively
C：Butler-Volmer equations describe electrochemical reaction
Wherein current density j:,i, accounting equation is as follows：
Overpotential ηiLiquid phase potential, which is subtracted, equal to solid phase potential subtracts equilibrium potential again i.e.：
Wherein equilibrium potential UiCalculating consider the influence of temperature change, i.e.,：
D：Ohm's law calculates electron charge conservation, i.e.,：
Absorbing boundary equation is：
L represents thickness in formula, and subscript n cc, n, sep, p represent negative current collector, negative pole, barrier film, positive pole respectively.
E：Concentrated solution theoretical description ionic charge conservation, i.e.,：
Absorbing boundary equation is：
F：Battery heat in charge and discharge process
Electrochemical reaction heat：
Ohm heat：
Polarize heat：
Equation zoning is distributed on the electrode of table 4
| A | B | C | D | F | |
| Negative current collector | Nothing | Nothing | Nothing | Have | Only Qohmic,i |
| Negative pole | Have | Have | Have | Have | Have |
| Barrier film | Nothing | Have | Nothing | Nothing | Only Qpolar,i |
| Positive pole | Have | Have | Have | Have | Have |
| Plus plate current-collecting body | Nothing | Nothing | Nothing | Have | Only Qohmic,i |
3) according to the parameter of table 1, three-dimensional batteries geometrical model is set up, the average heat of the electrode obtained is calculated with electrochemical model
Q is thermal source, and the thermal field that heat convection is loaded on 3-D geometric model obtains three-dimensional thermal model, carries out after mesh generation, calculates
Thermal field equation, and then obtain the change of temperature, and by the temperature Real-time Feedback into electrochemical model, realize electrochemical model and
The bidirectional couple of thermal model.
Thermal source Q accounting equation is：
The equation for describing heat convection is as follows：
ρ, C in formulap, K represent the density, specific heat capacity, thermal conductivity factor of electrode material respectively, and the relevant parameter of electrode material
As shown in table 5.Boundary condition is represented with Newtonian Cooling formula, i.e.,：
Wherein h is natural heat-exchange coefficient, value 7.17W/ (Km2)。
The thermal physical property parameter of the lithium ion battery electrode material of table 5
Verify the validity of electrochemical heat coupling model.In this example, using 2.0V as discharge cut-off voltage, in environment temperature
Under the conditions of 25 DEG C, natural cooling, battery is tested under 0.5C, 1C, 3C, 5C modeling discharge voltage profile and the same terms
Discharge curve is as shown in Figure 4.Release capacity calculation expression beI, t represent discharge current and discharge time respectively.
Compare phase in the same time, simulation calculates discharge voltage and the relative error of measurement result is 0.93% to the maximum.Battery 1C electric discharges terminate
The infrared thermal imaging figure at moment, i.e. temperature highest moment is as shown in Figure 5.Maximum temperature mould in battery discharge procedure under different multiplying
Intend result and experimental measurements contrast is as shown in table 6.Relative error is both less than 1% and shows that model is effective.
Maximum temperature analog result and experimental measurements in the discharge process of table 6
Increase optimization module in a model, optimization aim is turned to energy density (E) maximum.Object function in this example
For:
Wherein, V, M represent discharge voltage and electrode quality respectively.
Optimized algorithm selects BOBYQA, and Optimal Parameters are positive pole thickness (Lp) and negative pole thickness (Ln), according to client in capacity
With estimated under conditions of volume defining span be [30,110] and [20,60], unit be μm.
Qualifications are：Maximum temperature ＜ in three-dimensional thermal model
50℃。
Using thickness of electrode in table 2 as initial value, substitute into electrochemical heat coupling model, the value of calculating target function, feedback
Into BOBYQA algorithms, substituted into based on the next thickness of electrode design of BOBYQA algorithms selections in electrochemical heat coupling and calculate target
The value of function, then feed back in BOBYQA algorithms, next iteration is judged whether to according to whether object function maximizes, directly
Maximized to object function, and provide the thickness of electrode design of last time iteration, that is, the thickness of electrode design after optimizing.
BOBYQA algorithms provide optimal design：Positive pole thickness is 67.1 μm, and negative pole thickness is 38.4 μm, electrochemical heat coupling
The battery core energy density that matched moulds type calculates under thickness of electrode design is 164.80Wh/kg.
Using positive pole thickness as 67.1 μm, negative pole thickness is 38.4 μm, other design parameters and preparation technology with it is originally identical,
It is made under battery sample, 1C discharging conditions, the battery core mean energy density of battery sample is 164.34Wh/kg.Higher than most primary election
Select the battery core energy density (159.40Wh/kg) of battery specifications.
Claims (2)
1. a kind of method for optimizing battery electrode thickness, it is characterised in that：Comprise the following steps：1) maximized with energy density E
Or power density P is maximized as the optimization aim for this method, the wherein maximized object function of energy density isThe maximized object functions of power density P are：
Wherein, V, M represent discharge voltage and electrode quality respectively, and i is discharge current, and t is discharge time；According to the need of client
Ask, select battery specifications close therewith as the base batteries specification of this method, battery electrode thickness is carried out on this basis
Optimization；
2) electrode structural designs parameter, the kinetic parameter and thermal physical property parameter of electrode material of base batteries specification are obtained；
3) battery electrochemical thermal coupling model is set up, the electrochemical heat coupling model is an accurate two-dimentional electrochemical model and one
The coupling of individual three-dimensional thermal model；
1. accurate two-dimentional electrochemical model：The different parts of electrode are represented with line segment, one-dimensional geometrical model are set up, in one-dimensional geometry
The differential or partial differential equation of loading description discharge process, gridding is carried out by finite element theory to geometrical model on model
Afterwards, the differential or partial differential equation of description discharge process are calculated, the discharge voltage profile of the battery under a thickness of electrode is obtained；
2. three-dimensional thermal model：Three-dimensional batteries geometrical model is set up according to the three-dimensional dimension of base batteries specification, with electrochemical model
The average heat of electrode obtained by calculating is thermal source Q, and the thermal field of heat convection is loaded on 3-D geometric model, Three Dimensional Thermal is obtained
Model, is carried out after gridding by finite element theory, calculates thermal field equation, and then obtains the change of temperature, by this temperature change
Feed back in electrochemical model, realize the bidirectional couple of electrochemical model and thermal model；
4) verification step 3) in electrochemical heat coupling model validity, by modeling battery with different multiplying discharge voltage
And the maximum temperature of discharge process with minimum temperature result of calculation compared with experimental measurements under the same terms, it is relative by mistake
Poor ＜ 1%；
5) optimization aim is set according to customer demand, using thickness of electrode as optimized variable, under the conditions of the reasonable design limiting of electrode,
Establish the optional scope of thickness of electrode, by BOBYQA algorithms (quasi- two-dimensional linear border optimization) to step 3) in model carry out
Calculate, the battery electrode thickness after being optimized；It using the thickness of electrode of base batteries specification is initial that the BOBYQA algorithms, which are,
Value, is substituted into electrochemical heat coupling model, the value of calculating target function is fed back in BOBYQA algorithms, based on BOBYQA algorithms
Select next thickness of electrode design to substitute into the value of calculating target function in electrochemical heat coupling, then feed back to BOBYQA algorithms
In, next iteration is judged whether to according to whether object function maximizes, until object function is maximized, and provides last
The thickness of electrode design of an iteration, that is, the thickness of electrode design after optimizing.
2. the method for optimization battery electrode thickness according to claim 1, it is characterised in that：Step 5) in electrode rationally set
Counting qualifications includes negative pole theoretical capacity for 1.1~1.2 times of positive pole theoretical capacity, maximum temperature ＜ 50 in discharge process
℃；I.e.50 DEG C of maximum temperature T ＜ in three-dimensional batteries geometrical model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710209193.5A CN107145629A (en) | 2017-03-31 | 2017-03-31 | A kind of method for optimizing battery electrode thickness |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710209193.5A CN107145629A (en) | 2017-03-31 | 2017-03-31 | A kind of method for optimizing battery electrode thickness |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107145629A true CN107145629A (en) | 2017-09-08 |
Family
ID=59784187
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710209193.5A Pending CN107145629A (en) | 2017-03-31 | 2017-03-31 | A kind of method for optimizing battery electrode thickness |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107145629A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108595840A (en) * | 2018-04-25 | 2018-09-28 | 中国科学技术大学 | The method for establishing model of lithium ion battery, the optimization method of system and polar ear size |
| CN109446619A (en) * | 2018-10-18 | 2019-03-08 | 中国科学技术大学 | A kind of optimization method of lithium ion cell electrode design parameter |
| CN110232201A (en) * | 2019-04-02 | 2019-09-13 | 中南大学 | A kind of battery design method of multi-parameter synergistic effect |
| CN110580389A (en) * | 2019-08-29 | 2019-12-17 | 中南大学 | Method for optimizing battery tab |
| CN111062137A (en) * | 2019-12-23 | 2020-04-24 | 荣盛盟固利新能源科技有限公司 | Lithium ion battery performance prediction model and construction method and application thereof |
| CN111682255A (en) * | 2020-05-08 | 2020-09-18 | 深圳市鹏诚新能源科技有限公司 | Design method of battery |
| CN112713333A (en) * | 2020-12-31 | 2021-04-27 | 武汉理工大学 | Cooling plate structure of lithium ion square battery and cooling method thereof |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6016047A (en) * | 1996-11-21 | 2000-01-18 | U.S. Philips Corporation | Battery management system and battery simulator |
| CN104849675A (en) * | 2015-06-17 | 2015-08-19 | 哈尔滨工业大学 | Method for obtaining electrochemical and thermal coupling models of lithium ion 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 |
| CN105550452A (en) * | 2015-12-18 | 2016-05-04 | 中国科学院上海高等研究院 | Lithium ion battery P2D model parameter recognizing method based on heuristic algorithm |
-
2017
- 2017-03-31 CN CN201710209193.5A patent/CN107145629A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6016047A (en) * | 1996-11-21 | 2000-01-18 | U.S. Philips Corporation | Battery management system and battery simulator |
| 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 |
| CN104849675A (en) * | 2015-06-17 | 2015-08-19 | 哈尔滨工业大学 | Method for obtaining electrochemical and thermal coupling models of lithium ion battery |
| CN105550452A (en) * | 2015-12-18 | 2016-05-04 | 中国科学院上海高等研究院 | Lithium ion battery P2D model parameter recognizing method based on heuristic algorithm |
Non-Patent Citations (1)
| Title |
|---|
| 汤依伟等: "锂离子动力电池高倍率充放电过程中弛豫行为的仿真", 《物理学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108595840A (en) * | 2018-04-25 | 2018-09-28 | 中国科学技术大学 | The method for establishing model of lithium ion battery, the optimization method of system and polar ear size |
| CN109446619A (en) * | 2018-10-18 | 2019-03-08 | 中国科学技术大学 | A kind of optimization method of lithium ion cell electrode design parameter |
| CN109446619B (en) * | 2018-10-18 | 2022-12-23 | 中国科学技术大学 | Optimization method of design parameters of lithium ion battery electrode |
| CN110232201A (en) * | 2019-04-02 | 2019-09-13 | 中南大学 | A kind of battery design method of multi-parameter synergistic effect |
| CN110580389A (en) * | 2019-08-29 | 2019-12-17 | 中南大学 | Method for optimizing battery tab |
| CN111062137A (en) * | 2019-12-23 | 2020-04-24 | 荣盛盟固利新能源科技有限公司 | Lithium ion battery performance prediction model and construction method and application thereof |
| CN111062137B (en) * | 2019-12-23 | 2024-03-12 | 荣盛盟固利新能源科技有限公司 | Lithium ion battery performance prediction model, construction method and application thereof |
| CN111682255A (en) * | 2020-05-08 | 2020-09-18 | 深圳市鹏诚新能源科技有限公司 | Design method of battery |
| CN112713333A (en) * | 2020-12-31 | 2021-04-27 | 武汉理工大学 | Cooling plate structure of lithium ion square battery and cooling method thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107145629A (en) | A kind of method for optimizing battery electrode thickness | |
| CN110232201A (en) | A kind of battery design method of multi-parameter synergistic effect | |
| CN107145628B (en) | Method for predicting cycle life of lithium battery based on electrochemical-thermal coupling model | |
| He et al. | A 3D electrochemical-thermal coupled model for electrochemical and thermal analysis of pouch-type lithium-ion batteries | |
| CN109446619B (en) | Optimization method of design parameters of lithium ion battery electrode | |
| US10539621B2 (en) | Method and apparatus for identifying a battery model | |
| Golmon et al. | A design optimization methodology for Li+ batteries | |
| Li et al. | Physics-based CFD simulation of lithium-ion battery under the FUDS driving cycle | |
| CN111062137B (en) | Lithium ion battery performance prediction model, construction method and application thereof | |
| Yang et al. | An improved semi-empirical model for thermal analysis of lithium-ion batteries | |
| Romero-Becerril et al. | Comparison of discretization methods applied to the single-particle model of lithium-ion batteries | |
| Jiang et al. | Lattice Boltzmann simulation of ion and electron transport during the discharge process in a randomly reconstructed porous electrode of a lithium-ion battery | |
| CN109614754B (en) | Simulation method of three-dimensional simplified model of lithium ion battery | |
| Hellwig et al. | A multi-scale electrochemical and thermal model of a LiFePO4 battery | |
| CN109738806B (en) | Method, device and medium for simulating heat generation rate of battery | |
| CN109581064A (en) | Measure method, apparatus, the medium of battery low temperature direct resistance | |
| CN109033619A (en) | A kind of transient temperature model modelling approach of 18650 type lithium battery discharge cycles | |
| Wang et al. | An experimental and numerical examination on the thermal inertia of a cylindrical lithium-ion power battery | |
| Huang et al. | A co-simulation method based on coupled thermoelectric model for electrical and thermal behavior of the lithium-ion battery | |
| CN114547903A (en) | Method for predicting cycle life of lithium battery based on electrochemical-thermal coupling model | |
| Fang et al. | Performance simulation method and state of health estimation for lithium-ion batteries based on aging-effect coupling model | |
| An et al. | A comparison of model prediction from P2D and particle packing with experiment | |
| Tian et al. | Parallel-connected battery module modeling based on physical characteristics in multiple domains and heterogeneous characteristic analysis | |
| Qiu et al. | Understanding thermal and mechanical effects on lithium plating in lithium-ion batteries | |
| Du et al. | Optimization of LiMn 2 O 4 electrode properties in a gradient-and surrogate-based framework |
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 | ||
| AD01 | Patent right deemed abandoned | ||
| AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20210330 |