Next Article in Journal
Enhancement of Heavy-Duty Engines Performance and Reliability Using Cylinder Pressure Information
Previous Article in Journal
Potential of CBM as an Energy Vector in Active Mines and Abandoned Mines in Russia and Europe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Digital Twin of Li-Ion Battery Based on Generic MATLAB Model

Department of Electrical Engineering and Mechatronics, Faculty of Electrical Engineering and Informatics, Technical University of Kosice, Letna 9, 042 00 Kosice, Slovakia
*
Author to whom correspondence should be addressed.
Energies 2023, 16(3), 1194; https://doi.org/10.3390/en16031194
Submission received: 8 December 2022 / Revised: 9 January 2023 / Accepted: 17 January 2023 / Published: 21 January 2023
(This article belongs to the Section D2: Electrochem: Batteries, Fuel Cells, Capacitors)

Abstract

:
The paper describes the digital twin of a Li-ion battery cell based on the MATLAB/Simulink generic model. The digital twin is based on measured data for constant current/constant voltage charging and discharging cycles with State of Health (SoH) up to 79%, also including fast charging. Mathematical equations used for the digital twin are obtained by 3D data fitting of measured SoH, battery capacity, and battery cell current. The input to the proposed digital twin is only the measured battery cell current, and its output includes State of Charge (SoC), SoH, and battery cell voltage. The designed digital twin is tested and compared with MATLAB/Simulink generic model and battery cell measurements for constant discharging current and dynamically generated discharging current profile. The results show significant improvement in the generic MATLAB/Simulink model.

1. Introduction

With rapid development in electric vehicles and consumer electronics, smart grids, etc., energy storage systems are becoming an important sector. Electrical energy generation is shifting from predictable and stable sources such as power plants to less predictable and stable sources such as renewable energy sources. Carbon neutrality, for example, net-zero greenhouse gas emissions, which the EU aims to reach by 2050 [1], will place significant emphasis on highly reliable battery energy storage systems (BESS) and advanced smart grids. As the world’s second-largest producer of motor vehicles after China, the EU aims for zero-emission vehicles from 2035 [2]. All these applications require high-energy BESS to meet the demand for reliable systems used for a short time and highly efficient energy storage. The BESS are small in size, portable, and can be placed everywhere. Batteries are electrochemical energy storage systems; the energy stored in the battery cannot be accessed directly as in chemical storage systems such as hydrogen. This implies that it is not easy to assess the battery state-of-charge (SoC) instantaneously. The SoC is one of the most important parameters of the BESS. In general, the battery SoC is defined as the current capacity ratio to the battery’s nominal capacity. Moreover, the BESS must be safe during the whole life cycle as the battery performance gradually deteriorates. Lithium-ion (Li-ion) batteries are preferred over other battery chemistry in various applications because of their long life, small self-discharge, high energy density, no requirement for complete discharge (the memory effect), and lighter weight. The battery management system (BMS) is responsible for ensuring that the BESS is operated within a safe operating range. However, many accidents are being reported as a result of faults in the BESS [3]. In this regard, the current research of BMS is focused on precise battery State-of-Health (SoH) estimation, and Remaining-Useful-Life (RUL), which need to be monitored to prevent failure [3,4]. The SoH is defined as a ratio of the maximal instantaneous releasable capacity to the new battery’s capacity [5]. The RUL is a subjective estimate of the remaining number of months or years that the BESS is estimated to be able to function in accordance with its purpose before it will need to be replaced [4].
The basic property of a battery is its SoC, which needs to be precisely monitored by BMS. The existing methods for SoC estimation can be divided into three categories: direct, indirect, and adaptive [4]. Direct methods are based on measuring the physically accessible battery properties such as the battery voltage, or battery current [6,7,8]. Indirect methods are data-based and use a database of measured battery properties within a cycle [9,10]. Finally, adaptive methods use adaptive filters formed from Fuzzy logic [11], artificial neural networks [12], and the Kalman filter [12], which are based on adaptive algorithms for adjusting parameters. The direct methods combining coulomb counting and open circuit voltage measurements are the easiest to implement and provide reliable results in periodic battery cycling [6,7,8,9,10,11,12,13].
The SoH is critical for the proper operation of the BMS as it represents the current state of the battery within its lifetime. Four approaches are used for battery degradation modelling [10]: physical-based models, equivalent circuit models, machine learning models, and empirical and semi-empirical models. Empirical models are simple and easy to implement as they contain only mathematical formulas which try to describe complex electrochemical processes. This approach is intuitive for non-chemical engineers; however, the empirical models often over-simplify the complex behaviour of Li-ion batteries. Empirical models are based on curve fitting the relationship of various factors influencing the performance of the battery, such as cut-off charge voltage, temperature, SoH, DoD etc. Large data sets are necessary to decouple their mutual influence. Data sets for various charging (usually constant current followed by constant voltage) and discharging profiles with variable temperatures are available as open data [9,14,15,16,17]. They are used for developing ageing models, for example, for LiFePO4 with fast-charging up to 10C [18] or for LiCoO2 degradation modelling [19]. However, these datasets usually do not include fast charging and consider galvanostatic discharging.
MATLAB/Simulink has been used for many applications in estimations of SoC and SoH of a battery cell, including the Extended Kalman filter, artificial neural networks, neuro-fuzzy systems and various equivalent circuit models. In [20], the battery SoC estimation equivalent circuit provided in MATLAB has been modified by adding the 3-RC pairs in series with its internal resistance. In [21], authors propose a MATLAB function with the objective of providing a public tool that estimates the battery SoC and terminal voltage at different temperatures using a second-order resistor-capacitor (2RC) ECM along with an extended Kalman filter (EKF). In [22], authors propose a unified OCV–SOC model that is intended for both SOC estimation and SoH monitoring. In [23], a universal method for tuning the cell model against any standard dataset is proposed. A cell model of 1-RC pair was developed using MATLAB/Simulink. A new methodology is adopted to better estimate SoC by tuning the cell parameters. In [24], the authors propose a method for estimating the SoH of a Li-ion battery based on the assumption that limited data is obtained through a small number of reference performance tests (RPT). Performing frequent RPTs causes unnecessary degradation of the Li-ion battery. The authors design a neural network for SoH estimation using RPT-reduced experimental data proposed for achieving economic efficiency and mitigating the dispensable degradation due to RPT operation.
The paper describes the development of an empirical battery model based on a proprietary data set obtained by battery cycling for Samsung INR18650-25R cell (Figure 1b) and also includes fast charging and its influence on battery lifetime. The empirical battery model is developed in MATLAB/Simulink environment and is based on the generic battery model from MATLAB/Simulink. SoC estimation is based on coulomb counting, as the only measured variable which enters the model is the battery current. The initial capacity of the battery cell (corresponding to the initial SoH) is measured by a single charge and discharge cycle by a precise power analyser. The SoH model is based on data obtained from cell cycling up to SoH equal to 80% of the initially measured SoH. The data for empirical model creation, thus, contains information about the whole life cycle. The created model can be easily implemented for battery management system development. Comparison between measured data and the designed empirical model shows a good correlation even for randomly generated discharging profiles. Li-ion batteries are produced in various shapes (cylindrical, prismatic, button, pouch, etc.) (Figure 1a).
In Section 2, the basic charging and discharging parameters of Li-ion batteries regarding the lifetime of the battery are presented. Section 3 describes a methodology for collecting the data set used for model creation. Section 4 describes the test bench used for data collection during the cyclic test. In Section 5, the battery tests are described and obtained data sets are presented for SoH up to 79%. Section 6 presents the empirical model and its concept. Moreover, data fitting and the final model are described in Section 6. In Section 7, the comparison between measured data and created model summarises the properties of the model for randomly generated battery discharging profile.

2. Li-Ion Batteries

Li-ion battery cells have many advantages when compared to other battery chemistry. They offer relatively high capacities (188 mAh/g for NCA cathode cell) and high energy density (260 Wh/kg for NCA cathode cell) whilst being lightweight, accept fast charging, have very low self-discharge (0.35–2.5% per month) and have a long lifespan of the battery cell [26,27,28,29]. They also have relatively high nominal voltage (2.5 V–4.2 V) and no memory effect. Li-ion batteries can be operated in a wide temperature range from −20 °C to +60 °C
One of their disadvantages is that if we want to achieve a long service life for the battery cell, it is necessary to treat them according to certain requirements:
  • Do not exceed their maximal charging and discharging current
  • Maintain Constant Current–Constant Voltage (CC-CV) charging algorithm
  • Protect the cell from overcharging/undercharging
  • Protect the cell from over-temperature by proper thermal management.
Li-ion battery cells have a lifespan of approximately 300–2000 cycles, depending on the cathode chemistry, Depth of Discharge (DoD), and temperature [30]. As can be seen in Figure 1a, Li-ion batteries are produced in various shapes (cylindrical, prismatic, button, pouch, etc.) and in most cases have a nominal voltage of 3.6 V or 3.7 V [30,31,32].

3. Testing Method of Li-Ion Batteries

Batteries can be tested from several points of view, such as the electrical, mechanical or thermal properties. The main task was to retrieve models based on measured data from cycling testing accurately. This section focuses on the proper charging and discharging process of Li-ion batteries. In addition, this section also describes the protections that must be ensured during these processes. For this purpose, it was necessary to build test equipment that performs cyclic charging and discharging of battery cells. The performed tests were focused precisely on the electrical properties of the battery cell. These tests take into account the battery’s ability to deliver a certain voltage, current and capacity. Such testing takes place by repeatedly charging and discharging battery cells [33,34].

3.1. Charging

The optimal charging process for Li-ion cells is called CC/CV (Constant Current/Constant Voltage) and consists of two phases:
  • The constant current (CC) phase is the phase during which a regulated current with a constant value flows into the battery cell. A constant current in the range of 0.5–1 C is usually chosen (if the cell has a capacity of 2500 mAh 1C = 2.5 A), depending on the exact type of the battery cell. This current must be maintained by a converter with current control. Thus, the current is constant, and the voltage on the battery cell slowly begins to increase. This first phase, i.e., the constant current phase, is used until the moment when the battery cell’s maximum allowable voltage, determined by the manufacturer, appears on the battery cell terminals. In most cases, it is a value of around 4.2 V. After reaching this voltage value, the charging process continues with the next phase [35,36].
  • The constant voltage phase (CV) is followed by the CC phase and is the phase during which the battery cell is no longer able to receive a current with the value used in the first phase without its output voltage increases above the maximum allowed voltage value set by the manufacturer. Since the battery cell is not yet fully charged at the beginning of this phase, the charger creates a constant voltage applied to the cell with the maximum value determined by the manufacturer. In most cases, this value is around 4.2 V. Thanks to this, the battery will not exceed this voltage value at its terminals. However, a current will start flowing to the battery cell and will charge it up to 100%. The cell current is exponentially decreasing during the CV phase. The current value decreases until the minimum current value when the battery is considered fully charged and the charging process is terminated [35,36].
The overall charging process described above can be seen graphically in Figure 2. The length of the individual phases is significantly influenced by the exact type of battery and especially by the maximum current in the first phase. The boost charging method uses an additional higher current at the beginning of the charging cycle to shorten the charging time without compromising the cycle life. Moreover, the charging time is influenced by the value of the changing current in the classical CC-CV charging method [37,38,39].

3.2. Discharging

From the operation point of view, the discharge process of the battery cell is significantly simpler as it consists only of one phase. The discharge profile can be set to constant current (e.g., galvanostatic) and is terminated by reaching the cut-off voltage. However, in real applications, the discharging current is variable, and it is important to monitor the current during discharge, the minimum (cut-off) voltage and the temperature of the cell. The DoD significantly influences the cell cycle life as Li-ion cells prefer shallow discharge, and lowering DoD increases the cyclability of the battery cell [41,42].
Discharging with higher than permitted currents, under-discharging below minimum voltages and discharging the battery cell outside the working temperature range have a significant impact on its lifetime. Typical discharge characteristics for various current values can be seen in Figure 3. During discharging with low current, the cell remains relatively cold and discharging lasts longer. With higher discharging current also increase cell temperature and decreases discharging time. The discharging, which includes micro-cycling, for example, a small swing in the state of charge, can have a significant impact on the cycle life of the battery cell and its influence on the cycle life needs to be studied [42].

4. Device for Cyclic Testing of Battery Cells

The battery cell data, such as its voltage and current during charging and discharging, are necessary to create the empirical model of the cell. These data can be found in open access databases such as [14], but we have decided to obtain our own data for the modelled battery cell Samsung INR18650-25R. There are many different types of battery testers with different features and prices on the market, but we have decided to use our own tester whose design was part of a PhD thesis. The designed testing device for battery cell cycling ensures the charging and discharging of battery cells throughout its lifetime and is used for testing the chosen battery cell model. The designed device can also be extended by a thermal box with very precise temperature control in the future. The requirements for the test equipment are as follows:
  • Charging and discharging of various types of Li-ion cells
  • Possibility to charge/discharge up to 6 battery cells together in one group.
  • The possibility of setting the charging and discharge profile
  • Charging current max 5 A, discharging max 30 A
  • Current accuracy with a sampling rate
  • Protections: reverse voltage, undervoltage, overvoltage, short circuit, overheating
  • Monitoring of voltage, current, temperature and battery cell DoD
Based on these requirements, the concept of the entire test equipment was designed.

4.1. Concept

A concept with one master unit and several slave units was chosen for the development of the device for cycling testing. In this concept, every battery cell was measured by a dedicated device with its own microcontroller, which is able to work independently of the remaining test devices. The advantage of this solution is the robustness of the system when the devices will be able to work even in the event of failures on the remaining devices. Figure 4 shows a block diagram of the cycling battery cell tester. The blocks are divided into three groups with different colours.
The power parts and paths where a higher current is expected are shown in red. The sensing, control and communication parts of the device are marked by blue. The yellow colour shows the controllers. In the case of digital control, these regulators were implemented in a microcontroller. According to the requirements, we want to test at least six battery cells together in one group. It means that at least six test modules will be connected to the master processor equipped with the HMI interface [44,45].

4.2. Developed Test Device

The test device, which was developed based on requirements, is shown in Figure 5. The device is based on Cortex M4 single-core ARM architecture microcontroller from the STM32 family. The microcontroller provides a sufficient number of GPIO pins, PWM pins and pins for the AD converter. The microcontroller also offers a hardware interface for CAN bus connection.
The power stage of the device is based on two buck DC/DC converters connected in parallel. The first converter ensures the charging of the battery cell. The second DC/DC converter ensures discharging battery cell through 0.05 Ω resistor with rated power 100 W. Measured data are collected to the SD card.
The current version of the developed device has these options:
  • Constant current load test. Charging current to 5 A and discharging current to 20 A with auxiliary cooling.
  • Possibility to load discharging profile in .csv format, which can be modified later.
  • Possibility to cycle power charging/discharging cells to the defined SoH value. the SoH value can be set in a graphical user interface.
  • Possibility to simulate the overcharging, undercharging and fast charging.
  • Implemented protections: overvoltage, undervoltage, overheating and short circuit

5. Cell Tests

The Li-ion cell Samsung INR18650-25R was chosen as a reference for testing. The parameters of this battery cell are shown in Figure 6. The designed testing device for cycling tests uses the Nucleo STM32F303-K8 microcontroller with 12-bit ADC. The resolution for charging and discharging current was 4.8 mA and 13.3 mA, respectively. To keep the error low, an open-loop hall-effect current sensor was used for charging and discharging current measurements. The sampling rate of the ADC was set to 10 ksps. The testing device was used only to “wear out” the battery cells, and current sensors were needed to maintain constant charging and discharging current. Data for model creation and curve fitting (discharging curves and SoH measurement) was by a programmable DC load IT8518E and precision power analyser LMG500. The discharge process was terminated by reaching the cut-off voltage set in the programmable DC load. The accuracy of the current and voltage measurement for discharge curves and SoH measurement (Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12) was 0.045% and 0.05%, respectively. The battery cell was subjected to several hundreds of charge and discharge cycles, and then the SoH was measured.
The following series of tests were performed on this type of Li-ion cell with the aim of obtaining the maximum possible amount of data to create an accurate empirical model of the Li-ion cell:
  • Determining the maximum number of cycles which the battery cell can perform before its SoH drops to 80% at nominal charging and discharging current.
  • Determining the influence of discharging current changes to the battery cell ageing with nominal charging current and discharging current 5A, 10A, 20A to 80% SoH.
  • Determining the influence of change charging current to battery cell ageing with nominal discharging current and different charging current up to 80% SoH.
  • Fast charging with maximal current and subsequent discharging with current up to 20 A till the SoH reached 80%
  • Measuring of battery capacity with different charging and discharging currents and different SoH up to 80% SoH.
  • Measuring shapes of charging and discharging curves for different charging and discharging currents and different SoH up to 80% SoH.
Measured results are shown in the figures below. In Figure 7, discharging characteristics for SoH 100% reflect the ability to supply high current while preserving high capacity—it can be achieved by a small voltage drop on the internal resistance of the battery cell.
Figure 7 also shows a ripple in discharging part of the curve for current 1 A. This ripple does not apply to higher current curves.
SoH = Q m a x Q n o m · 100 [ % ]
where:
Q m a x —maximal capacity for a given cycle in which SoH is being calculated,
Q n o m —initial capacity of a new cell measured by one charge and discharge cycle.
One cycle consisted of a CC-CV charge with an initial current value of 0.4 C followed by a predefined constant current discharge (from 0.4 C to 8 C) until the cut-off voltage was reached. Then the next cycle was initiated without any rest period.
Characteristics in Figure 8 were measured at SoH = 94%. The effect of increasing internal resistance elements between curves for currents of 5 A and 10 A is clearly visible. It is also possible to see nonlinearities at less than the nominal current for this SoH. It can be seen that the cell at SoH = 94% does not achieve such excellent results, especially at currents higher than 10 A. At SoH = 94% and a current of 20 A, it is able to deliver only about 78% of the nominal capacity of the battery cell. As the SoH decreases, the increase in the losses on the internal elements of the battery cell, which results in lower voltage values, is present. The battery cell is able to deliver only a smaller capacity. In this case, it is only around 73% of the battery cell’s nominal capacity for the maximum permitted current.
With further reduction of SoH, this trend continues. Losses on the internal elements of the battery cell increase. It confirms Figure 10 for SoH 85% and Figure 11 for SoH 79%.
Figure 12 shows the changes in the shape of the curves at a current of 5 A; it is possible to see how the shape changes, especially in the end positions of the curves, i.e., in SoC = 100% and SoC = 0%. In addition to these measured curves, changes in the battery cell’s capacity based on the currents taken were also measured. A comparison of discharge capacities for different currents can be seen in Figure 13. The cell’s maximal capacity is significantly influenced by the drawn current.
Waveforms in Figure 14 and Figure 15 show the dependence of battery cell ageing depending on charge/discharge currents. The y-axes show the number of cell cycles which the cell is able to perform before reaching SoH = 80%. The value of the SoH was obtained by measuring one charge and discharge cycle with a power analyzer after a certain number of charge and discharge cycles and comparing it to the reference first cycle with SoH = 100%.
Measured data also allows us to describe capacity degradation due to a number of cycles and charging and discharging current. Waveforms in Figure 16 show dependency between capacity degradation and the number of cycles with constant charging current and different discharging currents.
Based on these measured data, all the necessary dependencies were obtained to accurately determine the value of the voltage during discharge at every moment of the cell’s lifetime.

6. Proposed Model

The concept of the battery cell model was designed depending on measured data and parts of the generic battery model from MATLAB. From MATLAB, a part for coulomb counting was used, and other parts of the proposed model were derived based on fitting functions of measured data. We took into account all relationships between these measured data, and the overall concept of the battery cell model is in Figure 17. The proposed model consists of eight blocks with defined relationships.
  • Current Sensor block, which measures the charging and discharging current of the battery cell model.
  • Coulomb Counting and SoC evaluation—these two blocks can be combined into one block. The Coulomb counting block integrates the value of current from the current sensor over time. Subsequently, the current SoC value is determined. These phenomena can be described by the following equations:
    i _ c o u n t = t i d t [ A h ]
    SoC = Q C _ a c t i _ c o u n t Q C _ a c t · 100 [ % ]
Internal connections of this block from MATLAB/Simulink are shown in Figure 18.
  • Average Current calculation block—the block processing the value from the current sensor and distinguishes whether it is a charging or discharging current. The output from this block are two values. The value of Idis_ave corresponds to the average discharge current in the previous cycle, and Ichar_ave corresponds to the charging current. Values of Idis_ave and Ichar_ave are calculated as average integrated current divided by the number of measured samples.
  • Block DoD Calculation—the block calculates and stores the depth of discharge of the cell from the last three measured cycles. It has 4 outputs, three of which give information about three consecutive DoD states from previous cycles. The fourth output is a coefficient determining the number of transitions between the charging and discharging process.
  • Actual capacity calculation block—the block recalculates the currently maximum available capacity of the battery cell based on SoH and current consumption. These values are then fed back to the SoC evaluation block. The calculations in the block are based on real measured data, which can be seen in Table 1.
Values of SoH were defined during measurements. A new cell from the manufacturer was tested with one charge and discharge cycle to determine the initial reference SoH. The SoH of the cell was considered equal to 100% after this first cycle. As initial capacity was considered, the capacity is given by the manufacturer in the datasheet. After that first cycle, all further measurements were taken, and the SoH of cells decreasing during these measurements was captured. Values of SoH for Table 1 were intended by decreasing capacity during measurements compared to the nominal capacity specified by the manufacturer.
All measured data were transferred to a 3D graph, where they were fitted by a MATLAB extension-Curve Fitter (Figure 19). By using this extension, we found the best-fitting functions for measured data. In this case, the founded fitting function is a polynomial of degree 2 when the fitting method used was a nonlinear least square.
All fitting functions mentioned above were selected as the best compromise between accuracy and simplicity for using these functions in the Li-ion battery model.
The final equation for actual cell capacity unify the format.
Q _ a c t = p 00 + p 10 · I + p 01 · SoH + p 20 · I 2 + p 11 · I · SoH + p 02 · SoH 2
Figure 20 shows the connection of the block based on (4). This block has 4 outputs. The first output gives us information about the maximum capacity at the actual load current. The second output represents nominal cell capacity at actual SoH. The remaining two outputs recalculate the capacitance value based on the temperature. Data for determining the temperature equations were obtained from the battery cell datasheet.
  • SoH evaluation block—the main goal of this block is to recalculate the SoH state during the simulation and adjust model parameters. The state of SoH changes during its ageing. The calculation of ageing in the model is based on the ageing coefficient ε , which assumes a value of zero at the beginning of the lifetime and 1 at its end—80% SoH for in our case.
    For determination of the actual SoH, we used [46]:
    ε ( n ) = ε ( n 1 ) + 0.5 N ( n 1 ) · 2 D o D ( n 2 ) + D o D ( n ) D o D ( n 1 )
    Coefficients DoD are obtained from the DoD Calculation block, while the last unknown in this block is the value of the coefficient N, which indicates the maximum estimated number of cycles at the actual testing current and the DoD. It can be described by the following equation [46]:
    N ( n ) = H D o D ( n ) 100 ξ · γ 01 H I d i s _ a v e ( n ) γ 11 · γ 02 H I c h _ a v e ( n ) γ 12
    where:
    -
    H—number of cycles at standard currents and DoD = 100% after which cell achieves 80% SoH
    -
    DoD(n)—DoD value from the previous cycle
    -
    ξ —influence DoD coefficient
    -
    γ 01 , γ 02 —coefficient for value to percentage
    -
    γ 11 —exponent for influence of discharge current
    -
    γ 12 —exponent for influence of charging current
    -
    I d i s _ a v e ( n ) —median of discharging current in previous discharging cycle
    -
    I c h _ a v e ( n ) —median of charging current in previous charging cycle
After transferring these equations into the model, the overall wiring of the ageing estimation block looks like in Figure 21.
  • Block Output voltage calculation—the output of this block is the battery cell voltage. Block is based on dependencies listed in (7).
    V b a t = f ( I , SoC , SoH , Q , Q n o m )
    where:
    -
    I—actual load current (A)
    -
    SoC—state of charge (%)
    -
    SoH—state of health (%)
    -
    Q/ Q n o m —actual/nominal capacity of the battery cell (Ah)
The battery cell voltage accordingly depends on load current, SoC, SoH, actual capacity and nominal capacity. All of these dependencies can be defined by one equation. We obtain this equation by fitting data from Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11. For fitting were tested different polynomial, exponential and Furier-based fitting functions. Some of these functions did not achieve accuracy in SoC > 95% and SoC < 5% areas. As a best final fitting function was selected exponential function of the third degree listed below:
a · e x p ( b · x ) + c · e x p ( d · x ) + e · e x p ( f · x )
This one equation describes every measured current waveform at all selected values of SoH. Subsequently, the individual coefficients of the equation were compared between different current waveforms for one selected SoH. Dependencies between these coefficients were fitted by a polynomial of degree 2, and final equations for individual coefficients we can see below the equation for Ubatt:
U b a t t = a · e x p ( b · S o C ) + c · e x p ( d · S o C ) + e · e x p ( f · S o C )
where
a = a 1 · I 2 + a 2 · I + a 3
b = b 1 · I 2 + b 2 · I + b 3
c = c 1 · I 2 + c 2 · I + c 3
d = d 1 · I 2 + d 2 · I + d 3
e = e 1 · I 2 + e 2 · I + e 3
The value of coefficient f is equal for all current waveforms with constant value f = 100. Therefore, the equation for coefficient f is not described. Equations of coefficients are second-degree polynomial functions Based on these equations of coefficients; it was possible to include another variable in the equations, namely the current. After fitting all the curves, we can obtain the following three-dimensional model of the battery cell for a given SoH, which is dependent on SoC and current (Figure 22 and Figure 23).
After calculating these coefficients for each SoH value, we find out that changes in these coefficients are linearly dependent on each other, and each coefficient is possible to describe with another equation in which SoH will appear as a variable. Equations for description coefficients dependent on SoH are first-degree polynomial functions. As a result, for each coefficient, we obtain an equation that depends on SoC, current and SoH.
a 1 = a 1 a · SoH + a 1 b ; a 2 = a 2 a · SoH + a 2 b
a 3 = a 3 a · SoH + a 3 b ; b 1 = b 1 a · SoH + b 1 b
b 2 = b 2 a · SoH + b 2 b ; b 3 = b 3 a · SoH + b 3 b
c 1 = c 1 a · SoH + c 1 b ; c 2 = c 2 a · SoH + c 2 b
c 3 = c 3 a · SoH + c 3 b ; d 1 = d 1 a · SoH + d 1 b
d 2 = d 2 a · SoH + d 2 b ; d 3 = d 3 a · SoH + d 3 b
e 1 = e 1 a · SoH + e 1 b ; d 2 = e 2 a · SoH + e 2 b
e 3 = e 3 a · SoH + e 3 b ; f = 100
A final simulation model of the battery cell is shown in Figure 24.
The main outputs of the block Output voltage calculation are outputs for controlling the supplies connected at the dynamic sub-model of the battery cell. These supplies can be seen in Figure 25 and provide the correct function of the overall battery model. In addition to these control outputs, the model also has several information or measurement outputs such as the current value of SoC, SoH, current nominal capacity and, of course, the terminal voltage of the simulated cell.

7. Comparison of Model and Real Battery Cell

After completing the digital twin, it was necessary to verify its accuracy by comparing it with the real battery cell. Real data from the battery cell was measured using a programmable DC electronic load (IT8518E from Itech, 60 V, 240 A, 6 kW) and a high-precision power analyzer (LMG500 from ZES Zimmer) which were also used for cycling tests. Several measurements took place with a constant value of current. These measurements were static, so the value of the current was constant. Results for different values of discharging current are in Figure 26a. Besides that, several dynamic measurements were performed with a variable current value to determine the accuracy of the model. Some of these measurements were selected for the results section to demonstrate the accuracy of the digital twin. In addition to the data from the real battery cell and the presented empirical model, a battery cell model available in MATLAB/Simulink, called the Generic Battery model, was also added to the comparison. This model was configured according to the manual in MATLAB. A comparison of dynamics of the real battery cell, the developed digital twin and the MATLAB model is in Figure 26b.
The digital twin utilises the Coulomb counting method for SoC estimation, which is prone to integration error due to the current measurement error and simplified integration. Moreover, this error will accumulate over time. Despite these disadvantages, Coulomb counting is still used due to its simplicity. The influence of the current measurement can be minimised by an open-loop hall-effect sensor, which is superior when compared to the shunt current sensing. A detailed discussion and analysis of the Coulomb counting approach for SoC estimation can be found in [47]. The SoH is estimated based on the DoD, and, thus, it contains the error generated by SoC estimation. It requires calibration in regular intervals. Another approach to improving the SoC is presented in [48]. The authors used the Kalman filter to correct the estimation of SoC by the Kalman filter through the feedback based on the relationship between SoC and OCV. The authors achieved the SoC estimation absolute error of 1.08%.

8. Conclusions

Battery cell models are a crucial part of any BMS. Successful implementation of a battery cell model into MATLAB/Simulink can help with the development of control methods for cell lifetime optimization. An empirical model of a battery cell can be easily generated from measured data obtained by the battery cell cycling by utilising the curve fitting functions. The empirical SoH model is based on data obtained from cell cycling up to SoH equal to 80% for the NR18650-25R cell. The data for empirical model creation, thus, contains information about the whole life cycle. Hence, the MATLAB/Simulink digital twin of the Li-ion cell from Samsung was developed based on the measured data. The data from battery cell cycling include fast charging, which is well known to accelerate the ageing process of the battery cell. However, it is not always included in publicly available battery cell data sets. The created digital twin provides information about SoC, SoH, cell voltage and cell capacity based only on the cell current. The model was successfully tested and compared to the Generic Simulink model of a Li-ion cell based on the data provided by the manufacturer. From the comparison, one can conclude that the static accuracy of the developed digital twin is better than the accuracy generic Simulink model. However, data used for digital twin creation from cycling tests are obtained by galvanostatic discharging. This leads to less accurate model performance for dynamic current changes. The digital twin could be improved by adding information to the charging and discharging cycles by implementing Hybrid Pulse Power Characterization (HPPC) and micro-cycling. It would enhance the data sets by the dynamical response of the battery cell, which would lead to more realistic parameter estimation. Moreover, the model could be extended to the battery pack. However, this would require a cell balancing technique to prevent the mismatch between cells, as the only parameter that is sensed is the battery current.

Author Contributions

Conceptualization, J.B., M.L. and F.D.; methodology, J.B.; software, J.B.; formal analysis, J.B.; investigation, J.B.; data curation, J.B.; writing—original draft preparation, M.P. and A.M.; writing—review and editing, M.P. and A.M.; visualization, A.M.; supervision, M.L.; project administration, F.D.; funding acquisition, M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Slovak Research and Development Agency under project APVV-18-0436. This work was supported by Slovak Research and Development Agency under project APVV-19-0210.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Special thanks are given to the support by T. Merva for developing data collecting software.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Council: Fit for 55. Available online: https://www.consilium.europa.eu/en/policies/green-deal/fit-for-55-the-eu-plan-for-a-green-transition (accessed on 16 November 2022).
  2. European Council: Infographic-Fit for 55: Why the EU Is Toughening CO2 Emission Standards for Cars and Vans. Available online: https://www.consilium.europa.eu/en/infographics/fit-for-55-emissions-cars-and-vans (accessed on 16 November 2022).
  3. Vermeer, W.; Mouli, G.R.C.; Bauer, P. A Comprehensive Review on the Characteristics and Modeling of Lithium-Ion Battery Aging. IEEE Trans. Transp. Electrif. 2022, 8, 2205–2232. [Google Scholar] [CrossRef]
  4. Jia, J.; Liang, J.; Shi, Y.; Wen, J.; Pang, X.; Zeng, J. SoH and RUL Prediction of Lithium-Ion Batteries Based on Gaussian Process Regression with Indirect Health Indicators. Energies 2020, 13, 375. [Google Scholar] [CrossRef] [Green Version]
  5. Watrin, N.; Blunier, B.; AMiraoui, A. Review of adaptive systems for lithium batteries State-of-Charge and State-of-Health estimation. In Proceedings of the IEEE Transportation Electrification Conference and Expo (ITEC), Dearborn, MI, USA, 18–20 June 2012. [Google Scholar] [CrossRef]
  6. Kim, J.; Kowal, J. Development of a Matlab/Simulink Model for Monitoring Cell State-of-Health and State-of-Charge via Impedance of Lithium-Ion Battery Cells. Batteries 2022, 8, 8. [Google Scholar] [CrossRef]
  7. Theiler, M.; Schneider, S.; Endisch, C. Kalman Filter Tuning Using Multi-Objective Genetic Algorithm for State and Parameter Estimation of Lithium-Ion Cells. Batteries 2022, 8, 104. [Google Scholar] [CrossRef]
  8. Mao, S.; Han, M.; Han, X.; Lu, L.; Feng, X.; Su, A.; Wang, D.; Chen, Z.; Lu, Y.; Ouyang, M. An Electrical–Thermal Coupling Model with Artificial Intelligence for State of Charge and Residual Available Energy Co-Estimation of LiFePO4 Battery System under Various Temperatures. Batteries 2022, 8, 140. [Google Scholar] [CrossRef]
  9. Deng, Z.; Hu, X.; Lin, X.; Che, Y.; Xu, L.; Guo, W. Data-driven state of charge estimation for lithium-ion battery packs based on Gaussian process regression. Energy 2020, 205, 118000. [Google Scholar] [CrossRef]
  10. Dickson, N.T.; Hannah, M.A.; Lipu Hossain, M.S.; Ker, P.J. State of Charge Estimation for Lithium-ion Batteries Using Model-Based and Data-Driven Methods: A review. IEEE Access 2022, 7, 136116–136136. [Google Scholar] [CrossRef]
  11. Du, J.; Liu, Z.; Wang, Y.; Wen, C. A Fuzzy Logic-based Model for Li-ion Battery with SoC and Temperature Effect. In Proceedings of the 11th IEEE International Conference on Control and Automation (ICCA), Taichung, Taiwan, 18–20 June 2014. [Google Scholar] [CrossRef]
  12. Kuchly, J.; Goussian, A.; Merveillaut, M.; Baghdadi, I.; Franger, S.; Nelson-Gruel, D.; Nouillant, C.; Chamaillard, Y. Li-ion battery SoC estimation method using a Neural Network trained with data generated by a P2D model. IFAC-PapersOnLine 2021, 54, 336–343. [Google Scholar] [CrossRef]
  13. Cui, Z.; Hu, W.; Zhang, G.; Zhang, Z.; Chen, Z. An extended Kalman filter based SoC estimation method for Li-ion battery. Energy Rep. 2022, 8, 81–87. [Google Scholar] [CrossRef]
  14. Birkl, C. Oxford Battery Degradation Dataset 1; University of Oxford: Oxford, UK, 2017; Available online: https://ora.ox.ac.uk/objects/uuid:03ba4b01-cfed-46d3-9b1a-7d4a7bdf6fac (accessed on 25 October 2022).
  15. Kollmeyer, P.; Vidal, C.; Naguib, M.; Skells, M. LG 18650HG2 Li-Ion Battery Data and Example Deep Neural Network xEV SoC Estimator Script, Mendeley Data, V3. 2020. Available online: https://data.mendeley.com/datasets/cp3473x7xv/3/ (accessed on 15 October 2022). [CrossRef]
  16. Kollmeyer, P. Panasonic 18650PF Li-ion Battery Data, Mendeley Data, V1. 2018. Available online: https://data.mendeley.com/datasets/wykht8y7tg/1/ (accessed on 15 October 2022). [CrossRef]
  17. Severson, K.A.; Attia, P.M.; Jin, N.; Perkins, N.; Jiang, B.; Yang, Z.; Chen, M.H.; Aykol, M.; Herring, P.K.; Fraggedakis, D.; et al. Data-driven prediction of battery cycle life before capacity degradation. Nat. Energy 2019, 4, 383–391. [Google Scholar] [CrossRef] [Green Version]
  18. Omar, N.; Monem, M.A.; Firouz, Y.; Salminen, J.; Smekens, J.; Hegazy, O.; Gaulous, H.; Mulder, G.; Van den Bossche, P.; Coosemans, T.; et al. Lithium iron phosphate based battery—Assessment of the aging parameters and development of cycle life model. Appl. Energy 2014, 113, 1575–1585. [Google Scholar] [CrossRef]
  19. Weiping, D.; Saurabh, S.; Michael, P. Accelerated cycle life testing and capacity degradation modeling of LiCoO2-graphite cells. J. Power Sources 2019, 453, 1575–1585. [Google Scholar] [CrossRef]
  20. Singh, K.V.; Bansal, H.O.; Singh, D. Hardware-in-the-loop Implementation of ANFIS based Adaptive SoC Estimation of Lithium-ion Battery for Hybrid Vehicle Applications. J. Energy Storage 2020, 27, 101124. [Google Scholar] [CrossRef]
  21. Khanum, F.; Louback, E.; Duperly, F.; Jenkins, C.; Kollmeyer, P.J.; Emadi, A. A Kalman Filter Based Battery State of Charge Estimation MATLAB Function. In Proceedings of the 2021 IEEE Transportation Electrification Conference and Expo (ITEC), Chicago, IL, USA, 21–25 June 2021; pp. 484–489. [Google Scholar] [CrossRef]
  22. Weng, C.; Sun, J.; Peng, H. A unified open-circuit-voltage model of lithium-ion batteries for state-of-charge estimation and state-of-health monitoring. J. Power Sources 2014, 258, 228–237. [Google Scholar] [CrossRef]
  23. Surya, S.; Saldanha, C.C.; Williamson, S. Novel Technique for Estimation of Cell Parameters Using MATLAB/Simulink. Electronics 2022, 11, 117. [Google Scholar] [CrossRef]
  24. Kim, J.; Chun, H.; Kim, M.; Yu, J.; Kim, K.; Kim, T.; Han, S. Data-Driven State of Health Estimation of Li-Ion Batteries With RPT-Reduced Experimental Data. IEEE Access 2019, 7, 106987–106997. [Google Scholar] [CrossRef]
  25. Guy, B. A History of the Battery. Available online: https://batteryguy.com/kb/knowledge-base/a-history-of-the-battery (accessed on 5 October 2022).
  26. Arm, M.; Axmann, P.; Bresser, D.; Copley, M.; Edström, K.; Ekberg, C.; Guyomard, D.; Lestriez, B.; Novák, P.; Petranikova, M.; et al. Lithium-ion batteries—Current state of the art and anticipated developments. J. Power Sources 2020, 479, 228708. [Google Scholar] [CrossRef]
  27. R-Smith, N.A.Z.; Gramse, G.; Moertelmaie, R.M.; Kasper, M.; Kienberge, R.F. Advanced Self-Discharge Measurements of Lithium-Ion Cells and Comparison to Modeling. In Proceedings of the IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Dubrovnik, Croatia, 25–28 May 2020; pp. 1–5. [Google Scholar] [CrossRef]
  28. Redondo-Iglesias, E.; Venet, P.; Pelissier, S. Measuring Reversible and Irreversible Capacity Losses on Lithium-Ion Batteries. In Proceedings of the 2016 IEEE Vehicle Power and Propulsion Conference (VPPC), Hangzhou, China, 17–20 October 2016; pp. 1–5. [Google Scholar] [CrossRef] [Green Version]
  29. Kraft, L.; Zünd, T.; Schreiner, D.; Wilhelm, R.; Günter, F.J.; Reinhart, G.; Gasteiger, H.A.; Jossen, A. Comparative Evaluation of LMR-NCM and NCA Cathode Active Materials in Multilayer Lithium-Ion Pouch Cells: Part II. Rate Capability, Long-Term Stability, and Thermal Behavior. J. Electrochem. Soc. 2021, 168, 020537. [Google Scholar] [CrossRef]
  30. Battery University. Bu-205: Types of Lithium-Ion. Available online: https://batteryuniversity.com/learn/article/types_of_lithium_ion (accessed on 7 October 2022).
  31. Battery University. What’s the Best Battery? Available online: https://batteryuniversity.com/learn/archive/whats_the_best_battery (accessed on 7 October 2022).
  32. Aditya, J.P.; Ferdowsi, M. Comparison of NiMh and li-ion batteries in automotive applications. In Proceedings of the 2008 IEEE Vehicle Power and Propulsion Conference, Harbin, China, 3–5 September 2008; pp. 1–6. [Google Scholar] [CrossRef]
  33. Meena, N.; Baharwani, V.; Sharma, D.; Sharma, A.; Choudhary, B.; Parmar, P.; Stephen, R.B. Charging and discharging characteristics of lead acid and li-ion batteries. In Proceedings of the 2014 Power and Energy Systems: Towards Sustainable Energy, Bangalore, India, 13–15 March 2014; pp. 1–3. [Google Scholar] [CrossRef]
  34. Wey, C.L.; Jui, P.C. A unitized charging and discharging smart battery management system. In Proceedings of the 2013 International Conference on Connected Vehicles and Expo (ICCVE), Las Vegas, NV, USA, 2–6 December 2013; pp. 903–909. [Google Scholar] [CrossRef]
  35. Khan, A.B.; Pham, V.L.; Nguyen, T.T.; Choi, W. Multistage constant-current charging method for li-ion batteries. In Proceedings of the 2016 IEEE Transportation Electrification Conference and Expo, Asia-Pacific (ITEC Asia-Pacific), Busan, Republic of Korea, 1–4 June 2016; pp. 381–385. [Google Scholar] [CrossRef]
  36. Luan, S.W.; Teng, J.H.; Lee, D.J.; Huang, Y.Q.; Sung, C.L. Charging/discharging monitoring and simulation platform for li-ion batteries. In Proceedings of the TENCON 2011-2011 IEEE Region 10 Conference, Bali, Indonesia, 21–24 November 2011; pp. 868–872. [Google Scholar] [CrossRef]
  37. Keil, P.; Jossen, A. Charging protocols for lithium-ion batteries and their impact on cycle life—An experimental study with different 18650 high-power cells. J. Energy Storage 2016, 6, 125–141. [Google Scholar] [CrossRef]
  38. Notten, P.H.L.; Op het Veld, J.H.G.; van Beek, J.R.G. Boostcharging Li-ion batteries: A challenging new charging concept. J. Power Sources 2005, 145, 89–94. [Google Scholar] [CrossRef]
  39. Tomaszewska, A.; Chu, Z.; Feng, X.; O’kane, S.; Liu, X.; Chen, J.; Ji, C.; Endler, E.; Li, R.; Liu, L.; et al. Lithium-ion battery fast charging: A review. eTransportation 2019, 1, 28. [Google Scholar] [CrossRef]
  40. Wan, H. High Efficiency dc-dc Converter for ev Battery Charger Using Hybrid Resonant and pwm Technique. Ph.D. Thesis, Virginia Tech, Blacksburg, VA, USA, 2012. Available online: https://vtechworks.lib.vt.edu/handle/10919/32343 (accessed on 5 September 2022).
  41. Simolka, M.; Heger, J.F.; Kaess, H.; Biswas, I.; Friedrich, K.A. Influence of cycling profile, depth of discharge and temperature on commercial LFP/C cell ageing: Post-mortem material analysis of structure, morphology and chemical composition. J. Appl. Electrochem. 2020, 1, 1101–1117. [Google Scholar] [CrossRef]
  42. Soto, A.; Alberto Berrueta, A.; Mateos, M.; Sanchis, P.; Ursúa, A. Impact of micro-cycles on the lifetime of lithium-ion batteries: An experimental study. J. Energy Storage 2022, 55, 105343. [Google Scholar] [CrossRef]
  43. PowerStream. Introduction of INR18650-25R. Available online: https://www.powerstream.com/p/INR18650-25R-datasheet.pdf (accessed on 8 February 2022).
  44. Biľanský, J.; Lacko, M. Design and simulation of cyclic battery tester. Power Electron. Drives 2020, 5, 13. [Google Scholar] [CrossRef]
  45. Biľanský, J.; Merva, T.; Ivan, J.; Marcinek, A.; Lacko, M. Cyclic tester of battery cells for electric vehicles. In Proceedings of the 2021 IEEE International Workshop of Electronics, Control, Measurement, Signals and Their Application to Mechatronics (ECMSM), Liberec, Czech Republic, 21–22 June 2021; pp. 1–7. [Google Scholar] [CrossRef]
  46. Motapon, S.N.; Lachance, E.; Dessaint, L.A.; AL-Haddad, K. A Generic Cycle Life Model for Lithium-Ion Batteries Based on Fatigue Theory and Equivalent Cycle Counting. IEEE Open J. Ind. Electron. Soc. 2020, 1, 207–217. [Google Scholar] [CrossRef]
  47. Movassagh, K.; Raihan, A.; Balasingam, B.; Pattipati, K. A Critical Look at Coulomb Counting Approach for State of Charge Estimation in Batteries. Energies 2019, 14, 4074. [Google Scholar] [CrossRef]
  48. Fang, L.; Li, J.; Peng, B. Online Estimation and Error Analysis of both SOC and SOH of Lithium-ion Battery based on DEKF Method. Energy Procedia 2019, 158, 3008–3013. [Google Scholar] [CrossRef]
Figure 1. (a) comparison of various Li-ion cells [25], (b) cell used for cycling testing.
Figure 1. (a) comparison of various Li-ion cells [25], (b) cell used for cycling testing.
Energies 16 01194 g001
Figure 2. Constant voltage/Constant current charging waveform [40]. Courtesy of CADEX.
Figure 2. Constant voltage/Constant current charging waveform [40]. Courtesy of CADEX.
Energies 16 01194 g002
Figure 3. Discharging characteristics for various discharge currents [43]. Courtesy of SAMSUNG.
Figure 3. Discharging characteristics for various discharge currents [43]. Courtesy of SAMSUNG.
Energies 16 01194 g003
Figure 4. Testing system topology.
Figure 4. Testing system topology.
Energies 16 01194 g004
Figure 5. Developed device for cycling battery cell testing.
Figure 5. Developed device for cycling battery cell testing.
Energies 16 01194 g005
Figure 6. Parameters of the Samsung INR18650-25R battery cell [43]. Courtesy of SAMSUNG.
Figure 6. Parameters of the Samsung INR18650-25R battery cell [43]. Courtesy of SAMSUNG.
Energies 16 01194 g006
Figure 7. Discharge characteristics for SoH = 100%.
Figure 7. Discharge characteristics for SoH = 100%.
Energies 16 01194 g007
Figure 8. Discharge characteristics for SoH = 94%.
Figure 8. Discharge characteristics for SoH = 94%.
Energies 16 01194 g008
Figure 9. Discharge characteristics for SoH = 89%.
Figure 9. Discharge characteristics for SoH = 89%.
Energies 16 01194 g009
Figure 10. Discharge characteristics for SoH = 85%.
Figure 10. Discharge characteristics for SoH = 85%.
Energies 16 01194 g010
Figure 11. Discharge characteristics for SoH = 79%.
Figure 11. Discharge characteristics for SoH = 79%.
Energies 16 01194 g011
Figure 12. Discharge characteristics for 5 A a SoH = 100%, 94%, 89%, 85%, 79%.
Figure 12. Discharge characteristics for 5 A a SoH = 100%, 94%, 89%, 85%, 79%.
Energies 16 01194 g012
Figure 13. Influence of discharge current and SoH on maximal capacity.
Figure 13. Influence of discharge current and SoH on maximal capacity.
Energies 16 01194 g013
Figure 14. Influence of the charging current on the maximal number of cycles before the battery cell reaches SoH = 80%.
Figure 14. Influence of the charging current on the maximal number of cycles before the battery cell reaches SoH = 80%.
Energies 16 01194 g014
Figure 15. Influence of the discharge current on the maximal number of cycles before the battery cell reaches SoH = 80%.
Figure 15. Influence of the discharge current on the maximal number of cycles before the battery cell reaches SoH = 80%.
Energies 16 01194 g015
Figure 16. Influence of the number of cycles and discharging current on capacity degradation.
Figure 16. Influence of the number of cycles and discharging current on capacity degradation.
Energies 16 01194 g016
Figure 17. Concept of the proposed model.
Figure 17. Concept of the proposed model.
Energies 16 01194 g017
Figure 18. Coulomb counting and SoC evaluation. Where Q_act parameter computing is based on the equation described below in the Actual capacity calculation block section, i is the actual sensed current.
Figure 18. Coulomb counting and SoC evaluation. Where Q_act parameter computing is based on the equation described below in the Actual capacity calculation block section, i is the actual sensed current.
Energies 16 01194 g018
Figure 19. Fitted date based on measurements from Table 1.
Figure 19. Fitted date based on measurements from Table 1.
Energies 16 01194 g019
Figure 20. Block for calculating capacity.
Figure 20. Block for calculating capacity.
Energies 16 01194 g020
Figure 21. Block for approximation cell ageing.
Figure 21. Block for approximation cell ageing.
Energies 16 01194 g021
Figure 22. Data fitted by the third-degree exponential function for SoH = 100%.
Figure 22. Data fitted by the third-degree exponential function for SoH = 100%.
Energies 16 01194 g022
Figure 23. Data fitted by the third-degree exponential function of for SoH = 89%.
Figure 23. Data fitted by the third-degree exponential function of for SoH = 89%.
Energies 16 01194 g023
Figure 24. Interconnection of an individual block within the model.
Figure 24. Interconnection of an individual block within the model.
Energies 16 01194 g024
Figure 25. Final digital twin based on an empirical model of the Li-ion cell.
Figure 25. Final digital twin based on an empirical model of the Li-ion cell.
Energies 16 01194 g025
Figure 26. Comparison of (a) static measurement and (b) dynamic measurements for SoH = 100%.
Figure 26. Comparison of (a) static measurement and (b) dynamic measurements for SoH = 100%.
Energies 16 01194 g026
Table 1. Battery cell capacity at different discharge current and different SoH.
Table 1. Battery cell capacity at different discharge current and different SoH.
SoH (%)10094898579
1 A2.5 Ah2.35 Ah2.22 Ah2.17 Ah2.05 Ah
2.5 A2.48 Ah2.30 Ah2.19 Ah2.13 Ah1.98 Ah
5 A2.45 Ah2.28 Ah2.17 Ah1.99 Ah1.88 Ah
10 A2.43 Ah2.12 Ah2.00 Ah1.92 Ah1.79 Ah
15 A2.42 Ah2.04 Ah1.96 Ah1.85 Ah1.74 Ah
20A2.38 Ah1.96 Ah1.85 Ah1.74 Ah1.64 Ah
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bilansky, J.; Lacko, M.; Pastor, M.; Marcinek, A.; Durovsky, F. Improved Digital Twin of Li-Ion Battery Based on Generic MATLAB Model. Energies 2023, 16, 1194. https://doi.org/10.3390/en16031194

AMA Style

Bilansky J, Lacko M, Pastor M, Marcinek A, Durovsky F. Improved Digital Twin of Li-Ion Battery Based on Generic MATLAB Model. Energies. 2023; 16(3):1194. https://doi.org/10.3390/en16031194

Chicago/Turabian Style

Bilansky, Juraj, Milan Lacko, Marek Pastor, Adrian Marcinek, and Frantisek Durovsky. 2023. "Improved Digital Twin of Li-Ion Battery Based on Generic MATLAB Model" Energies 16, no. 3: 1194. https://doi.org/10.3390/en16031194

APA Style

Bilansky, J., Lacko, M., Pastor, M., Marcinek, A., & Durovsky, F. (2023). Improved Digital Twin of Li-Ion Battery Based on Generic MATLAB Model. Energies, 16(3), 1194. https://doi.org/10.3390/en16031194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop