# Prognostic and Remaining Life Prediction of Electronic Device under Vibration Condition Based on CPSD of MPI

Authors: Ying Chen; Ning Tang; Zenghui Yuan
Journal: Shock and Vibration
 (2016)

---

## Abstract

Prognostic of electronic device under vibration condition can help to get information to assist in condition-based maintenance and reduce life-cycle cost. A prognostic and remaining life prediction method for electronic devices under random vibration condition is proposed. Vibration response is measured and monitored with acceleration sensor and OMA parameters, including vibration resonance frequency, especially first-order resonance frequency, and damping ratio is calculated with cross-power spectrum density (CPSD) method and modal parameter identification (MPI) algorithm. Steinberg vibration fatigue model which considers transmissibility factor is used to predict the remaining life of electronic component. Case study with a test board is carried out and remaining life is predicted. Results show that with this method the vibration response characteristic can be monitored and predicted.


---

## Body

## 1. Introduction

Prognostics and Health Monitoring (PHM) of electronic devices integrates sensor data with models that enable assessment of the degradation of a product from an expected normal operating condition and also assesses the future reliability of the product based on current and historic conditions [1]. The technique is intended for condition monitoring in high reliability applications where the knowledge of impending failure is critical and the risks in terms of loss-of-functionality are too high to bear. In avionics and aerospace electronics systems, electronic components, and dynamic loading are applied to electronic devices during normal operation, which includes high temperature and temperature cycling, vibration, and humidity. Data show that 20% [1] of electronic equipment failure is related to vibration and shock. PHM of electronic under vibration condition can help electronic devices to meet the safety, reliability, maintainability, and supportability requirements and reduce life-cycle cost.PHM of electronic devices under vibration condition has been studied by some researchers, which is mainly based on two kinds of method, data-based method and model-based method. Lall et al. developed a statistical method for prognostication of area-array electronics under shock and vibration loads during vibration testing, which is based on state space vectors from resistance spectroscopy measurements [2]. Kalman filter method is used to estimate the state space feature vector, extrapolate the future state, and predict the remaining useful life. Results will help to schedule downtime for maintenance and select the appropriate time to reorder a replacement part. Aiming at the phenomenon that the recent studies which predict impending failures of power electronic components’ packaging are not explicit, they proposed monitor system for ceramic area-array-package with resistance spectroscopy and high speed cameras, predicted with system-transfer function and statistical pattern recognition techniques during shock test [3]. Results show that there is degradation in statistical confidence values, and the impending failure may have been detected before failure of the electronic assembly.Compared with data-based prognostic method, model-based method provides the relation between failure and the characteristic of product, which helps to have the proper knowledge of failures and take measures to prevent them. Gu et al. proposed a model-based health monitoring and prognostics method for assessing the reliability of a printed circuit board (PCB) under random vibration by using strain gauges and accelerometer to measure the vibration response [1]. The Physics-of-Failure fatigue model and Miner accumulation rule are used to assess the remaining life of a test board. Modal analysis with FEA software is conducted to transfer global sensor information to local sensor information. Miner rule is used to accumulate the damage of different loading condition. An uncertainty analysis of prognostics of electronics subject to random vibration, which includes measurement uncertainty, parameter uncertainty, failure criteria uncertainty, and future usage uncertainty, is also conducted [4]. A case study is presented whereby prognostics with uncertainty are applied to an electronic circuit board subject to random vibration. Gucik-Derigny et al. presented a model-based prognosis approach for assessing the remaining useful life of a PCB subject to low frequency thermal fatigue damage and high frequency of vibration damage [5]. The proposed prognosis method was also divided into fast dynamics part which includes the estimation of nonmeasurable state and slowly drifting PCB damping parameter and slow dynamic part which includes estimation of low cycle fatigue BGA solder joint interconnection damage state and parameter. Simulation is done to test the efficiency of their approach.For model-based prognosis approach, modal parameters, such as frequency and damping, are associated with the prediction of remaining useful life, which will also change during the process of work. Recent research evaluates these parameters by simulation with FEA or CalcePWA software, which may add offline procedures and may need more experience of skilled person. Online measurement of these parameters and using them properly in the prediction model is critical. Operational Modal Analysis (OMA) technique can determine the inherent properties of a structure by measuring only the response of it without using an artificial excitation. Modal parameter identification is the technique to calculate the resonance frequencies, damping, and mode shapes of a structure. In OMA all modal parameters are to be determined without knowing the excitation forces. Therefore it is normally assumed that the excitation forces are Gaussian white noise [6]; recent research has extended the excitation to harmonic excitation [7], nonstationary excitation [8], and periodic excitations [9]. OMA has widely been used in performance assessment of structure whose properties would change with time [10–12], and recently the researchers have raised the attention on health diagnosis of structure and estimating the remaining service life in civil engineering [13] and oil industry [14].The research on physics model of vibration fatigue failure can be traced back to the 1970s. After many years of practical experience, Steinberg proposed Steinberg model applied to lifetime estimation of electronic devices working under sinusoidal or random vibration conditions [15, 16]. Although Manson model and other models appeared in this field later, Steinberg model is still wildly used in engineering for its obvious physical meaning [17]. Chesné and Deraemaeker studied the application of Steinberg model in tantalum capacitor [18]. Through experiment, they provided the S-N curve in different sinusoidal sweeping-frequency vibration conditions and compared the experimental result with Finite Element Analysis (FEA) simulation, with the aim of determining parameter values of Steinberg model. Steenackers et al. proposed some principles of electronic systems to withstand high vibration and shock conditions [19]. Lim analyzed PCB’s vibration with Steinberg model in CalcePWA [20]. Yang et al. studied the impact of transmissibility factor of Steinberg model on fatigue lifetime of electronic devices [21], which represent the coupling state of PCB and its case. Neglecting this factor will result in the inaccuracy of prediction.From the above discussion, previous model based prognostic methods of electronic device under vibration condition are limited, for example, the methods lack online measurement and calculation of model parameters. In this paper, an online health prognostic method for electronic devices under random vibration condition is proposed. Vibration response is measured and monitored with acceleration sensor under white noise random excitation, and OMA parameters, including vibration resonance frequency, especially first-order resonance frequency, and damping ratio, are calculated with cross-power spectrum density (CPSD) method and modal parameter identification (MPI) algorithm. Steinberg vibration fatigue model which considers transmissibility factor is used to predict the remaining life of electronic component.The remainder of the paper is organized as follows. Section2 presents the PHM flowchart with detail steps. Section 3 introduces the cross-power spectrum density (CPSD) method and modal parameter identification (MPI) algorithm. Section 4 discusses the Steinberg model with transmissibility factor. Section 5 illuminates the real case with a test board. Section 6 gives conclusions as well as directions for the future work.


## 2. PHM Methodology for Electronics under Vibration Loading

The proposed health prognostic method for electronics device under random vibration is shown in Figure1, in which either the vibration response is monitored and analyzed or the remaining life is predicted online.Figure 1
Proposed prognostic method of electronic devices under vibration condition.The first step is to select the sensor location on the PCB and prepare the monitoring devices. Since some complicated PCBs have hundreds of components, the space is limited; it is essential to know where to place the sensors. Sensors will be placed at the site where components are most likely to fail. The vibration result of the PCB under given random vibration condition could be obtained. For example, CalcePWA software [1] can identify the natural frequency of the PCB and the locations of critical components at certain vibration loading levels. Piezoelectric sensor is used to measure the time-domain acceleration response of PCB. After that the cross spectrum of signal from different sensor is calculated, results will be used to generate the CPSD of reference point and measurement point. With frequency domain, modal extraction method calculates the first-order resonant frequency and damping of the PCB, which are the main parameters needed to calculate the remaining life time with PoF models. During the process, monitor the vibration response signals and the modal parameters and make decisions when necessary.


## 3. CPSD Based Modal Parameter Identification

### 3.1. Cross-Power Spectrum Density Function

Piezoelectric acceleration sensor can measure the system response in time domain and Discrete Fourier Transform (DFT) should be carried out to obtain the power spectral density matrices that will contain all the frequency information. The relationship between the inputx(t) and the output y(t) in frequency domain can be written in the following form:(1)Gyyjω=Hjω∗GxxjωHjωT,where Gxx(jω) is the input power spectral density matrix that is constant in the case of a stationary zero mean white noise input. This constant will be called C in the rest of the mathematical derivation. Gyy(jω) is the output PSD matrix and H(jω) is the frequency response function (FRF) matrix. As seen in (1), the output Gyy will be highly sensitive to the input constant C. The rest of the equation derivations and single degree of freedom identification will provide relevant results, only by assuming that the input is effectively represented by a constant value (mean Gaussian). It is therefore important to realize how this input assumption will be crucial to the technique.The FRF matrix can be written in a typical partial fraction form (used in classical modal analysis), in terms of poles and residues where [Gxx(ω)] is the input power spectral density (PSD) matrix, [Gyy(ω)] is the output PSD matrix, and [H(ω)] is the frequency response function (FRF) matrix, and ∗ and superscript T denote complex conjugate and transpose, respectively. The FRF matrix can be written in a typical partial fraction form (used in classical Modal analysis), in terms of poles, λ, and residues, R, which is(2)Hω=∑k=1NRkjω-λk+Rk∗jω-λk∗,λk=-ξk+jωnk,where N is the total number of modes of interest, λk the pole of the kth mode, ξk the modal damping (decay constant), and ωnk the damped natural frequency of the kth mode. (3)ωdk=ω0k1-ςk2=ω0k1-σk2ω0k,ς k being the critical damping and ω0k the undamped natural frequency, both for mode k.The transfer function matrix [H] is symmetric and an element Hpq(jω) of this matrix is then written in terms of the component rkpq(jω) of the residue matrix as follows:(4)Hpqjω=∑k=1Nrkp,qjω-λk+rkp,q∗jω-λk∗.Using expression (1) for the matrix [Gyy(ω)] and the Heaviside partial fraction theorem for polynomial expansions, the following expression for the output PSD matrix is [Gyy(ω)] assuming the input is random in both time and space and has a zero mean white noise distribution; that is, [Gxx(ω)]=[C]:(5)Gyyjω=∑k=1NAkjω-λk+Ak∗jω-λk∗+Bk-jω-λk+Bk∗-jω-λk∗,where [Ak] is the kth residue matrix of the matrix [Gyy]. The matrix Gxx is assumed to be a constant value C, since the excitation signals are assumed to be uncorrelated zero mean white noise in all the measured DOFs.In (5), the form of Gyy(ω) is similar to frequency response function, the poles also including frequency and damping information. By replacing the traditional frequency response function with cross-power spectrum density function, it is possible to identify working modal parameters.Cross-power spectrum is the product of Fourier transforms of test response signal and the reference response signal. Its amplitude is the product of the amplitude of test response signal and the reference. Its phase is the difference of the phase of test response signal and the reference. The relations between FRF and cross-power spectrum and working modal parameter identification algorithm based on cross-power spectrum will be demonstrated at the next section.


### 3.2. Modal Parameter Identification Algorithm

For a system withm measurement point and p reference point, the cross-power spectrum density matrix is (6)Gjω=G11jωG12jω⋯G1PjωG21jωG22jω⋯G2Pjω⋮⋮⋮⋮GM1jωGM2jω⋯GMPjωM×Pand then (5) can be written as(7)Gyyjω=∑k=1NψkWkTjω-λk+ψk∗WkHjω-λk∗+WkψkT-jω-λk+Wk∗ψkH-jω-λk∗in which ψk, ψk∗, Wk, and Wk∗ are complex modal matrix.Lets=jω, Os=SI-Z-1W, and (6) can be written as(8)Gs=ϕSI-Z-1W=ϕOs.The time derivative of Gs is(9)Gs=SGs-R0=SGs-ϕW=SϕOs-ϕSI-ZOs=ϕZOswhich will get(10)GsGs=ϕϕZOsorGsSGs-0ϕW=ϕϕZOs.For Eigenvalue matrix and feature vector matrix, there must be a matrix, which will fulfill the Eigenvalue function of(11)Aϕ+ϕΖ=0.Another form of (9) is(12)A⋮IϕϕZ=0.Multiply Os to (10) and get(13)A⋮IϕϕZOs=A⋮IGsSGs-ϕW=0,(14)AGs+SGs-ϕW=0Consider S=jω; then (12) can be given by(15)AGjωk+jωkGjωk=ϕW,k=1,2,…,K.Let(16)D=Gjω1,Gjω2,…,Gjω1K,Ω=-diag⁡jω1I,jω2I,…,jωkIand then(17)A⋮-ϕWDI1,I2,…,Ik=DΩ.D and Ω are all functions of measured CPSD and ωi, which are known; the matrix A and ϕW can be solved by (14):(18)Modal  frequency:ωnk=Reλk2+Imλk2Dampingξk=-Reλkωnk,n=1,…,N.


## 3.1. Cross-Power Spectrum Density Function

Piezoelectric acceleration sensor can measure the system response in time domain and Discrete Fourier Transform (DFT) should be carried out to obtain the power spectral density matrices that will contain all the frequency information. The relationship between the inputx(t) and the output y(t) in frequency domain can be written in the following form:(1)Gyyjω=Hjω∗GxxjωHjωT,where Gxx(jω) is the input power spectral density matrix that is constant in the case of a stationary zero mean white noise input. This constant will be called C in the rest of the mathematical derivation. Gyy(jω) is the output PSD matrix and H(jω) is the frequency response function (FRF) matrix. As seen in (1), the output Gyy will be highly sensitive to the input constant C. The rest of the equation derivations and single degree of freedom identification will provide relevant results, only by assuming that the input is effectively represented by a constant value (mean Gaussian). It is therefore important to realize how this input assumption will be crucial to the technique.The FRF matrix can be written in a typical partial fraction form (used in classical modal analysis), in terms of poles and residues where [Gxx(ω)] is the input power spectral density (PSD) matrix, [Gyy(ω)] is the output PSD matrix, and [H(ω)] is the frequency response function (FRF) matrix, and ∗ and superscript T denote complex conjugate and transpose, respectively. The FRF matrix can be written in a typical partial fraction form (used in classical Modal analysis), in terms of poles, λ, and residues, R, which is(2)Hω=∑k=1NRkjω-λk+Rk∗jω-λk∗,λk=-ξk+jωnk,where N is the total number of modes of interest, λk the pole of the kth mode, ξk the modal damping (decay constant), and ωnk the damped natural frequency of the kth mode. (3)ωdk=ω0k1-ςk2=ω0k1-σk2ω0k,ς k being the critical damping and ω0k the undamped natural frequency, both for mode k.The transfer function matrix [H] is symmetric and an element Hpq(jω) of this matrix is then written in terms of the component rkpq(jω) of the residue matrix as follows:(4)Hpqjω=∑k=1Nrkp,qjω-λk+rkp,q∗jω-λk∗.Using expression (1) for the matrix [Gyy(ω)] and the Heaviside partial fraction theorem for polynomial expansions, the following expression for the output PSD matrix is [Gyy(ω)] assuming the input is random in both time and space and has a zero mean white noise distribution; that is, [Gxx(ω)]=[C]:(5)Gyyjω=∑k=1NAkjω-λk+Ak∗jω-λk∗+Bk-jω-λk+Bk∗-jω-λk∗,where [Ak] is the kth residue matrix of the matrix [Gyy]. The matrix Gxx is assumed to be a constant value C, since the excitation signals are assumed to be uncorrelated zero mean white noise in all the measured DOFs.In (5), the form of Gyy(ω) is similar to frequency response function, the poles also including frequency and damping information. By replacing the traditional frequency response function with cross-power spectrum density function, it is possible to identify working modal parameters.Cross-power spectrum is the product of Fourier transforms of test response signal and the reference response signal. Its amplitude is the product of the amplitude of test response signal and the reference. Its phase is the difference of the phase of test response signal and the reference. The relations between FRF and cross-power spectrum and working modal parameter identification algorithm based on cross-power spectrum will be demonstrated at the next section.


## 3.2. Modal Parameter Identification Algorithm

For a system withm measurement point and p reference point, the cross-power spectrum density matrix is (6)Gjω=G11jωG12jω⋯G1PjωG21jωG22jω⋯G2Pjω⋮⋮⋮⋮GM1jωGM2jω⋯GMPjωM×Pand then (5) can be written as(7)Gyyjω=∑k=1NψkWkTjω-λk+ψk∗WkHjω-λk∗+WkψkT-jω-λk+Wk∗ψkH-jω-λk∗in which ψk, ψk∗, Wk, and Wk∗ are complex modal matrix.Lets=jω, Os=SI-Z-1W, and (6) can be written as(8)Gs=ϕSI-Z-1W=ϕOs.The time derivative of Gs is(9)Gs=SGs-R0=SGs-ϕW=SϕOs-ϕSI-ZOs=ϕZOswhich will get(10)GsGs=ϕϕZOsorGsSGs-0ϕW=ϕϕZOs.For Eigenvalue matrix and feature vector matrix, there must be a matrix, which will fulfill the Eigenvalue function of(11)Aϕ+ϕΖ=0.Another form of (9) is(12)A⋮IϕϕZ=0.Multiply Os to (10) and get(13)A⋮IϕϕZOs=A⋮IGsSGs-ϕW=0,(14)AGs+SGs-ϕW=0Consider S=jω; then (12) can be given by(15)AGjωk+jωkGjωk=ϕW,k=1,2,…,K.Let(16)D=Gjω1,Gjω2,…,Gjω1K,Ω=-diag⁡jω1I,jω2I,…,jωkIand then(17)A⋮-ϕWDI1,I2,…,Ik=DΩ.D and Ω are all functions of measured CPSD and ωi, which are known; the matrix A and ϕW can be solved by (14):(18)Modal  frequency:ωnk=Reλk2+Imλk2Dampingξk=-Reλkωnk,n=1,…,N.


## 4. Remaining Useful Life Prediction

Steinberg model is used for remaining useful life prediction, which is derived from component fatigue characteristic fitting from fatigue experiments data [22], where (19)N1S1b=N2S2bin whichN1 and N2 are the stress cycle before fatigue. S1 and S2 are the stress when failure happens, and b is the fatigue index number related to the linear gradient of fatigue curve. For linear system, stress S is proportionate to displacement Z, and (19) can be written as(20)N1Z1b=N2Z2b,N2=N1Z1Z2b.Under vibration condition, the displacement of simply supported PCB with four sides can expressed as follows, where the fatigue lifetime of these component can reach1×107 stress cycles [22]: (21)Z1=0.00022BChRxyL,where B is the length of the circuit board edge parallel to the component, L is the length of the electronic component, h is the circuit board thickness, C is the constant for different types of electronic components (0.75<C<2.25, and refer to [23]), and Rxy is the relative position factor for the component mounted on the board. (22)Rxy=sin⁡πXasin⁡πYb,where X and Y are horizontal and longitudinal position of the component and a and b are length and width of PCB.Based on the vibration fatigue curve, incorporating the system dynamic response property analysis, Steinberg proposed a vibration fatigue model for electronic device. In his model, PCB can be approximated as a single degree of freedom system, when it vibrates under the fundamental resonance.When the input power spectral density (PSD) of random vibration is flat spectrum in resonance region, the mean square root acceleration response of a system is given by(23)Gout=π2PfnQRMS.The actual dynamic single amplitude displacement of PCB’s center is given by [15](24)Z2=9.8Goutfn2=12.3PQfn1.25,where P is the input PSD at resonant frequency, fn is resonant frequency, and Q is the transmissibility at resonant frequency, which is given by [15] (25)Q=1+2RΩRc21-RΩ22+2RΩRc21/2,(26)RΩ=ωωn,Rc=ξ1ξc,where ωn=K/m, being the first-order angular resonance frequency, K is rigidity, m is mass, ξ1 is the damping at first-order frequency, and(27)ξc=2Km.Thefn and damping ratio is measured and calculated using (18).


## 5. Case Study

### 5.1. The Test Board

A certain type of test board (shown in Figure2) is mainly used for ground automatic control systems. There are CPU microprocessor, memory and cache, sensor chips, counters, and various interfaces on the single board, and its function is like a microcomputer. There are many types of components on the test board including integrated circuits, chip capacitors, tantalum capacitors, metal-film resistors, inductors, connectors, crystals, switches, and PCB (including the solder joint, pad, PTH, vias, and metal interconnects). The material of PCB is FR-4; material of metal wires and plated through holes (PTH) is copper, and solder is 63Sn37Pb. The main types of package of the device are BGA, SOP, QFP, DIP, TSSOP, and so forth. And packaging materials are mainly plastic or ceramic. PCB is mounted through the holes on its four edges side, and it is considered that the installation site will not fail.Figure 2
Test board and acceleration sensor layout.


### 5.2. The Monitoring Devices

The monitoring system is shown in Figure3, which mainly includes the vibration table (LDS), acceleration sensor, vibration measurement instrument, and data collection and analysis software. The working frequency of LDS vibration table is 5–4000 Hz, and vibration table can vibrate horizontally and vertically. The payload weight of the vibration table is about 200 Kg. Vibration Signal Control System is the brain of LDS vibration table, which has the function of random and sinusoidal vibration, harmony frequency lingering, and shock wave controlling. The vibration signal should be amplified by a power amplifier before sent to the vibration table. Endevco 2222c piezoelectric sensor whose weight is only 2 grams is used to measure the acceleration response of PCB. Vibration measurement instrument is DP730 24-channel device which could analyze polytype of vibration and noise in real time with software. The monitoring system is composed of data collection module, modal parameter identification module, and RUL calculation module.Figure 3
Monitoring devices of electronic devices under vibration.The test board was mounted on vibration table, and its two edges were clamped, as shown in Figure3. Random vibration loading condition with a frequency from 10 to 2000 Hz is applied to the vibration table, which is shown in Figure 4. Vibration spectrum is given by LDS vibration control system and the signal is exported into a power amplifier to drive the vibration table. The vibration table would generate vibration in the vertical direction.Figure 4
Exaction vibration spectrum, a white noise spectrum.There are mainly four ways to connect sensors and the test object, by bolt, adhesives, wax, or permanent magnet. The Endevco 7703A-50 acceleration sensor in this test is with no magnetic base and requires high upper limited frequency. So we fix them via hard adhesives. Theoretically, more acceleration sensors can improve the measurement and prediction accuracy. But restricted to the space of PCB, the number of sensors is limited. This may be better if wireless and built-in sensors are used. In our case study, the layout of the test board is very compact and most of the components are on the front of the PCB. Modal analysis and random vibration analysis are carried out with CalcePWA. Results show the displacement response of this test board is well-distributed. Based on these pieces of information, we arranged eight acceleration sensors on the well-distributed eight locations where there is no component to avoid interaction between them.The broadband random vibration spectrum is set close to white noise and loaded on the vibration table along the normal direction of the PCB. Via the fixture, excitation signal will be transmitted from the vibration table to the electronic equipment.


### 5.3. Prognostic Data and Remaining Life Prediction

The acceleration data were collected by the DP730 vibration instrument. Meanwhile, the acceleration data was sent to the data collection module which sample and record the data in real time. The modal parameter identification module will identify working modal parameters including the natural frequency and damping of the test board in accordance with the cross-power spectrum. Figure4 is the exaction vibration spectrum, which is a white noise spectrum. Figure 5 shows the time-domain data of the measurement points 1 to 8, which express the dynamic response of the PCB. Choose point 1 as the reference point. Figure 6 is the cross-power spectrum density curve of points 2–8 with reference point 1.Figure 5
Time-domain dynamic response of measurement points 1 to 8.Figure 6
Cross-power spectrum density curve.The times series data could be processed using the Fourier transform method. Multiplying the reference point’s Fourier transform data by the conjugation of other measuring points, the cross-power spectrum could be obtained, which illustrate the correlation of the two-channel vibrate data. The cross-power spectrum density curves in Figure6 indicate the correlation between the monitoring point and reference point. PCB is mounted through the holes on its four edges side. Monitoring points 1–6 are located around the PCB, but points 7 and 8 are located in the center of the PCB. So the curves’ shape between monitoring points 2–6 and point 1 is consistent, and similarly, the curves’ shape between monitoring points 2–6 and point 1 is consistent. Due to the transmission delay of the vibration signal on the board, the first curves and the last two curves are different. The transmission delay of the vibration signal has an effect on the cross-power spectrum density curves. We can find that the vibration magnitude of output signal is more dependent on the input signal when the frequency component of input signal approach 1550 Hz or 1850 Hz. The resonance frequency is calculated based on the modal parameter identification algorithm according to the first six curves. The peaks in Figure 6 indicate that the vibration data of other measuring point has significant correlation with the reference point.Through the identification of modal parameters, the first three work mode shapes of the PCB board would be obtained. Modal analysis and random vibration analysis are carried out with CalcePWA. The modal analysis results of cross-power spectrum estimation and the modal analysis results of the CalcePWA are consistent.The remaining useful life of some important component predicted with (20)–(25) is presented in Table 1, in addition of the estimated first-order modal frequency, damping ratio at the beginning of the monitor (t=0), and when t=28 days, t=63 days, and t=121 days the data are also given. It can be seen from the results that, without consideration of measurement and computation deviation, the first-order modal frequency of the test board is not significantly changed; the damping will increase with the time.Table 1
Remaining life prediction results.
Component Time(days) First-order modal Frequency Damping Remaining life prediction(days) BGA64 t = 0 231.97 3.69 332.15 t = 28 227.41 3.73 298.73 t = 63 222.76 3.76 220.67 t = 121 220.32 3.78 183.32  QFP100 t = 0 231.97 3.69 384.29 t = 28 227.41 3.73 349.50 t = 63 222.76 3.76 281.25 t = 121 220.32 3.78 249.66  QFP48 t = 0 231.97 3.69 435.77 t = 28 227.41 3.73 405.41 t = 63 222.76 3.76 330.92 t = 121 220.32 3.78 387.45  SOP32 t = 0 231.97 3.69 592.86 t = 28 227.41 3.73 563.43 t = 63 222.76 3.76 478.32 t = 121 220.32 3.78 449.66Table1 shows that the prediction for four components. All but BGA64 of the components lasted more than one year. BGA64 and QFP100, which have large dimension and are located at the center of the test board, where the vibration response is strong than other components, have less lifetime.


## 5.1. The Test Board

A certain type of test board (shown in Figure2) is mainly used for ground automatic control systems. There are CPU microprocessor, memory and cache, sensor chips, counters, and various interfaces on the single board, and its function is like a microcomputer. There are many types of components on the test board including integrated circuits, chip capacitors, tantalum capacitors, metal-film resistors, inductors, connectors, crystals, switches, and PCB (including the solder joint, pad, PTH, vias, and metal interconnects). The material of PCB is FR-4; material of metal wires and plated through holes (PTH) is copper, and solder is 63Sn37Pb. The main types of package of the device are BGA, SOP, QFP, DIP, TSSOP, and so forth. And packaging materials are mainly plastic or ceramic. PCB is mounted through the holes on its four edges side, and it is considered that the installation site will not fail.Figure 2
Test board and acceleration sensor layout.


## 5.2. The Monitoring Devices

The monitoring system is shown in Figure3, which mainly includes the vibration table (LDS), acceleration sensor, vibration measurement instrument, and data collection and analysis software. The working frequency of LDS vibration table is 5–4000 Hz, and vibration table can vibrate horizontally and vertically. The payload weight of the vibration table is about 200 Kg. Vibration Signal Control System is the brain of LDS vibration table, which has the function of random and sinusoidal vibration, harmony frequency lingering, and shock wave controlling. The vibration signal should be amplified by a power amplifier before sent to the vibration table. Endevco 2222c piezoelectric sensor whose weight is only 2 grams is used to measure the acceleration response of PCB. Vibration measurement instrument is DP730 24-channel device which could analyze polytype of vibration and noise in real time with software. The monitoring system is composed of data collection module, modal parameter identification module, and RUL calculation module.Figure 3
Monitoring devices of electronic devices under vibration.The test board was mounted on vibration table, and its two edges were clamped, as shown in Figure3. Random vibration loading condition with a frequency from 10 to 2000 Hz is applied to the vibration table, which is shown in Figure 4. Vibration spectrum is given by LDS vibration control system and the signal is exported into a power amplifier to drive the vibration table. The vibration table would generate vibration in the vertical direction.Figure 4
Exaction vibration spectrum, a white noise spectrum.There are mainly four ways to connect sensors and the test object, by bolt, adhesives, wax, or permanent magnet. The Endevco 7703A-50 acceleration sensor in this test is with no magnetic base and requires high upper limited frequency. So we fix them via hard adhesives. Theoretically, more acceleration sensors can improve the measurement and prediction accuracy. But restricted to the space of PCB, the number of sensors is limited. This may be better if wireless and built-in sensors are used. In our case study, the layout of the test board is very compact and most of the components are on the front of the PCB. Modal analysis and random vibration analysis are carried out with CalcePWA. Results show the displacement response of this test board is well-distributed. Based on these pieces of information, we arranged eight acceleration sensors on the well-distributed eight locations where there is no component to avoid interaction between them.The broadband random vibration spectrum is set close to white noise and loaded on the vibration table along the normal direction of the PCB. Via the fixture, excitation signal will be transmitted from the vibration table to the electronic equipment.


## 5.3. Prognostic Data and Remaining Life Prediction

The acceleration data were collected by the DP730 vibration instrument. Meanwhile, the acceleration data was sent to the data collection module which sample and record the data in real time. The modal parameter identification module will identify working modal parameters including the natural frequency and damping of the test board in accordance with the cross-power spectrum. Figure4 is the exaction vibration spectrum, which is a white noise spectrum. Figure 5 shows the time-domain data of the measurement points 1 to 8, which express the dynamic response of the PCB. Choose point 1 as the reference point. Figure 6 is the cross-power spectrum density curve of points 2–8 with reference point 1.Figure 5
Time-domain dynamic response of measurement points 1 to 8.Figure 6
Cross-power spectrum density curve.The times series data could be processed using the Fourier transform method. Multiplying the reference point’s Fourier transform data by the conjugation of other measuring points, the cross-power spectrum could be obtained, which illustrate the correlation of the two-channel vibrate data. The cross-power spectrum density curves in Figure6 indicate the correlation between the monitoring point and reference point. PCB is mounted through the holes on its four edges side. Monitoring points 1–6 are located around the PCB, but points 7 and 8 are located in the center of the PCB. So the curves’ shape between monitoring points 2–6 and point 1 is consistent, and similarly, the curves’ shape between monitoring points 2–6 and point 1 is consistent. Due to the transmission delay of the vibration signal on the board, the first curves and the last two curves are different. The transmission delay of the vibration signal has an effect on the cross-power spectrum density curves. We can find that the vibration magnitude of output signal is more dependent on the input signal when the frequency component of input signal approach 1550 Hz or 1850 Hz. The resonance frequency is calculated based on the modal parameter identification algorithm according to the first six curves. The peaks in Figure 6 indicate that the vibration data of other measuring point has significant correlation with the reference point.Through the identification of modal parameters, the first three work mode shapes of the PCB board would be obtained. Modal analysis and random vibration analysis are carried out with CalcePWA. The modal analysis results of cross-power spectrum estimation and the modal analysis results of the CalcePWA are consistent.The remaining useful life of some important component predicted with (20)–(25) is presented in Table 1, in addition of the estimated first-order modal frequency, damping ratio at the beginning of the monitor (t=0), and when t=28 days, t=63 days, and t=121 days the data are also given. It can be seen from the results that, without consideration of measurement and computation deviation, the first-order modal frequency of the test board is not significantly changed; the damping will increase with the time.Table 1
Remaining life prediction results.
Component Time(days) First-order modal Frequency Damping Remaining life prediction(days) BGA64 t = 0 231.97 3.69 332.15 t = 28 227.41 3.73 298.73 t = 63 222.76 3.76 220.67 t = 121 220.32 3.78 183.32  QFP100 t = 0 231.97 3.69 384.29 t = 28 227.41 3.73 349.50 t = 63 222.76 3.76 281.25 t = 121 220.32 3.78 249.66  QFP48 t = 0 231.97 3.69 435.77 t = 28 227.41 3.73 405.41 t = 63 222.76 3.76 330.92 t = 121 220.32 3.78 387.45  SOP32 t = 0 231.97 3.69 592.86 t = 28 227.41 3.73 563.43 t = 63 222.76 3.76 478.32 t = 121 220.32 3.78 449.66Table1 shows that the prediction for four components. All but BGA64 of the components lasted more than one year. BGA64 and QFP100, which have large dimension and are located at the center of the test board, where the vibration response is strong than other components, have less lifetime.


## 6. Discussion and Conclusion

This paper presents an online health prognostic and monitoring method for assessing the reliability of electronic devices under random vibration condition, using acceleration sensor’s information measured on the place where there are some critical components and the acceleration sensor could be placed perfectly. Based on cross-power spectrum density (CPSD) method and modal parameter identification (MPI) algorithm, the relationship between modal parameter and the remaining life is established. In the case study, we can assess the reliability of electronic devices without modal analysis and other simulation using software. As a result, this paper proved that an accelerometer is suitable for PHM and makes the field application of PHM much easier, since its installation and associated signal conditioning circuitry are simpler, cheaper, and more reliable.Future work should focus on increasing PHM accuracy, by using wireless sensor systems to enhance prediction capability and accuracy by developing or application of algorithms for other vibration excitation, for example, nonwhite noise spectrum.



---

*Source: 2650569-2016-11-23.xml*