SciELO - Scientific Electronic Library Online

 
vol.14 número13Recent Developments in Damage Identification of Structures Using Data MiningAssessing the Notch Sensitivity of High Strength Wires índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

Compartilhar


Latin American Journal of Solids and Structures

versão impressa ISSN 1679-7817versão On-line ISSN 1679-7825

Lat. Am. j. solids struct. vol.14 no.13 Rio de Janeiro  2017

http://dx.doi.org/10.1590/1679-78253163 

Articles

An Enhanced Random Vibration and Fatigue Model for Printed Circuit Boards

Bruno de Castro Braza  * 

Flávio Luiz de Silva Bussamrab 

a Space System Division, Instituto Nacional de Pesquisas Espaciais, São José dos Campos, Brazil. Email: bruno.braz@inpe.br

b Aerospace Engineering Division, Instituto Tecnológico de Aeronáutica, São José dos Campos, Brazil. Email: flaviobu@ita.br

Abstract

Aerospace vehicles are mostly exposed to random vibration loads during its operational lifetime. These harsh conditions excites vibration responses in the vehicles printed circuit boards, what can cause failure on mission functionality due to fatigue damage of electronic components. A novel analytical model to evaluate the useful life of embedded electronic components (capacitors, chips, oscillators etc.) mounted on Printed Circuit Boards (PCB) is presented. The fatigue damage predictions are calculated by the relative displacement between the PCB and the component, the lead stiffness, as well the natural vibration modes of the PCB and the component itself. Statistical methods are used for fatigue cycle counting. The model is applied to experimental fatigue tests of PCBs available on literature. The analytical results are of the same magnitude order of the experimental findings.

Keywords: Random vibration; fatigue; PCB

1 INTRODUCTION

Aerospace vehicles (airplanes, launcher vehicles and satellites) are exposed to harsh random vibrating environments. These mechanical stresses occur during qualification and acceptance tests, launching and orbit injection (spacecraft) and during continuous flight operations (airplanes). The aerospace vehicles shall be designed to withstand these conditions, providing a safe environment to the onboard equipment. Electronic equipment is constituted mainly of Printed Circuit Boards (PCB) and a support structure, which shall ensure a proper mechanical, thermal and electrical environment to the embedded electronic components (EC), as capacitors, resistors, chips, oscillators etc. The design of aerospace equipment shall be capable to predict the time to failure of the electronic components. Random vibration is the common environment of aerospace vehicles, then, fatigue calculations are treated statistically.

Wong et al. (2000), Wu (2009) and Grieu et al. (2008) presented stress calculations at PCBs and ECs based on accurate Finite Element Method (FEM) modeling. The methodology consisted of meshing the entire PCB and the electronic component parts, as leads and body. This approach seems to be too time consuming when applied to PCBs with several components. Cifuentes (1994), Yang et al. (2000), Amy et al. (2010, 2006) and Sayles and Stoumbos (2015) proposed simplified FEM models for calculation of PCB’s eigenvalues and eigenvectors. The mounted components on the board were modeled by updating the mass and stiffness properties of the PCB.

Experimental vibration tests were presented by Singal et al. (1992), Amy et al. (2010, 2009) and Yang et al. (2002). Veprik and Barbitsky (2000), Veprik (2003) and Esser and Huston (2003) presented experimental data for damped systems. The experimental approach validates only the tested or similar setups, so additional tests are necessary for different configurations.

Silva and Gonçalves (2013) proposed an analytical model for lead stress prediction. The leads were modeled as beams, the EC as a 6-degree-of-freedom rigid bodies and the PCB as a simply supported plate. These authors did not apply a fatigue model for component life evaluation. Steinberg (2000) proposed a semi-analytical method for estimating the useful life of electronic components mounted on PCBs under random vibrations. The PCB was modeled as a beam. The beam displacements were taken into account, but the rotations not. The fatigue damage prediction was performed for an electronic component (EC) in the center of the board, and the results were approximated elsewhere. Steinberg’s method is one of the methods most used in aerospace industry.

The objective of this work is to present a novel accurate analytical model with low computational cost for fatigue life prediction of electronic components mounted on PCBs, which typically fails at the leads. The loads at the leads are calculated by the relative displacements and rotations between the PCB and the EC when the system is excited. The PCB is modeled as a thin plate, considering a uniform smearing of electrical component masses. The EC is modeled as: a) a rigid body and b) a system of two perpendicular beams (3-degree-of-freedom system). The fatigue life is evaluated by Dirlik (1985) statistical fatigue damage model.

The present work enhances the well-known semi-analytical method developed by Steinberg (2000) by proposing a plate model (instead of a beam) to the PCB, adding the contributions of the PCB rotations, the EC body flexibility and the EC first natural frequency. A similar model seems not to be available in the literature. Also, the present work adopt the Dirlik’s random fatigue damage counting method (instead of the 3-band technique), one of the methods most used in industry for fatigue calculations. Steinberg’s method compensates the lack of these features by the use of the experimental constant C st , defined in his book for each type of electronic component. The Steinberg’s method is given by

N=20106{(0.0002225.41.5BPCBCsttPCBrL)1Zmax_rms}bb, (1)

where, N is the cycles to failure, B PCB is the PCB length, L is the EC length parallel to B PCB , r is a factor given by a sine function to account the EC position on the PCB (r = 1 in the PCB’s center, r < 1 in the other regions), C st is a constant for different types of ECs, t PCB is the PCB thickness and Z max_ rms is the root mean square of the PCB maximum displacement, all given in millimeters.

2 ELECTRONIC COMPONENT LEADS FATIGUE

Electronic components can be simply approximated by three different parts. The main part is the component body, usually covered by a case made of plastic or ceramic. The leads connect the component to the PCB circuitry. The third part is the solder joints, which attach the leads to the board. The typical failure of electronic components under vibration is the fatigue of leads and/or solder joints due to the relative displacements between the component and the PCB.

The lead is assumed as a beam with bending, axial and torsional stiffness. The lead geometry is modeled with the polygonal form exhibited in Fig. 1. For simpler cases some quotas can be set to zero.

Figure 1: Simplified (polygonal) lead geometry. 

The lead fatigue life model is presented below. Section 2.1 describes the lead stiffness calculations, followed by the lead stretch model definition on Section 2.2. The formulation is extended to account for the component body flexibility (Sections 2.3 and 2.4) and PCB bending rotations (Sections 2.5 and 2.6). At last, the stresses (Section 2.7) and the fatigue damage (Section 2.8) complete the model.

2.1 Lead Stiffness

Figure 2a shows the free-body diagram of a lead when the component moves perpendicular to the PCB plane (lead longitudinal direction). The model considers that the leads are clamped on the component body (extremity 1), which is so far adopted rigid, and the deformed PCB applies an external force P at extremity 2 (solder joint). Figures 2b and 2c present the free-body diagrams when the PCB curvature introduces bending and torsional moments, M and T, respectively, at point 2.

Figure 2: Lead free-body diagrams: (a) concentrated lead P; (b) bending moment M; (c) torsional moment T

Based on Euler-Bernoulli beam theory, the lead stiffness in the lead longitudinal direction due to a force P is given by

kl={(c33EIz'2+b2c2EIz'2+b1c2EIz'1+(a3+3a2c+3ac2)3EIz'1(a+c)2dEIz'1)+ (b1+dEA1+b2EA2)}1, (2)

where E is the elasticity modulus of the lead, I z’i and A i are respectively the moment of inertia parallel to z’ axis taken at the section geometric center and the cross section area of the part i, i =1,2 (b 1 , b 2 ), as depicted in Fig. 1 (lateral view 2).

Similarly, the bending and torsional stiffness are given by, respectively,

kb=(cEIz'2+b2EIz'2+b1EIz'1+aEIz'1dEIz'1)1, (3)

kt=(cGJ2+b2EIx'2+b1EIx'1+aGJ1dEIx'1)1, (4)

where G is the shear modulus of elasticity, J i and I x’i are respectively the polar moment of inertia and the moment of inertia parallel to x’ axis taken at the section geometric center of the part i, i =1,2 (b 1 , b 2 ).

2.2 Lead Stretch

The lead stretches are due to the PCB displacements and the EC acceleration. The PCB is modeled as a simply supported rectangular plate of dimensions ā × b¯ . The x- y coordinate system origin is located on the PCB left lower corner, and the x c - y c EC local coordinate system is located on component center, as depicted in Fig. 3.

Figure 3: PCB with an electronic component (EC). 

The leads local coordinate system x’- y’ is perpendicular to the plane x- y of the PCB: x’ lays on the x- y plane and y’ is perpendicular to this plane. The W mn natural frequencies of the PCB and its Z mn displacements are given by (Blevins 2001)

Wmn=π2[(ma¯)2+(nb¯)2]EPCBtPCB312ρPCB(1υPCB2), (5)

Zmn(x,y,W)=2Bmn(W)ρa¯b¯sin(mπxa¯)sin(nπyb¯), (6)

where m and n are the number of half waves in x and y directions, respectively. E PCB , t PCB , ρ PCB and υ PCB are the PCB elasticity modulus, thickness, density and Poisson’s ratio. The PCB density ρ PCB is the ratio between the total mass of the system (PCB plus ECs) and the volume of the PCB. B mn (W) determines the vibration magnitude. For a uniform acceleration distribution α over the PCB, B mn is given by

Bmn(W)=Gmn(W)α(W), (7)

where

Gmn(W)=2ρabπ2mn(1cos(πm))(1cos(πn))Wmn2W2+2iξWmnW. (8)

So far, it is assumed that the PCB and the leads are both flexible, all the displacements are small (senθ ~ θ), and the electronic component body is a rigid plate. Then the transverse (lead longitudinal) displacements can be written as

Zc,mn(x,y)=Z0+(xe¯)θc+(yf¯)φc,e¯c¯2xe¯+c¯2 , f¯d¯2yf¯+d¯2 (9)

where Z 0 stands for displacement in z direction at the center of the EC (x=e¯,y=f¯) and θ c (or ϕ c ) is the rotation around x (or y) axis, as displayed in Fig. 4. These values are determined considering a uniform distribution of springs connecting the component to the board and the equilibrium equations (details at Appendix A).

Figure 4: Kinematics of the PCB and EC: (a) underformed position; (b) displacement and rotations. 

Now, it is possible to determine the δ cg,mn lead stretch in the position (x,y) and in the W frequency due to the PCB vibration as

δcg,mn(x,y,W)=Zc,mnZmn=δmn(x,y,W)Bmn, (10)

where δ mn is a transfer function from B mn to δ cg_mn .

In addition to the stretches due to the PCB displacement, the acceleration itself induces a mass/spring/damper behavior on the component, modeled as a 1-degree-of-freedom system with total stiffness calculated as the sum of all n leads k l (Eq. 2). The component natural frequency, vibration transfer function and stretch due to the acceleration are given respectively by

Wc,long=nklmc, (11)

Gc,long(W)=1Wc,long2W2+2iξc,longWc,longW, (12)

δca,mn(x,y,W)=[Zmn(x,y)Gmn(W)(W2)+1]Gc,long(W)α(W), (13)

where m c is the total mass of the electronic component (EC), ξ c,long is the lead damper coefficient, adopted as 0.05 (as usual for some metal alloys) and i is the complex number.

2.3 Component Body Flexibility to Displacements

In the previous Section, the geometric stretch δ cg,mn and the acceleration stretch δ ca,mn are calculated assuming the component to be a rigid body. However, the component body is flexible and absorbs part of the strain, what relieves the stresses on the leads. The component flexibility calculation is not a trivial task, because the lead itself is a flexible structure and the loads are dependent on the component position and on the PCB vibration mode. The EC is supposed to be a rectangle of side lengths c¯×d¯ . Each lead, located on the EC boundary at distances l x and l y from the EC center, applies a concentrated load P, as shown in Fig. 5. The EC is modeled as a sum of two distinct beams clamped in the component center, where R x and R y are the fractions of the load applied on each direction:

Rx=lxlx+ly,Ry=lylx+ly. (14)

Figure 5: Double-beam model for EC. 

Therefore, the total component body displacement (δ c ) at the lead is the sum of the displacements due to the R x P force (δ c,x ), and due to the R y P force (δ c,y ), say

δc=δc,x+δc,y, (15a)

δc,x=c¯3RxP24EIx,c,δc,y=d¯3RyP24EIy,c, (15b)

where E is the elasticity modulus, I x,c and I y,c are the inertia moment around x c and y c axes (Fig. 5) and t c the component body thickness. The EC inertia moments are calculated considering an internal void of t c /2:

Ix,c=d¯[tc3(tc/2)3]12,Iy,c=c¯[tc3(tc/2)3]12. (16)

The force applied on the component body equals to the force applied on the leads, and the total stretch is the sum of the component stretch and the lead stretch (δ lead,x , δ lead,y ). For equilibrium of EC

RxP=24EIx,cc¯3δc,x=klnxδlead,x (17)

then

δc,x=c¯3klnx24EIx,cδlead,x,δc,y=d¯3klny24EIy,cδlead,y (18)

where n x and n y are the number of leads on each side of length d¯ and c¯ , respectively. Thus, the total stretch δ can be calculated as the sum of leads and component body stretches

δ=(c¯3klnx24EIx,c+1)δlead,x+(d¯3klny24EIy,c+1)δlead,y (19)

Assuming a linear relation between lead stretches δ lead,x and δ lead,y , with R x + R y = 1, the total stretch and consequently the fat c lead stretch reduction factor can be written as

δlead,x=Rxδlead,δlead,y=Ryδlead (20)

δ=fatc(x,y)1δlead, (21)

where

fatc(x,y)=(Rxc¯3klnx24EIx,c+Ryd¯3klny24EIy,c+1)1. (22)

2.4 Lead Total Stretch

The δ c,long total lead stretch can be calculated for all PCB vibration modes from Eqs. 10, 13 and 22

δc,long=m,n=1,1m,n=φ,φ[δcg,mn+δca,mn]fatc=m,n=1,1m,n=φ,φ{δmn(x,y,W)Gmn(W) +[Zmn(x,y)Gmn(W)(W2)+1φ2]Gc,long(W)}fatc(x,y)α(W). (23)

Note that the term 1 is divided by φ 2 in order to not account φ 2 times the α input acceleration.

2.5 Lead Rotation

The lead rotations are caused by the PCB bending curvature during the vibrations (θ PCB or ϕ PCB ) minus the component body rotation (θ c or ϕ c ) and minus the rotation due to the lead longitudinal stretch (θ 2 - see Appendix A). The PCB rotations, which induce torsional and bending moments, are given by the derivative of Eq. 6 by x and y. The rotations necessary to describe the problem are depicted in Fig. 4 and summarized at Appendix A.

2.6 Component Body Flexibility to Rotations

The component body flexibility to rotations and moments shall also be applied to the model. The flexibility factor due to rotations differs from the flexibility factor due to forces by the difference in beam stiffness for these two distinct types of loads. The cantilever beam stiffness to moments is

kbeam,moment=EIL, (24)

then, the flexibility factors for bending and torsional moments assume a slightly different format

fatcm1(x,y)=(Rx(x,y)c¯kbnx2EIx,c+Ry(x,y)d¯kbny2EIy,c+1)1fatc(x,y), (25)

fatcm2(x,y)=(Rx(x,y)c¯ktnx2EIx,c+Ry(x,y)d¯ktny2EIy,c+1)1fatc(x,y), (26)

where fat c is the flexibility reduction factor for force (Eq. 22), fat cm1 is the flexibility reduction factor for lead bending rotations (Fig. 2b), and fat cm2 is the flexibility reduction factor for lead torsional rotations (Fig. 2c). In the same way fat c multiplies the lead total stretch (Eq. 23), fat cm1 multiplies the lead total bending rotation around y and x, and fat cm2 multiplies the lead total torsional rotation around y and x, as presented on Appendix A (Eqs A.8-11). The inclusion of fat c in Eqs. 25 and 26 is a simplified way to account for component body additional rotation due to the longitudinal stretches.

2.7 Lead Stress

The stretches and rotations found on the previous Sections are now used to evaluate the lead stress. The von Mises stress is directly proportional to the input acceleration

σVM(x,y,W)=GVM(x,y,W)α(W), (27)

where G VM is the transfer function from input acceleration α to the von Mises stress σ VM . For random vibrations, the variables are described in Power Spectral Density (PSD) (Lalanne 2002), resulting in

σVM,PSD(x,y,W)=|GVM(x,y,W)|2αPSD(W). (28)

2.8 Fatigue Damage

The fatigue damage model used on this work is the well-known Miner’s rule (Lalanne 2002), where the fatigue damage D is calculated as the relation between the number n i of performed cycles in a prescribed stress, to the number N i of cycles necessary to produce a failure in the same stress level, as described in Eq. 29. The failure occurs when n i =N i .

D=niNi. (29)

The fatigue damage in random vibrations, once obtained the stress response in PSD (Eq. 28), is calculated in this work by Dirlik’s method (Dirlik 1985), one of the most used in industry and literature. The damage, evaluated with a probability distribution function developed statistically through the use of several random signals and with the Rainflow cycle counting method (ASTM 2011), is calculated as

D=np+T(σVM,rms)bbCb[G1QbbΓ(1+b)+(2)bbΓ(1+bb2)(G2|R|bb+G3)], (30)

where T is the time, σ VM,rms is the root mean square of the von Mises stress, Γ is the Gamma function, b b and C b are the S-N variables given by Basquin relation (Nσbb=Cb) , and G 1 , G 2 , G 3 , R and Q are variables defined on Appendix B.

All stresses are directly proportional to the input accelerations. It is possible to affirm that for each PCB vibration mode, the relations between shear and normal stresses are constant. When describing the stresses in PSD, these relations are always constant (in statistical sense) for each frequency band, apart the PCB vibration modes, because all vibration modes are summed (Eq. 23). Finally, as the root mean square of the power spectral densities represents a statistical average of the stresses, it is possible to assume that the relations between shear and normal stresses are in average constant. Therefore, the von Mises root mean square (rms) stress can be used to approximate the fatigue damage calculations in random vibrations, as presented on Eq. 30.

3 RESULTS

The present model is applied to PCB’s experimental fatigue tests available on the literature. The results are also compared to Steinberg’s method predictions. The experimental results findings discussed in this Section include most of the types of leads, assemblies and solder joints. All materials, dimensions and properties used for each case are presented on Appendix C.

In each system, the model is applied with four distinct configurations in order to clarify the fatigue damage contributors:

  1. Complete / 9 modes: the PCB is modeled as a plate with 9 natural vibration modes, and the ECs as double beams with 1 natural mode;

  2. Acceleration off: the PCB is modeled as a plate with 9 modes, and the ECs as double beams without the contribution of its natural frequency, expressed by δ ca,mn (Eq. 13);

  3. Moment off: the PCB is modeled as a plate with 9 natural modes, and the ECs as double beams with 1 natural mode and not accounting the influence of PCB curvature;

  4. Complete / 1 mode: the PCB is modeled as a plate with only 1 natural mode, and the ECs as double beams with 1 natural mode.

3.1 Plastic Leaded Chip Carrier (PLCC) and Leadless Ceramic Chip Carrier (LCCC)

Liguore and Followell (1995) performed random vibration tests in 8 PCBs with 9 electronic components, 3 of the type PLCC with 68 J-leads and the others of the type LCCC with 32, 68 and 84 pins, as described in Fig. 6. The epoxy fiberglass PCBs (E-Glass) with 152.4 × 152.4 × 1.45 mm3 were mounted in a shaker with the side close to the components U1, U2 and U3 free. As the analytical model uses a 4-side simply supported plate, the PCB was mirrored about the free edge. The PCBs were vibrated in a band frequency of 50 Hz around the first natural frequency of the board with a constant PSD. The applied levels and PCB’s amplifications varied in order to produce results between 30 to 80 Grms.

Figure 6: Test setup - PLCC and LCCC (Liguore and Followell 1995). 

In the present model, the PLCC leads of copper are modeled with the following geometric parameters: a = 0.08 mm, b 1 = 1.02 mm, b 2 = 1.02 mm, t = 0.25 mm, l 1 = 0.735 mm and l 2 = 0.43 mm. The leads made of Pb solder (LCCC) have more or less conical shape and the rupture often appears in the middle of the cone. Thus the LCCC lead is modeled with the b 1 = b 2 = 0.85/2 mm parts and with a constant diameter given by the average (l 1 = l 2 = 1.27 mm).

The results of the present model are exposed in Tables 1 to 5 and compared with experimental results from Liguore and Followell (1995) and semi-analytical beam model from Steinberg (2000). Liguore and Followell (1995) presented their results only for component type. No information about component position was given. So, the results of components U1, U2 and U3 are analyzed altogether. All the present model results are of the same magnitude order of the experimental ones.

Table 1: Cycles to failure for PLCC with 30.6 Grms, 68 J-Lead. 

Model U1 and U3 components U2 component Fail Point (Fig. 1)
Complete / 9 modes 4.79×106 14.73×106 a
Acceleration off 5.14×106 16.47×106 a
Moment off 5.22×106 16.00×106 a
Complete / 1 mode 4.78×106 14.73×106 a
Steinberg (2000) 1.96×106 0.21×106 not available
Liguore and Followell (1995) 1×106 to 9×106 not available

Table 2: Cycles to failure for PLCC with 80.7 Grms, 68 J-Lead. 

Model U1 and U3 components U2 component Fail Point (Fig. 1)
Complete / 9 modes 0.86×104 2.65×104 a
Acceleration off 0.92×104 2.96×104 a
Moment off 0.94×104 2.87×104 a
Complete / 1 mode 0.86×104 2.65×104 a
Steinberg (2000) 0.37×104 0.038×104 not available
Liguore and Followell (1995) 2×104 to 10×104 not available

Table 3: Cycles to failure for LCCC with 30.6 Grms, 32 pins. 

Model U5 and U6 components Fail Point (Fig. 1)
Complete / 9 modes 5.19×105 solder
Acceleration off 5.44×105 solder
Moment off 11.52×105 solder
Complete / 1 mode 5.19×105 solder
Steinberg (2000) 0.90×105 not available
Liguore and Followell (1995) 10×105 to 40×105 solder

Table 4: Cycles to failure for LCCC with 30.6 Grms, 68 pins. 

Model U4 and U7 components Fail Point (Fig. 1)
Complete / 9 modes 1.87×105 solder
Acceleration off 1.90×105 solder
Moment off 2.95×105 solder
Complete / 1 mode 1.87×105 solder
Steinberg (2000) 1.17×105 not available
Liguore and Followell (1995) 2×105 to 20×105 solder

Table 5: Cycles to failure for LCCC with 30.6 Grms, 84 pins. 

Model U8 and U9 components Fail Point (Fig. 1)
Complete / 9 modes 5.94×105 solder
Acceleration off 6.00×105 solder
Moment off 8.52×105 solder
Complete / 1 mode 5.95×105 solder
Steinberg (2000) 0.63×105 not available
Liguore and Followell (1995) 2.0×105 to 10×105 solder

The influences of acceleration stretch are lower than 11% in the fatigue life, explained by the ECs calculated natural frequencies 200 times higher than the maximum vibration frequency of 175 Hz. The component with the lower calculated natural frequency is the 68 J-lead (U1, U2 and U3) with 31123 Hz. Also, as the geometrical stretch (Eq. 10) is the main damage contributor for this case, the regions in the diagonal of the board have lower fatigue lives.

The PLCCs have longer fatigue lives than the LCCCs due to the difference in stiffness. Both components have close fat c factors (Eq. 22). However the LCCC components have ceramic bodies and short leads made of solder, which are stiffer and produce more stresses.

Steinberg’s results for all EC but U2 are of the same magnitude order of experimental results since the geometrical stretch is the most damaging contributor in these components. However, it is conservative for the EC in the PCB center region (U2).

3.2 Plastic Ball Grid Array (PBGA)

The Plastic Ball Grid Array (PBGA) components are constituted by plastic body and spherical leads of solder. YU et al (2011) presented random vibration results for PBGA with SAC 305 and SAC 405 Pb free solder joints (lead). The tests were performed in FR4 PCBs of dimensions 100 × 50 × 0.6 mm3 attached to the shaker through 6 screws, as described in Fig. 7. Two tests were performed for each solder configuration, one with 0.15 g2Hz-1, another with 0.25 g2Hz-1 (PSD), both in 40 to 1000 Hz frequency range.

Figure 7: Test setup - PBGA (Yu et al. 2011). 

The present model results are obtained by assuming the leads as a cylinder of 2 parts, b 1 and b 2 . The first part has the average diameter of l 1 = 0.266 mm and the total length of the solder joint/lead b 1 = 0.2 mm. The second part has a negligible length b 2 = 0 and diameter of l 2 = 0.22 mm (the smallest diameter of the solder ball), used for stress calculations.

The results are exposed in Tables 6 to 9. As expected from solder leads, the high stiffness of the connection makes the geometric stretch (Eq. 10) the main contributor to the fatigue damage. The components calculated natural frequency is 606320 Hz, which is far beyond the tested frequency band. Therefore, the acceleration stretch contribution was minimized. It is also possible to observe the good adherence of the present predictions.

Table 6: Minutes to failure for PBGA SAC 305 with 0.25 g2Hz-1

Model x= ā/4; y= b¯ /2 (Fig. 3) Fail Point (Fig. 1)
Complete / 9 modes 1.27×102 solder
Acceleration off 1.33×102 solder
Moment off 1.71×102 solder
Complete / 1 mode 1.35×102 solder
Steinberg (2000) 0.15×102 not available
YU et al (2011) 0.1×102 to 2.20×102 solder
average: 1.13×102

Table 7: Minutes to failure for PBGA SAC 305 with 0.15 g2Hz-1

Model x= ā/4; y= b¯ /2 (Fig. 3) Fail Point (Fig. 1)
Complete / 9 modes 7.61×102 solder
Acceleration off 7.99×102 solder
Moment off 10.24×102 solder
Complete / 1 mode 8.08×102 solder
Steinberg (2000) 0.92×102 not available
YU et al (2011) 1.0×102 to 4.0×102 solder
average: 3.95×102

Table 8: Minutes to failure for PBGA SAC 405 with 0.25 g2Hz-1

Model x= ā/4; y= b¯ /2 (Fig. 3) Fail Point (Fig. 1)
Complete / 9 modes 1.41×102 solder
Acceleration off 1.46×102 solder
Moment off 1.72×102 solder
Complete / 1 mode 1.48×102 solder
Steinberg (2000) 0.62×102 not available
YU et al (2011) average: 1.02×102 solder

Table 9: Minutes to failure for PBGA SAC 405 with 0.15 g2Hz-1

Model x= ā/4; y= b¯ /2 (Fig. 3) Fail Point (Fig. 1)
Complete / 9 modes 4.77×102 solder
Acceleration off 4.93×102 solder
Moment off 5.80×102 solder
Complete / 1 mode 4.99×102 solder
Steinberg (2000) 2.09×102 not available
YU et al (2011) average: 2.82×102 solder

In this case, Steinberg’s method presents results in conformance due to the high contribution of the geometrical stretch, the only structural feature modeled by Steinberg.

3.3 Plastic Dual in-line Package (PDIP) and Tantalum Capacitor

The Plastic Dual in-line Package (PDIP) components are composed of a plastic body with 14 leads mounted on 2 parallel sides. The Tantalum capacitors are cylindrical components with 2 leads. Genç (2006) tested 6 PCBs, 3 fully populated with 24 PDIPs (Fig. 8a) and other 3 fully populated with 30 capacitors (Fig. 8a). Both components were mounted on the boards by through hole solder joints. All the PCBs are of 233 × 15 × 1 mm3 dimensions, attached to the shaker as shown on Fig. 8.

Figure 8: Test setup (Genç 2006): a) PDIP; b) Tantalum Capacitor. 

The geometric parameters of the PDIP’s leads of CDA 194 (copper) are: b 1 = 1.715 mm; b 2 = 0.551 mm; l 1 = 1.65 mm; l 2 = 0.5 mm and t = 0.304 mm. The model of the capacitor’s leads of nickel has: a = 3.5 mm, b = 4.51 mm and diameter of 0.64 mm.

Both tests were executed through vibration steps of 60 minutes each in the frequency range of 20 to 2000 Hz with constant PSD. Step 1 level is 2.02×10-3 g2Hz-1 and step 14 level is 6.694×10-3 g2Hz-1. The PSD ratio between steps is close to 1.56. Genç (2006) has not achieved failures for the PDIDs. For these components, the test started at step 3 and lasted until 2.5 minutes of step 14, when the shaker aborted. Then, any result higher than 662.5 minutes for the PDIPs lifetime is in agreement with the experimental tests. The capacitors test started at step 1 and failed mainly in steps 3 and 4.

The comparisons between experimental and analytical times to failure are described in Tables 10 and 11. The results presented by Genç (2006) were not discriminated by EC positions on the board. Therefore, it was only possible to compare the ranges of times to failure.

Table 10: Failure step (minutes to failure) for PDIP. 

Model x= ā/2; y= b¯ /2 (Fig. 3) x= ā/10; y= b¯ /10 (Fig. 3) x= ā/4; y= b¯ /4 (Fig. 3) Fail Point (Fig. 1)
Complete / 9 modes >step 14 (8.43×102) step 14 (6.64×102) step 14 (7.05×102) no fail
Acceleration off >step 14 (8.70×102) step 14 (6.64×102) step 14 (7.06×102) no fail
Moment off >step 14 (8.60×102) step 14 (6.64×102) step 14 (7.08×102) no fail
Complete / 1 mode >step 14 (9.65×102) step 14 (6.67×102) step 14 (7.12×102) no fail
Steinberg (2000) step 10 (5.43×102) step 13 (7.27×102) step 11 (6.31×102) not available
Genç (2006) step 3 up to 2.5 min. of step 14 with no failure no fail
total time: 6.62×102

Table 11: Failure step (minutes to failure) for tantalum capacitor. 

Model x= ā/2; y= b¯ /2 (Fig. 3) x= ā/10; y= b¯ /10 (Fig. 3) x= ā/4; y= b¯ /4 (Fig. 3) Fail Point (Fig. 1)
Complete / 9 modes step 3 (1.52×102) step 6 (3.07×102) step 3 (1.59×102) a
Acceleration off step 9 (5.22×102) >14 step step 12 (6.62×102) a
Moment off step 2 (1.01×102) step 5 (2.53×102) step >14 (9.90×102) a
Complete / 1 mode step 4 (2.07×102) step 4 (1.89×102) step 7 (4.06×102) a
Steinberg (2000) step 11 (5.17×102) step 13 (7.20×102) step 11 (6.12×102) not available
Genç (2006) first fail: between step 3 and 4 (1.52×102 and 2.04×102) a
eleventh fail: step 6 (3.50×102)

The present analyses performed step 14 without a time constraint, so it could surpass 60 minutes. The PDIPs calculated natural frequency of 55730 Hz is much higher than the vibrated band of 20 to 2000Hz. Thus, the acceleration stretch almost did not contribute in fatigue. The capacitors calculated natural frequency of 990 Hz is inside the vibrated band of 20 to 2000Hz. Therefore, the acceleration stretch (Eq. 13) is the most damaging contributor, and the middle region of the PCB is the most critical one for the capacitors. The present results are of the same magnitude order of the experimental ones.

Steinberg’s method showed to be conservative for the PDIDs, and optimistic for the capacitors. In the case of the capacitors, Steinberg’s method failed in describing the time to failure by not considering the EC natural frequency. This is evident when comparing the similarity between the acceleration off model to Steinberg’s results.

4 CONCLUSIONS

A novel analytical model is presented to predict the useful life of electronic components under random vibration. The complete/ 9 modes model presented the closest approximation to the experimental results. As this is the most complete model, and seemed to be always more conservative than the others, this is the recommended configuration. The acceleration off model excludes the acceleration stretch term (δ ca,mn ). Its fatigue lives estimation varied from 0%, for stiff lead attachments, to 316% higher than Complete/9 modes model, for the most flexible one (Tantalum). The moment off configuration proved the importance of the PCB bending effects, since it presented in average life estimation values 35% higher than the recommended model, and the increase in fatigue life estimation varied between 0% and 522%. The complete/ 1 mode model results depend on the test frequency range: for a large range, this model is 4% to 14% higher in life to failure estimation than the recommended model, with a maximum of 155% for the Tantalum capacitor.

The lifetime predictions for the full present model are of the same magnitude order of the experimental findings. The obtained results are more accurate than the Steinberg’s method and could predict reasonable results for all electronic components studied. Steinberg’s method failed mainly in components where EC’s natural frequency is relatively low (Tantalum Capacitors). The proposed method calculates the distributed stress along the leads, what allows a workaround in case of fatigue problems.

The objective of creating an accurate model with low computational cost has been achieved. The practical applications show that the model is suitable for design of aerospace vehicle electronic embedded systems, or any other vehicle exposed to dynamic loads.

References

Amy, R. A., Aglietti, G. S., Richardson G., (2010). Accuracy of simplified printed circuit board finite element models, Elsevier Microelectronics Reliability 50: 86-97. [ Links ]

Amy, A., Aglietti, G. S., Richardson, G., (2010). Board-level vibration failure criteria for printed circuit assemblies: An experimental approach, IEEE Transactions on electronics packaging manufacturing 33, No. 4: 303-311. [ Links ]

Amy, A., Aglietti, G. S., Richardson, G., (2009). Reliability analysis of electronic equipment subjected to shock and vibration - A review, Shock and Vibration 16: 45-59. [ Links ]

Amy, R. A., Aglietti, G. S., Richardson, G., (2006). Simplified modeling of printed circuit boards for spacecraft applications, 57th International Astronautical Congress. [ Links ]

ASTM (American Society for testing and materials) (2011). E1049-85: Standard practices for cycle counting in fatigue analysis. [ Links ]

Basaran, C., Jiang, J., (2002). Measuring intrinsic elastic modulus of Pb/Sn solder alloys, Mechanics of Materials 34: 349-362. [ Links ]

Blevins, R. D., (2001). Formulas for natural frequency and mode shape, Krieger Publishing Company (New York). [ Links ]

Cifuentes, A. O., (1994). Estimating the dynamic behavior of printed circuit boards, IEEE Transactions on components, packaging, and manufacturing technology-part B: Advanced packaging 17: 69-75. [ Links ]

Dirlik, T., (1985). Application of computers in fatigue analysis, Ph.D. Thesis, University of Warwick, Coventry, UK. [ Links ]

Esser, B., Huston, D., (2003). Active mass damping of electronic circuit boards, Elsevier Journal of Sound and Vibration. [ Links ]

Fairchild Semiconductor (2001). 68-Lead plastic leaded chip carrier (PLCC): package dimensions. [ Links ]

Genç, C., (2006). Mechanical fatigue and life estimation analysis of printed circuit board components, Thesis of Science Master Degree, Middle East Technical University, Beirut, Lebanon. [ Links ]

Grieu, M., Massiot, G., Maire, O., Chaillot, A., Bienvenu, C. M. Y., Renard, J., (2008). Durability Modeling of a BGA Component under Random Vibration, Mechanical and Multiphysics Simulation and Experiments in Micro-Electronics and Micro-Systems, IEEE 9th. Int. Conf. on Thermal, EuroSimE. [ Links ]

Lalanne, C., (2002). Mechanical vibration & shock: Random vibration, Hermes Penton Ltd (London). [ Links ]

Liguore, S., Followell, D., (1995). Vibration fatigue of surface mount technology (SMT) solder joints, IEEE Proceedings Annual Reliability and Maintainability Symposium. [ Links ]

Mrnik, M., Slavic, J., Boltezar, M., (2012). Frequency-domain methods for a vibration-fatigue-life estimation - Application to real data, International Journal of Fatigue 47: 8-17. [ Links ]

Sayles, A. A., Stoumbos, T., (2015). Evaluating the dynamic behavior and analytically predicted displacements of printed circuit boards (PCBs) using the “smeared-mass” & fine mesh approach, 56th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference. [ Links ]

Singal, R. K., Gorman, D. J., Forgues, S. A., (1992). A comprehensive analytical solution for free vibration of rectangular plates with classical edge conditions: Experimental Verification, Experimental Mechanics: 21-23. [ Links ]

Silva, G. H. C., Gonçalves, P. J. P., (2013). A model for computing vibration induced stresses of electronic components in a general flexible mounting, Elsevier Journal of Sound and Vibration 332: 5192-5206. [ Links ]

Steinberg, D. S., (2000). Vibration analysis for electronic equipment, John Wiley Sons (New York). [ Links ]

Sumitomo Bakelite (2000). Environmentally friendly epoxy molding compound for semiconductor: Sumikon EME-G series. [ Links ]

Texas Instruments (2001). SMJ320C25, SMJ320C25-50: digital signal processor. [ Links ]

Veprik, A. M., Barbitsky, V. I., (2000). Vibration protection of sensitive electronic equipment from harsh harmonic vibration, Elsevier Journal of Sound and Vibration 238: 19-30. [ Links ]

Veprik, A. M., (2003). Vibration protection of critical components of electronic equipment in harsh environmental conditions, Elsevier Journal of Sound and Vibration 259: 161-175. [ Links ]

Wong, T. E., Reed, B. A., Cohen, H. M., Chu, D. W., (2000). Development of BGA solder joint vibration fatigue life prediction model, IEEE. [ Links ]

Wu, M. L., (2009). Vibration-induced fatigue life estimation of ball grid array packaging, Journal Micromechanics and Microengineering 19. [ Links ]

Yang, Q. J., Wang, Z. P., Lim, G. H., Pang, J. H. L., Yap, F. F., Lin, R. M., (2002). Reliability of PBGA assemblies under out-of-plane, IEEE Transactions on components and packaging technologies 25, No. 2: 293-300. [ Links ]

Yang, Q. J., Pang, H. L. J., Wang, Z. P., Lim, G. H., Yap, F. F., Lin, R. M., (2000). Vibration reliability characterization of PBGA assemblies, Pergamon Microelectronics Reliability 40: 1097-1107. [ Links ]

Yu, D., Al-Yafawi, A., Nguyen, T. T., Park, S., Chung, S., (2011). High-cycle fatigue life prediction for Pb-free BGA under random vibration loading, Elsevier Microelectronics Reliability 51: 649-656. [ Links ]

APPENDIX A: ELECTRONIC COMPONENT KINEMATICS

The electronic component movement due to PCB displacements can be expressed as

Zc,mn(x,y)=Z0+(xe¯)θc+(yf¯)ϕc (A.1)

where the variables Z 0 , θ c and ϕ c can be evaluated according to the following equations:

Z0=4a¯b¯π2mnc¯d¯2Bmn(W)ρa¯b¯sen(nπf¯b¯)sen(nπd¯2b¯)sen(mπe¯a¯)sen(mπc¯2a¯) (A.2)

θc=12d¯c¯32b¯πmn2Bmn(W)ρa¯b¯sen(nπf¯b¯)sen(nπd¯2b¯){2a¯2π2msen(mπc¯2a¯)sen(mπe¯a¯) a¯π[c¯cos(mπe¯a¯)cos(mπc¯2a¯)2e¯sen(mπe¯a¯)sen(mπc¯2a¯)]2a¯e¯πsen(mπe¯a¯)sen(mπc¯2a¯)} (A.3)

ϕc=12c¯d¯32a¯πmn2Bmn(W)ρa¯b¯sen(nπe¯a¯)sen(nπc¯2a¯){2b¯2π2msen(mπd¯2b¯)sen(mπf¯b¯) b¯π[d¯cos(mπf¯b¯)cos(mπd¯2b¯)2f¯sen(mπf¯b¯)sen(md¯2b¯)]2b¯f¯πsen(mπf¯b¯)sen(mπd¯2b¯)} (A.4)

The PCB curvature and the lead total rotations are calculated by:

  • - PCB rotation around y axis

θPCB(x,y,W)=Zmn(x,y,W)x=2mπBmn(W)a¯ρa¯b¯cos(mπxa¯)sin(nπyb¯) (A.5)

  • - PCB rotation around x axis

ϕPCB(x,y,W)=Zmn(x,y,W)y=2nπBmn(W)bρa¯b¯sin(mπxa¯)cos(nπyb¯) (A.6)

  • - Lead lateral 1 (bending) rotation due to the P force in the extremity 2 (longitudinal stretch)

θ2=P(c22EIz'2+b2cEIz'2+b1cEIz'1+(a2+2ac)2EIz'1(a+c)dEIz'1) (A.7)

  • - Lead lateral 1 (bending) rotation around y axis

θlead,lat1(x,y,W)=[θPCB(x,y,W)θc(W)θ2(x,y,W)]fatcm1(x,y) (A.8)

  • - Lead lateral 2 (torsion) rotation around y axis

θlead,lat2(x,y,W)=[θPCB(x,y,W)θc(W)]fatcm2(x,y) (A.9)

  • - Lead lateral 1 (bending) rotation around x axis

ϕlead,lat1(x,y,W)=[ϕPCB(x,y,W)ϕc(W)θ2(x,y,W)]fatcm1(x,y) (A.10)

  • - Lead lateral 2 (torsion) rotation around x axis

ϕlead,lat2(x,y,W)=[ϕPCB(x,y,W)ϕc(W)]fatcm2(x,y) (A.11)

APPENDIX B: SUMMARY OF DIRLIK FATIGUE DAMAGE METHOD

The Dirlik fatigue damage method variables are resumed below (Mrnik 2012):

mi=12π0+WiσVM(W)dW (B.1)

α2=m2m0m4 (B.2)

xm=m1m0(m2m4)12 (B.3)

G1=2(xmα22)1+α22 (B.4)

R=α2xmG121α2G1+G12 (B.5)

G2=1α2+G1+G121+R (B.6)

G3=1G1G2 (B.7)

Q=1.25(α2G3G2R)G1 (B.8)

APPENDIX C: SUMMARY OF PCB AND EC PROPERTIES

The PCB and EC properties used on the calculations are presented on Tables C.1 and C.2. The properties were obtained preferably in the reference paper (Liguore and Followell (1995) / YU et al (2011) / Genç (2006)). Otherwise the data was taken from components datasheet or from Steinberg (2000).

Table C.1: PCB properties. 

PCB properties ā (mm) b¯ (mm) Material t PCB (mm) E PCB (Mpa) ν PCB
PLCC / LCCC 304.8(1) 152.4(1) E-Glass(1) 1.45(1) 80000 0.3
PBGA 85(2) 45(2) FR4(2) 0.6(2) 25000(2) 0.28(2)
PDIP/ Capacitor 318(3) 233(3) not available 2.5* 17000(3) 0.3(3)

(References: 1 - Liguore and Followell (1995) ; 2 - YU et al (2011) ; 3 - Genç (2006))

Table C.2: EC properties. 

EC properties mass (g) lead material E lead (Mpa) body material E body (Mpa) a (mm) b 1 (mm) b 2 (mm) c (mm)
PLCC 10* Copper(1) 128000 Plastic(1) 26000(5) 0.08(7) 1.53(7) 1.53(7) 0(7)
LCCC 32 10* Pb solder(1) 44000(6) Ceramic(1) 275700(4) 0(8) 0.305(8) 0.305(8) 0(8)
LCCC 68 10* Pb solder(1) 44000(6) Ceramic(1) 275700(4) 0(8) 0.425(8) 0.425(8) 0(8)
LCCC 84 10* Pb solder(1) 44000(6) Ceramic(1) 275700(4) 0(8) 0.425(8) 0.425(8) 0(8)
PBGA 0.21(2) SAC 305(2) 51000(2) Plastic(2) 20900(2) 0(2) 0.2(2) 0.01(2) 0(2)
PBGA 0.21(2) SAC 405(2) 51000(2) Plastic(2) 20900(2) 0(2) 0.2(2) 0.01(2) 0(2)
PDIP 1.08(3) CDA194(3) 121000(3) Plastic(3) 14500(3) 0.457(3) 1.715(3) 0.551(3) 0(3)
Cap. Tant. 6.04(3) Nickel(3) 207000(3) - 186000(3) 3.5(3) 4.5(3) 0(3) 0(3)

(References: 1 - Liguore and Followell (1995) ; 2 - YU et al (2011) ; 3 - Genç (2006) ; 4 - Steinberg (2000) ; 5 - Sumitomo Bakelite (2000); 6 - Basaran and Jiang (2002) ; 7 - Fairchild Semiconductor (2001) ; 8 - Texas Instruments (2001))

Table C.2: EC properties (continuation). 

EC prop. d (mm) l 1 (mm) l 2 (mm) t leads (mm) c¯ (mm) d¯ (mm) t c (mm) C b (MPa) b b
PLCC 0(7) 0.735(7) 0.43(7) 0.25(7) 25(7) 25(7) 4.125(7) 1.57e19(4) 6.5(4)
LCCC 32 0(8) 1.27(8) 1.27(8) cylinder 13(8) 13(8) 1.83(8) 4.032e9(4) 4(4)
LCCC 68 0(8) 1.27(8) 1.27(8) cylinder 24(8) 24(8) 2.55(8) 4.032e9(4) 4(4)
LCCC 84 0(8) 1.27(8) 1.27(8) cylinder 29(8) 29(8) 2.55(8) 4.032e9(4) 4(4)
PBGA 0(2) 0.266(2) 0.22(2) cylinder(2) 10.5(2) 10.5(2) 1(2) 4.454e12(2) 7(2)
PBGA 0(2) 0.266(2) 0.22(2) cylinder(2) 10.5(2) 10.5(2) 1(2) 2.534e10(2) 4.75(2)
PDIP 0(3) 1.65(3) 1.65(3) 0.3045(3) 19(3) 9(3) 3.43(3) 4.345e20(4) 6.5(4)
Cap. Tant. 0(3) 0.64(3) 0.64(3) cylinder(3) 9(3) 19(3) 9(3) 9.051e17(4) 6.5(4)

(References: 1 - Liguore and Followell (1995) ; 2 - YU et al (2011) ; 3 - Genç (2006) ; 4 - Steinberg (2000) ; 5 - Sumitomo Bakelite (2000); 6 - Basaran and Jiang (2002) ; 7 - Fairchild Semiconductor (2001) ; 8 - Texas Instruments (2001))

* Property not available. Value calculated in order to fit the PCB first natural frequency.

Received: June 13, 2016; Revised: April 10, 2017; Accepted: September 05, 2017

* Corresponding author

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License