Determination of discharge current equation parameters of ESD using genetic algorithms G.P. Fotis, I.F. Gonos and I.A. Stathopulos In the literature there is an absence of an accurate equation describing the current of the electrostatic discharge (ESD) phenomenon. Reported, is a method that is a genetic algorithm, which having as input data current measurements from ESD generators optimises the parameters of the discharge current’s equation.

Introduction: The IEC 61000-4-2 Standard [1] has an aberration between the typical waveform of the output current of the electrostatic discharge (ESD) generator and the discharge current that the ESD generator of the Standard produces in reality. This fact has as a result computer simulations for the circuit defined in the Standard to insert an error in the calculated voltages and currents. This problem can be solved by introducing either a new circuit of the ESD generator or an ESD current source, where the produced current waveform is a function of a number of parameters. In the work reported in this Letter the second way has been followed and the minimisation of the parameters of the current discharge equation is obtained, using genetic algorithms (GA). The GA is applied on four different types of equations. A well-known equation, which does not correspond to the discharge current, but is used in the further analysis for the application of the GA, is the equation of the lightning current given by the following equation:   ð1Þ iðtÞ ¼ i0  et=t1  et=t2

Application of GA: The GA has been used in the past and produces excellent results for computation of parameters of the earth structure [5]. A careful selection of the experimental data has to be made. In this application, the use of the GA does not require the use of the whole measured data; this would be not only a time consuming procedure, but does not give more accurate solutions than using properly selected measured data and applying a greater number of parents and iterations. For computation of the parameters for each equation the minimisation of the function Fg is necessary. Fg is given by the following equation:   N I m  I c  P i i ð5Þ Fg ¼ Iim i¼1 where Iim is the ith measured value of the discharge current. Iic is the computed value of the discharge current for the unknown parameters of (1) – (4) The waveforms’ points have been selected as follows: from 0 to 5 ns all the measured points (51 points) have been used, while for 5 to 100 ns the stepwidth was equal to 3(1 þ round(exp( j=N)) (92 points), where j is the jth point of the measured data of 1000 points. This was done in order that the GA takes all the measured points at the first 5 ns, and in order that the initial nanoseconds of the discharge current be treated as the most important part of the waveform owing to the fact that most of the radiation and the largest values of the current derivative occur during this period. Results: In Fig. 1 common graphs of the experimental data of the discharge current and the discharge current for the optimised parameter values for (1) – (4) are depicted.

where i0 is current in amperes, and t1, t2 are time constants in nanoseconds. A first approximate equation of the discharge current for commercial simulators was first introduced by [2] using a double exponential function: ð2Þ

where i1, i2 are current in amperes, and t1, t2 are time constants in nanoseconds. The reference waveform for the discharge current according to [3] is: 2

2

iðtÞ ¼ A  eðtt1 =s1 Þ þ B  t  eðtt2 =s2 Þ

discharge current, A

iðtÞ ¼ i1  et=t1  i2  et=t2

15

GA result for (1) GA result for (2) GA result for (3) GA result for (4) experimental data

10

5

ð3Þ

(tt1=s1)2

The factors A  e and B  t  e(tt2=s2) represent the narrow and broad Gaussian, respectively. In [4] the referred waveform is given by the formula:

0

n

1=n ðt=t1 Þ  e1=t2 ½t1 ðnt2 =t1 Þ t iðtÞ ¼ i1  1 þ ðt=t1 Þn   1=n ðt=t3 Þn 1=t4 t3 ðnt4 =t3 Þ t þ i2  ne 1 þ ðt=t3 Þ

ð4Þ

where i1, i2 are currents in amperes, T1, T2, T3, T4 are time constants in nanoseconds and n signifies how many times the equation can be differentiated with respect to time. Here we assume that n ¼ 3.

Experimental setup: The measurement system used was in accordance with the Standard [1] and provides high fidelity data. The current for a charging voltage of 4 kV was measured by a four-channel Tektronix oscilloscope model TDS 7254B. An ESD generator, model NSG-438 of Schaffner, was producing contact discharges and it was grounded to the earth via a ground strap. To measure the current a resistive load known as the Pellegrini target (MD 101 of Schaffner) was placed in the centre of a grounded metal plane with dimensions 1.5  1.5 m. The Pellegrini target was connected to the oscilloscope by a HF coaxial cable and an attenuator. The measurements were conducted in an anechoic chamber so that the measurement system was unaffected by the surrounding equipment and the cables were set away from the discharge point. To minimise the uncertainty of the position of the ground strap into the GA application the ground strap was at a distance 1 m from the target as the Standard defines and the loop was as large as possible.

0

20

40

60 time, ns

80

100

Fig. 1 Curve comparison between experimental data of discharge current and discharge current for optimised parameter values of (1)–(4)

The optimised parameters for each equation are as follows: for (1) i0 ¼ 7.03 A, t1 ¼ 91.80 ns, t2 ¼ 0.12 ns and Fg ¼ 28.24; for (2) i1 ¼ 21.17 A, i2 ¼ 15.93 A, t1 ¼ 53.97 ns, t2 ¼ 21.71 ns and Fg ¼ 27.21; for (3) A ¼ 4.95 A, B ¼ 0.27 A, t1 ¼ 5.18 ns, t2 ¼ 1.62 ns, s1 ¼ 9.78, s2 ¼ 54.72 and Fg ¼ 28.87; and for (4) i1 ¼ 17.46 A, i2 ¼ 7.81 A, t1 ¼ 0.75 ns, t2 ¼ 0.82 ns, t3 ¼ 3.43 ns, t4 ¼ 68.70 ns and Fg ¼ 22.91. Comparing the error (Fg) for each equation it can be concluded that the equations can be sorted as follows: (4) ! (2) ! (1) ! (3), with (4) giving the best result. Comparing the curves of Fig. 1 it is obvious that the equation which has the best fitting to the experimental data is (4). This is the most suitable of all the examined equations since it simulates the discharge current in the best way. The second more suitable equation, (2), cannot simulate the first peak of the discharge current, however it can calculate accurately the parameters of the double exponential function. Equation (3) has a shape similar to the experimental results, but inserts a higher error. Conclusions: A methodology based on a GA is proposed to calculate the parameters of the discharge current, produced by an ESD generator. The calculated discharge current is very close to the current that is measured. Evaluation of the equations and the sampling rates

ELECTRONICS LETTERS 6th July 2006 Vol. 42 No. 14

has shown that they can be sorted as follows: (4) ! (2) ! (1) ! (3). Equation (4) has the best behaviour. Therefore, a current source, which produces the ESD current of (4), is better for use than the ESD generator circuit described by IEC 61000-4-2. This equation could be included in the next revision of the Standard.

References

Acknowledgment: G.P. Fotis is supported by a PhD study scholarship from the State Scholarships Foundation of Greece.

3

1 2

4

# The Institution of Engineering and Technology 2006 13 March 2006 Electronics Letters online no: 20060767 doi: 10.1049/el:20060767

5

G.P. Fotis, I.F. Gonos and I.A. Stathopulos (School of Electrical and Computer Engineering, High Voltage Laboratory 9, Iroon Politechniou Str., 15780 Zografou, Athens, Greece)

International Standard IEC 61000-4-2: ‘Electromagnetic Compatibility (EMC), Part 4: Testing and measurement techniques, Section 2: Electrostatic discharge immunity test – Basic EMC Publication’, 1995 Cerri, G., Leo, R., and Primiani, V.M.: ‘ESD indirect coupling modelling’, IEEE Trans. Electromagn. Compat., 1996, 38, (3), pp. 274–281 Berghe, S.V., and Zutter, D.: ‘Study of ESD signal entry through coaxial cable shields’, J. Electrost., 1998, 44, pp. 135–148 Wang, K., Pommerenke, D., Chundru, R., Doren, T.V., Drewniak, J.L., and Shashindranath, A.: ‘Numerical modeling of electrostatic discharge generators’, IEEE Trans. Electromagn. Compat., 2003, 45, (2), pp. 258–270 Gonos, I.F., and Stathopulos, I.A.: ‘Estimation of multi-layer soil parameters using genetic algorithms’, IEEE Trans. Power Deliv., 2005, 20, (1), pp. 100–106

E-mail: [email protected]

ELECTRONICS LETTERS 6th July 2006 Vol. 42 No. 14

## Determination of discharge current equation parameters ... - IEEE Xplore

Jul 6, 2006 - input data current measurements from ESD generators optimises the ... computer simulations for the circuit defined in the Standard to insert.

#### Recommend Documents

Determination of discharge current equation parameters ... - IEEE Xplore
Jul 6, 2006 - Page 1 ... discharge (ESD) generator and the discharge current that the ESD generator of the Standard produces in reality. This fact has as a ...

Postbreakdown Current in MOS Structures: Extraction of ... - IEEE Xplore
series resistance. Nevertheless, as important as having a good model for the leakage current is the availability of a suitable parameters extraction strategy, and ...

IEEE Photonics Technology - IEEE Xplore
AbstractâDue to the high beam divergence of standard laser diodes (LDs), these are not suitable for wavelength-selective feed- back without extra optical ...

wright layout - IEEE Xplore
tive specifications for voice over asynchronous transfer mode (VoATM) [2], voice over IP. (VoIP), and voice over frame relay (VoFR) [3]. Much has been written ...

Device Ensembles - IEEE Xplore
Dec 2, 2004 - time, the computer and consumer electronics indus- tries are defining ... tered on data synchronization between desktops and personal digital ...

wright layout - IEEE Xplore
ACCEPTED FROM OPEN CALL. INTRODUCTION. Two trends motivate this article: first, the growth of telecommunications industry interest in the implementation ...

Modeling Current-Mode-Controlled Power Stages for ... - IEEE Xplore
Abstract- This paper continues the efforts on behaviorally modeling multiple-module inter-connected power supply systems and presents detailed derivations of ...

Evolutionary Computation, IEEE Transactions on - IEEE Xplore
search strategy to a great number of habitats and prey distributions. We propose to synthesize a similar search strategy for the massively multimodal problems of ...

I iJl! - IEEE Xplore
Email: [email protected]. Abstract: A ... consumptions are 8.3mA and 1.lmA for WCDMA mode .... 8.3mA from a 1.5V supply under WCDMA mode and.

determination of dynamical parameters of the human ...
The 3rd International Conference on ... and Virtual Engineeringâ³ ... used an install form of a plate (platform) Kistler, an electrical signal amplifier, two DAQ ...

Gigabit DSL - IEEE Xplore
(DSL) technology based on MIMO transmission methods finds that symmetric data rates of more than 1 Gbps are achievable over four twisted pairs (category 3) ...

Significance of Parameters of the Conic Equation ...
THE CONIC EQUATION HOUGH TRANSFORM FOR IMAGE ANALYSIS. 1. Significance of ... for each set of points, it has a relatively low time complexity. The advantage of the ... building identification from satellite images also uses the linear Hough ... left

IEEE CIS Social Media - IEEE Xplore
Feb 2, 2012 - interact (e.g., talk with microphones/ headsets, listen to presentations, ask questions, etc.) with other avatars virtu- ally located in the same ...

Grammatical evolution - Evolutionary Computation, IEEE ... - IEEE Xplore
definition are used in a genotype-to-phenotype mapping process to a program. ... evolutionary process on the actual programs, but rather on vari- able-length ...

SITAR - IEEE Xplore
SITAR: A Scalable Intrusion-Tolerant Architecture for Distributed Services. â. Feiyi Wang, Frank Jou. Advanced Network Research Group. MCNC. Research Triangle Park, NC. Email: {fwang2,jou}@mcnc.org. Fengmin Gong. Intrusion Detection Technology Divi

striegel layout - IEEE Xplore
tant events can occur: group dynamics, network dynamics ... network topology due to link/node failures/addi- ... article we examine various issues and solutions.

Digital Fabrication - IEEE Xplore
we use on a daily basis are created by professional design- ers, mass-produced at factories, and then transported, through a complex distribution network, to ...

Iv~~~~~~~~W - IEEE Xplore
P. Arena, L. Fortuna, G. Vagliasindi. DIEES - Dipartimento di Ingegneria Elettrica, Elettronica e dei Sistemi. Facolta di Ingegneria - Universita degli Studi di Catania. Viale A. Doria, 6. 95125 Catania, Italy [email protected]. ABSTRACT. The no

Device Ensembles - IEEE Xplore
Dec 2, 2004 - Device. Ensembles. Notebook computers, cell phones, PDAs, digital cameras, music players, handheld games, set-top boxes, camcorders, and.

Fountain codes - IEEE Xplore
7 Richardson, T., Shokrollahi, M.A., and Urbanke, R.: 'Design of capacity-approaching irregular low-density parity check codes', IEEE. Trans. Inf. Theory, 2001 ...

Quantitative determination of tip parameters in ...
the solution of a rigid dielectric problem, b the stress field is calculated using .... contrary to the behavior anticipated in contact mode imag- ing. To complement ...

Multipath Matching Pursuit - IEEE Xplore
AbstractâIn this paper, we propose an algorithm referred to as multipath matching pursuit (MMP) that investigates multiple promising candidates to recover ...

Accurate determination of saturation parameters for ...
determined from the cw data and the excited-state ab- sorption cross section is ..... when the beam is focused near the center of the crystal. This suggests that ...