420 IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 153 HARMONIC DISTORTION ASSESSMENT BY AREA BASED APPROACH AT SINGLE PHASING OF AN INDUCTION MOTOR SURAJIT CHATTOPADHYAY AND SAMARJIT SENGUPTA 1 Dept of Electrical Engineering, MCKV Institute of Engineering, (West Bengal University of Technology), 243 GT Road (N), Howrah, West Bengal, India 2 Dept of Applied Physics, Prof. and Head, University of Calcutta, APC Road 92, Kolkata, 700009, West Bengal, India; surajitchattopadhyay@gmail.com; samarsgp@rediffmail.com ABSTRACT: This paper presents harmonic distortion assessment by area based approach at single phasing of an induction motor. This has been achieved by assessing voltage and current signals using area based approach. Some specific reference signals have been defined, after which, the real power system data are plotted with this reference signal and areas. Formed by the real power system data with the reference signal, the reference signal and areas have been calculated, where the contributions of fundamental waveform and the harmonic components in real and reactive powers are being assessed separately. Single phasing is done on induction machine as well as the total harmonic distortion factors have been calculated afterwards. Moreover, the significant change is observed in harmonic distortion due to single phasing. ABSTRAK: Kertas kerja ini membentangkan penilaian herotan harmonik menggunakan kaedah keluasan kawasan pada pemfasaan tunggal motor aruhan. Menggunakan kaedah keluasan kawasan, penilaian terhasil dengan memantau isyarat arus dan voltan. Sesetengah isyarat rujukan tertentu dikenal pasti, di mana, data sistem kuasa sebenar diplotkan berdasarkan isyarat rujukan ini. Kawasan kemudiannya dibentuk dengan adanya data sistem kuasa sebenar dengan menggunakan pengiraan isyarat rujukan. Pengiraan ini memberikan bentuk gelombang asas dan komponen harmonik sebenar di mana kuasa reaktif ditentukan secara berasingan. Pemfasaan tunggal ditentukan menggunakan mesin aruhan dan faktor jumlah herotan harmonik diambil kira .Perubahan yang ketara dikenal pasti dalam herotan harmonik yang disebabkan oleh pemfasaan tunggal. KEYWORDS: Area Based Approach; Harmonics; Single-Phasing. 1. INTRODUCTION Maintaining the near sinusoidal waveform of power distribution bus voltages and currents at rated magnitude and frequency is important for good quality of electric power in respect of harmonics. Harmonics are generated from various types of sources like actual loads, equipment and components, subsystems of transmission and distribution systems, etc. With the increasing use of solid state control devices, power converter, nonlinear loads like switch-mode-power-supplier and phase control rectifiers in distribution and transmission networks, the amount of such disturbances is also increasing. Different type of faults occurred in electrical machine also contribute in harmonic distortion. It leads IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 154 catastrophic consequences such as long production downtimes, mal-function of devices and shortened equipment life. Use of induction motor has increased and it demands an early and accurate determination of motor fault for its reliable operation. Single phasing is one common fault often observed in industrial application. Research work is going on to study the sources and effects of faults. A detail review on measurement of harmonic contribution in impedances and different issues of their practical implementation has been done [1]. Modeling is done for large load areas considering harmonics in distributions networks [2]. Modeling is also done for assessing harmonics and inter-harmonic distortion in multi-converter system [3]. Measurement technique is introduced for determination of harmonic/ inter-harmonics of time varying frequencies [4]. Different passive models are introduced for motor control [5-11]. A critical impedance based method has been developed for identifying harmonic sources [5]. Passivity-based control system is developed for Euler-Lagrange systems [7]. A lot of experimental investigations have been done on the passivity based induction motor control [8-10]. Modeling of adjustable speed drives for power system harmonic analysis [11] has been done. For this purpose different electric quantities under sinusoidal, non sinusoidal, balanced or unbalanced conditions [12] have been defined. For motor control and fault diagnosis, analysis of air-gap flux, current and vibration signals as a function of the combination of static and dynamic air-gap eccentricity in 3-phase induction motors [13] has been done. For fault diagnosis, use of motor current signature analysis (MCSA) [14] [15] has become very popular. Condition monitoring of induction motor is also being done by MCSA [15]. A detail review of induction motor signature analysis as a medium for Fault detection has been provided in [16]. Some online diagnosis of induction motors using MCSA [14, 17, 18] have been introduced. Wavelet and PSD are being used as fault detection technique [19]. Induction motor stator faults diagnosis has been achieved by a current Concordia pattern-based fuzzy decision system [20]. Fault detection and diagnosis in an induction machine drive is done using pattern recognition approach based on Concordia stator mean current vector [21]. Monitoring and diagnosis of induction motor electrical faults using a current Park's vector pattern approach [22-25]. Current Concordia based assessment of single phasing of an induction motor by feature pattern extraction method and radar analysis [26] has been introduced. Assessment of crawling of an induction motor has also been done by stator current Concordia analysis [27]. In this paper, stator current of induction motor at single phasing are assessed by area based approach. Some specific reference signals having sinusoidal nature and different frequencies have been defined. Areas formed by the real power system data with the reference signal have been calculated. Mathematical relations have been established between those areas with amplitude and phase angles of voltage and current wave forms wherefrom contributions of fundamental waveform and harmonic components have been assessed separately. Active power, reactive power and total harmonic distortion factors have been measured. Distortion factors of obtained from normal and fault condition have been compared. 2. AREA BASED APPROACH FOR HARMONIC ASSESSMENT 2.1 Power in Terms of Area Consider v(t) be the voltage and i(t) be the current in a power system. In presence of harmonics, these waveforms can be written as follows IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 155 )sin()( ,.....3,2,1 n nn n tnVtv V ϕω −= ∑ = (1) )(sin)( ,......3,2,1 n nn n tnIti I θω −= ∑ = (2) where, n denotes order of harmonics, nV and nI be the highest order of harmonics present in voltage waveform and current waveform respectively. Fundamental frequency (ω) is assumed to be constant in single complete cycle. If these waveforms are plotted in voltage – current plane (v-i plane), covered area can be written as follows ∑∫ == − n nnn iv TOTAL InVkvdiA φsin 1 (3) where, ωθϕφ 1.01 =−= kandnnn From (1), iv TOTAL n nn n nn AKnQInV − = == ∑∑ 1 ,......3,2,1 2 1 )(sin 2 1 φ (4) where, 1 1 1 k K = If voltage and current are plotted in vi-t plane, covered area can be written as follows ∑∫ ∑∑∫ =−−== == − n nnnn n nn n n tvi TOTAL IVkdttnItnVvidtA φθωϕω cos)(sin)sin( 2 ,......3,2,1,......3,2,1 (5) Where, tconsaiskand nnn tan 2 θϕφ −= From (5), tvi TOTAL n nn n nn AKPIV − = == ∑∑ 2 ,......3,2,1 2 1 )(cos 2 1 φ (6) where, 2 2 1 k K = Fundamental and harmonic components of same order, which are present in both voltage and current waveform, contribute to the active and reactive power (4, 6). Equation (3) shows that area in v-i plane contains the information of sine of angle between the voltage and current wave form whereas, (5) shows that the area under the curve formed by vi in vi-t plane gives the information of cosine of the phase angle difference between voltage and current waveform. In the proposed scheme, a reference signal (M) of sinusoidal nature has been defined and voltage (v) and current (i) signals have been plotted with the reference signal (M) in v-M, vM-t, i-M and iM-t planes respectively, areas have been calculated wherefrom amplitude and phase angles of voltage and current signals of fundamental as well as harmonic components have been isolated. IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 156 2.2 Reference Signals For assessment of distortion factor, contribution of fundamental and harmonic components should be isolated. To extract powers contributed only by fundamental waveform a reference signal of unity amplitude is defined with the following properties ttM 11 sin)( ω= (7) Here ω1 is frequency of the reference signal. The frequency of the reference signal ω1 should be equal to the fundamental frequency ω of the real signal collected from the power system, which is initially an unknown quantity. To select the value of ω1, first take, ωω ∆±= 50 1 , where, ω∆ is a small incremental value. If this signal is plotted against any one of the signals which have been derived from the real power system data, then the starting and ending point will not meet each other. There will be a finite gap between these two points as shown in Fig. 1. Fig. 1 Curves formed by M1 and voltage having unequal fundamental frequency. If ω1 is changed gradually by changing ω∆ in such a way that the starting and ending points come closer to each other then, it can be said that ω1 is approaching towards ω. When starting and ending points come at minimum distance from each other, ideally zero, then it can be said that ω1= ω. Then the reference signal will become ttM ωsin)(1 = (8) It should be noted that the reference signal does not consist of harmonic components which may be in the real data. As a result of this, area formed by reference and real signals will not depend on any other frequency except fundamental. Real voltage and current signals are plotted with the reference signal (8) and the areas covered by the curve in one complete cycle have been calculated wherefrom amplitude and phase angle of fundamental component, have been derived for assessment of contribution of IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 157 fundamental component in distortion. Similarly, reference signal for assessment of harmonic component is defined as tntM n ωsin)( = (9) 2.3 Fundamental Component in Phase R Covered by vR and M1 waveforms in vR –M1 plane can be written as follows 1111 sin RR Mv R VkA ϕ= − Mv RRR AKVor − = 1111 sin, ϕ (10) Covered area in VR M1 -t can be written as follows: 1121 cos RR tvM R VkA ϕ= − tvM RRR AKVor − = 1211 cos, ϕ (11) From (10) and (11), ( ) ( )212 2 111 tvM R Mv RR AKAKV −− += (12)         = − − − tvM R Mv R R AK AK 12 111 1 tanϕ (13) Covered area in iR – M1 plane can be written as follows: 1111 sin RR Mi R IkA θ= − Mi RRR AKIor − = 1111 sin, θ (14) Covered area in iR M1 -t plane can be written as follows: 1121 cos RR tiM R IkA θ= − tiM RRR AKIor − = 1211 cos, θ (15) From (14) and (15), ( ) ( )212 2 111 tiM R Mi RR AKAKI −− += (16)         = − − − tiM R Mi R R AK AK 12 111 1 tanθ (17) 2.4 Harmonic Component in Phase R Covered area by curve vR and Mn waveforms in vR –Mn can be written as follows: RnRn Mv Rn VkA ϕsin1= − Mv RnRnRn AKVor − = 1sin, ϕ (18) Covered area in VR Mn -t can be written as follows: RnRn tvM Rn VkA ϕcos2= − tvM RnRnRn AKVor − = 2cos, ϕ (19) From (18) and (19), ( ) ( )22 2 1 tvM Rn Mv RnRn AKAKV −− += (20) IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 158         = − − − tvM Rn Mv Rn Rn AK AK 2 11 tanϕ (21) Covered area in iR – Mn plane can be written as follows RnRn Mi Rn IkA θsin1= − Mi RnRnRn AKIor − = 1sin, θ (22) Covered area formed in iR Mn -t plane can be written as follows RnRn tiM Rn IkA θcos2= − tiM RnRnRn AKIor − = 2cos, θ (23) From (22) and (23), ( ) ( )22 2 1 tiM Rn Mi RnRn AKAKI −− += (24)         = − − − tiM Rn Mi Rn Rn AK AK 2 11 tanθ (25) 2.5 Harmonic Distortion in Phase R Harmonics distortion is expressed in terms of total harmonic distortion factor. For phase R, using (12), (16), (20) and (24) distortion factors can be written as follows: ( ) ( ) ( ) ( )212 2 11 ,....4,3,2 2 2 2 1 1 .,.........3,2 2 )( tvM R Mv R n tvM Rn Mv Rn R n Rn RV AKAK AKAK V V THD −− = −− = + + == ∑∑ (26) ( ) ( ) ( ) ( )212 2 11 ,....4,3,2 2 2 2 1 1 ,........3,2 2 )( tiM R Mi R n tiM Rn Mi Rn R n Rn RI AKAK AKAK I I THD −− = −− = + + == ∑∑ (27) 3. EXPERIMENTATION Block diagram for experimentation is shown in Fig. 2. An induction motor (Ventwell Corporation make, 415 volts, 3 phase, 50 Hz, Squirrel cage, 1 HP, 1425 rpm) is fed by balanced supply and is run at no-load. The line currents are stepped down with current transformers (CT) and voltages are stepped down by potential transformer. They are then sampled and digitized and captured thorough data acquisition system (DAS). Each phase is connected through separate switch. Suddenly one phase (in Fig. 2 it is phase B) is switched off and single phasing is done. Voltage and current signals of healthy phases are observed for analysis. IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 159 Fig. 2: Block diagram for experimentation. 4. RESULTS The motor is first run at normal condition and harmonics distortions are measured. Then, three cases single phasing are done by switching off phase R, phase Y and phase B respectively. For each case harmonics distortions of healthy phase are measured using area based approach. The result is compared with the normal value as presented in Table 1. Significant rise of harmonics distortions are observed. Table 1:Distortions factors at normal and fault condition. Condition Current distortion in healthy phase Voltage distortion in healthy phase Normal 1.1586 1.1592 Single phasing at phase R 16.4655 8.3450 Single phasing at phase Y 16.4620 8.2550 Single phasing at phase B 16.5685 8.6575 Accuracy increases with the increase of sampling rate of the data acquisition system. The same technique may be extended in Clarke plane which will reduce the computational effort. 5. CONCLUSION In this approach fundamental frequency (ω) has been assumed constant only for one complete cycle and may change in next cycle, which is more advantageous in real system analysis than other methods where (ω) is to be constant throughout the process. Small error of the method strengthens the acceptability of the above methods. Single phasing has shown significant rise of harmonics distortion that has been detected by area based approach. This is done by capturing voltage and current signals. Some specific reference signals are defined, then, real power system data are plotted with this reference signal and areas are then calculated wherefrom contributions of fundamental waveform and harmonic components in real and reactive powers have been assessed separately. Single phasing is done on induction machine in each phase and harmonic distortion factors are calculated. Significant change of distortion is observed due to single phasing. The application area may be extended in Clarke plane for less computational effort. IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 160 REFERENCES [1] W. Xu, E. Ahmed, and X. Zhang, "Measurement of Harmonic Impedances: Practical Implementation Issues and Solutions", IEEE Transaction. on Power Delivery, 17.1 (2002):210-16. [2] R. Lamedica, A. Prodenzi, E. Tironi, D. Zaninelli, “A model of large load areas for harmonic studies in distributions networks”, IEEE Transactions on Power Delivery, 12.1 (1997) [3] R. Carbone, D. Menniti, R. E. Morison. and A. Testa, “Harmonics and inter-harmonic distortion modeling in multi-converter system,” IEEE Transaction on Power Delivery, 10.3(1995):1685-92. [4] M. K. Ghartemani, M. R. Iravani, “Measurement of Harmonic/ Inter-Harmonics of Time Varying Frequencies”, IEEE Transaction on Power Delivery, 20.1(2005):23-31. [5] C. Li, W. Xu, T. Tayjasanant, “A critical impedance based method for identifying harmonic sources”, IEEE Transaction. On Power Delivery,. 19.2(2004). [6] C. Li, W. Xu and T. Tayjasanant, "Interharmonics: Basic Concepts and Techniques for Their Detection and Measurement", Electric Power Systems Research,. 66.1(2003):39-48. [7] R. Ortega, A. Loria, P. J. Nickasson, and H. Sira-Ramyrez. “Passivity-based control of Euler-Lagrange systems”, Communications and Control Engineering. Springer-Verlag, Berlin, (1998). [8] C. Cecati. “Experimental investigation of the passivity based induction motor control”. Industrial Electronics, 1998. [9] . Gokdere, Levent U., and Marwan A. Simaan. " “A passivity-based method for induction motor control”.. Industrial Electronics, IEEE Transactions on 44.5 (1997): 688-95. [10] P. J. Nicklasson, R. Ortega, and G. Espinosa. “Passivity based control of a class of blondel- park transformable machines”. Automatic Control, IEEE Transactions on 42.5 (1997): 629- 647. [11] W. Xu, H.W. Dommel, and M. B. Hughes, "Modeling of Adjustable Speed Drives for Power System Harmonic Analysis", IEEE Transaction. on Power Delivery,. 14.2(1999):595-601. [12] “Definitions for the Measurement of Electric Quantities under sinusoidal, non sinusoidal, balanced or unbalanced conditions”, IEEE Std. 1459-2000, Jan. 2000D. [13] G. Dorrel, W. T. Thomson and S. Roach, “Analysis of airgap flux, current and vibration signals as a function of the combination of static and dynamic airgap eccentricity in 3-phase Induction motors”, IEEE Transaction. on Industry. Applications ,. 33(1997):24-34. [14] Jee-Hoon Jung, Jong-Jae Lee and Bong-Hwan Kwon, “Online Diagnosis of Induction Motors using MCSA”, IEEE Transactions on Industrial Electronics, 53.6(2006):1842-52. [15] Neelam Mehla and Ratna Dahiya, “An Approach of Condition Monitoring of Induction Motor Using MCSA”, International Journal of Systems Applications, Engineering and Development 1.1(2007):13-17. [16] Benbouzid, M.E.H., “A review of Induction motor Signature analysis as a medium for Fault detection”, IEEE Transaction Industrial. Electron., 47.5(2000):984-993. [17] Jordi Cusido, Luis Romeral, Juan A. Ortego, Javier A. Rosero and Antonio Garefa Espinosa, “Fault Detection in Induction Machines Using Power Spectral Density in Wavelet Decomposition”, IEEE Transactions on Industrial Electronics, 55.2 (2008):633-643. [18] Jee-Hoon Jung, Jong-Jae Lee and Bong-Hwan Kwon, “Online Diagnosis of Induction Motors using MCSA”, IEEE Transactions on Industrial Electronics, 53.6(2006):1842-52. [19] J. Cusido, A. Jornet, L. Romeral, J.A. Ortega and A. Garcia, “Wavelet and PSD as a Fault Detection Technique”, ITMC 2006-Instrumentation and Measurement Technology Conference, Proceeding of IEEE, Sorrento, Italy, 24-27 April 2006.. 1397(2006). [20] Fatiha Zidani, Mohamed El Hachemi Benbouzid, Demba Diallo, and Mohamed Saïd Naït- Saïd, “Induction Motor Stator Faults Diagnosis by a Current Concordia Pattern-Based Fuzzy Decision System”, IEEE Transaction on Energy Conversion, 18.4(2003):469-75. [21] Demba Diallo, Mohamed El Hachemi Benbouzid, Denis Hamad, and Xavier Pierre, “Fault Detection and Diagnosis in an Induction Machine Drive: A Pattern Recognition Approach IIUM Engineering Journal, Vol. 14, No. 2, 2013 Chattopadhyay and Sengupta 161 Based on Concordia Stator Mean Current Vector”, IEEE Transaction on Energy Conversion,. 29.3(2005):512-19. [22] M. E. H. Benbouzid “Monitoring and diagnosis of Induction motor electrical faults using a current Park's vector pattern approach”, Proceeding. of IEEE International. Conference. On Electrical Machines and Drives, Seattle, WA, (1999): 275 – 77. [23] H. Nejjari “Monitoring and diagnosis of induction motors electrical faults using a current Park’s vector pattern learning approach”, IEEE Transaction. Industrial. Application., 36.3(2000):730-35. [24] S. Chattopadhyay, M. Mitra, and S. Sengupta: “Harmonic Analysis in a Three-Phase System using Park Transformation Technique”, International Journal on Modeling, Measurement and Control of General Physics and Electrical Applications, AMSE, . 80.3 (2007):42-58. [25] Chattopadhyay, Surajit, Madhuchhanda Mitra, and Samarjit Sengupta “Electric Power Quality”, Springer, 2011. [26] S Chattopadhyay, S Karmakar, M Mitra, S Sengupta, “Symmetrical Components and Current Concordia based Assessment of Single Phasing of an Induction Motor by Feature Pattern Extraction Method and Radar Analysis”, International Journal of Electrical Power and Energy Systems,. 37.1(2012):43-49. [27] S Chattopadhyay, S Karmakar, M Mitra, S Sengupta, “Assessment of crawling of an induction motor by stator current Concordia analysis”, IET Electronic. Letters. -- . 48.14(2012):841-42.