CN112883638B - Online estimation method for temperature distribution of super capacitor module - Google Patents

Online estimation method for temperature distribution of super capacitor module Download PDF

Info

Publication number
CN112883638B
CN112883638B CN202110140366.9A CN202110140366A CN112883638B CN 112883638 B CN112883638 B CN 112883638B CN 202110140366 A CN202110140366 A CN 202110140366A CN 112883638 B CN112883638 B CN 112883638B
Authority
CN
China
Prior art keywords
temperature
module
monomers
neural network
loss
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.)
Active
Application number
CN202110140366.9A
Other languages
Chinese (zh)
Other versions
CN112883638A (en
Inventor
韦莉
张婉婷
吴铭
闫孟迪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN202110140366.9A priority Critical patent/CN112883638B/en
Publication of CN112883638A publication Critical patent/CN112883638A/en
Application granted granted Critical
Publication of CN112883638B publication Critical patent/CN112883638B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/048Activation functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/13Energy storage using capacitors

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

The invention relates to an online estimation method for temperature distribution of a super capacitor module, which can realize online estimation of the temperature of all monomers in the module by measuring the temperature of a few monomers in the module. The method avoids using a large number of temperature sensors, not only greatly reduces the hardware cost, but also solves the problem that the temperature sensors are not easy to install in partial areas inside the module. Meanwhile, the algorithm occupies less computing resources, is suitable for online estimation, and can meet higher computing precision. The method has important significance for heat management optimization and safe operation of the super capacitor.

Description

Online estimation method for temperature distribution of super capacitor module
Technical Field
The invention relates to an online estimation method for temperature distribution of a super capacitor module.
Background
The super capacitor is used as a novel energy storage device, and has wide application prospect in the fields of renewable clean energy power generation, energy-saving operation of electric automobiles, urban rails, ships and the like. Along with the work of the super capacitor, the accumulation of the internal temperature of the super capacitor can continuously raise the temperature of the monomer, influence the electrical characteristics of the super capacitor, cause the performance to be reduced, accelerate the aging of the super capacitor and further shorten the service life of the super capacitor. The thermal runaway of the whole system can be caused after the temperature of the super capacitor is too high or is out of control, so that the online estimation of the temperature of the super capacitor is particularly important, the management of the super capacitor is more accurate and efficient, and the safety of an energy storage system is improved.
With the wide application of the super capacitor energy storage system, the temperature of each super capacitor monomer needs to be monitored on line in more and more application occasions. However, in practical applications, the temperature sensor is limited by the problems of cost, installation and the like, and often only detects the temperature of a few monomers in the module, and cannot monitor the temperature of each monomer. In addition, finite element thermal simulation software is generally used for researching the temperature distribution of the module, but the method is complex in model, large in calculation amount and long in calculation time, and is not beneficial to online operation.
Disclosure of Invention
Aiming at the defects of the conventional supercapacitor module temperature estimation method, the invention provides the supercapacitor module temperature online estimation method which can effectively utilize a small amount of monomer temperature to estimate all the monomer temperatures in the module and has a simple calculation process. The method combines the characteristics of monomer arrangement in the module, optimizes and selects the installation position of the temperature sensor, considers the correlation among monomers of the super capacitor module and establishes a neural network model for online temperature estimation. The method is realized by the following steps as shown in figure 1:
the method comprises the following steps: the method comprises the following steps of (1) taking the temperature correlation of monomers in a super capacitor module into consideration, and partitioning the module by combining the arrangement characteristics of the monomers;
step two: selecting a tested monomer according to the subareas for testing to obtain temperature and current data of the selected monomer;
step three: constructing a neural network model, designing a training algorithm, training and optimizing the neural network model to obtain an estimation model, and inputting measured data into the temperature online estimation model for calculation;
step four: and outputting the estimated value of each monomer temperature in the module.
In the third step, the neural network model is trained by adopting an algorithm based on Levenberg-Marquardt, and the training process is as follows:
(1) initializing neural network parameters θ 0 Updating the step length lambda and the step length reduction coefficient beta according to the maximum iteration number K, wherein the current iteration number K is 0, the training error limit epsilon and the early-stop initial value P are 0, and early-stop tolerance coefficient P is obtained;
(2) selecting Sigmoid function as activation function
Figure BDA0002928538390000021
Computational neural network forward propagation
Figure BDA0002928538390000022
Wherein,
Figure BDA0002928538390000023
value, K, representing the ith neuron of the mth layer m-1 Represents the number of neurons in the m-1 layer,
Figure BDA0002928538390000024
represents the weight of the jth input neuron of the (m-1) th layer to the ith neuron of the m layer,
Figure BDA0002928538390000025
is the bias term for the ith neuron in the mth layer. Selecting a loss function of
Figure BDA0002928538390000026
Wherein,
Figure BDA0002928538390000027
is the square of the error of sample i, and l is the number of samples; calculating the loss k If loss k If ∈ or K ═ K or P ═ P, then step (5) is carried out;
(3) according to the formula theta-lambda [ J T (θ)J(θ)+μI] -1 J T (theta) e (theta) calculating network parameters, wherein mu is a constant larger than zero, I is an identity matrix, and J is a Jacobian matrix; update parameter derivation
Figure BDA0002928538390000031
And calculating the loss value loss of the model after updating the one-step parameter k+1 (ii) a If loss k+1 <loss k If yes, the parameters are updated, and
Figure BDA0002928538390000032
and p is 0, k is k +1, and λ is λ/β, and the process goes to step (2). Otherwise, turning to the step (4);
(4) making λ ═ λ β, p ═ p +1, and going to step (2);
(5) and outputting the temperature values estimated by the T1, T2, T4, T5, T7, T8 and T9.
Compared with the prior art, the invention has the following beneficial effects:
the temperature of all the monomers in the super capacitor module can be effectively estimated through a small number of temperature measurement values, and the number and the installation cost of the temperature sensors are greatly reduced. The method has simple process, is easy to realize on line and has good engineering application prospect.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic layout diagram of a super capacitor module according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a neural network model according to an embodiment of the present invention.
Detailed Description
The method comprises the following specific implementation processes:
the module partitioning method in the first step is described by taking the super capacitor module shown in fig. 2 as an example. Under the condition of no forced heat dissipation, in consideration of the symmetry of the module, 9 sections of monomers of half of the module are selected as research objects and respectively marked as T1-T9. By performing 30A, 50A, 70A and 90A constant current experiments on the super capacitor module, the temperature distribution of the tested monomer can be mainly divided into three areas. The temperature distribution result of each monomer and the arrangement mode of the supercapacitor module can be analyzed, and the monomer temperature rising rate is higher when the monomer is closer to the inner arrangement, and the monomer temperature rising rate is higher when the monomer is more in contact with the rest of the monomers.
In an example, as shown:
the temperatures of T1 and T3 were the lowest, and the number of contact surfaces with the remaining monomers was 2, so the first region was defined.
The temperature of the T2, T6, T7, T9 and T4 monomers is similar, the number of contact surfaces with the rest monomers is 3, and the monomers can be divided into a second area.
The temperatures of T5 and T8 were the highest, and the number of contact surfaces with the remaining monomers was 4, which was defined as the third region.
The selection method of the test monomer in the step two comprises the following steps: and selecting a reasonable number of monomers for testing according to the requirements on algorithm precision and occupied computing resources, and recommending one monomer to be selected in each temperature partition. In case of limited computational resources, at least one monomer may be selected. In this example, T3 and T6 were selected as input monomers.
And the module temperature online estimation model in the third step is realized based on a neural network. And taking the temperature and the current I of the input monomers T3 and T6 selected in the step two as the input of the module temperature online estimation model, wherein the output result is the temperature estimation values of the rest monomers T1, T2, T4, T5, T7, T8 and T9 in the module, and the precision of the temperature estimation values is related to the structure of the neural network model and the number of training samples. And determining the structure of the neural network by off-line training of the neural network model to obtain the temperature on-line estimation model. The neural network model built by the example is shown in fig. 3, and comprises 3 inputs, 7 outputs and 2 successive hidden layers in the middle, and the number of neurons in the hidden layers is 28 and 14 in sequence.
The neural network model is trained by adopting an algorithm based on Levenberg-Marquardt, and the Levenberg-Marquardt algorithm utilizes approximate second-order gradient information of errors, so that the method is more accurate and efficient in theory compared with a traditional gradient descent method. The training process is as follows:
(1) initializing neural network parameters θ 0 The maximum iteration number K, the updating step length lambda, the step length reduction coefficient beta, the current iteration number K is 0, the training error limit epsilon, the early-stop initial value P is 0, and the early-stop endurance coefficient P is stopped.
(2) Selecting a Sigmoid function as the activation function
Figure BDA0002928538390000041
Computational neural network forward propagation
Figure BDA0002928538390000051
Wherein,
Figure BDA0002928538390000052
represents the m-th layerValue of the ith neuron, K m-1 Represents the number of neurons in the m-1 layer,
Figure BDA0002928538390000053
represents the weight of the jth input neuron of the (m-1) th layer to the ith neuron of the m-th layer,
Figure BDA0002928538390000054
is the bias term for the ith neuron in the mth layer. Selecting a loss function of
Figure BDA0002928538390000055
Wherein,
Figure BDA0002928538390000056
is the square of the error for sample i and l is the number of samples. Calculating the loss k If loss k If ∈ or K ═ K or P ═ P, then step (5) is executed.
(3) According to the formula theta-lambda [ J T (θ)J(θ)+μI] -1 J T (θ) e (θ) calculating network parameters, where μ is a constant greater than zero, I is an identity matrix, and J is a Jacobian matrix. Update parameter derivation
Figure BDA0002928538390000057
And calculating the loss value loss of the model after updating the one-step parameter k+1 . If loss k+1 <loss k If yes, the parameters are updated, and
Figure BDA0002928538390000058
and p is 0, k is k +1, and λ is λ/β, and the process goes to step (2). Otherwise, the step (4) is carried out.
(4) And (3) enabling the step (2) to be shifted to the step (2), wherein the step (lambda) is equal to the step (beta) and the step (p) is equal to the step (p + 1).
(5) And outputting the temperature values estimated by the T1, the T2, the T4, the T5, the T7, the T8 and the T9.
The present invention is only a preferred embodiment of the present invention, and it should be noted that those skilled in the art can make various corresponding changes and modifications according to the present invention without departing from the spirit and the essence of the present invention, but those changes and modifications should fall within the protection scope of the appended claims.

Claims (2)

1. An online estimation method for temperature distribution of a super capacitor module is characterized by comprising the following steps:
the method comprises the following steps: the method comprises the following steps of (1) taking the temperature correlation of monomers in a super capacitor module into consideration, and partitioning the module by combining the arrangement characteristics of the monomers;
under the condition of no forced heat dissipation, considering the symmetry of the module, selecting 9 sections of monomers of a half of the module as research objects which are respectively marked as T1-T9;
the temperature of T1 and T3 is the lowest, and the number of contact surfaces with the rest monomers is 2, so that the first region can be divided;
the temperature of the T2, T6, T7, T9 and T4 monomers is similar, the number of contact surfaces with the rest monomers is 3, and the monomers can be divided into a second area;
the temperature of T5 and T8 is the highest, the number of contact surfaces with the rest monomers is 4, and the three-dimensional space can be divided into a third area;
step two: selecting a tested monomer according to the subareas for testing to obtain temperature and current data of the selected monomer;
the selection method of the test monomer in the step two comprises the following steps: selecting a reasonable number of monomers for testing according to the requirements on algorithm precision and occupied computing resources, and recommending one monomer to be selected in each temperature partition; in case of limited computational resources, at least one monomer may be selected; selecting T3 and T6 as input monomers;
step three: constructing a neural network model, designing a training algorithm, training, optimizing and establishing the neural network model to obtain a module temperature online estimation model, inputting measured data, and calculating in the module temperature online estimation model;
the module temperature online estimation model in the third step is realized based on a neural network; taking the temperature and current I of the input monomers T3 and T6 selected in the step two as the input of the module temperature online estimation model, wherein the output result is the temperature estimation values of the rest monomers T1, T2, T4, T5, T7, T8 and T9 in the module, and the precision of the temperature estimation values is related to the structure of the neural network model and the number of training samples; determining the structure of a neural network by off-line training of a neural network model to obtain an on-line temperature estimation model; the built neural network model comprises 3 inputs, 7 outputs and a middle 2-layer hidden layer, wherein the number of neurons in the hidden layer is 28 and 14 in sequence;
step four: and outputting the estimated value of each monomer temperature in the module.
2. The online estimation method for the temperature distribution of the supercapacitor module group according to claim 1, characterized in that in the third step, the neural network model is trained by using an algorithm based on Levenberg-Marquardt, and the training process is as follows:
(1) initializing neural network parameters θ 0 Updating the step length lambda and the step length reduction coefficient beta according to the maximum iteration number K, wherein the current iteration number K is 0, the training error limit epsilon and the early-stop initial value P are 0, and early-stop tolerance coefficient P is obtained;
(2) selecting Sigmoid function as activation function
Figure FDA0003677804410000021
Computational neural network forward propagation
Figure FDA0003677804410000022
Wherein,
Figure FDA0003677804410000023
value, K, representing the ith neuron of the mth layer m-1 Represents the number of neurons in the m-1 layer,
Figure FDA0003677804410000024
represents the weight of the jth input neuron of the (m-1) th layer to the ith neuron of the m-th layer,
Figure FDA0003677804410000025
a bias term for the ith neuron of the mth layer; selecting a loss function of
Figure FDA0003677804410000026
Wherein,
Figure FDA0003677804410000027
is the square of the error of sample i, and l is the number of samples; calculating the loss k If loss k If ∈ or K ═ K or P ═ P, then step (5) is carried out;
(3) according to the formula theta-lambda [ J T (θ)J(θ)+μI] -1 J T (theta) e (theta) calculating network parameters, wherein mu is a constant larger than zero, I is an identity matrix, and J is a Jacobian matrix; update parameter derivation
Figure FDA0003677804410000031
And calculating the loss value loss of the model after updating the one-step parameter k+1 (ii) a If loss k+1 <loss k If yes, the parameters are updated, and
Figure FDA0003677804410000032
p is 0, k is k +1, λ is λ/β and go to step (2); otherwise, turning to the step (4);
(4) making λ ═ λ β and p ═ p +1, and going to step (2);
(5) and outputting the temperature values estimated by the T1, the T2, the T4, the T5, the T7, the T8 and the T9.
CN202110140366.9A 2021-02-02 2021-02-02 Online estimation method for temperature distribution of super capacitor module Active CN112883638B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110140366.9A CN112883638B (en) 2021-02-02 2021-02-02 Online estimation method for temperature distribution of super capacitor module

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110140366.9A CN112883638B (en) 2021-02-02 2021-02-02 Online estimation method for temperature distribution of super capacitor module

Publications (2)

Publication Number Publication Date
CN112883638A CN112883638A (en) 2021-06-01
CN112883638B true CN112883638B (en) 2022-08-09

Family

ID=76052500

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110140366.9A Active CN112883638B (en) 2021-02-02 2021-02-02 Online estimation method for temperature distribution of super capacitor module

Country Status (1)

Country Link
CN (1) CN112883638B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3096389A1 (en) * 2015-05-18 2016-11-23 Nokia Technologies Oy An apparatus and associated methods for electrical storage
CN111474487A (en) * 2020-04-13 2020-07-31 重庆大学 Battery state of charge-internal temperature joint online estimation method
WO2020178616A1 (en) * 2019-03-02 2020-09-10 Qatar University Carbon nanotubes decorated with carbon nanospheres

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103675694B (en) * 2012-09-26 2016-11-02 同济大学 The online Subarea detecting system of a kind of fuel cell pile and detection method
CN109088129B (en) * 2018-08-28 2024-02-13 华霆(合肥)动力技术有限公司 Heat abstractor and battery module
CN109239489A (en) * 2018-08-29 2019-01-18 中车株洲电力机车有限公司 A kind of health status on-line analysis, device and medium for energy-storage system
CN110008557A (en) * 2019-03-28 2019-07-12 江苏理工学院 Stable and convergent supercapacitor service life On-line Estimation method
CN110474119B (en) * 2019-08-13 2024-06-11 许昌学院 Lithium ion battery pack thermal management system and thermal management method thereof
CN110687452B (en) * 2019-09-05 2022-05-20 南京理工大学 Lithium battery capacity online prediction method based on K-means clustering and Elman neural network

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3096389A1 (en) * 2015-05-18 2016-11-23 Nokia Technologies Oy An apparatus and associated methods for electrical storage
WO2020178616A1 (en) * 2019-03-02 2020-09-10 Qatar University Carbon nanotubes decorated with carbon nanospheres
CN111474487A (en) * 2020-04-13 2020-07-31 重庆大学 Battery state of charge-internal temperature joint online estimation method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
不同空间结构的有轨电车超级电容热行为;戴朝华等;《西南交通大学学报》;20201031;第55卷(第05期);正文第1-4节,第920-926页 *
电动汽车锂离子动力电池健康状态估计方法研究;黄业伟;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20150615(第06期);正文第4.1-4.8节,第31-41页 *
锂离子超级电容组高可靠性数据采集系统设计与研究;肖璆;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20130715(第07期);正文第2.1、4.2.2节,第10-11、50-56页 *

Also Published As

Publication number Publication date
CN112883638A (en) 2021-06-01

Similar Documents

Publication Publication Date Title
Lipu et al. Intelligent algorithms and control strategies for battery management system in electric vehicles: Progress, challenges and future outlook
Ren et al. A data-driven auto-CNN-LSTM prediction model for lithium-ion battery remaining useful life
Ge et al. A review on state of health estimations and remaining useful life prognostics of lithium-ion batteries
Haris et al. Early and robust remaining useful life prediction of supercapacitors using BOHB optimized Deep Belief Network
Tian et al. Feature fusion-based inconsistency evaluation for battery pack: Improved Gaussian mixture model
CN106295003B (en) Lithium battery life prediction method based on degradation track coordinate reconstruction and multiple linear regression
CN115618732B (en) Nuclear reactor digital twin key parameter autonomous optimization data inversion method
CN114548368A (en) Modeling method, prediction model and prediction method of lithium battery temperature field prediction model based on multilayer nuclear overrun learning machine
Wang et al. Health diagnosis for lithium-ion battery by combining partial incremental capacity and deep belief network during insufficient discharge profile
CN111191824B (en) Power battery capacity attenuation prediction method and system
CN115407207A (en) Online prediction method and system for lithium ion battery charge state
Li et al. Remaining useful life prediction of lithium battery using convolutional neural network with optimized parameters
CN115455842A (en) Method for estimating state of charge of super capacitor based on variable-temperature model weighted fusion
Cho et al. An LSTM-PINN hybrid method to estimate lithium-ion battery pack temperature
CN114186437A (en) Multi-physical-field coupling degradation model order reduction method for power system reliability simulation analysis
CN112731183A (en) Lithium ion battery life prediction method based on improved ELM
CN113688465B (en) Aircraft structural strength digital twin method based on combination of load and state
CN115994486A (en) Residual life prediction method of IGBT power module based on depth time sequence migration characteristic
Lin et al. Algorithm of BPNN‐UKF based on a fusion model for SOC estimation in lithium‐ion batteries
CN112883638B (en) Online estimation method for temperature distribution of super capacitor module
Ding et al. Jaya-based long short-term memory neural network for structural damage identification with consideration of measurement uncertainties
CN107301266B (en) LOC estimation method and system for lithium iron phosphate battery
CN112016241B (en) Soft measurement modeling method based on mLASSO-MLP model
CN117725859A (en) Fatigue crack growth prediction method and system for aircraft fatigue damage key part
Li et al. Machine learning algorithm based battery modeling and management method: A Cyber-Physical System perspective

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
GR01 Patent grant
GR01 Patent grant