Abstract

The reliability and economy of dc-link capacitors are important concerns in multiphase drive systems. Due to the parallel connection of several converters, the dc-link capacitors are subjected to a higher RMS current, and the root mean square (RMS) current of dc-link capacitors is an important reference standard to determine its lifetime, cost, and volume. In this paper, the RMS current of dc-link capacitor is calculated by using the dual Fourier integral method and the effect of carrier interleave is studied. Meanwhile, the modulation ratio, harmonic sidebands, and switching frequency are also considered. In order to optimize the reliability and economy of the multiphase drive system, a Cotes method combined with carrier-phase shifting technology (CPST) for calculating RMS current of the dc-link capacitor is proposed. The proposed method can provide optimization guidance for the design of dc-link capacitors. Finally, the analytical and experiment results are compared with the existing methods, and the experimental results verify the effectiveness of the proposed method.

1. Introduction

Nowadays, the multiphase machine is gaining more and more attention because of the advantages of low power density, flexible control, and high fault tolerance [1]. In multiphase machine drives, the multi-inverters connected on a common capacitor structure. Inevitably larger current ripples are generated in the dc-link capacitors. Naturally, the RMS current on the capacitors will significantly increase and have an effect on its thermal management, size design, and lifetime [24]. Therefore, it is inevitable to design a reliable and economical capacitor in a multiphase drive system.

On the one hand, the reduction of the RMS value is a very important aspect. Some researchers have focused on improving system stability and reducing the size of dc-link capacitors by reducing the current RMS value [59]. Taking the six-phase machine drives as an example, carrier interleaving between upper and lower inverters will reduce the dc-link current ripple [1014]. If the interleaving angle reaches 90 degrees, the total current amplitude is reduced by nearly 50% [14]. Obviously, the RMS current of capacitors has sharply changed. Therefore, CPST can reduce the size and temperature rise in the capacitor design benchmark.

On the other hand, the accuracy of RMS current calculation needs careful attention, which influences the capacitance selection and lifetime monitoring [15]. Actually, the inverter uses different control strategies and the dc-link capacitor will bear different harmonic current [16, 17]. The calculation method of the RMS value of capacitor current has been reported in the literature, and it could be classified as two types: time domain analysis and harmonic spectrum analysis.

In time domain, three combined parameters (modulation depth, current amplitude, and phase angle) decided the RMS current together. A simple analytical expression is used to calculate the RMS current, while it possibly exists about 5% deviation compared to its actual value [18].

In frequency domain, the capacitor RMS current is calculated according to the harmonic frequency analysis, in which the harmonic components are obtained with a double Fourier integration [17, 19, 20]. This method is clearly simple and easy to accomplish [19]. Besides, the effects on the RMS current under different modulation methods are analyzed in detail, and its analysis results were verified [17]. A polynomial interpolation integral method was used to calculate the RMS current, which could eliminate the external integral term and thus simplify the analytical modeling [20, 21]. The calculation methods have low accuracy and complicated processes in the existing literature, and there is still no paper to calculate the RMS current considering the CPST effect.

Therefore, a calculation method of RMS current is proposed for multiphase machine drives with CPST in this paper. This is sufficient to reflect the contribution of reliability and economy to the multiphase drive system [22]. Firstly, the effect of converter modulation on dc-link capacitor current is analyzed. Secondly, the double Fourier integration is obtained through the frequency spectrum analysis of capacitor current. In order to reduce the integration calculation process and improve the accuracy, an improved polynomial interpolation quadrature method (Cotes) is proposed. Finally, the RT-lab platform is used to verify the effectiveness of the proposed method under various operating conditions.

The rest of this article is organized as follows. Section 2 introduces the system and the impact of converter’s modulation on capacitor current. In Section 3, an improved polynomial interpolation method using double Fourier analysis is proposed. Section 4 discusses the effect of CPST on multiphase machine drive. Section 5 gives the simulation results of the proposed method. Section 6 concludes this article.

2. Basic Theoretical Analysis

Figure 1 shows the topology of the six-phase PMSM drives with two inverters, which are connected to a common capacitor in parallel way. The uncontrolled rectifier is used to provide the power supply. The following theoretical analysis in this paper is based on this topology, and the six-phase machine is assumed to be working under a steady state. In addition, the inverter adopts the SVPWM control method.

2.1. Rectifier Output Current Analysis

The switching rule of diode rectifier bridge is natural commutation. In a switch cycle, the opening sequence of the rectifier diode is D1⟶D6⟶D2⟶D4⟶D3⟶D5. Therefore, the output current is composed of the dc component and the sum of harmonic component as follows:

Only the ac harmonic components can flow through the dc-link capacitor. Therefore, the capacitor RMS current is mainly composed of the harmonic component of rectifier output.

2.2. Inverter Input Current Analysis

The dc-link current consists of the input current of the two inverters, including the ac component and dc component. It can be represented by

The dc component is supplied directly by the rectifier output, and the ac component is affected by the switching strategy of the inverter. The ac component consists of the sum of the six-phase switching currents. The switching current of each phase can be represented by the multiplication of the ac phase current and switching function T(t). Therefore, the ac component can be written as

Due to the switching behavior of the converters, the AC current can be represented by just fundamental frequency component aswhere and represent the individual current amplitudes in the respective phases. The fundamental frequency is defined as ωf1 = 2πff1 or ωf2 = 2πff2, where ff1 and ff2 are the current frequency of the upper and lower inverter ( = 0, 1, 2). φ1 and φ2 are the phase angles of the upper and lower inverter. τ is the phase shift angle. If the pulse width of the switch signal T1(t) is normalized to the switching period, equation (5) shows the time function in first and second sectors.

Equation (5) is based on the time definitions of the switching states for SVM (space vector modulation) as explained in [16, 21]. The variable M is the modulation index according towhere is the reference voltage vector. The switch signal of switch T7(t) is normalized to the switching period as

The angles φ of the normalized and the fundamental wave shift β are defined in Figure 2. The angle between the current time function of the U phase and fundamental wave time functions α1(t) in Figure 2 is φ1. In addition, the fundamental wave time functions α1(t) and α2(t) of the upper and lower inverter control strategies are shifted by an angle β.

3. Analytic Spectrum Calculation

Based on the topology in Figure 1, the RMS value of dc-link capacitor current can be affected by the ac component of both rectifier output and inverter input [15]. can be calculated as follows:

It can be seen that the RMS value of the inverter input current affects the value. Obviously, the value is the sum of multi-inverters. For six-phase machine drives, can be represented as

3.1. Description of the Current Spectrum as an Integral Function

According to Parseval’s theorem, the energy of a time signal can be equal to the energy of its spectrum. Therefore, the time signal can be expressed by spectrum amplitude . As explained in [20], the upper-inverter spectrum results from the sum of the spectra of the three phase currents according to

Furthermore, the geometric wall model method can be used to describe the switching current spectrum of the dual inverter. The method of the geometric wall model is used in [12]. Different from the analysis of switching current in time domain, the spectrum of switching current can be expressed by a double Fourier integral. Using this method, the spectrum of the inverter switch current is listed as follows:where , , and can be represented by (4), yu, yd, xu, and xd are the integration limits of the inner and outer integral, m is the carrier index variable, and n is the baseband index variable. The upper and lower limits of integration are shown in Table 1. y = , x = , and with m, n.

As explained in Figure 2, the normalized switch functions α1 and α7 are shifted by the angle β. Similar to the representation of , the spectrum of the lower inverter can be derived with (11) and the results are as follows:where , , and can be represented by (4). y = , x = , and , with , h ∈ N.

3.2. Simplified Analytical Modeling for the Calculation of the Spectrum

The expressions given in equations (11)-(12) are elementary integrals which involves many unsolved integrals. In order to obtain a simplified calculation expression, a Cotes method is carried out for the expression (11)-(12). The Cotes method is based on the idea of interpolating the original function [22]. For example, the Newton 3/8 formulas are numerical quadrature formulas for the approximate determination of an antiderivative of a function. Within the considered area, the locations of the required sampling points are in equidistant intervals. For the integral of the resulting function, an exact solution can be generated, and equations (11)-(12) are the double Fourier summation integral, which has a huge and tedious amount of calculation. Using the proposed Cotes interpolation quadrature method can eliminate a lot of the integration calculation process, and it also has a high calculation accuracy.

However, the number of sampling points in Newton rule is relatively small, resulting in a low accuracy of approximate results. Therefore, in order to obtain an approximate solution with higher accuracy, the Cotes method is proposed in this paper to simplify the calculation formula (11)-(12). The results based on the Cotes rule are presented below.

Firstly, the inner integral term in (11)-(12) is solved, as follows:where (U, V, W) corresponds to τ = (1, 2, 3). Then, the results of applying the Cotes rule are as follows:

Here, C(y) represents the function approximated to the original function. The coefficient represents the weights, which can be found in Table 2. The coefficient describes the location of the steps and is defined as

Hence, k ∈ {1, 2, 3, 4, 5}. The variable represents the normalized location of the steps, which can also be found in Table 2. It can be seen that this solution does not contain any integral or differential expressions. Therefore, the Cotes method has more simplified calculation accuracy and the process.

Similarly, can be described. The sum of to then defines the total spectrum . An approximate description of the total spectrum takes into account the first six harmonics and their sidebands. In the case of odd harmonics, the sidebands are around ±3ff, and for even harmonics, they occur around ±6ff. Equation (16) can be obtained from this rule.

4. Influence of CPST on the RMS Value of Capacitive

Compared to the analysis method in the time domain, it is easier to explain the principle of interleaving in the frequency domain. According to equations (11)-(12), the fundamental wave displacement β is not considered, and the difference between and is determined by G. As shown in the following equation,

Obviously, the is maximum when the switching cycle clock signal without interleaving. is going to decrease by selecting an appropriate displacement angle G within [0-]. According to (8) and (17), can be reduced by CPST. The switching cycle shift is the angle between the control signals from the upper inverter and lower inverter [23], as shown in Figure 3.

5. Verify the Validity of the Proposed Method

In order to verify the effectiveness of the proposed method, we used the RT-lab experimental platform for verification. All the experiment data are obtained by using the RT-LAB experiment platform, and the RT-LAB experiment platform is shown in Figure 4. The RT-LAB can implement the hardware in the loop simulation (HILS) for the six-phase machine drive system. The model of the DSP controller is TMS320F2812, which runs the algorithm, and RT-LAB (OP5600) is used to construct other parts of the system such as the six-phase machine and dual inverter.

The system parameters are shown in Table 3. The switching frequency and control frequency are both 4 kHz. The phase currents are considered as ideal cosine functions. The amplitude of the input three-phase voltage is 326 V, and the machine operation frequency is 50 Hz.

Figure 5 shows the dc-link capacitor current before and after the CPST method. To make it a fair comparison, the parameters of the drive system are strictly the same. It can be seen from Figure 5 that the capacitor current ripple is significantly reduced using the CPST method. Besides, the amplitude of current ripple decreases obviously. Therefore, Figure 5 shows that CPST has an inhibitory effect on the RMS value of the dc-link capacitor current.

5.1. Verification of the Effect of CPST by the Experiment

As described in previous sections, the RMS value of dc-link capacitor current could be reduced by the CPST. As can be seen from equations (11)-(12), is highly correlated with variables M, G, and β. Figure 6 shows the spectrum analysis of under different conditions of the modulation index M and the displacement angle G = 0 and G = 90°, and the variable β is considered to be the constant zero.

In order to verify the effect of CPST, the experiments for dual-three phase machine are performed, and the system models were loaded into RT-lab, respectively, to observe the effect of dc-link capacitor current. Figure 6 shows that the results of the capacitor current FFT analysis using CPST at three different modulation ratios 1, 0.8, and 0.5. The results show that the proposed CPST method could reduce the 2 fs harmonics effectively at various modulation ratios, and the diminution of 2 fs harmonics current is the main reason of dc-link capacitor RMS current reduction.

5.2. Verification of the Approximate Effect of the Cotes Rule by Calculation

In order to verify the simplification of the Cotes method, a comparison experiment of the calculation speed is established, such as the Cotes method and MATLAB double Fourier calculation method. Equations (11)-(12) can be computed in two different ways. The programs of the two calculation methods are executed by the DSP controller (TMS320F2812), and the execution time is recorded by inserting breakpoints. Figure 7 shows the comparison result of calculation speed. The program execution time of the Cotes method is about 5 ms, and the MATLAB FFT method is about 8 ms. From the results, we can see that the Cotes method is about 30% faster. Therefore, the proposed Cotes method is more simplified.

In order to verify the effectiveness of the proposed method, the results of two approximate methods are compared, considering the carrier index m. Equations (11)-(12) define the original function, which are applied to the Newton 3/8 method proposed in literature [20] and the Cotes method proposed in this paper. Then, the standard integral calculation method in [17], the approximate results of the Cotes method, and Newton 3/8 method are compared.

In order to verify the accuracy of the proposed Cotes method, three different methods (Newton 3/8, FFT analysis, and the proposed method) are used to calculate the RMS value of dc-link capacitor current. Figure 8 shows the RMS value of the dc-link capacitor current as a function of M for different 1–6 times the switching frequency (after CPST). The black line represents the calculation results of the method proposed in [17], the red line represents the calculation results of the method proposed in [20], and the blue line represents the calculation results of the method proposed in this paper. The high precision of the Cotes method is proved by the graphic results. Comparing the three methods, the Cotes method proposed in this paper has a better approximation degree than the MATLAB FFT analysis method. Therefore, the higher calculation accuracy and simplified calculation process of the Cotes method are confirmed.

Figure 9 shows the calculation error of comparison results between the Newton 3/8 rule and Cotes rule calculated spectrum for the switch frequency. The calculated results are 7.021 A and 7.101 A, respectively, under two different methods, and the proposed Cotes method is more similar to the MATLAB FFT calculation results. Therefore, the error of the Cotes method is less than that of the Newton 3/8 rule method.

6. Conclusions

In this paper, a new calculation expression of the RMS value of dc-link capacitor current is given in the six-phase machine drive. Compared with existing methods, a simplified form of calculation is obtained. This method calculates the integral by using the idea of polynomial interpolation, which not only reduces the calculation process but also improves the accuracy of calculation results. Therefore, the simplified calculation method of the capacitor current RMS value is combined with the CPST method, which lays a foundation for the optimal design of dc-link capacitor cost, volume, and service life.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Key Projects of Scientific Research of the Hunan Provincial Education Department (22A0603), major special projects of the Changsha Science and Technology Plan (kq2105001), and the National Natural Science Foundation of China (51907061)