Review of Methodologies for Structural Integrity Evaluation of Power Modules

Miyazaki, Noriyuki; Shishido, Nobuyuki; Hayama, Yutaka

doi:10.1115/1.4048038

Abstract

This paper reviews the previous research on the methodologies for evaluating structural integrity of wire bonds and die-attachments in power modules. Under power module operation, these parts are subjected to repeated temperature variations which induce repeated thermal stress due to the mismatch in coefficients of thermal expansion (CTE) of the constituent materials. Thus, thermal fatigue phenomena are critical issues for the structural integrity of power modules. In the present paper, we also deal with the evaluation methodologies for thermal fatigue in the temperatures over 200 $° C$ ⁠, which are expected operational temperatures for wide bandgap semiconductor power modules. The failure models based on the temperature range $Δ T$ widely used in the power electronics community are critically reviewed from a mechanical engineering viewpoint. Detailed discussion is given concerning the superiority of failure models based on the physical quantities such as the inelastic strain range $Δ ε_{in}$ ⁠, the inelastic strain energy density range $Δ W_{in}$ ⁠, and the nonlinear fracture mechanics parameter range $Δ T$ * over the conventional $Δ T$ -based failure models. It is also pointed out that the distributed state concept (DSC) approaches based on the unified constitutive modeling and the unified mechanics theory are promising for evaluating the structural integrity of power modules. Two kinds of test methods, a power cycling test (PCT) and a thermal cycling test (TCT), are discussed in the relation to evaluating the lifetimes of wire-liftoff and die attach cracking.

1 Introduction

Among various kinds of energy, the ratio of electric power consumption to total energy consumption increases year by year, because electric power is safe, clean, and easy to use. Power modules are utilized for electric power control and play a key role in efficient energy conversion. They are embedded in various systems such as electric power transmission systems, railway systems, automobiles, factory automation systems, home appliances, and so on. Thus, research on highly reliable and long-life power modules contributes to realizing energy-saving and sustainable society.

Let us consider an insulated gate bipolar transistor (IGBT) shown in Fig. 1. Key locations for the structural integrity are the connection between a chip and a bonding wire, a wire-bonding location, and the connection between a chip and a copper substrate by means of a die attachment material. Under power module operation, these connections are subjected to repeated temperature cycles which induce repeated variations of thermal stress due to the mismatch of coefficients of thermal expansion (CTE) of materials consisting of a power module. Thus, thermal fatigue phenomena are critical issues for the structural integrity of power modules. They are wire-liftoff, in which a bonding wire is delaminated from a chip surface, heel cracking, in which the heel part of a bonding wire is fractured, and die attach cracking, in which a crack occurs and propagates in a die attach layer. In the present review paper, we deal with these thermal fatigue phenomena. Power modules using wide bandgap semiconductors such as SiC and GaN have been developed in recent years. Their operating temperatures are expected to be over 200 °C. In this paper, we also deal with the thermal fatigue phenomena at elevated temperature over 200 °C.

Fig. 1

View large Download slide

Schematic of IGBT power device

Historically, the academic disciplines concerning structural integrity evaluation have advanced in the safety design and operation of large-scaled structures used as civil engineering structures, hull structures, aircraft structures, pressure vessels and piping, etc. Especially, fracture mechanics and fatigue strength have greatly advanced. Electronic packaging can be regarded as a small-scaled or miniature structure which has sources of stress concentration such as interfaces between dissimilar materials. Thus, methodologies for the structural integrity evaluation of large-scaled structures have been necessarily applied to electronic packaging [1]. In 1986, a Japanese automobile company decided to recall faulty automobiles causing sudden start. The reason for such an incident was the fatigue failure of solder joints in integrating circuits (ICs) of an autodrive computer. This trouble attracted attention to the structural integrity of electronic packaging. Research on the structural integrity of IC packaging has been actively performed since 1980s. On the other hand, such research for power devices has been performed since 1990s [2,3], and several review papers have been published [4–8]. According to the authors' opinion, main concern of the structural integrity of power modules is condition monitoring and remaining life estimation of in-service power devices, as shown in Refs. [5], [7], and [8].

Fewer power modules are utilized in electric power transmission systems and railway vehicles than in the automobiles and home electric appliances, and their structural integrity can be confirmed by in-service condition monitoring and periodic inspection. On the other hand, an enormous number of power modules are utilized in the automobiles and home electric appliances without proper inspection. We cannot deny that a malfunction of electronic devices including power modules in automobiles may cause a serious accident as described before, and recall of an enormous number of automobiles may be possibly caused. It is necessary to confirm the structural integrity of electronic devices in the design phase to avoid their malfunction during in-service period. The CAE-based stress analysis, fracture mechanics, and strength of materials have been applied to securing not only the structural integrity of large-scaled structures but also that of IC packaging structures.

In the present paper, we review the previous research on the methodologies for the thermal fatigue evaluation of power modules. The review is made based on a mechanical engineering viewpoint such as fracture mechanics and strength of materials. In addition, we will propose the research policy for securing the structural integrity of power modules in the design phase.

2 Failure of Bonding Wire

Here, we deal with heel cracking and wire-liftoff as the failure modes of a bonding wire.

2.1 Heel Cracking.

When a wire is connected to different components with different CTEs, it experiences relative motion during temperature variation. Under the repeated temperature cycles, high local stress is repeatedly generated at the heel part of a bonding wire, and heel cracking is caused as a fatigue phenomenon. In this case, the number of cycles to fatigue failure (fatigue life or lifetime)

N_{f}

is given by the Coffin–Manson law [9,10] as follows [11–13]:

N_{f} = C {(Δ ε_{p})}^{p}

(1)

where

Δ ε_{p}

is the plastic strain range per cycle, and C and p are material constants. At elevated temperature, creep strain is induced in addition to plastic strain, the Coffin–Manson law, Eq. , is modified by changing

Δ ε_{p}

into the inelastic strain range per cycle

Δ ε_{in}

given by

Δ ε_{in} = Δ ε_{p} + Δ ε_{c}

(2)

where $Δ ε_{c}$ is the creep strain range per cycle. $Δ ε_{p}$ and $Δ ε_{c}$ can be obtained by the inelastic finite element analysis considering both plastic and creep strains.

Another methodology for predicting the lifetime of heel cracking is an energy approach [14]. Celnikier et al. propose the following equation for the lifetime estimation of an aluminum ribbon [14]:

N_{f} (I) = w_{p}^{cr} / Δ w_{p} (I)

(3)

where

{Δ w}_{p} (I)

is the plastic energy density per cycle under current I, and

w_{p}^{cr}

is the limit of the plastic energy density for heel cracking determined from

N_{f} (50 A) = 250, 000

cycles. In Ref. [14], the lifetimes of heel cracking obtained from Eq. are compared with those obtained from the Coffin–Manson law, Eq. . Although the authors of Ref. [14] conclude that the results obtained from both approaches are in good agreement, the agreement is not good from a quantitative viewpoint. Such discrepancy results from the method for determining

w_{p}^{cr}

⁠. If Eq. is valid for an arbitrary value of I, the following failure model should be true for an arbitrary value of I:

N_{f} (I) = C {({Δ w}_{p} (I))}^{- 1}, (C : constant)

(4)

As shown in Eq. , the failure model based on

{Δ w}_{p} (I)

should be generally the following form:

N_{f} (I) = C {({Δ w}_{p} (I))}^{p}

(5)

There is no logical necessity of p = –1 for an arbitrary value of I.

References [4] and [5] state that the heel cracking of a bonding wire rarely occurs in advanced IGBT multichip modules, although it can be observed mainly after long endurance tests in the case where the ultrasonic bonding process is not optimized, and conclude that the main failure mode of a bonding wire is wire-liftoff. According to the recent experimental research [15], wire-liftoff always occurs in the IGBTs molded by silicone-gel, while, in the IGBTs molded by resin, wire-liftoff occurs for large temperature range (⁠ $Δ T = 80 ° C, 100 ° C$ ⁠) and both wire-liftoff and heel cracking is possible to occur for small temperature range (⁠ $Δ T = 60 ° C$ ⁠). Resin constrains the movement of a bonding wire more strongly than silicone-gel, so that the stress and strain in a bonding wire are greatly different between the silicone-gel molded IGBTs and the resin molded IGBTs. This is the reason why the failure mode of a bonding wire is different between these two types of IGBTs. Heel cracking has been recognized again when a molding material is changed from silicone-gel to resin. For the resin-molded IGBTs, similar heel cracking is observed by Choi et al. [16] and Zeng et al. [17].

The residual stresses induced in a bonding wire during the wire-bonding process has an impact on the lifetime of heel cracking. The scientific knowledge of the impact of residual stresses on fatigue crack propagation rate was gained by investigating the midair explosion accident of Comet I, the first commercial jet airliner, over the Mediterranean Sea, which occurred in 1954 [18]. The following scientific knowledge was gained from the accident investigation. The fatigue crack propagation rate accelerates or decelerates according as the residual stress component perpendicular to the crack propagation direction is tensile or compressive. It is therefore expected that crack initiation lifetime becomes longer under residual compressive stress and shorter under residual tensile stress. According to Ref. [14], the residual stress induced in a bonding wire during the wire-bonding process has an impact on reducing heel cracking lifetime. This implies that the region around heel cracking is dominated by residual tensile stress. The wire-bonding process is composed of a prepress step, in which a wire is attached and pressed on the surface of a chip/substrate and an ultrasonic step, in which ultrasonic vibration is added to the bonding part of a wire, keeping a prepress load constant. We need information on the residual stress state around the location where heel cracking occurs to discuss the impact of residual stress on the lifetime of heel cracking. It seems that the authors of Ref. [14] did not carry out the stress analysis of the wire-bonding process to obtain a residual stress state around the location where heel cracking occurs. Conclusively sufficient information is not provided in Ref. [14] to discuss the relation between the residual stress induced during the wire-bonding process and the lifetime of heel cracking. Examples of the stress analysis and the stress measurement during the wire-bonding process are given in Refs. [19] and [20].

2.2 Wire-Liftoff

2.2.1 Conventional Failure Models.

As described in the Sec. 2.1, wire-liftoff is a main failure mode for aluminum wires of IGBTs, and a lot of papers have been published on this failure mode. In these papers, thermal fatigue life has been related to a temperature ranges. The fatigue life

N_{f}

in a wire-liftoff phenomenon has been represented by two parameters, the temperature range per cycle

Δ T

and the mean temperature T_m [3]. In a mathematical expression, it is given by

N_{f} = C_{1} {(Δ T)}^{p_{1}}

(6)

where $p_{1}$ is a material constant, and C₁ is a function of T_m. Moreover, it is expected that the temperature swing duration $t_{Δ T}$ has an impact on the lifetime prediction [21], and C₁ is a function of both T_m and $t_{Δ T}$ ⁠.

From a mechanical engineering viewpoint, the fatigue life

N_{f}

is given by the Coffin–Manson law as

N_{f} = C_{2} {(Δ ε_{p})}^{p_{2}}

(7)

If the elastic strain range per cycle

Δ ε_{e}

is very small compared with the plastic strain range per cycle

Δ ε_{p}

⁠, Eq. is equivalent to the Coffin–Manson law under a pure thermal fatigue condition without mechanical loading [3]. As pointed out by Zhang [22], Eq. gives erroneous results when

Δ ε_{e}

cannot be neglected in comparison with

Δ ε_{p}

⁠, that is, in the case of high cycle fatigue. At elevated temperature, creep strain is induced in addition to plastic strain. In such a case, the Coffin–Manson law, Eq. , is modified using the inelastic strain range per cycle

Δ ε_{in}

(=

Δ ε_{p} + Δ ε_{c})

as follows [23]:

N_{f} = C_{3} {(Δ ε_{in})}^{p_{3}}

(8)

where $p_{3}$ and C₃ are material constants. In the thermal fatigue of a bonding wire subjected to repeated temperature variations, Eq. (8) should be used if the effect of creep strain cannot be neglected. Unlike $Δ T$ in Eq. (6), $Δ ε_{in}$ in Eq. (8) is a physical quantity difficult to measure experimentally, but its accurate value can be obtained from the finite element analysis considering both plastic and creep deformations.

The $Δ T$ -based failure models have been developed and used exclusively to estimate the wire-liftoff lifetime for power modules by European researchers. That is not a consensus for researchers engaged in the reliability studies of IC electronic packaging, in which Eq. (8) is exclusively utilized. Why does there such difference exist? In the present authors' opinion, main concern is directed to monitoring power modules under operation in the power electronics community, while the lifetime estimations of IC packages in the design phase are of crucial importance in the IC community. An enormous number of IC packages are supplied to the market after a short R&D period and their services periods are relatively short compared with power modules'. From a viewpoint of economic risk management, the lifetime estimation of IC packages in the design phase is more important than their in-service monitoring. $Δ T$ is easy to measure for monitoring power modules under operation, and various failure models are proposed based on Eq. (6), as shown in Ref. [22]:

Coffin–Manson model [24]
$N_{f} = A {Δ T}^{- n}$
(9)
General Coffin–Manson model [25]
$N_{f} = {A (Δ T - Δ T_{0})}^{- n}$
(10)
Modified Coffin-Manson model [3]
$N_{f} = A {Δ T}^{- n} \exp (E_{a} / k_{b} T_{m})$
(11)
General modified Coffin–Manson model [22]
$N_{f} = {A (Δ T - Δ T_{0})}^{- n} exp (E_{a} / k_{b} T_{m})$
(12)
Bayerer's model [26]
$N_{f} = A {Δ T}^{- n} \exp (β_{2} / T_{min}) t_{on}^{β_{3}} I^{β_{4}} V^{β_{5}} D^{β_{6}}$
(13)
General Bayerer's model [22]:
$N_{f} = {A (Δ T - Δ T_{0})}^{- n} exp (β_{2} / T_{min}) t_{on}^{β_{3}} I^{β_{4}} V^{β_{5}} D^{β_{6}}$
(14)

Here, k_b is the Boltzmann constant. In these failure models shown in Eqs. (10)–(14), additional parameters are included to accommodate Eq. (9) to various operating conditions. Such additional parameters are the temperature range $Δ T_{0}$ corresponding to the elastic strain range, the mean absolute temperature T_m, the minimum absolute temperature $T_{min}$ ⁠, the power-on-time t_on, the current per wire bond I, the blocking voltage of the chip V, and the diameter of bonding wire D. The additional parameters are needed, because the effects of plastic and creep strains on the lifetime of power modules are not represented only by the temperature range $Δ T$ ⁠. The constants A, n, E_a, $β_{2}$ ⁠, $β_{3}$ ⁠, $β_{4}, β_{5}$ ⁠, and $β_{6}$ are those determined from experimental data. As the number of additional parameters increases, we need more experimental data to determine the constants relevant to the additional parameters.

In the failure model based on $Δ ε_{in} (= Δ ε_{p} + Δ ε_{c})$ ⁠, Eq. (8), the effect of the temperature parameter $Δ T_{0}$ is taken into account, because the elastic strain range is excluded in $Δ ε_{in}$ ⁠. Other temperature parameters T_m and T_min are parameters relevant to both plastic and creep strains. The time parameter t_on is a parameter relevant to creep strain. In several studies on die attach cracking [27–29] described in Sec. 3 of this paper, the thermal cycling frequency f is introduced into the failure model based on $Δ T$ as a time parameter to take account of the effect of creep strain on the lifetime of power modules. These time parameters have an impact on creep strain. If we perform thermoplastic-creep analysis using a temperature distribution obtained from electrothermal analysis, the effects of these time factors are included in $Δ ε_{in} .$ The effects of other parameters I, V, and D are also considered in the same way as the time factors. Consequently, the effects of all parameters included in the failure models based on $Δ T$ are represented by $Δ ε_{in}$ ⁠, so that no additional parameter is required, when the failure model based on $Δ ε_{in}$ ⁠, Eq. (8), is used for estimating the lifetime of power modules. On the other hand, it must be considered what type of the failure model shown in Ref. [22] is selected for lifetime estimation, when the failure model based on $Δ T$ ⁠, Eq. (6), is used. It is concluded that the failure model based on $Δ ε_{in}$ is superior to those based on $Δ T$ ⁠, because the former being simpler than the latter, the number of constants to be determined by experimental data is less in the former model than in the latter ones. As a result, the former require a smaller number of experimental data than the latter to formulate failure models. Furthermore, as pointed out in Ref. [8], the values of the constants relevant to the parameters in the failure models based on $Δ T$ ⁠, Eqs. (9)–(14), are determined for a specific power module, and they would need to be recalibrated for another power module type. On the other hand, such recalibration is not required for the failure model based on $Δ ε_{in}$ ⁠, Eq. (8), regardless of power module dimensions if power modules are composed of the same materials.

Wire-liftoff is regarded as a crack propagation phenomenon due to thermal fatigue, and the strain intensity factor is applied to a wire-liftoff phenomenon [30]. The crack propagation rate per cycle is given by the following Paris law:

\frac{d a}{d N} = C {(Δ K_{ε eq})}^{n}

(15)

where a, N, and $Δ K_{ε eq}$ denote the crack length, the number of cycle, and the equivalent strain intensity factor range considering both the mode I (opening mode) and the mode II (in-plane shearing mode), respectively, and C and n are material constants. The remaining length of a jointed section is calculated based on Eq. (15) after constant temperature cycles, and compared with experimental data. Both results are relatively in good agreement. Although this methodology is reasonable from strength of materials, this cannot be applied to the failure at elevated temperature, in which plastic and creep deformations are significant.

2.2.2 Failure Models at Elevated Temperatures.

Power devices using wide bandgap semiconductors such as SiC and GaN are expected to operate over 200 $° C$ ⁠. Matsunaga and Uegai [31] report that aluminum wire-liftoff lifetime becomes saturated in the wire–bond junction temperatures more than $200 ° C$ as shown in Fig. 2. The similar phenomenon is observed by Yamada et al. [32] and Agyakwa et al. [33]. In the lifetime estimation based on Eq. (6), much more additional temperature parameters would be necessary to represent the lifetime saturation phenomenon. The inelastic strain range $Δ ε_{in}$ increases monotonically with increase in temperature because of increase in creep strain and increase in plastic strain due to the reduction of the yield stress of an aluminum wire. According to the lifetime estimation based on Eq. (8), the lifetime of aluminum wire-liftoff decreases monotonically with increase in temperature, and Eq. (8) does not explain the saturation phenomenon of the aluminum wire-liftoff lifetime shown in Fig. 2. Yang et al. [34] perform a study to explain this phenomenon. They propose a damage model of the interface between an aluminum wire and a silicon die which includes the effect of temperature and time-dependent material properties. The proposed model takes account of the damage removal due to recrystallization at high temperature as well as the damage build-up. This model has five parameters. The values of the constants relevant to the parameters are determined so as to represent experimental results of a particular system [33]. The time variations of interface damage and crack length can be obtained by solving simultaneous ordinary differential equations. It is not certain that the values of the constants determined so as to represent a particular system has generality, that is, they are applicable to a system with the same material as a particular system but different shape and dimensions.

Fig. 2

View large Download slide

Relationship between normalized number of cycles N_fand maximum temperature T_max (experimental results are adapted from Ref. [31])

Hereafter, much attention is paid to the failure models of wire-liftoff applicable to the wire-bond junction temperature over $200 ° C$ ⁠, in which the aluminum wire-liftoff lifetime becomes saturated. Any failure models described above are not suitable for representing such a saturation phenomenon. Above $200 ° C$ large stress relaxation occurs in an aluminum wire due to the large creep deformation and the large plastic deformation resulted from drastic reduction of the yield stress. The impact of stress relaxation on the wire-liftoff lifetime should be taken into account in the failure models.

2.2.2.1 Inelastic strain energy density range $Δ W_{in}$ ⁠.

The inelastic strain energy density range per cycle

Δ W_{in}

is the area of a cyclic stress–inelastic strain hysteresis loop, and denoted by the following integral:

Δ W_{in} = \oint_{C}^{} σ_{i j} d ε_{i j}^{i n}

(16)

where

σ_{i j}

and

ε_{i j}^{i n}

are the stress tensor and the inelastic strain tensor, respectively, and C denotes an integral path corresponding to a hysteresis loop of a stress–inelastic strain curve. The failure model is given by

N_{f} = C_{4} {(Δ W_{in})}^{p_{4}}

(17)

where $p_{4}$ and $C_{4}$ are material constants.

Equation (17) is proposed by Morrow [35] as an energy model for low cycle fatigue. This model contains not only inelastic strain but also stress, so that $Δ W_{in}$ is a physical quantity taking account of stress relaxation. Kanda et al. [36] show that Eq. (17) is equivalent to Eq. (8) in the case where a steady-state stress–inelastic stain loop is obtained under a cyclic loading condition. The energy model, Eq. (17), is applied to the lifetime prediction of solder joints in electronics packaging [29,37–39], and also to a wire-liftoff lifetime [40].

2.2.2.2 Nonlinear fracture mechanics parameter T* range $Δ T$ *.

Other possible physical quantities useful to the lifetime prediction of thermal fatigue at elevated temperatures are fracture mechanics parameters. They fall into two categories. One is a linear fracture mechanics parameter such as the stress intensity factor K. The other is a nonlinear fracture mechanics parameter such as the J-integral. At elevated temperature, large inelastic strain occurs due to the reduction of yield stress and creep deformation, so that nonlinear fracture mechanics parameters should be utilized to discuss the fracture at elevated temperatures. The J-integral proposed by Rice [41] is given by

J = \int_{Γ_{0}}^{} [W n_{1} - t_{i} u_{i, 1}] d Γ

(18)

where W, n_i, and t_i denote the stain energy density, the outward normal vector and the traction vector, respectively, and W and t_i are defined by

W = \int_{0}^{ε_{i j}} σ_{i j} d ε_{i j}

(19)

t_{i} = σ_{i j} n_{j}

(20)

$Γ_{0}$ is an arbitrary integral path enclosing the crack tip, as shown in Fig. 3. The physical meaning of the J-integral is the energy flux flowing into a fracture process zone in the vicinity of a crack tip. Under the condition of the deformation theory of plasticity, the J-integral has path-independence which enables calculating the accurate value of the J-integral by using far-field values of $σ_{i j}$ and $ε_{i j}$ around a crack tip. This feature is important because a crack tip is a singular point where $σ_{i j}$ and $ε_{i j}$ diverge.

Fig. 3

View large Download slide

J-integral path and T*-integral path

Afterward the J-integral is extended to the J-integral range

Δ J

for the application of fatigue crack growth under cyclic loading [42–45].

Δ J

is defined as follows:

Δ J = \int_{Γ_{0}}^{} [Δ W n_{1} - {Δ t}_{i} {Δ u}_{i, 1}] d Γ

(21)

where

Δ W = \int_{0}^{Δ ε_{i j}} Δ σ_{i j} d Δ ε_{i j}

(22)

In Eq. ,

Δ σ_{i j}

and

Δ ε_{i j}

are the stress range and the strain range, respectively, and defined as the difference between a current state and a reference state on the loading process in cyclic loading

Δ σ_{i j} = σ_{i j} - σ_{i j}^{(0)}

(23)

Δ ε_{i j} = ε_{i j} - ε_{i j}^{(0)}

(24)

The superscript (0) denotes a reference state. Similarly, the ranges of traction

{Δ t}_{i}

and displacement

Δ u_{i}

are defined by

Δ t_{i} = t_{i} - t_{i}^{(0)}

(25)

Δ u_{i} = u_{i} - u_{i}^{(0)}

(26)

According to Lamba [42], Wüthrich [44], and Tanaka [45], the J-integral range $Δ J$ has the same path-independence as the J-integral, if a crack is open and the strain energy density range $Δ W$ is given by a single valued function of $Δ ε_{i j}$ ⁠. Kubo et al. [46] has a detailed discussion on the path-independence of $Δ J$ when crack closure occurs during a load cycle.

The J-integral defined by Eq. loses path-independence for thermal stress problems, creep deformation problems, flow theory of plasticity problems, and dynamic problems. For such problems, Atluri et al. [47,48] propose a path-independent nonlinear fracture mechanics parameter T*-integral defined by

\begin{matrix} T * = \int_{Γ_{ε}}^{} [W n_{1} - t_{i} u_{i, 1}] d Γ \\ = \int_{Γ_{0}}^{} [W n_{1} - t_{i} u_{i, 1}] d Γ - \int_{A_{0} - A_{ε}}^{} [W_{, 1} - {(σ_{i j} u_{i, 1})}_{, j}] d A \end{matrix}

(27)

where the integral paths $Γ_{ε}$ and $Γ_{0}$ ⁠, and the domains $A_{0}$ and $A_{ε}$ are shown in Fig. 3. Here, $Γ_{0}$ is an arbitrary integral path enclosing a crack tip, and $Γ_{ε}$ is an integral path enclosing a fracture process zone existing near a crack tip. For a problem subjected to the deformation theory of plasticity, in which no unloading occurs, the T*-integral is identical to the J-integral, because the second term of the right-hand side in Eq. (27) becomes equal to zero. Thus, the T*-integral is generalization of the J-integral for various loading cases.

For thermal fatigue problems under repeated temperature cycles, in which both thermoelastic–plastic strain and creep strain occur, the authors of Ref. [40] newly propose the T*-integral range

Δ

T* defined in the same way as

Δ

J

\begin{matrix} Δ T^{*} = \int_{Γ_{ε}}^{} [Δ W n_{1} - {Δ t}_{i} {Δ u}_{i, 1}] d Γ \\ = \int_{Γ_{0}}^{} [Δ W n_{1} - {Δ t}_{i} Δ u_{i, 1}] d Γ - \int_{A_{0} - A_{ε}}^{} [{Δ W}_{, 1} - {(Δ σ_{i j} Δ u_{i, 1})}_{, j}] d A \end{matrix}

(28)

This expression is valid in the case where a crack is open. When a crack is closed,

Δ

T* in Eq. is modified by adding the correction term expressing a crack closure effect. For more details refer to Ref. [40]. The failure model based on

Δ

T* are given by

N_{f} = C_{5} {(Δ T^{*})}^{p_{5}}

(29)

2.2.2.3 Application of failure models at elevated temperatures.

The results of application of the failure models at elevated temperature, Eqs. (17) and (29), are shown in Ref. [40], where the wire-liftoff lifetimes predicted by these models are compared with the experimental data of Matsunaga and Uegai [31]. Figure 4 shows a finite element model of the wire bonding structure. The part of wire bonding is subjected to temperature cycles with the initial temperature of 25 $° C$ ⁠, the period of 2 s and the fixed temperature range $Δ T$ of $80 ° C$ ⁠. Finite element thermoelastic–plastic-creep analyses were performed for five cases of the maximum temperature T_max, that is, T_max = 105, 150, 200, 250, and 300 $° C$ ⁠.

Fig. 4

View large Download slide

Finite element model of wire bonding structure and boundary conditions

In the practical application of

Δ W_{in}

given by Eq. to the finite element analysis, we employ the area averaged inelastic strain energy density range

\bar{Δ W_{in}}

defined by

\bar{Δ W_{in}} = \sum_{k = 1}^{N} {{(Δ W_{in})}_{k} \times A_{k}} / \sum_{k = 1}^{N} A_{k}

(30)

where N and $A_{k}$ indicate the number of finite elements contained in the selected domain and the area of element k, respectively. In the lifetime calculation, $Δ W_{in}$ in Eq. (17) should be replaced by $\bar{Δ W_{in}}$ ⁠.

The exact calculation of $Δ$ T* can be done by using Eq. (28). There is no function of calculating the exact $Δ$ T* values in a commercial finite element program such as marc, ansys, etc. Thus, the following approximate calculation method for evaluating $Δ$ T* is employed to apply the commercial finite element program MSC.MARC2015 to the study on a wire-liftoff problem:

Smelser and Gurtin [49] show that the J-integral defined by Eq. (18) is valid for a bimaterial body if a crack is located along a bond line. Thus, the J-integral values can be calculated for several integral paths using a standard function of the general purpose finite element program MSC MARC2015.
The near-field J-integral value at the integral path $Γ_{ε}$ in the vicinity of a crack tip is obtained by using the extrapolation of the far-field J-integral values calculated in the process (1). As shown in Eq. (27), the near-field J-integral value is equivalent to the T*-integral value.
The near-field $Δ J$ -value is calculated from
$Δ J = J_{max} - J_{min}$
(31)
where $J_{max}$ and $J_{min}$ are the maximum and the minimum of near-field J-integral values during per loading cycle.
According to Eq. (28), the near-field $Δ J$ -value is regarded as $Δ$ T*-value.

Theoretically the J-integral range $Δ J$ defined by Eq. (31) is incorrect [50] but such a J-integral range is used in a lot of papers [51–56]^, because fatigue crack propagation rate is well correlated by it. Moreover, Ref. [57] shows that the $Δ J$ -value calculated from Eq. (31) agrees reasonably well with that calculated from the exact definition of $Δ J$ ⁠, Eq. (21), for a thermal fatigue problem.

The area average of the inelastic strain energy density range $\bar{Δ W_{in}}$ defined by Eq. (30) are calculated from the finite elements located on the extended line of a crack from the crack tip X = 0 to X_ave. Figure 5 shows the variation of $\bar{Δ W_{in}}$ with distance from crack tip. It is clearly found that $\bar{Δ W_{in}}$ depends on the domain selected for the calculation of $\bar{Δ W_{in}}$ ⁠. Thus the domain selection cannot be uniquely determined. The similar situation occurs when the inelastic strain range $Δ ε_{in}$ is utilized in the lifetime estimation.

Fig. 5

Variation of area average of strain energy density range ΔWin¯ with distance from crack tip for Tmax=200 °C and ΔT=80 °C

View large Download slide

Variation of area average of strain energy density range $\bar{Δ W_{in}}$ with distance from crack tip for $T_{max} = 200 ° C$ and $Δ T = 80 ° C$

The variations of the normalized number of cycle to failure with the maximum temperature $T_{max}$ are shown in Fig. 6, where the results obtained from failure models based on $Δ ε_{in}$ ⁠, $\bar{Δ W_{in}}$ ⁠, and $Δ$ T* are compared with the experimental data Matsunaga and Uegai [31]. In comparison with the experimental data, the failure model based on $Δ$ T* provides the best result.

Fig. 6

Relationship between normalized number of cycles Nf and maximum temperature Tmax. The predictions obtained from failure models based on ΔT*, ΔWin¯, and Δεin are compared with experimental data.

View large Download slide

Relationship between normalized number of cycles N_f and maximum temperature T_max. The predictions obtained from failure models based on $Δ$ T*, $\bar{Δ W_{in}}$ ⁠, and $Δ ε_{in}$ are compared with experimental data.

3 Die Attach Cracking

A power device chip is joined to a substrate by die attach. A die attach layer is subjected to repeated thermal stress due to the CTE mismatch between a semiconductor chip and a substrate, when temperature varies repeatedly. Such repeated thermal stress induces the thermal fatigue of a die attach layer, which results in die attach cracking. Although several review papers were published on the methodologies for estimating the thermal fatigue lifetime of solder joints in electronics packaging and the constitutive equations of solder materials [58–61], we pay attention to the papers on the lifetime estimation of die attach cracking in power modules, and review them critically.

Almost all the previous papers on die attach cracking deal with crack propagating in the horizontal direction of a die attach layer, and research on such die attach cracking is reviewed in Secs. 3.1 and 3.2. The die attach cracking propagating in the thickness or vertical direction of a die attach layer has been recently observed in power modules operated under elevated temperatures. Such a type of die attach cracking is also treated in Sec. 3.3.

Total lifetime of die attach cracking is generally the sum of crack initiation lifetime and crack propagation lifetime. Crack propagation lifetime can be neglected compared with total lifetime in the cases of wire-liftoff and the fatigue failure of solder bumps in flip-chip packaging, because the dimension of a wire jointed on a chip surface and that of solder bumps in flip-chip packaging are at most a few millimeters. On the other hand, the dimension of a jointed part is a few centimeters in the case where the surface of a power device chip fully adheres to a substrate. Thus, the ratio of crack propagation lifetime to total lifetime becomes large. Here, we describe the lifetime estimation methodologies of die attach cracking without and with consideration of crack propagation.

3.1 Lifetime Estimation Methodologies Without Consideration of Crack Growth

3.1.1 Lifetime Estimation Formulae Using Temperature Range.

Inspired by the paper of Norris and Landzberg [27], the following estimation formulae using the temperature range $Δ T$ are proposed:

Refs. [27] and [62]:
$N_{f} = A f^{m} {Δ T}^{- n} e xp (E_{a} / k_{b} T_{max})$
(32)
Refs. [62] and [63]:
$N_{f} = A f^{m} {Δ T}^{- n} e xp (E_{a} / k_{b} T_{max}) {[corr (Δ T)]}^{- 1 / c}$
(33)
Ref. [64]:
$N_{f} = A {(1 / t_{hot})}^{m} {Δ T}^{- n} e xp (E_{a} / k_{b} T_{max})$
(34)

where A, m, n, and

E_{a}

are the material constants determined from experimental data,

k_{b}

denotes the Boltzmann constant, and

T_{max}

is the maximum absolute temperature. Equation is proposed for Pb–Sn eutectic solder alloys. Equation is extended to the application of Sn–Ag–Cu lead-free solder alloys, and Eqs. and are such modified versions of Eq. . In Eqs. and , f denotes the frequency of temperature cycle, and

[Corr (Δ T)]

is a correction term given by

corr (Δ T) = a \ln (Δ T) + b

(35)

where a and b are package and solder material specific constants provided in Ref. [65], and c denotes a constant depending on temperature profile [65]. The parameter $t_{hot}$ in Eq. (34) is the time per cycle, during which the component is hot. As described in Sec. 2.2, additional parameters are added to the lifetime estimation formulae based on $Δ T$ to express complex phenomena.

We can define the acceleration factor AF as follows:

A F = \frac{{(N_{f})}_{2}}{{(N_{f})}_{1}}

(36)

where the subscripts 1 and 2 indicate separate states. Based on this equation, the lifetime in the service condition of a power module

{(N_{f})}_{2}

can be estimated from that of an acceleration test

{(N_{f})}_{1}

⁠. The expression of AF for Eq. is given by

A F = \frac{{(N_{f})}_{2}}{{(N_{f})}_{1}} = {(\frac{f_{2}}{f_{1}})}^{m} {(\frac{{Δ T}_{1}}{{Δ T}_{2}})}^{n} e xp [(E_{a} / k_{b}) {1 / {(T_{max})}_{2} - {1 / (T_{max})}_{1}}]

(37)

In deriving the above equation, we utilize the fact that the values of A, m, n, and $E_{a}$ are identical both in the state 1 and in the state 2. The formulae similar to Eq. (37) can be obtained for Eqs. (33) and (34).

3.1.2 Lifetime Estimation Formula Using Inelastic Strain Range.

The lifetime predictions can be done based on the Coffin–Manson law, Eq. , that is,

N_{f} = C_{3} {(Δ ε_{in})}^{p_{3}}

(38)

Being different from the temperature range $Δ T$ ⁠, the inelastic strain range $Δ ε_{in}$ is difficult to measure, but its accurate value can be obtained from inelastic finite element analysis.

The lifetime estimation formulae based on $Δ T$ ⁠, Eqs. (32)–(34), include the maximum temperature $T_{max}$ ⁠, the temperature cycle frequency f, and the hot dwell time per temperature cycle $t_{hot}$ in addition to $Δ T$ to accommodate these formulae to complex operating conditions. The constants relevant to these parameters, A, m, n, and $E_{a}$ ⁠, are determined so as to fit experimental data. All the additional parameters $T_{max}$ ⁠, f, and $t_{hot}$ are related to plastic and creep strains, so that the inelastic strain range $Δ ε_{in}$ obtained from inelastic finite element analysis automatically includes the effects of $T_{max}$ ⁠, f, and $t_{hot}$ ⁠. Consequently, no additional parameter is needed in the failure model based on the inelastic strain range $Δ ε_{in}$ even under complex operating conditions.

The lifetime estimation formula based on $Δ ε_{in}$ ⁠, Eq. (38), has been widely used for solder joints in IC packaging since 1990s [23,66–70], but the applications of the failure model based on $Δ ε_{in}$ ⁠, Eq. (38), to the lifetime of die attach cracking in power modules [71–74] are fewer than those of the failure models based on $Δ T$ ⁠, Eqs. (32)–(34).

3.1.3 Lifetime Estimation Formula Using Inelastic Strain Energy Density Range.

Although the estimation formula using the inelastic strain energy density range is shown in the wire-liftoff lifetime, this formula can be also applied to predicting the die attach cracking lifetime. The formula is given again as follows:

N_{f} = C_{4} {(Δ W_{in})}^{p_{4}}

(39)

The above formula has been applied to the lifetime predictions of small bumps in IC packages [37,68,70,75–78], in which crack propagation lifetime can be neglected in comparison with total lifetime. In Refs. [29], [39], [79], and [80], the formula, Eq. (39), is applied to the lifetime prediction of die attach cracking in power modules. The size of a die attach layer is a few centimeters and much larger than the wire bonding size, so that the estimation formula for die attach cracking lifetime, Eq. (39), is usually used together with the formula for crack propagation lifetime as described later.

In the practical application of Eq. ,

Δ W_{in}

is replaced by the area/volume-averaged inelastic strain energy density range

\bar{Δ W_{in}}

given by Eq. . In such a case, we must consider what domain is selected as the inelastic strain energy density range concerned with thermal fatigue failure. Fatigue failure generally occurs at a stress concentration part. Thus, we should select a domain including a stress concentration part, but we cannot determine such a domain uniquely. Nakajima et al. [81] deal with a crack in a homogeneous material, assuming that stress and strain around a crack tip are governed by the HRR-fields [82,83]. They consider a square domain with the side length of L_area, the center of which is located at the crack tip. Let

Δ W_{in}

be the inelastic strain energy density range averaged in the square domain with L_area. The following relationship holds, because

{Δ W}_{in}

and L_area are inversely proportional to each other

{Δ W}_{in - c} = Δ W_{in} {\cdot L}_{area}

(40)

where $Δ W_{in - c}$ is a proportional constant independent of $L_{area}$ ⁠. It is therefore possible to discuss die attach cracking lifetime using the new parameter ${Δ W}_{in - c}$ independent of the domain selected for the calculation of $Δ W_{in}$ ⁠. Equation (40) is true in the cases of plastic deformation under the deformation theory of plasticity and steady-state creep deformation. Otherwise ${Δ W}_{in - c}$ is domain-dependent, for example, in the cases of plastic deformation with unloading and transient creep deformation where the stress and strain around a crack tip is not governed by the HRR-fields. Also the stress singularity at a notch and a corner is not represented by HRR-fields, so that ${Δ W}_{in - c}$ is domain-dependent.

3.2 Lifetime Estimation Methodologies With Consideration of Crack Propagation.

Fatigue crack propagation rate, that is, crack propagation per cycle da/dN is written by the Paris law [84]

\frac{d a}{d N} = C_{1} {(Δ K)}^{n_{1}}

(41)

where

C_{1}

and

n_{1}

are material constants, and

Δ K

denotes the stress intensity factor range per cycle. Equation is valid under the condition that small-scale yielding is satisfied around a crack tip. A crack in a ductile material does not satisfy this condition. In such a case, the original Paris law Eq. should be modified using the inelastic strain range

Δ ε_{in}

defined by Eq. , the inelastic strain energy density range

Δ W_{in}

defined by Eq. , and the J-integral range

Δ

J defined by Eq. . They are given by

\frac{d a}{d N} = C_{2} {(Δ ε_{in})}^{n_{2}}

(42)

\frac{d a}{d N} = C_{3} {(Δ W_{in})}^{n_{3}}

(43)

\frac{d a}{d N} = C_{4} {(Δ J)}^{n_{4}}

(44)

As

Δ

J loses path-independence in the cases of plastic deformation with unloading and creep deformation, the T*-integral range

Δ

T* should be used instead of

Δ

J

\frac{d a}{d N} = C_{5} {(Δ T^{*})}^{n_{5}}

(45)

The lifetime estimation formulae shown in Sec. 3.1 basically provide fatigue crack initiation lifetime and good predictions of the fatigue lifetime for small-sized solder joints such as microbumps. The die attach size of power modules is so large that crack propagation lifetime cannot be neglected. Let

N_{0}

and

N_{1}

be fatigue crack initiation lifetime and crack propagation lifetime, respectively, then total lifetime N is given by

N = N_{0} + N_{1} = N_{0} + \frac{L_{c}}{(d a / d N)}

(46)

where $L_{c}$ denotes the critical crack length, at which a rapid increase of temperature or a rapid variation of electric quantities occurs.

3.2.1 Lifetime Estimation Methodology Using Inelastic Strain Range.

In this methodology the formulae Eqs. and are, respectively, used for crack initiation lifetime and for crack propagation lifetime.

{Formula for crack initiation lifetime : N}_{0} = α_{1} {(Δ ε_{in})}^{β_{1}}

(47)

Formula for crack propagation rate : \frac{d a}{d N} = γ_{1} {(Δ ε_{in})}^{δ_{1}}

(48)

The point at issue for the application of the above formulae is as follows. As shown in Sec. 2.2.2.3, a fatigue crack arises at a stress concentration region. There is arbitrariness in selecting the evaluation point of $Δ ε_{in}$ ⁠, because a crack tip is a stress and strain singular point after crack initiation. In the original Paris law shown in Eq. (41), the stress intensity factor range $Δ K$ is utilized to avoid such arbitrariness. The stress intensity factor $K$ is one of the indicator representing the intensity of stress singularity. The $J$ -integral and T*-integral, respectively, given in Eqs. (44) and (45) are also the similar indicators to the stress intensity factor K, and determine the intensity of stress singularity uniquely, if they have path-independent nature.

Déplanque et al. [85] applied this methodology to the fatigue lifetime estimation of the die attach cracking for a Sn–Pb eutectic solder alloy and a Sn–Ag–Cu lead free solder alloy in a power-MOS transistor. They used the accumulated inelastic (creep) strain averaged from a crack tip (x = 0) to a certain point (x = u) along the extended line of a crack. Although the definition of accumulated inelastic (creep) strain is not clear in this paper, it is equivalent to the inelastic (creep) strain range $Δ ε_{in}$ used in Eqs. (47) and (48) if it represents the accumulated inelastic (creep) strain per temperature cycle. Other researchers also use Eq. (48) to obtain the crack propagation lifetime of a power module [86–88]. They divide a crack propagation path into multiple subsections, each of which is a side of the finite element allocated along a crack path, and the lifetime of crack propagation that makes advance of each subsection is obtained from Eq. (48). It is not clear whether crack initiation lifetime is considered in this methodology. If it is not considered, the present methodology provides shorter lifetime.

Japanese researchers propose a methodology for simulating not only fatigue crack initiation but also a crack propagation path [89,90]. In this methodology, a damage parameter D is obtained both from the Coffin–Manson law and the linear cumulative damage rule (Miner's rule) based on the finite element stress analysis taking account of both plastic and creep strains, and the elements satisfying $D = 1$ are removed by reducing the Young's modulus at a small value compared with an undamaged material. This methodology is utilized to predict the crack propagation path of QFP solder joints [89] and BGA solder joints [90] under thermal fatigue conditions, and also applied to the crack propagation in a die attach layers of a IGBT power module [91,92]. This methodology can deal with the propagation of a two-dimensional crack in a three-dimensional structure, although the methodology employed in Refs. [85–88] deals with only that of a one-dimensional crack in a two-dimensional structure.

3.2.2 Lifetime Estimation Methodology Using Inelastic Strain Energy Density Range.

In this methodology, the formulae Eqs. and are, respectively, used for crack initiation lifetime and for crack propagation lifetime

Formula of crack initiation lifetime : N_{0} = α_{2} {(Δ W_{in})}^{β_{2}}

(49)

Formula of crack propagation rate : \frac{d a}{d N} = γ_{2} {(Δ W_{in})}^{δ_{2}}

(50)

This methodology is proposed by Darveaux for evaluating the lifetime of die attach cracking in IC packages [93,94], and is used by other authors not only for the estimation of die attach cracking lifetime of IC packages [95,96] but also for that of power modules [97–99]. As described in Sec. 2.2.2.3 and 3.1.3, in the practical application of $Δ W_{in}$ to Eqs. (49) and (50), Darveaux proposed the use of volume averaged $Δ W_{in}$ ⁠, ${(Δ W_{in})}_{ave}$ ⁠, which is basically equivalent to $\bar{Δ W_{in}}$ given by Eq. (30) described in 2.2.2.3. Then what region is chosen to calculate ${(Δ W_{in})}_{ave}$ ? Any answers have never been provided for this question by the authors using the inelastic strain energy density range, except that the region where the highest stress is expected should be always included.

Kariya et al. [100] obtain the constants

γ_{2}

and

δ_{2}

in a formula of crack growth rate Eq. for sintered Ag nanoparticles which is expected as a die attachment material of next-generation power modules. In Ref. [100], they use the volume-averaged

Δ W_{in}

⁠, the value of which depends on the domain selected for the calculation. So they propose the use of

{Δ W}_{in - c}

defined by Eq. instead of volume-averaged

Δ W_{in}

⁠, because

{Δ W}_{in - c}

is independent of the domain selected for the calculation of

Δ W_{in}

and uniquely determined [101]. In this case, fatigue crack propagation rate is given by

\frac{d a}{d N} = γ_{2} {(Δ W_{in - c})}^{δ_{2}}

(51)

They perform fatigue crack propagation tests under a displacement-controlled pulsating tensile–tensile mode using single-edge notched specimens made of sintered Ag nanoparticles and determine the constants $γ_{2}$ and $δ_{2}$ in Eq. (51) from experimental data. As described in Sec. 3.1.3, the use of the formula Eq. (51) is limited to plastic deformation without unloading obeying the deformation theory of plasticity and steady-state creep deformation.

3.2.3 Lifetime Estimation Methodology Using Fracture Mechanics Parameters.

The formulae of crack propagation rate using fracture mechanics parameters are Eqs. (41), (44), and (45). The inelastic strain range $Δ ε_{in}$ and the inelastic strain energy density range $Δ W_{in}$ depend on the point or domain selected for their estimations and they are not determined uniquely. On the other hand, the stress and strain around a crack tip are uniquely determined by fracture mechanics parameters, and $Δ$ J and $Δ$ T* defined by path integral have an advantage of numerical analysis because the values of $Δ$ J and $Δ$ T* can be evaluated from the far-field mechanical quantities if they have path-independence. The parameter $Δ W_{in - c}$ proposed by Kariya et al. can be also regarded as one of the fracture mechanics parameters.

Thermal fatigue lifetime analysis of a power module is performed using the original Paris law, Eq. (41), as a formula of crack propagation lifetime together with the Coffin–Manson law, Eq. (38), as a formula of crack initiation lifetime [102]. According to the results of the analysis, almost all the total lifetime is occupied by the crack initiation lifetime, and the crack propagation lifetime is only 3.6% of the total lifetime. The die attachment material considered in the analysis is a 95Pb–5Sn solder alloy. The use of the original Paris law written by the stress intensity factor range $Δ$ K is not valid because a 95Pb–5Sn solder alloy is ductile and the small-scale yielding condition is not satisfied, and provide shorter crack propagation life time.

Pao performs pioneering research on the thermal fatigue life prediction of solder joints in electronic packaging [103]. He uses the

Δ

J-based equation, Eq. , as a formula of crack propagation rate due to thermal fatigue and the C*-based equation [104] as a formula of crack propagation rate due to creep, that is

Formula of crack propagation rate due to thermal fatigue : \frac{d a}{d N} = γ_{3} {(Δ J)}^{δ_{3}}

(52)

Formula of crack propagation rate due to creep : \frac{d a}{d t} = γ_{3^{'}} {(C^{*})}^{δ_{3^{'}}}

(53)

where t denotes time, and the total lifetime $N_{f}$ is determined from the number of cycles, at which a crack length reaches the critical one $L_{c}$ ⁠. Although crack initiation life time is not considered in the analysis, the lifetime obtained from the analysis agrees well with experimental data. The reason for good agreement would be due to the fact that the critical crack length $L_{c}$ of one half of solder joint width, a large crack length, is assumed in the analysis.

Kariya and coworkers [105] carry out fatigue crack propagation tests of single-edge notch specimens made of sintered Ag nanoparticles under displacement-controlled tensile–tensile loading condition and the isothermal temperature condition of 140

° C

⁠. They employ the

Δ

J-based relationship as a formula of fatigue crack propagation rate:

\frac{d a}{d N} = γ_{3} {(Δ J)}^{δ_{3}}

(54)

and the constants $γ_{3}$ and $δ_{3}$ are determined from experimental data. In applying Eq. (54), $Δ J$ value is calculated from a simplified evaluation formula proposed by Asada et al. [106]. If thermal and creep deformations are taken into account, $Δ$ J loses path-independence and cannot be uniquely determined like $Δ ε_{in}$ and $Δ W_{in}$ ⁠.

The fracture mechanics parameters T* and

Δ

T* described in detail in Sec. 2.2.2.2 are path-independent and determined uniquely under thermal and creep deformations; thus valid fracture mechanics parameters for such deformations. The present authors of this paper recommend that the following formula be used to represent fatigue crack propagation rate under elevated temperature condition, in which thermal and creep deformations are significant in addition to plastic deformation:

\frac{d a}{d N} = γ_{4} {(Δ T^{*})}^{δ_{4}}

(55)

A crack must be assumed in the methodologies for the lifetime estimation using fracture mechanics parameters such as $Δ T^{*}$ ⁠. Thus, an important problem in applying a fracture mechanics parameter to lifetime estimation is how to estimate fatigue crack initiation lifetime. The following are solutions of this problem:

Fatigue crack initiation lifetime is estimated by the methodologies without assuming a crack, that is, Eqs. (47) and (49).
Fatigue crack initiation lifetime $N_{0}$ is obtained from the following formula:
$N_{0} = α_{4} {(Δ T^{*})}^{β_{4}}$
(56)

by assuming a short-length crack. In this case, a problem what length of a crack should be regarded as a short-length crack is left unresolved. One of the method to resolve this problem is to carry out sensitivity analysis of $N_{0}$ by varying a crack length. We regard a crack with the length, below which the sensitivity of crack length to $N_{0}$ becomes low, as a short-length crack.

3.2.4 Other Lifetime Estimation Methodologies

3.2.4.1 Continuum damage mechanics.

The continuum damage mechanics deals with the failure process from the generation and coalescence of microdefects such as voids to macrocrack propagation within the framework of the continuum mechanics by introducing a damage variable D ranging from zero to one. The values of D equal to zero and one represent no damage and full damage, respectively. D is usually defined as the ratio of the damaged area to the initial area. The damage variable D has a direct influence on stress and an indirect influence on constitutive equations of materials via stress. A detailed description of the continuum damage mechanics is given in the monograph published by Murakami [107]. In this approach, it is possible to deal with a failure problem of a structure without assuming a crack explicitly.

The applications of the continuum damage mechanics to the failure of solder joints are presented in Refs. [108] and [109]. In these papers, the damage levels of solder joints are obtained for simple test systems and loading conditions. The damage mechanics is applied to the die attach cracking of power modules [110,111]. These studies use the approach similar to that employed by Yang et al. in the estimation of wire-liftoff lifetime [34]. In accordance with Yang et al.'s approach, the factors impacting on the variation of D are expressed by the functions of the temperature T, the inelastic strain $ε_{in}$ ⁠, the damage variable D, and the distance from a crack tip x, and the constants included in each function are determined from the experimental data on the relationship between the crack length and the number of thermal cycles shown in Ref. [112]. The time variations of damage variable and crack length can be obtained by solving simultaneous ordinary differential equations. In the continuum damage mechanics, the damage variable D is incorporated into the framework of the continuum mechanics and included in the governing equations of the continuum mechanics. On the other hand, in the approach used in Refs. [110] and [111], the selections of the functions for representing the damage variable D are empirical rather than theoretical. So it would not be suitable to classify this approach as the continuum damage mechanics.

Desai et al. [113–116] apply another type of the continuum damage mechanics approach. Their approach is named the distributed state concept (DSC) approach based on the unified constitutive modeling and takes account of elastic, plastic and creep strains, microcracking, damage and degradation or softening, and postpeak stiffening under thermomechanical loading for materials. The basic assumption of the DSC approach is as follows: the deforming material element is the mixture of material parts in two reference state, relative intact (RI) and fully adjusted (FA) states. Under this assumption, the observed stress tensor,

σ_{i j}^{a}

⁠, is given by

σ_{i j}^{a} = (1 - D) σ_{i j}^{i} + D σ_{i j}^{c}

(57)

where the superscripts a, i, and c denote the observed, RI and FA states, respectively, and D is a scalar disturbance function or the ratio of the change in disorder parameter to the original reference state, which is regarded as degradation/damage metric. In Refs. [113–116], a degradation function D is empirically obtained by curve fitting to test data. This approach is applied to the thermal fatigue problems of Sn–Pb solder joints of electronic packaging. A good review on the DSC approach based on the unified constitutive modeling is given by Desai [117]. Furthermore Basaran and Yan [118] and Basaran and Tang [119] propose the unified mechanics theory, that is, the unification of Newton's universal laws of motion and laws of thermodynamics, and derive the expression of D from the entropy, a measure of disorder in the system. We can obtain the expression of D theoretically without test data in the unified mechanics theory. The DSC approach based on the unified mechanics theory is applied to low cycle fatigue [118], thermal fatigue [119], dynamic loading with harmonic vibration [120], and concurrent thermal cycle and dynamic loading [121]. Sophisticated derivations of both analytical and computational three-dimensional models based on the unified mechanics theory have been recently provided for predicting damage and fatigue life by a research team led by Basaran and coworkers [122].

We can easily incorporate the DSC approach into nonlinear finite element formulation [115,120], and obtain the variation of D in the process of fatigue using a nonlinear finite element analysis. Although the DSC approaches based on the unified constitutive modeling and the unified mechanics theory have not yet been utilized to estimate the lifetime of die attach cracking in power modules, they are the promising approaches. We are expecting that these approaches are utilized to structural integrity evaluation of power modules.

3.2.4.2 Cohesive zone model.

The cohesive zone model (CZM) is a model for dealing with a crack propagation phenomenon in ductile materials. In this model, cohesive force acts on crack surfaces even after the separation of crack surfaces occurs, and a crack gradually propagates under the resistance of the cohesive force. The region of crack surfaces on which cohesive force acts is called a cohesive zone. It is assumed that the crack surfaces completely separate when cohesive force becomes zero. Detailed description of the CZM is given in Ref. [123]. The model can deal with the failure of a structure without assuming a crack with a crack tip having stress and strain singularity. The CZM is applied not only to the thermal fatigue of microsolder bumps [124] but also that of a die attach layer in a power module [125].

3.3 Lifetime Estimation Methodologies for Die Attach Failure in the Thickness Direction.

Demands of the high temperature operation of power modules have been increasing in recent years. Die attach failure in the thickness direction of a die attach layer (hereafter abbreviated as “vertical failure”) has been observed in such an environment [126–128]. In Ref. [126], the vertical failure of a Sn–Cu die-attach solder joint is observed, when power cycling tests (PCTs) are performed under the temperature conditions of (⁠ $T_{j \min} = 50 ° C, T_{j \max} = 150 ° C)$ and (⁠ $T_{j \min} = 50 ° C, T_{j \max} = 175 ° C$ ⁠), where $T_{j \min}$ and $T_{j \max}$ denote the minimum junction temperature and the maximum junction temperature, respectively. In Ref. [127], three kinds of die-attach solder alloys, Pb–5Sn–1.5Ag, Sn–7Cu, and Sn–3Cu–10Sb, are tested under the power cycling condition of $T_{j \max} = 175 ° C$ and $Δ T =$ 150 $° C$ ⁠. As a result, the dominant die attach failure mode of the Pb-based solder alloy, Pb–5Sn–1.5Ag, is the failure of the horizontal direction in a die attach layer (hereafter abbreviated as “horizontal failure”), while that of the Sn-based solders alloys, Sn–7Cu and Sn–3Cu–10Sb, is vertical failure. In Ref. [128], the thermal cycling tests (TCTs) ranging from 50 $° C$ to 175 $° C$ are performed for four kinds of test specimens consisting of a Si-die, a Sn–Ag–Cu–Sb die attach layer and four kinds of substrates with different CTEs. The test results are summarized as follows:

The horizontal failure is mainly observed in the test specimens consisting of the substrate with large CTE, that is, in the test specimens with the large CTE mismatch between a Si-die and a substrate.
The vertical failure is mainly observed in the test specimens consisting of the substrate with small CTE, that is, in the test specimens with the small CTE mismatch between a Si-die and a substrate.

A study on the vertical failure mechanism of a Sn-based solder layer is also performed based on the microstructure images of solder layers [129]. According to this study, fatigue cracks are originated from the intermetallic compounds such as Cu₆Sn₅ around $β$ -Sn dendrite boundaries. They propagate along the dendrite boundaries/the grain boundaries, and then, the vertical direction in a solder layer under the normal biaxial-stress condition caused by cyclic temperature variation. The failure mechanism mentioned above is confirmed by a finite element simulation.

Now, let us show the lifetime estimation methodologies for vertical failure. Harubeppu et al. [130] propose a lifetime estimation formula, assuming that voids are generated and grown on the grain boundaries by creep deformation at high temperature, and further that the grain boundary voids are connected and form intergranular cracks by the tensile stress at low temperature. They define the creep damage per cycle $D_{c}$ as accumulated creep strain per cycle and the fatigue damage per cycle $D_{f}$ as a power of temperature range $Δ T$ ⁠. The creep damage is caused only by creep deformation, so that $D_{c}$ is independent of $D_{f}$ ⁠. On the other hand, the fatigue damage promotes the coalescence of grain boundary voids generated by creep damage, so that the contribution of fatigue damage to the total damage is small when $D_{c}$ is small and becomes large when $D_{c}$ increases. Considering these facts, they propose the total damage $D$ per cycle as $D = D_{c} (1 + D_{f}) .$ Then, the lifetime $N_{f}$ is given by $N_{f} \propto 1 / D$ ⁠. The proposed lifetime estimation formula contains three material properties characterizing creep strain rate, three material properties characterizing thermal fatigue, and one coefficient dependent on the structure of a power module, which is the so-called fitting parameter, determined from the lifetime experiments of a specific power module. Conclusively, the lifetime estimation formula proposed by Harubeppu et al. [130] cannot be used as the lifetime estimation in the design phase of a power module, because the lifetime of a power module is not estimated only by material properties.

Studies on the creep damage caused by repeated thermal stress due to temperature cycle have been performed in the structural reliability studies of the equipment of chemical plants and nuclear plants. The linear cumulative damage rule [131] and the strain range-partitioning method [132,133] are representative methodologies for dealing with fatigue-creep interaction problems that both fatigue damage and creep damage have to be simultaneously taken into consideration. In the former methodology, the summation of time fractions is used as a measure of creep damage proposed by Robinson [134], and that of cycle fractions is used as a measure of fatigue damage called Miner's rule. Total damage is defined as the summation of both damages. The continuum damage mechanics provides more reasonable creep damage than that proposed by Robinson [135]. In the latter methodology, we consider two strain modes, that is, plastic flow and creep, and consider four combinations of plasticity and creep depending on tensile and compression as follows:

$Δ ε_{p p}$ ⁠: tensile plastic strain reversed by compressive plastic strain.
$Δ ε_{c p}$ ⁠: tensile creep strain reversed by compressive plastic strain.
$Δ ε_{p c}$ ⁠: tensile plastic strain reversed by compressive creep strain.
$Δ ε_{c c}$ ⁠: tensile plastic creep reversed by compressive creep strain.

Then, we divide total inelastic strain range into the above four stain range, and obtain the fraction of each strain range

F_{i j}

(⁠

i j = p p, c p, p c c c)

⁠. Based on the partitioned strain range–lifetime relationships

Δ ε_{i j} - N_{f, i j} (i j = p p, c p, p c, c c)

obtained from the experiments, we can estimate the lifetime

N_{f}

as

\frac{F_{p p}}{N_{f, p p}} + \frac{F_{c p}}{N_{f, c p}} + \frac{F_{p c}}{N_{f, p c}} + \frac{F_{c c}}{N_{f, c c}} + \frac{1}{N_{f}}

(58)

The lifetime of the vertical failure in the die attach of power modules can be estimated using the above methodologies, but such a study is not yet performed.

4 Power Cycling Test and Thermal Cycling Test

Both a power cycling test (PCT) or a thermal cycling test (TCT) have been utilized to the structural integrity evaluation of power modules. A PCT is done by setting test conditions such as voltage, current, cooling method, and ambient temperature of a test specimen. Thus, a test specimen is subjected to temperature cycling due to the joule heating of a power module. In a TCT, a test specimen is subjected to temperature cycling in a temperature cycle chamber. A temperature distribution is heterogeneous for a PCT and homogeneous for a TCT. As shown in Refs. [3] and [26], wire-liftoff lifetimes are usually obtained from PCTs. In such cases, the junction temperature is evaluated through the collector–emitter voltage V_CE as described in IEC 60747 standard [136]. The junction temperature obtained from V_CE represents the average chip temperature and does not represent the temperature at the bonding part of a wire. A wire-liftoff phenomenon is a local failure of a bonding wire caused by the CTE mismatch of a bonding wire and a semiconductor die. The local temperature at the bonding part contributes to the wire-liftoff. In some cases, there exists much difference between such a local temperature and the average chip temperature obtained from V_CE, as shown in Ref. [137]. The temperature of the bonding part in a wire measured by infrared cameras or thermocouples, or the temperatures calculated using electrothermal simulation tools such as ansys [136] should be used instead of the temperatures obtained from V_CE to formulate the wire-liftoff failure models shown in Eqs. (9)–(14). When Eq. (8) is used for evaluating the wire-liftoff lifetime, $Δ ε_{in}$ is an indicator representing a local inelastic strain of the bonding part in a wire which directly influences a local failure mode called wire-liftoff. Other physical quantities such as the inelastic strain energy density range $Δ W_{in}$ ⁠, the nonlinear fracture mechanics parameter J-integral range $Δ$ J and the nonlinear fracture mechanics parameter T*-integral range $Δ$ T* are similar to $Δ ε_{in}$ ⁠. From a scientific view point, the temperature range $Δ T$ obtained from V_CE in a PCT is not a valid parameter for describing wire-liftoff lifetime. If we insist on using the $Δ T$ -based wire-liftoff failure models, we should perform TCTs instead of PCTs or should use the local temperature of a bonding part measured experimentally or calculated by CAE tools. Better failure models for predicting wire-liftoff lifetime are those based on $Δ ε_{in}$ ⁠, $Δ W_{in}$ ⁠, $Δ$ J, and $Δ$ T*.

Now, let us consider the horizontal failure of a die attach layer. All the lifetime estimation formulae given in Eqs. – are obtained from TCTs. A temperature distribution is homogeneous in a TCT, so that the temperature of a substrate is the same as that of a semiconductor die. The driving force of the horizontal failure of a die attach layer in a TCT is the shear stress caused by a global mismatch in the CTEs of a semiconductor die and a substrate,

Δ

CTE, and given by

Driving force \propto (α_{D} - α_{S}) (T_{max} - T_{min}) = Δ CTE \times Δ T

(59)

where

α_{D}

and

α_{S}

denote the CTEs of a semiconductor die and a substrate, respectively, and

T_{max}

and

T_{min}

denote the maximum temperature and the minimum temperature in a temperature cycle, respectively. In a PCT, there is a temperature gap between a semiconductor die and a substrate. Let

Δ_{max}

and

Δ_{min}

be the temperature gaps at the temperatures of

T_{max}

and

T_{min}

⁠, respectively. The driving force of the horizontal failure of a die attach layer in a PCT is given by

\begin{matrix} Driving force \propto (α_{D} - α_{S}) (T_{max} - T_{min}) + α_{S} (Δ_{max} - Δ_{min}) \\ = Δ CTE \times Δ T + α_{S} (Δ_{max} - Δ_{min}) \end{matrix}

(60)

The temperature gaps $Δ_{max}$ and $Δ_{min}$ are expected to be very small, because a die attach layer is very thin compared with a semiconductor die and a substrate. Thus, the second term can be neglected in Eq. (60), and the driving force of the horizontal failure of a die attach layer is almost the same both in a TCT and a PCT. It is concluded that the failure models obtained from TCTs, Eqs. (32)–(34), can be applied to the lifetime estimation of power modules under operational conditions.

5 Summaries of Previous Research

In the present paper, we review the methodologies for estimating the thermal fatigue lifetimes of bonding wires and die attach layers in power modules. The following are summaries of the present review:

The purely empirical formulae based on the temperature range $Δ T$ have been proposed to estimate a wire-liftoff lifetime and a die attach cracking lifetime mainly by European and American researchers. They are Eqs. (9)–(14) for wire-liftoff lifetime and Eqs. (32)–(34) for die attach cracking lifetime. In these formulae, the lifetime $N_{f}$ is expressed by the product of physical quantities influencing $N_{f}$ ⁠, and the constants relevant to the physical quantities are determined so as to fit experimental data, and a lot of experimental data are required to determine these constants for complex operating conditions of power modules. Furthermore, the use of these formula is limited to a specific power module, from which experimental data are obtained to formulate the equations for lifetime estimation.
The lifetime estimation formulae for wire-liftoff and die attach cracking based on the inelastic strain range $Δ ε_{in}$ ⁠, the inelastic strain energy density range $Δ W_{in}$ ⁠, the nonlinear fracture mechanics parameter J-integral range $Δ$ J and the nonlinear fracture mechanics parameter T*-integral range $Δ$ T* are simpler than those based on the temperature range $Δ T$ ⁠, because the above-mentioned quantities are physical quantities directly influencing thermal fatigue. This is because the lifetime estimation formula based on $Δ ε_{in}$ ⁠, for example, contains all the temperature factors and the time factors of that based on the temperature range $Δ T$ ⁠. Other physical quantities except for the temperature range $Δ T$ are similar to $Δ ε_{in}$ ⁠.
Large stress relaxation occurs due to large plastic deformation resulted from reduction of yield stress and creep deformation. Thus, the lifetime estimation formula based on $Δ ε_{in}$ cannot represent the phenomena peculiar to the high temperature operation of power modules. The impact of stress relaxation on the thermal fatigue lifetime should be taken into account in the formulae. The formulae based on $Δ W_{in}$ ⁠, $Δ$ J, and $Δ$ T* are effective in the case of the high temperature operation of power modules. Among them, $Δ W_{in}$ depends on the domain selected for their estimation and is not determined uniquely. As for the fracture mechanics parameters such as $Δ$ J and $Δ$ T*, their values are uniquely determined if they have path-independence. Under the condition of thermal and creep deformation, $Δ$ J is path-dependent, while $Δ$ T* is path-independent. Thus, the use of the formula based on $Δ$ T* is recommended to estimate the thermal fatigue lifetime at elevated temperature.
In wire-liftoff lifetime, crack propagation lifetime can be neglected, because a wire-bonding length is very short. On the other hand, the die attach size of power modules is so large compared with a wire-bonding length that crack propagation lifetime should be considered in the lifetime estimation of die attach cracking by using the fatigue crack propagation rate given by the Paris law.
The lifetime estimation methodologies based on the continuum damage mechanics and the cohesive zone model are able to deal with both crack initiation and crack propagation without assuming a crack with a crack tip having stress and strain singularity. Among the methodologies based on the continuum damage mechanics, Desai and Basaran et al. proposed the DSC approaches based on the unified constitutive modeling and the unified mechanics theory, and successfully applied these approaches to lifetime estimations of IC packaging. We are expecting that these approaches would be powerful methodologies for structural integrity evaluation of power modules, especially for lifetime estimation of die attach crack.
Commercial finite element codes such as marc, ansys, etc., do not have the standard functions for the continuum damage mechanics and the cohesive zone model. So great effort is required to add these functions to the commercial finite element codes or to develop new special-purpose finite element programs. The rainflow method [138–140] has been frequently utilized as a cycle count algorithm to predict the lifetime for complex load versus time histories. The lifetime estimation methodologies based on the continuum damage mechanics and the cohesive zone model are ill-matched to the rainflow method, because they are different from the Coffin–Manson type failure models. Thus, these methodologies need to calculate a damage variable D step by step. An enormous amount of computational time would be required for such a calculation.
New mechanical reliability problems of power modules are becoming apparent due to the change of manufacturing technology and operating environment for power modules. For example, heel cracking is becoming apparent as a failure mode of a bonding wire, when a molding material of power modules changes from silicone-gel to hard resin. Die attach failure in the thickness direction is becoming obvious, when a substrate with a low CTE is used to cope with high operating temperature. It becomes important to understand why such a failure mode is becoming apparent, and to develop a methodology for the lifetime estimation of such a failure mode.
The use of the failure models based on the temperature range $Δ T$ evaluated through the collector-emitter voltage V_CE in a power cycling test, which is described in IEC 60747 standard, is unsuitable for evaluating the lifetime of a local failure such as wire-liftoff. This is because the temperature obtained from V_CE is the average temperature of a Si/semiconductor die and not the local temperature at a failure location. The local temperature at the bonding part contributes to wire-liftoff. In summary, although the formulae for the wire-liftoff lifetimes obtained from power cycling tests based on IEC 60747 standard are useful for estimating the residual lifetimes of in-service power modules and for predicting the lifetimes of existing power modules, it is inappropriate to use them for estimating the wire-liftoff lifetimes in the design phase, because the temperatures determined based on IEC 60747 standard do not represent the local temperatures at the wire-bonding parts. The wire-liftoff lifetime data obtained from thermal cycling tests should be used for predicting the lifetimes in the design phase.

6 Proposals for Future Research

We should make a paradigm shift from estimating residual lifetimes/monitoring of in-service power modules using the purely empirical failure models such as $Δ T$ -based failure models to the prior assessment of lifetime in the design phase of power modules using the physically/mechanically based failure models such as $Δ ε_{in}$ -, $Δ W_{in}$ - and $Δ$ T*-based failure models or the DSC approaches based on the unified constitutive modeling and the unified mechanics theory. As described in Sec. 1, this is because enormous numbers of power modules will be utilized from now on in various fields, and, in such a circumstance, the malfunction and recall of power modules would be enormously unfavorable impact on our daily lives. So power device manufactures should provide reliable power modules to the market. From such a viewpoint, the present authors emphasize the need of the paradigm shift mentioned above, and propose the research policy for securing the structural integrity of power modules in the design phase, as follows:

Material testing is performed to obtain the material properties of a bond wire and a die attachment material. The following material characteristics are necessary:
- elastoplastic characteristics,
- creep characteristics,
- fatigue characteristics of the Coffin–Manson law type, and
- fatigue crack propagation characteristics of the Paris law type.
The physical quantity representing (c) and (d) should be a fracture mechanics parameter that can be uniquely determined. Among fracture mechanics parameters, a nonlinear fracture mechanics parameter should be selected, because ductile materials are utilized as a bonding wire/ribbon and die attach. We recommend the nonlinear fracture mechanics parameter T* range $Δ$ T*, because it has path-independence for various loading conditions including creep deformation as well as thermo-elastoplastic deformation.
We propose that mechanical loading tests are used to obtain the fatigue characteristics of the Coffin–Manson law type, and the fatigue crack propagation characteristics of the Paris law type as shown in Ref. [105] instead of time-consuming thermal cycle tests and power cycle tests. These test data are useful for the lifetime estimation of crack initiation and crack propagation in a die attach layer of a power module. The mechanical loading tests are also proposed for the estimation of wire-liftoff lifetime [141,142].
The lifetime estimation of power modules in design phase can be done by making full use of CAE tools such as finite element software together with the inelastic material properties of the materials consisting of a power module obtained from material testing and using the following procedure:

Step 1: Based on circuit design, the temperature distribution in a power module under the operating condition is obtained from electrothermal analysis.
Step 2: Based on the temperature distribution in a power module obtained in step 1 and the design dimensions of a power module, the physical quantities such as displacements, strains, and stresses are obtained from the thermo-elastoplastic creep analysis using the elastoplastic and creep characteristics of a bond wire and a die attachment material.
Step 3: The physical quantities relevant to fatigue characteristics and fatigue crack propagation characteristics such as $Δ ε_{in}$ ⁠, $Δ W_{in}$ ⁠, $Δ$ J, $Δ$ T*, etc., are obtained from those obtained in step 2.
Step 4: Based on the result of step 3, wire-liftoff lifetime and die attach cracking lifetime are, respectively, obtained from the fatigue characteristics of the Coffin–Manson law type and the fatigue crack propagation characteristics of the Paris law type.

The lifetime estimation in the design phase of a power module is performed according to the procedure mentioned above, and the lifetime obtained is checked whether or not it satisfies the design goal. If it does not satisfy the design goal, design change is necessary to secure the structural integrity of a power module.

Funding Data

JSPS Grant-in-Aid for Scientific Research (Grant No. JP18K03863; Funder ID: 10.13039/501100001691).

References

1.

Wild

,

R. N.

,

1972

, “

Fatigue Properties of Solder Joints

,”

Weld. J.

,

51

(

11

), pp.

521s

–

526s

.https://app.aws.org/wj/supplement/WJ_1972_11_s521.pdf

2.

Wu

,

W.

,

Held

,

M.

,

Jacob

,

P.

,

Scacco

,

P.

, and

Birolini

,

A.

,

1995

, “

Investigation on the Long Term Reliability of Power IGBT Modules

,”

Proceedings of 1995 International Symposium on Power Semiconductor Devices and ICs

, Yokohama, Japan, May 23–25, pp.

443

–

448

.10.1109/ISPSD.1995.515079

3.

Held

,

M.

,

Jacob

,

P.

,

Nicoletti

,

G.

,

Scacco

,

P.

, and

Poech

,

M.-H.

,

1999

, “

Fast Power Cycling Test for Insulated Gate Bipolar Transistor Modules in Traction Application

,”

Int. J. Electron.

,

86

(

10

), pp.

1193

–

1204

.10.1080/002072199132743

4.

Ciappa

,

M.

,

2002

, “

Selected Failure Mechanisms of Modern Power Modules

,”

Microelectron. Reliab.

,

42

(

4–5

), pp.

653

–

667

.10.1016/S0026-2714(02)00042-2

5.

Yang

,

S.

,

Xiang

,

D.

,

Bryant

,

A.

,

Mawby

,

P.

,

Ran

,

L.

, and

Tavner

,

P.

,

2010

, “

Condition Monitoring for Device Reliability in Power Electronic Converters: A Review

,”

IEEE Trans. Power Electron.

,

25

(

11

), pp.

2734

–

2752

.10.1109/TPEL.2010.2049377

6.

Yang

,

L.

,

Agyakwa

,

P. A.

, and

Johnson

,

C. M.

,

2013

, “

Physics-of-Failure Lifetime Prediction Models for Wire Bond Interconnects in Power Electronic Modules

,”

IEEE Trans. Device Mater. Reliab.

,

13

(

1

), pp.

9

–

17

.10.1109/TDMR.2012.2235836

7.

Moeini

,

R.

,

Tricoli

,

P.

,

Hemida

,

H.

, and

Baniotopoulos

,

C.

,

2018

, “

Increasing the Reliability of Wind Turbines Using Condition Monitoring of Semiconductor Devices: A Review

,”

IET Renewable. Power Gener.

,

12

(

2

), pp.

182

–

189

.10.1049/iet-rpg.2017.0477

8.

Hanif

,

A.

,

Yu

,

Y.

,

DeVoto

,

D.

, and

Khan

,

F.

,

2019

, “

A Comprehensive Review Toward the State-of-Art in Failure and Lifetime Predictions of Power Electronic Devices

,”

IEEE Trans. Power Electron.

,

34

(

5

), pp.

4729

–

4746

.10.1109/TPEL.2018.2860587

9.

Coffin

,

L. F.

,

1954

, “

A Study of the Effects of Cyclic Thermal Stresses in a Ductile Metal

,”

Trans. ASME

,

76

(

6

), pp.

931

–

950

.https://www.osti.gov/biblio/4363016-study-effects-cyclic-thermal-stresses-ductile-metal

10.

Manson

,

S. S.

,

1954

, “

Behavior of Materials Under Conditions of Thermal Stress

,” NACA, Washington DC, Report No.

NACA-TN-1170

.https://ntrs.nasa.gov/citations/19930092197

11.

Ramminger

,

S.

,

Seliger

,

N.

, and

Wachutka

,

G.

,

2000

, “

Reliability Model for Al Wire Bonds Subjected to Heel Crack Failures

,”

Microelectron. Reliab.

,

40

(

8–10

), pp.

1521

–

1525

.10.1016/S0026-2714(00)00139-6

12.

Merkle

,

L.

,

Kaden

,

T.

,

Sonner

,

M.

,

Gademann

,

A.

,

Turki

,

J.

,

Dresbach

,

C.

, and

Petzold

,

M.

,

2008

, “

Mechanical Fatigue Properties of Heavy Aluminum Wire Bonds for Power Applications

,”

Proceedings of Second Electronics Systemintegration Technology Conference

, Greenwich, UK, Sept. 1–4, pp.

1363

–

1367

.10.1109/ESTC.2008.4684554

13.

Merkle

,

L.

,

Sonner

,

M.

, and

Petzold

,

M.

,

2014

, “

Lifetime Prediction of Thick Aluminum Wire Bonds for Mechanical Cyclic Loads

,”

Microelectron. Reliab.

,

54

(

2

), pp.

417

–

424

.10.1016/j.microrel.2013.10.009

14.

Celnikier

,

Y.

,

Benabou

,

L.

,

Dupont

,

L.

, and

Coquery

,

G.

,

2011

, “

Investigation of the Heel Crack Mechanism in Al Connections for Power Electronics Modules

,”

Microelectron. Reliab.

,

51

(

5

), pp.

965

–

974

.10.1016/j.microrel.2011.01.001

15.

Nakajima

,

D.

,

Motoyama

,

K.

,

Masuhara

,

T.

,

Tamura

,

S.

, and

Fujimoto

,

K.

,

2019

, “

The Failure Mode of Wire Bond in Transfer Molded Power Modules

,”

Proceedings of 25th Symposium on Microjoining and Assembly Technology in Elecctronics (Mate 2019)

, Yokohama, Japan, Jan. 29–30, pp.

123

–

126

.

16.

Choi

,

U. –M.

,

Blaabjerg

,

F.

, and

Jørgensen

,

S.

,

2017

, “

Study on Effect of Junction Temperature Swing Duration on Lifetime of Transfer Molded Power IGBT Modules

,”

IEEE Trans. Power Electron.

,

32

(

8

), pp.

6434

–

6443

.10.1109/TPEL.2016.2618917

17.

Zeng

,

G.

,

Borucki

,

L.

,

Wenzel

,

O.

,

Schilling

,

O.

, and

Lutz

,

J.

,

2018

, “

First Results of Development of a Lifetime Model for Transfer Molded Discrete Power Devices

,” International Exhibition and Conference for Power Electronics, Intelligent Motion, Renewable Energy and Energy Management (

PCIM Europe 2018

), Nuremberg, Germany, June 5–7, pp.

706

–

713

.https://ieeexplore.ieee.org/document/8402917

18.

Kobayashi

,

H.

, and

Terada

,

H.

, “

Mid-Air Explosion of Comet I over the Mediterranean Sea, Failure Knowledge Data Base/100 Selected Cases

,” Association for Study of Failure, Tokyo, Japan.http://www.shippai.org/fkd/en/cfen/CB1071012.html

19.

Ikeda

,

T.

,

Miyazaki

,

N.

,

Kudo

,

K.

,

Arita

,

K.

, and

Yakiyama

,

H.

,

1999

, “

Failure Estimation of Semiconductor Chip During Wire Bonding Process

,”

ASME J. Electron. Packag.

,

121

(

2

), pp.

85

–

91

.10.1115/1.2792672

20.

Hizukuri

,

M.

, and

Asano

,

T.

,

2000

, “

Measurement of Dynamic Strain During Ultrasonic Au Bump Formation on Si Chip

,”

Jpn. J. Appl. Phys.

,

39

(

Pt. 1 4B

), pp.

2478

–

2482

.10.1143/JJAP.39.2478

21.

Choi

,

U.-M.

,

Ma

,

K.

, and

Blaabjerg

,

F.

,

2018

, “

Validation of Lifetime Prediction of IGBT modules Based on Linear Damage Accumulation by Means of Superimposed Power Cycling Tests

,”

IEEE Trans. Ind. Electron.

,

65

(

4

), pp.

3520

–

3529

.10.1109/TIE.2017.2752142

22.

Zhang

,

Y.

,

Wang

,

H.

,

Wang

,

Z.

,

Yang

,

Y.

, and

Blaabjerg

,

F.

,

2017

, “

Impact of Lifetime Model Selections on the Reliability Prediction of IGBT Modules in Modular Multilevel Converters

,”

Proceedings of 2017 IEEE Energy Conversion Congress and Exposition

, Cincinnati, OH, Oct. 1–5, pp.

4202

–

4207

.10.1109/ECCE.2017.8096728

23.

Yu

,

Q.

,

Shiratori

,

M.

,

Kaneko

,

S.

,

Ishihara

,

T.

, and

Wang

,

S.

,

1998

, “

Analytical and Experimental Hybrid Study on Thermal Fatigue Strength of Electronic Solder Joints (2nd Report, Evaluation by Isothermal Mechanical Fatigue Test)

,”

Trans. JSME (Ser. A)

,

64

(

619

), pp.

558

–

563

.10.1299/kikaia.64.558

24.

Manson

,

S. S.

,

1966

,

Thermal Stress and Low Cycle Fatigue

,

McGraw-Hill

,

New York

.

25.

Wang

,

H.

,

Ma

,

K.

, and

Blaabjerg

,

F.

,

2012

, “

Design for Reliability of Power Electronic Systems

,” Proceedings of 38th Annual Conference of the IEEE Industrial Electronics Society (

IECIN 2012

), Montreal, QC, Canada, Oct. 25–28, pp.

33

–

44

.10.1109/IECON.2012.6388833

26.

Bayerer

,

R.

,

Herrmann

,

T.

,

Licht

,

T.

,

Lutz

,

J.

, and

Feller

,

M.

,

2008

, “

Model for Power Cycling Lifetime of IGBT Modules—Various Factors Influencing Lifetime

,”

Proceedings of Fifth International Conference on Integrated Power Systems

, Nuremberg, Germany, Mar. 11–13.https://ieeexplore.ieee.org/document/5755669

27.

Norris

,

K. C.

, and., and

Landzberg

,

A. H.

,

1969

, “

Reliability of Controlled Collapse Interconnections

,”

IBM J. Res. Dev.

,

13

(

3

), pp.

266

–

271

.10.1147/rd.133.0266

28.

Kovacevic-Badstuebner

,

I. F.

,

Kolar

,

J. W.

, and

Shilling

,

U.

,

2015

, “

Modelling for the Lifetime Prediction of Power Semiconductor Modules

,”

Reliability of Power Electronic Converter System

,

IET

,

London

, pp.

103

–

140

.10.1049/PBPO080E_ch5

29.

Yang

,

X.

,

Lin

,

Z.

,

Ding

,

J.

, and

Long

,

Z.

,

2019

, “

Lifetime Prediction of IGBT Modules in Suspension Choppers of Medium/Low-Speed Maglev Train Using an Energy-Based Approach

,”

IEEE Trans. Power Electron.

,

34

(

1

), pp.

738

–

747

.10.1109/TPEL.2018.2812732

30.

Sasaki

,

K.

, and

Iwasa

,

N.

,

2008

, “

Thermal and Structural Simulation Techniques for Estimating Fatigue Life of an IGBT Model

,”

Proceedings of 20th International Symposium on Power Semiconductor Devices and IC's

, Orlando, FL, May 18–22, pp.

181

–

184

.10.1109/ISPSD.2008.4538928

31.

Matsunaga

,

T.

, and

Uegai

,

Y.

,

2006

, “

Thermal Fatigue Life Evaluation of Aluminum Wire Bonds

,”

Proceedings of First Electronics System Integration Technology Conference

, Dresden, Germany, Sept. 5–7, pp.

726

–

731

.10.1109/ESTC.2006.280092

32.

Yamada

,

Y.

,

Takaku

,

Y.

,

Yagi

,

Y.

,

Nakagawa

,

I.

,

Atsumi

,

T.

,

Shirai

,

M.

,

Ohnuma

,

I.

, and

Ishida

,

K.

,

2007

, “

Reliability of Wire-Bonding and Solder Joint for High Temperature Operation of Power Semiconductor Device

,”

Microelectron. Reliab.

,

47

(

12

), pp.

2147

–

2151

.10.1016/j.microrel.2007.07.102

33.

Agyakwa

,

P. A.

,

Corfield

,

M. R.

,

Yang

,

L.

,

Li

,

J. F.

,

Marques

,

V. M. F.

, and

Johnson

,

C. M.

,

2011

, “

Microstructural Evolution of Ultrasonically Bonded High Purity Al Wire During Extended Range Thermal Cycling

,”

Microelectron. Reliab.

,

51

(

2

), pp.

406

–

415

.10.1016/j.microrel.2010.08.018

34.

Yang

,

L.

,

Agyakwa

,

P. A.

, and

Johnson

,

C. M.

,

2011

, “

A Time-Domain Physics-of-Failure Model for the Lifetime Prediction of Wire Bond Connection

,”

Microelectron. Reliab.

,

51

(

9–11

), pp.

1882

–

1886

.10.1016/j.microrel.2011.07.052

35.

Morrow

,

J.

,

1965

, “

Cyclic Plastic Strain Energy and Fatigue of Metals

,”

Internal Friction, Damping, and Cyclic Plasticity

, B. Kazan, ed., ASTM International, West Conshohocken, PA, pp.

45

–

85s

.10.1520/STP43764S

36.

Kanda

,

Y.

,

Kariya

,

Y.

, and

Oto

,

Y.

,

2012

, “

Influence of Cyclic Strain-Hardening Exponent on Fatigue Ductility Exponent for a Sn–Ag–Cu Micro-Solder Joint

,”

J. Electron. Mater.

,

41

(

3

), pp.

580

–

587

.10.1007/s11664-011-1830-7

37.

Clech

,

J.-P.

,

1997

, “

Solder Reliability Solutions: A PC-Based Design-for-Reliability Tool

,”

Solder. Surf. Mount Technol.

,

9

(

2

), pp.

45

–

54

.10.1108/09540919710800638

38.

Ciappa

,

M.

,

Carbognani

,

F.

,

Cova

,

P.

, and

Fichtner

,

W.

,

2002

, “

A Novel Thermomechanics-Based Lifetime Prediction Model for Cycle Fatigue Failure Mechanisms in Power Semiconductors

,”

Microelectron. Reliab.

,

42

(

9–11

), pp.

1653

–

1658

.10.1016/S0026-2714(02)00206-8

39.

Riedel

,

G. J.

,

Schmidt

,

R.

,

Liu

,

C.

,

Beyer

,

H.

, and

Alaperä

,

I.

,

2012

, “

Reliability of Large Area Solder Joints Within IGBT Modules: Numerical Modeling and Experimental Results

,”

Proceedings of Seventh International Conference on Integrated Power Electronics Systems

, Nuremberg, Germany, Mar. 6–8, Paper No. 06.4.https://ieeexplore.ieee.org/document/6170637

40.

Shishido

,

N.

,

Hayama

,

Y.

,

Morooka

,

W.

,

Hagihara

,

S.

, and

Miyazaki

,

N.

,

2019

, “

Application of Nonlinear Fracture Mechanics Parameter to Predicting Wire-Liftoff Lifetime of Power Module at Elevated Temperature

,”

IEEE J. Emerg. Sel. Top. Power Electron.

,

7

(

3

), pp.

1604

–

1614

.10.1109/JESTPE.2019.2914244

41.

Rice

,

J. R.

,

1968

, “

A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks

,”

ASME J. Appl. Mech.

,

35

(

2

), pp.

379

–

386

.10.1115/1.3601206

42.

Lamba

,

H. S.

,

1975

, “

The J-Integral Applied to Cyclic Loading

,”

Eng. Fract. Mech.

,

7

(

4

), pp.

693

–

703

.10.1016/0013-7944(75)90025-9

43.

Dowling

,

N. E.

, and

Begley

,

J. A.

,

1976

, “

Fatigue Crack Growth During Gross Plasticity and the J-Integral

,”

Mechanics of Crack Growth

,

J.

Rice

and

P.

Paris

, ed.,

ASTM International

,

West Conshohocken, PA

, pp.

82

–

103

.10.1520/STP33940S

44.

Wüthrich

,

C.

,

1982

, “

The Extension of the J-Integral Concept to Fatigue Cracks

,”

Int. J. Fract.

,

20

(

2

), pp.

R35

–

R37

.10.1007/BF01141264

45.

Tanaka

,

K.

,

1983

, “

The Cyclic J-Integral as a Criterion for Fatigue Crack Growth

,”

Int. J. Fract.

,

22

(

2

), pp.

91

–

104

.10.1007/BF00942715

46.

Kubo

,

S.

,

Yafuso

,

T.

,

Nohara

,

M.

,

Ishimaru

,

T.

, and

Ohji

,

K.

,

1989

, “

Investigation on Path-Integral Expression of the J-Integral Range Using Numerical Simulations of Fatigue Crack Growth

,”

JSME Int. J. (Ser. I)

,

32

(

2

), pp.

237

–

244

.10.1299/jsmea1988.32.2_237

47.

Atluri

,

S. N.

,

1982

, “

Path-Independent Integrals in Finite Elasticity and Inelasticity, With Body Forces, Inertia, and Arbitrary Crack-Face Conditions

,”

Engng. Fract. Mech.

,

16

(

3

), pp.

341

–

364

.10.1016/0013-7944(82)90113-8

48.

Brust

,

F. W.

,

McGowan

,

J. J.

, and

Atluri

,

S. N.

,

1986

, “

A Combined Numerical/Experimental Study of Ductile Crack Growth After a Large Unloading, Using T*, J, and CTOA Criteria

,”

Eng. Fract. Mech.

,

23

(

3

), pp.

537

–

550

.10.1016/0013-7944(86)90161-X

49.

Smelser

,

R. E.

, and

Gurtin

,

M. E.

,

1977

, “

On the J-Integral for Bi-Material Bodies

,”

Int. J. Fract.

,

13

(

3

), pp.

382

–

384

.10.1007/BF00040155

50.

Banks-Sills

,

L.

, and

Volpert

,

Y.

,

1991

, “

Application of the Cyclic J-Integral to Fatigue Crack Propagation of Al 2024-T351

,”

Engng. Fract. Mech.

,

40

(

2

), pp.

355

–

370

.10.1016/0013-7944(91)90270-B

51.

Gasiak

,

G.

, and

Rozumek

,

D.

,

2004

, “

ΔJ-Integral Range Estimation for Fatigue Crack Growth Rate Description

,”

Int. J. Fatig

,

26

(

2

), pp.

135

–

140

. pp.10.1016/S0142-1123(03)00111-7

52.

Shahani

,

A. R.

,

Moayeri Kashani

,

H.

,

Rastegar

,

M.

, and

Botshekanan Dehkordi

,

M.

,

2009

, “

A Unified Model for the Fatigue Crack Growth Rate in Variable Stress Ratio

,”

Fatig. Fract. Eng. Mater. Struct.

,

32

(

2

), pp.

105

–

118

.10.1111/j.1460-2695.2008.01315.x

53.

Takeda

,

T.

,

Shindo

,

Y.

, and

Narita

,

F.

,

2011

, “

Vacuum Crack Growth Behavior of Austenitic Stainless Steel Under Fatigue Loading

,”

Streng. Mater.

,

43

(

5

), pp.

532

–

536

.10.1007/s11223-011-9324-7

54.

Božić

,

Ž.

,

Mlikota

,

M.

, and., and

Schmauder

,

S.

,

2011

, “

Application of the

Δ K

⁠,

Δ J

and

Δ

CTOD Parameters in Fatigue Crack Growth Modelling

,”

Tech. Gazette

,

18

(

3

), pp.

459

–

466

.https://core.ac.uk/download/pdf/14440251.pdf

55.

Ktari

,

A.

,

Baccar

,

M.

,

Shah

,

M.

,

Haddar

,

N.

,

Ayedi

,

H. F.

, and

Rezai-Aria

,

F.

,

2014

, “

A Crack Propagation Criterion Based on

Δ

CTOD Measured With 2D-Digital Image Correlation Technique

,”

Fatig. Fract. Eng. Mater. Struct.

,

37

(

6

), pp.

682

–

694

.10.1111/ffe.12153

56.

Li

,

L.

,

Yang

,

Y. H.

,

Xu

,

Z.

,

Chen

,

G.

, and

Chen

,

X.

,

2014

, “

Fatigue Crack Growth Law of API X80 Pipeline Steel Under Various Stress Ratios Based on J-Integral

,”

Fatig. Fract. Eng. Mater. Struct.

,

37

(

10

), pp.

1124

–

1135

.10.1111/ffe.12193

57.

Furuhashi

,

I.

, and

Wakai

,

T.

,

1995

, “

Revisions of Fracture Mechanics Parameters Analysis Code CANIS-J(2D)

,” JAEA, Tokai-mura, Japan, Report No. PNC-TN9410

95

–

080

.

58.

Ridout

,

S.

, and

Bailey

,

C.

,

2007

, “

Review of Methods to Predict Solder Joint Reliability Under Thermo-Mechanical Cycling

,”

Fatig. Fract. Eng. Mater. Struct.

,

30

(

5

), pp.

400

–

412

.10.1111/j.1460-2695.2006.01065.x

59.

Wong

,

E. H.

,

van Driel

,

W. D.

,

Dasgupta

,

A.

, and

Pecht

,

M.

,

2016

, “

Creep Fatigue Models of Solder Joints: A Critical Review

,”

Microelectron. Reliab.

,

59

(

1

), pp.

1

–

12

.10.1016/j.microrel.2016.01.013

60.

Yao

,

Y.

,

Long

,

X.

, and

Keer

,

L. M.

,

2017

, “

A Review of Recent Research on the Mechanical Behavior of Lead-Free Solders

,”

ASME Appl. Mech. Rev.

,

69

(

4

), Paper No. 040802.10.1115/1.4037462

61.

Chen

,

G.

,

Zhao

,

X.

, and

Wu

,

H.

,

2017

, “

A Critical Review of Constitutive Models for Solders in Electronic Packaging

,”

Adv. Mech. Eng.

,

9

(

8

), pp.

1

–

21

.10.1177/1687814017714976

62.

Huang

,

X.

,

Wu

,

W.-F.

, and

Chou

,

P.-L.

,

2012

, “

Fatigue Life and Reliability Prediction of Electronic Packages Under Thermal Cycling Conditions Through FEM Analysis and Acceleration Models

,”

Proceedings of 14th International Conference on Electronic Materials and Packaging

, Macao, China, Dec.

13

–

16

.10.1109/EMAP.2012.6507914

63.

Salmela

,

O.

,

2007

, “

Acceleration Factors for Lead-Free Solder Materials

,”

IEEE Trans. Comp. Packag. Technol.

,

30

(

4

), pp.

700

–

707

.10.1109/TCAPT.2007.900076

64.

Dauksher

,

W.

,

2008

, “

A Second-Level SAC Solder-Joint Fatigue-Life Prediction Methodology

,”

IEEE Trans. Device Mater. Reliab.

,

8

(

1

), pp.

168

–

173

.10.1109/TDMR.2007.912253

65.

Salmela

,

O.

,

Putaala

,

J.

,

Nousiainen

,

O.

,

Uusimäki

,

A.

,

Särkkä

,

J.

, and

Tammenmaa

,

M.

,

2016

, “

Multipurpose Lead-Free Reliability Prediction Model

,”

Proceedings of 2016 Pan Pacific Microelectronics Symposium

, Big Island, HI, Jan.

25

–

28

.10.1109/PanPacific.2016.7428389

66.

Darveaux

,

R.

, and

Banerji

,

K.

,

1991

, “

Fatigue Analysis of Flip Chip Assemblies Using Thermal Stress Simulations and a Coffin-Manson Relation

,”

Proceedings of 41st Electronic Components and Technology Conference

, Atlanta, GA, May 11–16, pp.

797

–

805

.10.1109/ECTC.1991.163971

67.

Pang

,

J. H. L.

,

Tan

,

T.-I.

, and

Sitaraman

,

S. K.

,

1998

, “

Thermo-Mechanical Analysis of Solder Joints Fatigue and Creep in a Flip Chip on Board Package Subjected to Temperature Cycling Loading

,”

Proceedings of 48th Electronic Components and Technology Conference

, Seattle, WA, May 25–28, pp.

878

–

883

.10.1109/ECTC.1998.678811

68.

Schubert

,

A.

,

Dudek

,

R.

,

Auerswald

,

E.

,

Gollhardt

,

A.

,

Michel

,

B.

, and

Reichl

,

H.

,

2003

, “

Fatigue Life Models for SnAgCu and SnPb Solder Joints Evaluated by Experiments and Simulation

,”

Proceedings of 53rd Electronic Components and Technology Conference

, New Orleans, LA, May 27–30, pp.

603

–

610

.10.1109/ECTC.2003.1216343

69.

Kim

,

I. H.

,

Park

,

T. S.

,

Yang

,

S. Y.

, and

Lee

,

S. B.

,

2005

, “

A Comparative Study of the Fatigue Behavior of SnAgCu and SnPb Solder Joints

,”

Key Engng. Mater.

,

297–300

, pp.

831

–

836

.10.4028/www.scientific.net/KEM.297-300.831

70.

Karppinen

,

J. S.

,

Li

,

J.

,

Mattila

,

T. T.

, and

Paulasto-Kröckel

,

M.

,

2010

, “

Thermomechanical Reliability Characterization of a Handheld Product in Accelerated Tests and Use Environment

,”

Microelectron. Reliab.

,

50

(

12

), pp.

1994

–

2000

.10.1016/j.microrel.2010.07.011

71.

Hansen

,

P.

, and

McCluskey

,

P.

,

2007

, “

Failure Models in Power Device Interconnects

,”

Proceedings of 2007 European Conference on Power Electronics and Applications

, Aaborg, Denmark, Sept.

2

–

5

.10.1109/EPE.2007.4417406

72.

O'Keefe

,

M.

, and

Vlahinos

,

A.

,

2009

, “

Impacts of Cooling Technology on Solder Fatigue for Power Modules in Electric Traction Drive Vehicles

,”

Proceedings of 2009 IEEE Vehicle Power and Propulsion Conference

, Dearborn, MI, Sept. 7–10, pp.

1182

–

1188

. 10.1109/VPPC.2009.5289717

73.

Barbagallo

,

C.

,

Malgioglio

,

G. L.

,

Petrone

,

G.

, and

Cammarata

,

G.

,

2017

, “

Thermal Fatigue Life Evaluation of SnAgCu Solder Joints in a Multi-Chip Power Module

,”

IOP Conf. Ser.: J. Phys.: Conf. Ser.

,

841

, Paper No. 012014.10.1088/1742-6596/841/1/012014

74.

Cavallaro

,

D.

,

Greco

,

R.

, and

Bazzano

,

G.

,

2018

, “

Effect of Solder Material Thickness on Power MOSFT Reliability by Electro-Thermo-Mechanical Simulations

,”

Microelectron. Reliab.

,

88–90

, pp.

1168

–

1171

.10.1016/j.microrel.2018.07.082

75.

Darveaux

,

R.

,

1997

, “

Solder Joint Fatigue Life Model

,”

Proceedings of 1997 TMS Annual Meeting

, Orlando, FL, Feb. 10–13, pp.

213

–

218

.

76.

Amagai

,

M.

, and

Nakao

,

M.

,

1998

, “

Ball Grid Array (BGA) Packages With the Copper Core Solder Balls

,”

Proceedings of 48th Electronic Components and Technology Conference

, Seattle, WA, May 25–28, pp.

692

–

701

.10.1109/ECTC.1998.678773

77.

Amagai

,

M.

,

1999

, “

Characterization of Chip Scale Packaging Materials

,”

Microelectron. Reliab.

,

39

(

9

), pp.

1365

–

1377

.10.1016/S0026-2714(99)00059-1

78.

Zahn

,

B. A.

,

2003

, “

Solder Joint Fatigue Life Model Methodology for 63Sn37Pb and 95.5Sn4Ag0.5Cu Materials

,”

Proceedings of 53th Electronic Components and Technology Conference

, New Orleans, LA, May 27–28, pp.

83

–

94

.10.1109/ECTC.2003.1216261

79.

Thébaud

,

J.-M.

,

Woirgard

,

E.

,

Zardini

,

C.

,

Azzopardi

,

S.

,

Briat

,

O.

, and

Vinassa

,

J.-M.

,

2003

, “

Strategy for Designing Accelerated Aging Tests to Evaluate IGBT Power Modules Lifetime in Real Operation Mode

,”

IEEE Trans. Comp. Packag. Technol.

,

26

(

2

), pp.

429

–

438

.10.1109/TCAPT.2003.815112

80.

Hu

,

B.

,

Gonzalez

,

J. O.

,

Ran

,

L.

,

Ren

,

H.

,

Zeng

,

Z.

,

Lai

,

W.

,

Gao

,

B.

,

Alatise

,

O.

,

Lu

,

H.

,

Bailey

,

C.

, and

Mawby

,

P.

,

2017

, “

Failure and Reliability Analysis of a SiC Power Module Based on Stress Comparison to a Si Device

,”

IEEE Trans. Device Mater. Reliab.

,

17

(

4

), pp.

727

–

737

.10.1109/TDMR.2017.2766692

81.

Nakajima

,

Y.

,

Ono

,

K.

, and

Kariya

,

Y.

,

2019

, “

Evaluation of Fatigue Crack Propagation of Sn-5.0Sb/Cu Joint Using Inelastic Strain Energy Density

,”

Mater. Trans.

,

60

(

6

), pp.

876

–

881

.10.2320/matertrans.MH201813

82.

Hutchinson

,

J. W.

,

1968

, “

Singular Behaviour at the End of a Tensile Crack in a Hardening Material

,”

J. Mech. Phys. Solids

,

16

(

1

), pp.

13

–

31

.10.1016/0022-5096(68)90014-8

83.

Rice

,

J. R.

, and

Rosengren

,

G. F.

,

1968

, “

Plane Strain Deformation Near a Crack Tip in a Power Law Hardening Material

,”

J. Mech. Phys. Solids

,

16

(

1

), pp.

1

–

12

.10.1016/0022-5096(68)90013-6

84.

Paris

,

P.

, and

Erdogan

,

F.

,

1963

, “

A Critical Analysis of Crack Propagation Laws

,”

ASME J. Basic Eng.

,

85

(

4

), pp.

528

–

533

.10.1115/1.3656900

85.

Déplanque

,

S.

,

Nüchter

,

W.

,

Wunderle

,

B.

,

Schacht

,

R.

, and

Michel

,

B.

,

2006

, “

Lifetime Prediction of SnPb and SnAgCu Solder Joints of Chips on Copper Substrate Based on Crack Propagation FE-Analysis

,” Proceedings of Seventh International Conference on Thermal, Mechanical and Multiphysics Simulation and Experiments in Micro-Electronics and Micro-Systems (

EuroSimuE 2006

), Como, Italy, Apr. 24–26.10.1109/ESIME.2006.1643976

86.

Lu

,

H.

,

Tilford

,

T.

,

Bailey

,

C.

, and

Newcombe

,

D. R.

,

2007

, “

Lifetime Prediction for Power Electronics Module Substrate Mount-Down Solder Interconnect

,”

Proceedings of 2007 International Symposium on High Density Packaging and Microsystem Integration

, Shanghai, China, June

26

–

28

.10.1109/HDP.2007.4283558

87.

Yin

,

C. Y.

,

Lu

,

H.

,

Musallam

,

M.

,

Bailey

,

C.

, and

Johnson

,

C. M.

,

2010

, “

In-Service Reliability Assessment of Solder Interconnect in Power Electronics Module

,”

Proceedings of 2010 Prognostics and System Health Management Conference

, Macao, China, Jan. 12–14, Paper No. MU3095.10.1109/PHM.2010.5413346

88.

Kostandyan

,

E. E.

, and

Sørensen

,

J. D.

,

2011

, “

Reliability Assessment of Solder Joints in Power Electronic Modules by Crack Damage Model for Wind Turbine Applications

,”

Energies

,

4

(

12

), pp.

2236

–

2248

.10.3390/en4122236

89.

Mukai

,

M.

,

Hirohata

,

K.

,

Takahashi

,

H.

,

Kawakami

,

T.

, and

Takahashi

,

K.

,

2005

, “

Damage Path Simulation of Solder Joints in QFP

,”

ASME

Paper No. IPACK2005-73297.10.1115/IPACK2005-73297

90.

Tanie

,

H.

,

Terasaki

,

T.

, and

Naka

,

Y.

,

2005

, “

A New Method for Evaluating Fatigue Life of Micro-Solder Joints in Semiconductor Structures

,”

ASME

Paper No. IPACK2005-73331.10.1115/IPACK2005-73331

91.

Anzawa

,

T.

,

Yu

,

Q.

,

Shibutani

,

T.

, and

Shiratori

,

M.

,

2007

, “

Reliability Evaluation for Power Electronics Devices Using Electrical Thermal and Mechanical Analysis

,”

Proceedings of Nineth Electronics Packaging Technology Conference

, Singapore, Dec. 10–12, pp.

94

–

99

.10.1109/EPTC.2007.4469833

92.

Shinohara

,

K.

, and

Yu

,

Q.

,

2010

, “

Reliability Evaluation of Power Semiconductor Devices Using Coupled Analysis Simulation

,”

Proceedings of 12th IEEE Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems

, Las Vegas, NV, June

2

–

5

.10.1109/ITHERM.2010.5501267

93.

Darveaux

,

R.

,

2000

, “

Effect of Simulation Methodology on Solder Joint Crack Growth Correlation

,”

Proceedings of 50th Electronic Components and Technology Conference

, Las Vegas, NV, May 21–24, pp.

1048

–

1058

.10.1109/ECTC.2000.853299

94.

Darveaux

,

R.

,

2002

, “

Effect of Simulation Methodology on Solder Joint Crack Growth Correlation and Fatigue Life Prediction

,”

ASME J. Electron. Packag.

,

124

(

3

), pp.

147

–

154

.10.1115/1.1413764

95.

Zhang

,

L.

,

Sitaraman

,

R.

,

Patwardhan

,

V.

,

Nguyen

,

L.

, and

Kelkar

,

N.

,

2003

, “

Solder Joint Reliability Model With Modified Darveaux's Equations for the Micro SMD Wafer Level-Chip Scale Package Family

,”

Proceedings of 53rd Electronic Components and Technology Conference

, New Orleans, LA, May 27–30, pp.

572

–

577

.10.1109/ECTC.2003.1216338

96.

Hossain

,

M. M.

,

Jagarkal

,

S. G.

,

Agonafer

,

D.

,

Lulu

,

M.

, and

Reh

,

S.

,

2007

, “

Design Optimization and Reliability of PWR Level Electronic Package

,”

ASME J. Electron. Packag.

,

129

(

1

), pp.

9

–

18

.10.1115/1.2429704

97.

Bai

,

J. G.

,

Calata

,

J. N.

, and

Lu

,

G.-Q.

,

2006

, “

Discussion on the Reliability Issues of Solder-Bump and Direct-Solder Bonded Power Device Packages Having Double-Sided Cooling Capability

,”

ASME J. Electron. Packag.

,

128

(

3

), pp.

208

–

214

.10.1115/1.2229218

98.

Xie

,

X.

,

Bi

,

X.

, and

Li

,

G.

,

2009

, “

Thermal-Mechanical Fatigue Reliability of PbSnAg Solder Layer of Die Attachment for Power Electronic Devices

,”

Proceedings of Tenth International Conference on Electronic Packaging Technology and High Density Packaging

, Beijing, China, Aug. 10–13, pp.

1181

–

1185

.10.1109/ICEPT.2009.5270615

99.

Chen

,

H.-C.

,

Guo

,

S.-W.

, and

Liao

,

H.-K.

,

2018

, “

The Solder Life Prediction Model of Power Module Under Thermal Cycling Test (TCT)

,”

Proceedings of 13th International Microsystems, Packaging, Assembly and Circuit Technology Conference

, Taipei, Taiwan, Oct. 24–26, pp.

104

–

107

.10.1109/IMPACT.2018.8625778

100.

Kimura

,

R.

,

Kariya

,

Y.

,

Mizumura

,

N.

, and

Sasaki

,

K.

,

2018

, “

Effect of Sintering Temperature on Fatigue Crack Propagation Rate of Sintered Ag Nanoparticles

,”

Mater. Trans.

,

59

(

4

), pp.

612

–

619

.10.2320/matertrans.M2017392

101.

Sato

,

T.

,

Kariya

,

Y.

,

Takahashi

,

H.

,

Nakamura

,

T.

, and

Aiko

,

Y.

,

2019

, “

Evaluation of Fatigue Crack Propagation Behavior of Pressurized Sintered Ag Nanoparticles and Its Application to Thermal Fatigue Life Prediction of Sintered Joint

,”

Mater. Trans.

,

60

(

6

), pp.

850

–

857

.10.2320/matertrans.MH201802

102.

Sundararajan

,

R.

,

McCluskey

,

P.

, and

Azarm

,

S.

,

1998

, “

Semi Analytic Model for Thermal Fatigue Failure of Die Attach in Power Electronic Building Blocks

,”

Proceedings Fourth International High Temperature Electronic Conference

, Albuquerque, NM, June 14–18, pp.

94

–

102

.10.1109/HITEC.1998.676768

103.

Pao

,

Y.-H.

,

1992

, “

A Fracture Mechanics Approach to Thermal Fatigue Prediction of Solder Joint

,”

IEEE Trans. Comp., Hybrids Manuf. Technol.

,

15

(

4

), pp.

559

–

570

.10.1109/33.159886

104.

Riedel

,

H.

,

1981

, “

Creep Deformation at Crack Tips in Elastic-Viscoplastic Solids

,”

J. Mech. Phys. Solids

,

29

(

1

), pp.

35

–

49

.10.1016/0022-5096(81)90014-4

105.

Shioda

,

R.

,

Kariya

,

Y.

,

Mizumura

,

N.

, and

Sasaki

,

K.

,

2017

, “

Low-Cycle Fatigue Life and Fatigue Propagation of Sintered Ag Nanoparticles

,”

ASME J. Electron. Mater.

,

46

(

2

), pp.

1155

–

1162

.10.1007/s11664-016-5068-2

106.

Asada

,

Y.

,

Shimakawa

,

T.

,

Kitagawa

,

M.

,

Kodaira

,

T.

,

Wada

,

Y.

, and

Asayama

,

T.

,

1992

, “

Analytical Evaluation Method of J-Integral in Creep-Fatigue Fracture for Type 304 Stainless Steel

,”

Nucl. Engng. Des.

,

133

(

3

), pp.

361

–

367

.10.1016/0029-5493(92)90162-O

107.

Murakami

,

S.

,

2012

,

Continuum Damage Mechanics: Approach to the Analysis of Damage and Fracture

,

Springer-Verlag

, Dordrecht/Heidelberg/London/New York.

108.

Xiao

,

H.

,

Li

,

X.

,

Liu

,

N.

, and

Yan

,

Y.

,

2011

, “

A Damage Model for SnAgCu Solder Under Thermal Cycling

,”

Proceedings of 12th International Conference on Electronic Packaging Technology and High Density Packaging

, Shanghai, China, Aug. 8–11, pp.

772

–

776

.10.1109/ICEPT.2011.6066946

109.

Yao

,

Y.

,

He

,

X.

,

Keer

,

L. M.

, and

Fine

,

M. E.

,

2015

, “

A Continuum Damage Mechanics-Based Unified Creep and Plasticity Model for Solder Materials

,”

Acta Mater.

,

83

, pp.

160

–

168

.10.1016/j.actamat.2014.09.051

110.

Rajaguru

,

P.

,

Lu

,

H.

, and

Bailey

,

C.

,

2015

, “

A Time Dependent Damage Indicator Model for Sn3.5Ag Solder Layer in Power Electronic Module

,”

Microelectron. Reliab

,.

55

(

11

), pp.

2371

–

2381

.10.1016/j.microrel.2015.07.047

111.

Lai

,

W.

,

Chen

,

M.

,

Ran

,

L.

,

Xu

,

S.

,

Pan

,

L.-M.

,

Alatise

,

O.

, and

Mawby

,

P.

,

2018

, “

Study on Lifetime Prediction Considering Fatigue Accumulative Effect for Die-Attach Solder Layer in an IGBT Module

,”

IEEJ Trans. Electr. Electron. Eng.

,

13

(

4

), pp.

613

–

621

.10.1002/tee.22607

112.

Darveaux

,

R.

,

Enayet

,

S.

,

Reichman

,

C.

,

Berry

,

C. J.

, and

Zafar

,

N.

,

2011

, “

Crack Initiation and Growth in WLCSP Solder Joints

,”

Proceedings of the 61st Electronics Components Technology Conference

, Lake Buena Vista, FL, May 31–June 3, pp.

940

–

953

.10.1109/ECTC.2011.5898624

113.

Desai

,

C. S.

, and

Toth

,

J.

,

1996

, “

Distributed State Constitutive Modelling Based on Stress-Strain and Nondestructive Behavior

,”

Int. J. Solids Struct.

,

33

(

11

), pp.

1619

–

1650

.10.1016/0020-7683(95)00115-8

114.

Desai

,

C. S.

,

Basaran

,

C.

, and

Zhang

,

W.

,

1997

, “

Numerical Algorithms and Mesh Dependence in the Distributed State Concept

,”

Int. J. Num. Meth. Eng.

,

40

(

16

), pp.

3059

–

3083

.10.1002/(SICI)1097-0207(19970830)40:16<3059::AID-NME182>3.0.CO;2-S

115.

Desai

,

C. S.

,

Basaran

,

C.

,

Dishongh

,

T.

, and

Prince

,

J. L.

,

1998

, “

Thermomechanical Analysis in Electronic Packaging With Unified Constitutive Model for Materials and Joints

,”

IEEE Trans. Comp. Packag. Technol. Part B

,

21

(

1

), pp.

87

–

97

.10.1109/96.659511

116.

Basaran

,

C.

,

Desai

,

C. S.

, and

Kundu

,

T.

,

1998

, “

Thermomechanical Finite Element Analysis of Problems in Electronic Packaging Using the Distributed State Concept—Part 1: Theory and Formulation

,”

ASME J. Electron. Packag.

,

120

(

1

), pp.

41

–

47

.10.1115/1.2792284

117.

Desai

,

C. S.

,

2016

, “

Distributed State Concept as Unified Constitutive Modeling Approach

,”

J. Rock Mech. Geotech. Eng.

,

8

(

3

), pp.

277

–

293

.10.1016/j.jrmge.2016.01.003

118.

Basaran

,

C.

, and

Yan

,

C.-Y.

,

1998

, “

A Thermodynamic Framework for Damage Mechanics of Solder Joints

,”

ASME J. Electron. Packag.

,

120

(

4

), pp.

379

–

384

.10.1115/1.2792650

119.

Basaran

,

C.

, and

Tang

,

H.

,

2002

, “

Implementation of a Thermodynamic Framework for Damage Mechanics of Solder Interconnects in Microelectronic Packaging

,”

Int. J. Damage Mech.

,

11

(

1

), pp.

87

–

108

.10.1106/105678902022259

120.

Basaran

,

C.

, and

Chandaroy

,

R.

,

1998

, “

Mechanics of Pb40/Sn60 Near-Eutectic Solder Alloys Subjected to Vibrations

,”

Appl. Math. Model.

,

22

(

8

), pp.

601

–

627

.10.1016/S0307-904X(98)10059-8

121.

Chandaroy

,

R.

, and

Basaran

,

C.

,

1999

, “

Damage Mechanics of Surface Mount Technology Solder Joints Under Concurrent Thermal and Dynamic Loading

,”

ASME J. Electron. Packag.

,

121

(

2

), pp.

61

–

68

.10.1115/1.2792669

122.

Bin Jamal M

,

N.

,

Kumar

,

A.

,

Lakshmana Rao

,

C.

, and

Basaran

,

C.

,

2019

, “

Low Cycle Fatigue Life Prediction Using Unified Mechanics Theory in Ti–6Al–4V Alloys

,”

Entropy

,

22

(

1

), p.

24

.10.3390/e22010024

123.

Elices

,

M.

,

Guinea

,

G. V.

,

Gómez

,

J.

, and

Planas

,

J.

,

2002

, “

The Cohesive Zone Model: Advantages, Limitations and Challenges

,”

Eng. Fract. Mech.

,

69

(

2

), pp.

138

–

163

.10.1016/S0013-7944(01)00083-2

124.

Bhate

,

D.

,

Chan

,

D.

,

Subbarayan

,

G.

, and

Nguyen

,

L.

,

2008

, “

A Nonlinear Fracture Mechanics Approach to Modeling Fatigue Crack Growth in Solder Joints

,”

ASME J. Electron. Packag.

,

130

(

2

), p.

021003

.10.1115/1.2840057

125.

Benabou

,

L.

,

Sun

,

Z.

, and

Dahoo

,

P. R.

,

2013

, “

A Thermo-Mechanical Cohesive Zone Model for Solder Joint Lifetime Prediction

,”

Int. J. Fatig.

,

49

, pp.

18

–

30

.10.1016/j.ijfatigue.2012.12.008

126.

Miyazaki

,

T.

, and

Ikeda

,

O.

,

2015

, “

Development of Sn-Cu Based Solder for Power Modules

,”

J. Smart Process.

,

4

(

4

), pp.

184

–

189

.10.7791/jspmee.4.184

127.

Miyazaki

,

T.

,

Ikeda

,

O.

,

Kushima

,

T.

, and

Kawase

,

D.

,

2019

, “

Reliability and Failure Modes in Power Cycling Tests for Solder Bonding

,”

Proceedings of 25th Symposium on Microjoining and Assembly Technology in Electronics (Mate 2019)

, Yokohama, Japan, Jan. 29–30, pp.

337

–

340

.

128.

Tanaka

,

Y.

,

Fukumoto

,

A.

,

Endo

,

K.

,

Taya

,

M.

,

Yamazaki

,

K.

, and

Nishikawa

,

K.

,

2016

, “

Evaluation of Vertical Degradation in Lead-Free Solder

,”

Proceedings of 22nd Symposium on Microjoining and Assembly Technology in Electronics (Mate 2016)

, Yokohama, Japan, Feb. 2–3, pp.

135

–

140

.

129.

Sugimoto

,

H.

,

Kariya

,

Y.

,

Hanada

,

R.

,

Ito

,

Y.

,

Yokoyama

,

Y.

, and

Soda

,

S.

,

2019

, “

Fracture Analysis of Vertical Direction Fracture in Die Attach Joint for Power Semiconductor Device

,”

Proceedings of 32nd JSME Computational Division Conference

, Kawagoe, Japan, Sept. 16–18, Paper No. 186.10.23919/LTB-3D.2019.8735334

130.

Harubeppu

,

Y.

,

Tanie

,

H.

,

Ikeda

,

O.

,

Miyazaki

,

T.

,

Kawase

,

D.

,

Morita

,

T.

, and

Sasaki

,

K.

,

2019

, “

Study on Vertical Crack Mechanism of Solder

,”

Proceedings of JSME M&M2019 Conference

, Fukuoka, Japan, Nov. 2–4, Paper No. OS0312.10.1299/jsmemm.2019.OS0312

131.

Taira

,

S.

,

1962

, “

Lifetime of Structures Subjected to Varying Load and Temperature

,”

Creep in Structures

,

Springer-Verlag

, Berlin/Goettingen/Heidelberg, pp.

96

–

124

.

132.

Halford

,

G. R.

,

Hirschberg

,

M. H.

, and

Manson

,

S. S.

,

1973

, “

Temperature Effects on the Strain Range Partitioning Approach for Creep Fatigue Analysis

,”

Fatigue at Elevated Temperatures

, A. E. Carden, A. J. McEvily and C. H. Wells, ed.,

ASTM International

,

West Conshohocken, PA

, pp.

658

–

669

.10.1520/STP38877S

133.

Manson

,

S. S.

,

1973

, “

The Challenge to Unify Treatment of High Temperature Fatigue–a Partisan Proposal Based on Strainrange Partitioning

,”

Fatigue at Elevated Temperatures

, A. E. Carden, A. J. McEvily, and C. H. Wells, ed.,

ASTM International

,

West Conshohocken, PA

, pp.

744

–

782

.10.1520/STP38885S

134.

Robinson

,

E. L.

,

1952

, “

Effect of Temperature Variation on the Long-Time Rupture Strength of Steels

,”

Trans. ASME

,

74

(

5

), pp.

777

–

781

.

135.

Miyazaki

,

N.

,

Aihara

,

Y.

,

Hagihara

,

S.

, and

Munakata

,

T.

,

1991

, “

Finite Element Creep Damage Analysis of a Tube Using Damage Mechanics

,”

Kagaku Kogaku Ronbunshu

,

17

(

1

), pp.

187

–

193

.10.1252/kakoronbunshu.17.187

136.

SEMIKRON

,

2014

, “

Application Note AN1404: Thermal Resistance of IGBT Modules–Specification and Modelling

,” SEMIKRON International GmbH, Nuremberg, Germany.

137.

Schwabe

,

C.

,

Seidel

,

P.

, and

Lutz

,

J.

,

2019

, “

Power Cycling Capability of Silicon Low-Voltage MOSFETs Under Different Operation Conditions

,”

Proceedings of 31th International Symposium on Power Semiconductor Devices and ICs

, Shanghai, China, May 19–23, pp.

495

–

498

.10.1109/ISPSD.2019.8757631

138.

Endo

,

T.

,

1988

, “

Review on Life Prediction for Complex Load Versus Time Histories (Critical Review on Rainflow Algorithm)

,”

Trans. JSME. (Ser. A)

,

54

(

501

), pp.

869

–

874

.10.1299/kikaia.54.869

139.

Khosrovaneh

,

A. K.

, and

Dowling

,

N. E.

,

1990

, “

Fatigue Loading History Reconstruction Based on the Rainflow Technique

,”

Int. J. Fatig.

,

12

(

2

), pp.

99

–

106

.10.1016/0142-1123(90)90679-9

140.

Samavatian

,

V.

,

Iman-Eini

,

H.

, and

Avenas

,

Y.

,

2018

, “

An Efficient Online Time-Temperature-Dependent Creep-Fatigue Rainflow Counting Algorithm

,”

Int. J. Fatig.

,

116

, pp.

284

–

292

.10.1016/j.ijfatigue.2018.06.037

141.

Khatibi

,

G.

,

Lederer

,

M.

,

Weiss

,

B.

,

Licht

,

T.

,

Bernardi

,

J.

, and

Danninger

,

H.

,

2010

, “

Accelerated Mechanical Fatigue Testing and Lifetime of Interconnects in Microelectronics

,”

Procedia Eng.

,

2

(

1

), pp.

511

–

519

.10.1016/j.proeng.2010.03.055

142.

Czerny

,

B.

, and

Khatibi

,

G.

,

2018

, “

Cyclic Robustness of Heavy Wire Bonds: Al, AlMg, Cu and CucorAl

,”

Microelectron. Reliab.

,

88–90,

pp.

745

–

751

.10.1016/j.microrel.2018.07.003

Review of Methodologies for Structural Integrity Evaluation of Power Modules

Abstract

1 Introduction

2 Failure of Bonding Wire

2.1 Heel Cracking.

2.2 Wire-Liftoff

2.2.1 Conventional Failure Models.

2.2.2 Failure Models at Elevated Temperatures.

2.2.2.1 Inelastic strain energy density range $Δ W_{in}$ ⁠.

2.2.2.2 Nonlinear fracture mechanics parameter T* range $Δ T$ *.

2.2.2.3 Application of failure models at elevated temperatures.

3 Die Attach Cracking

3.1 Lifetime Estimation Methodologies Without Consideration of Crack Growth

3.1.1 Lifetime Estimation Formulae Using Temperature Range.

3.1.2 Lifetime Estimation Formula Using Inelastic Strain Range.

3.1.3 Lifetime Estimation Formula Using Inelastic Strain Energy Density Range.

3.2 Lifetime Estimation Methodologies With Consideration of Crack Propagation.

3.2.1 Lifetime Estimation Methodology Using Inelastic Strain Range.

3.2.2 Lifetime Estimation Methodology Using Inelastic Strain Energy Density Range.

3.2.3 Lifetime Estimation Methodology Using Fracture Mechanics Parameters.

3.2.4 Other Lifetime Estimation Methodologies

3.2.4.1 Continuum damage mechanics.

3.2.4.2 Cohesive zone model.

3.3 Lifetime Estimation Methodologies for Die Attach Failure in the Thickness Direction.

4 Power Cycling Test and Thermal Cycling Test

5 Summaries of Previous Research

6 Proposals for Future Research

Funding Data

References

Get Email Alerts

Cited By

Journals

Conference Proceedings

eBooks

Resources

Opportunities

Review of Methodologies for Structural Integrity Evaluation of Power Modules

Abstract

1 Introduction

2 Failure of Bonding Wire

2.1 Heel Cracking.

2.2 Wire-Liftoff

2.2.1 Conventional Failure Models.

2.2.2 Failure Models at Elevated Temperatures.

2.2.2.1 Inelastic strain energy density range ΔWin⁠.

2.2.2.2 Nonlinear fracture mechanics parameter T* range ΔT*.

2.2.2.3 Application of failure models at elevated temperatures.

3 Die Attach Cracking

3.1 Lifetime Estimation Methodologies Without Consideration of Crack Growth

3.1.1 Lifetime Estimation Formulae Using Temperature Range.

3.1.2 Lifetime Estimation Formula Using Inelastic Strain Range.

3.1.3 Lifetime Estimation Formula Using Inelastic Strain Energy Density Range.

3.2 Lifetime Estimation Methodologies With Consideration of Crack Propagation.

3.2.1 Lifetime Estimation Methodology Using Inelastic Strain Range.

3.2.2 Lifetime Estimation Methodology Using Inelastic Strain Energy Density Range.

3.2.3 Lifetime Estimation Methodology Using Fracture Mechanics Parameters.

3.2.4 Other Lifetime Estimation Methodologies

3.2.4.1 Continuum damage mechanics.

3.2.4.2 Cohesive zone model.

3.3 Lifetime Estimation Methodologies for Die Attach Failure in the Thickness Direction.

4 Power Cycling Test and Thermal Cycling Test

5 Summaries of Previous Research

6 Proposals for Future Research

Funding Data

References

Get Email Alerts

Cited By

Related Articles

Related Proceedings Papers

Related Chapters

Journals

Conference Proceedings

eBooks

Resources

Opportunities

This Feature Is Available To Subscribers Only

2.2.2.1 Inelastic strain energy density range $Δ W_{in}$ ⁠.

2.2.2.2 Nonlinear fracture mechanics parameter T* range $Δ T$ *.