CN116205053B - Method and device for solving simulation output working condition sequence - Google Patents
Method and device for solving simulation output working condition sequence Download PDFInfo
- Publication number
- CN116205053B CN116205053B CN202310051806.2A CN202310051806A CN116205053B CN 116205053 B CN116205053 B CN 116205053B CN 202310051806 A CN202310051806 A CN 202310051806A CN 116205053 B CN116205053 B CN 116205053B
- Authority
- CN
- China
- Prior art keywords
- working condition
- sequence
- coupling
- full
- condition sequence
- 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
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 202
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000010168 coupling process Methods 0.000 claims abstract description 171
- 238000005859 coupling reaction Methods 0.000 claims abstract description 171
- 238000004364 calculation method Methods 0.000 claims description 47
- 230000008859 change Effects 0.000 abstract description 8
- 238000012546 transfer Methods 0.000 description 24
- 239000007790 solid phase Substances 0.000 description 23
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 22
- 229910001416 lithium ion Inorganic materials 0.000 description 22
- 239000007791 liquid phase Substances 0.000 description 17
- 230000008878 coupling Effects 0.000 description 10
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 5
- 229910052744 lithium Inorganic materials 0.000 description 5
- 238000009792 diffusion process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003487 electrochemical reaction Methods 0.000 description 2
- 239000003792 electrolyte Substances 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000012798 spherical particle Substances 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 238000006479 redox reaction Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Secondary Cells (AREA)
Abstract
The application discloses a method and a device for solving an analog output working condition sequence, wherein the method comprises the steps of collecting a control input working condition sequence of a battery, receiving an input electrochemical parameter of the battery, and inputting the electrochemical parameter and the control input working condition sequence to an electrochemical solver; performing full-coupling simulation on the control input working condition sequence within a first preset time by using the electrochemical solver to obtain a pre-simulation output working condition sequence; calculating the derivative of the fitted pre-simulation output working condition sequence along with time; when the derivative is larger than a preset value, the electrochemical solver performs full-coupling simulation to obtain a full-coupling simulation output working condition sequence; and when the derivative is smaller than a preset value, the electrochemical solver performs half-coupling simulation to obtain a half-coupling simulation output working condition sequence. The method carries out half-coupling solution on the field with slow time change of the control input working condition sequence, carries out full-coupling solution on the field with fast time change, reduces the calculated amount of the electrochemical model solution, and ensures the accuracy of the calculated result.
Description
Technical Field
The application relates to the field of battery parameter identification, in particular to a method for solving a simulation output working condition sequence.
Background
There are two models for state modeling of lithium batteries at present, one is to consider the battery as a connection of abstract circuit elements, and this method is called equivalent circuit model. Another approach is to describe the physicochemical process in a lithium battery using mutually coupled partial differential equations, which model is an electrochemical model. The numerical simulation of the equivalent circuit model is simple, but the accuracy is low; the electrochemical model has high precision and predictability, and can help to know the real situation inside the lithium battery, but the numerical simulation is complex, and tens of electrochemical parameters which have real physical and chemical significance and represent the state of the internal material of the lithium battery are needed.
In solving the full-order electrochemical model, there are two general different methods for obtaining a process of simulating the gradual change of the battery working condition with time: one is to perform full-coupling solution by an electrochemical solver: all fields of the involved coupling are put together for solving, namely, PDE equation sets of electrons and lithium ions are calculated simultaneously; the second is to perform half-coupling solution by an electrochemical solver: only lithium ion PDE equations were calculated, including: solid phase mass transfer, liquid phase mass transfer, solid phase potential, liquid phase potential calculation formulas. The former is stable in numerical solution relative to the latter, and has high numerical precision, but the calculation cost is high, specifically, the former is not easy to crash a program, and the numerical precision is high, but if fields with quicker changes and fields with slower changes in all fields respectively occupy half, the time required for half-coupling solution is about one tenth of that for full-coupling solution.
For most battery conditions, half-coupling can be solved, but for too extreme battery conditions, situations where full-coupling can be solved stably but half-coupling cannot be solved often occur. At this time, the calculation cost of the full-coupling solution is too great, which can be unacceptable.
Disclosure of Invention
In order to solve the problem that full-coupling simulation is large in calculated amount and semi-coupling simulation is easy to distort under the working condition of an extreme battery, the application provides a method and a device for solving a simulation output working condition sequence.
Specifically, the technical scheme of the application is as follows:
in a first aspect, the present application provides a method for solving a sequence of simulated output conditions, comprising:
collecting a control input working condition sequence of a battery, receiving an electrochemical parameter of the battery input from the outside, and inputting the electrochemical parameter and the control input working condition sequence to an electrochemical solver;
according to the electrochemical parameters and the control input working condition sequence within a first preset time, the electrochemical solver performs full-coupling simulation to obtain a pre-simulation output working condition sequence;
fitting the pre-simulation output working condition sequence to obtain a pre-simulation output working condition curve;
calculating the derivative of the pre-simulation output working condition curve along with time;
when the derivative is larger than a preset value, the electrochemical solver continues full-coupling simulation to obtain a full-coupling simulation output working condition sequence;
and when the derivative is smaller than the preset value, the electrochemical solver performs half-coupling simulation to obtain a half-coupling simulation output working condition sequence.
According to the method, the battery control input working condition sequence is classified, and different calculation modes are used for solving, so that a simulation working condition value is obtained, the calculated amount is reduced, and the accuracy of a calculation result is guaranteed.
In some embodiments of the method for solving a simulated output condition sequence, according to the control input condition sequence within a first preset time, the electrochemical solver performs full-coupling simulation to obtain a pre-simulated output condition sequence, including:
taking the first preset time as the window length of a sliding time window, and performing full-coupling simulation for a plurality of times by using the second preset time as the simulation interval time according to the control input working condition sequence in the time range of the sliding time window while sliding the sliding time window along the time increasing direction to obtain a plurality of initial simulation output working condition sequences;
and calculating the average value of the initial simulation output working condition sequences, and taking the average value as a pre-simulation output working condition sequence.
According to the embodiment, the time range of the control input working condition sequence is changed through the sliding time window, so that the value of the control input working condition sequence is more uniform.
In some embodiments of the method of solving the sequence of analog output conditions,
the control input working condition sequence comprises a voltage time sequence or a current time sequence or a power time sequence or a temperature time sequence, and the control input working condition sequence is determined by a control mode of the battery.
The present embodiment provides a definition of a sequence of control input conditions to the electrochemical solver.
In some embodiments of the method for solving the analog output working condition sequence, when the derivative is greater than a preset value, the electrochemical solver continues to perform full-coupling simulation, and after obtaining a full-coupling analog working condition value, the method further includes:
acquiring an actual working condition sequence of the battery, and calculating working condition loss values of the actual working condition sequence and the full-coupling simulation output working condition sequence through a loss function;
when the working condition loss value is larger than a preset working condition value, the full-coupling simulation output working condition sequence is distorted, and error reporting information is output;
and outputting the full-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
The embodiment provides that after full-coupling simulation is performed, the loss function is used for comparing the full-coupling simulation output working condition sequence with the actual working condition sequence difference value, and whether the result obtained by full-coupling simulation is distorted or not is analyzed.
In some embodiments for solving the analog output working condition sequence, after the electrochemical solver performs half-coupling simulation when the derivative is smaller than the preset value, the method further includes:
collecting an actual working condition value of the battery, and calculating working condition loss values of the actual working condition value and the semi-coupling simulation working condition value through a loss function;
when the working condition loss value is larger than a preset working condition value, the half-coupling simulation working condition value is distorted, and full-coupling simulation is carried out;
and outputting the semi-coupling simulation working condition value when the working condition loss value is smaller than the preset working condition value.
The embodiment provides that after half-coupling simulation is performed, the difference value of the half-coupling simulation output working condition sequence and the actual working condition sequence is compared through a loss function, and whether a result obtained by the half-coupling simulation is distorted or not is analyzed.
In a second aspect, the present application provides an apparatus for solving a sequence of simulated output conditions, comprising:
the acquisition unit is used for acquiring a control input working condition sequence of the battery, receiving an electrochemical parameter of the battery input from the outside, and inputting the electrochemical parameter and the control input working condition sequence into the calculation unit;
the calculation unit comprises a full-coupling calculation subunit, wherein the full-coupling calculation subunit is used for carrying out full-coupling simulation by using an electrochemical solver according to the electrochemical parameters and the control input working condition sequence within a first preset time to obtain a pre-simulation output working condition sequence;
the computing unit is also used for fitting the pre-simulation output working condition sequence and computing the derivative of the pre-simulation output working condition curve obtained by fitting along with time;
the full-coupling calculation subunit is further configured to, when the derivative is greater than a preset value, continue full-coupling simulation by the electrochemical solver to obtain a full-coupling simulation output working condition sequence;
the calculation unit further comprises a half-coupling calculation subunit, and the half-coupling calculation subunit is used for performing half-coupling simulation on the electrochemical solver to obtain a half-coupling simulation output working condition sequence when the derivative is smaller than the preset value.
In some embodiments of the apparatus for solving a sequence of analog output conditions,
the calculation unit is further configured to perform full-coupling simulation for several times with a second preset time as a simulation interval time according to the control input condition sequence in the time range of the sliding time window while the sliding time window slides along the time increasing direction, so as to obtain several initial simulation output condition sequences; and calculating the average value of the initial simulation output working condition sequences, and taking the average value as the pre-simulation output working condition sequence.
In some embodiments of the apparatus for solving a sequence of analog output conditions,
the control input working condition sequence comprises a voltage time sequence or a current time sequence or a power time sequence or a temperature time sequence, and the control input working condition sequence is determined by a control mode of the battery.
In some embodiments of the apparatus for solving a sequence of analog output conditions,
the acquisition unit is also used for acquiring the actual working condition sequence of the battery and sending the actual working condition sequence to the calculation unit;
the calculation unit is further used for calculating working condition loss values of the actual working condition sequence and the full-coupling simulation output working condition sequence through a loss function when the derivative is larger than the preset value; when the working condition loss value is larger than a preset working condition value, the full-coupling simulation output working condition sequence is distorted, and error reporting information is output; and outputting the full-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
In some embodiments of the apparatus for solving a sequence of analog output conditions,
the calculation unit is further used for calculating working condition loss values of the actual working condition sequence and the semi-coupling simulation output working condition sequence through a loss function when the derivative is smaller than the preset value; when the working condition loss value is larger than a preset working condition value, the half-coupling simulation output working condition sequence is distorted, and full-coupling simulation is performed; and outputting the semi-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
Compared with the prior art, the application has at least one of the following beneficial effects:
1. according to the method, the battery control input working condition sequence is classified, one of the full-coupling simulation mode and the half-coupling simulation mode is selected to solve the simulation output working condition sequence of the battery according to the classification result, when the battery control working condition changes fast, the full-coupling simulation mode is used for solving, and when the battery control working condition changes slowly, the half-coupling simulation mode is used for solving, so that the simulation output working condition sequence is obtained. The full-coupling simulation mode is used for simultaneously calculating PDE equation sets of electrons and lithium ions at all points in space, and the half-coupling simulation mode is used for only calculating PDE equation sets of lithium ions. The method reduces the calculated amount and ensures the accuracy of the calculated result.
2. The application uses the sliding time window to select sequences with different time periods for the control input working condition sequence as the input of the electrochemical solver, so that the input control input working condition sequence can more represent the current control working condition, and further a more accurate pre-simulation output working condition sequence is obtained.
3. According to the application, the difference value between the actual working condition sequence and the obtained full-coupling analog output working condition sequence is calculated through the loss function, if the difference value accords with the expectation, the full-coupling analog output working condition sequence is normal, the full-coupling analog output working condition sequence is output, and if the difference value does not accord with the expectation, the error reporting information is output, so that the accuracy of the obtained full-coupling analog output working condition sequence is further ensured.
4. According to the application, the difference value between the actual working condition sequence and the obtained half-coupling simulation output working condition sequence is calculated through the loss function, if the difference value accords with the expectations, the half-coupling simulation output working condition sequence is normal, the half-coupling simulation output working condition sequence is output, and if the difference value does not accord with the expectations, the calculation is performed by switching to full-coupling simulation, so that the accuracy of the obtained half-coupling simulation output working condition sequence is further ensured.
Drawings
The above features, technical features, advantages and implementation of the present application will be further described in the following description of preferred embodiments with reference to the accompanying drawings in a clear and easily understood manner.
FIG. 1 is a flow chart of one embodiment of a method of solving a sequence of simulated output conditions of the present application;
FIG. 2 is a flow chart of one embodiment of a method of solving a sequence of simulated output conditions of the present application;
FIG. 3 is a flow chart of one embodiment of a method of solving a sequence of simulated output conditions of the present application;
FIG. 4 is a system block diagram of one embodiment of an apparatus for solving a sequence of simulated output conditions of the present application.
Reference numerals illustrate: 10- -an acquisition unit; 20- -a calculation unit; 21- -a fully coupled computation subunit; 22-semi-coupled computation subunit.
Detailed Description
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the following description will explain the specific embodiments of the present application with reference to the accompanying drawings. It is evident that the drawings in the following description are only examples of the application, from which other drawings and other embodiments can be obtained by a person skilled in the art without inventive effort.
For simplicity of the drawing, only the parts relevant to the application are schematically shown in each drawing, and they do not represent the actual structure thereof as a product. Additionally, in order to simplify the drawing for ease of understanding, components having the same structure or function in some of the drawings are shown schematically with only one of them, or only one of them is labeled. Herein, "a" means not only "only this one" but also "more than one" case.
It should be further understood that the term "and/or" as used in the present specification and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes such combinations.
In this context, it should be noted that the terms "mounted," "connected," and "connected" are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected, unless explicitly stated or limited otherwise; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present application will be understood in specific cases by those of ordinary skill in the art.
In addition, in the description of the present application, the terms "first," "second," and the like are used merely to distinguish between descriptions and are not to be construed as indicating or implying relative importance.
In one embodiment, referring to FIG. 1 of the drawings, the present application provides a method for solving a sequence of simulated output conditions.
S110, collecting a control input working condition sequence of the battery, receiving an externally input electrochemical parameter of the battery, and inputting the electrochemical parameter and the control input working condition sequence to an electrochemical solver.
Step S110, the battery control working conditions are collected, and the battery control working conditions and electrochemical parameters which need to be judged whether to be accurate or not are input into an electrochemical solver. The control input condition sequence includes a voltage time sequence or a current time sequence or a power time sequence or a temperature time sequence, which is determined by a control mode of the battery, for example: the battery is controlled by voltage, and the control input working condition sequence is a voltage time sequence.
And S120, performing full-coupling simulation by the electrochemical solver according to the control input working condition sequence within the first preset time to obtain a pre-simulation output working condition sequence.
Step S120, performing full-coupling simulation for a certain period of time, and simultaneously calculating PDE equation sets of electrons and lithium ions, wherein the PDE equation sets comprise solid-phase mass transfer, liquid-phase mass transfer, solid-phase potential and liquid-phase potential, the solid-phase mass transfer and the liquid-phase mass transfer represent the movement of the lithium ions, the solid-phase potential and the liquid-phase potential represent the movement of electrons, and the movement of the lithium ions and the movement of the electrons can be seen from the following formula to be coupled together through jn;
wherein η± (x, t) is the battery overpotential, ++>Is positive and negative solid phase potential->Is the liquid phase potential, U ± Is an open-circuit voltage of the positive electrode and the negative electrode>Is electrochemical lithium ion flux,/->Is the solid phase surface lithium ion concentration, F is Faraday constant,>is SEI film equivalent resistance;
and the electrochemical reaction kinetics process between the solid phase and the liquid phase meets the Butler-Volmer equation:
wherein R is a gas constant, T is a thermodynamic temperature, a + Charge transfer coefficient in negative (anode) direction, a - Is the charge transfer coefficient in the positive (cathode) direction;
where ao is an electrochemical reaction factor, ac is a redox reaction factor, c e (x, t) is a function of the liquid phase lithium ion concentration as a function of space time coordinates.
The full coupling simulation is a PDE equation set for simultaneously calculating electrons and lithium ions at all points in space, wherein the PDE equation set comprises a solid phase mass transfer calculation formula, a liquid phase mass transfer calculation formula, a solid phase potential calculation formula and a liquid phase potential calculation formula. The electrochemical solver obtains a pre-simulation output working condition sequence by utilizing solid phase mass transfer, liquid phase mass transfer, solid phase potential and liquid phase mass transfer which are obtained through full-coupling simulation.
The solid mass transfer calculation formula in the PDE equation set is:
the boundary conditions of solid phase mass transfer are:
the concentration difference of the spherical centers of the active particles is 0;
the solid phase diffusion amount of the surface of the active particles is Li + Entering electrolyte flux;
li in electrolyte at time 0 + The concentration is uniform.
Boundary conditions of solid phase mass transfer and Fick's second law can be derived simultaneously
In calculating solid phase mass transfer, radial dimension r of lithium ion/electron in spherical particle spherical coordinate system and ion concentration of lithium ion/electron at (x, r, t)Diffusion coefficient of lithium ion/electron->Substituting a solid phase mass transfer calculation formula, and substituting related parameters of lithium ions and electrons into a liquid phase mass transfer and a solid phase mass transferAnd finally obtaining solid phase mass transfer, liquid phase mass transfer, solid phase potential and liquid phase potential according to the calculation formulas of the phase potential and the liquid phase mass transfer.
S130, fitting the pre-simulation output working condition sequence to obtain a pre-simulation output working condition curve.
And S130, fitting the discrete working condition sequences into a continuous working condition curve, and conveniently deriving the following step S140 by using the working condition curve.
S140, calculating the derivative of the pre-simulation output working condition curve along with time.
And S150, when the derivative is larger than a preset value, the electrochemical solver continues full-coupling simulation to obtain a full-coupling simulation output working condition sequence.
Step S150, when the battery simulation working condition changes rapidly, full-coupling simulation is used, the calculated amount of the full-coupling simulation is large, the numerical accuracy obtained by simulation is high, and the distortion of the calculation result is avoided;
and S160, when the derivative is smaller than a preset value, the electrochemical solver performs half-coupling simulation to obtain a half-coupling simulation output working condition sequence.
Step S160, when the battery simulation working condition changes rapidly, half-coupling simulation is used, the calculated amount of the half-coupling simulation is small, the default electronic movement speed can correspond to the movement of the upper lithium ion, and only relevant parameters of the lithium ion are substituted into the PDE equation set to calculate, so that relevant parameters of the electronic are not calculated.
In the embodiment, when the change of the working condition of the simulated battery is slow, half-coupling simulation is utilized: the electrochemical solver only calculates the motion of the lithium ions, ignoring the motion of the electrons, for example: in calculating the solid phase mass transfer, only the radial dimension r of lithium ions in the spherical particle coordinate system and the ion concentration of lithium ions at (x, r, t) are calculatedDiffusion coefficient of lithium ion->Substituting the solid phase mass transfer calculation formula, so that the solving time of the electrochemical solver is reduced; in simulating battery working condition changesFast, full coupling simulation is utilized: the electrochemical solver calculates the movement of lithium ions and electrons simultaneously, and a more accurate simulation output working condition sequence is obtained, so that the accuracy of the result of the electrochemical solver is ensured.
The embodiment provides a method for solving a simulated output working condition sequence based on the foregoing embodiment, where the electrochemical solver performs full-coupling simulation according to a control input working condition sequence within a first preset time to obtain a pre-simulated output working condition sequence, and the method includes:
taking the first preset time as the window length of a sliding time window, sliding the sliding time window along the time increasing direction, and simultaneously, performing full-coupling simulation for a plurality of times by using the second preset time as the simulation interval time according to the control input working condition sequence in the time range of the sliding time window by the electrochemical solver to obtain a plurality of initial simulation output working condition sequences;
and calculating the average value of a plurality of initial simulation output working condition sequences, and taking the average value as a pre-simulation output working condition sequence.
In the embodiment, the time range of the control input working condition sequence is changed through the sliding time window, so that the value of the control input working condition sequence is wider and more uniform.
Based on the foregoing embodiment, referring to fig. 2 of the specification, this embodiment provides a method for solving a sequence of analog output conditions, in the step S150: when the derivative is larger than a preset value, the electrochemical solver continues full-coupling simulation, and after obtaining a full-coupling simulation output working condition sequence, the method further comprises the following steps:
s151, acquiring an actual working condition sequence of the battery, and calculating working condition loss values of the actual working condition sequence and the full-coupling simulation output working condition sequence through a loss function;
s152, judging whether the working condition loss value is larger than a preset working condition value;
s152, when the working condition loss value is larger than a preset working condition value, fully coupling simulation output working condition sequence distortion is carried out, and error reporting information is output;
and S153, outputting a full-coupling simulation output working condition sequence when the working condition loss value is smaller than a preset working condition value.
The embodiment provides that after full-coupling simulation is performed, the difference value of the full-coupling simulation output working condition sequence and the actual working condition sequence is compared through the loss function, whether the result obtained by full-coupling simulation is distorted is analyzed, the accuracy of the obtained result is further ensured, and the distorted full-coupling simulation output working condition sequence is filtered, for example: using a mean square error loss functionCalculated, wherein MSE is a working condition loss value, n is the number of sampling time points, i is a positive integer, U 1，i Is the value of the full-coupling analog output working condition sequence of the battery at the ith time point, U 2，i Is the value of the actual battery operating mode sequence of the battery at the ith time point.
Based on the foregoing embodiment, referring to fig. 3 of the specification, the present embodiment provides a method for solving a sequence of analog output conditions, where in step S160, when the derivative is smaller than a preset value, the electrochemical solver performs half-coupling simulation to obtain a sequence of half-coupling analog output conditions, and then further includes:
s161, acquiring an actual working condition sequence of the battery, and calculating working condition loss values of the actual working condition sequence and the semi-coupling simulation output working condition sequence through a loss function;
s162, judging whether the working condition loss value is larger than a preset working condition value;
s163, when the working condition loss value is larger than a preset working condition value, the half-coupling simulation output working condition sequence is distorted, and full-coupling simulation is carried out;
and S164, outputting a semi-coupling simulation output working condition sequence when the working condition loss value is smaller than a preset working condition value.
The embodiment provides that after half coupling simulation is performed, the difference value of the half coupling simulation output working condition sequence and the actual working condition sequence is compared through the loss function, whether the result obtained by half coupling simulation is distorted is analyzed, when the result is distorted, more complex full coupling is adopted for operation, and the situation that the working condition of a battery is extremely extreme is avoided, for example: when the temperature is too low and the current is too large, the situation that full coupling can be solved stably but half coupling cannot be solved occurs.
In one embodiment, referring to fig. 4 of the specification, the apparatus for solving a sequence of analog output conditions provided by the present application includes:
the acquisition unit 10 is configured to acquire a control input condition sequence of the battery, receive an externally input electrochemical parameter of the battery, and input the electrochemical parameter and the control input condition sequence to the calculation unit 20;
the calculation unit 20 is connected with the acquisition unit 10 and comprises a full-coupling calculation subunit 21, and the full-coupling calculation subunit 21 is used for carrying out full-coupling simulation by using an electrochemical solver according to the electrochemical parameters and the control input working condition sequence within a first preset time to obtain a pre-simulation output working condition sequence;
the calculating unit 20 is further configured to fit the pre-simulation output working condition sequence, and calculate a derivative of the pre-simulation output working condition curve obtained by the fitting along with time;
the full-coupling calculation subunit 21 is further configured to, when the derivative is greater than a preset value, continue full-coupling simulation by the electrochemical solver to obtain a full-coupling simulation output working condition sequence;
the computing unit 20 further includes a half-coupling computing subunit 22, where the half-coupling computing subunit 22 is configured to perform half-coupling simulation by the electrochemical solver when the derivative is less than a preset value, so as to obtain a half-coupling simulation output working condition sequence.
In this embodiment, the acquisition unit 10 collects the battery control input working condition sequence, and the calculation unit 20 selects half-coupling simulation or full-coupling simulation according to the working condition change speed of the collected control input working condition sequence, so as to quickly solve the simulation output working condition sequence. Under the extreme condition of the guarantee working condition, the full-coupling simulation is used, the simulation output working condition sequence is accurately solved, and under the condition of small working condition change, the semi-coupling simulation is used, the movement of electrons is ignored, the solution of a related equation set of lithium ions is carried out, so that the calculated amount is saved.
The present embodiment provides a device for solving a sequence of analog output conditions based on the foregoing embodiment,
the calculating unit 20 is further configured to perform full-coupling simulation for several times with the second preset time as a simulation interval time according to the control input working condition sequence in the time range of the sliding time window while the sliding time window slides along the time increasing direction with the first preset time as the window length of the sliding time window, so as to obtain several initial simulation output working condition sequences; and calculating the average value of a plurality of initial simulation output working condition sequences, and taking the average value as a pre-simulation output working condition sequence.
In the embodiment, the time range of the control input working condition sequence is changed through the sliding time window, so that the value range of the control input working condition sequence is wider and more uniform.
The embodiment provides a device for solving the sequence of analog output working conditions based on the previous embodiment,
the acquisition unit 10 is further used for acquiring an actual working condition sequence of the battery and sending the actual working condition sequence to the calculation unit;
the calculating unit 20 is further configured to calculate, when the derivative is greater than a preset value, a working condition loss value of the actual working condition sequence and the full-coupling analog output working condition sequence through a loss function; when the working condition loss value is larger than a preset working condition value, the full-coupling simulation output working condition sequence is distorted, and error reporting information is output; and outputting the full-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
After full-coupling simulation is used for the pre-simulation output working condition sequence with rapid working condition change, working condition loss values between the simulation output working condition sequence and the actual working condition sequence are obtained through analysis, when the working condition loss values are large, data are distorted, and when the working condition loss values are small, the data are normal, and the data are output.
The embodiment provides a device for solving the sequence of analog output working conditions based on the previous embodiment,
the calculating unit 20 is further configured to calculate, when the derivative is smaller than a preset value, a working condition loss value of the actual working condition sequence and the semi-coupling analog output working condition sequence through a loss function; when the working condition loss value is larger than a preset working condition value, the half-coupling simulation output working condition sequence is distorted, and full-coupling simulation is performed; and outputting a semi-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
In the embodiment, after half-coupling simulation is used for the pre-simulation output working condition sequence with slower working condition change, working condition loss values between the obtained simulation output working condition sequence and the actual working condition sequence are analyzed, when the working condition loss values are larger, data distortion is calculated by using full-coupling simulation, and when the working condition loss values are smaller, data are normal and data are output.
It should be noted that the above embodiments can be freely combined as needed. The foregoing is merely a preferred embodiment of the present application and it should be noted that modifications and adaptations to those skilled in the art may be made without departing from the principles of the present application, which are intended to be comprehended within the scope of the present application.
Claims (8)
1. A method for solving a sequence of analog output conditions, comprising:
collecting a control input working condition sequence of a battery, receiving an electrochemical parameter of the battery input from the outside, and inputting the electrochemical parameter and the control input working condition sequence to an electrochemical solver;
according to the electrochemical parameters and the control input working condition sequence within a first preset time, the electrochemical solver performs full-coupling simulation to obtain a pre-simulation output working condition sequence;
fitting the pre-simulation output working condition sequence to obtain a pre-simulation output working condition curve;
calculating the derivative of the pre-simulation output working condition curve along with time;
when the derivative is larger than a preset value, the electrochemical solver continues full-coupling simulation to obtain a full-coupling simulation output working condition sequence;
when the derivative is smaller than the preset value, the electrochemical solver performs half-coupling simulation to obtain a half-coupling simulation output working condition sequence;
the electrochemical solver performs full-coupling simulation according to the control input working condition sequence within a first preset time to obtain a pre-simulation output working condition sequence, and the method comprises the following steps:
taking the first preset time as the window length of a sliding time window, and performing full-coupling simulation for a plurality of times by using the second preset time as the simulation interval time according to the control input working condition sequence in the time range of the sliding time window while sliding the sliding time window along the time increasing direction to obtain a plurality of initial simulation output working condition sequences;
and calculating the average value of the initial simulation output working condition sequences, and taking the average value as a pre-simulation output working condition sequence.
2. A method of solving a sequence of analog output conditions according to claim 1, wherein:
the control input working condition sequence comprises a voltage time sequence or a current time sequence or a power time sequence or a temperature time sequence, and the control input working condition sequence is determined by a control mode of the battery.
3. The method of any one of claims 1-2, wherein after said electrochemical solver continues full-coupling simulation when said derivative is greater than a predetermined value to obtain a full-coupling simulated output condition sequence, further comprising:
acquiring an actual working condition sequence of the battery, and calculating working condition loss values of the actual working condition sequence and the full-coupling simulation output working condition sequence through a loss function;
when the working condition loss value is larger than a preset working condition value, the full-coupling simulation output working condition sequence is distorted, and error reporting information is output;
and outputting the full-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
4. The method of any one of claims 1-2, wherein after performing a semi-coupled simulation by the electrochemical solver when the derivative is less than the predetermined value to obtain a semi-coupled simulated output condition sequence, further comprising:
acquiring an actual working condition sequence of the battery, and calculating working condition loss values of the actual working condition sequence and the semi-coupling simulation output working condition sequence through a loss function;
when the working condition loss value is larger than a preset working condition value, the half-coupling simulation output working condition sequence is distorted, and full-coupling simulation is performed;
and outputting the semi-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
5. An apparatus for solving a sequence of analog output conditions, comprising:
the acquisition unit is used for acquiring a control input working condition sequence of the battery, receiving an electrochemical parameter of the battery input from the outside, and inputting the electrochemical parameter and the control input working condition sequence into the calculation unit;
the calculation unit is connected with the acquisition unit and comprises a full-coupling calculation subunit, and the full-coupling calculation subunit is used for carrying out full-coupling simulation by utilizing an electrochemical solver according to the electrochemical parameters and the control input working condition sequence within a first preset time to obtain a pre-simulation output working condition sequence;
the computing unit is also used for fitting the pre-simulation output working condition sequence and computing the derivative of the pre-simulation output working condition curve obtained by fitting along with time;
the full-coupling calculation subunit is further configured to, when the derivative is greater than a preset value, continue full-coupling simulation by the electrochemical solver to obtain a full-coupling simulation output working condition sequence;
the calculation unit further comprises a half-coupling calculation subunit, wherein the half-coupling calculation subunit is used for performing half-coupling simulation on the electrochemical solver to obtain a half-coupling simulation output working condition sequence when the derivative is smaller than the preset value;
the calculation unit is further configured to perform full-coupling simulation for several times with a second preset time as a simulation interval time according to the control input condition sequence in the time range of the sliding time window while the sliding time window slides along the time increasing direction, so as to obtain several initial simulation output condition sequences; and calculating the average value of the initial simulation output working condition sequences, and taking the average value as the pre-simulation output working condition sequence.
6. The apparatus for solving a sequence of analog output conditions of claim 5, wherein:
the control input working condition sequence comprises a voltage time sequence or a current time sequence or a power time sequence or a temperature time sequence, and the control input working condition sequence is determined by a control mode of the battery.
7. An apparatus for solving a sequence of analog output conditions according to any one of claims 5-6,
the acquisition unit is also used for acquiring the actual working condition sequence of the battery and sending the actual working condition sequence to the calculation unit;
the calculation unit is further used for calculating working condition loss values of the actual working condition sequence and the full-coupling simulation output working condition sequence through a loss function when the derivative is larger than the preset value; when the working condition loss value is larger than a preset working condition value, the full-coupling simulation output working condition sequence is distorted, and error reporting information is output; and outputting the full-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
8. The apparatus for solving a series of analog output conditions according to claim 7,
the calculation unit is further used for calculating working condition loss values of the actual working condition sequence and the semi-coupling simulation output working condition sequence through a loss function when the derivative is smaller than the preset value; when the working condition loss value is larger than a preset working condition value, the half-coupling simulation output working condition sequence is distorted, and full-coupling simulation is performed; and outputting the semi-coupling simulation output working condition sequence when the working condition loss value is smaller than the preset working condition value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310051806.2A CN116205053B (en) | 2023-02-02 | 2023-02-02 | Method and device for solving simulation output working condition sequence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310051806.2A CN116205053B (en) | 2023-02-02 | 2023-02-02 | Method and device for solving simulation output working condition sequence |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116205053A CN116205053A (en) | 2023-06-02 |
| CN116205053B true CN116205053B (en) | 2023-12-01 |
Family
ID=86518398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310051806.2A Active CN116205053B (en) | 2023-02-02 | 2023-02-02 | Method and device for solving simulation output working condition sequence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116205053B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110741267A (en) * | 2018-02-01 | 2020-01-31 | 株式会社Lg化学 | Method for estimating parameters of equivalent circuit model for battery and battery management system |
| CN113281655A (en) * | 2021-05-20 | 2021-08-20 | 中南大学 | Predictive control method and device for internal heating of power battery in low-temperature environment |
| CN115062474A (en) * | 2022-06-17 | 2022-09-16 | 上海玫克生储能科技有限公司 | Working condition simulation method and system of lithium battery, terminal device and storage medium |
| CN115116554A (en) * | 2022-06-29 | 2022-09-27 | 上海玫克生储能科技有限公司 | Electric field decoupling method and device of electrochemical model |
| CN115512777A (en) * | 2022-09-23 | 2022-12-23 | 上海玫克生储能科技有限公司 | Electrochemical model parameter identification method and system based on capacity change rate |
-
2023
- 2023-02-02 CN CN202310051806.2A patent/CN116205053B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110741267A (en) * | 2018-02-01 | 2020-01-31 | 株式会社Lg化学 | Method for estimating parameters of equivalent circuit model for battery and battery management system |
| CN113281655A (en) * | 2021-05-20 | 2021-08-20 | 中南大学 | Predictive control method and device for internal heating of power battery in low-temperature environment |
| CN115062474A (en) * | 2022-06-17 | 2022-09-16 | 上海玫克生储能科技有限公司 | Working condition simulation method and system of lithium battery, terminal device and storage medium |
| CN115116554A (en) * | 2022-06-29 | 2022-09-27 | 上海玫克生储能科技有限公司 | Electric field decoupling method and device of electrochemical model |
| CN115512777A (en) * | 2022-09-23 | 2022-12-23 | 上海玫克生储能科技有限公司 | Electrochemical model parameter identification method and system based on capacity change rate |
Non-Patent Citations (2)
| Title |
|---|
| Enhanced Heat Transfer in Porous Media for Reducing Degradation of the of Lithium-Ion Batteries;Masoud Jabbari 等;《The Electrochemical Society》;摘要 * |
| 磷酸铁锂电池放电过程的电化学—热全耦合模型数值模拟研究;廖连生;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》(第12期);C042-1399 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116205053A (en) | 2023-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109143083B (en) | data-driven lithium ion battery lithium analysis diagnosis method for electric vehicle | |
| CN113702845B (en) | Retired lithium battery core parameter evaluation method and equipment | |
| CN113176503B (en) | Full SOC range lithium ion battery equivalent model based on electrochemical process | |
| CN109991548A (en) | A kind of OCV-SOC method of calibration experiment, battery equivalent model parameter identification method and SOC estimation method | |
| CN112182890A (en) | Lithium ion battery electrochemical model for low-temperature application | |
| CN110442901B (en) | Electrochemical simplified model of lithium ion battery and parameter acquisition method thereof | |
| CN110888070A (en) | Battery temperature estimation method, device, equipment and medium | |
| Mu et al. | Modeling, evaluation, and state estimation for batteries | |
| CN114755582A (en) | Lithium ion battery health state estimation method when environmental temperatures are different | |
| Li et al. | Linear parameter-varying modeling and identification of lithium-ion battery for control-oriented applications | |
| Jiang et al. | State-of-charge estimation of lithium-ion batteries from a gas-liquid dynamics model including the direct temperature input | |
| CN115116555B (en) | Electric field decoupling method and device of electrochemical model based on targeting method | |
| Wang et al. | State of health estimation of lithium-ion battery in wide temperature range via temperature-aging coupling mechanism analysis | |
| CN116205053B (en) | Method and device for solving simulation output working condition sequence | |
| CN114912327B (en) | Method and device for calculating current density in battery electrochemical model and storage medium | |
| Xin et al. | A novel state of charge estimation method for ternary lithium batteries based on system function and extended kalman filter | |
| Guo et al. | Nonlinear state-variable method for solving physics-based li-ion cell model with high-frequency inputs | |
| CN114865117A (en) | Lithium ion battery electrode lithium embedding amount detection method and device and battery management system | |
| Gao et al. | Determination of half-cell open-circuit potential curve of silicon-graphite in a physics-based model for lithium-ion batteries | |
| Ma et al. | State of charge and state of health estimation based on dual nonlinear adaptive observer and hysteresis model of lithium-ion battery | |
| CN113642153B (en) | Reference electrode measurement error calibration method and computer equipment | |
| CN116165547A (en) | Method, system, device and storage medium for calculating and correcting voltage of lithium battery | |
| CN116125314A (en) | Lithium battery state of charge estimation method based on electrochemical model and particle filtering | |
| Kerdja et al. | Microscopy imaging based numerical model of Li-ion batteries electrode: a parametric study through a design of experiment approach | |
| Chen et al. | Research on state-of-charge estimation of lithium-ion batteries based on an improved gas-liquid dynamics model |
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 | ||
| CP03 | Change of name, title or address |
Address after: Room 1101, No. 2, Lane 288, Qianfan Road, Xinqiao Town, Songjiang District, Shanghai 201612 Patentee after: Shanghai Meikesheng Energy Technology Co.,Ltd. Country or region after: China Address before: Room 1101, No. 2, Lane 288, Qianfan Road, Xinqiao Town, Songjiang District, Shanghai Patentee before: Shanghai Meike Energy Storage Technology Co.,Ltd. Country or region before: China |