METHOD AND SYSTEM FOR BUILDING THERMAL MODEL OF POWER LITHIUM-ION BATTERY BASED ON ELECTROCHEMICAL MECHANISM

Abstract

A method and system for building a thermal model of a power lithium-ion battery based on an electrochemical mechanism. The method includes: discretizing a second order partial differential heat conduction equation of a power lithium-ion battery according to a finite differential method, thereby building a thermal model of the power lithium-ion battery; carrying out a dynamic working condition test by using a cylindrical power lithium-ion battery selected as an object, thereby acquiring experimental data such as a temperature, a current, a voltage, and a temperature of a surface of the battery; identifying an electrochemical parameter of the power lithium-ion battery according to an optimal parameter algorithm by using test data acquired in a dynamic working condition, thereby building a thermal model of the power lithium-ion battery; and verifying accuracy of the thermal model of the power lithium-ion battery by using test data acquired in another dynamic working condition.

Claims

1. A method for building a thermal model of a power lithium-ion battery based on an electrochemical mechanism, comprising: step S1: discretizing a second order partial differential heat conduction equation of a power lithium-ion battery according to a finite differential method, thereby building a thermal model of the power lithium-ion battery; step S2: carrying out a dynamic working condition test by using a cylindrical power lithium-ion battery selected as an object, thereby acquiring test data such as a temperature, a current, a voltage, and a temperature of a surface of the battery; step S3: identifying an electrochemical parameter of the power lithium-ion battery according to an optimal parameter algorithm by using test data acquired in a dynamic working condition, thereby building a thermal model of the power lithium-ion battery; and step S4: verifying accuracy of the thermal model of the power lithium-ion battery by using test data acquired in another dynamic working condition; step S5: acquiring a temperature of a specific power lithium-ion battery by using a verified thermal model of the power lithium-ion battery, comparing the acquired temperature with a predetermined threshold, and controlling a switch of the specific power lithium-ion battery to turn off according to a comparison result.

2. The method according to claim 1, wherein step S1 comprises: step S1.1: discretizing a second order partial differential import equation of a cylindrical power lithium-ion battery according to the finite differential method, thereby building a thermal model of a one-dimensional state space of the cylindrical power lithium-ion battery; and step S1.2: determining an impact of a temperature on an electrochemical parameter of the battery according to an Arrhenius equation, thereby establishing a coupling relationship of the temperature to the electrochemical parameter of the battery.

3. The method according to claim 2, wherein step S1.1 specifically comprises: step S1.1.1: assuming that a temperature distribution of the cylindrical power lithium-ion battery complies with the following one-dimensional unsteady-state heat conduction equation of cylindrical coordinates: $ρ c \frac{\partial T}{\partial t} = \frac{λ}{r} \frac{\partial T}{\partial r} + λ \frac{\partial^{2} T}{\partial r^{2}} + q;$ building a boundary condition: ${\begin{matrix} {- λ \frac{\partial T}{\partial r} .Math.}_{r_{0}} = h_{0} (T_{0} - T_{1}) = 0 \\ {- λ \frac{\partial T}{\partial r} .Math.}_{R} = h (T_{R} - T_{amb}) \end{matrix};$ building an initial condition: T(t)=T.sub.amb; and building a supplementary condition: T.sub.1(t)≈T.sub.0(t), wherein T.sub.0 denotes a temperature of a thin air layer that is of a hollow portion of the battery and that is closest to an inner-most layer of the battery; T.sub.1 denotes a temperature of the inner-most layer of the battery; T.sub.R denotes a temperature of a surface of the battery; h.sub.0 denotes a convection-diffusion coefficient of the thin air layer; h denotes a convection-diffusion coefficient of the surface of the battery; $a = \frac{λ}{ρ c}$ denotes a thermal diffusivity; ρ denotes a density of the power lithium-ion battery; c denotes a specific heat capacity J/(kg.Math.° C.) of the power lithium-ion battery; and λ denotes a radial thermal conductivity; step S1.1.2: respectively approximating a first-order partial differential equation and a second-order partial differential equation by using a backward difference method and a central difference method, thereby discretizing the second order partial differential heat conduction equation: $\frac{\partial T}{\partial r} \approx \frac{T_{k + 1} - T_{k - 1}}{2 Δ r}, \frac{\partial^{2} T}{\partial r^{2}} \approx \frac{T_{k + 1} - 2 T_{k} + T_{k - 1}}{{(Δ r)}^{2}};$ step S1.1.3: representing the one-dimensional unsteady-state heat conduction equation of the cylindrical power lithium-ion battery as follows: ${\dot{T}}_{k} = \frac{a_{k}}{r_{k}} \frac{T_{k + 1} - T_{k - 1}}{2 Δ r} + a_{k} \frac{T_{k + 1} - 2 T_{k} + T_{k - 1}}{{(Δ r)}^{2}} + \frac{a_{k}}{λ} q (t), k = 1, 2, 3 .Math., M - 1,$ wherein $Δ r = \frac{R - r_{0}}{M}; r_{k} = r_{0} + k Δ r;$ M denotes a quantity of walls of the cylindrical power lithium-ion battery; and r.sub.M=R; and when $b_{k} = \frac{a_{k}}{2 r_{k} Δ r} and d_{k} = \frac{2 a_{k}}{{(Δ r)}^{2}},$ ${\dot{T}}_{k} = (b_{k} + d_{k}) T_{k + 1} - 2 d_{k} T_{k} + (- b_{k} + d_{k}) T_{k - 1} + \frac{a_{k}}{λ} q (t);$ and representing the temperature T.sub.M of the surface of the battery as follows: $T_{k = M} = \frac{λ_{M}}{h Δ r_{M} + λ_{M}} T_{k = M - 1} + \frac{h Δ r_{M}}{h Δ r_{M} + λ_{M}} T_{amb},$ wherein q(t)=q.sub.p+q.sub.r; the equation is solved in step S1.1.4; T.sub.M=T.sub.R; T.sub.M denotes a temperature of an outer layer of the battery; λ.sub.M denotes a thermal conductivity of an outer-layer material; λr.sub.M denotes a thickness of the outer-layer material; and q(t) denotes a heating power of the battery at moment t; step S1.1.4: according to a polarization phenomenon and a current heating effect in charging and discharging processes of the battery, building the following calculation formula for a heat release rate of polarization heat and ohmic heat in a use process of the battery:
q.sub.p=I.sup.2R.sub.act+I.sup.2R.sub.ohm=I.sup.2R.sub.t, wherein I denotes a current of the battery; R.sub.act denotes a polarization internal resistance of the battery; R.sub.ohm denotes an ohmic internal resistance of the battery; R.sub.t denotes a total internal resistance of the battery; and q.sub.p denotes a heat release rate of polarization heat and ohmic heat in a discharging process of the battery; and combining calculation steps to acquire the following total heating power:
q=q.sub.r+q.sub.p, wherein q denotes a total heating power; q.sub.r denotes a heat release rate of reaction heat of the battery; and q.sub.p denotes the heat release rate of the polarization heat and ohmic heat in the discharging process of the battery; and step S1.1.5: determining a final model output according to an actual design demand for a system; carrying out approximation according to the finite differential method; iteratively updating an electrochemical parameter of the battery according to an Arrhenius equation; and calculating a system output temperature at a current moment; defining a system state equation: {dot over (x)}=Ax+Bu; and defining a system output equation: y=Cx+Du, wherein A and B denote system matrices: a system state is denoted by x=(T.sub.1, T.sub.2, . . . , T.sub.i, . . . , T.sub.M-1).sup.T; a system input is denoted by q; and system output y denotes a temperature of an (M−1).sup.th layer of the battery.

4. The method according to claim 2, wherein step S1.2 specifically comprises: based on the relationship between a temperature and an electrochemical parameter of the battery, determining an impact of a temperature on a parameter of the battery according to the following Arrhenius equation: $ψ = ψ_{ref} \exp [\frac{E_{a}^{ψ}}{R} (\frac{1}{T_{ref}} - \frac{1}{T})],$ wherein ψ.sub.ref is a general variable denoting a diffusion coefficient of a substance, an electrical conductivity of an electrolyte, an exchange current density of an electrode reaction electrode reaction, or the like; subscript ref denotes a value at a reference temperature; and E.sub.a.sup.ψ denotes an activation energy.

5. The method according to claim 3, wherein in the step S1.1.5, each system matrix is denoted as follows: $A = (- (b_{1} + d_{1}), (b_{1} + d_{1}), 0, .Math., 0; - (b_{2} + d_{2}), - 2 d_{2}, (b_{2} + d_{2}), .Math., 0 : .Math.; 0, 0, 0, .Math., - (b_{M - 1} + d_{M - 1}), (\frac{λ_{M - 1} (b_{M - 1} + d_{M - 1})}{h Δ r + λ_{M - 1}} - 2 d_{M - 1}))$ $B (i, 1) = ρ_{i} c_{i};$ $C = {(0, 0, .Math., 0, 1)}^{1}; and$ $D = 0,$ wherein a temperature denoted by system output y may be determined by adjusting a location of 1 in matrix C.

6. The method according to claim 1, wherein step S2 comprises: step S2.1: selecting a power lithium-ion battery to be tested; and adhering thermocouples to the power lithium-ion battery according to a layout solution; step S2.2: allowing the battery to stand still in a 25° C. incubator for 2 h; step S2.3: charging the battery in a constant-current and constant-voltage manner to a fully charged state, namely, SOC=100%; discharging the battery at a speed of C/3 to SOC=95%; and allowing the battery to stand still for 2 h; step S2.4: loading a dynamic working condition (UDDS) to the battery by an appropriate proportion till SOC of the battery decreases to about 5%; step S2.5: recording data such as a current, a voltage, an ambient temperature, and a surface temperature in the working condition; step S2.6: repeating steps S2.2 to S2.5 at the same ambient temperature; and acquiring test data at the temperature in dynamic working conditions such as FUDS and UDDS; and step S2.7: changing the temperature of the incubator to 5° C., 10° C., and 35° C.; repeating steps S2.2 to 52.6; and acquiring test data at each of the temperatures in the dynamic working conditions.

7. The method according to claim 1, wherein the optimal parameter algorithm in step S3 is a least square method.

8. A system for building a thermal model of a power lithium-ion battery based on an electrochemical mechanism, comprising: module M1 configured to: discretize a second order partial differential heat conduction equation of a power lithium-ion battery according to a finite differential method, thereby building a thermal model of the power lithium-ion battery; a module M2 configured to: carry out a dynamic working condition test by using a cylindrical power lithium-ion battery selected as an object, thereby acquiring test data such as a temperature, a current, a voltage, and a temperature of a surface of the battery; a module M3 configured to: identify an electrochemical parameter of the power lithium-ion battery according to an optimal parameter algorithm by using test data acquired in a dynamic working condition, thereby building a thermal model of the power lithium-ion battery; and a module M4 configured to: verify accuracy of the thermal model of the power lithium-ion battery by using test data acquired in another dynamic working condition, acquire a temperature of a specific power lithium-ion battery by using a verified thermal model of the power lithium-ion battery, compare the acquired temperature with a predetermined threshold, and control a switch of the specific power lithium-ion battery to turn off, according to a comparison result.

9. The system according to claim 8, wherein the module M1 comprises: a module M1.1 configured to: discretize a second order partial differential import equation of a cylindrical power lithium-ion battery according to the finite differential method, thereby building a thermal model of a one-dimensional state space of the cylindrical power lithium-ion battery, and a module M1.2 configured to: determine an impact of a temperature on an electrochemical parameter of the battery according to an Arrhenius equation, thereby establishing a coupling relationship of the temperature to the electrochemical parameter of the battery.

10. The system according to claim 9, wherein the module M1.1 specifically comprises: a module M1.1.1 configured to: assume that a temperature distribution of the cylindrical power lithium-ion battery complies with the following one-dimensional unsteady-state heat conduction equation of cylindrical coordinates: $ρ c \frac{\partial T}{\partial t} = \frac{λ}{r} \frac{\partial T}{\partial r} + λ \frac{\partial^{2} T}{\partial r^{2}} + q;$ building a boundary condition: ${\begin{matrix} {- λ \frac{\partial T}{\partial r} .Math.}_{r_{0}} = h_{0} (T_{0} - T_{1}) = 0 \\ {- λ \frac{\partial T}{\partial r} .Math.}_{R} = h (T_{R} - T_{amb}) \end{matrix};$ building an initial condition: T(t)=T.sub.amb; and building a supplementary condition: T.sub.1(t)≈T.sub.0(t), wherein T.sub.0 denotes a temperature of a thin air layer that is of a hollow portion of the battery and that is closest to an inner-most layer of the battery; T.sub.1 and T.sub.R respectively denote a temperature of the inner-most layer of the battery and a temperature of a surface of the battery; h.sub.0 and h respectively denote a convection-diffusion coefficient of the thin air layer and that of the surface of the battery; $a = \frac{λ}{ρ c}$ denotes a thermal diffusivity; ρ denotes a density of the power lithium-ion battery; c denotes a specific heat capacity J/(kg.Math.° C.) of the power lithium-ion battery; and λ denotes a radial thermal conductivity; a module M1.1.2 configured to: respectively approximate a first-order partial differential equation and a second-order partial differential equation by using a backward difference method and a central difference method, thereby discretizing the second order partial differential heat conduction equation: $\frac{\partial T}{\partial r} \approx \frac{T_{k + 1} - T_{k - 1}}{2 Δ r}, \frac{\partial^{2} T}{\partial r^{2}} \approx \frac{T_{k + 1} - 2 T_{k} + T_{k - 1}}{{(Δ r)}^{2}};$ module M1.1.3 configured to: represent the one-dimensional unsteady-state heat conduction equation of the cylindrical power lithium-ion battery as follows: ${\dot{T}}_{k} = \frac{a_{k}}{r_{k}} \frac{T_{k + 1} - T_{k - 1}}{2 Δ r} + a_{k} \frac{T_{k + 1} - 2 T_{k} + T_{k - 1}}{{(Δ r)}^{2}} + \frac{a_{k}}{λ} q (t), k = 1, 2, 3 .Math., M - 1,$ wherein $Δ r = \frac{R - r_{0}}{M}; r_{k} = r_{0} + k Δ r;$ M denotes a quantity of walls of the cylindrical power lithium-ion battery; and r.sub.M=R; and when $b_{k} = \frac{a_{k}}{2 r_{k} Δ r} and d_{k} = \frac{2 a_{k}}{{(Δ r)}^{2}},$ ${\dot{T}}_{k} = (b_{k} + d_{k}) T_{k + 1} - 2 d_{k} T_{k} + (- b_{k} + d_{k}) T_{k - 1} + \frac{a_{k}}{λ} q (t);$ and represent the temperature T.sub.M of the surface of the battery as follows: $T_{k = M} = \frac{λ_{M}}{h Δ r_{M} + λ_{M}} T_{k = M - 1} + \frac{h Δ r_{M}}{h Δ r_{M} + λ_{M}} T_{amb},$ wherein q(t)=q.sub.p+q.sub.r; the equation is solved by using module M1.1.4; T.sub.M=T.sub.R; T.sub.M denotes a temperature of an outer layer of the battery; λ.sub.M denotes a thermal conductivity of an outer-layer material; and Arm denotes a thickness of the outer-layer material; a module M1.1.4 configured to: according to a polarization phenomenon and a current heating effect in charging and discharging processes of the battery, build the following calculation formula for a heat release rate of polarization heat and ohmic heat in a use process of the battery:
q.sub.p=I.sup.2R.sub.act+I.sup.2R.sub.ohm=I.sup.2R.sub.t, wherein I denotes a current of the battery; R.sub.act denotes a polarization internal resistance of the battery; R.sub.ohm denotes an ohmic internal resistance of the battery; R.sub.t denotes a total internal resistance of the battery; and q.sub.p denotes a heat release rate of polarization heat and ohmic heat in a discharging process of the battery; and combine calculation steps to acquire the following total heating power:
q=q.sub.r+q.sub.p, wherein q denotes a total heating power; q.sub.r denotes a heat release rate of reaction heat of the battery; and q.sub.p denotes the heat release rate of the polarization heat and ohmic heat in the discharging process of the battery; and a module M1.1.5 configured to: determine a final model output according to an actual design demand for a system; carry out approximation according to the finite differential method; iteratively update an electrochemical parameter of the battery according to an Arrhenius equation; and calculate a system output temperature at a current moment; define a system state equation: {dot over (x)}=AX+Bu; and define a system output equation: y=Cx+Du, wherein A and B denote system matrices; a system state is denoted by x=(T.sub.1, T.sub.2, . . . , T.sub.i, . . . , T.sub.M-1).sup.T; a system input is denoted by q and is a heat generation rate per unit volume; and system output y denotes a temperature of an (M−1).sup.th layer of the battery.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0064] To describe the technical solutions in the embodiments of the present disclosure or in the prior art more clearly, the accompanying drawings required for the embodiments are briefly described below. Apparently, the accompanying drawings in the following descriptions show merely some embodiments of the present disclosure, and those of ordinary skill in the art may still derive other accompanying drawings from these accompanying drawings without creative efforts.

[0065] FIG. 1 is a flowchart of building and verifying a thermal model according to an embodiment of the present disclosure;

[0066] FIG. 2 is a flowchart of details according to an embodiment of the present disclosure;

[0067] FIG. 3 is a schematic diagram of a building process of a thermal model an embodiment of the present disclosure; and

[0068] FIG. 4 is a layered diagram of an internal structural of a cylindrical power lithium-ion battery according to an embodiment of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0069] The technical solutions of the embodiments of the present disclosure are clearly and completely described below with reference to the accompanying drawings in the embodiments of the present disclosure. Apparently, the described embodiments are merely some rather than all of the embodiments of the present disclosure. All other embodiments acquired by those of ordinary skill in the art based on the embodiments of the present disclosure without creative efforts shall fall within the protection scope of the present disclosure.

[0070] The present disclosure provides a method and system for building a thermal model of a power lithium-ion battery based on an electrochemical mechanism, to verify accuracy of a heat generation model of a battery by estimating a temperature of a surface of the battery, thereby improving reliability and safety of a battery pack.

[0071] To make the above-mentioned objective, features, and advantages of the present disclosure clearer and more comprehensible, the present disclosure will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0072] An embodiment of the present disclosure provides a method for building a thermal model of a power lithium-ion battery based on an electrochemical mechanism. Referring to FIG. 1 and FIG. 2, the method specifically includes the following steps.

[0073] Step S1: Discretize a second order partial differential heat conduction equation of a power lithium-ion battery according to a finite differential method, thereby building a thermal model of the power lithium-ion battery. Step S1 specifically includes the following steps.

[0074] Discretize a second order partial differential import equation of a cylindrical power lithium-ion battery according to the finite differential method, thereby building a thermal model of a one-dimensional state space of the cylindrical power lithium-ion battery. For details, referring to FIG. 3 and FIG. 4. A left diagram in FIG. 4 is a top view of an internal structural of a cylindrical power lithium-ion battery; and a right diagram in FIG. 4 is a front view of the internal structural of the cylindrical power lithium-ion battery.

[0075] Assume that a temperature distribution of the cylindrical power lithium-ion battery complies with the following one-dimensional unsteady-state heat conduction equation of cylindrical coordinates:

[00021] $ρ c \frac{\partial T}{\partial t} = \frac{λ}{r} \frac{\partial T}{\partial r} + λ \frac{\partial^{2} T}{\partial r^{2}} + q, namely, ρ c \frac{\partial T}{\partial t} = \frac{λ}{r} \frac{\partial T}{\partial r} + \frac{\partial^{2} T}{\partial^{2} r} + q (t) .$

T denotes a temperature of the battery; t denotes a current moment; r denotes an electrode radius of a layer; q denotes a heat generation rate per unit volume; and q(t) denotes a heat generation rate per unit volume at the current moment.

[0076] Build a boundary condition:

[00022] ${\begin{matrix} - λ \frac{\partial T}{\partial r} {.Math.}_{r_{0}} = h_{0} (T_{0} - T_{1}) = 0 \\ - λ \frac{\partial T}{\partial r} {.Math.}_{R} = h (T_{R} - T_{amb}) \end{matrix} .$

r.sub.0 denotes a radius of a hollow portion of the battery; R denotes a diameter of the surface of the battery; and T.sub.amb denotes an ambient temperature.

[0077] Build an initial condition: T(t)=T.sub.amb.

[0078] Build a supplementary condition: T.sub.1(t)≈T.sub.0(t).

T.sub.0 denotes a temperature of a thin air layer that is of a hollow portion of the battery and that is closest to an inner-most layer of the battery; T.sub.1 and T.sub.R respectively denote a temperature of the inner-most layer of the battery and a temperature of a surface of the battery; h.sub.0 and h respectively denote a convection-diffusion coefficient of the thin air layer and that of the surface of the battery;

[00023] $a = \frac{λ}{ρ c}$

denotes a thermal diffusivity; ρ denotes a density of the power lithium-ion battery; c denotes a specific heat capacity J/(kg.Math.° C.) of the power lithium-ion battery; and λ denotes a radial thermal conductivity.

[0079] Respectively approximate a first-order partial differential equation and a second-order partial differential equation by using a backward difference method and a central difference method, thereby discretizing the second order partial differential heat conduction equation:

[00024] $\frac{\partial T}{\partial r} \approx \frac{T_{k + 1} - T_{k - 1}}{2 Δ r}, \frac{\partial^{2} r}{\partial r^{2}} \approx \frac{T_{k + 1} - 2 T_{k} + T_{k - 1}}{{(Δ r)}^{2}} .$

[0080] Represent the one-dimensional unsteady-state heat conduction equation of the cylindrical power lithium-ion battery as follows:

[00025] ${\dot{T}}_{k} = \frac{a_{k}}{r_{k}} \frac{T_{k + 1} - T_{k - 1}}{2 Δ r} + a_{k} \frac{T_{k + 1} - 2 T_{k} + T_{k - 1}}{{(Δ r)}^{2}} + \frac{a_{k}}{λ} q (t), k = 1, 2, 3 .Math., M - 1. Δ r = \frac{R - r_{0}}{M}; r_{k} = r_{0} + k Δ r;$

M denotes a quantity of walls of the cylindrical power lithium-ion battery; and r.sub.M=R.

[0081] When

[00026] $b_{k} = \frac{a_{k}}{2 r_{k} Δ r} and d_{k} = \frac{2 a_{k}}{{(Δ r)}^{2}},$

[00027] ${\dot{T}}_{k} = (b_{k} + d_{k}) T_{k + 1} - 2 d_{k} T_{k} + (- b_{k} + d_{k}) T_{k - 1} + \frac{a_{k}}{λ} q (t) .$

[0082] Represent the temperature T.sub.M of the surface of the battery as follows:

[00028] $T_{k = M} = \frac{λ_{M}}{h Δ r_{M} + λ_{M}} T_{k = M - 1} + \frac{h Δ r_{M}}{h Δ r_{M} + λ_{M}} T_{amb} .$

[0083] q(t)=q.sub.p+q.sub.r; the equation is solved in step S1.1.4; T.sub.M=T.sub.R; T.sub.M denotes a temperature of an outer layer of the battery; λ.sub.M denotes a thermal conductivity of an outer-layer material; and Δr.sub.M denotes a thickness of the outer-layer material.

[0084] According to a polarization phenomenon and a current heating effect in charging and discharging processes of the battery, build the following calculation formula for a heat release rate of polarization heat and ohmic heat in a use process of the battery:

q.sub.p=I.sup.2R.sub.act+I.sup.2R.sub.ohm=I.sup.2R.sub.t.

I denotes a current of the battery; R.sub.act denotes a polarization internal resistance of the battery; R.sub.ohm denotes an ohmic internal resistance of the battery; R.sub.t denotes a total internal resistance of the battery; and q.sub.p denotes a heat release rate of polarization heat and ohmic heat in a discharging process of the battery.

[0085] Combine calculation steps to acquire the following total heating power:

q=q.sub.r+q.sub.p.

q denotes a total heating power; q.sub.r denotes a heat release rate of reaction heat of the battery; and q.sub.p denotes the heat release rate of the polarization heat and ohmic heat in the discharging process of the battery.

[0086] Determine a final model output according to an actual design demand for a system; carry out approximation according to the finite differential method; iteratively update an electrochemical parameter of the battery according to an Arrhenius equation; and calculate a system output temperature at a current moment.

[0087] Define a system state equation: {dot over (x)}=Ax+Bu.

[0088] Define a system output equation: y=Cx+Du.

u denotes a system input; both C and D denote system matrices; a.sub.1 denotes a coefficient of a thermal diffusivity of a first inner layer of the cylinder; λ.sub.1 denotes a thermal conductivity of the first inner layer of the cylinder; and T denotes a temperature of the battery.

[0089] A and B denote system matrices; a system state is denoted by x=(T.sub.1, T.sub.2, . . . , T.sub.i, . . . , T.sub.M-1).sup.T; a system input is denoted by q and is a heat generation rate per unit volume; and system output y denotes a temperature of an (M−1).sup.th layer of the battery.

[0090] Specifically, each system matrix is denoted as follows:

[00029] $A = (- (b_{1} + d_{1}), (b_{1} + d_{1}), 0, .Math., 0; - (b_{2} + d_{2}), - 2 d_{2}, (b_{2} + d_{2}), .Math., 0; .Math.; 0, 0, 0, .Math., - (b_{M - 1} + d_{M - 1}), (\frac{λ_{M - 1} (b_{M - 1} + d_{M - 1})}{h Δ r + λ_{M - 1}} - 2 d_{M - 1})); B (i, 1) = ρ_{i} c_{i}; C = {(0, 0, .Math., 0, 1)}^{T}; and D = 0.$

custom-character denotes an axial thermal conductivity of an (M−1).sup.th layer; |B(i,1) denotes a system matrix expression of an i.sup.th layer; i denotes an expression of a layer number; ρ.sub.i denotes a density of a battery material of the i.sup.th layer; c.sub.i| denotes a battery specific heat capacity of the i.sup.th layer; and a temperature denoted by system output y may be determined by adjusting a location of 1 in matrix C.

[0091] Secondly, an impact of a temperature on an electrochemical parameter of the battery is determined according to an Arrhenius equation, thereby establishing a coupling relationship of the temperature to the electrochemical parameter of the battery.

[0092] An impact of a temperature on a parameter of the battery is determined based on the relationship between a temperature and an electrochemical parameter of the battery according to the following Arrhenius equation:

[00030] $ψ = ψ_{ref} \exp [\frac{E_{a}^{ψ}}{R} (\frac{1}{T_{ref}} - \frac{1}{T})] .$

ψ.sub.ref is a general variable denoting a diffusion coefficient of a substance, an electrical conductivity of an electrolyte, an exchange current density of an electrode reaction electrode reaction, or the like; subscript ref denotes a value at a reference temperature; E.sub.Ω.sup.ψ denotes an activation energy; ψ denotes a diffusion coefficient; and T.sub.ref denotes a reference temperature.

[0093] The thermal model of the one-dimensional state space of the cylindrical power lithium-ion battery is built based on finite difference and discretization.

[0094] Step S2: Carry out a dynamic working condition test by using a cylindrical power lithium-ion battery selected as an object, thereby acquiring experimental data such as a temperature, a current, a voltage, and a temperature of a surface of the battery. Step S2 specifically includes the following steps:

[0095] step S2.1: selecting a power lithium-ion battery to be tested; and adhering thermocouples to the power lithium-ion battery according to a layout solution;

[0096] step S2.2: allowing the battery to stand still in a 25° C. incubator for 2 h;

[0097] step S2.3: charging the battery in a constant-current and constant-voltage manner to a fully charged state, namely, SOC=100%; discharging the battery at a speed of C/3 to SOC=95%; and allowing the battery to stand still for 2 h;

[0098] step S2.4: loading a dynamic working condition (UDDS) to the battery by an appropriate proportion till SOC of the battery decreases to about 5%;

[0099] step S2.5: recording data such as a current, a voltage, an ambient temperature, and a surface temperature in the working condition;

[0100] step S2.6: repeating steps S2.2 to S2.5 at the same ambient temperature; and acquiring test data at the temperature in dynamic working conditions such as FUDS and UDDS; and

[0101] step S2.7: changing the temperature of the incubator to 5° C., 10° C., and 35° C.; repeating steps S2.2 to S2.6; and acquiring test data at 5° C., 10° C. and 35° C. in the dynamic working conditions.

[0102] Step S3: Identify an electrochemical parameter of the power lithium-ion battery according to an optimal parameter algorithm by using test data acquired in a dynamic working condition, thereby building a thermal model of the power lithium-ion battery, where the optimal parameter algorithm in step S3 is a least square method. This step is not limited to the optimal parameter algorithm.

[0103] Step S4: Verify accuracy of the thermal model of the power lithium-ion battery by using test data acquired in another dynamic working condition.

[0104] In practical application, the temperature of each layer of a specific power lithium-ion battery can be obtained by using the verified thermal model of the power lithium-ion battery. It is determined whether the power lithium-ion battery operates in a safety state by comparing the obtained temperature with a predetermined threshold. If the power lithium-ion battery does not operate in the safety state, a warning such as audio and text warnings is issued, to instruct the user to turn off a switch of the power lithium-ion battery. Alternatively, a controlling signal can be sent to a switching control circuit of the power lithium-ion battery, which can turn off the power lithium-ion battery in response to the controlling signal.

[0105] The thermal model of the power lithium-ion battery according to embodiments of the present disclosure can be carried on hardware devices such as computers or embedded in Battery Management System (BMS) for implementation. BMS is an important component for connecting on-board battery and battery electric vehicle. Functions of the BMS can include: voltage collection of a battery, temperature collection of a battery, current detection of a battery pack, SOC measurement of a battery/a battery pack, assessment of state of health (SOH) of a battery pack, insulation detection and leakage protection, thermal management control and communication, key data recording, fault analysis on a battery and online alarm, etc. When the thermal model of the power lithium-ion battery is embedded in the BMS, the BMS acquires the temperature of each layer of the power lithium-ion battery in real time, and determine whether the power lithium-ion battery operates in a safety state by comparing the key data with the predetermined threshold. If the power lithium-ion battery does not operate in the safety state, an audible alarm can be alarmed, or a text prompt can be displayed on a display, to instruct the user to perform relevant operations or turn off the switch of the power lithium-ion battery.

[0106] In conclusion, an embodiment of the present disclosure provides a method for building a thermal model of a power lithium-ion battery based on an electrochemical mechanism. According to this method, a heat generation rate of a battery is acquired by analyzing an electrochemical mechanism of the battery; and accuracy of a heat generation model of the battery is verified by estimating a temperature of a surface of the battery. Therefore, convenience is brought for state calculation and fault diagnosis of a BMS, which facilitates implementation of a function of a thermal management system for the battery, thereby improving reliability and safety of a battery pack.

[0107] Those skilled in the art are aware that in addition to being realized by using pure computer-readable program code, the system and each apparatus, module, and unit thereof provided in the present disclosure can realize a same program in a form of a logic gate, a switch, an application-specific integrated circuit, a programmable logic controller, or an embedded microcontroller by performing logic programming on the method steps. Therefore, the system and each apparatus, module, and unit thereof provided in the present disclosure can be regarded as a kind of hardware component. The apparatus, module, and unit included therein for realizing each function can also be regarded as a structure in the hardware component; and the apparatus, module, and unit for realizing each function can also be regarded as a software module for implementing the method or a structure in the hardware component.

[0108] The specific embodiments of the present disclosure are described above. It should be understood that the present disclosure is not limited to the above specific implementations, and a person skilled in the art can make various variations or modifications within the scope of the claims without affecting the essence of the present disclosure. The embodiments in the present disclosure and features in the embodiments may be freely combined with each other in a non-conflicting manner.

[0109] Each embodiment of the present specification is described in a progressive manner, each embodiment focuses on the difference from other embodiments, and the same and similar parts between the examples may refer to each other.

[0110] In this specification, some specific embodiments are used for illustration of the principles and implementations of the present disclosure. The description of the foregoing embodiments is used to help illustrate the method of the present disclosure and the core ideas thereof. In addition, those of ordinary skill in the art can make various modifications in terms of specific implementations and the scope of application in accordance with the ideas of the present disclosure. In conclusion, the content of the present description shall not be construed as limitations to the present disclosure.

METHOD AND SYSTEM FOR BUILDING THERMAL MODEL OF POWER LITHIUM-ION BATTERY BASED ON ELECTROCHEMICAL MECHANISM

Inventors

Cpc classification

Classification Explorer

G06F17/13

PHYSICS

Classification Explorer

G06F2119/08

PHYSICS

Classification Explorer

Y02E60/10

GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS

Classification Explorer

H01M10/486

ELECTRICITY

Classification Explorer

G06F30/23

PHYSICS

Classification Explorer

G01R31/367

PHYSICS

Classification Explorer

G06F17/11

PHYSICS

International classification

Classification Explorer

G06F30/23

PHYSICS

Classification Explorer

G06F17/11

PHYSICS

Abstract

Claims

Description