Multi-Cell-to-Multi-Cell Battery Equalization in Series Battery Packs Based on Variable Duty Cycle

: Batteries are widely used in our lives, but the inevitable inconsistencies in series-connected battery packs will seriously impact their energy utilization, cycle life and even jeopardize their safety in use. This paper proposes a balancing topology structure combining Buck-Boost circuit and switch array to reduce this inconsistency. This structure can realize multi-cell-to-multi-cell (MC2MC) battery balancing by controlling the switch array and having a fast balancing speed, easy expansion and few magnetic components. Then, the operation principle of the proposed balancing topology is analyzed, and the simulation model is veriﬁed. In addition, the effects of switching frequency and voltage difference on the equalization effect are further analyzed. The results show that the higher the switching frequency, the lower the time efﬁciency, but the higher the energy efﬁciency. The voltage difference signiﬁcantly impacts the duty cycle, so it is absolutely necessary to introduce a variable duty cycle in the multi-cell-to-multi-cell equalization. Finally, eight series batteries are selected for simulation veriﬁcation. The simulation results show that, compared with any-cell-to-any-cell (AC2AC) equalization, the time efﬁciency of multi-cell-to-multi-cell equalization is improved considerably, the energy efﬁciency is improved slightly, and the variance of the completed equalization is reduced, demonstrating the excellent performance of multi-cell-to-multi-cell equalization. inductor L switching frequency f the discharging and loop resistance R R the above Equation (4) duty cycle of D = 0.506. calculated of inductor above parameter Equations (1) (2). The simulated value of inductor the the equivalent The


Introduction
Considering energy security, environmental protection and technological progress, industries such as the global power grid, manufacturing and transportation have gradually turned to low energy consumption, low pollution and low carbon [1][2][3]. Lithium-ion batteries are widely used in electric vehicles, mobile communications, aerospace and energy storage systems [4][5][6] due to their high energy density, low self-discharge rate, long cycle life, small size, and lightweight [7]. In order to meet the high voltage, high power and large capacity requirements of electric vehicles and energy storage systems, batteries are often used in groups in series or parallel [8]. However, due to the complexity of lithium-ion battery manufacturing technology and production process, the capacity, voltage, internal resistance and other performance parameters of lithium-ion batteries produced in the same batch will still be inconsistent [9]. This inconsistency will be further deepened during use and storage [10]. According to the "bucket theory", the maximum charging capacity of the battery pack is limited by the highest SOC battery, and the maximum discharging capacity is limited by the smallest SOC battery. This makes the actual usable capacity of the battery pack smaller than the theoretical maximum usable capacity [11]. This phenomenon persists throughout the life cycle of the battery pack.
In order to reduce the energy difference between the individual cells of the battery pack, the battery balancing technology came into being. With the refinement of battery modeling theory and the improvement of SOC estimation methods, the accuracy of battery models is becoming higher and higher [12,13], which provides a solid basis for research in areas such as battery equalization. The battery balancing system mainly includes balancing topology and control strategy [14]. The balancing topology is related to the dissipation or transfer of the balancing energy. Improving the topology can effectively improve the speed and efficiency of battery balancing, reduce circuit losses and lower balancing costs. There are various forms of balancing topology. According to whether the energy is dissipated in the balancing process, it can be divided into passive balancing (also called dissipative balancing) and active balancing (also called non-dissipative balancing) [15]. According to the energy transfer path in the equalization process, it can be divided into adjacent cell-to-cell (AC2C), pack-to-cell (P2C), cell-to-pack (C2P), any-cell-to-any-cell (AC2AC) and multi-cell-to-multi-cell (MC2MC) equalization methods [16].
Passive equalization refers to the technology of using parallel resistance to consume the batteries with more energy or higher voltage in the battery pack until the energy of each battery tends to be consistent [17]. Although passive equalization is easy to design and implement, it has the disadvantages of serious resistance heating effect, long equalization time and low energy utilization rate [18].
Pascual et al. [20] first proposed a multi-switched capacitor structure, but its balancing energy could only be transmitted between adjacent batteries through capacitors. Lee et al. [21] improved the switched capacitor equalization scheme and used zero-current switched capacitors to improve the equalization efficiency. Daowd et al. [22] proposed a single switched capacitor structure with a single capacitor as the energy storage element, and equalization energy could be transferred between any two cells. However, the capacitorbased balancing topology can only ensure the same voltage between batteries but cannot ensure the same energy. The balancing effect is poor when the battery voltage difference is not obvious. Thus, it is only suitable for occasions that do not require high equalization speed and accuracy.
In order to ensure the reliability of the balancing circuit under high voltage and high current conditions, transformer-based isolated balancing topologies are gaining more and more attention. It uses the transformer as the energy transfer carrier and transfers energy through mutual inductance. Li et al. [23] proposed a coaxial multi-winding transformer topology, but it could only realize the one-way transfer of balancing energy. Kutkut et al. [24] proposed an improved scheme for multi-secondary transformers, which could reduce the number of transformer secondary sides by half. Liu et al. [25] combined the switch array with the forward transformer and reduced conduction losses through synchronous rectification. But the magnetic reset circuit designed to prevent transformer core saturation complicated the topology and reduced reliability. Shang et al. [26] introduced forward transformation and flyback transformation into a multi-winding transformer simultaneously. The magnetic energy in the transformer at the end of each cycle was eliminated by the flyback transformation, which effectively solved the mismatch problem of the multiwinding transformer. In order to optimize the equalizer structure, Shang et al. [27] reduced the number of transformer windings by half by introducing a coupled half-bridge converter based on Reference [26], but it complicated the control process. Transformer-based equalization topologies are bulky, expensive, and difficult to expand; it is only suitable for occasions with few series batteries.
Using DC/DC converter circuits, such as the Cuk circuit, Boost circuit and Buck-Boost circuit to achieve equalization is called converter-based battery equalization, which also uses capacitors, inductors or inductors or transformers as energy storage elements. Daowd et al. [28] compared a variety of active and passive equalization schemes in MAT-LAB/Simulink. They pointed out that the DC/DC converter equalization had high speed and high-efficiency advantages. Das et al. [29] proposed a distributed bidirectional active equalization structure based on a Buck-Boost converter. This structure used inductors as energy storage elements, but equalization energy could only be transferred between adjacent batteries in sequence. In [30], the first and last batteries were also connected with a basic Buck-Boost topology unit so that the batteries in series could form a circular energy loop, which increased the flexibility of balancing energy transfer. Wu et al. [31] proposed a hierarchical equalization topology based on the Buck-Boost converter. The equalization process was divided into multiple equalization processes of different layers according to the equalization energy flow. The equalization processes between different layers did not affect each other and could occur simultaneously to effectively speed up the equalization speed. Wang et al. [32] proposed an equalization topology that could construct multiple Buck-Boost converters between adjacent cells or adjacent cell clusters. The equalization path of this structure was flexible, and different equalization processes could occur simultaneously, so output power could be significantly improved by forming cell groups. However, the deficiency of reference [30][31][32] is that the equalization energy can only be transmitted successively between adjacent batteries (clusters). Chen et al. [33] proposed a bidirectional active equalization method based on a single inductor, which allowed equalization energy to be directly transferred from any cell to any cell. To optimize the equalization process, the influence of various factors on the equalization current and efficiency is then analyzed. Wu et al. [34] adopted a balancing topology combining a Cuk equalizer and double-layer selector switches, which could also achieve energy transfer between any battery. Although this topology reduced the number of inductors and capacitors exponentially, it required many switches. In order to equalize multiple cells simultaneously in one equalization process, Shang et al. [35] proposed a direct multi-cell-to-multi-cell equalization topology based on the matrix LC converter or the half-bridge LC converter to realize the direct energy transfer from continuous high-energy unit group to continuous low-energy unit group. However, this topology had a complex structure and required many switches. Converter-based equalization topology has low loss, high speed and good scalability, and it is one of the key development directions of active equalization.
In summary, the capacitor-based equalization topology is only suitable for occasions with low equalization speed and accuracy requirements. The inductance-based equalization topology is often combined with transformers and converters. The transformer-based equalization topology has a large volume, high cost and is not easy to expand. Therefore, this paper designs the balancing topology based on the Buck-Boost converter. The topology structure is simplified by introducing the switch array while realizing MC2MC equalization. According to the analysis results of the factors affecting the equalization effect, an equalization control strategy based on variable duty ratio is proposed. Finally, the simulation verification of eight series battery packs is carried out.
The structure of this paper is as follows. The second part introduces the topology and operating principle of MC2MC equalization based on the Buck-Boost circuit. The third part analyzes the factors that affect the equalization effect. The fourth part uses Matlab/Simulink to model and simulate the proposed topology and control strategy based on a variable duty cycle. Finally, the fifth section gives the relevant conclusions.

Balancing Topology
This paper uses the Buck-Boost circuit combined with the switching array to design the equalization topology, which enables the equalization of the multi-cell-to-multi-cell energy transfer by controlling the switching array. The overall balancing topology is shown in Figure 1. The main structure of this topology is the switching array and the balancing main circuit. The switch array consists of positive switch S11, S21, S31~Sn1 and negative switch S12, S22, S32~Sn2. All positive switches are connected together to form a positive bus, and all negative switches are connected to form a negative bus. The main balancing circuit adopts the non-isolated Buck-Boost circuit, and its circuit elements include inductance L and switching devices Q1, Q2, Q3 and Q4. The circuit operates under continuous current mode (CCM) to ensure that more energy is transferred in one switching cycle. energy transfer by controlling the switching array. The overall balancing topology is shown in Figure 1. The main structure of this topology is the switching array and the balancing main circuit. The switch array consists of positive switch S11, S21, S31~Sn1 and negative switch S12, S22, S32~Sn2. All positive switches are connected together to form a positive bus, and all negative switches are connected to form a negative bus. The main balancing circuit adopts the non-isolated Buck-Boost circuit, and its circuit elements include inductance L and switching devices Q1, Q2, Q3 and Q4. The circuit operates under continuous current mode (CCM) to ensure that more energy is transferred in one switching cycle. The proposed balancing topology has the following advantages.
(1) Multiple batteries can be equalized simultaneously in one equalization process, which greatly reduces the equalization time.
(2) Only additional positive and negative switches are needed when the number of series batteries increases. Thus, the circuit complexity remains unchanged; that is, the circuit is easy to expand. (3) The number of magnetic components, such as inductance and capacitance, is greatly reduced.

Principle of Operation
The proposed topology has three typical working modes, namely any-cell-to-any-cell equalization (Mode 1), multi-cell-to-any-cell equalization (Mode 2) and multi-cell-tomulti-cell equalization (Mode 3). This paper takes eight series batteries as examples to analyze their working principle.
Mode 1: Assuming that the source battery B1 balances with the target battery B5, as shown in Figure 2, this process is divided into two steps.
1. Source battery B1 discharging: The discharging loop is composed of source battery B1, positive and negative switches S11, S12, MOS transistors Q1, Q2 and inductor L. The current flows along the red arrow. At this time, the electric energy of the source battery B1 is converted into the magnetic energy of the inductor L.  The proposed balancing topology has the following advantages.
(1) Multiple batteries can be equalized simultaneously in one equalization process, which greatly reduces the equalization time.
(2) Only additional positive and negative switches are needed when the number of series batteries increases. Thus, the circuit complexity remains unchanged; that is, the circuit is easy to expand. (3) The number of magnetic components, such as inductance and capacitance, is greatly reduced.

Principle of Operation
The proposed topology has three typical working modes, namely any-cell-to-any-cell equalization (Mode 1), multi-cell-to-any-cell equalization (Mode 2) and multi-cell-to-multicell equalization (Mode 3). This paper takes eight series batteries as examples to analyze their working principle.
Mode 1: Assuming that the source battery B1 balances with the target battery B5, as shown in Figure 2, this process is divided into two steps. Mode 2: Assuming that adjacent batteries B2 and B3 form a high-energy battery cluster to balance the target battery B5. As shown in Figure 3, this process is divided into the following two steps.  Source battery B1 discharging: The discharging loop is composed of source battery B1, positive and negative switches S11, S12, MOS transistors Q1, Q2 and inductor L. The current flows along the red arrow. At this time, the electric energy of the source battery B1 is converted into the magnetic energy of the inductor L.

2.
Target battery B5 charging: The charging loop consists of target battery B5, positive and negative switches S51, S52, MOS transistors Q3, Q4 and inductor L. The current flows along the blue arrow. At this time, the magnetic energy stored in the inductor L is converted into electric energy in the target battery B5.
Mode 2: Assuming that adjacent batteries B2 and B3 form a high-energy battery cluster to balance the target battery B5. As shown in Figure 3, this process is divided into the following two steps. Mode 2: Assuming that adjacent batteries B2 and B3 form a high-energy battery cluster to balance the target battery B5. As shown in Figure 3, this process is divided into the following two steps.
1. Battery clusters B2, B3 discharging: The discharging loop is composed of battery clusters B2 and B3, positive and negative switches S21 and S32, MOS transistors Q1, Q2 and inductor L. The current flows along the red arrow. At this time, the electrical energy of the high-energy battery cluster B2 and B3 is converted into the magnetic energy of the inductor L.  Mode 3: Assuming that high-energy battery clusters B2 and B3 balance with lowenergy battery clusters B6 and B7. As shown in Figure 4, this process is divided into the following two steps.
1. Battery clusters B2, B3 discharging: The discharging loop is composed of battery clusters B2 and B3, positive and negative switches S21, S32, MOS transistors Q1, Q2 and inductor L. The current flows along the red arrow. At this time, the electrical energy of the high-energy battery cluster B2 and B3 is converted into the magnetic energy of the inductor L. 2. Battery clusters B6, B7 charging: The charging loop consists of battery clusters B6 and B7, positive and negative switches S61, S72, MOS transistors Q3, Q4 and inductance L. The current flows along the blue arrow. At this time, the magnetic energy stored in the inductor L is converted into electrical energy in the low-energy battery cluster B6 and B7.

1.
Battery clusters B2, B3 discharging: The discharging loop is composed of battery clusters B2 and B3, positive and negative switches S21 and S32, MOS transistors Q1, Q2 and inductor L. The current flows along the red arrow. At this time, the electrical energy of the high-energy battery cluster B2 and B3 is converted into the magnetic energy of the inductor L.

2.
Target battery B5 charging: The charging loop consists of target battery B5, positive and negative switches S51, S52, MOS transistors Q3, Q4 and inductor L. The current flows along the blue arrow. At this time, the magnetic energy stored in the inductor L is converted into electric energy in the target battery B5.
Mode 3: Assuming that high-energy battery clusters B2 and B3 balance with lowenergy battery clusters B6 and B7. As shown in Figure 4, this process is divided into the following two steps.

1.
Battery clusters B2, B3 discharging: The discharging loop is composed of battery clusters B2 and B3, positive and negative switches S21, S32, MOS transistors Q1, Q2 and inductor L. The current flows along the red arrow. At this time, the electrical energy of the high-energy battery cluster B2 and B3 is converted into the magnetic energy of the inductor L.

2.
Battery clusters B6, B7 charging: The charging loop consists of battery clusters B6 and B7, positive and negative switches S61, S72, MOS transistors Q3, Q4 and inductance L. The current flows along the blue arrow. At this time, the magnetic energy stored in the inductor L is converted into electrical energy in the low-energy battery cluster B6 and B7.

Simulation Model Verification
The accuracy of the established Matlab/Simulink simulation model is verified by the equivalent circuit model of the proposed balancing topology. A simplified equivalent model of the main circuit is shown in Figure 5. Figure 6 shows its switch transistor gate drive signal. For ease of analysis, we assume that the switch is ideal.  Assuming that batteries Cell1 and Cell2 are the source and target batteries, respectively. The equalization process is divided into the following two steps.
Step 1 [0~DT]: At this time, the switch transistor Q1 is turned on, and Q2 is turned off. Source battery Cell1 charges inductor L. The red arrow shows the current path. The expression of the discharging current ( ) of the source battery Cell1, that is, the charging current ( ) of the inductor L is given by

Simulation Model Verification
The accuracy of the established Matlab/Simulink simulation model is verified by the equivalent circuit model of the proposed balancing topology. A simplified equivalent model of the main circuit is shown in Figure 5. Figure 6 shows its switch transistor gate drive signal. For ease of analysis, we assume that the switch is ideal.

Simulation Model Verification
The accuracy of the established Matlab/Simulink simulation model is verifie equivalent circuit model of the proposed balancing topology. A simplified eq model of the main circuit is shown in Figure 5. Figure 6 shows its switch transi drive signal. For ease of analysis, we assume that the switch is ideal.  Assuming that batteries Cell1 and Cell2 are the source and target batteries tively. The equalization process is divided into the following two steps.
Step 1 [0~DT]: At this time, the switch transistor Q1 is turned on, and Q2 i off. Source battery Cell1 charges inductor L. The red arrow shows the current p expression of the discharging current ( ) of the source battery Cell1, that is, th ing current ( ) of the inductor L is given by Step 2 [DT~T]: At this time, the switch transistor Q2 is turned on, and Q1 i

Simulation Model Verification
The accuracy of the established Matlab/Simulink simulation model is verified equivalent circuit model of the proposed balancing topology. A simplified equ model of the main circuit is shown in Figure 5. Figure 6 shows its switch transist drive signal. For ease of analysis, we assume that the switch is ideal.  Assuming that batteries Cell1 and Cell2 are the source and target batteries, tively. The equalization process is divided into the following two steps.
Step 1 [0~DT]: At this time, the switch transistor Q1 is turned on, and Q2 is off. Source battery Cell1 charges inductor L. The red arrow shows the current pa expression of the discharging current ( ) of the source battery Cell1, that is, the ing current ( ) of the inductor L is given by Step 2 [DT~T]: At this time, the switch transistor Q2 is turned on, and Q1 is off. Inductor L charges the target battery Cell2. The blue arrow shows the curren Assuming that batteries Cell1 and Cell2 are the source and target batteries, respectively. The equalization process is divided into the following two steps.
Step 1 [0~DT]: At this time, the switch transistor Q1 is turned on, and Q2 is turned off. Source battery Cell1 charges inductor L. The red arrow shows the current path. The expression of the discharging current i 1 (t) of the source battery Cell1, that is, the charging current i L (t) of the inductor L is given by Step 2 [DT~T]: At this time, the switch transistor Q2 is turned on, and Q1 is turned off. Inductor L charges the target battery Cell2. The blue arrow shows the current path. The expression of the charging current i 2 (t) of the target battery Cell2, that is, the discharging current i L (t) of the inductor L is given by where, E 1 is the voltage of the source battery Cell1 (V), E 2 is the voltage of the target battery Cell2 (V), R Σ1 is the total resistance of the discharging loop (Ω), R Σ2 is the total resistance of the charging loop (Ω) (both including the on-resistance of the MOSFET, the equivalent resistance of the inductance, the equivalent resistance of the wire and the internal resistance of the battery), L is the inductance value of the inductor L (H), D is the duty cycle of the PWM control signal and T is the switching period of the switching transistor (s). T = 1 f , where, f is switching frequency (Hz). Since the circuit works in continuous current mode, in order to avoid magnetic saturation of the inductor, the charging current i 2 (t), that is, the inductor's current i L (t) must be zero at the end of each cycle. This can be expressed by the following equation.
This can be simplified to (4).
Assuming that the source battery voltage E 1 = 3.45 V, the target battery voltage E 2 = 3.2 V, the inductor L = 100uH, the switching frequency f = 10 KHz and the discharging and charging loop resistance R Σ1 = R Σ2 = 0.2 Ω. Substituting the above parameters into Equation (4) gives a duty cycle of D = 0.506.
The calculated value of inductor current can be obtained by substituting the above parameter data into Equations (1) and (2). The simulated value of inductor current can be obtained by substituting into the Simulink simulation model of the equivalent circuit. The results are represented in Figure 7.
where, is the voltage of the source battery Cell1 (V), is the voltag battery Cell2 (V), is the total resistance of the discharging loop (Ω), resistance of the charging loop (Ω) (both including the on-resistance of the equivalent resistance of the inductance, the equivalent resistance of the w ternal resistance of the battery), is the inductance value of the inductor L duty cycle of the PWM control signal and is the switching period of the s sistor (s).
Since the circuit works in continuous current mode, in order to avoid m ration of the inductor, the charging current ( ), that is, the inductor's curre be zero at the end of each cycle. This can be expressed by the following equ This can be simplified to (4). (1 ) Assuming that the source battery voltage  It can be seen from Figure 7 that the inductive current changes approxim during battery equalization. The calculation results show that the maximum rent is generated at = 5.06 × 10 , which is 1.6603 A. The simulation show current at = 5.1 × 10 , which is 1.7143 A. The relative error between them It can be seen from Figure 7 that the inductive current changes approximately linearly during battery equalization. The calculation results show that the maximum inductor current is generated at t = 5.06 × 10 −5 , which is 1.6603 A. The simulation shows a maximum current at t = 5.1 × 10 −5 , which is 1.7143 A. The relative error between them is 3.25%, but this is due to the simulation time step setting. In practice, the maximum inductive current should occur at t = DT and be less than the simulated maximum value. After that, the inductor current decreases and decreases to zero at the end of the cycle. The results show the simulation model can truly express the equalization process.

Analysis of Factors Affecting the Equalization Effect
Equalization speed and efficiency are two basic indexes to judge the quality of an equalization system. When a certain amount of charge needs to be transferred to the target cell, the larger the average equalization current, the faster the equalization speed, the lower the heat loss caused by the loop resistance and the higher the equalization efficiency. From Equations (1) and (2), we can find that the average values of the source battery discharging current i 1 (t) and the target battery charging current i 2 (t) over the entire period T and in the respective discharging cycle DT and charging cycle (1 − D)T are given by where, i 1 and i 2 are the average values of the discharging current i 1 (t) and the charging current i 2 (t) over the entire period T, respectively. i 1 and i 2 are the average values of the discharging current i 1 (t) and the charging current i 2 (t) in respective discharging cycle DT and charging cycle (1 − D)T, respectively. E 1 and E 2 are the voltages of the source battery Cell1 and the target battery Cell2. R Σ1 and R Σ2 are the total resistance of the discharging and charging circuits. L is the inductance value, D is the duty cycle and T is the switching period. Heat loss rate due to equivalent resistances R Σ1 and R Σ2 is given by where, W 1 is the energy released by the source battery in one switching cycle. W R Σ1 and W R Σ2 are the heat dissipation energy caused by the equivalent resistors R Σ1 and R Σ2 in the discharging loop and the charging loop, respectively. η re represents the percentage of heat loss energy caused by all equivalent resistors in the energy released by the source battery. From Equations (5)-(12), we can know that the average equalization current and heat loss rate are affected by the source battery and target battery voltages E 1 and E 2 , duty cycle D, switching frequency f and charging and discharging loop resistors R Σ1 and R Σ2 . According to practical experience, the larger the equivalent resistance, the lower the time efficiency and energy efficiency. In addition, since the circuit operates in continuous current mode, the duty cycle D is determined by Equation (4). Therefore, the influence of the switching frequency f and the battery voltages E 1 and E 2 on the equalization effect is analyzed below. Assuming that the source battery voltage E 1 = 3.45 V, the target battery voltage E 2 = 3.2 V, the inductor L = 100 uH and the discharging and charging loop resistance R Σ1 = R Σ2 = 0.2 Ω. The relationship between the duty ratio D and the switching frequency f can be obtained by Equation (4), as shown in Figure 8, and the values are shown in Table 1 Assuming that the source battery voltage = 3.45 V, the target b = 3.2 V, the inductor = 100uH and the discharging and charging lo = = 0.2 Ω. The relationship between the duty ratio and the quency can be obtained by Equation (4), as shown in Figure 8, and the val in Table 1.  It can be found from Figure 8 and Table 1 that the duty cycle decreas crease of switching frequency. In addition, as the switching frequency incr of change of the duty cycle gradually decreases to zero. That is, the duty c tends to a certain limit value.

Influence of Switching Frequency on Average Equalizing Current
Assuming that the source battery voltage = 3.45V, the target battery 3.2V, the inductor = 100uH and the discharging and charging loop res = 0.2Ω. By substituting the switching frequency and its correspond obtained in Section 3.1.1 into Equations (5) and (6) can obtain the variatio age equalization current with the switching frequency, as shown in Figure   Figure 8. Relationship between duty cycle and switching frequency. It can be found from Figure 8 and Table 1 that the duty cycle decreases with the increase of switching frequency. In addition, as the switching frequency increases, the rate of change of the duty cycle gradually decreases to zero. That is, the duty cycle gradually tends to a certain limit value.

Influence of Switching Frequency on Average Equalizing Current
Assuming that the source battery voltage E 1 = 3.45 V, the target battery voltage E 2 = 3.2 V, the inductor L = 100 uH and the discharging and charging loop resistance R Σ1 = R Σ2 = 0.2 Ω. By substituting the switching frequency f and its corresponding duty cycle D obtained in Section 3.1.1 into Equations (5) and (6) can obtain the variation of the average equalization current with the switching frequency, as shown in Figure 9.
Assuming that the source battery voltage = 3.45V, the target battery 3.2V, the inductor = 100uH and the discharging and charging loop res = 0.2Ω. By substituting the switching frequency and its correspond obtained in Section 3.1.1 into Equations (5) and (6) can obtain the variatio age equalization current with the switching frequency, as shown in Figure   Figure 9. Variation of average equalizing current with switching frequency.  Figure 9 shows the average values i 1 and i 2 of the discharging current i 1 (t) and the charging current i 2 (t) in the whole period T decreased with the increase of the switching frequency f . For comparison at the same scale, the equation for the amount of charge q 2 transferred to the target battery Cell2 per unit time can be derived as Equation (13). The amount of charge transferred to the target cell Cell2 per unit time at different switching frequencies is obtained by substituting the corresponding parameter data, as shown in Figure 10.
nergies 2022, 15, x FOR PEER REVIEW Figure 9 shows the average values and of the discharging cur the charging current ( ) in the whole period decreased with the in switching frequency . For comparison at the same scale, the equation for charge transferred to the target battery Cell2 per unit time can be derive (13). The amount of charge transferred to the target cell Cell2 per unit tim switching frequencies is obtained by substituting the corresponding para shown in Figure 10. (1 ) Figure 10 shows as the switching frequency increases, the amount of c time transferred to the target cell Cell2 decreases and tends to a certain lim from the point of view of increasing the equalization current and improvin ciency, the switching frequency should be chosen as a smaller value.  Figure 10 shows as the switching frequency increases, the amount of charge per unit time transferred to the target cell Cell2 decreases and tends to a certain limit. Therefore, from the point of view of increasing the equalization current and improving the time efficiency, the switching frequency should be chosen as a smaller value.

Influence of Switching Frequency on Heat Loss Rate
Assuming that the source battery voltage E 1 = 3.45 V, the target battery voltage E 2 = 3.2 V, the inductor L = 100 uH and the discharging and charging loop resistance R Σ1 = R Σ2 = 0.2 Ω. By substituting the switching frequency f and its corresponding duty cycle D obtained in Section 3.1.1 into Equations (9)- (12) can obtain the variation of the heat loss rate with the switching frequency, as shown in Figure 11. from the point of view of increasing the equalization current and ciency, the switching frequency should be chosen as a smaller va 3.1.3. Influence of Switching Frequency on Heat Loss Rate Assuming that the source battery voltage = 3.45V, the tar 3.2V, the inductor = 100uH and the discharging and chargin = 0.2Ω. By substituting the switching frequency and its co obtained in Section 3.1.1 into Equations (9)- (12) can obtain the rate with the switching frequency, as shown in Figure 11.  It can be seen from Figure 11 that the higher the switching frequency, the lower the heat loss rate, and it tends to a constant value. Therefore, selecting a medium frequency switch is beneficial to improving the equalization circuit's energy efficiency while reducing the hardware cost. Given that later studies show that equalization time can be significantly reduced by MC2MC equalization, the switching frequency f = 10 KHz is chosen in this paper to reduce hardware costs and improve the energy efficiency of the equalization circuit.

Influence of Voltage Difference on Duty Cycle
The voltage difference ∆E refers to the difference between the source battery voltage E 1 and the target battery voltage E 2 during the equalization process. It is given by Compared with realizing the AC2AC equalization path, the MC2MC equalization method will produce a large voltage difference between the source and target batteries. This voltage difference will change in the equalization process. Assuming that the target battery voltage E 2 = 3.2 V, the inductor L = 100 uH, the switching frequency f = 10 KHz and the discharging and charging loop resistance R Σ1 = R Σ2 = 0.2 Ω. The relationship between the duty cycle D and the voltage difference ∆E can be obtained by Equation (4), as shown in Figure 12, and the partial values are shown in Table 2.
This voltage difference will change in the equalization process. Assuming battery voltage = 3.2 V, the inductor = 100 uH, the switching frequenc and the discharging and charging loop resistance = = 0.2 Ω. The re tween the duty cycle and the voltage difference Δ can be obtained by as shown in Figure 12, and the partial values are shown in Table 2.  It can be found from Figure 12 and Table 2 that the duty cycle decreases difference increases. When one battery is balancing to one battery (Δ = 0), is 0.525. When two batteries are balancing to one battery (Δ = 3.2), the duty When three batteries are balancing to one battery (Δ = 6.4), the duty cycle i fore, in the MC2MC equalization system, the change in the duty cycle change in the voltage difference cannot be ignored; it is necessary to introd duty ratio.

Influence of Voltage Difference on Average Balancing Current
Assuming that the target battery voltage = 3.2V, the inductor switching frequency = 10KHz and the discharging and charging loop res  It can be found from Figure 12 and Table 2 that the duty cycle decreases as the voltage difference increases. When one battery is balancing to one battery (∆E = 0), the duty cycle is 0.525. When two batteries are balancing to one battery (∆E = 3.2), the duty cycle is 0.356. When three batteries are balancing to one battery (∆E = 6.4), the duty cycle is 0.269. Therefore, in the MC2MC equalization system, the change in the duty cycle caused by the change in the voltage difference cannot be ignored; it is necessary to introduce a variable duty ratio.

Influence of Voltage Difference on Average Balancing Current
Assuming that the target battery voltage E 2 = 3.2 V, the inductor L = 100 uH, the switching frequency f = 10 KHz and the discharging and charging loop resistance R Σ1 = R Σ2 = 0.2 Ω. By substituting the voltage difference ∆E and its corresponding duty cycle D obtained in Section 3.2.1 into Equations (7) and (8) can obtain the variation of the average balancing current with the voltage difference, as shown in Figure 13. Figure 13 shows the average values i 1 and i 2 of discharging current i 1 (t) and charging current i 2 (t) in their respective discharging and charging cycles increase with the increase of voltage difference. Based on the transfer of the same amount of charge to the target battery Cell2, the rate of improvement of the equalization time by increasing the voltage difference can be calculated. It is given by where, i 2 (1), D 1 and X 1 respectively represent the average value, duty cycle and required number of cycles of the charging current i 2 (t) under a certain voltage difference. i 2 (2), D 2 and X 2 respectively represent the average value, duty cycle and required number of cycles of the charging current i 2 (t) under another voltage difference. = 0.2Ω. By substituting the voltage difference Δ and its corresponding du obtained in Section 3.2.1 into Equations (7) and (8) can obtain the variation of t age balancing current with the voltage difference, as shown in Figure 13.  (1) (1 ) where, ′ (1), and respectively represent the average value, duty cycle quired number of cycles of the charging current ( ) under a certain voltage dif ′ (2), and respectively represent the average value, duty cycle and r number of cycles of the charging current ( ) under another voltage difference.
(1) The time improvement rate of the two-to-one battery equalization than the one battery equalization.
Substituting ′ (1) = 0.7849 , = 0.525 as well as ′ (2) = 1.0759 , into Equation (15) gives the following result. 2 1 =0.538 X X ⋅ It shows that under the condition of transferring the same amount of charg target battery, using two batteries to balance one battery, the time efficiency is im by 46.2% compared to using only one battery to balance one battery.
(2) The time improvement rate of the three-to-one battery equalization than the one battery equalization.
Substituting ′ (1)=1.0759, = 0.356 as well as ′ (2) = 1.2244, = 0.2 Equation (15) gives the following result. 2 1 =0.7742 X X ⋅ It shows that under the condition of transferring the same amount of charg target battery, using three batteries to balance one battery, the time efficiency is im by 22.58% compared to using two batteries to balance one battery.
(3) The time improvement rate of the three-to-one battery equalization than the one battery equalization. (1) The time improvement rate of the two-to-one battery equalization than the one-to-one battery equalization.
It shows that under the condition of transferring the same amount of charge to the target battery, using two batteries to balance one battery, the time efficiency is improved by 46.2% compared to using only one battery to balance one battery.
(2) The time improvement rate of the three-to-one battery equalization than the two-toone battery equalization.
It shows that under the condition of transferring the same amount of charge to the target battery, using three batteries to balance one battery, the time efficiency is improved by 22.58% compared to using two batteries to balance one battery.
(3) The time improvement rate of the three-to-one battery equalization than the one-toone battery equalization.
It shows that under the condition of transferring the same amount of charge to the target battery, using three batteries to balance one battery, the time efficiency is improved by 58.35% compared to using only one battery to balance one battery.
To further illustrate the tendency of the voltage difference to enhance the equalization time, the amount of charge Q 2 received by the target cell in one cycle is derived for different voltage differences. It is given by Equation (16). By substituting the parameter data into the above equation, the relationship between the amount of charge Q 2 received by the target battery, and the voltage difference in a cycle can be obtained in Figure 14. 16) It shows that under the condition of transferring the same amou target battery, using three batteries to balance one battery, the time eff by 58.35% compared to using only one battery to balance one battery.
To further illustrate the tendency of the voltage difference to en tion time, the amount of charge received by the target cell in one different voltage differences. It is given by Equation (16). By substitu data into the above equation, the relationship between the amount of by the target battery, and the voltage difference in a cycle can be obta (1 ) (1 ) Figure 14 shows that as the voltage difference increases, the am ceived by the target battery in one period also increases accordin MC2MC equalization method can effectively reduce the equalization the time efficiency.

Influence of Voltage Difference on Heat Loss Rate
Assuming that the target battery voltage = 3.2 V, the induc switching frequency = 10 KHz and the discharging and charging lo = 0.2 Ω. Substituting the voltage difference and its corresponding in Section 3.2.1 into Equations (9) to (12) can obtain the variation of th the voltage difference, as shown in Figure 15.  Figure 14 shows that as the voltage difference increases, the amount of charge received by the target battery in one period also increases accordingly. Therefore, the MC2MC equalization method can effectively reduce the equalization time and improve the time efficiency.

Influence of Voltage Difference on Heat Loss Rate
Assuming that the target battery voltage E 2 = 3.2 V, the inductor L = 100 uH, the switching frequency f = 10 KHz and the discharging and charging loop resistance R Σ1 = R Σ2 = 0.2 Ω. Substituting the voltage difference and its corresponding duty cycle obtained in Section 3.2.1 into Equations (9) to (12) can obtain the variation of the heat loss rate with the voltage difference, as shown in Figure 15. It can be seen from Figure 15 that the heat loss rate and the vo positively correlated. The heat loss rate of one-to-one battery equ The heat loss rate of two-to-one battery equalization is 12.54% It can be seen from Figure 15 that the heat loss rate and the voltage difference are positively correlated. The heat loss rate η re1 of one-to-one battery equalization is 12.47%. The heat loss rate η re2 of two-to-one battery equalization is 12.54%. The heat loss rate η re3 of three-to-one battery equalization is 12.56%. The increase in heat loss rate is small, so the equalization method of MC2MC can effectively shorten the equalization time. The increase in heat loss is small, which has the dual advantages of high time efficiency and high energy efficiency.

Balancing Control Strategy
In this paper, eight series batteries are selected for equalization simulation verification. The simulation model of this topology is established in MATLAB/Simulink2021a, including battery module, switching array, Buck-Boost equalization unit, switching control module and duty cycle control module based on voltage difference. The duty cycle control module determines the duty cycle in the MC2MC equalization process. SOC is selected as the equalization variable. Firstly, the initial charge state of each cell is used to determine whether a multi-cell cluster can be formed in the series battery pack. If the battery cluster cannot be formed, the AC2AC balancing strategy is performed, as shown in Figure 16a; otherwise, the MC2MC balancing strategy is performed, as shown in Figure 16b. In order to prevent the starting condition of the equalization circuit from being quickly satisfied again after the equalization of the battery pack is completed, the equalization start threshold Th_1 and the equalization termination threshold Th_2 are respectively set. In this paper, the start threshold Th_1 is 3, and the termination threshold Th_2 is 1.

Multi-Cell-to-Any-Cell Simulation Verification
Assuming that the SOC values of the eight series batteries are 42. 3, 43.8, 42, 46.5, 44.5, 45.8, 46.2, and 42.5, respectively, and the energy state distribution diagram of the battery pack is shown in Figure 17. If the SOC value of a battery meets SOC j ≥ ∑ 8 i=1 SOC i 8 + Th 2 2 , this battery is a high-energy battery, which is represented by yellow. If the SOC value of a battery meets SOC j ≤ ∑ 8 i=1 SOC i 8 − Th 2 2 , this battery is a low-energy battery, which is represented by red. The remaining batteries are normal batteries, indicated in green.
From Figure 17, it can be seen that two adjacent batteries, Cell6 and Cell7, can form a high-energy cell cluster. Figure 18a shows the simulation results of the AC2AC balancing strategy, which shows that the highest SOC battery is balancing with the lowest SOC battery. The equalization time is 662s, and the SOC of each battery after equalization is 43. 49, 43.79, 43.49, 44.49, 44.49, 44.49, 44.49 and 43.49, with an average of 44.0275 and the variance of 0.2223. Figure 18b shows the simulation results of the MC2AC equalization strategy. It can be seen that, at time 0 ∼ t1, the high-energy battery cluster formed by Cell6 and Cell7 balances the lowest SOC battery. After time t1, since Cell6 has finished equalizing, Cell6 and Cell7 can no longer form the high-energy battery cluster, and the subsequent equalization is transferred to the AC2AC equalization process, where t1 is 180.228 s. The equalization time is 557.126 s, and the SOC of each battery after equalization is 43.54, 43. 8, 43.54, 44.54, 44.49, 44.54, 44.54 and 43.54, with an average of 44.06625 and the variance of 0.2193. The comparison of the two control strategies is shown in Table 3.
As can be seen from the above, compared with the AC2AC equalization strategy, MC2AC equalization has the following advantages.   From Figure 17, it can be seen that two adjacent batteries, Cell6 and Cell7, c a high-energy cell cluster. Figure 18a shows the simulation results of the AC2AC ing strategy, which shows that the highest SOC battery is balancing with the low battery.  Figure 18b shows the simulation results of the MC2AC equa strategy. It can be seen that, at time 0~ 1, the high-energy battery cluster for Cell6 and Cell7 balances the lowest SOC battery. After time 1, since Cell6 has equalizing, Cell6 and Cell7 can no longer form the high-energy battery cluster, subsequent equalization is transferred to the AC2AC equalization process, whe    From Figure 17, it can be seen that two adjacent batteries, Cell6 and Cell7, can form a high-energy cell cluster. Figure 18a shows the simulation results of the AC2AC balancing strategy, which shows that the highest SOC battery is balancing with the lowest SOC battery. The equalization time is 662s, and the SOC of each battery after equalization is 43. 49 Figure 19. Among them, yellow represents the high-energy battery, green represents the normal battery, and red represents the low-energy battery.

Multi-Cell-to-Multi-Cell Simulation Verification
Assuming that the SOC values of the eight series batteries are 43.8, 42.3 44.5, 45.8, 46.2 and 42, respectively, the energy state distribution diagram of pack is shown in Figure 19. Among them, yellow represents the high-ener green represents the normal battery, and red represents the low-energy batter From Figure 19, it can be seen that two adjacent batteries, Cell6 and Cell a high-energy cell cluster, and Cell2 and Cell3 can form a low-energy cell clu 20a shows the simulation results of the AC2AC balancing strategy, which sho highest SOC battery is balancing with the lowest SOC battery. The equalizat 662 s, and the SOC of each battery after equalization is 43. 8  From Figure 19, it can be seen that two adjacent batteries, Cell6 and Cell7, can form a high-energy cell cluster, and Cell2 and Cell3 can form a low-energy cell cluster. Figure 20a shows the simulation results of the AC2AC balancing strategy, which shows that the highest SOC battery is balancing with the lowest SOC battery. The equalization time is 662 s, and the SOC of each battery after equalization is 43. 8 Figure 20b shows the simulation results of the MC2MC equalization strategy. It can be seen that, at time 0 ∼ t1, the high-energy battery cluster formed by Cell6 and Cell7 balances the low-energy battery cluster formed by Cell2 and Cell3. At time t1 ∼ t2, since Cell6 has finished equalizing, Cell6 and Cell7 can no longer form the high-energy battery cluster. The subsequent equalization is transferred to the AC2MC balancing process; the highest SOC battery balances the low-energy battery cluster formed by Cell2 and Cell3. After time t2, since Cell3 has finished equalizing, Cell2 and Cell3 can no longer form the low-energy battery cluster, and the subsequent equalization is transferred to the AC2AC equalization process, where the equalization turning points t1 and t2 are 88.037 s and 157.458 s respectively. The equalization time is 467.379 s, and the SOC of each battery after equalization is 43.8, 43.54, 43.63, 44.54, 44.49, 44.54, 44.54 and 43.54, with an average of 44.0775, and the variance of 0.2084. The comparison of the two equalization strategies is shown in Table 4.    As can be seen from the above, compared with the AC2AC equalization strategy, MC2MC equalization has the following advantages. (1) It can effectively shorten the equalization time and increase the time efficiency by 29.4%. (2) The variance is reduced by 6.04%, indicating that the equalization results are less discrete and the results are more stable. (3) Due to the shortened equalization time, the energy efficiency is not reduced but slightly increased by 0.1107%.

Conclusions
This paper proposes an active balancing scheme for battery packs to achieve multicell-to-multi-cell balancing. The operation principle of the balancing topology combining the Buck-Boost circuit and switch array is analyzed in detail, and the accuracy of the simulation model is verified. In addition, to explore the influence of various factors on the equalization effect, the switching frequency and voltage difference are selected for analysis and calculation. The results show that the larger the switching frequency is, the smaller the average equalization current and heat loss rate are. Therefore, considering the time efficiency and energy efficiency, the switching frequency should be about 10 KHz. Voltage difference significantly affects the duty cycle, so it is absolutely necessary to introduce a variable duty cycle in MC2MC equalization. Using three batteries and two batteries to balance one battery results in a 58.35% and 46.2% improvement in time efficiency over using only one battery to balance one battery, respectively. Finally, eight batteries are selected to build the simulation model. Under the condition that each battery's initial state of charge is the same, but the location is different, compared with any-cell-to-any-cell (AC2AC) equalization, the time efficiency of multi-cell-to-multi-cell equalization is increased by 15.842% and 29.4%, respectively. The variance is reduced by 1.35% and 6.04%, respectively, and because the equalization time is shortened, the energy efficiency is slightly improved.