Next Article in Journal
Recognition of American Sign Language Gestures in a Virtual Reality Using Leap Motion
Next Article in Special Issue
A Stochastic Deterioration Process Based Approach for Micro Switches Remaining Useful Life Estimation
Previous Article in Journal
Wireless Channel Models for Over-the-Sea Communication: A Comparative Study
Previous Article in Special Issue
Latent Leakage Fault Identification and Diagnosis Based on Multi-Source Information Fusion Method for Key Pneumatic Units in Chinese Standard Electric Multiple Units (EMU) Braking System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fault Prediction Model of High-Power Switching Device in Urban Railway Traction Converter with Bi-Directional Fatigue Data and Weighted LSM

1
School of Electrical Engineering, Beijing Jiaotong University, Beijing 100044, China
2
Beijing Electrical Engineering Technology Research Center, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(3), 444; https://doi.org/10.3390/app9030444
Submission received: 23 December 2018 / Revised: 23 January 2019 / Accepted: 23 January 2019 / Published: 28 January 2019
(This article belongs to the Special Issue Fault Detection and Diagnosis in Mechatronics Systems)

Abstract

:
The switching device is relatively weakest in the traction converter, and this paper aims at the fault prediction of it. Firstly, the mathematical distribution is analyzed based on the results that were obtained in electro thermal simulation and a single-directional accelerated fatigue test. Then, the accelerated fatigue test with bi-directional fatigue current is proposed, the data from which reflects the accelerating effect from FWD on the device aging process. The analytical model of fatigue process is fitted with the data that were obtained in the test. In order to shorten the test time consumption, we propose a weighted least squares method (LSM) to fit the failure data. Finally, the prediction model is presented with the consideration of fatigue signature and Arrhenius temperature factor.

1. Introduction

The urban rail transit system is a reliable and environmentally friendly way of transport; it generates low pollution, has high passenger capacity and is not prone to delays. The traction converter (TC) in the EMU provides motive and electric braking power through traction motors. TC’s switching IGBTs are high-power devices that are prone to failure compared to other components [1]. Switching device fault accounts for more than 10% of the total TC failures. In most cases, fatigue causes the failure of the switching device [2]. If the fatigue state of the switching device can be predicted and manually maintained or replaced in the device [3], the train failure caused by the failure of the switching device can be reduced.
The switching device in the TC is usually a package module composed of an IGBT and an FWD. The fatigue of the switching device is mainly manifested by package failure, including two main forms: Wire liftoffs and solder layer cracking [4,5]. Wire liftoffs and solder layer cracking are caused by shear thermal stresses that are produced by temperature drift inside the device. In the intermittent operation of the device when it is turned on, heat is generated by current; when it is turned off, i.e., the current is cut off, the heat is dissipated. In such a way, the device is heated and cooled periodically. The internal temperature of the device (not only the IGBT junction temperature Tj, however also the temperature of each single layer) would increase and decrease accordingly. One switching device consists of multiple layers (as shown in Figure 1) with different coefficients of thermal expansion (CTE) of each layer.
The relationship between thermal stress and the change of CTE is shown in Equation (1).
S = f ( Δ C T E , Δ T )
S is the thermal stress, ΔCTE is the difference of CTE and ΔT is junction temperature variation. Equation (1) shows that the thermal stress is caused by the difference of CTE and ΔT.
Such temperature change is caused by internal thermal resistance (such as Rth-c_1 and Rth-c_2 in Figure 1) and heat flow, as is shown in Equation (2)
Δ T = P R t h
In Equation (2), P is the heating power. Rth is the thermal resistance.
Temperature changes causes shear thermal stress, while shear thermal stress causes liftoffs and cracks [6]. Due to the encapsulation, internal liftoffs and cracks cannot be observed directly; several sensitive parameters have been considered to represent the fatigue characteristics instead, such as the IGBT saturation voltage drop (UCES) [7], SOA [8], thermal resistance between IGBT junction and device substrate (Rth [9], threshold voltage of the IGBT gate [10], etc. These methods are based on the one-to-one relationship between fatigue signatures and fatigue levels. However, these parameters are difficult to obtain reliably in the situation of high power and high EMI (Electromagnetic Interference), such as converters.
In other cases, liftoffs and cracks can be estimated by experimental data where a certain number of thermal cycles and Temperature drift (obtained from the mission profile) are applied to the IGBT under test; we call it AFT (accelerated fatigue test). Then, the model is built based on fatigue observations from the AFT. The model that is established by the AFT is the fitting relationship between the fixed working conditions and the fatigue life of switching devices.
An AFT tests a DUT with a predetermined time-effect force and obtains data showing the fatigue characteristics of it. However, the existing fatigue model requires stable and fixed operating conditions because these models were established under the AFT results of fixed experimental conditions [4,11,12]. Therefore, the existing model cannot be fully applied to the changing conditions of the actual operating conditions, and this is the actual application of the TC. For example, in practical TC operation, the switching device is affected by the changing ambient temperature due to the movement of the EMU, while the ambient temperature remains constant in AFT [13]; in addition, the switching frequency of the switching device (i.e., The number of currents turn-ons and turn-offs in one second) and currents are also constantly changing as the speed of the EMU or the passenger load changes, however the conventional AFT only considers the constant switching frequency and the constant test current [14].
In this paper, it shows that Weibull distribution is more suitable to describe the cumulative failure rate of high-power switching device fatigue based on the data from conventional AFT. Then, a novel bidirectional aging approach is proposed to show the aging effect from FWD inside the switching device. In the process of AFT, the experiment time is often too long under low current. Therefore, a weighted LSM (least squares method) is proposed to complete the piecewise fitting of BAFT (bidirectional AFT) experimental data, which improves the accuracy of modeling. Through the aging model that is proposed in this paper, the service life of switching devices under actual working conditions can be predicted and the fault trend can be evaluated online when the load current of switching devices is known.

2. Fatigue Characteristics of High-Power IGBT Model

2.1. Analysis of the Results of Conventional Single-Directional Aging Experiments

The analytic model of the IGBT is the analytical relationship between the sensitive factor in the IGBT and the number of failures (i.e., the lifetime). The aging data is fitted correspondingly, and the expression is intuitive. In order to describe the relationship between the cumulative failure rate F and the number of experiments (i.e., the number of cycles), it is necessary to select an appropriate life distribution. If the distribution model is properly selected, the F under various working conditions could be obtained with more accuracy [15]. Common life distribution models include normal distribution, lognormal distribution, Weibull distribution and log-Weibull distribution. The choosing of distribution model is the key to switching device fault reconfiguration.
The normal distribution, also known as the Gauss distribution, is the most common and widely used distribution of the probability distribution of all random phenomena. It can be used to describe many natural phenomena and various physical properties. In practical applications, there are many experimental data that can be fitted with a normal distribution as its approximate distribution [16]. The lognormal distribution is a typical model that is derived from a normal distribution. When the logarithm of the expiration time “ln (t)” obeys the normal distribution, t obeys logarithmic normal distribution. Logarithmic transformation can reduce the larger number to a smaller number. Lognormal distribution is widely used in fatigue life research [17]. The Weibull distribution is a continuous distribution that is widely used in reliability, and its advantage lies in its strong adaptability to various types of test data. It is widely used in reliability engineering, especially for the distribution of wear and tear failure of electromechanical products. Since it can easily infer its distribution parameters using probability values, it is widely used in data processing of various life tests [15]. The log-Weibull distribution is a special correlation distribution of the Weibull distribution which is obtained by taking the logarithm of x and obeying the Weibull distribution. The effect is similar to the lognormal distribution, and the larger number can be reduced by logarithmic transformation [18].
Table 1 is the equation for the distribution function and density function of the four distributions.
Because the design life of IGBTs is generally more than 20 years, aging data is generally obtained by means of accelerated fatigue test (AFT) [19].
In Figure 2, G1 is auxiliary GTO, I1 is a fatigue current and I2 is used to calculate Tj. R and C form an absorption branch. If G1 is turned on, I1 flows through the IGBT part. If G1 is turned off, G2 provides the freewheeling route of I1. C1 and C2 produce I1 and I2, respectively. During the experiment, aging current I1 only flows through IGBT. So, it can be called SAFT.
During AFT, the ambient temperature is approximately 25 °C. I2 is 1.5 A, which is no more than 0.1% of the DUT rated current (1500 A) as a rule of thumb. The cycle time is 12 s with a duty cycle of 48%. During the AFT, DUT is placed in the actual TC power module and is then connected to the test platform. In AFT, F is calculated from the observed UCES [7], as shown in Equation (3).
F = U CESat U CESat 0 U CESat 0
In Equation (3), UCES0 is the initial saturation voltage of a new DUT, and UCES is the current value. When UCES rises to 1.2 times that of UCES0, F is defined to be 1.
With the data from SAFT, the corresponding F is calculated. Four different distribution models are fitted by the LSM method. As for the normal distribution, it gives:
F ( x ) = 1 35500 2 π x e ( t 309254 ) 2 2 × 35500 2 d t , μ = 309254 , σ = 35500
As for the lognormal distribution, it gives:
F ( x ) = 1 0.127 2 π ln x e ( t 12.665 ) 2 2 × 0.127 2 d t , μ = 12.665 , σ = 0.127
As for the Weibull distribution, it gives:
F ( x ) = 1 e ( x 317254 ) 12.98 , λ = 317254 ,   k = 12.98
As for the Logarithmic Weibull distribution, it gives:
F ( x ) = 1 e ( ln x 165 ) 12.667 , λ = 165 ,   k = 12.667
In order to select the most suitable one of the four distributions, we plot the curve of the experimental data and four distributions, and the density functions are also plotted. The results are shown in Figure 3 and Figure 4. It should be noted that in Figure 3 and Figure 4, the curve of the Weibull distribution coincides with the curve of the Logarithmic Weibull distribution. In Figure 3 and Figure 4, the horizontal axis is the number of cycle times, and the vertical axis is cumulative failure rate, which is calculated according to DUT saturation voltage drop [13,17].
From Figure 4, it seems that the Weibull distribution is closer to the practical fatigue process (hence the fault transition process) compared with the normal distribution, because the latter one shows a larger error. What is more, the lognormal distribution is basically consistent with normal distribution; hence the same fitting accuracy is obtained. However, the transient component of failure density in the case of the lognormal distribution function is too high; therefore, it is not suitable to describe the failure density of high-power IGBTs. The logarithmic Weibull distribution has the same characterization accuracy as the Weibull distribution, however its calculation process is much more complicated. Finally, the Weibull distribution is the most suitable one.

2.2. Effect of FWD Actions on Device Aging

It is not enough to consider IGBT only in the fatigue model. In the device, current can flow in two directions, as shown in Figure 5, where C and E are the external terminals of the device. In Figure 5a, current flows through IGBT. In Figure 5b, current flows through FWD. In AFT, the DUT is heated by the test current flowing through it. However, conventional AFT only offers the test current in one direction, i.e., the direction shown in Figure 5a [13,14,15], considering that IGBT is more fragile. In conventional AFT (hence SAFT, single-directional AFT), the aging of IGBT is equal to the aging of DUT. However, the current through FWD also generates heat, hence FWD actions also cause DUT to age. SAFT should be improved to bidirectional AFT (BAFT) where the test current flows through both the IGBT and FWD. The data that is obtained by BAFT accurately shows the aging effects caused by IGBTs and FWD. Such an effect is shown more clearly in junction temperature variation, as shown in the simulation of the TC operation below.
Figure 6 shows TC topology, where Q11~Q32 are IGBTs, D11~D32 are FWDs, QB, RB and DB constitute the brake energy consumption branch. Q11~Q32 and D11~D32 are usually packaged in six modules onto the heat sink. QB and DB is not used frequency, therefore only the fatigue of Q11~Q32 and D11~D32 are considered here.
The internal temperature of the IGBTs and FWDs under the given operating conditions was obtained in an electro-thermal simulation of the TC model. In the simulation, the EMU went through two operating cycles. The simulation results are shown in Figure 7. The horizontal axis in Figure 7 is time, and the vertical axis is junction temperature. The switching device used is 1500 A/3300V, the switching frequency is 57 Hz~1000 Hz and the ambient temperature is 25 °C.
In Figure 7, the maximum variation of FWD junction temperature (ΔTj2) is 56.3 °C, and that of IGBT (ΔTj1) is 39.7 °C. Junction temperature variation severely affects device lifetime [20,21]. According to the experimental data of LESIT [22] and CIPS08 [23], the junction temperature change and the average junction temperature are the key factors affecting the aging process. Table 2 shows the maximum values of ΔTj2 andΔTj1 under different load conditions [24]. ΔTj2 and ΔTj1 are affected by various factors, such as power loss of IGBT and FWD, heat dissipation condition and thermal resistance. In Table 2, ΔTj2 and ΔTj1 increase as the vehicle load increases (which appears in the form of load current, i.e., the fatigue current), which means that there is certain one-to-one correspondence between ΔTj1, ΔTj2 and the load current.
In Figure 7, it should be noted that ΔTj1 is higher in the traction phase than ΔTj2 and lower in the braking phase. What is more, ΔTj2 grows much faster than ΔTj1 in the high-speed braking phase. This stems from the fact that more TC power output is required in the high-speed braking phase than in the traction phase, and such braking power mainly flows through the FWD instead of the IGBT [6]. Therefore, the role of FWD should not be ignored during the fatigue of the switching device.
As stated above, the IGBT ages faster with the FWD actions, i.e., the fatigue process is accelerated. If the aging state of the DUT is still represented by the aging state of the IGBT (as is the case in the existing references [2,3,4]), this acceleration effect of the FWD action must be considered.

2.3. BAFT and the Analysis of BAFT Results

The fatigue model reveals the relationship between the main factors affecting fatigue accumulation and the degree of fatigue. The load condition of the vehicle directly affects the fatigue of the switching device. The TC current output reflects the load level and this current flow completely through the switching device. Therefore, there is a one-to-one correspondence between the TC current output and the fatigue level of the switching device, and there is a one-to-one correspondence between the device current and the fatigue level and between the device current and the device lifetime. The correspondence is shown in Equation (8).
N f = f ( I eq )
Nf is the number of cycles the device can withstand before failure. We denote Nf as the lifetime of the device; Ieq is the equivalent current amplitude of the device.
The acceleration effect of the FWD action means that the life of the device is shortened, as shown in Equation (9).
N f 2 = n N f 1
In Equation (7), Nf2 is the practical lifetime with FWD action, while Nf1 is the lifetime of only considering IGBT action; n is a coefficient indicating an acceleration effect.
However, n is difficult to obtain through experimental testing, so other ways need to be considered to express the aging acceleration effect of FWD on power modules. For better applicability, we present the Equation (8) in its current form, as shown in Equation (10).
I eq 2 = i G 2 + μ 2 i D 2
In Equation (10), iG is the current of the IGBT, iD is the current of FWD and μ is the acceleration factor. The fatigue of semiconductor devices is caused by the accumulation of thermal stress effects, so the current is in the form of RMS (Root-Mean-Square), and there is a square operation.
In order to obtain μ, we propose a new AFT (BAFT) with bi-directional fatigue current. The BAFT test platform is shown in Figure 8.
In Figure 8, the device functions are similar to those in the single-directional accelerated fatigue test platform of Figure 2, except that when G2 and G4 are turned on, I1 flows through the FWD portion inside the DUT. In BAFT, the number of cycles of IGBT and FWD is equal, so only half of the total number of actions is considered.
For comparison, a bi-directional experiment was performed at I1 of 1500 A. Based on the data obtained from the bi-directional experiment, the cumulative failure rate curve of BAFT is shown together with the cumulative failure rate curve of SAFT at 1500 A in Figure 9. In Figure 9, the horizontal axis is the number of cycle times and the vertical axis is cumulative failure rate. It can be seen that the bi-directional experiment significantly accelerates the aging of the entire switching device.
Since the fatigue process of the switching device follows the Weibull distribution, the Weibull distribution is fitted by the least squares method according to the bi-directional aging experimental data under 1500 A, and the mathematical analysis relationship between the BAFT and the number of cycles x can be fitted to:
F B ( x ) = 1 e ( x 283769 ) 15.05
The result curve of SAFT is fitted as:
F S ( x ) = 1 e ( x 317254 ) 12.98
The mean of the Weibull distribution:
E = λ Γ ( 1 k + 1 )
In Equation (13), Γ is a gamma function, and λ and k are parameters of the Weibull distribution.
To derive μ, we define an intermediate variable m which is derived from the quotient between the mean of FB (x) and FS (x), as shown in Equation (14).
m = λ 2 Γ ( 1 + 1 k 2 ) λ 1 Γ ( 1 + 1 k 1 )
In Equation (14), λ1 = 283769, k1 = 15.05, λ2 = 317254, k2 = 12.98. Using these values, m can be calculated as 1.1125.
From the difference between BAFT and SAFT, it is given that:
{ I eq = q i G i G = i D
Consider iG = iD = 1500A and m = 1.1125. Equation (10) shows that FWD accelerates the aging speed by 23.77%.

2.4. Fault Model Fitted by Weighted LSM

However, it is not enough to analyze the experimental data of bi-directional aging under a single current. In order to obtain a more consistent distribution function with the practical trend of failure rate, this paper will analyze the bi-directional aging under different current values. Take the bi-directional 900A experiment as an example. The data collected from the experiment are plotted in Figure 10. In Figure 10, the horizontal axis is the number of cycle times and the vertical axis is cumulative failure rate.
It should be noted that when the test current decreases, the aging speed decreases sharply [25]. It is too difficult to obtain complete aging data through experiments. Therefore, only the part with a failure rate of 0.4 is used for analysis.
Through the previous analysis, Weibull distribution is selected to fit the failure rate distribution. The least squares method can be used to estimate the parameters of non-linear regression analysis. The traditional methods of least squares are direct search method, lattice search method, Gauss-Newton method, Newton-Raphson method, etc. [26]. Direct search method and lattice search method are inefficient and seldom used in practical engineering [27]. Gauss-Newton method is also called linearization method. Newton-Raphson method is an improvement of Gauss-Newton method. However, in some cases, these two methods may require much iteration to converge, and in a few cases, they may not converge at all [26].
To a certain extent, the problem of the above method can be avoided by using a genetic algorithm for least squares. GA (genetic algorithm) is an adaptive stochastic search method that is based on natural selection and biogenetic theory [28]. Firstly, the solution of the problem to be solved is regarded as a number of the individuals in the population. The individual in the population is coded and expressed as a string. Then, according to the fitness of the individuals, the selection, crossing and variation operations between the individuals are simulated to generate new solutions. In the process of population cyclic evolution, individuals gradually evolve to the state of the approximate optimal solution; therefore, the optimal solution of the problem can be obtained. The flow chart is shown in Figure 11.
The most important part of GA is fitness function. When GA simulates the selection of nature, the adaptation of each generation of population to conditions corresponds to a necessary parameter, which determines how much probability the individuals in this population participate in the heredity. The function of evaluating this probability is called fitness function. The suitability of fitness function selection will affect the whole GA. Selecting a different fitness function will result in different convergence and accuracy of the whole algorithm [29].
Aiming at the fitting of Weibull distribution in this paper, the purpose of our fitting is to find the appropriate values of λ and k, so we take λ and k as variables in the initial population.
The fitness function selected is:
f i t n e s s = i = 1 n ( 1 e ( x i λ ) k Y i ) 2
In Equation (16), xi is the number of experimental cycles and Yi is the failure rate corresponding to the number i experimental cycle in the actual experimental data.
It should be noted that GA is to solve the fitness function to a minimum value, and the corresponding variables of the minimum value as the output of the optimal solution. The concept of LSM is also used in the construction of fitness function. The purpose of the fitness function is to minimize the absolute value of the difference between the output value of the function to be fitted and the actual data. After summing the squares of the difference between each output value and the practical data, we can calculate a function that is closest to the curve that was formed by the practical data.
Because of the existence of crossing and variation operations, GA has randomness in the process of finding the optimal solution, which can avoid the algorithm falling into local optimum to a certain extent. The search for the optimal solution by GA is general. The fitting degree of λ and k fitted by GA should be the closest to all data points, and the fitted function should be the closest to the actual curve.
Taking λ and k as input variables and Equation (16) as fitness function, the algorithm is simulated by the GA toolbox of MATLAB. After 200 generations of the calculation, the optimal value is λ = 925830, k = 5.7629. The relationship curve between the optimal solution of each generation and the number of iterations of the algorithm is shown in Figure 12. In Figure 12, the horizontal axis is the number of generation and the vertical axis is fitness value. It can be seen that the algorithm converges basically after about 10 generations.
Using GA to calculate this can get better fitting accuracy and convergence, however when facing a large number of data in IGBT acceleration test, GA takes too much time, which is the biggest problem of genetic algorithm. Therefore, we introduce a weighted least squares method.
The weighted algorithm is helpful to select the most important region of a set of data. Because the aging process is too slow and the data that are obtained are less when the current is small, in order to describe the average level of bi-directional aging of IGBT under low current (900 A), the LSM method of weighted residual is introduced in this paper. The LSM method of weighted residuals can conveniently determine which interval in the whole definition domain is the most concerned.
The standard LSM is to form a hyper plane in two dimensions (single independent variable, single dependent variable or SISO) by some functions. The sum of the distances between the hyper plane and the points in the sample space is the smallest. This hyper plane can be defined as s p a n { φ 1 ( x ) , φ 2 ( x ) ... , φ n ( x ) } ; the equation at any point on the extended hyper plane can be defined as:
f ( x ) = a 1 φ 1 ( x ) + a 2 φ 2 ( x ) ... + a n φ n ( x )
In fact,
y i = f ( x i ) , ( x i , y i ) C
In Equation (17), C is the set of sample spaces.
In weighted LSM, we define residual as the sum of the distances between points and the hyper planes.
Ψ = i = 1 m f ( x i ) y i 2
In the improved LSM, the definition
Ψ = i = 1 m α i f ( x i ) y i 2 , i = 1 m α i = 1
where α i is defined as the quotient between the number of abscissa x i ( X i , X i + Δ X ) of sample point ( x i , y i ) and the total number of sample points.
Δ X is usually determined by experiment or experience. For IGBT aging failure rate fitting, it is actually the resolution of junction temperature fluctuation (the relative relationship between failure rate and junction temperature fluctuation) or the time resolution (the relationship between failure rate and time).
Ψ a h = a h i = 1 m α i f ( x i ) y i 2 Ψ a h = a h { i = 1 m α i [ ( j = 1 n a j φ j ( x i ) ) 2 ] } 2 a h { i = 1 m [ y i j = 1 n a j φ j ( x i ) ] }
Ψ a h = 0 j = 1 n { a j i = 1 m [ α i φ h ( x i ) φ j ( x i ) ] } = [ i = 1 m y i α i φ h ( x i ) ]
If defined j = 1 n [ α i φ h ( x i ) φ j ( x i ) ] = < α φ j , φ h > , i = 1 m y i α i φ h ( x i ) = < α f , φ h >
This is actually the definition of product in vector space.
Ψ a h = 0 j = 1 n [ a j < α φ j , φ h > ] = < α f , φ h >
{ Ψ a 1 = 0 Ψ a 2 = 0 ... Ψ a n = 0 ( < α φ 1 , φ 1 > < α φ 2 , φ 1 > ... < α φ n , φ 1 > < α φ 1 , φ 2 > < α φ 2 , φ 2 > ... < α φ n , φ 2 > ... ... ... ... < α φ 1 , φ n > < α φ 2 , φ n > ... < α φ n , φ n > ) ( a 1 a 2 ... a n ) = ( < α f , φ 1 > < α f , φ 2 > ... < α f , φ n > )
Note: Here A can be given according to the normal distribution.
According to the weighted LSM theory, we select different weights and obtain three curves, as shown in Figure 13. In Figure 13, the horizontal axis is the number of cycle times and the vertical axis is cumulative failure rate.
Curve 1 is obtained by maximizing the weight of the interval [7 × 105, 8 × 105] of experiments’ number, parameter λ = 868359, k = 8.99, mean value calculated by Equation (5), E1 = 822222.
Curve 2 is obtained by maximizing the weight of interval [4 × 105, 7 × 105]. The parameters λ = 918264, k = 5.96, E2 = 851580.
Curve 3 is obtained by maximizing the weight of interval [1 × 105, 5 × 105]. The parameters λ = 987945, k = 4.02, E3 = 895730. Compared with the experimental data curves, the three curves fit the original data most accurately in the interval with the largest weights.
From the point of view of the relationship between failure rate and time, in fact, what we need to focus on is not only the accuracy of characterization of failure rate in the whole life cycle of products, however also the accuracy of the most product failure processes. Accurately grasping the change of the latter is more conducive to accurately predicting the state of aging damage of components. According to this idea, observing the experimental data curve under 900 A, the failure rate rises rapidly in the interval where the number of tests is [4.5 × 105, 7.5 × 105], so we think that curve 2 fits the aging situation under 900 A most accurately. By fitting the data with this weighted least squares method, the function corresponding to the above curve 2 is the closest to the actual distribution function.
In order to validate the distribution function fitted by weighted LSM, we compare the fitting results with those of GA. The function curve obtained is compared with the curve that was fitted by weighted LSM as shown in Figure 14, where the horizontal axis is the number of cycle times and the vertical axis is cumulative failure rate.
According to Figure 14, the curve of the function that was obtained by GA and curve 2 almost coincides in the interval of cycle number [4 × 105, 8 × 105], which proves that the weighted LSM can more accurately fit the distribution function. Compared with the GA, the computational complexity of weighted LSM is greatly reduced and the computation time is shortened by 27% compared with the GA. So, the weighted LSM is more suitable for fitting the distribution.

2.5. Analytical Fatigue Model

Different BAFTs are based on different I1, as shown in Figure 15. The horizontal axis in Figure 15 is the number of cycle times, and the vertical axis is cumulative failure rate. Under I1 at 1200 A and 900 A, the mathematical model between cumulative failure rate F and cycle number x is as follows:
F B 1200 ( x ) = 1 e ( x 508785 ) 8.01
F B 900 ( x ) = 1 e ( x 918264 ) 5.96
Obtained from Equation (12)
F B 1500 ( x ) = 1 e ( x 317254 ) 12.98
In Equations (21)–(23), FB1500(x), FB1200(x) and FB900(x) are the cumulative failure efficiencies of I1 at 1500A, 1200A and 900A, respectively, and their mean values are 304890, 479170 and 851580. The corresponding equipment service life Nf can be expressed in terms of the mean value.
In addition, the ambient temperature Ta should also be considered [30,31]. According to Equation (9), the Nf expression considering Ta is obtained.
N f = f ( I eq ) × g ( T a )
In Equation (24), g is specific functions.
Finally, after the fitting calculation, Equation (9) is expressed as:
N f = ( 1.17 I 1 2 3697.11 I 1 + 3.23 × 10 6 ) × e k ( 1 T a 1 25 )
In Equation (25), k is an Arrhenius coefficient. According to Arrhenius’s law (a temperature increase of 10 °C means that the service life is shortened by half), k is calculated to be 60.65. The experimental and fitting curves are shown in Figure 16.
Considering the RMS calculation scheme, the amplitude is 2 times of the effective value, so:
I 1 = 2 I eq
From Equation (26):
N f = ( 2.34 I eq 2 5228.35 I eq + 3.23 × 10 6 ) × e 60.65 ( 1 T a 1 25 )
The data obtained in the experiment and the curve that was obtained by mathematical fitting is shown in Figure 16. The vertical axis in Figure 16 is the cycle counts that a DUT could survive before failure. The horizontal axis is the amplitude value of the test current. As shown in Figure 16, the obtained curve conforms to the obtained data very well.
By Equation (27), when Ieq and Ta are known, the service life or failure rate of the IGBT can be predicted. Equation (27) has been subjected to actual field data. The example we consider here is a metro line application in Shenzhen, a southern-China city, with the annual temperature average of 23 °C. The Ieq has been recorded to be 489 A in a single year; therefore, Nf is calculated to be 2,197,460. Considering that a train normally operates 204 cycles per day and normally 300 days per year, the train is capable of operating around 36 years under the same working conditions.
According to statistics offered by metro operating agency, the initial failure rate of a new module is 157 FIT. Given that they hold that the failure rate rise of 200% means that the switching devices in a batch have failed, the statistical service life of the applied power module in the specific metro line is around 47 years (offered by the agency). This means that the approach that was proposed has a prediction error of around 23%, which is within project tolerance (normally 20%~30%). It should be noted that the lifetime predicted by our model in Equation (27) is more likely to be true and more convincing, because a metro train is designed to operate for no more than 30 years, and the 47-year lifetime prediction does not take into consideration the acceleration of the aging process, as is shown by a bathtub curve.

3. Conclusions

The object of the fault model in this paper is the fault model of the high-power switching device that was applied to rail transit. Based on a new test method using bidirectional current aging, we propose a method for extracting the acceleration factor of the freewheeling diode in the switching device for the entire module. Then, the weighting factor is introduced in the traditional least squares method, which improves the processing efficiency of the test results and the fitting accuracy of the model. Finally, an IGBT fault prediction model was presented by using the aging current and ambient temperature as sensitive factor.

Author Contributions

L.W. and S.L. designed the experiment, S.L. and L.W. performed the experiment, R.Q. did the data analysis, C.X. contributed analysis tools, L.W. wrote the paper.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities, grant number 2018JBM060.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations and Nomenclature

FWDFree Wheeling Diode
EMUElectric Multiple Unit
IGBTInsulated Gate Bipolar Transistor
TCTraction Converter
CTECoefficient of Thermal Expansion
DUTDevice Under Test
AFTAccelerated Fatigue Test
GTOGate-Turn-Off thyristor
SAFTSingle-directional Accelerated Fatigue Test
BAFTBi-directional Accelerated Fatigue Test
EMIElectromagnetic Interference
GAGenetic Algorithm
SThermal stress
RthThermal resistance
Rth-c_1Thermal resistance between the IGBT die and device case
Rth-c_2Thermal resistance between the FWD die and device case
ΔCTEDifference of CTE
ΔTJunction temperature variation
PHeating power
FCumulative failure rate
I1, I2the currents used in SAFT and BAFT
UCES0Initial saturation voltage drop of the device
UCESSaturation voltage drop of the device
i1, i2, i3Current outputs of TC
UdcDC voltage input of TC
TjJunction temperature
Tj1Junction temperature of IGBT
Tj2Junction temperature of FWD
ΔTj1Junction temperature variation of IGBT
ΔTj2Junction temperature variation of FWD
IeqEquivalent current amplitude through the device or fatigue current
NfThe service life of the device
n, mCoefficients to show the accelerating effect from FWD actions on device
μAcceleration factor
iGIGBT current
iDFWD current
FB(x), FS(x), FB1500(x), FB1200(x), FB900(x)the function relationship between F and the number of cycle times x under various conditions
λ1, k1, λ2, k2Parameters in Weibull distribution
TaAmbient temperature

References

  1. Enrico, P.; Riccardo, T.; Pietro, C. Fault Current Detection and Dangerous Voltages in DC Urban Rail Traction Systems. IEEE Trans. Ind. Appl. 2017, 53, 4109–4115. [Google Scholar] [Green Version]
  2. Shen, C.C.; Hefner, A.R.; Berning, D.W.; Bernstein, J.B. Failure dynamics of the IGBT during turn-off for unclamped inductive loading conditions. IEEE Trans. Ind. Appl. 2000, 36, 614–624. [Google Scholar] [CrossRef]
  3. Liu, Y.; Xiang, D.; Fu, Y. Prognosis of Chip-loss Failure in High-power IGBT Module by Self-testing. IEEE Ind. Electron. Soc. 2016, 6836–6840. [Google Scholar] [CrossRef]
  4. De Vega Angel, R.; Ghimire, P. Test setup for accelerated test of High Power IGBT modules with online monitoring of V-ce and V-f voltage during converter operation. In Proceedings of the 2014 International Power Electronics Conference, Hiroshima, Japan, 18–21 May 2014; pp. 2547–2553. [Google Scholar]
  5. Thebaud, J.M.; Woirgard, E.; Zardini, C. Strategy for designing accelerated aging tests to evaluate IGBT power modules lifetime in real operation mode. IEEE Trans. Compon. Packag. Technol. 2003, 26, 429–438. [Google Scholar] [CrossRef]
  6. Li, L.; Xu, Y.; Li, Z. The effect of electro-thermal parameters on IGBT junction temperature with the aging of module. Microelectron. Reliab. 2016, 66, 58–63. [Google Scholar] [CrossRef]
  7. Mei Sundaramoorthy, V.K.; Bianda, E.; Bloch, R. A study on IGBT junction temperature (Tj) online estimation using gate-emitter voltage (Vge) at turn-off. Microelectron. Reliab. 2014, 54, 2423–2431. [Google Scholar] [CrossRef]
  8. Coquery, G.; Lallemand, R. Failure criteria for long term Accelerated Power Cycling Test linked to electrical turn off SOA on IGBT module. A 4000 hours test on 1200A-3300V module with AlSiC base plate. Microelectron. Reliab. 2000, 40, 1665–1670. [Google Scholar] [CrossRef]
  9. Li, H.; Wang, Y.; Ashle, P. Modeling, Analysis and Implementation of Junction-to-Case Thermal Resistance Measurement for Press-Pack IGBT Modules. In Proceedings of the 19th International Conference on Electronic Packaging Technology (Icept), Shanghai, China, 8–11 August 2018; pp. 1605–1610. [Google Scholar]
  10. Ding, J.; Zhang, P.; Li, J. Fatigue Life Prediction of IGBT Module for Metro Vehicle Traction Converter Based on Traction Calculation. In Proceedings of the 2015 IEEE 11th International Conference on Power Electronics and Drive Systems, Sydney, Australia, 9–12 June 2015; pp. 1116–1121. [Google Scholar]
  11. Barlini, D.; Ciappa, M.; Castellazzi, A. New technique for the measurement of the static and of the transient junction temperature in IGBT devices under operating conditions. Microelectron. Reliab. 2006, 46, 1772–1777. [Google Scholar] [CrossRef]
  12. Wuest, F.; Trampert, S.; Janzen, S. Comparison of temperature sensitive electrical parameter based methods for junction temperature determination during accelerated aging of power electronics. Microelectron. Reliab. 2018, 88–90, 534–539. [Google Scholar] [CrossRef]
  13. Daniel, A.; Ibanez, M.; Ainhoa, G.; Jose, M.E. Analysis of the Results of Accelerated Aging Tests in Insulated Gate Bipolar Transistors. IEEE Trans. Power Electron. 2016, 31, 7953–7962. [Google Scholar]
  14. Wang, L.; Xu, J.; Wang, G. Lifetime estimation of IGBT modules for MMC-HVDC application. Microelectron. Reliab. 2018, 82, 90–99. [Google Scholar] [CrossRef]
  15. Cao, J.; Cheng, K. Common life distribution. In Introduction to Reliability Mathematics; Higher Education Press: Beijing, China, 2006; pp. 7–37. [Google Scholar]
  16. Liu, W. Several Probability Distributions in Reliability Theory. In Mechanical Reliability Design; Tsinghua University Press: Beijing, China, 1996; pp. 50–91. [Google Scholar]
  17. Ui-Min, C.; Soren, J.; Frede, B. Advanced Accelerated Power Cycling Test for Reliability Investigation of Power Device Modules. IEEE Trans. Power Electron. 2016, 32, 8371–8386. [Google Scholar]
  18. Jiang, R.; Wang, T. Log-Weibull Distribution as A Lifetime Distribution. In Proceedings of the 2013 International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering, Chengdu, China, 15–18 July 2013; pp. 813–816. [Google Scholar]
  19. Dornic, N.; Ibrahim, A.; Khatir, Z.; Tran, S.H.; Ousten, J.P. Analysis of the degradation mechanisms occurring in the topside interconnections of IGBT power devices during power cycling. Microelectron. Reliab. 2018, 12, 462–469. [Google Scholar] [CrossRef]
  20. Song, W.; Ma, J.; Zhou, L.; Feng, X. Deadbeat Predictive Power Control of Single-Phase Three-Level Neutral-Point-Clamped Converters Using Space-Vector Modulation for Electric Railway Traction. IEEE Trans. Power Electron. 2016, 31, 721–732. [Google Scholar] [CrossRef]
  21. Tian, B.; Qiao, W.; Wang, Z.; Tanya, G.; Hudgins Jerry, L. Monitoring IGBT’s Health Condition via Junction Temperature Variations. In Proceedings of the 2014 IEEE Applied Power Electronics Conference and Exposition (APEC), Fort Worth, TX, USA, 16–20 March 2014; pp. 2550–2555. [Google Scholar]
  22. Miao, Y.; Lei, W.; Li, S.; Lv, X. Influence of Inverter DC Voltage on the Reliability of IGBT. In Proceedings of the 6th International Conference on Computer-Aided Design, Manufacturing, Modeling and Simulation (CDMMS 2018), Busan, Korea, 4–15 April 2018. [Google Scholar]
  23. Josef, L. IGBT-Modules: Design for Reliability. In Proceedings of the 13th European Conference on Power Electronics and Applications, EPE’09, Barcelona, Spain, 8–10 September 2009; pp. 6045–6047. [Google Scholar]
  24. Shi, H.; Wu, P.; Luo, R.; Zeng, J. Estimation of the damping effects of suspension systems on railway vehicles using wedge tests. Proc. Inst. Mech. Eng. Part F J. Rail Rapid Transit 2016, 203, 392–406. [Google Scholar] [CrossRef]
  25. Vanessa, S.; Forest, F.; Huselstein, J.-J. Ageing and Failure Modes of IGBT Modules in High-Temperature Power Cycling. IEEE Trans. Ind. Electron. 2011, 58, 4931–4941. [Google Scholar]
  26. Neter, J.; Wasserman, W.; Kutner, M.H. Nonlinear Regression. In Applied Linear Regression Models; China Statistics Press: Beijing, China, 1990; pp. 518–529. [Google Scholar]
  27. Bieniasz, L.K.; Speiser, B. Use of sensitivity analysis methods in the modelling of electrochemical transients—Part 3. Statistical error uncertainty propagation in simulation and in nonlinear least-squares parameter estimation. J. Electroanal. Anal. Chem. 1998, 458, 209–229. [Google Scholar] [CrossRef]
  28. Li, Y.; Zhao, L.; Zhou, S. Review of Genetic Algorithm. Adv. Mater. Res. 2011, 179–180, 365–367. [Google Scholar] [CrossRef]
  29. Zhu, Y. Evolving dynamic fitness measures for genetic programming. Expert Syst. Appl. 2018, 109, 162–187. [Google Scholar]
  30. Held, M.; Jacob, P.; Nicoletti, G.; Scacco, P.; Poech, M.H. Fast power cycling test of IGBT modules in traction application. In Proceedings of the Second International Conference on Power Electronics and Drive Systems, Singapore, 26–29 May 1997; pp. 425–430. [Google Scholar]
  31. Bayerer, R.; Herrmann, T.; Licht, T.; Lutz, J.; Feller, M. Model for Power Cycling lifetime of IGBT Modules—Various factors influencing lifetime. In Proceedings of the 5th International Conference on Integrated Power Electronics Systems, Nuremberg, Germany, 11–13 March 2008; pp. 1–6. [Google Scholar]
Figure 1. Internal structure of the switchgear.
Figure 1. Internal structure of the switchgear.
Applsci 09 00444 g001
Figure 2. Schematic diagram of single-directional accelerated fatigue test.
Figure 2. Schematic diagram of single-directional accelerated fatigue test.
Applsci 09 00444 g002
Figure 3. Comparison between data curves of single-direction 1500 A fatigue test and four distribution fitting results.
Figure 3. Comparison between data curves of single-direction 1500 A fatigue test and four distribution fitting results.
Applsci 09 00444 g003
Figure 4. Comparison of the density function of Weibull distribution, Logarithmic Weibull distribution and normal distribution, lognormal distribution under the single-direction 1500 A fatigue test.
Figure 4. Comparison of the density function of Weibull distribution, Logarithmic Weibull distribution and normal distribution, lognormal distribution under the single-direction 1500 A fatigue test.
Applsci 09 00444 g004
Figure 5. IGBT and FWD current flow differences. (a) IGBT current flow; (b) FWD current flow.
Figure 5. IGBT and FWD current flow differences. (a) IGBT current flow; (b) FWD current flow.
Applsci 09 00444 g005
Figure 6. TC (traction converter) circuit topology.
Figure 6. TC (traction converter) circuit topology.
Applsci 09 00444 g006
Figure 7. Junction Temperature Simulation of Q32 and D32 under AW0.
Figure 7. Junction Temperature Simulation of Q32 and D32 under AW0.
Applsci 09 00444 g007
Figure 8. Schematic diagram of a bi-directional accelerated fatigue test.
Figure 8. Schematic diagram of a bi-directional accelerated fatigue test.
Applsci 09 00444 g008
Figure 9. BAFT (bidirectional AFT) and SAFT curves obtained under 1500 A.
Figure 9. BAFT (bidirectional AFT) and SAFT curves obtained under 1500 A.
Applsci 09 00444 g009
Figure 10. Curve of bi-directional aging under 900A.
Figure 10. Curve of bi-directional aging under 900A.
Applsci 09 00444 g010
Figure 11. Flow chart of genetic algorithm.
Figure 11. Flow chart of genetic algorithm.
Applsci 09 00444 g011
Figure 12. Optimal Solution Curve of Each Generation.
Figure 12. Optimal Solution Curve of Each Generation.
Applsci 09 00444 g012
Figure 13. Weibull distribution function curves fitted under different weights.
Figure 13. Weibull distribution function curves fitted under different weights.
Applsci 09 00444 g013
Figure 14. Comparison of Weibull distribution fitted by GA and curves fitted by weighted LSM (least squares method).
Figure 14. Comparison of Weibull distribution fitted by GA and curves fitted by weighted LSM (least squares method).
Applsci 09 00444 g014
Figure 15. Curves obtained under different BAFT aging currents.
Figure 15. Curves obtained under different BAFT aging currents.
Applsci 09 00444 g015
Figure 16. Experimental results of BAFT and fitting curves.
Figure 16. Experimental results of BAFT and fitting curves.
Applsci 09 00444 g016
Table 1. Distribution functions and density functions of the four distributions.
Table 1. Distribution functions and density functions of the four distributions.
NameDistribution FunctionDensity Function
Normal distribution F ( x ) = 1 σ 2 π x e ( t μ ) 2 2 σ 2 d t f ( x , μ , σ ) = { 1 σ 2 π e ( x μ ) 2 2 σ 2 , x > 0 0 , x 0
Lognormal distribution F ( x ) = 1 σ 2 π ln x e ( t μ ) 2 2 σ 2 d t f ( x , μ , σ ) = { 1 σ x 2 π e ( ln x μ ) 2 2 σ 2 , x > 0 0 , x 0
Weibull distribution F ( x ) = { 1 e ( x λ ) k , x > 0 0 , x 0 f ( x , k , λ ) = { k λ ( x λ ) k 1 e ( x λ ) k , x > 0 0 , x 0
Logarithmic Weibull distribution F ( x ) = { 1 e ( ln x λ ) k , x > 0 0 , x 0 f ( x , k , λ ) = { k λ x ( x λ ) k 1 e ( ln x λ ) k , x > 0 0 , x 0
Table 2. ΔTj of Q32 and D32 under different load conditions.
Table 2. ΔTj of Q32 and D32 under different load conditions.
Load ConditionΔTj1ΔTj2
AW039.7 °C56.3 °C
AW252.6 °C72.1 °C
AW360.9 °C87.2 °C

Share and Cite

MDPI and ACS Style

Wang, L.; Liu, S.; Qiu, R.; Xu, C. Fault Prediction Model of High-Power Switching Device in Urban Railway Traction Converter with Bi-Directional Fatigue Data and Weighted LSM. Appl. Sci. 2019, 9, 444. https://doi.org/10.3390/app9030444

AMA Style

Wang L, Liu S, Qiu R, Xu C. Fault Prediction Model of High-Power Switching Device in Urban Railway Traction Converter with Bi-Directional Fatigue Data and Weighted LSM. Applied Sciences. 2019; 9(3):444. https://doi.org/10.3390/app9030444

Chicago/Turabian Style

Wang, Lei, Shenyi Liu, Ruichang Qiu, and Chunmei Xu. 2019. "Fault Prediction Model of High-Power Switching Device in Urban Railway Traction Converter with Bi-Directional Fatigue Data and Weighted LSM" Applied Sciences 9, no. 3: 444. https://doi.org/10.3390/app9030444

APA Style

Wang, L., Liu, S., Qiu, R., & Xu, C. (2019). Fault Prediction Model of High-Power Switching Device in Urban Railway Traction Converter with Bi-Directional Fatigue Data and Weighted LSM. Applied Sciences, 9(3), 444. https://doi.org/10.3390/app9030444

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