ASRUKF-based multi-point internal temperature online estimation method for full-life wide-temperature square lithium battery

Abstract

The invention relates to the technical field of lithium ion battery state estimation, and discloses an ASRUKF-based multi-point internal temperature online estimation method for a full-life wide-temperature square lithium battery. According to the method, a lithium battery electrothermal coupling model is constructed based on a second-order RC equivalent circuit model and a three-heat-source heat production model in combination with a two-state lumped heat transfer model; model parameters are identified, so that precise modeling of internal temperature estimation of the lithium ion battery under the full-life wide temperature is realized; establishing a multi-point temperature heat transfer model, discretizing representative region temperature to obtain a system state space equation and a measurement equation, and realizing multi-point temperature estimation of the square battery; and the internal temperatures of different points of the square lithium battery are estimated by fusing an adaptive square root unscented Kalman filtering algorithm, and the estimated internal temperatures are used as feedback quantity to update circuit model parameters, so that coupling of an electric heating model and online estimation of the internal temperatures of multiple points are realized.

Description

Technical Field

[0001] The invention relates to the technical field of lithium-ion battery state estimation, and in particular to an ASRUKF-based full-life wide-temperature square lithium battery multi-point internal temperature online estimation method. Background Art

[0002] With the development of electric vehicle technology and the expansion of the market, the safety and reliability of lithium-ion batteries as the internal power components of electric vehicles have attracted much attention. As a power source, lithium batteries often need to work under complex driving conditions, and the battery temperature state changes accordingly. Accurately estimating the internal temperature of lithium batteries under different working environments is crucial to ensure the safe use of batteries.

[0003] In practical applications, the internal temperature state of battery cells can only be obtained through built-in sensors, but installing sensors in batteries used in groups on vehicles will undoubtedly increase costs and the risk of loss of control. At present, the research approaches for estimating the internal state of batteries mainly include: establishing a finite element model based on simulation software to simulate and analyze the heat generation and heat transfer mechanism inside lithium batteries. This method cannot fully consider the influence of various external factors in practical applications, and the results are too idealized; although the data-driven research method is more realistic, it requires a large amount of historical data to support model training; although the partial differential equation based on the chemical reaction mechanism of lithium batteries can fully reflect the reaction kinetics behavior inside the battery, it requires a large amount of partial differential calculations, has strong parameter dependence, and the model is complex and has a large amount of calculation; based on the resistance and capacitance characteristics of lithium batteries, an equivalent circuit model is established, and the internal charge and discharge characteristics of lithium batteries are simulated by a combination of electronic components, and the internal temperature of the battery is estimated by coupling with the thermal model. The physical meaning is clear, the model is simple, and the result is highly accurate. Among them, the square battery has a large area, and the heating conditions in different areas inside the battery are different. Excessive local temperature may cause thermal runaway of the battery starting from a point.

[0004] At present, there is still a lack of research on multi-point temperature estimation of square batteries. It is very meaningful to establish an electrothermal coupling model that conforms to realistic application scenarios and covers a wide temperature operating range within the entire cycle life cycle of the battery to accurately estimate the internal temperature of the battery and enable online application. Summary of the invention

[0005] 1. Technical issues to be resolved

[0006] In view of the shortcomings of the prior art, the present invention provides an online estimation method for multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF, which has the advantages of small calculation amount, high precision and strong robustness, and solves the above technical problems.

[0007] (II) Technical solution

[0008] To achieve the above object, the present invention provides the following technical solution: an online estimation method for multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF, comprising the following steps:

[0009] S1: Obtain the charge and discharge data of lithium-ion single cells under full life and wide temperature conditions. The specific steps are as follows: perform cyclic charge and discharge aging tests on different batteries of the same model, fully charge the battery using constant current and constant voltage, let it stand for 30 seconds, and then discharge it to the battery cut-off voltage. Repeat the above process, discharge the battery to 95%, 85%, and 80% of the nominal capacity, respectively, and then perform working condition experiments within the temperature range of -20℃ to 60℃, and record the battery's current, voltage, representative area, and lug temperature;

[0010] S2. Establish an equivalent circuit model and a thermal model, wherein the equivalent circuit model is a second-order RC equivalent circuit model, and the battery thermal model includes a heat generation model and a heat transfer model. The heat generation model uses a simplified three-heat source heat generation model to represent the heat generation inside the battery, and uses a two-state lumped parameter thermal model as a heat transfer model to represent the heat transfer inside the battery. The second-order RC equivalent circuit model, the simplified three-heat source heat generation model, and the two-state lumped parameter thermal model are discretized. After discretization, a lithium-ion battery internal temperature estimation model is established to obtain the state space equation and measurement equation of the system;

[0011] S3. Model parameter identification: Using the experimental data under HPPC working conditions, the relationship curve between SOC and OCV is obtained through 6th-order polynomial fitting. The fitting expression is as follows:

[0012]

[0013] in, represents the polynomial fitting parameters, represents the sum of the entire sixth-order polynomial, express Different powers of ;

[0014] Under different aging conditions at a wide temperature, based on the battery charge and discharge experimental data obtained by S1, the FFRLS algorithm is used to perform parameter identification on the second-order RC equivalent circuit and parameter thermal model established in S2. The identification process is as follows:

[0015] Calculation error:

[0016]

[0017] Calculate the gain vector:

[0018]

[0019] Update the parameter vector:

[0020]

[0021] Update the covariance matrix:

[0022]

[0023] in, represents the prediction error at the current moment, Represents the system output at the current moment, is the system regression vector at the current moment, for The transposition of represents the covariance at the previous moment, is the parameter vector at the current moment, represents the parameter vector at the previous moment, is the gain vector at the current moment, is the system covariance matrix at the current moment, For the forgetting factor;

[0024] S4, online estimation of the internal temperature of different points of the square lithium battery based on the ASRUKF algorithm: Based on the system state space equation and measurement equation obtained in S2, the terminal voltage and open circuit voltage output by the second-order RC equivalent circuit model are substituted into the three-heat source heat generation model to calculate the battery heat generation power, and then the two-state lumped parameter thermal model is combined with the ASRUKF algorithm to estimate the internal temperature of multiple points of the lithium battery;

[0025] S5. Coupling of equivalent circuit model and thermal model and updating of model parameters: RC parameters in the second-order RC equivalent circuit model , , , , Thermal parameters in heat transfer models , It will change with the temperature of the battery, and update the parameters identified based on S3 under different temperature gradients.

[0026] As a preferred technical solution of the present invention, the expression of the second-order RC equivalent circuit model in step S2 is as follows:

[0027]

[0028]

[0029] in, express The open circuit voltage at time express The terminal voltage at the moment, express The current at the moment, , express The polarization voltage at time Represents the battery's ohmic internal resistance, , They represent polarization resistance, , They represent polarization capacitance, and They represent the time constants, represents a natural constant, represents the sampling time interval, and the second-order RC equivalent circuit model is discretized using the Euler method, and the expression obtained is as follows:

[0030]

[0031] in, is the current battery open circuit voltage, is the battery terminal voltage at the current moment, is the battery current at the current moment, Indicates the battery current at the previous moment, The battery's ohmic internal resistance, , They are the polarization resistance at the current moment , and polarization capacitance , The corresponding terminal voltage, , are the polarization resistance at the previous moment , and polarization capacitance , The corresponding terminal voltage, and is the time constant, Represents the current state quantity, Indicates the state quantity at the previous moment;

[0032] The simplified three-heat source heat generation model is expressed as follows:

[0033]

[0034] in, express The heat generation power of the battery cell at the moment, express The open circuit voltage at time express The terminal voltage at the moment, express The current at the moment, , Respectively The heat generation power of the two tabs at any given moment, express The internal temperature of the battery at all times, express right The partial guide, , , Represent the internal resistance of the battery cell, positive electrode and negative electrode respectively. Discretize the simplified three-heat source heat generation model, and the expression obtained is as follows:

[0035]

[0036] in, Indicates the heat generation power of the battery cell at the current moment. is the battery current at the current moment, is the current battery open circuit voltage, is the battery terminal voltage at the current moment, , Respectively represent the heat generation power of the two tabs at the current moment;

[0037] The two-state lumped parameter thermal model expression is as follows:

[0038]

[0039] in, , , , Respectively The internal temperature of the battery at four measurement points at the moment, Represents the internal resistance between the battery surface and the interior, Indicates the internal resistance of the battery surface shell, express The temperature of the positive electrode ear at all times, express The temperature of the negative electrode ear at the moment, , , , Respectively The surface temperature of the four measurement points on the battery surface at the moment, the expression after discretization of the two-state lumped parameter thermal model is as follows:

[0040]

[0041] in, , , , They represent the estimated internal temperature of the four measurement points of the battery at the current moment. , , , Respectively represent the surface temperature of the four measurement points of the battery at the current moment, Represents the internal resistance between the battery surface and the interior, Indicates the internal resistance of the battery surface shell;

[0042] The internal temperature estimation model of lithium-ion batteries is established, and the state space equation and measurement equation of the system are obtained as follows:

[0043]

[0044] in, represents the current state vector, represents the state transfer matrix; is the system control matrix; is the control quantity at the previous moment; Measure the state vector at the current moment; is the measurement noise matrix at the previous moment White Gaussian measurement noise with covariance Represents the transformation matrix from state quantity to observation quantity.

[0045] As a preferred technical solution of the present invention, the expression of the state transfer matrix is ​​as follows:

[0046]

[0047] The expression of the current state vector is as follows:

[0048]

[0049] The expression of the system control matrix is ​​as follows:

[0050]

[0051] The expression of the control quantity at the current moment is as follows:

[0052]

[0053] in, , , , Respectively The internal temperature of the battery at four measurement points at the moment, Represents the internal resistance between the battery surface and the interior, Indicates the internal resistance of the battery surface shell, express The temperature of the positive electrode ear at all times, express The temperature of the negative electrode ear at the moment, , , , Respectively The surface temperature of the four measurement points on the battery surface at the moment, , They respectively represent the heat generation power of the two pole ears at the current moment.

[0054] As a preferred technical solution of the present invention, step S4 specifically includes the following steps:

[0055] S4.1. Initialization of state function. The specific expression is as follows:

[0056]

[0057] in, represents the state vector at time 0, represents the expected value of the state vector, represents the mathematical expectation of the calculated parameters, represents the initial value of the error covariance matrix, the superscript represents transpose, Express Perform Cholesky decomposition, represents the initial value of the square root, represents the initial value of process noise, represents the initial value of the measurement noise, , They represent the process noise and measurement noise at time 0 respectively;

[0058] S4.2, Sigma point weight selection, the specific expression is as follows:

[0059]

[0060] in, represents the scale factor, represents the mean weight of the 0th Sigma point, represents the covariance weight of the 0th Sigma point, Indicates the serial number of other Sigma points, Represents the mean weight of other Sigma points, represents the covariance weight of other Sigma points, Represents the state parameter, represents the state variable dimension, represents the secondary scaling parameter;

[0061] The expression for constructing the Sigma point is as follows:

[0062]

[0063] in, represents the state vector at the previous moment, Represents the square root value at the current moment, represents the scale factor, represents the vector formed by the Sigma points at the previous moment, Represents the state variable dimension;

[0064] S4.3, time update, the expression for calculating the state estimate and the square root error covariance is as follows:

[0065]

[0066] in, Represents the Sigma point of the current state estimate, represents the current estimated covariance matrix, Express Perform QR decomposition, represents the state transition function of the system, represents the updated covariance matrix, express The function is used to update the covariance matrix. represents the process noise at the current moment;

[0067] S4.4, Sigma point resampling, the expression is as follows:

[0068]

[0069] Predict the observed values ​​of the system:

[0070]

[0071] in, represents the estimated observed value, Represents the observation function, and the state Sigma point Mapped to the observation space, represents the current state estimate, Represents the residual at the current moment;

[0072] Estimation of measurement noise:

[0073]

[0074] in, represents the measurement noise covariance matrix, represents the updated measurement noise covariance matrix, represents the further updated measurement noise covariance matrix, represents the final updated measurement noise covariance matrix, represents the scalar gain factor, , represents the forgetting factor, represents the forgetting factor at the next moment, represents the predicted observation value vector obtained through Sigma point propagation, Indicates the construction of a diagonal matrix;

[0075] Updates to the variance matrix and error covariance:

[0076]

[0077] in, is the updated observation error covariance matrix, is the further updated observation error covariance matrix, is the cross covariance matrix between the state quantity and the observation quantity;

[0078] Calculate the Kalman gain matrix and update the estimator and variance:

[0079]

[0080] Update of process noise:

[0081]

[0082] in, represents the process noise covariance matrix at the previous moment, represents the updated process noise covariance matrix, represents the further updated process noise covariance matrix, Represents the final updated process noise covariance matrix.

[0083] Compared with the prior art, the present invention provides an online estimation method for multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF, which has the following beneficial effects:

[0084] 1. The present invention adopts a three-heat source heat generation model through the heat generation model, which characterizes the heat generation of the tab and the inside of the battery respectively, and can fully describe the different characteristics of heat generation in different areas of the square lithium battery during the charging and discharging process. It is coupled with the second-order RC equivalent circuit model to accurately describe the resistance-capacitance characteristics and heat generation characteristics of the battery. The model has low calculation complexity, which can improve the calculation efficiency while ensuring accuracy. The ASRUKF algorithm is integrated to estimate the multi-point internal temperature of the square lithium battery. The adaptive mechanism and square root filtering method are combined to achieve more stable temperature estimation in high noise or nonlinear systems.

[0085] 2. The present invention considers the influence of external temperature on the parameters of the electrothermal model under different aging degrees of the battery, adopts battery circuit parameters that match the current temperature state of the battery under full cycle life to calculate the terminal voltage and open circuit voltage, considers the influence of different working environment temperatures of the lithium battery on the model parameters, uses the currently estimated internal temperature as feedback to update the circuit model parameters online, and dynamically adjusts the model parameters to adapt to the changes in the actual working environment temperature, thereby improving the accuracy and real-time performance of temperature estimation. Considering the large area of ​​square lithium batteries and the uneven internal heat generation, combined with experimental analysis, representative heating areas are modeled, and a heat generation and heat transfer model from the tabs to the inside of the battery is established to achieve online estimation of the multi-point internal temperature of square lithium batteries. BRIEF DESCRIPTION OF THE DRAWINGS

[0086] Figure 1 It is a schematic diagram of the process of the present invention;

[0087] Figure 2 It is a schematic diagram of the second-order RC equivalent circuit model of the present invention;

[0088] Figure 3 This is a schematic diagram of the three-heat source heat generation model of the present invention. DETAILED DESCRIPTION

[0089] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0090] See also Figure 1-3 , an online estimation method for multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF, comprising the following steps:

[0091] S1: Obtain the charge and discharge data of lithium-ion single cells under full life and wide temperature conditions. The specific steps are: perform cyclic charge and discharge aging tests on different batteries of the same model, fully charge the battery using constant current and constant voltage, let it stand for 30 seconds, and then discharge it to the battery cut-off voltage. Repeat the above process, discharge the battery to 95%, 85%, and 80% of the nominal capacity, respectively, and then perform working condition experiments within the temperature range of -20℃ to 60℃, record the battery's current, voltage, representative area, and lug temperature. The experimental data are recorded as follows:

[0092] Table 1 Experimental record data

[0093]

[0094] S2. Establish an equivalent circuit model and a thermal model, wherein the equivalent circuit model is a second-order RC equivalent circuit model, and the battery thermal model includes a heat generation model and a heat transfer model. The heat generation model uses a simplified three-heat source heat generation model to represent the heat generation inside the battery, and uses a two-state lumped parameter thermal model as a heat transfer model to represent the heat transfer inside the battery. The second-order RC equivalent circuit model, the simplified three-heat source heat generation model, and the two-state lumped parameter thermal model are discretized. After discretization, a lithium-ion battery internal temperature estimation model is established to obtain the state space equation and measurement equation of the system. The expression of the second-order RC equivalent circuit model in step S2 is as follows:

[0095]

[0096]

[0097] in, express The open circuit voltage at time express The terminal voltage at the moment, express The current at the moment, , express The polarization voltage at time Represents the battery's ohmic internal resistance, , They represent polarization resistance, , They represent polarization capacitance, and They represent the time constants, represents a natural constant, represents the sampling time interval, and the second-order RC equivalent circuit model is discretized using the Euler method, and the expression obtained is as follows:

[0098]

[0099] in, is the current battery open circuit voltage, is the battery terminal voltage at the current moment, is the battery current at the current moment, Indicates the battery current at the previous moment, The battery's ohmic internal resistance, , They are the polarization resistance at the current moment , and polarization capacitance , The corresponding terminal voltage, , are the polarization resistance at the previous moment , and polarization capacitance , The corresponding terminal voltage, and is the time constant, Represents the current state quantity, Indicates the state quantity at the previous moment;

[0100] The simplified three-heat source heat generation model is expressed as follows:

[0101]

[0102] in, express The heat generation power of the battery cell at the moment, express The open circuit voltage at time express The terminal voltage at the moment, express The current at the moment, , Respectively The heat generation power of the two tabs at any given moment, express The internal temperature of the battery at all times, express right The partial guide, , , Represent the internal resistance of the battery cell, positive electrode and negative electrode respectively. Discretize the simplified three-heat source heat generation model, and the expression obtained is as follows:

[0103]

[0104] in, Indicates the heat generation power of the battery cell at the current moment. is the battery current at the current moment, is the current battery open circuit voltage, is the battery terminal voltage at the current moment, , Respectively represent the heat generation power of the two tabs at the current moment;

[0105] The two-state lumped parameter thermal model expression is as follows:

[0106]

[0107] in, , , , Respectively The internal temperature of the battery at four measurement points at the moment, Represents the internal resistance between the battery surface and the interior, Indicates the internal resistance of the battery surface shell, express The temperature of the positive electrode ear at all times, express The temperature of the negative electrode ear at the moment, , , , Respectively The surface temperature of the four measurement points on the battery surface at the moment, the expression after discretization of the two-state lumped parameter thermal model is as follows:

[0108]

[0109] in, , , , They represent the estimated internal temperature of the four measurement points of the battery at the current moment. , , , Respectively represent the surface temperature of the four measurement points of the battery at the current moment, Represents the internal resistance between the battery surface and the interior, Indicates the internal resistance of the battery surface shell;

[0110] The internal temperature estimation model of lithium-ion batteries is established, and the state space equation and measurement equation of the system are obtained as follows:

[0111]

[0112] in, represents the current state vector, represents the state transfer matrix; is the system control matrix; is the control quantity at the previous moment; Measure the state vector at the current moment; is the measurement noise matrix at the previous moment White Gaussian measurement noise with covariance Represents the conversion matrix from state quantity to observation quantity;

[0113] The expression of the state transfer matrix is ​​as follows:

[0114]

[0115] The expression of the current state vector is as follows:

[0116]

[0117] The expression of the system control matrix is ​​as follows:

[0118]

[0119] The expression of the control quantity at the current moment is as follows:

[0120]

[0121] in, , , , Respectively The internal temperature of the battery at four measurement points at the moment, Represents the internal resistance between the battery surface and the interior, Indicates the internal resistance of the battery surface shell, express The temperature of the positive electrode ear at all times, express The temperature of the negative electrode ear at the moment, , , , Respectively The surface temperature of the four measurement points on the battery surface at the moment, , Respectively represent the heat generation power of the two tabs at the current moment;

[0122] S3. Model parameter identification: Using the experimental data under HPPC working conditions, the relationship curve between SOC and OCV is obtained through sixth-order polynomial fitting. The fitting expression is as follows:

[0123]

[0124] in, represents the polynomial fitting parameters, represents the sum of the entire sixth-order polynomial, express Different powers of ;

[0125] Under different aging conditions at a wide temperature, based on the battery charge and discharge experimental data obtained by S1, the FFRLS algorithm is used to perform parameter identification on the second-order RC equivalent circuit and parameter thermal model established in S2. The identification process is as follows:

[0126] Calculation error:

[0127]

[0128] Calculate the gain vector:

[0129]

[0130] Update the parameter vector:

[0131]

[0132] Update the covariance matrix:

[0133]

[0134] in, represents the prediction error at the current moment, Represents the system output at the current moment, is the system regression vector at the current moment, for The transposition of represents the covariance at the previous moment, is the parameter vector at the current moment, represents the parameter vector at the previous moment, is the gain vector at the current moment, is the system covariance matrix at the current moment, For the forgetting factor;

[0135] S4, online estimation of the internal temperature of different points of the square lithium battery based on the ASRUKF algorithm: Based on the system state space equation and measurement equation obtained in S2, the terminal voltage and open circuit voltage output by the second-order RC equivalent circuit model are substituted into the three-heat source heat generation model to calculate the battery heat generation power, and then the two-state lumped parameter thermal model is combined with the ASRUKF algorithm to estimate the internal temperature of multiple points of the lithium battery;

[0136] Step S4 specifically includes the following steps:

[0137] S4.1. Initialization of state function. The specific expression is as follows:

[0138]

[0139] in, represents the state vector at time 0, represents the expected value of the state vector, represents the mathematical expectation of the calculated parameters, represents the initial value of the error covariance matrix, the superscript represents transpose, Express Perform Cholesky decomposition, represents the initial value of the square root, represents the initial value of process noise, represents the initial value of the measurement noise, , They represent the process noise and measurement noise at time 0 respectively;

[0140] S4.2, Sigma point weight selection, the specific expression is as follows:

[0141]

[0142] in, represents the scale factor, represents the mean weight of the 0th Sigma point, represents the covariance weight of the 0th Sigma point, Indicates the serial number of other Sigma points, Represents the mean weight of other Sigma points, represents the covariance weight of other Sigma points, Represents the state parameter, represents the state variable dimension, represents the secondary scaling parameter;

[0143] The expression for constructing the Sigma point is as follows:

[0144]

[0145] in, represents the state vector at the previous moment, Represents the square root value at the current moment, represents the scale factor, represents the vector formed by the Sigma points at the previous moment, Represents the state variable dimension;

[0146] S4.3, time update, the expression for calculating the state estimate and the square root error covariance is as follows:

[0147]

[0148] in, Represents the Sigma point of the current state estimate, represents the current estimated covariance matrix, Express Perform QR decomposition, represents the state transition function of the system, represents the updated covariance matrix, express The function is used to update the covariance matrix. represents the process noise at the current moment;

[0149] S4.4, Sigma point resampling, the expression is as follows:

[0150]

[0151] Predict the observed values ​​of the system:

[0152]

[0153] in, represents the estimated observed value, Represents the observation function, and the state Sigma point Mapped to the observation space, represents the current state estimate, Represents the residual at the current moment;

[0154] Estimation of measurement noise:

[0155]

[0156] in, represents the measurement noise covariance matrix, represents the updated measurement noise covariance matrix, represents the further updated measurement noise covariance matrix, represents the final updated measurement noise covariance matrix, represents the scalar gain factor, , represents the forgetting factor, represents the forgetting factor at the next moment, represents the predicted observation value vector obtained through Sigma point propagation, Indicates the construction of a diagonal matrix;

[0157] Updates to the variance matrix and error covariance:

[0158]

[0159] in, is the updated observation error covariance matrix, is the further updated observation error covariance matrix, is the cross covariance matrix between the state quantity and the observation quantity;

[0160] Calculate the Kalman gain matrix and update the estimator and variance:

[0161]

[0162] Update of process noise:

[0163]

[0164] in, represents the process noise covariance matrix at the previous moment, represents the updated process noise covariance matrix, represents the further updated process noise covariance matrix, represents the final updated process noise covariance matrix;

[0165] S5. Coupling of equivalent circuit model and thermal model and updating of model parameters. During operation, the internal temperature of the battery will change with the change of current and voltage, and the physical parameters of the battery will also change accordingly. The resistance and capacitance parameters in the second-order RC equivalent circuit model , , , , Thermal parameters in heat transfer models , It needs to be updated as the temperature changes, so that the current, voltage and other parameters in the current state can be passed to the thermal model to calculate the heat generation power, and the new battery current temperature state T can be estimated. The estimated result of T is fed back to the equivalent circuit model to update the resistance and capacitance parameters, thereby realizing the coupling of the electrical model and the thermal model.

[0166] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

Claims

1. An online estimation method for multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF, characterized by: The following steps are involved: S1: Obtain the charge and discharge data of lithium-ion single cells under full life and wide temperature conditions. The specific steps are as follows: perform cyclic charge and discharge aging tests on different batteries of the same model, fully charge the battery using constant current and constant voltage, let it stand for 30 seconds, and then discharge it to the battery cut-off voltage. Repeat the above process, discharge the battery to 95%, 85%, and 80% of the nominal capacity, respectively, and then perform working condition experiments within the temperature range of -20℃ to 60℃, and record the battery's current, voltage, representative area, and lug temperature; S2. Establish an equivalent circuit model and a thermal model, wherein the equivalent circuit model is a second-order RC equivalent circuit model, and the battery thermal model includes a heat generation model and a heat transfer model. The heat generation model uses a simplified three-heat source heat generation model to represent the heat generation inside the battery, and uses a two-state lumped parameter thermal model as a heat transfer model to represent the heat transfer inside the battery. The second-order RC equivalent circuit model, the simplified three-heat source heat generation model, and the two-state lumped parameter thermal model are discretized. After discretization, a lithium-ion battery internal temperature estimation model is established to obtain the state space equation and measurement equation of the system; S3. Model parameter identification: Using the experimental data under HPPC working conditions, the relationship curve between SOC and OCV is obtained through sixth-order polynomial fitting. The fitting expression is as follows: ； in, represents the polynomial fitting parameters, represents the sum of the entire sixth-order polynomial, express Different powers of ; Under different aging conditions at a wide temperature, based on the battery charge and discharge experimental data obtained by S1, the FFRLS algorithm is used to perform parameter identification on the second-order RC equivalent circuit and parameter thermal model established in S2. The identification process is as follows: Calculation error: ； Calculate the gain vector: ； Update the parameter vector: ； Update the covariance matrix: ； in, represents the prediction error at the current moment, Represents the system output at the current moment, is the system regression vector at the current moment, for The transposition of represents the covariance at the previous moment, is the parameter vector at the current moment, represents the parameter vector at the previous moment, is the gain vector at the current moment, is the system covariance matrix at the current moment, For the forgetting factor; S4, online estimation of the internal temperature of different points of the square lithium battery based on the ASRUKF algorithm: Based on the system state space equation and measurement equation obtained in S2, the terminal voltage and open circuit voltage output by the second-order RC equivalent circuit model are substituted into the three-heat source heat generation model to calculate the battery heat generation power, and then the two-state lumped parameter thermal model is combined with the ASRUKF algorithm to estimate the internal temperature of multiple points of the lithium battery; S5. Coupling of equivalent circuit model and thermal model and updating of model parameters: RC parameters in the second-order RC equivalent circuit model , , , , Thermal parameters in heat transfer models , It will change with the temperature of the battery, and update the parameters identified based on S3 under different temperature gradients.

2. The method for online estimation of multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF according to claim 1 is characterized in that: The expression of the second-order RC equivalent circuit model in step S2 is as follows: ； ； in, express The open circuit voltage at time express The terminal voltage at the moment, express The current at the moment, , express The polarization voltage at time Represents the battery's ohmic internal resistance, , They represent polarization resistance, , They represent polarization capacitance, and They represent the time constants, represents a natural constant, represents the sampling time interval, and the second-order RC equivalent circuit model is discretized using the Euler method, and the expression obtained is as follows: ； in, is the current battery open circuit voltage, is the battery terminal voltage at the current moment, is the battery current at the current moment, Indicates the battery current at the previous moment, The battery's ohmic internal resistance, , They are the polarization resistance at the current moment , and polarization capacitance , The corresponding terminal voltage, , are the polarization resistance at the previous moment , and polarization capacitance , The corresponding terminal voltage, and is the time constant, Represents the current state quantity, Indicates the state quantity at the previous moment; The simplified three-heat source heat generation model is expressed as follows: ； in, express The heat generation power of the battery cell at the moment, express The open circuit voltage at time express The terminal voltage at the moment, express The current at the moment, , Respectively The heat generation power of the two tabs at any given moment, express The internal temperature of the battery at all times, express right The partial guide, , , Represent the internal resistance of the battery cell, positive electrode and negative electrode respectively. Discretize the simplified three-heat source heat generation model, and the expression obtained is as follows: ； in, Indicates the heat generation power of the battery cell at the current moment. is the battery current at the current moment, is the current battery open circuit voltage, is the battery terminal voltage at the current moment, , They respectively represent the heat generation power of the two tabs at the current moment; The two-state lumped parameter thermal model expression is as follows: ； in, , , , Respectively The internal temperature of the battery at four measurement points at the moment, Represents the internal resistance between the battery surface and the interior. Indicates the internal resistance of the battery surface shell, express The temperature of the positive electrode ear at all times, express The temperature of the negative electrode ear at the moment, , , , Respectively The surface temperature of the four measurement points on the battery surface at the moment, the expression after discretization of the two-state lumped parameter thermal model is as follows: ； in, , , , They represent the estimated internal temperature of the four measurement points of the battery at the current moment. , , , Respectively represent the surface temperature of the four measurement points of the battery at the current moment, Represents the internal resistance between the battery surface and the interior. Indicates the internal resistance of the battery surface shell; The internal temperature estimation model of lithium-ion batteries is established, and the state space equation and measurement equation of the system are obtained as follows: ； in, represents the current state vector, represents the state transfer matrix, is the system control matrix, is the control quantity at the previous moment, is the measured state vector at the current moment, is the measurement noise matrix at the previous moment The covariance of the white Gaussian measurement noise, Represents the transformation matrix from state quantity to observation quantity.

3. The method for online estimation of multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF according to claim 2 is characterized in that: The expression of the state transfer matrix is ​​as follows: ； The expression of the current state vector is as follows: ； The expression of the system control matrix is ​​as follows: ； The expression of the control quantity at the current moment is as follows: ； in, , , , Respectively The internal temperature of the battery at four measurement points at the moment, Represents the internal resistance between the battery surface and the interior, Indicates the internal resistance of the battery surface shell, express The temperature of the positive electrode ear at all times, express The temperature of the negative electrode ear at the moment, , , , Respectively The surface temperature of the four measurement points on the battery surface at the moment, , They respectively represent the heat generation power of the two pole ears at the current moment.

4. The method for online estimation of multi-point internal temperature of full-life wide-temperature square lithium battery based on ASRUKF according to claim 1, characterized in that: The step S4 specifically comprises the following steps: S4.

1. Initialization of state function. The specific expression is as follows: ； in, represents the state vector at time 0, represents the expected value of the state vector, represents the mathematical expectation of the calculated parameters, represents the initial value of the error covariance matrix, the superscript represents transpose, Express Perform Cholesky decomposition, represents the initial value of the square root, represents the initial value of process noise, represents the initial value of the measurement noise, , They represent the process noise and measurement noise at time 0 respectively; S4.2, Sigma point weight selection, the specific expression is as follows: ； in, represents the scale factor, represents the mean weight of the 0th Sigma point, represents the covariance weight of the 0th Sigma point, Indicates the serial number of other Sigma points, Represents the mean weight of other Sigma points, represents the covariance weight of other Sigma points, Represents the state parameter, represents the state variable dimension, represents the secondary scaling parameter; The expression for constructing the Sigma point is as follows: ； in, represents the state vector at the previous moment, Represents the square root value at the current moment, represents the scale factor, represents the vector formed by the Sigma points at the previous moment, Represents the state variable dimension; S4.3, time update, the expression for calculating the state estimate and the square root error covariance is as follows: ； in, Represents the Sigma point of the current state estimate, represents the current estimated covariance matrix, Express Perform QR decomposition, represents the state transition function of the system, represents the updated covariance matrix, express The function is used to update the covariance matrix. represents the process noise at the current moment; S4.4, Sigma point resampling, the expression is as follows: ； Predict the observed values ​​of the system: ； in, represents the estimated observed value, Represents the observation function, and the state Sigma point Mapped to the observation space, represents the current state estimate, Represents the residual at the current moment; Estimation of measurement noise: ； in, represents the measurement noise covariance matrix, represents the updated measurement noise covariance matrix, represents the further updated measurement noise covariance matrix, represents the final updated measurement noise covariance matrix, represents the scalar gain factor, , represents the forgetting factor, represents the forgetting factor at the next moment, represents the predicted observation value vector obtained through Sigma point propagation, Indicates the construction of a diagonal matrix; Updates to the variance matrix and error covariance: ； in, is the updated observation error covariance matrix, is the further updated observation error covariance matrix, is the cross covariance matrix between the state quantity and the observation quantity; Calculate the Kalman gain matrix and update the estimator and variance: ； Update of process noise: ； in, represents the process noise covariance matrix at the previous moment, represents the updated process noise covariance matrix, represents the further updated process noise covariance matrix, Represents the final updated process noise covariance matrix.

Images