A DETAILED STUDY OF THE HIGH ORDER SERIAL RESONANT INVERTER FOR INDUCTION HEATING

Similar documents
LC transfer of energy between the driving source and the circuit will be a maximum.

RESONANCE SERIES RESONANT CIRCUITS. 5/2007 Enzo Paterno 1

Power efficiency and optimum load formulas on RF rectifiers featuring flow-angle equations

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Conventional Paper-I (a) Explain the concept of gradient. Determine the gradient of the given field: ( )

15 B1 1. Figure 1. At what speed would the car have to travel for resonant oscillations to occur? Comment on your answer.

Dynamic Performances of Self-Excited Induction Generator Feeding Different Static Loads

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

Duality between Statical and Kinematical Engineering Systems

Experiment I Voltage Variation and Control

Analysis of high speed machining center spindle dynamic unit structure performance Yuan guowei

Steady State and Transient Performance Analysis of Three Phase Induction Machine using MATLAB Simulations

Physics 313 Practice Test Page 1. University Physics III Practice Test II

ME 3600 Control Systems Frequency Domain Analysis

J. Electrical Systems 1-3 (2005): Regular paper

Swissmetro: design methods for ironless linear transformer

Current, Resistance and

TECHNICAL SIMULATIONS AND MEASUREMENTS OF CAGE-SWITCHING AT THE INDUCTION MACHINE. Tulbure Adrian*, Risteiu Mircea**

Reading Assignment. Problem Description for Homework #9. Read Chapters 29 and 30.

[2007] IEEE. Reprinted, with permission, from [Jiaxin Chen, Jianguo Zhu, Youguang Guo, A 2-D nonlinear FEA tool embedded in Matlab/Simulink

Chapter 5 Force and Motion

Electrodynamic Forces in Steel Strip during Induction Heating

Chapter 5 Force and Motion

Three-dimensional Quantum Cellular Neural Network and Its Application to Image Processing *

STABILITY AND PARAMETER SENSITIVITY ANALYSES OF AN INDUCTION MOTOR

EKT 345 MICROWAVE ENGINEERING CHAPTER 2: PLANAR TRANSMISSION LINES

Tutorial Exercises: Central Forces

Chapter 33 Alternating Current

2.5 The Quarter-Wave Transformer

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS

Phys 1215, First Test. September 20, minutes Name:

Lecture 2 Date:

A new approach in classical electrodynamics to protect principle of causality

Bifurcation Analysis for the Delay Logistic Equation with Two Delays

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

Using Laplace Transform to Evaluate Improper Integrals Chii-Huei Yu

EXAM NMR (8N090) November , am

International Journal of Mathematical Archive-3(12), 2012, Available online through ISSN

Analysis and Optimization of a Special Type of Dielectric Loaded Resonant Cavity for Mobile Communication Filters

AP Physics C: Electricity and Magnetism 2003 Scoring Guidelines

4. Electrodynamic fields

Electricity Revision ELECTRICITY REVISION KEY CONCEPTS TERMINOLOGY & DEFINITION. Physical Sciences X-Sheets

A Three-Dimensional Magnetic Force Solution Between Axially-Polarized Permanent-Magnet Cylinders for Different Magnetic Arrangements

DonnishJournals

ASTR415: Problem Set #6

(Sample 3) Exam 1 - Physics Patel SPRING 1998 FORM CODE - A (solution key at end of exam)

ITI Introduction to Computing II

A method for solving dynamic problems for cylindrical domains

Physics 2020, Spring 2005 Lab 5 page 1 of 8. Lab 5. Magnetism

Basic Bridge Circuits

Explosive Contagion in Networks (Supplementary Information)

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

Interaction of Feedforward and Feedback Streams in Visual Cortex in a Firing-Rate Model of Columnar Computations. ( r)

Physics NYB problem set 5 solution

(a) Unde zeo-bias conditions, thee ae no lled states on one side of the junction which ae at the same enegy as the empty allowed states on the othe si

COUPLED MODELS OF ROLLING, SLIDING AND WHIRLING FRICTION

Physics 221 Lecture 41 Nonlinear Absorption and Refraction

Qualifying Examination Electricity and Magnetism Solutions January 12, 2006

HOW TO TEACH THE FUNDAMENTALS OF INFORMATION SCIENCE, CODING, DECODING AND NUMBER SYSTEMS?

Rotor Flux Estimation of Induction Motors Using Sliding-Mode Observer

Exam 3, vers Physics Spring, 2003

Computers and Mathematics with Applications

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS

A Simple Method to Control of Variable Speed Wind Generation System Coupled with Squirrel Cage Induction Generator

Gravitational Radiation from Oscillating Gravitational Dipole

EM-2. 1 Coulomb s law, electric field, potential field, superposition q. Electric field of a point charge (1)

Calculation of Temperature Rise in Dry-type Air-core Reactors Using Strong Coupling of Fluid-Temperature Field

This brief note explains why the Michel-Levy colour chart for birefringence looks like this...

M.M.A HIGHER SECONDARY SCHOOL - PAPPANADU 12 TH PHYSICS - HALF YEARLY EXAM ANSWER KEY

MSE 561, Atomic Modeling in Material Science Assignment 1

SEARCH-BASED DYNAMIC IDENTIFICATION OF INDUCTION MOTORS

APPLICATION OF MAC IN THE FREQUENCY DOMAIN

Graphs of Sine and Cosine Functions

11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.

Contact impedance of grounded and capacitive electrodes

Research on the Method of Identifying the Indicator Diagram of the Sucker Rod Pumping Unit

A scaling-up methodology for co-rotating twin-screw extruders

Proceedings of the 11th WSEAS International Conference on Automatic Control, Modelling and Simulation

Ferrite Loss Measurement and Models in Half Bridge and Full Bridge Waveforms

Functions Defined on Fuzzy Real Numbers According to Zadeh s Extension

Galilean Transformation vs E&M y. Historical Perspective. Chapter 2 Lecture 2 PHYS Special Relativity. Sep. 1, y K K O.

4/18/2005. Statistical Learning Theory

Conservative Averaging Method and its Application for One Heat Conduction Problem

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Measurement of Electrical Properties of Rapeseed Seeds with LCR Meter Good Will 8211

Improvement in Accuracy for Design of Multidielectric Layers Microstrip Patch Antenna

Compactly Supported Radial Basis Functions

EFFECTS OF FRINGING FIELDS ON SINGLE PARTICLE DYNAMICS. M. Bassetti and C. Biscari INFN-LNF, CP 13, Frascati (RM), Italy

Fields and Waves I Spring 2005 Homework 8. Due: 3 May 2005

Aalborg Universitet. Load Estimation from Natural input Modal Analysis Aenlle, Manuel López; Brincker, Rune; Canteli, Alfonso Fernández

Numerical Integration

Physics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!

Physics: Work & Energy Beyond Earth Guided Inquiry

FI 2201 Electromagnetism

EKT 356 MICROWAVE COMMUNICATIONS CHAPTER 2: PLANAR TRANSMISSION LINES

FREE SPACE OPTICS - MEASUREMENT OF TRANSMISSION QUALITY LINK PARAMETERS

Phys-272 Lecture 18. Mutual Inductance Self-Inductance R-L Circuits

A High Performance Rotor-Flux Oriented Indirect Vector Controlled Induction Motor Drive using Fuzzy Logic Control

CHAPTER 2 DERIVATION OF STATE EQUATIONS AND PARAMETER DETERMINATION OF AN IPM MACHINE. 2.1 Derivation of Machine Equations

Phys 222 Sp 2009 Exam 1, Wed 18 Feb, 8-9:30pm Closed Book, Calculators allowed Each question is worth one point, answer all questions

Transcription:

ELECTRONICS 005 1 3 Septembe, Sozopol, BULGARIA A DETAILED STUDY OF THE HIGH ORDER SERIAL RESONANT INVERTER FOR INDUCTION HEATING Evgeniy Ivanov Popov, Liliya Ivanova Pindeva, Elisaveta Histova Mileva, Boislav Nikolaev Tsankovski Depatment of Powe Electonics, Technical Univesity of Sofia, Kl. Ohidski Bul. No.8, 1000 Sofia, Bulgaia, phone: 359 96564, e-mail: epopov@tu-sofia.bg The eal high ode powe cicuit of the seial esonant invete without fee- wheeling diodes widely applied fo induction heating and melting is analyzed assuming a sine-wave shape of the load voltage. The load cicuit is eplaced by seially connected equivalent esistance and eactance. The paametes of the invete cicuit, especially the load voltage, the voltage acoss the powe semiconducto devices and in the thyisto vaiant the cicuit tun-off time ae detemined. A MATLAB pogam calculates and gaphically displays the fequency contolling chaacteistics. A diect steady state simulation of the invete opeation confims the esults. The equiements towads the closed loop automatic system of fequency contol can be easily defined. Keywods: seial invete, esonant invete, high-ode cicuit, chaacteistics 1. INTRODUCTION The seial esonant invete without fee - wheeling diodes applies widely fo induction heating and melting due to its convenience fo high fequency pefomance. The invete may be constucted with seially connected fully contollable devices (tansistos o GTO thyistos) and diodes. Fo cheape and not so high fequency applications when thee is a negative voltage acoss the devices ight away afte thei conduction inteval the powe cicuit can be built with thyistos. The seial esonant invete [1] has not been fully analyzed in the pactical case when the invete cicuit is of high ode. That s why this aticle is aimed at investigation of the high ode powe cicuit of the seial esonant invete fo induction heating and obtaining the fequency contol chaacteistics of the main paametes. That will allow the powe cicuit to be coectly designed taking into account the lage vaiations of the load and the estictions of the powe cicuit. The equiements towads the contolling system can also be defined. Simulations and possibly expeiments should confim the esults.. MAIN ASSUMPTIONS A half bidge cicuit (Fig. 1) is unde study but it can be easily conveted into the bidge one. The powe losses in the invete ae neglected. The commutation of the powe semiconducto devices is instantaneous. A paallel equivalent cicuit epesents the induction heate. The quality facto of the load cicuit is sufficiently high that the voltage acoss the load has a close to the sine wave shape. A combination of the method of tansitoy values and the Fouie tansfomation is applied in the analysis.

ELECTRONICS 005 1 3 Septembe, Sozopol, BULGARIA 3. REPRESENTATION OF THE LOAD CIRCUIT A paallel equivalent cicuit containing the active esistance R l and the inductance L l that ae fequency dependent epesents the inducto heate: n (1) R ( ) 1 l R n 1 () ( ) L l L, whee is the esonant fequency of the load cicuit R l, Ll, Cl 1 (3) Cl L R, L ae the inducto paametes fo the esonant fequency. π f is the contolling fequency. The values of n 1 and n ae given in []. If the quality facto of the load cicuit R l, Ll, Cl at esonance (index ) C (4) Q l R L is high enough the load voltage is almost sine wave and the load cicuit can be epesented by a seial equivalent cicuit containing the paametes R eq, X eq (Fig.). n ( ) 1 Req (5) R (n 1) 1 ( ) 1 n Q ( 1 n ( n1 ) Q ( ) ( X eq (6) R (n 1) 1 ( ) 1 n Q ( Below the esonance of the load cicuit the equivalent eactance X eq is in fact a seial equivalent inductance 1 n ( ) n1 Q ( ) ( 1) Leq (7) L (n 1) 1 ( ) 1 n Q ( Theefoe the seial invete cicuit consists of R Req, L L eq LS, C CS. Above the esonance of the load cicuit the equivalent eactance X eq is a seial equivalent capacitance

ELECTRONICS 005 1 3 Septembe, Sozopol, BULGARIA (n 1 ( ) 1 n 1) Q ( ) C eq (8) C 1 n l (n ) 1 Q ( ) ( 1) CeqCS and the seial invete cicuit consists of R Req, L LS, C, Ceq CS α C /C eq. Then the voltage acoss the invete capacito C is divided between the capacitos C S and C eq in accodance with the theoy of electical engineeing [3]. Gaphical epesentations of the functions (5), (6), (7) and (8) ae given in Fig., Fig. 3, Fig. 4 and Fig. 5 espectively. 4. PARAMETERS OF THE SERIAL INVERTER CIRCUIT The whole fouth ode invete cicuit is educed to a second ode RLC esonant ( R eq < L / C ) cicuit. This invete is analyzed in [4], [5], [1], [6]. 5. CALCULATION OF THE LOAD VOLTAGE The RMS of the load voltage can be detemined fom its hamonic components as follows (11) V l V l ( m) m 1,3,5,... Each hamonic component of the load voltage is c( m) (1) V l( m) Inom.. R eq( m) X eq( m) whee R eq(m) and X eq(m) ae calculated fom (5) and (6), exchanging / with m /. I nom Vd / Req fo natual commutation ( 0 ) and Inom Vd /( 0L) fo foced commutation ( > 0 ) whee 0 1/( LC) Req /(4L ) is the inheent fequency of the seial invete cicuit and πf is the contolling fequency. The nomalized peak value of the m hamonic component (m1,3,5, ) of the invete cuent c ( m) a ( m) b ( m) is given in [1] fo natual commutation and in [7] fo foced commutation. 6. VOLTAGE ACROSS THE POWER SEMICONDUCTOR DEVICES At natual commutation duing the inteuption inteval t X1 0 to tp π / π / 0 the voltage of the tuned-off devices is c(1) Vd (13) VVS 0.5Vd m ( k 0.5)(1 α ) Vd ± sin( tx1 tp ϕ(1) ϕeq ) cosϕeq

ELECTRONICS 005 1 3 Septembe, Sozopol, BULGARIA whee the uppe sign coesponds to the just tuned-off devices, ϕ ( 1) actg( a(1) / b(1) ), ϕ eq actg( X eq / Req ). When the othe device conducts t X 0 to π / 0 the voltage of the tuned-off device is δ t δ V 0.5 (1 ) [0.5 VS Vd α Vd ke X ( sin0t X cos0t X )] 0 (14) c(1) Vd sin( tx ϕ(1) ϕeq ) cosϕeq R eq 1 whee δ is the damping coefficient and k is the coefficient of L π δ 1 e 0 hesitation. Expessions (13) and (14) ae applied fo calculation of the maximal device voltage and in the thyisto vaiant the cicuit tun-off time. Fo foced commutation the necessay expession is given in [6]. 7. FREQUENCY CHARACTERISTICS A MATLAB pogam pocesses all the mathematical infomation descibing the steady-state opeation of the invete togethe with its load cicuit in the allowed fequency ange. The fequency chaacteistics of the invete ae obtained and gaphically displayed in Fig. 6. These paticula chaacteistics coespond to a pactically implemented invete with the following data: V d 500 V; L S 45 µh; C S 84 µf; R 0.4739599 Ω, f 083 Hz, L 8.8717 µh, C l,657.88 µf,n 1 0, n 1, f1600-600 Hz. The fist set of gaphics shows: two times the aveage input cuent of the invete I d [A] (solid line), two times the load voltage V l [V] (dashed line), ten times the cicuit tun-off time of the thyistos 10 t q.c. [µs] (dotted line). The second set of gaphics shows: the peak seial capacito voltage V CSm [V] (dash-dotted line), the maximal device voltage V VSm (solid line). The same MATLAB pogam can calculate fequency chaacteistics fo diffeent input data. Othe paametes of the cicuit can also be calculated and displayed. The calculation is moe pecise if the shape of the load voltage is close to the sine-wave fom. That is checked fo each paticula fequency. The stability and efficiency of the invete can be studied fom the chaacteistics. In geneal thee is a minimum of the input cuent and powe, load voltage, seial capacito voltage and device voltage aound the esonance of the load cicuit. If the paametes of the inducto heate vay duing the induction heating pocess it is advisable to maintain almost constant input cuent (powe) in a slight capacitive deangement of the load cicuit R l, L l, C l, whee the opeation is stable, by execising an influence on the contolling fequency. The calculated slope of the fequency chaacteistic of the input cuent is used fo detemining the paametes of the closed loop automatic contol system [1].

ELECTRONICS 005 1 3 Septembe, Sozopol, BULGARIA 8. DIRECT SIMULATION OF THE STEADY STATE The steady state is simulated by a MATLAB pogam based on the method descibed in [8] that is expeimentally confimed. The esults fom the fequency chaacteistics and fom the simulation ae in good ageement accoding to Table 1. That confims the coectness of the whole study. The detailed simulation esults of the powe invete ae gaphically displayed (Fig. 7) fo f195 Hz. The fist set of diagams shows: the seial capacito voltage v CS [V] (solid line), invete cuent i CS [A] (dashed line) and twice the load voltage v l [V] (dotted line). The second set of diagams shows: a quate of the load inductance cuent i Ll [A] (dotted line) and the voltage acoss a device v VS1 [V] (solid line). The simulation esults fo the same cicuit at f694,4 Hz displaying an inefficient but possible opeation in an inductive deangement of the load geneally with the thid hamonic component of the invete cuent (a case that is difficult to be studied analytically) ae given in Fig. 8. Table 1. f[hz] 694,4 183 018 083 Recom. 195 394 Study Sim. F.ch. Sim. F.ch. Sim. F.ch. Sim. F.ch. Sim. F.ch. Sim. Load Ind. Ind. Ind. Res. Cap. Cap I d [A] 169,5 400,7 388,7 00 196,5 178, 180,1 199,8 08, 400,6 416,3 V l [V] 00,4 304,4 303,5 15,4 15,8 01,7 06,6 11,9,1 30, 314,1 t q.c. [µs] 604 113,6 15,9 61 68,3 63,8 7,5 8,9 89,9 97,3 99,6 V CSm [V] 1456 1308 170 590 580 509 515 54 565 996 1034 V VSm [V] 197 134 134 866 875 863 886 1003 1045 1636 1689 9. CONCLUSION A detailed study of the eal seial esonant invete without fee - wheeling diodes fo induction heating is pefomed. The steady state paametes and chaacteistics ae detemined. The equiements fo the contol system can be defined. 10. REFERENCES [1] Popov E., Automatic Contol of Seial Thyisto Invetes with Vaying Load, Ph.D. Thesis, VMEI, Sofia, 1978 (in Bulg.). [] E. I. Bekovitch, G. V. Ivenskiy, Yu. S. Yoffe, A T. Matchak, V. V. Mogun, Highe fequency thyisto convetes fo electo technology applications, Leningad, Enegoatomizdat, 1983 (in Russian), pp. 8-13. [3] G. I. Atabekov, Fundamentals of cicuit theoy, Moscow, Enegy, 1969 (in Russian), pp.131-13, 100. [4] Nachev N., G. Maleev, Powe Electonics, Technica, Sofia, 1979 (in Bulg.). [5] Kusteva A., Investigation of the system seial thyisto invete induction heate, Ph.D. Thesis, VMEI, Sofia, 1974 (in Bulg.). [6] Popov E., Investgation and Analysis of Seial Thyisto RLC Invetes, Jounal Electical Engineeing and Electonics EE,vol.7-9, 003 (in Bulg.), pp.39-43. [7] Popov E., A. Popova, Study of Seial Thyisto RLC Invetes. Nomalized Examples, Poceedings of the Eleventh National Scientific and Applied Science Confeence with Intenational Paticipation ELECTRONICS ET 00, Sozopol, Bulgaia, Septembe 00 (in Bulg.). [8] Popov E., Diect detemination of the steady-state mode in autonomous invetes, Scientific Jounal Electical Engineeing and Electonics EE, vol. 11-1, 004, Sofia, Bulgaia, (In Bulg.).

ELECTRONICS 005 1-3 Septembe, Sozopol, BULGARIA i d VS1 V d / L l R l i L S1 o V d / C l C S L S v l VS Fig. 1. Fig. 5. Fig.. Fig. 6. Fig. 3. Fig. 7. Fig. 4. Fig. 8.