Citation
Impact of electric vehicle loads on electric power distribution systems

Material Information

Title:
Impact of electric vehicle loads on electric power distribution systems
Creator:
Buchh, Tariq Aslam
Publication Date:
Language:
English
Physical Description:
vii, 100 leaves : ill. ; 29 cm.

Subjects

Subjects / Keywords:
Batteries ( jstor )
Customers ( jstor )
Electric current ( jstor )
Electric potential ( jstor )
Electric power ( jstor )
Harmonics ( jstor )
Recordings ( jstor )
Simulations ( jstor )
Transformers ( jstor )
Waveforms ( jstor )
Dissertations, Academic -- Electrical and Computer Engineering -- UF ( lcsh )
Electric power consumption -- Forecasting ( lcsh )
Electric power distribution ( lcsh )
Electric vehicles ( lcsh )
Electrical and Computer Engineering thesis, Ph. D ( lcsh )
Genre:
bibliography ( marcgt )
non-fiction ( marcgt )

Notes

Thesis:
Thesis (Ph. D.)--University of Florida, 1998.
Bibliography:
Includes bibliographical references (leaves 92-99).
General Note:
Typescript.
General Note:
Vita.
Statement of Responsibility:
by Tariq Aslam Buchh.

Record Information

Source Institution:
University of Florida
Holding Location:
University of Florida
Rights Management:
Copyright [name of dissertation author]. Permission granted to the University of Florida to digitize, archive and distribute this item for non-profit research and educational purposes. Any reuse of this item in excess of fair use or other copyright exemptions requires permission of the copyright holder.
Resource Identifier:
029532204 ( ALEPH )
41284352 ( OCLC )

Downloads

This item has the following downloads:


Full Text









IMPACT OF ELECTRIC VEHICLE LOADS ON ELECTRIC POWER
DISTRIBUTION SYSTEMS













By

TARIQ ASLAM BUCHH


A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL
OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT
OF THE REQUIREMENTS FOR THE DEGREE OF
DOCTOR OF PHILOSOPHY


UNIVERSITY OF FLORIDA


1998














In the Name of God, the Most Beneficent, the Most Merciful.

All the praises and thanks be to God, the Cherisher and Sustainer of the Worlds.
The Most Beneficent, the Most Merciful.
The Only Master of the Day of Recompense.
You (Alone) we worship and You (Alone) we ask for help.
Guide us to the Straight Way.
The way of those on whom You have bestowed Your Grace and not of those who
earned your anger, nor of those who went astray.


Translation of Opening Chapter of the Holy Ouran













ACKNOWLEDGEMENTS


I acknowledge the conviction of those who believe in the Truth and have the

courage to stand for it. Having said that, I would like to mention the names of those who

have been there when I needed them. A word of appreciation for my adviser Dr. Alex

Domijan Jr. for agreeing to chair my advisory committee and to Dr. Dennis Carroll, Dr.

Khai D Ngo, Dr. Herman Lam and Dr. Barney Capehart for being on the committee.

My mother, who has been a constant source of inspiration and drive, my father,

who is an embodiment of patience, my wife, Falaknaz, who sacrificed a part of her

dreams and aspirations and is taking excellent care of our kids, Shaima and Ziad, my

loving brothers, Mehmood and Khalid and my best friend, Ashfaq. There has also been a

huge treasure of friends and well-wishers who were there to support and pray for me

whenever the going got tough. A special word of thanks for Dr. Song for his constructive

criticism and suggestions.

Florida Power and Light Company, provided the funding and technical support

and I appreciate it.















TABLE OF CONTENTS



ACKNOWLEDGEMENTS.. ...............................................................iii

ABSTRACT ............................ .......................... ..................vi

CHAPTERS

1 INTRODUCTION ................................... ..................... ..... .......I
1.1. Background................................................. 1
1.2. Project Justification .................................................. 3
1.3. Objectives and Goals..................................... ....... ... .................. 4
1.4. M ajor Contributions ...................................................................... 4
1.4. Research Phases....................................................... 5
2 LITERATURE SURVEY............................... ...... ............................. 7

3 POWER QUALITY AND HARMONIC DISTORTION ..................................... 10
3.1. Harm onics........................................................... 11
3.2. Estimation of Harmonic Distortion............................. ...................... 12
3.2.1. C rest Factor ..................................................... ...................... 12
3.2.2. Percentage Total Harmonic Distortion (THD).................. 13
3.2.3. The K-Factor ............................................. 13
3.3. Transform er Derating .................................................................... 14
3.3.1. Transformer Capability Equivalent Calculation ...................... 15
3.3.2. Effect of Current Harmonics on the Transformer Neutral Current ....... 17
4 FIELD DATA RECORDING, REPRESENTATION AND OBSERVATION........... 19
4.1. Introduction. .................................... ................ ......... .................... 19
4.2. General Considerations In Data Acquisition And Processing...................... 20
4.2.1. Data Collection...................................................... 21
4.2.2. Data Recording........................................................ 22
4.2.3. Data Preparation..................................................... 22
4.2.4. D ata Q ualification ............................................... .............. 23
4.2. Data Analysis..................................................... 24
4.3. Field Data M monitoring ............................................................. 26
4.3.1. Data Acquisition Methodology...................... ....... 26
4.3.2. Data M monitoring Instrument......................... .......... ... .................... 26
4.3.3. Data M management ........................................................... 27


iv








5 DATA REPRESENTATION AND ANALYSIS..................... ....... 30
5.1. Introduction ............................................................ 30
5.2. Data Trends and Characteristics.......................... ........... ...................... 32
5.3. Data Observation ............................................... .................. 33
5.4. Load Characteristics ......................................................... 34
6 LABORATORY TESTING OF THE EV CHARGER........................................... 47
6.1. Testing Procedure and Results .............................................. ....... 47
6.2. Observation of Laboratory Results................................... 48
7 COMPUTER SIMULATION.................................................. 61
7.1. Component Models............................... .......................... 61
7.1.1. System M odel........................................................... ..................... 6 1
7.1.2. Electric Vehicle Charger And Battery Models .................................. 64
7.1.3. Transformer Model .......................................... 64
7.1.4. Domestic Load M odel.................................. .............. 66
7.2. Mathematical Model Of The EV Charger....................... 66
6.3. Simulation Results................................. .......................68
7.4. Simulation of Short Charging Time Chargers....................................... 69
7.4.1. Simulation Methodology For Fast Chargers..................................... 70
8 COMPARISON OF FIELD, LABORATORY AND SIMULATION RESULTS....... 73
8.1. The Field D ata ............................................................ ........... 74
8.2. Laboratory Test....................................... ............74
8.3. Simulation Results...................... ........ ........ .......... 75
8.4. Comparison.................................................. 75
9 DERATING OF DISTRIBUTION TRANSFORMERS..................................... 83
9.1 M ethod of Calculation................................................... 83
9.1.1 ANSI/IEEE C 57.110 M ethod............................... .................... 83
9.1.2 Modified ANSI/IEEE C 57.110 Method.................................. ..... 85
9.2. Transformer Derating with Fast Time Chargers................................... 85
9.3. Conclusion.......................... .... ........................ 86
10 CONCLUSIONS.................................................. 89
10.1 .Summary and Conclusions................... ................... 89
10.2. Scope For Further Research ................... .................... 90

REFERENCES................................ .. ..........92

BIOGRAPHICAL SKETCH................... ............... 100













Abstract of Dissertation Presented to the Graduate School
of the University of Florida in Partial Fulfillment of the
Requirements for the Degree of Doctor of Philosophy



IMPACT OF ELECTRIC VEHICLE LOADS ON ELECTRIC POWER
DISTRIBUTION SYSTEMS

By

Tariq Aslam Buchh

August 1998


Chairman: Dr. Alexander Domijan, Jr.
Major Department: Electrical and Computer Engineering

The voltages and currents on four distribution sites (FPL_1 to FPL_4) have been

monitored. The load behavior with and without the EV load has been recorded. By

performing laboratory tests on the GMC EVI charger charging a fully discharged battery,

the characteristics of the charger have been determined. The harmonic spectrum of the

transformer current charging a fully discharged battery has been determined. The snap

shots of the current and voltage waveforms have been recorded and studied. Based on the

laboratory data a computer model of the charger and the distribution system has been

developed. Comparing the simulation results to the laboratory results and the field data

the simulation results have been validated. The simulated currents are roughly 4 to 10%

away from the field recorded currents, in their magnitude and distortion. Based on the

field data, laboratory results and the simulation results the derating of the distribution








transformer has been calculated. The two methods suggested for the calculation of

derating have been presented and compared. The maximum and minimum limits of

derating have been suggested. Using the first method the maximum derating is 20% and

the minimum value is 99%. Using the modified method the derating remains close to 99

% throughout the charging cycle. (It must be noted here, that a derating to 20% means

that the transformer can safely handle only 20% of its full rated load. Likewise a derating

to 99% means that the transformer can handle 99% of its full rated capacity. Thus, a

transformer derated to 20% is highly derated and a transformer derated to 99% is slightly

derated.) After establishing the reliability of the simulation, by comparing the simulation

results with the field data and laboratory results, a charging scenario with shorter

charging times (45 minutes and 1.2 hours) has been simulated. The distribution

transformer derating corresponding to these short-charging times has been calculated and

presented. Furthermore, recommended scope for further research has been suggested

which among other things, involves validating the derating formula for distribution

transformers.


VII













CHAPTER 1
INTRODUCTION


1.1. Background


From the state of being a concept vehicle, a possible means of future

transportation and a dream, the Electric Vehicle (EV) has become a reality. But even after

the first major commercial launch of the EVI by General Motor Corporation, the future

of EVs is still under a cloud of uncertainty and some unanswered questions.

The popularity of EVs as a future mode of transport can be attributed partly to

exhausting reserves of natural energy resources and environmental considerations. These

and many other factors have boosted the research and development of the EV technology

and a possibility of it being accepted as the popular mode of ground transport in the not

so distant future. There has been a great boost to EV research because of government and

industrial funding. This research has been mainly aimed at improving the battery

technology, charging techniques, weight, comfort and reliability. For the

environmentalist, the EVs' noiseless operation compliments its pollution free

performance, and, for the customer, the EV is a good performance machine with the

downside of short battery life, long recharging time and high initial cost.

The EVI manufactured by General Motors Corporation is powered by a Delco Valve-

Regulated Lead Acid battery pack, which contains 26 12-volt modules, capable of

carrying 16.2-kilowatt hours of energy. The batteries can be recharged from 0 to a full








state of charge in approximately 3 hours using a 240-volt inductive Magne charge

charger and about 15 hours from a 120-volt portable charger. Inductive charging is

achieved without metal-to-metal contact. The alternating current powers a device that

produces a fluctuating magnetic field in a paddle that is inserted into a port at the front of

the car. This induces an alternating emf in a winding inside the car. This alternating emf

is then rectified and fed into the batteries. The inductive charging system is inherently

safe even under hazardous climatic conditions, since the windings are safely encased in a

plastic cover. Typically this system can handle power levels from 1.5 to 25 kW with

overall efficiency of 90% while the power transfer frequency is 40-350 kHz [1]. The

specifications of the EV, as published by General Motors [PrEView Drive program

literature], are summarized as follows:

1. 137 H.P. three-phase ac induction motor.

2. 16.8 kWh lead-acid battery.

3. Inductive coupled standard charger (220 volts, 30 Amps)

4. IGBT power inverter module.

5. Range: 90 miles highway, 70 miles EPA city.

6. Acceleration from 0 to 60 mph in 8.5 seconds.

7. Electronically regulated speed of 80 mph.

8. Charging time from 85% depth of discharge is 2 to 3 hours.

9. Curb weight is 2970 pounds.

While the popularity and widespread use of EVs on one hand will take ground

transportation into a new era, on the other hand it can put additional and unprecedented strains

on power supply and distribution systems.








For power systems, the EV loads pose twofold problems. These are:

1. Increase in power demand

2. Decrease in power quality.

The power distribution system looks at the EV as an additional load of about 7kw/EV

(based on the GMC EV ), which was almost certainly not foreseen when most of the present

distribution systems were designed more than a few decades ago. The fact that the EV charger

load is a non-linear load with power electronic devices makes matters worse, for it is a source

of harmonics and distortion in the distribution system. These harmonics which get injected

into the system have a significant impact on the distribution system in terms of losses and

overall performance of the power system. The non-linearities and harmonics caused by EV

chargers can adversely affect other sensitive loads that may be connected to the distribution

system. The power quality can become a major concern if the battery and charger technology

is developed such as to put an extra strain on the distribution system. For example, research is

now being aimed at developing batteries and their chargers which would make it possible to

charge a completely exhausted battery in a very short period of time (15 minutes).




1.2. Project Justification




Existing distribution systems may not have the capability to cater to the new

needs of customers if the use of EVs becomes widespread. The electric power utility

companies need to be aware of and prepared to face such load scenarios. The utility

companies need to first have a reliable assessment of the possible future load scenarios

with EV loads. After this, they need to have a ready design of the new distribution system








equipment in terms of their sizes and ratings and also in terms of the configuration or the

layout of the entire distribution system. Thus, in case of an increase in the use of EVs by

the customers, the utilities will find themselves well prepared and thus continue to

reliably supply quality power to their customers.




1.3. Objectives and Goals




The aim of this project is to provide Florida Power and Light Company with an

assessment of the possible load effects with EV loads. The outline of the objectives of the

project is given below.

* To study the normal load behavior for a typical distribution system feeding two to ten
houses.

* To study the load behavior when one of the houses in the above case owns an Electric
Vehicle.

* To perform laboratory tests to determine the behavior of the General Motor EVI
charger, in terms of its current requirements and harmonic effect on the input line.

* Based on the above data to prepare a computer model of the distribution system
feeding loads consisting of EV loads.

* Prepare the mathematical model for various components and the whole system.

* Perform basic statistical analysis of the recorded data.

* Use computer simulations to study various possible load scenarios.



1.4. Maior Contributions


Major contributions of this dissertation work are itemized below.









Actual field data was recorded from distribution transformer sites. This data
contains information about the load behavior and characteristics at these
particular sites. It has been for the first time, that such data has been made
available, for load scenario comparison between loads with and without the
EV loads.

Harmonic analysis has been performed on the field-recorded data, to identify
the signature that the EV load leaves on the system. This has provided with an
understanding of the EV load behavior when put on the system along with
other domestic loads.

Actual laboratory testing of the EV charger, charging into a fully discharged
battery has been carried out. For the first time, GMEVI, charger has been
tested when operated in isolation from any other external disturbance or
interference. The testing procedure has been developed and applied to find the
current, power, voltage, power factor and harmonic characteristics of the EV
charger.

Computer simulation of the distribution system feeding EV loads has been
carried out. The simulation results have been compared with the actual field
and laboratory results. Thus, establishing the accuracy of the simulation
results.

Fast charging time EV chargers have been simulated and their expected
behavior has been studied. Although these chargers are not commercially
available, however, these simulation results have provided a basis for
understanding the expected behavior of such chargers once such a technology
becomes available.

Generalized load profile on a distribution transformer site has been developed.

An IEEE/ANSI method for the calculation of the derating of transformer
feeding non-sinusoidal loads has been modified. The proposed method yields
more accurate and realistic results.


1.4. Research Phases



The objectives of the thesis outlined in 1.3 have been achieved by dividing the

whole research work into various phases. These phases are outlined as follows and later


on in the section they are explained in more detail.








These phases are:


PHASE I DATA PHASE:


This phase consists of:

The recording and analysis of the field load data with and without the EV
being a part of the load.


The laboratory testing of the characteristics of the EV charger.




PHASE I. MODELLING PHASE:


This phase consists of:

The physical model of the various components of the distribution system.




PHASE 1 COMPUTERSIMULATIONPHASE:


This phase consists of:

Computer simulation of the power distribution system and the analysis of
simulation results.
Create and study of different load scenarios with multiple EV loads.
Create and study load scenarios with EV loads of variable charging times.



PHASE IV DESIGN PHASE:


This phase consists of:

Calculation of expected derating of a distribution system transformer
based on the simulation results for:
Normal charging time chargers.
Short charging time chargers. (Quick Chargers)













CHAPTER 2
LITERATURE SURVEY


The issue of power quality is as old as the issue of electric power itself But due to

changing load conditions and customer needs, the importance and significance of supplying

high quality power has changed with time. In the last decade, due to the increase in the use of

power electronic equipment, power quialr, has attracted the attention of the utilities and

customers alike. The utilities look to supply quality power to its customers at a price

dependent on the amount of distortion that the customer is causing. The customers on the

other hand want a reliable, sag free power to be delivered to them at a highly competitive

price.


The study and research in the field of electric power quality can be traced to 60's [1]

and 70's [2,3]. The fact though is that such studies were sparse. Power quality during that era

was just a limit applied to fluctuations of frequency and voltage, to the voltage unbalance,

voltage transients and flickers and power cuts [1,2]. Owing to the absence of vast amount of

non-linear loads, power quality and harmonics was a factor more of an academic importance

rather than a serious practical concern. The harmonics in the power systems were mostly

caused due to load switching and telephone interference. During the 80's with the first power

electronic revolution, the use of power electronic equipment was on an increase. Most of the

power electronic loads are switching loads and are highly non-linear; these loads thus cause a








lot of harmonics to be injected into the system. This prompted the academia to study the

power system harmonics in a greater detail, this included devising methods for the actual

measurement of system harmonics, making case studies and estimating their effect on

transformers, feeders and lines [4,5,6]. In the mid 80s [7] the increase in the presence of non-

linear power electronic loads was obvious and a further increase was anticipated. This increase

was owing to first the change in the technology e.g. faster and more efficient power electronic

devices, better converter topologies and converter control methodologies, second the advent of

fast microprocessors and their application to power electronic converter drives. The

microprocessor controlled drives were more efficient and reliable, had the flexibility of speed

control over a wide range, were easy to control and were compact in size. The need for

technical understanding of the nature of electric disturbances caused by these non-linear loads

was felt. Thus, in the face of changing load scenarios, there was a need for power quality

education that included definitions, grounding practices, power line dynamics, case experience

and trouble shooting [8]. In and after the late 80s a sharp increase in the published literature on

power quality and harmonics can be found. The number of these publications is so large that it

is impossible to list all of them, thus a selection of these papers has been listed here. Decrease

in power quality that not too long ago was a mere speculation, was now viewed as a serious

threat to the normal operation of power systems [10]. Research was thus directed towards

development of harmonic mitigation techniques [11] as well as development of power

electronic converters which would cause minimum harmonic distortion. Papers

recommending different methods for the estimation and monitoring of power system

harmonics were published [12,13,14]. References 15 25 include papers on the study and

characterization of specific non-linear loads and their effect on power systems. The studied








loads consist of fluorescent lighting loads, electric drive loads, railroad loads and some other

types of non-linear loads. In addition to these papers, scores of publications dealing with the

study of the effect of harmonics on various power system equipment can be found in the

literature [25 30].


The last few years have seen a phenomenal increase in the use of computers by the

power customers. The computer loads, which are highly non-linear [31], although, exist

significantly as commercial loads but the number of home computers is also increasing every

day. [32]. Power quality is no longer an issue that is linked with just reliability, efficiency,

voltage sags and flickers; today power quality is understood literally as the quality of the

actual current and voltage waveforms. The power customer today has become very sensitive

to the quality of power in terms of waveforms that the utility is providing. This sensitivity is

owing to the changing needs of the customers, for instance in some cases the equipment

powered from the line is highly sensitive to harmonics. Also, both commercial and domestic

customers are more aware of harmonics and the impact they have on the performance of

equipment and utility bill.


One of the very sources of harmonics in the power system is used to assess and

remedy the situation. Computer modeling and simulations are being used for the study of

various non-linear loads and the impact they have on power systems.[33-37]. In the midst of

the talk and concern about increasing disturbance and harmonics, yet another type of non-

linear load has appeared on the horizon. This is the electric vehicle load. The electric vehicle

load has attracted the attention of the researchers and the utilities alike.













CHAPTER 3
POWER QUALITY AND HARMONIC DISTORTION




Power quality, in simple words can be defined as the amount of distortion in the

system current and/or voltage waveforms. The power systems and equipment are designed for

pure sinusoidal current and voltage waveforms. Nowadays, widespread use of computers,

florescent lights, Adjustable Speed Drives and many other power electronic equipment by

customers cause a considerable amount of distortion of the voltage and current waveforms.

The distortion is caused due to the inherently non-linear characteristics of these types of loads.

In the presence of this increasing level of non-linear loads, the utility companies need to

continue catering to the customer needs, which is supplying quality power with a high degree

of reliability. This could necessitate new and stricter distortion guidelines for power

customers, as well as a variable rate structure being developed. Also, utility engineers need to

deal with the issue by analyzing the distortion and plan to control it by employing methods to

either minimize the distortion by using some mitigation techniques or derate the existing

system equipment to better cope with the steadily increasing non-linear load scenario.

Electrical energy consumption in the United States in the year 1992 was of the order of 600 *

109 kWh/year [38]. With the advances in power electronic converters, computers and many

other non-linear loads, more and more of the proportion of this power would be used by such

non-linear, distortion causing loads.








3.1. Harmonics




The non-linear characteristics of the loads cause the distortion of the current

and/or voltage waveforms in a power system. These current and voltage waveforms are

periodic in nature but they are not pure sinusoids. The distorted waveforms can be

expressed in the form of a Fourier series. Using the Fourier series expansion, a distorted

periodic wave can be expressed as the sum of an infinite number of sine waves, whose

frequencies are an integral multiple of the fundamental frequency.

The distorted periodic waveform f (t) can be expressed as:


f(t) = ao + ,[a,, cos(ncot) + b ni ...... i
Where:

ao is the dc offset of the waveform
a,, is the amplitude of the cos term
b,, is the amplitude of the sin term
n is the order of the harmonic


To study the effect of the distorted voltage and current waveforms on the power

system, the principle of superposition is used. This means that the effect of each

constituent (harmonic) of the distorted wave is studied separately and then the overall

effect is found by adding the effect of individual harmonics. The fundamental constituent

of the distorted waveform has a frequency of 60 Hz, the second harmonic has a frequency

of 120 Hz, the third 180 Hz, and so forth. The equipment used in the power system is

designed for a supply frequency of 60 Hz i.e. the fundamental frequency. The other

constituent waveforms (harmonics) certainly have an undesirable impact on the

equipment and the overall power system and its performance.










3.2. Estimation of Harmonic Distortion




Before assessing the qualitative and quantitative effect of harmonics on the power

system, some kind of estimate or quantitative expression of harmonic distortion needs to

be defined. Three methods of estimating harmonic load content are practiced [41]. The

Crest Factor and The Total Harmonic Distortion (THD) are the most common methods.

The third method is the K-Factor method.



3.2.1. Crest Factor





The Crest Factor gives a very simple estimate of the harmonic content of the

waveform. It is defined as the ratio of the peak of the wave to its Root Mean Square (RMS)

value. A perfect sine wave by definition will thus have a crest factor of 1.414.



Crest Factor = Peak Value of the Waveform
RMS Value of the Waveform


If the Crest Factor of the waveform (voltage or current) is other than 1.414, it

indicates some kind of distortion in the wave. This technique has a limitation, it does not

provide information about the constituent harmonic frequencies. Knowledge of the

harmonic frequencies and their respective amplitudes is important to study the effect of

wave distortion on the system.








3.2.2. Percentage Total Harmonic Distortion (THD)



Total Harmonic Distortion (THD) is defined as the ratio of the sum of the RMS

values of the harmonics to the RMS value of the fundamental.




THD = 2


The above expression shows that THD provides with a measure of the distortion

of a waveform. Like the Crest Factor, the THD is also limited in that it does not provide

any information about the frequencies and amplitudes of individual harmonics. While

computing the THD for a wave, the relative phases of the various harmonics are also not

considered. Nevertheless both Crest Factor and the THD are very useful to establish the

presence or absence of waveform distortion.



3.2.3. The K-Factor





The K-Factor while providing an estimate of distortion in a waveform takes into

account the harmonic frequencies and amplitudes. This fact makes it the most accurate

available method of estimation of non-linear load harmonics in a system. This is the method

used for the calculation of derating of dry-type power distribution transformers [42].

Reference 4 proposes an expansion of the K-Factor method. The author [41] proposes to use

the K-Factor as a weighted sum. Each type of load is considered separately and then all the K-

Factors are combined to create a composite harmonic fingerprint.








The K-Factor is defined as:


K = I/,(pu)2h'
h=1



3.3. Transformer Delating




Power transformers and distribution transformers are designed for sinusoidal

currents and voltages. Name plate ratings of transformers hold for purely sinusoidal

supply and loads. If the current gets non-sinusoidal, for instance because of the use of

switching power supplies, solid state inverters, converters, Adjustable Speed Drives etc.,

the indicated ratings of transformers no longer hold. A high harmonic content in the

output current waveform effects the transformer ratings. Load harmonic currents add a

disproportionate amount of additional transformer heating, particularly in the form of

winding eddy current losses and stray losses which are supply frequency dependent.

Several methods are currently used in the industry to assess the derating of transformers

for a harmonic application. The most widely used method is IEEE/ANSI C57.110. The

IEEE/ANSI method was developed for power transformers and does not specifically

address other applications such as distribution and dry type transformers. On the other

hand, the K factor method can be applied to only dry type transformers. In general, the

application of the IEEE/ANSI method provides conservative results when applied to a

distribution transformer. This standard outlines the following method for calculating the

transformer harmonic capability.








3.3.1. Transformer Capability Equivalent Calculation.




In this method a per unit value of non-sinusoidal current is calculated which will

result in load losses in the transformer just equal to the highest loss region for 60 Hz rated

operation. The winding eddy current loss under rated conditions at the point of maximum

loss density is 15% of the local 12R losses. Taking the fundamental component of current

as 1 per unit and taking the measured per unit values of the current harmonics up to the

SI 1h harmonic the derating of a transformer is calculated.

The maximum loss density of the local 12R losses is assumed as 15%. Solving

Equation 8 in ANSI/IEEE C 57.110 for the permissible current.

Following the procedure in the standard:


1.15

I + f x 0.15



Using the values shown in Table. I. and obtaining the amplitude of each harmonic

from Figure 4, the maximum permissible non-sinusoidal load current is given by:




113.089 0.501
Hu.15
4.75

Hence the ratio of transformer rating with and without the harmonic loads = .501
















Table. I Harmonic Content of the Sample Waveform


HARMONIC PER UNIT CURRENT (FH)

(H) FH2 FH2H2

1 1 1 1

2 .16 .0256 .102

3 1.3 1.69 15.21

4 .267 .0712 1.139

5 1 1 25

6 .101 .010 0.36

7 .68 .462 22.638

9 .54 .291 23.57

11 .447 .199 24.07



Total 4.75 113.089








3.3.2. Effect of Current Harmonics on the Neutral Current.





The previous sections of this chapter discuss the effect of harmonics on the

distribution transformer. However there is some more concern for harmonic loads on multi-

grounded wye systems [43,44]. The current harmonics in a Y-Y distribution transformer

cause a high neutral current and an elevated emf of the neutral.


In an ideal case when the loads on the transformer are balanced, the neutral current

which is the vector sum of the three phase currents is zero.


If:

I'A = Asincot
1B = Asin(t +2400)
Ic = Asin(cot+1200)

where A is the amplitude 'A 1, and Ic are the currents in phase A, B and C respectively


The neutral current is the vector sum of the three phase currents and is given by:


IN =IA +IB+IC

'N = Asin at + A sin(w t + 2400)+ A sin(w t + 1200)
.'. IN =0

Now if the load current is rich in the third harmonic, the third harmonic component


with amplitude T can be written as:

I'A = T sin(3a)
1B =T sin(3 t+ 2400)
Ic = Tsin(3ot +1200)








And as before


'N =IA+IB+IC

NI = T sin(3o t) + T sin(3co t + 2400)+ T sin(3o t +1200)
.'. I, =3Tsin(3wt)


The net neutral current because of all the other tripplen harmonics will yield similar

results. As reported by Mansoor et al. [44], in a power system, even with load diversity, there

is very little cancellation of the third harmonic. These high neutral currents cause the

following problem [45]:


1. The heavy neutral currents cause the over loading of feeder transformers. They
can cause further distortion of the voltage because they push the transformer into
saturation. The saturated transformer then becomes a source of harmonics.

2. When the neutral conductor carries harmonic currents, additional heat is generated
and the current carrying capacity of the feeder is reduced.

3. The neutral conductor for transformers is generally sized same as the phase
conductor, but due to the tripplen harmonic currents the neutral conductor current
can exceed its rated value thus overloading the neutral that can lead to equipment
failure.

4. The balanced tripplen currents flow as circulating currents in the transformer delta
primary winding. As a result, more current flows in the winding of the transformer
than is detected by the transformer primary circuit over current protection device.
This can result in the overloading of the transformer.

These high neutral harmonic currents cause a voltage differential between the neutral

and the ground since the ground wire at harmonic frequencies can cause significant neutral

conductor voltage drop. This is known as Common Mode Noise Voltage.













CHAPTER 4
FIELD DATA RECORDING, REPRESENTATION AND OBSERVATION.


4.1. Introduction.



The effect of the EV loads on the power distribution system would depend on a

number of factors. Once the EVs become popular one of the factors that is going to be of

significance, is the usage pattern of the EV by its users. The usage pattern would

obviously vary from individual to individual, it would depend on the lifestyles and

specific needs of the customer, e.g. nature of his job, weather conditions, age of the user,

so on and so forth. To generalize the usage pattern and come up with an accurate

generalized model is surely a nontrivial problem. Based on field-recorded data and

knowledge of the EV load behavior, a typical load profile of a distribution node feeding

EV loads can be constructed. Even though this load model won't be accurate but it will

contain enough information to enable the utilities to design the distribution transformers,

so that they are able to supply power reliably and economically. The usage pattern of

users on the same distribution transformer shall determine the load diversity factor, the

expected peak load on the transformer and the overall derating of the transformer.

The objective of the field data monitoring thus is threefold. First to get the load

profile during two separate periods viz. when none of the customers are using an EV

and when one of the customer is using an EV. Second to try and identify the signature or

effect of the EV load when it is connected to the system. This effect is in terms of








increase in power demand and possible decrease in power quality. Lastly, this field

recorded data can form the basis for configuring system computer model used for

carrying out computer simulation study.

The objectives of field data monitoring have been achieved by monitoring the

field data (system voltages, currents, power factor etc.) first for a period of time (two

weeks) when the EV load is not a part of the loads on the distribution transformer. This is

followed by a two-week period then when one of the customers is using the EV. The data

obtained during the above two load scenarios has been analyzed. The comparison of these

load scenarios has been used to establish a generalized load scenario when one of the

customer on a distribution node is using an EV.



4.2. General Considerations In Data Acquisition And Processing



Appropriate techniques for the acquisition and processing of random data are

heavily dependent upon the physical phenomenon represented by the data and the desired

engineering goals of the processing. In broad terms, however, the required operations

may be divided into five primary categories as follows:

a. Data collection

b. Data recording (including transmission)

c. Data preparation

d. Data qualification

e. Data analysis

Each of these categories involves a number of sequential steps. The purpose of

this section is to summarize basic considerations associated with each of these key steps.








4.2.1. Data Collection



The primary element in data collection is the instrumentation transducer. In general

terms, a transducer is any device which translates power from one form to another. In an

engineering context, this usually means the translation of some measure of a physical

phenomenon of interest into an analog signal with a calibrated relationship between the

input and output quantities. This translation may involve up to three basic operations; (a)

a mechanical conversion of the physical quantity of interest into an intermediate

mechanical quantity, (b) a pick off step which converts the intermediate mechanical

quantity into an intermediate electrical quantity, and (c) an electrical conversion into a

final electrical quantity, usually voltage. Some transducers may combine any two or all

three of the operations, depending upon the physical quantity being measured and the

specific nature of the transducer.

For example, a thermocouple, which is widely used as a temperature transducer,

converts a difference in temperature directly to a difference in voltage without

intermediate steps. On the other hand, a resistance thermometer, another commonly used

temperature transducer, first converts temperature to a change in resistance, and then

converts the resistance change to an electrical voltage change.

Ideally, the above operations would be accomplished without distortion or

modification of the time history of the physical quantity being measured. In other words,

if the input time history is x(t) and the output time history is y(t), the perfect transducer

would provide an analog output, y(t) = c x(t), where c is a simple calibration constant.

Unfortunately, this ideal situation is difficult to achieve in practice. Gain and phase

modifications as well as distortion producing non-linearities are often inherent in the








transducer operations. This fact makes the transducer a potential source of error in any

data acquisition and processing program.



4.2.2. Data Recording



For some applications, it is possible to perform all desired data processing directly

on the transducer signals in real time. For most applications, however, this is not

practical and some form of storage (and perhaps remote transmission) of the transducer

signals will be required. The most desirable and convenient type of data storage system

is the magnetic recorder or a computer hard drive memory. Although other types of

recorders could be used, the magnetic recorder has the advantage of being able to store

large quantities of data, and to reproduce them in electrical form. The most desirable

way to transmit data signals is through electrical lines or telephone lines. There are

obvious situations where this is not feasible; for example, retrieving data from a

spacecraft in earth orbit. For such cases, radio transmission (telemetry) of the transducer

signals is usually required.



4.2.3. Data Preparation

The next key phase in data acquisition and processing is the preparation of raw data

for detailed analysis. The raw data are usually supplied from the recorder in the form of

voltage time histories (or as direct analog signals from the transducers, if the data are being

analyzed on-line). A number of operations are needed at this point to make the voltage time

histories suitable for detailed analysis. The first of these operations is generally classified as

data editing.








Data editing refers to those pre-analysis procedures which are designed to detect

and eliminate spurious and/or degraded data signals which might have resulted from

acquisition and recording problems such as excessive noise, signal dropout, loss of signal

due to transducer malfunctions. Editing can often be accomplished through visual

inspection of the data time history signals by a talented analyst. In more elaborate data

acquisition and processing systems, a specific instrument for quick-look evaluation might

be employed. Real-time spectrum analyzers are popular for this application. It should be

noted that the data editing step is more critical for the case of digital processing than for

analog processing. This is true because once the data have been converted to a digital

format, it is often difficult to detect even the most obvious errors in the original signal.

For the case of analog data processing, data preparation beyond editing usually

includes only conversion into engineering units (calibration).



4.2.4. Data Qualification



The correct procedures for analyzing random data, as well as interpreting the

analyzed results, are strongly influenced by certain basic characteristics which may or

may not be exhibited by the data. The three most important of these basic characteristics

are the stationarity of the data, the presence of periodicities in the data, and the normality

of the data. Stationarity is of concern because the analysis procedures required for non-

stationary data are generally more complicated than those which are appropriate for

stationary data. Periodicities in the data should at least be identified to avoid erroneous

interpretations of later results. The validity of an assumption that the data (excluding

periodicities) have a Gaussian probability density function should be investigated since








the normality assumption is vital to many analytical applications for random data.

Qualification of sampled data in terms of these basic characteristics is indicated as a

separate operation to be performed prior to detailed data analysis. In practice, however, it

is often accomplished as an integral part of the data analysis phase. Practical

considerations and procedures for such qualification will now be discussed.



4.2. Data Analysis



The procedures for analyzing the properties of random data may be divided logically

into two categories: the procedure for analyzing individual sample records, and the procedure

for analyzing a collection of sample records given the properties of the individual records.

Applicable data analysis procedures for these two categories are now outlined.


MEANAND MEAN SQUARE VALUE ANALYSIS.


The first step is a mean and mean square value (or variance) measurement. This step

is almost universally performed for one or more of three sound reasons. First, since the mean

and mean square values are the basic measures of central tendency and dispersion, their

calculation is generally required for even the most rudimentary applications. Second, the

calculation of short time averaged mean and mean square value estimates provides a basis for

evaluating the stationarity of the data. Third, mean and mean square value estimates can be

extracted from other descriptive properties (probability density plots, correlograms, and/or

power spectra) which might be measured later. The comparison of directly measured mean

and mean square values estimates to the corresponding estimates extracted from other








analyses provides an excellent method for checking the data analysis equipment or computer

programs for correct operation.


POWER SPECTRAL DENSITYANALYSIS.

Perhaps the most important single descriptive characteristic of stationary random

data is the power spectral density function, which defines the frequency composition of

the data. For constant parameter linear physical systems, the output power spectrum is

equal to the input power spectrum multiplied by the square of the gain factor of the

system. Thus power spectra measurements can yield information concerning the

dynamic characteristics of the system. To be more general, the mean square value of the

data in any frequency range of concern is determined by the area under the power

spectrum bounded by the limits of that frequency range. Obviously, the measurement of

power spectra data, will be valuable for many analysis objectives. Secondary

applications include its use for the detection of periodicities and as an intermediate step in

the calculation of autocorrelation functions.

PROBABILITY DENSITYANALYSIS.

The last fundamental analysis included in the procedure is probability density

analysis. Probability density analysis is often omitted from a data analysis procedure

because of the tendency to assume that all random phenomena are normally distributed.

In some cases, however, random data may deviate substantially from the Gaussian form.

If such deviations are detected by a test for normality, then the probability density

function of the data must be measured to establish the actual probabilistic characteristics

of the data. Furthermore, a probability density function estimate is sometimes used as a

basis for a normality test.








4.3. Field Data Monitoring

4.3.1. Data Acquisition Methodology.




The data was recorded in the form of one cycle snap shots of voltage, current and

power factor. These waveforms were acquired once every one minute interval. The

voltage and current waveforms at the distribution transformer output terminals were

recorded and stored temporarily in the local hard drive of the data acquisition system.

The stored data was then down loaded to a personal computer over a telephone line using

modems. The distribution nodes under study were feeding typical domestic loads at 110

volts and EV loads at 220 volts. At each of these locations the distribution transformer

has typical domestic loads of 2 to 10 homes. At each location one EV was provided to

one of the customers for a period of two weeks. The customer was asked to put the EV to

normal use. The data at each location was recorded during time intervals that can be

classified into three states. The first state is the period (about two weeks) when there was

no EV load. The second state is the period (about 2 weeks) when one of the houses on the

distribution transformer was using an EV. The third state was again when there was no

EV load. Thus the field data collected represents the load conditions before; during and

after the EV load was connected to the system. Figure 4.1 shows the data acquisition

system layout.

4.3.2. Data Monitoring Instrument.


The field data was recorded using several BMI models 8020 Plus PQNodeTM data

acquisition instruments. The BMI model 8020 Plus PQNodeTM monitors power quality

phenomena in electric power systems. It combines the capabilities of instruments that







monitor continuous quantities, (such as voltage level, load level, power factor, and

harmonic distortion), with the capabilities of instruments that record disturbances

(impulses, wave shape faults, swells and sags, outages, cold load pickups). Monitoring

continuous quantities requires that the monitored signals be sampled at a regular basis,

while monitoring disturbances requires that the signals be continuously sampled, and

recorded only if the signals exceed specified values. The BMI 8020 Plus PQNodeTM

continuously samples up to four voltages and four current three phases and neutral. The

instrument periodically records waveforms, RMS levels, frequency and temperature so

those trends can be calculated. The PQNodeTM is designed to interface with a personal

computer via telephone line, and is equipped with an internal modem for this purpose.



4.3.3. Data Management



The data acquired by the BMI data acquisition system was stored by the BMI unit

in the hard drive of the device in the form of ASCII files. Each set of data was stored as a

4 K file. The actual information is in only 1 kilobyte of the 4 K space. After downloading

the data files to the personal computer in the Power Quality and Power Electronics Lab,

at The University of Florida, the data needed to be stored safely with adequate backups

and at the same time consuming optimum space. The data also needs to be easily

restorable whenever needed for analysis purposes. Using data acquisition and data storing

software, the data was archived in the computer hard drive. The archived data occupies

only one fourth of the space since the files are crunched together. The data was also

backed up on tape drives connected to the computer. In this way it was made certain that

the precious data is not lost in the event of a hard drive failure. For the purpose of








analysis, only the required amount of data is restored and after it is processed it is deleted

from the hard drive and a fresh batch of data to be analyzed and processed is restored

from the archives.

For each batch of data the THD trend is plotted. Also plotted are the voltage and

current trends for all the harmonics from the second through to the fiftieth harmonic

component. In addition to the harmonic trends, steady state trends are also plotted; these

are the real power, apparent power, RMS current and power factor trends for the two

phases between which the EV charger is connected. Each of the trend waveforms is

displayed one by one.

After examining the harmonic and steady state trends for the selected days,

typical real time steady state voltage and current waveforms at certain instances of time

were also plotted. It may be mentioned here that the BMI instrument acquires the

distribution node data in the form of real time waveforms of voltages and currents. Using

these waveforms which are acquired during every fixed period of time, the computer

software uses them to create various harmonic and steady state trends.












































Figure 4.1. Data Acquisition System Layout













CHAPTER 5
DATA REPRESENTATION AND ANALYSIS


5.1. Introduction



This section deals with the graphical representation of the recorded data. A visual

inspection of the data trends is enough to show that the data represents stationary random

process. It is assumed that all random phenomena are normally distributed. [data book] In

some cases, however, random data may deviate substantially from Gausian form. Looking at

the data trends included in this chapter, the properties of the recorded data can be described at

any instant of time by computing the average values over a collection of sample functions that

define the random process. The stationary random data represented here falls in the category

of ergodic random data. The properties of ergodic random process can be determined by

performing time averages over a single sample function. Fortunately, in practice, random data

representing stationary physical phenomena are generally ergodic. It is for this reason that the

properties of stationary random phenomena can be measured properly, in most cases, from a

single observed time history record.

The recorded data is presented in the form of various trends and waveforms. Each data

trend for a particular site with and without the EV load are presented together to facilitate

comparison between the two load scenarios namely when the EV load is the part of the load

and when the EV is not a part of the load. The data recorded using the BMI data acquisition








units was synthesized using the software PASSTM. Broadly speaking two types of data trends

were created and displayed graphically, viz. the steady state trends and the harmonic trends.

For the purpose of studying the impact of Electric Vehicle loads on the distribution system,

only the following trends have been selected:

1. The Steady State real power trend.

2. The Load Current Total Harmonic Distortion (THD) trend.

3. The various harmonic trends of the load current.

Although there are many other steady state and harmonic trends that were

displayed using the PASSTM software, only the above mentioned three trends have been

used for data analysis for the following reasons:

The EV charger load is a real power load, which consumes a peak power of 7

kilowatts (based on GMI charger). The 7-kilowatt power is consumed when a fully

discharged EVI battery is connected to the charger. The power decreases as the charge .

on the battery increases. So one of the impressions that the EV charger load causes in the

system is an increase in real power. Thus the real power trends for the distribution node

with and without the EV load were created. This made it possible to compare the load

scenario with and with out the EV load. From the laboratory testing of the EV charger, it

was found that the EV load causes a distortion of the input line current. The General

Motors EVI is known to cause very little input voltage distortion [2].

The data represented in this chapter can be broadly classified into two categories

viz.:

I. When the distribution transformer is feeding normal domestic loads.

II. When the distribution transformer is feeding normal domestic loads that
also include Electric Vehicle Load.








Representing the data in this manner has facilitated the comparison of the load

scenarios with and without the EV load.



5.2. Data Trends and Characteristics




Figure 5.1 through 5.7 are sample data trends created from the actual monitored field

data. Each of these figures consist of two parts, viz. part a and part b. Parts a are the data trends

without the EV loads, parts b are the data trends when the EV load is a part of the loads on the

transformers. Following is the summary of highlights of each of the figures 5.1 to 5.7.


Figure 5.1: This figure represents the Total Harmonic Distortion trends of the line

voltage. The voltage THDs for both cases are very close. Maximum voltage THD variation

between "a" and "b" is 2%. The presence of EV loads does not have a significant effect on the

input voltage THD. This fact is verified by the laboratory testing of the EV charger.


Figure 5.2: Unlike the voltage, the current THD trends shown in figure 5.2, are

sensitive to EV loads. Comparing 5.2 a and 5.2 b, it can be seen that the current THD trend for

two cases is not the same. When EV was a part of the domestic load, we see in 5.2 b, the

current THD between the hours of 2:15 and 4:40 is very high. The current THD during this

time varies between 80% to 40%. For the remaining part of the day, the current THD for both

cases is close to each other.


Figure 5.3: These figures show the real power trends. The real power trend, when EV

is a part of the load is not much different from the case, when EV is not a part of the load. This

could be due to the fact that the customers do not fully discharge the batteries before








recharging them. Thus each time the battery charger load is put on the system, it draws only a

fraction of its maximum rated power of 7 kilowatts.


Figure 5.4: These show the power factor trends. The EV load does not effect the

power factor. Thus we see, in both parts of this figure the power factor trend stays close to .8.


Figure 5.5: These figures represent the third harmonic current trends. As in the case of

the overall current THD trends, the third harmonic current trend between the hours of 2:15 and

4:00, show a high, for the case when EV is part of the load. The third harmonic component

touches a high of 1.5 amperes during this period. This again indicates a possible presence of

EV loads during these hours.


Figure 5.6 and 5.7: These are the fifth and the seventh harmonic current trends. The

conclusions drawn from these trends are similar to those drawn from figure 5.5.





5.3. Data Observation




Observing the data trends represented in figures 5.1 through 5.7, the following

conclusions can be drawn:


I. The EV load does not have a distorting effect on the input line voltage. Thus
the input line voltage THD is not effected by the presence of EV loads.

II. The EV loads cause a considerable distortion of the line current. The line
current THD is particularly high when the amplitude of the fundamental of the
line current is small. Thus the current distortion is maximum close to the end
of the charging cycle, when the charger draws a small current.








III. The EV users generally would not be expected to fully discharge the batteries
before recharging them. Thus the EV loads may not have a significant effect
on the expected peak load on the distribution transformers.

IV. The EV battery charger operates at a high power factor. Thus the reactive
power drawn by the charger is low.

V. The EV charger injects a significant amount of lower odd current harmonics
into the power line.




5.4. Load Characteristics




This section contains the computed maximum and the average daily power output of

the distribution transformer. The daily average power is the average of the load over a period

of two weeks. The load average has been found, first for the period weeksk) when EV was

not a part of the loads and next for the period of 2 weeks, when EV was a part of the load.

Also, for these two two-week periods the maximum load profile has been computed. The

figure X.8 a, has been derived by computing the load average for a 2-week period when EV

was not a part of the loads. The figure 5.9b is the 2-week load average when EV was a part of

the loads. Comparing, 5.8a and 5.8b, we see that there is an increase in the average power

consumed, when the loads contain the EV load. This increase is particularly visible, during the

early morning hours, this reinforces the previous conclusion about the times at which the EV

load was connected to the system.


Figure 5.9 a and b, show the maximum power recorded during the two two-week

periods. Plot in figure 5.9a is the maxima of the daily load curves of the two weeks period

when EV was not a part of the loads. Plot 5.9b is the maxima of the daily load curves of the








two-week period when EV was a part of the loads. It must be noted, the EV loads although do

not effect the peak loads, nevertheless being an additional load it does affect the average daily

power consumed.


Figures 5.10 a and 5.10 b, show the snap shots of the recorded line current waveforms.

The line current snap shot of figure 5.10a was recorded when EV was not expected to be a

part of the loads and 5.1 la was recorded, when EV was expected to be a part of the loads.

From the Fast Fourier Transforms (FFTs) of waveforms in figure 5.10a, shown in figure

5.1 Ob, we see that the line current does not have a high content of the lower odd harmonics

(viz. 3rd, 5th and the 7"). This is not the case with figure 5.1 lb, where these harmonics are

significant.


Based on the preceding discussion, the load profile of Figure 5.9, qualifies for a

generalized load profile of loads on a distribution transformer, for that particular site.









Voltage THD Trend: No EV Loads


Time of Day



(a)







Voltage THD with EV Loads


0 0 0 (N (0 S mo 0 V 2

Time of Day




(b)


Figure 5.1 Voltage Total Harmonic Distortion Trends

















Current THD Trend: No EV Load


o 20

I 15

25


U. I


0


'I U, vi 10 CO



Time of Day




(a)




Current THD Trend: EV Load


100


80 -
60

(b)20
O. Z 0




Time of Day



(b)


Figure 5.2. Line Current Total Harmonic Distortion Trends








Real Power Trend: No EV Load



D 4
3 3
co 2
0

77o 7 -4 N N o o in 9 o Sn n p

Time of Day



(a)





Real Power Trend: EV Load

10

6 ---

I g
0 0
t 2


'N Ip .O. ., ; z b d,

Time of Day





(b)


Figure 5.3. Real Power Trends














Power Factor Trend: No EV


12


0.6
0.2
a0 0

0-~ s _b 0 l .- -. <. ,:. ---<- IV N.- N' 4-N
&u ^ n' 41> 1' 41. 41 <1 41 41 ,* 4'-1 441 41. (? >^ n i

Time of Day



(a)





Power Factor Trend: EV Load


*5 1
0.8
0.6
0.4
_026 -------------- -- ---------- I
0.2
M 04


1 aS az o~o -i**''oOO.. 41 .; 41 41 41 \0 i 41

Time of Day




(b)


Figure 5.4. Power Factor Trends















Third Harmonic Current Trend: No EV Load


E 5
< 4


-- 3? N 'N- N N N
2- j--- ---- -




C ) C- N V V N '


-Ji



N C CN T t C C CO N rN 0


Time of Day


(a)







Third Harmonic Current Trend: EV Load


4
3 -----------------

2







Time of Day


-n +,- A

'2 2 2


Figure 5.5. Third Harmonic Current Trends


BC:
8=













Fifth Harmonic Current Trend :No EV Load

0- 25
E 2
C 15

1

0
25 --------------------------.--.-----.------.-------,--





Time of Day




(a)








Fifth Harmonic Current Trend :EV Load

m 2
< 1.5

0 0





Time of Day




(b)


Figure 5.6. Fifth Harmonic Current Trends













Seventh Harmonic Trend : No EV Load


S 12
E- 1 2 ----------------------------
E 08







Time of Day



(a)







Seventh Harmonic Current Trend : EV Load

12

0
S 0.6 -- --
q 04
` <, 8; 1,( I' 0*^.^*.^**{*' *^ *^ ^y^^^^ '.


























Time of Day



(b)
.E 0 .6 ,-------- -----










(b)


Figure 5.7. Seventh Harmonic Current Trend












Average Power: No EV Load



6 1 1



C 6 I... .i .$ A Ji._ ..i-M ., ...A ,---
(a 0

'S. ^ y *y O{* (6 6 N6 C8. 3. 1 *w

Time of Day




(a)





Average Power: EV Load




8 &:
o> 1 2r---------------------------------

. 0 .
a-o


11P '?1* o"1 o11P 11 N 11% 14" NIP le 11 6N ^'10
Time of Day




(b)

Figure 5.8. Average Power Computed on the Basis of Two Week Daily Load


















Max Power Demand: No EV Loads


12 ---- --
10
8
6-


4 -PP 0'-- P *s -.P P- -- -
2 -------------------------------------



Time of Day



(a)






Max Power Demand : EV Load



25
20 -
c 15

10 -
0

0

4' 4' "''0 S,> ,0 0 "* e 0 4 ,0 e* ",?-l .

Time of Day





(b)

Figure 5.9. Maximum Power Demand Computed over Two Week Period














radwClirventrwawemi wi NEV Lo[ad


5 10 15
noe ( 5

(a)






FFrOftie wad curritWaeirm wj WNo EV ad



TD :E
S: 3.310773
RMSi : 0


06 --- --- ----- RMS 32.6343
AsJM : 395032
TF 130.916
04- ---r :427332





00 1 1- 1 I I II 1. I I ll II l i. II ... II I I I
120 720 1320 1920 2520 3120



(b)

Figure 5.10. a) Line Current Waveform Snap Shot with No EV Load
b) Fast Fourier Transform of the Current Waveform of "a"


4.




lo~---- -
-10 _


I
L -<

v


L LJ.IT

CF : 1.61574
-FFr: 112339


v
|C'


U













Fe i Recorded b Jrnentwil EV Loads


Hr, I.,'1.-
Ag,: 533581
Abs: 10.7422
RMS: 6/3585
CF : 175933
FF : 114432


10
Treiusk


(a)


Hrmon m ii Coentofte Fe-iRecrGedLad Current
r- -: 6
Rrd 153EI
T-H :12278
TD :0
TD :12S78
RMSh 197152
-- RMSi: 0
RMSft: 15A95
RMS :15A95

T : 356.55


--t -.. -. _- ___







720 1320 1920 2520 3120




(b)


Figure 5.11.
a.) Field Recorded Current Waveform.
b.) FFT of the Field Recorded Waveform in Figure a.


ISO


125


| 100

S075-


050-


025


000
120


i : iL, ru tff-l-c


-













CHAPTER 6
LABORATORY TESTING OF THE EV CHARGER



6.1. Testing Procedure and Results



The EVI charger is a 240 volt split single-phase device. In the laboratory

testing of the charger, a nearly fully discharged (5% of full charge) battery was connected

to the charger. The voltages and currents on both input lines were measured; the

equipment set up for the laboratory testing is shown in figure 6.1. Electrical

measurements were taken using the power monitor, automatically once every one-minute

interval for the entire charging cycle. The fully discharged battery was connected to the

charger at 10:17 AM. The battery reached its almost full charge at 2 PM, a charging time

of a little less than 3 hours. The EV1 charger when charging a fully discharged battery

draws a current of 30 amps at 240 volts and is also a source of significant harmonics as

shown in the following laboratory results. These results have been presented in the form

of waveforms and graphs that depict the electrical characteristics of the EVI charger.

The first graph, Figure 6.2 shows the power consumed by the EVI charger over

the entire charging cycle. The load varied from the maximum of 7kw at the beginning of

the cycle to about Ikw when the battery is nearly completely charged. During the entire

charging cycle the input voltage remained unaffected as far as either the magnitude or the

harmonic distortion is concerned (Figure 6.3).








The input current however varied significantly in percentage harmonic content.

Figures 6.4, 6.5, 6.6 and 6.7 show the current waveforms recorded at specified instants of

time. It can be observed that as the amplitude of the current decreases the waveform distortion

becomes more and more visible. Figure 6.9 shows the RMS values and figure 5.10 the Total

Harmonic Distortion (THD) of the input currents for phases a and b, over the entire charging

cycle. The figure 6.11 shows the harmonic content of the input current. It was seen that the

amplitude of various harmonics throughout the charging cycle does not vary significantly.

Since the amplitude of the fundamental component decreases with time and the THD, which

is a percent measure of harmonics to the fundamental shows a significant increase with time.



6.2. Observation of Laboratory Results



The EV charger charging into a fully discharged battery draws a line current of 30

Amps at 240 Volts. This corresponds to a power consumption of approximately 7 kW at time

10:17 as shown in Figure 6.2. The voltage waveform as shown in Figure 6.3., stays distortion

free throughout the charging cycle. The input line current distortion at the beginning of the

charging cycle is less visible (Figure 6.4). This is because the amplitude of the line current is

high and its THD which is expressed as a percent of the fundamental RMS current, is low.

Note that in Figure 6.10 the THD at time 10:18 is just about 2 % for both line currents laand

Ib. As the charge on the battery increases and the amplitude of the line current decreases, the

current distortion becomes more and more visible (Figure 6.5 to 6.8). This is verified by

Figure 6.10 where in at the end of the charging cycle, the THD of the line current climbs to

almost 30 %. Although the THD of the line current increases as the charge on the battery

increases, the actual harmonic content of the line current stays fairly constant. Figure 6.11





49


shows the line current harmonic content up to the 50h harmnic. The most significant

harmonics are the 3rd, 5h and the 7t.The even harmonics are almost not present.













Line 1
Line 2
Neutral I


Currents


Power Monitor
Digital Oscilloscope
Dynamic
Signal Analyzer
Power Harmonics
Analyzer


EV
Charger

240VAC
Single Phase


Voltages


Figure 6.1. Experimental Setup for EV Charger Testing


N
O











Charger Power Output

8000
m 6000
1 4000
2000



-1 i-1 r-i I--I --i 1--I -I --


Figure 6.2. Power Curve for the EVI Charger.



















200

150

-too i- / : ; : /


50




S- \50

U3 -100-




-200
0 5 10 15 20 25
Time (ms)




Figure 6.3. Typical Charger Output Voltage Waveform
















60



40


20 \






-20-


-40 ,



-60
0 5 10
Time (ms)


Figure 6.4. Charger Input Current Waveform Sampled at 11:19


\;
\i


-i


\


\o














40 -


30


20


10 -


0

S-10 -
-io

-20




-40
0 5 10 15 20 25
Time (ms)




Figure 6.5. Charger Input Current Waveform Sampled at 11:40
















I 1



T


I T










10
Ti me


(mS)


Figure 6.6. Charger Input Current Waveform Sampled at 12:07


Tj me


\ /






5 20 25






































Figure 6.7. Charger Input Current Waveform Sampled at 12:36












1 .0 1 1 111


0.5 --




0.0




-0.5




-1.0
0


<4


10 15
Time (ms)


Figure 6.8. Charger Input Current Waveform Sampled at 13:53


20


i


\I,



















Current RMS

35


30 -- 11 RMS
--i2RMS
25 --


20

15 -

10 -






10:18 10:42 11:02 11:22 11:42 12:02 12:22 12:42 13:02 13:22 13:42





Figure 6.9. The Charger Current RMS Variation over the Charging
Cycle



















Current THD


35

30

2la THD
25 lb THD

10 20




15

0



10:18 10:42 11:02 11:22 11:42 12:02 12:22 12:42 13:02 13:22 13:42



Figure 6.10. The Charger Current THD Variation Over the Charging Cycle.



















Ia Max Harmonics


0 0.6


0.4


02



0\ a 'r aO M 0 CM,1 'D 00 0 N *' D 0 0 0 CM x, c) 0
o 'N M M M M ( 't 't I
Harmonic



Figure 6.11. The Charger Input Current FFT.













CHAPTER 7
COMPUTER SIMULATION


This section contains the results of computer simulations of the power distribution

system with domestic and Electric Vehicle (EV) loads. The domestic load included in

these simulations is arbitrary. The EV model used is based on the laboratory data

obtained by testing the actual charging characteristics of an EV charger used by the GMC

EVI. The aim of this simulation is to see the effect of EVs connected to the distribution

system. The transformer current with EVs present in the system is simulated. The FFT of

the transformer current with and without EVs are shown and different harmonic

components are found to study the effect of these harmonics on the distribution

transformer in terms of the derating of transformers.

The power distribution system with an EV charger load was simulated using

SABERTM. This section includes the simulation results and discussion about the

simulation model.



7.1. Component Models


7.1.1. System Model


The system model consists of the distribution transformer, domestic load and the

EV charger and battery load. Figure 7.1 shows the layout of the simulated power

distribution system. The distribution transformer has typical domestic loads connected at








































Figure 7.1. Distribution System Model.











































Figure 7.2. Simplified Model of the Simulated EV Charger








120 volts between the center tap and the outer terminals. The outer terminals, which are

at 240 volts, are connected to a full bridge rectifier. The bridge rectifier is used to model

the inductive charger used by the EV to charge its batteries.



7.1.2. Electric Vehicle Charger And Battery Models




The EV charger is modeled on the basis of laboratory data obtained during the

charging of a fully discharged EV battery connected to charger. The full bridge rectifier

represents the charger. The input to the rectifier is 240 V from the distribution

transformer. The EV charger and battery model is shown in figure 7.2. The battery has

been modeled as a series combination of a resistor and a capacitor. The resistance

represents the loss component of the battery. It is presumed that the charging pattern of

the capacitor shall closely represent the charging pattern of the EV battery. The values of

the resistor and the capacitor can thus be chosen to achieve the desired charging current

and charging time. In the simulation while the full wave bridge rectifier represents the

EV load, a composite current source has been used to represent the harmonic effect of the

EV charger. The composite current source is injecting harmonics of the same amplitudes

and respective orders as shown in figure 3.2. The actual waveform of the harmonic

current injected by the composite current source is shown in Figure 7.3.



7.1.3. Transformer Model

The transformer used for the simulation is a gap less three winding transformer

with one primary and two secondary windings. The transformer in the simulation has the

following specifications.



























3 i m ------ -. _________ --- _








-3 Hroict ha re






Figure 7.3. The Harmonic Current Injected by the Charger








These specifications were provided by FPL.


7.1.4. Domestic Load Model




The domestic load has been modeled as a combination of some lighting and heating

load represented by resistive load. The load has been modeled such as to exhibit the

characteristics similar to the actual load characteristics as determined from the field data.

Figure 7.4 shows a typical domestic load model used in simulation.



7.2. Mathematical Model Of The EV Charger.


The EV charger can be represented by a single-phase bridge rectifier with an R-C

load, as shown in figure 7.1 can represent the EV charger. A current source model is used for

the injected harmonics when there is very little voltage distortion [9]. Since the EV charger is


KVA RATING 37

Turns Ratio 7260, 120/240

Material of the Core Grain oriented silicon steel

Thickness of Lamination .009 inches

Net cross sectional area of the core 28.9 sq. inches

Path length 25.74 inches

Magnetization Saturation 20.0 Kilogauss



























Figure 7.4. Domestic Loads Model


known to cause very little voltage distortion, the harmonics injected by the charger into the

line can be modeled by a composite current source. The behavior of the charger can thus be

determined by a second order system.


LC +RC-dvc + v =v (1)
dt2 dt



where v, the transformer secondary voltage

vc the capacitor voltage

L the equivalent battery inductance

C the equivalent battery capacitance


R the equivalent battery internal resistance








The net harmonic current from the composite current source can be expressed as:

50
iH Imk sin kcol (2)
k=2

where it, is the net harmonic current

k is the order of the harmonic

Ik is the amplitude of the k th harmonic current

The SABERTM software package was used to simulate the system, both transient

and steady state operating conditions, including RMS values, THD etc can be obtained

using SABERTM. In the battery model the L is small and hence can be neglected.

Therefore the charging time constant (T) for the battery can be given by:



T = RC (3)




7.3. Simulation Results




In the computer simulation two types of loads have been considered, viz. the

lighting loads and the EV load. The simulation results with two houses and a single EV

load are presented. The lighting load current is taken as a pure sinusoidal current. The

current harmonics injected into the mains by the EV charger are shown in figure 7.3. This

is based on the experimental results obtained by testing the EVI charger in the

laboratory. The harmonic components of the charger input currents (figure 7.5) are

supplied by the composite current source. The charger input current in the simulation is








controlled by varying R. We know from the lab test data that the maximum current that

the charger draws (for fully discharged battery) is approximately 30 Amps.

By modeling the GM EVI charger as a combination of a full wave bridge rectifier

with a resistive and inductive load and a composite current source, simulation results very

close to the lab results were obtained. A detailed discussion on the comparison of the

laboratory results, the simulation results and the field-recorded data has been included in

Chapter 7.


The simulation runs of the distribution system made it possible to check the various

load scenarios in the distribution system. The system design can now be checked for an

arbitrary number of EVs connected to it. This data obtained from simulations is important and

makes it possible to optimally design/redesign future distribution systems. This includes

replacing and/or adding transformer units in the existing systems in order to make them able to

withstand the new load conditions caused by the EV load and the harmonic distortion

associated with it.


Based on the field and laboratory data the EV charger has been successfully modeled.

The simulation results are close to the laboratory results obtained during the actual charging of

the EV battery. The waveforms on the input side of the transformer have a low order current

harmonics.


7.4. Simulation of Short Charging Time Chargers


The computer simulations were used to check scenarios when the charging time

for charging the battery is much shorter. The GM EV1 charger takes about three and a

half hours to charge a fully discharged battery (from about 15% of full charge). The long








charging time is a drawback with the EV and research is being done to reduce this

charging time substantially. The fast charging time chargers are not commercially

available and their behavior in terms of harmonic effect is not known. The computer

simulations have been used to provide an estimate of the effect of the fast charging EV on

the distribution system. The results of these simulations are based on the researcher

intuition and experience of working with highly non-linear loads. The fast charger study

is aimed at providing an idea about how such a technology, when developed can effect

the distribution system. These results need to be verified once the fast charging

technology becomes available.



7.4.1. Simulation Methodology For Fast Chargers



The charger model for the fast charging charger is essentially the same as that for

the normal charging time charger. The difference however is only in the values of the

resistor and the capacitor that represents the battery. Also the amplitude of the harmonic

currents has been increased proportional to the decrease in the charging time. Thus the

fast charging time charger has also been simulated as a combination of a full-wave bridge

rectifier and a harmonic current injection source. As mentioned earlier, these simulation

results are based on intuition and experience of working with non-linear loads. These

results can only be verified once fast charging technology becomes available. Figures 7.6

through 7.8 show the current waveforms and their respective FFTs for simulated

scenarios where the charging time is reduced as indicated in the figures. The distribution

transformer derating corresponding to this charging scenario is shown in Section 8.2.



































Figure 7.5 Transformer Current for the Fast Charging Charger with Charging
Time of 45 minutes.








Trmuuformar Cun Lr


40.0








o1ask 1~k e.75k 1.0K I ItK 1.Ek



Figure 7.6 FFT of the Transformer Current for the Fast Charging Charger with
Charging Time of 45 minutes.


(A): I(S
Trmfanner Curuwmn













0.o




































Figure 7.7. Transformer Current for the Fast Charging Charger with Charging
Time of 1.2 Hours.


30.0.


ZoA


0.6k 1.0k
f"m)


Mlg(A) : f(Hz)
FFT of the Trunsformw Curnw














1.5k 2.0k


Figure 7.8. FFT of the Transformer Current for the Fast Charging Charger with Charging
Time of 1.2 hours.


(A) : t()
Transonner Cumnrt













CHAPTER 8
COMPARISON OF FIELD, LABORATORY AND SIMULATION RESULTS


In this section, the field data and laboratory test results are compared with those

done by computer simulation to present the effect of EV loads on the power distribution

transformer.

The aim of this project is to study the impact of the EV load on the power distribution

system. In summary, this goal has been accomplished by recording the field data with and

without the EV load, recording the behavior of the EV charger by the laboratory testing of the

GM EVl charger charging a fully exhausted battery. The data acquired from the field provides

the information about the characteristics of the distribution loads when EV load is a part of the

load and when there is no EV load. The laboratory testing of the charger provides the input

current characteristics of the charger and the line current harmonic content of the charger.

Using the data from the laboratory test, the computer simulation model of the EV

charger has been built. The details of the computer model have been discussed in Section

8.1. Thus the laboratory testing of the charger forms the basis on which the computer

model is built. The charger computer model has been incorporated into the overall model

of the distribution system. The distribution system consists of the transformers, feeders,

different types of normal domestic loads and the EV loads.

First of all, in this section, the accuracy of the computer model of the EV charger

will be established. This has been achieved by comparing the laboratory test results of the








EV charger with the computer simulation results. Second, the similarity of the field-

recorded data with the simulation results will be discussed.



8.1. The Field Data



The field data was recorded in the form of one-cycle snapshots of current and

voltage as well as power factor at the distribution transformers at four different locations

in South Florida and then transmitted to a computer over the telephone line. The detailed

description about the data monitoring and acquiring methodology are mentioned in

Section 4.2.



8.2. Laboratory Test




For the laboratory test the EV charger is selected as a 240 volt, split single-phase

device as described in the previous section. The voltage and current input into the charger

were monitored and recorded over the full charging cycle for the EV, which was nearly

completely discharged to a level of 5% of full charge. Electrical measurements were

taken using a power monitor automatically every minute over the entire charging period

of about three hours.

From the measurement it was found that the current supplied to the EV charger

varies throughout the charging cycle. The charger current variation is shown in Figure

5.9. The actual amplitude of each of the harmonics in the charger current does not vary

significantly during the charging cycle.








8.3. Simulation Results



The power distribution system consisting of the feeder, distribution transformer,

the EV charger, together with domestic loads (mainly the lighting lamps) was simulated

using SABERIM. The models of various components used for simulation are based on

both field data recorded and the laboratory data. To simulate variation of input current

and distortion at different instants over its charging period, an adjustable resistor, as well

as a composite current source which is based on the FFT analysis of the laboratory test

data, is introduced in the loop. By varying the resistance, the identical situations to those

in the laboratory test can be established. Therefore, current waveform snapshots over a

whole cycle of 60 Hz are obtained and then compared with those recorded in the

laboratory test.



8.4. Comparison




The field data, laboratory test data and simulation results are compared in terms of

the waveform snapshots with the same magnitude of charger input current for all three

cases viz. the field data, the laboratory test data and the simulation results. The aim is to

establish the reliability of the laboratory test data; the simulation model and simulation

results by comparing them to each other and the field recorded data. Some of the

waveforms that have been included in previous chapters have been put in this section for

convenience of comparison.















FeI Recorded Currentwiti EV Load


Mi: -302734
Ag: 1527
Ats: 302734
RMS: 181012
CF : 167246
IF:1 85


5 10 15 20



(a)

Ha ont ATf flkide oftie F bRAecordeaurrent

Fi :179963

TD :0
To :10.79
RMh : 194302
.-- ----- ---- ------- RMSi :0
RMSih: 1:0
RMS : 18009
P511M 27.771
T :289335
S:53723







'. I..'.' 3120



(b)


Figure 8.1.
a.) Field Recorded Current

b.) Waveform FFT of the Field Recorded Waveform.


^^ L^TbnitBTiE

20 *- -
<

10--V---


^ o-----
4


125-





L-075-
U
050-


025-











M34 'fZl
Mr. 1l".C
Ag: 533581

CF : 1.7533
FF :114432


10
Tile5


(a)





Hanmormt rc ntoftie iefiRecordi Qarent --


TM :12B78
S TD :0
T :123278
RMs: 197152


RMS :1SA95

-- 5


r 5518.64



25 -



120 720 1320 1920 2520 3120





(b)


Figure 8.2.
a.) Field Recorded Current Waveform.
b.) FFT of the Field Recorded Waveform in Figure 8.3.


-," _____




C 0
ia r^ l /a~




L.


150


125


.07


* 075-


0
0


0









40














0 10 1I 20 0i
Time (mS)


Figure 8.3. Laboratory Test Transformer Current at 11:40 AM



















Time (ms)

Figure 8.4. Laboratory Test Transformer Current at 12:36 PM






















Figure 8.5. FFT of the Current Waveform of Figure 8.6.

























muld a Chrgr RCur wth corr-pong wFFT
(correpondnig to Lab Te CuMrent 12:0)
Mwg(A): ;f(H

2Z




0.5k 1.0k Aik 2.0k 2.k 0.0k


SO IKhorNt-xjn2)





0Oj2 0028 0.06 45 004
A*m


Figure 8.6. The Simulated Current Waveform and the Corresponding FFT


















The SkAuld Current Wvetornm nd Nb FT
(Conoepondnngto 11:40 urent Wefnm of the Lab TeM)
1. ___~_________ MCgA) t:b)





S | I I
O.k Ik 1.6k 8k 2.k Mk

~~~____~~_____(A) ,t(*)
210 |l(horwttn)





0.0 0m01 0. 0.4 0.04






Figure 8.7. The Simulated Current Waveform and the Corresponding FFT








The waveform of figure 8.lshows the field recorded current for the distribution

transformer feeding two homes (FPL_1I). This is the distribution transformer current when the

load on the transformer contains EV load. Figure 8.2 is the FFT of the waveform in Figure

8.1. The observation after comparing the current waveform of Figure 8.1 and its FFT in figure

8.2 with the laboratory recorded waveform in Figure 8.5 and its corresponding FFT in Figure

8.6 are presented in a tabular form in Table II. Both the currents (of Figure 8.1 and 8.5) are

then compared to the waveform of Figure 8.8. Figure 8.8 shows the simulated transformer

current. Also compared are the FFTs of the three waveforms. (Figure 8.2, 8.4 and 8.8).





Table. II



Laboratory Test Simulation Results Field Data

(Figure 8.5 and 8.7) (Figure 8.8) (Figure 8.1 and

8.2)

Peak Value of 30 Amps 29 Amps 27.82 Amps
Current
RMS Value of 19.8 Amps 20.2 Amps 18.1 Amps
Current
Crest Factor 1.515 1.435 1.537

THD 14.1 11.2 12.8



Holding the field-recorded data as reference the relative percent values of the

laboratory and simulation results is shown below in Table III.










Table. III


Looking at the waveforms figure 8.1 through figure 8.8 and Tables II & III, we

note that the field recorded and the simulated current contain the domestic load current.

The laboratory test waveform however is only the charger current. In the simulation the

domestic load is taken as linear load (Figure 6.4) while in reality (Field recorded) the

domestic loads may contain some non-linear loads. Thus the THD of the field current is

higher than the THD of the simulated current. The laboratory test current consists of only

the charger current.


Laboratory Test Simulation Results

(Figure 8.5 and 8.7) (Figure 8.8)

Peak Value of Current 7.8% 4.24%

RMS Value of Current 9.39% 11.6%

Crest Factor -1.43% -6.636%

THD 10.15% -12.5%













CHAPTER 9
DERATING OF DISTRIBUTION TRANSFORMERS

9.1 Method of Calculation


9.1.1 ANSI/IEEE C 57.110 Method


The method used for calculating the derating of the distribution transformer is the

ANSI/IEEE C 57.110 (see Section2.3.1). Using the method described in Section 2.3.1, the

derating has been calculated as a function of the current drawn by the charger. The derating

curve is shown in Figure 9.1. In the calculations the derating is calculated taking the charger

current as the base. So fh in equation 2.1 is calculated dividing each harmonic component by

the amplitude of the fundamental charger current. Thus, in this method the value of fh for the

fundamental component is always 1. This is the procedure employed by ANSI/IEEE C57.110

method. From the curve of Figure 9.1. we get the minimum and maximum limits of derating.





Table IV. Maximum and Minimum Derating Using ANSI/IEEE C57.110












Derating of Distribution Transformer


.-6 .4p SP3 .<^ .<^ .61 SP .(, 0(Z OZ .Q .Q .(


Time of Charging

Figure 9.1. The Transformer Derating Curve for the Complete Charging Cycle








9.1.2 Modified ANSI/IEEE C 57.110 Method





The ANSI/IEEE C57.110 method calculates the derating of the transformer

based on the charger fundamental current. In reality the transformer in a distribution

system is supplying other loads. Also from figure 9.1 it can be seen that the

transformer derating at 13:53 hours is as low as 20%. This calculation is done taking

the charger fundamental as base, and from Figure 5.9 we see that the charger current at

13:53 hours is just about 5 Amperes. This would form only a small percentage of the

total load on a distribution transformer of 37 KVA rating. Thus, although the derating

based on C57.110 method is very high the effect of the small current on the

transformer would be much less significant. In the modified ANSI/IEEE C57.110

method the calculation of transformer derating is done based on the maximum current

carrying capacity of the transformer. Since the transformer derating should be a

concern when the load on the transformer is full or close to full load this method

calculates the derating based on full load current of the transformer. The derating of a

37 KVA transformer using the modified method is shown in Figure 9.2.





92. Transformer Derating with Fast Time Chargers





Based on the simulation results of section 6.4, the modified ANSI/IEEE

C57.110 method has been used to calculate the derating of distribution transformers.








The derating has been calculated based on the 37 KVA distribution transformer whose

specifications are given in Section 6.1.3. Figures 9.3 and 9.4 show the transformer

derating for the two simulated load scenarios.





9.3. Conclusion





The derating of the distribution transformer based on ANSI/IEEE C57.110 method is

significant but it does not consider the current carrying capacity of the transformer when

calculating the derating. The modified approach calculates the derating of the transformer

based on the full load current of the transformer. It is obvious that based on this method we see

that the GM EVI charger does not significantly derate the transformer. One thing that may be

noted here is that the ANSI/IEEE C57.110 method is generally applied to find the derating of

large power transformer and gives highly conservative results for distribution transformers. It

is strongly recommended that the equation for derating of distribution transformers be

experimentally determined with KVA ratings of concern to FPL. (See Chapter 10). The

derating for the fast charging time chargers as seen from figures 9.3 and 9.4 is much higher. It

is expected that the fast chargers will be more prevalent than slow chargers if consumer

preference is not modified by the utility industry.












Derating of Distribution transformer


S-T- of Charging
Time of Charging


Figure 9.2. Transformer Derating using the Modified ANSI/IEEE Method


1.0001
1
0.9999
0.9998
0.9997
0.9996


)) (( i









Derating of Transformer





11%


)erated
Value
89%


*Derated Value





a.




Derating of Transformer


)erated
Value
59%


gDerated Value




b.


Figure 9.4. Derating Corresponding Short Charging Times a.) 1.2 Hours b.) 45 minutes













CHAPTER 10
CONCLUSIONS

10.1. Summary and Conclusions.


The objectives of the project outlined in section 1.3 have been achieved. The

behavior of four distribution nodes (FPL_1I to FPL_4) has been monitored. The load

behavior with and without the EV load has been recorded. Volume II of the final report

contains samples of the recorded data with an executive summary. By performing the

laboratory test on the GMC EV 1 charger charging a fully discharged battery, the charging

characteristics of the charger has been determined. The charger takes about 3 hours to

charge a nearly fully discharged (15% charge) battery to full (100%) charge. The

harmonic spectrum of the transformer current charging the discharged battery has been

determined. The snap shots of the current waveforms have been recorded and studied.

Some of these waveforms have been included in this report. Based on the laboratory data

a computer model of the charger and the distribution system has been developed.

Comparing the simulation results to the laboratory results and the field data has validated

the simulation results. The simulated currents are roughly 4 to 10% away from the field

recorded currents, in their magnitude and distortion. Based on the field data, laboratory

results and the simulation results the derating of the distribution transformer has been

calculated. Based on the field data, laboratory results and the simulation results the

derating of the distribution transformer has been calculated. The two methods suggested








for the calculation of derating have been presented and compared. The maximum and

minimum limits of derating have been suggested. Using the first method the maximum

derated value is 80% and the minimum value is 1%. Using the modified method the

derating remains close to 99 % throughout the charging cycle. The transformer derating

results corresponding to fast charging (45 minutes and 1.2 hours) chargers has been

presented.

A Microsoft Excel based program has been provided to FPL. This program shall

enable the user to find the derating of the transformer using the harmonic content

information of the load on the transformer.



10.2. Scope For Further Research



The transformer derating has been calculated using the IEEE /ANSI

method. The formula is an empirical formula. To get a highly accurate measure of

the derating of transformers, laboratory tests should be performed. These tests can

be performed in the Power Laboratory in the Electrical and Computer

Engineering Department at the University of Florida. Using the real time program

R4TM, the computer can generate low voltage controlled harmonic signals. These

signals will be then amplified using the power amplifier in the laboratory and fed

to the transformer under test. Temperature sensors placed selectively on and

inside the transformers will sense the rise in the core and winding temperatures.

The change in the temperature rise with a change in harmonic content from pure

sinusoidal to very high harmonic content can be gauged. The transformers can be





91


tested for different harmonic scenarios and loads. Thus an exact formula for the

derating of transformers can be evolved. As the harmonic generation can easily be

controlled using the computer, affect of individual harmonics can also be studied.

This information can be particularly useful to FPL for the design of harmonic

filters that may need to be installed in places where the load is highly sensitive to

harmonics. Also, at sites where the load is known to generate a significant amount

of a particular harmonic, these tests can provide a strong basis for calculating the

derating of equipment and thus the design of appropriate transformers.














REFERENCES


1. J. Pages, "On the quality of public electricity supply" Revue Francaise de I'Electricite,
vol.41, no.223, p. 14-21, 1968.

2. C. H. Hilger, "What is the quality of electric power?" Elektroteknikeren, vol.68,
no.19, p. 418-22, October 7, 1972.

3. H. Heikkila, "Harmonics in a power system" Saehkoe, vol.49, no.1, p. 27-30, January
1976.

4. G. L. Johnson, A. N. Mortensen, "A power system digital harmonic analyzer" IEEE
Trans. Instrum. Meas. (USA), IEEE Transactions on Instrumentation and
Measurement, vol.37, no.4, p. 537-40, 1988.

5. C. K. Duffey, R. P. Stratford, "Update of harmonic standard IEEE-519: IEEE
Recommended Practices and Requirements for Harmonic Control in Electric Power
Systems. Industrial Applications Society 35th Annual Petroleum and Chemical
Industry Conference pp. 249-55, 1988.

6. P. Meynaud, G. Hensman, J. L. Javerzac, L. Lagostena, M. J. Lahtinen, F. Lionetto,
W. Monbauer, P. Herager, A. Robert, J. Smid, G. Allen, C. Britten, A. Emanuel,
K. Murotani, P. Ribeiro, P. Wright, Yu Zhelezko, "Equipment producing harmonics
and conditions governing their connection to the mains power supply." Electra,
no.123, p. 20-37, March 1989.

7. T. S. Key, R. J. Ferraro, R. W. Taylor, "Electric power for the 21st Century: changing
loads-changing needs." Power Quality '89. Official Proceedings of the First
International Conference, p. 613, 58-73, 1989.

8. R. J. Gilleskie, "Power quality and emerging electrotechnologies their relationship,
impact and suggested precautions." Power Quality '89. Official Proceedings of the
First International Conference, p. 613, 108-23, 1989.

9. W.E. Kazibwe, R. J. Ringlee, G. W. Woodzell, H. M. Sendaula, "Power quality: a
review." IEEE Computer Applications in Power, vol.3, no.1, p. 39, January 1990.

10. M. S. Padbidri, D. G. Pantulu, "Studies on harmonic pollution in power systems."
Power Systems for the Year 2000 and Beyond. Proceedings of the Sixth National
Power Systems Conference, p. xvi+704, 31-5, 1990.








11. A. Domijan, Jr., E. Embriz-Santander, "Harmonic mitigation techniques for the
improvement of power quality of adjustable speed drives (ASDs)." APEC '90. Fifth
Annual Applied Power Electronics Conference and Exposition. Conference
Proceedings. p. 839, 96-105, 1990.

12. D. Ollivier, "Simulation of harmonic disturbances (in power systems)." Bulletin de la
Direction des Etudes et Recherches, Serie B, no.1, p. 1-31, 1991.

13. J. K. Phipps, J. P. Nelson, P. K. Sen, "Power quality and harmonic distortion on
distribution systems." 1991 Rural Electric Power Conference. Papers Presented at the
35th Annual Conference. p. xi+164, Dearborn, MI, USA; 28-30 April 1991.

14. D. A. Van Den Eijkel, D. A. Marshall, "The simulation and measurement of system
harmonic impedance." CIGRE. Proceedings of the 33rd Session. International
Conference on Large High Voltage Electric Systems, p. 2 vol. 2772, 36-302/1-4
vol.2. Paris, France; 26 Aug.-I Sept. 1990.

15. H. M. Beides, G. T. Heydt, "Power system harmonics estimation and monitoring."
Electric Machines and Power Systems, vol.20, no.2, p. 93-102. March-April 1992.

16. B. M. Hughes, J. S. Chan, D. 0. Koval, "Distribution customer power quality
experience." Conference Record of the 1991 IEEE Industry Applications Society
Annual Meeting. p. 2 vol. 1910, 1556-63 vol.2.Dearbom, MI, USA; 28 Sept.-4 Oct.
1991.

17. R. R. Verderber, 0. C. Morse, W. R. A.bir, "Harmonics from compact fluorescent
lamps." Conference Record of the 1991 IEEE Industry Applications Society Annual
Meeting. p. 2 vol. 1910, 1853-8 vol.2. Dearborn, MI, USA; 28 Sept.-4 Oct. 1991.

18. R. Christiansen, "Effects of high levels of harmonics from lighting equipment and
systems." Conference Record of the 1991 IEEE Industry Applications Society Annual
Meeting (Cat. No.91CH3077-5), p. 2 vol. 1910, 1859-62 vol.2. Dearborn, MI, USA;
28 Sept.-4 Oct. 1991.

19. K. Johnson, R. Zavadil, "Assessing the impacts of nonlinear loads on power quality
in commercial buildings-An overview." Conference Record of the 1991 IEEE
Industry Applications Society Annual Meeting (Cat. No.91CH3077-5), p. 2 vol. 1910,
1863-9 vol.2.Dearbom, MI, USA; 28 Sept.-4 Oct. 1991.

20. R. H. Kitchin, S. Etminan, "Computer simulation of harmonic waveform distortion,
and comparison with measured results." IEE Colloquium on 'Simulation of Power
Systems', p. 62, 4/1-7 London, UK; 2 Dec. 1992.

21. R. S. Wallace, "The harmonic impact of variable speed air conditioners on residential
power distribution." APEC '92. Seventh Annual Applied Power Electronics




Full Text
xml version 1.0 encoding UTF-8
REPORT xmlns http:www.fcla.edudlsmddaitss xmlns:xsi http:www.w3.org2001XMLSchema-instance xsi:schemaLocation http:www.fcla.edudlsmddaitssdaitssReport.xsd
INGEST IEID EDKGIPLV4_XZ3UL5 INGEST_TIME 2013-02-07T16:44:33Z PACKAGE AA00012982_00001
AGREEMENT_INFO ACCOUNT UF PROJECT UFDC
FILES



PAGE 1

IMPACT OF ELECTRIC VEHICLE LOADS ON ELECTRIC POWER DISTRIBUTION SYSTEMS By TARIQ ASLAM BUCHH A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY UNIVERSITY OF FLORIDA 1998

PAGE 2

In tlie Name of God, tlie Most Beneficent, tlie Most Merciful. All the praises and thanks be to God, the Cherisher and Sustainer of the Worlds. The Most Beneficent, the Most Merciful. The Only Master of the Day of Recompense. You (Alone) we ,vorship and You (Alone) we ask for help. Guide us to the Straight Way. The way of those on whom You have bestowed Your Grace and not of those who earned your anger, nor of those who went astray. Translation of Opening Chapter o_f the Ho(y Qi,ran

PAGE 3

ACKNOWLEDGEMENTS I I acknowledge the conviction of those who believe in the Truth and have the courage to stand for it. Having said that, I would like to mention the names of those who have been there when I needed them. A word of appreciation for my adviser Dr Alex Domijan Jr. for agreeing to chair my advisory committee and to Dr. Dennis Carroll, Dr. Khai D Ngo, Dr. Herman Lam and Dr. Barney Capehart for being on the committee. My motl1er, who has been a constant source of inspiration and drive, my father, who is an embodiment of patience my wife, Falaknaz who sacrificed a part of her dreams and aspirations and is taking excellent care of our kids, Shaima and Ziad, my loving brothers, Mehmood and Khalid and my best friend, Ashfaq. There has also been a huge treasure of friends and well-wishers who were there to support and pray for me whenever the going got tough. A special word of thanks for Dr. Song for his constructive criticism and suggestions. Florida Power and Light Company, provided the funding and technical support and I appreciate it ... 111

PAGE 4

TABLE OF CON T ENTS A CKN O WLEDG EMENTS ....... .. .... . . ... .. .... ... .. ... .. .... . .. . .... ........ ....... iii ABST:RACT ....... . .... .... ...... ... .. . ... . .... ... .... .. .. ................. .. .. .. .... .. . ..... .. vi CHAPTERS I I N TROD UCT I ON .. . . ...... . ..... ...... . . .. .......... .. .. ..... .... .. ... . . ... . ...... . 1 1 1 Background ........... ... ..... .. .. .. .. . . .. ... .. ..... ... ..... .......... ...... ..... . . ... ... .. ..... . .. 1 1 .2. Project Justification ................................................... .. .. ... ......... ............. ...... 3 1 .3. Objectives and Goals ..... .. . .. . . . . .. . .. . .. .. ..... .. .. ... .. .. . .. .... .... . ...... .... ... . . .. 4 1.4. Major Con tributi ons ....... .. .... .. .... . .... . . .. .. . ................ .. ............. . .... ..... ... . 4 I 4 Re searc h Pha ses . .. . .. .. .. .... .. ....... . .. ..... ........ .. . ..... .... ........... . .. ................. 5 2 LITERA. TURE SURVEY .. .... ...... . ..... ... .. . . .... .. .. .. .. . . . .... .. .. .. ... ... .. .. ... ... . ............... 7 3 POWER QUALITY AND HARMONIC DISTORTION .. ..... .. .......... .. ... . .... ... ... 10 3.1. Harmon ics .. . . ........ ..... .. . ... .... .. ........ . .. .. .. .... .. .... .. .................................... 11 3 .2 Estimation of H arrnonic Di stortio n . . .... .. .. ..... .. . .. .. ...... .. . ........ ...... .. .. ..... 1 2 3.2 1. C r est Factor .. . .. ...... ... .. .. .... . .... ......... . .. .. ............................................ 12 3 .2.2. Percenta ge Total Harmonic Di s tortion (THD) .. ....... . . .. .. . . . . .... ........ 13 3.2.3. The K-factor .. ...... ....... .... . ... .. .... .. . .. ... ................. ......... . ........ .. ........ 13 3.3. Transforrner Der ati n g .. .... .. .... ....... . ........ .. .. ............ .. .................... . . ....... 14 3 .3 1 Tran sfo rmer Capability E quival e nt Ca lculation ........ . ....... ... .. .. .... ... 15 3.3.2. Effect of Current Harmonics on the Transforrner Neutral Current ....... 17 4 FIELD DATA RECORDING, REPRESENTATION AND OBSERVATION . ......... 19 4.1. Introducti on . ..... . .. . ......... . .. .. . ...... .............. .. .. . ..... . .... ...... .... .. ....... . .. .. ..... 19 4.2. General Cons ideration s In Dat a Acquisition And Pro cessing ................ ........ 20 4.2.1. Data Collection ......... .. .. .. ................... ...... .. .......... . . . .. . ... ............. 21 4.2.2. Data Recording .... .. .. . .... .. ... .. ....... .... ... . . .... .. .. ... .. ... .. ... . ... ............. 22 4 .2 3. Data Preparation . .......... ... ........... .. .. .. ......................................... ... ... 22 4. 2 .4 Data Qualification .. ... .. .... . ... ... . ... . .. ...... . . . ............ . .. ....... ......... 23 4 .2. Dat a Analysis . .. .. . .. .. ..... . . .. . .. . . ............................................................... 24 4.3. Field Data Monitoring .. .. .. ....... . ... .... .. ........ . .. . .... .. . ..... ... ... ... . ..... . .. ........ 26 4 .3 .1. D ata Acquisition Methodology .. ..... . ... . .......... . ............ . . ... .. ... .... 26 4 3 2. D ata Monit or ing In st rum e nt ....................... .......................... . . ..... .. . ... 26 4 .3.3 Data Management ....... .. . ...... . ... ... ..... . .... .. . . .. . . . ..... .. .................... 27 l V

PAGE 5

5 DATA REPRESENTATION AND ANALYSIS .... .......... .. . .. . ... ........ . . ... .. ........ 30 5 1. Introdu ction . .. ... ............ .. ....... .. .. ... ... . .. ...... .. ......... .. ..... ..... .. ..... ... ...... .. .. .. .. .. .. . . . ... ..... ... 30 5.2. Data Trend s a nd Characte ri stics .... . ..... . .. . .. ...... .. . .. . .. .... . .. .. . ... ........... 32 5.3. Data Ob se rvation ............ ..... . ....... .. . ... .. ..... ........ . .. .. ....... . .. ........... ................. . ... 33 5 4. Load Character istic s ............... .. . .. . .... ...... ....................... .......... . . ......... . 3 4 6 LABO RA TORY TESTING OF THE EV CHARGER .. ........... .. .... . .. .... . . .. ........ 4 7 6 .1 Te st in g Pr oce dure and Re s ults ........... .. ...... .. .......... . .... .. ............ . ......... .. .. 4 7 6.2. Observation of Laboratory Re s ults ..... . . ..... .. .. . .... ......... .......... ..... ..... .... 48 7 COMPUTER SIMULATION .... .. .. .... .................................. .. . .. . . .. . .. .... ..... .. ...... 6 1 7 1. Component Model s . .. ... ............. .. ... .. . . . ........ .. . ....... .................................... 61 7.1.1 S ys t em Model ..... .......... .. .. . . .... .... .... .. . .. .. . .. ....... ........ .. . . ... .. . ........ 61 7 .1.2. Electric Vehicle C h arge r And Battery Models ..... ... . .. .. . . ... .. ..... ... 64 7 1. 3. Transfo rm e r Model .... .. . .. . .. . . . . .. .. .. .. ... ..... ..................... . .. .. ...... 64 7.1.4 D o me s ti c Load Model .. .. ...... . .. . .. .. .. .. .... . .. ................... ...... .. . .. ... 66 7 .2. Mathematical Model Of The EV Charger .... .. . .. ........................ . ..... .. . . 66 6.3. Simulation Results ... .. .. . . ... ...... .......... .. .. .... ............ . .. . ... . ... .. .... . ... .. .. ...... 68 7.4 Simulation of Short C h arg in g Ti m e Chargers .. .... . .... ........... ..... . . . ... .. : . 69 7 .4.1. Simulation Methodology For Fast C har gers . .... . . ... ..... .... ... . .......... 70 8 COMPARISON OF FIELD LABORATORY AND SIMULATION RESULTS ....... 73 8 I The Field Data .... ... .. . .. ........... ....... . .. .... ..... .. .. .. ........ ............ . . .............. 74 8 .2 Laborator y Te s t .... .. .... .. .... .. . ... . .. . . .. . ..... .. .. ...... ............ .... .... . . .. ... ... . . 7 4 8 3. Simulati on Re s ults .. ....... ................... . . .. ..... .... .. ...... .. . .. . . .. .. .... ..... ......... 75 8 4. Comparison ...... .. ........... ..... . .. .. . ... . .... ..... .. . .. . .. .... . ..... . .... .. ....... .. .. ....... 75 9 DERA TING OF DISTRIBUTION TRANSFORMERS . .. . ........ ......... .. ... ... . .. ....... 83 9 1 Metl1od of Calc ul atio n .. .. ............. . .. . ... .. .... .. .. . .. .... .. . ............... ........... .. .. 83 9 .1.1 A N SI/IEEE C 57 110 Method .... . ........ . .. .... .......... ..... ... . . ... . .. .. .. .. 83 9.1.2 Modified ANSI/IEEE C 57.110 Method .. . . . .... . . .. . . . .. ..... . ... .. .......... 85 9.2. Transformer Derating w ith Fa s t Time C h argers . .. .. . .. .... ......... .. . . . .. ... .. . .. 85 9.3. Conclusion ..... .. ... .. .. ................ .. . .... . . .. .. ........ .... . . .. .. . . .. .. .. .. . .. . ......... 86 10 CONCLUSIONS ....... ..... . ..... .. ....... .. ......... .... .. . .. .... ... .................. ...... . ....... ... .. . . 89 10. I .Summary and Co nclu s i ons ....... .. ... .. ... .. ... . .. .. ... ... .. ........ ...... .. 89 10 .2. Scope For Further Research .. .. . . ... ........ . .. .. .. ....... .. .. . . .... .. .. .... .. . ... .. ... 90 REFERENCES .... ... . .......... ......... .. . . ... ..... ....... . .... .... ........ .... . .. 92 BIOGRAPHICAL SKETCH ... .... ... ... .. ............... ... ..... ....... .. .. .. ......... 100 V

PAGE 6

Abstract of Dissertation Presented to the Graduate School of the University of Florida in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy IMPACT OF ELECTRIC VEHICLE LOADS ON ELECTRIC POWER DISTRIBUTION SYSTEMS By Tariq Aslam Buchh August 1998 Chair1r1an: Dr. Alexander Domijan, Jr. Major Department: Electrical and Computer Engineering The voltages and current s 011 four distribution sites (FPL l to FPL_ 4) have been monitored. The load behavior with and without the EV load has been recorded By performing laboratory tests on the GMC EVl charger charging a fully discharged battery, the characteristics of tl1e charger have been deter1nined. The ha1"111onic spectrum of the transformer current charging a fully discharged battery has been deter1nined. The snap shots of the current and voltage wavefor1ns have been recorded and studied. Based on the laboratory data a computer model of the charger and the distribution system has been developed. Comparing the simulation results to the laboratory rest1lts and the field data the simulation results have been validated. The simulated currents are roughly 4 to 10% away from the field recorded currents, in their magnitude and distortion Based on the field data, laboratory results and the simulation results the derating of the distribution VI

PAGE 7

transformer has been calculated. Tl1e two methods suggested for the calculation of derating have been presented and compared. The maximum and minimum limits of derating have been suggested. Using the first method the maximum derating is 20% and the minimum value is 99%. Using the modified method the derating remains close to 99 % throughout the charging cycle. (It must be noted here, that a derating to 20o/o means that the transformer can safely handle only 20o/o of its full rated load. Likewise a derating to 99% means that the transformer can handle 99% of its full rated capacity. Thus, a transformer derated to 20% is highly derated and a transformer derated to 99% is slightly derated.) After establishing the reliability of the simulation, by comparing the simulation results with the field data and laboratory results, a charging scenario with shorter charging times (45 minutes and 1.2 hours) has been simulated The distribution transformer derating corresponding to these short-charging times has been calculated and presented. Furtherrnore, recommended scope for further research has been suggested which among other things, involves validating the derating formula for distribution transformers. . Vil

PAGE 8

CHAPTER 1 INTRODUCTION 1 1 Background From the state of being a concept vehicle, a possible means of future transportation and a dream, the Electric Vehicle (EV) has become a reality. But even after the fust major commercial launch of the EV 1 by General Motor Corporation, the future of EV s is still under a cloud of uncertainty and some unanswered questions. The popularity of EV s as a future mode of transport can be attributed partly to exhausting reserves of natural energy resources and environmental considerations. These and many other factors have boosted the research and development of the EV technology and a possibility of it being accepted as the popular mode of ground transport in the not so distant future. There has been a great boost to EV research because of government and industrial funding. This research has been mainly aimed at improving the battery technology, charging techniques, weight, comf art and reliability. For the environmentalist, the EV s' noiseless operation compliments its pollution free pertormance, and, for the customer, the EV is a good perforrr1ance machine with the downside of short battery life, long recharging time and high initial cost. The EV 1 manufactured by General Motors Corporation is powered by a Delco Valve Regulated Lead Acid battery pack, which contains 26 12-volt modules, capable of carrying 16.2-kilowatt hours of energy. The batteries can be recharged from O to a full 1

PAGE 9

2 state of charge in approximately 3 hours using a 240-volt inductive Magne charge charger and about 15 hours from a 120-volt portable charger. Inductive charging is achieved without metal-to metal contact. The alternating current powers a device that produces a fluctuating magnetic field in a paddle that is inserted into a port at the front of the car. This induces an alternating emf in a winding inside the car. This alternating emf is then rectified and fed into the batteries. The inductive charging system is inherently safe even under hazardous climatic conditions, since the windings are safely encased in a plastic cover. Typically this system can handle power levels from 1.5 to 25 kW with overall efficiency of 90% while the power transfer frequency is 40-350 kHz [l ]. The specifications of the EV, as published by General Motors [PrEView Drive pro gram literature], are summarized as follows: 1. 137 H.P. three-phase ac induction motor. 2. 16.8 kWh lead-acid battery. 3. Inductive coupled standard charger (220 volts, 30 Amps) 4. IGBT power inverter module. 5. Range: 90 miles highway, 70 miles EPA city. 6. Acceleration from O to 60 mph in 8.5 seconds. 7. Electronically regulated speed of 80 mph. 8. Charging time from 85% depth of discharge is 2 to 3 hours. 9. Curb weight is 2970 pounds. While the popularity and widespread use of EV s on one hand will take ground transportation into a new era, on the other hand it can put additional and unprecedented strains on power supply and distribution systems.

PAGE 10

3 For power systems, the EV loads pose twofold problems. These are: 1. Increase in power demand 2. Decrease in power quality. The power distribution system looks at the EV as an additional load of about 7kw/EV (based on the GMC EV 1 ), which was almost certainly not foreseen when most of the present distribution systems were designed more than a few decades ago. The fact that the EV charger load is a non-linear load with power electronic devices makes matters worse, for it is a source of harmonics and distortion in the distribution system. These harmonics which get injected into the system have a significant impact on the distribution system in terms of losses and overall perfor1nance of the power system. The non-linearities and harmonics caused by EV chargers can adversely affect other sensitive l oads that may be connected to the distribution system. The power quality can become a major concern if the battery and charger technology is developed such as to put an extra strain on the distribution system. For example, research is now being aimed at developing batteries and their chargers which would make it possible to charge a completely exhausted battery in a very short period of time (15 minutes). 1.2. Proiect Justification < Existing distribut i on systems may not have the capability to cater to the new needs of customers if the use of EV s becomes widespread. The electric power utility companies need to be aware of and prepared to face such load scenarios. The utility companies need to first have a reliable assessment of the possible future load scenarios with EV loads. After this, they need to have a ready design of the new distribution system

PAGE 11

4 equipment in terrns of their sizes and ratings and also in terrns of the configuration or the layout of the entire distribution system. Thus, in case of an increase in the use of EV s by the customers, the utilities will find themselves well prepared and thus continue to reliably supply quality power to their customers. 1.3. Objectives and Goals The aim of tl1is project is to provide Florida Power and Light Company with an assessment of the possible load effects w i th EV loads. The outline of the objectives of the project is given below To stu dy the nortnal load behavior for a typical distribution system feeding two to ten houses To study the load behavior when one of the houses in the above case owns an Electric Vehicle. To perfortn laboratory tests to deter1nine the behavior of the General Motor EVl charger, in terms of its current requirements and harmonic effect on the input line. Based on the above data to prepare a computer model of the distribution system feeding loads consisting of EV loads. Prepare the mathematical model for various components and the whole system. Perform basic statistical analysis of the recorded data. Use computer simulations to study various possible load scenarios. 1.4. Maior Contributions Major contributions of this dissertation work are itemized below.

PAGE 12

5 Actual field data was recorded from distribution transfor1ner sites. This data contains information about the load behavior and characteristics at these particular sites. It has been for the first time, that such data has been made available, for load scenario comparison between loads with and without the EV loads. Hrumonic analysis has been performed on the field -r ecorded data, to identify the signature that the EV load leaves on the system. This has provided with an understanding of the EV load behavior when put on the system along with other domestic l oads Actual laboratory testing of the EV charger, charging into a fully discharged battery has been carried out For the first time, GMEV 1, charger has been tested when operated in isolation from any other external disturbance or interference. The testing procedure has been developed and applied to find the current, power, voltage power factor and hartnonic characteristics of the EV charger Computer simulation of the distribution system feeding EV loads has been carried out. The simulation results have been con1pared with the actual field and laboratory results. Thus, establishing the accuracy of the simulation results Fast charging time EV chargers have been simulated and their expected behavior has been studied. Although these chargers are not commercially available, however, these s i mulation results have provided a basis for understanding the expected behavior of such chargers once such a technology becomes available. Generalized load profile on a distribution transfor1ner site has been developed. An IEEE / ANSI method for the calculation of the derating of transformer feeding non sinusoidal loads has been modified. The proposed method yields more accurate and realistic results 1.4. Research Phases The objectives of the thesis outlined in 1.3 have been achieved by dividing the whole research work into various phases. These phases are outlined as follows and later on in the section they are explained in more detail.

PAGE 13

6 These phases are: PHASE 1 DATA PHASE: This phase consists of: The recording and analysis of the field load data with and without the EV being a part of the load. The laboratory testing of the characteristics of the EV charger. PHASE II. MODELLING PHASE: This phase consists of: The physical model of the various components of the distribution system. PHASE III. COMPUI'ERSIMULATION PHASE: This phase consists of: Computer simulation of the power distribution system and the analysis of simulation results. Create and study of different load scenarios with multiple EV loads. Create and study load scenarios with EV loads of variable charging times. PHASE IV. DESIGN PHASE: This phase consists of: Calculation of expected derating of a distribution system transfor1ner based on the simulation results for: Nor1nal charging time chargers. Short charging time chargers. (Quick Chargers)

PAGE 14

CHAPTER2 LITERATURE SURVEY The issue of power quality is as o ld as the issue of electric power itself But due to changing load conditions and customer needs, the importance and s ignificance of supplying high quality power has changed with time. In the last decade, due to the increase in the use of power electronic equipment, power quality ha s attracted the attention of the utilitie s and customers alike. The utilities look to supply quality power to its customers at a price dependent on the amount of distortion that the customer is causing. The customers on the other hand want a reliable, sag free power to be delivered to them at a highly competitive pnce. The study and research in the field of electric power quality can be traced to 60's [1] and 70's [2,3] The fact though is that such studies were sparse Power quality during that era was just a limit applied to fluctuations of frequency and voltage, to the voltage unbalance voltage transient s and flickers and power cuts [1,2]. Owing to the absence of vast amount of non-linear loads, power quality and hannonics was a factor more of an academic importance rather than a serious practical concern. The harmonic s in the power systems were mostly caused due to load switching and telepl1one interference. During the 80' s with the fust power electronic revolution the use of power electronic equipment was on an increase. Most of the power electronic loads are switching loads and are highly non-linear ; these loads thus cause a 7

PAGE 15

8 lot of harmonics to be injected into the system. This prompted the academia to study the power system harmonic s in a greater detail, this included devising methods for the actual measurement of syste m harmonic s, making case studies and estimating their effect on transformers, feeders and lines [4,5,6]. In the mid 80s [7] the increase in the presence of non linear power electronic lo a ds was obvious and a further increase was anticipated This increase was owing to fir st the change in the technology e.g. faster and more efficient power electronic devices better converter topologies and converter control methodologies, second the advent of fast microproces so rs and their application to power electronic converter drive s The microprocessor controlled drives were more efficient and reliable, had the flexibility of speed control over a wide range, were ea sy to control and were compact in size. The need for technical under s tanding of the nature of electric di s turbances caused by these non-linear loads was felt Thus, in the face of changing load scenarios, there was a need for power quality education that included definitions grounding practice s, power line dynamic s, case experience and trouble shooting [8] In and after the late 80 s a sharp increase in the published literature on power quality and harmonics can be found. The number of these publications i s so large that it i s impossible to list all of them thus a selection of the se papers has been listed here Decrease in power quality that not too long ago was a mere s peculation, was now viewed as a serious threat to the normal operation of power systems [1 OJ. Research was thus directed toward s development of harmonic mitigation techniques [11] as well as de v elopment of p owe r electronic converters which would cause minimum harmonic di sto rtion Papers recommending different method s for the estimation and monitoring of power system harmonics were publi she d [12 13 14] References 15 25 include papers on the study and characterization of s pecific non linear loads and their effect on po\.\i er systems. The studied

PAGE 16

9 loads consist of fluorescent lighting loads, electric drive loads, railroad loads and some other types of non-linear loads. In addition to these papers, scores of publications dealing with the study of the effect of harmonics on various power system equipment can be found in the literature [25 30]. The last few years have seen a phenomenal increase in the use of computers by the power customers. The computer loads, wluch are highly non-linear [31 ], although exist significantly as commercial loads but the ntirober of home computers is also increasing every day. [32]. Power quality is no longer an issue that is linked with just reliability efficiency voltage sags and flickers; today power quality is understood literally as the quality of the actual current and voltage wavefo11ns. The power customer today has become very sensitive to the quality of po,;ver in terms of waveforms that the utility is providing. This sensitivity is owing to the changing needs of the customers, for instance in some cases the equipment powered from the line is highly sensitive to harmonics. Also, both commercial and domestic customers are more aware of harmonics and the impact they have on the performance of equipment and utility bill One of the very sources of harmonics in the power system is used to assess and remedy the situation. Computer modeling and simulations are being used for the study of various non-linear loads and the impact they have on power systems.[33-37]. In the midst of the talk and concern about increasing disturbance and harmonics, yet another type of non linear load has appeared on the horizon. This is the electric vehicle load. The electric vehicle load has attracted the attention of the researchers and the utilities alike.

PAGE 17

CHAPTER3 POWER QUALITY AND HARMONIC DISTORTION Power quality, in simple words can be defmed as the amount of distortion in the system current and/or voltage waveforms. The power systems and equipment are designed for pure sinusoidal current and voltage waveforms Nowadays, widespread use of computers, florescent lights, Adjustable Speed Drives and many other power electronic equipment by customers cause a considerable amount of distortion of the voltage and current wavefor1ns. The distortion is caused due to the inherently non-linear characteristics of these types of loads. In the presence of this increasing level of non-linear loads, the utility companies need to continue catering to the customer needs, which is supplying quality pov~1er with a high degree of reliability. This could necessitate new and stricter distortion guidelines for power customers, as well as a variable rate structure being developed. Also, utility engineers need to deal with the issue by analyzing the distortion and plan to control it by employing methods to either minimize the distortion by using some mitigation techniques or derate the existing system equipment to better cope with the steadily increasing non linear load scenario. Electrical energy consumption in the United States in the year 1992 \\ : as of the order of 600 10 9 kWh/year (38]. With the advances in power electronic converters, computers and many other non-linear loads, more and more of the proportion of this power \.\ : ould be used by such non-linear, distortion causing loads. 10

PAGE 18

l l 3 .1. Hannonics The non-linear characteristics of the loads cause the distortion of the current and / or voltage waveforms in a power system. These current and voltage waveforms are periodic in nature but they are not pure si11usoids. The distorted waveforrns can be expressed in the form of a Fourier series. Using the Fourier series expansion, a distorted periodic wave can be expressed as the sum of an infinite number of sine waves, whose frequencies are an integral multiple of the fundamental frequency Where: The distorted periodic waveform f (t) can be expressed as: f(t) = a 0 + L :_ 1 [a,, cos(ncvt) + b,, sin(nmt)] a 0 is the de offset of the waveform a,, is the amplitude of the cos term b,, is the amplitude of the sin term n is the order of the harmonic To study the effect of tl1e distorted voltage and current wavefor1ns on the power system, the principle of superposition is used. This means that the effect of each constituent (harmonic) of the distorted wave is studied separately and then the overall effect is found by adding the effect of individual harmonics. The fundamental constituent of the distorted waveform has a frequency of 60 Hz, the second har1nonic has a frequency of 120 Hz, the third 180 Hz, and so forth. The equipment used in the povver system is designed for a supply frequency of 60 Hz i e. the fundamental frequency. The other constituent wavefor1ns (harmonics) certainly have an undesirable impact on the equipment and the overall power system and its perfonnance.

PAGE 19

12 3 .2. Estimation of Harmonic Distortion Before assessing the qualitative and quantitative effect of har1r1onics on the p owe r syste m, so me kind of estimate or quantitative expression of harmonic distortion needs to be defined. Three methods of estimating harmonic load content are practiced [ 41] The Crest Factor and The Total Harmonic Distortion (THD) are the mo s t common method s. The third method is the K-Factor method. 3 .2 1. Crest Factor The Crest Factor gives a very simple estimate of the hannonic content of the waveform It i s defin e d as the ratio of the peak of the wave to its Root Mean Square (RMS) value. A perfect sine wave by definition will thus have a crest factor of 1.414 C t F t Peak Value of th e Waveform res ac or = _________ ,;;___ RMS Value of the Waveform If the Crest Factor of the wavefor1n (voltage or c11rrent) is other than 1.414 it indicates some kind of distortion in the wave. This technique has a limitati o n it doe s not provide inf or1nation about the constituent harmonic frequencies. Knowledge of the har1nonic frequencie s and their re s pective amplitude s is important to study the effect of wave di s tortion on the sys tem

PAGE 20

3.2.2. Percentage Total Hannonic Distortion ('I'HD) Total Harmonic Distortion (THD) is defined as the ratio of the sum of the RMS values of the harmonics to the RMS value of the fundamental '' co 2 ~h= 2 (f,,) THD = -=---I, 13 The above expression shows that THD provides with a measure of the distortion of a waveform. Like the Crest Factor, the THD is also limited in that it does not provide any information about the frequencies and amplitudes of individual harmonics. While computing the THD for a wave, the relative phases of the various harmonics are also not considered. Nevertheless both Crest Factor and the THD are very useful to establish the presence or absence of waveforrn distortion. 3.2.3. The K-factor The K-factor while providing an estimate of distortion in a wavefo11n takes into account the harmonic frequencies and amplitudes. This fact makes it the most accurate available method of estimation of non-linear load harmonics in a system. This is the method used for the calculation of derating of dry-type power distribution transformers [42]. Reference 4 proposes an expansion of the K-factor method. The author (41] proposes to use the K-factor as a weighted sum. Each type of load is considered separately and then all the Factors are combined to create a composite harmonic fingerprint.

PAGE 21

14 The K-factor is defined as: aJ K = L 11, (pi,) 2 h 2 h=l 3 .3. Transformer Derating Power transformers and distribution tran sfor 111er s are designed for s inu soidal ct1rrents and voltages. Name plate ratings of transformers hold for purely sinusoidal supply and loads If the current gets non-sinusoidal, for in stance becau se of the use of switching power supp lies so lid state inverter s, converters, Adjustable Speed Drives etc., the indicated ratings of tran sformers no longer hold. A high harmonic content in the output current wavef or1r1 effects the transfonner ratings. Load har1nonic currents add a disproportionate amount of additional transformer heating particularly in the for1n of winding eddy current losse s and stray losses which are supply frequency dependent. Several methods are currently used in the industry to assess the derating of transformers for a har1nonic application. The most widely used method is IEEE / ANSI C57. l 10. The IEEE / ANSI method was developed for power transfortners and does not specifically address other applications such as di str ibution and dry type transfor1ners. On the ot h er hand, the K factor method can be applied to only dry type tran sf ortners. In general the application of the IEEE / ANSI method provides conservative re sults when applied to a distribution transf orme r. Thi s standard outlines the following method for calculating the transformer harmonic capability.

PAGE 22

15 3.3.1. Transformer Capability Equivalent Calculation. In this method a per unit value of non-sinusoidal cunent is calculated which will result in load losse s in the transforrner just eqt1al to the highest lo ss region for 60 Hz rated operation. The winding eddy cu1Tent loss under rated conditions at the point of maximum loss density is 15 % of the local I 2 R lo sses Taking the fundamental component of current as 1 per unit and taking the measured per unit values of the current harmonic s up to the 11 th harmonic the derating of a transformer is calculated. The maximum l oss den s ity of the local I 2 R los ses i s a ss umed as 15 %. Solving Equation 8 in ANSI/IEEE C 57 .11 0 for the permissible current. Following the procedure in the standard: I max = l 15 l + X O. 15 U s ing the values shown in Table. I. and obtaining the amplitude of each harmonic from Figure 4 the maximum permi ss ible non-sinu so idal load current is given b y: I 1 15 O O mnx = 113.089 = .S l l +--x. 15 4.75 Hence the ratio of transformer rating with and without the harmonic load s = .50 l

PAGE 23

16 Tab l e. I Harmonic Content of the Sa mpl e Waveform HARMONIC PER UNIT CURRENT (F H) ( H ) FH 2 FH2H2 1 1 1 1 2 .16 .0256 .102 3 1. 3 1 69 15 .2 1 4 .267 0712 1.139 5 1 1 25 6 101 .0 10 0.36 7 .68 .462 22.638 9 .54 291 23.57 1 1 .447 .199 24.07 Total 4.75 113.089

PAGE 24

17 3 3.2. Effect of Current Harmonics on the Neutral Current. The previous sect i ons of this chapter discuss the effect of l1armonics on the distribution transformer. However there is some more concern for harmonic loads on multigrounded wye systems [43 44]. The current harmonics in a Y Y distribution transformer cat 1 se a high neutral current and an elevated emf of the neutral. In an ideal case when the loads on the transfonner are balanced, the neutral current which is the vector sum of the three phase currents is zero. If: I A= Asincvt l a = A sin(cv t + 240 ) l e = Asin(cv t + 120 ) wh e re A is the amplitude I A, 1 8 and l e are the currents in phase A, Band C respectively The neutral current is the vector sum of the three phase currents and is given by: IN = A sin cv t + A sin( cv t + 240 ) + A sin( cv t + 12 0 ) ... IN = 0 Now if the load current is rich in the third harmonic, the third harmonic component with amplitude T can be written as: I A= Tsin(3cv) l a = T sin(3cv t + 240 ) l e= Tsin(3cvt+l20 )

PAGE 25

And as before IN= T sin(3m t) + T sin(3m t + 240 ) + T sin ( 3 m t + 120 ) :. IN = 3T sin ( 3m t) 18 The net neutral current becau s e of all the other tripplen hrun1onics will yield similar results. As reported by Mansoor et al [44] in a power system, even with load diversity there is very little ca.11cellation of the third hannonic. The s e high neutral current s cause the following problem [ 45]: 1. The heavy neutral currents cause the over loading of feeder transformers . They ca.11 cau s e further distortion of the voltage because they push the transfonr1er into saturation. The saturated transformer then becomes a source of harmonics. 2. W11en the neutral conductor carries harmonic currents additional heat is generated and the current carrying capacity of the feeder is reduced. 3. The neutral conductor for tra.11sf armers is generally sized same as the phase conductor, but due to the tripplen harmonic currents the neutral conductor current ca.11 exceed its rated value thus overloading the neutral that can lead to equipment failure. 4. The balanced tripplen currents flow as circulating currents in the transforrner delta p1imary winding. As a result more current flows in the \,\,1t1ding of the transfo11ner tha.11 is detected by the transf onner primary circuit over current protection device. This can result in the overloading of the transformer. These high neutral hannonic currents cause a voltage differential between the neutral and the ground since the ground wire at hannonic frequencies ca.11 cause significa.11t neutral conductor voltage drop. This is known as Common Mode Noise Voltage.

PAGE 26

CHAPTER4 FIELD DATA RECORDING REPRESENTATION AND OBSERVATION. 4 .1. Introduction The effect of the EV loads on the power distribution system would depend on a number of factors. Once the EV s ~ecome popular one of the factors that is going to be of significance, is the usage pattern of the EV by its users. The usage pattern would obviously vary from individual to individual it would depend on the lifestyles and specific needs of the customer, e.g. nature of his job, weather conditions, age of the user, so on and so forth. To generalize the usage pattern and come up with an accurate generalized model is surely a nontrivial problem Based on field-recorded data and knowledge of the EV load behavior, a typical load profile of a distribution node feeding EV loads can be constructed Even though this load model won't be accurate but it will contain enough information to enable the utilities to design the distribution transforrr1ers, so that they are able to supply power reliably and econo1nically. The usage pattern of users on the same distribution transformer shall determine the load diversity factor, the expected peak load on the transfor1ner and the overall derating of the trru1sformer The objective of the field data monitoring tl1us is threefold First to get the load profile during two separate periods viz. when none of the customers are using an EV and when one of the customer is using an EV. Second to try and identify the signature or effect of the EV load when it is connected to the system. This effect is in terrns of 19

PAGE 27

20 increase in power demand and possible decrease in power quality. Lastly, this field recorded data can form the basis for configuring system computer model used for carrying out computer simulation study. The objectives of field data monitoring l1ave been achieved by monitoring the field data (system voltages, currents, power factor etc.) first for a period of time (two weeks) when the EV load is not a part of the loads on the distribution transformer This is followed by a two week period then when one of the customers is using the EV. The data obtained during the above two l oad scenarios has been analyzed. The comparison of these l oad scenarios ha s been used to estab l ish a generalized l oad scenario when one of the customer on a distribut i on node is using an EV. 4 .2. General Considerations In Data Acquisition And Processing Appropriate techniques for the acquisition and processing of random data are heavily dependent upon the physical phenomenon represented by the data and the desired engineering goals of the processing In broad terms, however, the required operations may be divided into five primary categor i es as fol l ows: a. Data collection b. Data recording (inc l uding transmission) c. Data preparation d. Data qualification e. Data analysis Each of these categories involves a number of sequent i al steps. The purpose of this sect i on is to summarize basic considerations associated with each of these key steps.

PAGE 28

21 4.2.1. Data Collection Tl1e primary element in data collection is the instrumentation transducer. In general terms, a transducer is any device which translates power from one form to another. In an engineering context, this usually means the translation of some measure of a physical phenomenon of interest into an analog signal with a calibrated relationship between the input and output quantities. This translation may involve up to three basic operations; (a) a mechanical conversion of the physical quantity of interest into an inter1nediate mechanical quantity, (b) a pick off step which converts the intermediate mechanical quantity into an intermediate electrical quantity, and (c) an electrical conversion into a fmal electrical quantity, usually voltage. Some transducers may combine a11y two or all three of the operations, depending upon the physical quantity being measured and the specific nature of the transducer. For example, a ther111ocouple, which is widely used as a temperature transducer, converts a difference in temperature directly to a difference in voltage without intermediate steps. On the other hand, a resistance therrr1ometer, another commonly used temperature transducer, first converts temperature to a change in resistance, and then converts the resistance change to an electrical voltage change. Ideally, the above operations would be accomplished without distortion or modification of the time history of the physical quantity being measured. In other words, if the input time history is x(t) and the output time history is y(t), the perfect transducer would provide an analog output, y(t) = c x(t), where c is a simple calibration constant. Unfortunately, this ideal situation is difficult to achieve in practice. Gain and phase modifications as well as distortion producing non-linearities are often inherent in the

PAGE 29

22 transducer operations. This fact makes the transducer a potential source of error in any data acqt1isition and processing program. i 4.2.2. Data Recording For some applications, it is possible to perform all desired data processing directly on the transducer signals in real time. For most applications, however, this is not practical and some form of storage (and perhaps remote transmission) of the transducer signals will be required. The most desirable and convenient type of data storage system is the magnetic recorder or a computer hard drive memory. Although other typ es of recorders could be used, the magnetic recorder has the advantage of being able to store large quantities of data, and to reproduce them in electrical forrn The most desirable way to transmit data signals is through electrical lines or telephone lines. There are obvious situations where this is not feasible; for example, retrieving data fron1 a spacecraft in earth orbit. For such cases, radio transmission (telemetry) of the transducer signa l s is usually required. 4 .2.3. Data Preparation The next key pha se in data acquisition and processing is the preparation of raw data for detailed analysis The raw data are usually supplied from the recorder in the for1n of voltage time histories ( or as direct analog signals from the transducers, if the data are being analyzed on-line). A number of operations are needed at this point to make the voltage time histories suitable for detailed analysis. The first of these operations is generally classified as data editing.

PAGE 30

23 Data editing refers to those pre-analysis procedures which are designed to detect and eliminate spurious and/or degraded data signals wl1ich might have resulted from acquisition and recording problems such as excessive noise, signal dropout, loss of signal due to transducer malfunctions. Editing can often be accomplished through visual inspection of the data time history signals by a talented analyst. In more elaborate data acquisition and processing systems a specific instrument for quick-look evalttation might be employed. Real-time spectrum analyzers are popt1lar for this application. It should be noted that the data editing step is more critical for the case of digital processing than for analog processing. This is true because 011ce the data have been converted to a digital format, it is often difficult to detect even the most obvious errors in the original signal. For the case of analog data processing, data preparation beyond editing usually includes only conversion into engineering units ( calibration). 4.2.4. Data Qualification The correct procedures for analyzing random data, as well as interpreting the analyzed results, are strongly influenced by certain basic characteristics which may or may not be exhibited by the data. The three most important of these basic characteristics are the stationarity of the data, the presence of periodicities in the data, and the normality of the data. Stationarity is of concern because the analysis procedures required for non stationary data are generally more complicated than those which are appropriate for stationary data. Periodicities in the data should at least be identified to avoid erroneous interpretations of later results. The validity of an assumption that the data ( excluding periodicities) have a Gaussian probability density function should be investigated since

PAGE 31

24 I the normality assumption is vital to many analytical applications for random data. Qualification of sampled data in terms of these basic characteristics is indicated as a separate operation to be performed prior to detailed data analysis. In practice however, it is often accomplished as an integral part of the data analysis phase. Practical considerations and procedu res for such qualification will now be discussed 4.2. Data Analysis The procedures for analyzing the properties of random data may be divided logically into two categories: the procedure for analyzing individual sample records, and the procedure for analyzing a collection of sample records given the properties of the individual records. Applicable data analysis procedures for these two categories are now outlined MEAN AND MEAN SQUARE VALUE ANALYSIS. The first step is a mean and mean sqt1are value ( or variance) measurement. 1bis step is almost universally performed for one or more of three sound reasons. First, since the mean and mean square values are the basic measures of central tendency and dispersion, their calculation is generally required for even the most rudimentary applications. Second, the calculation of short time averaged mean and mean square value estimates provides a basis for evaluating the stationarity of the data. Third, mean and mean square value estimates can be extracted from other descriptive properties (probability density plots, correlograms, and/or power spectra) which might be measured later The comparison of directly measured mean and mean square values estimates to the corresponding estimates extracted from other

PAGE 32

25 ... analyses provides an excellent method for checking the data analysis equipment or computer programs for correct operation. POWER SPECTRAL DENSITY ANALYSIS. Perl1aps the most important single descriptive characteristic of stationary random data is the power spectral density function which defines the frequency composition of the data. For constant parameter linear physical systems, the output power spectrum is equal to the inpt1t power spectrum multiplied by the square of the gain factor of the system. Thus power spectra measurements can yield infor111ation concerning the dynamic characteri s tics of the system. To be more general, the mean square value of the data in any frequency range of concern is determined by the area under the power spectrum bounded by the limits of that frequency range. Obviously, the measurement of power spectra data, will be valuable for many analysis objectives. Secondary applications include its use for the detection of periodicities and as an intermediate step in the calculation of autocorrelation functions. PROBABILITY DENSITY ANALYSIS. The last fundamental analysis included in the procedure is probability density analysis. Probability density analysis is often omitted from a data analysi s procedure because of the tendency to assume that all random phenomena are normally distributed. In some cases however, random data may deviate substantially from the Gaussian form. If such deviations are detected by a test for normality, then the probability density function of the data must be measured to establish the actual probabilistic characteristics of the data. Furthermore a probability density function estimate is sometimes used as a basis for a normality test.

PAGE 33

26 4.3. Field Data Monitoring 4 .3 .1 Data Acquisition Methodology. The data was recorded in the form of one cycle snap shots of voltage, current and power factor. These waveforrns were acquired once every one minute interval. The voltage and current waveforms at the distribution transformer output terrninals were recorded and stored temporarily in the local hard drive of the data acquisition system. The stored data was then down loaded to a per so nal computer over a telephone line using modems. The distribution nodes under study were feeding typical domestic loads at 110 volts and EV load s at 220 volts. At each of these locations the distribution transforrner has typical domestic loads of 2 to 10 homes. At each location one EV was provided to one of the customers for a period of two weeks. The customer was asked to put the EV to normal use The data at each location was recorded during time intervals that can be classified into three states. The first state i s the period (about two weeks) when there was no EV load. The second state is the period ( about 2 weeks) when one of the houses on the distribution transforrner was using an EV. The third state was again when there was no EV load. Thus the field data collected represents the load conditions before; during and after the EV load was connected to the system. Figure 4.1 shows the data acquisition system layout. 4.3 .2. Data Monitoring Instrument. The field data was recorded using several BMI models 8020 Plus PQNode data acquisition instruments. The BMI model 8020 Plus PQNode monitors power quality phenomena in electric power systems. It combines the capabilities of instruments that

PAGE 34

27 monitor continuous quantities, (such as voltage level, load level, power factor, and harmonic distortion), with the capabilities of instruments that record disturbances (impulses, wave shape faults, swells and sags, outages, cold load pickups). Monitoring continuous quantities requires that the monitored signals be sampled at a regular basis, wh i le monitoring disturbances requires that the signals be continuously sampled, and recorded only if the signals exceed specified values. The BMI 8020 Plus PQNode continuously samples up to four voltages and four current three phases and neutral. The instrument periodically records wavefor1ns, RMS levels, frequency and temperature so those trends can be ca l culated. The PQNode TM is designed to interface with a personal computer via telephone line, and is equipped with an inten1al modem for this purpose. 4 .3 .3. Data Management The data acquired by the BMI data acquisition system was stored by the BMI unit in the hard drive of the device in the for1n of ASCII files. Each set of data was stored as a 4 K file. The actual information is in only I kilobyte of the 4 K space. After downloading the data files to the personal computer in the Power Quality and Power Electronics Lab, at The University of Florida, the data needed to be stored safely with adequate backups and at the same time consuming optimum space. The data also needs to be easily restorable whenever needed for analysis purposes Using data acquisition and data storing software, tl1e data was archived in the computer hard drive. The archived data occupies only one fourth of the space since the files are crunched together. The data was also backed up on tape drives connected to the computer. In this way it vvas made certain that the precious data is not l ost in the event of a hard drive failure. For the purpose of

PAGE 35

28 analysis, only the required amount of data is restored a11d after it is processed it is deleted from the hard drive and a fresh batch of data to be analyzed and processed is restored I from the archives. For each batch of data the THD trend is plotted. Also plotted are the voltage and current trends for all the harmonics from the second through to the fiftieth harmonic component. In addition to the harmonic trends, steady state trends are also plotted; these are the real power, apparent po1vver, RMS current and power factor trends for the two phases between which the EV charger is connected. Each of the trend wavefo11ns is displayed one by one . After exan1ining the harmonic and steady state trends for the selected days, typical real time steady state voltage and current waveforms at certain instru1ces of time were also plotted It may be mentioned here that the BMI instrument acquires the distribution node data in the form of real time waveforms of voltages and currents. Using these waveforms which are acquired during every fixed period of time, the computer software uses them to create various har111onic and steady state trends.

PAGE 36

L. L (/) E ... I.. 0 > VI 0 C N tG I.. ,IC en 0 .... 0 +' ;:, > J:J 0 I.. ., N Cl> ,C L Data Acquisition Unit Mair1 Unit u ~ Hard Telephone Line Drive ......... Aod~n ....... ---------~ -I I Personal C0n1puter Domestic Do m estic oo m es~ l c Loidi l 08C!! L~ads Electric (/J ... 0 > Vehicle 0 Load o~mest 1 c Domest i c "d' Do m~st i e Load! L o a d s LO~ds N .. F i g ur e 4 1 D ata Acq u isition System Layout 29

PAGE 37

CHAPTERS DATA REPRESENTATION AND ANALYSIS 5 .1. Introduction This section deals with the graphical representation of the recorded data A visual inspection of the data trends is enough to show that the data represents stationary random process. It is assumed that all random phenomena are normally distributed. [data book] In some cases, however, random data may deviate substantially from Gausian forrn. Looking at the data trends included in this chapter, the properties of the recorded data can be described at any instant of time by computing the average vales over a collection of sample functions that defme the random process The stationary random data represented here falls in the category of ergodic random data The properties of ergodic random process can be determined by performing time averages over a single sample function. Fortunately in practice, random data representing stationary physical phenomena are generally ergodic. It is for this reason that the properties of stationary random phenomena can be measured properly, in most cases, from a single observed time history record. The recorded data is presented in the forrn of various trends and waveforms. Each data trend for a particular site with and without the EV load are presented together to facilitate comparison between the two load scenarios namely when the EV load is the part of the load and when the EV is not a part of the load. The data recorded using the BMI data acquisition 30

PAGE 38

31 units was synthesized using the software PASS. Broadly speaking two types of data trends were created and displayed graphically, viz. the steady state trends and the harmonic trends. For the purpose of studying tl1e impact of Electric Vehicle loads on the distribution system, only the following trends have been selected: 1. The Steady State real power trend. 2. The Load Current Total Harrr1onic Distortion (THD) trend. 3. The various harmonic trends of the load current. Although there are many other steady state and har111onic trends that were displayed using the PASS software, only the above mentioned three trends have been used for data analy s is for the following reasons: The EV charger load is a real power load which consumes a peak power of 7 kilowatts (based on GMl charger). The 7-kilowatt power is consumed when a fully discharged EVI battery is connected to the charger. The power decreases as the charge on the battery increases. So one of the impressions that the EV charger load causes in the system is an increa se in real power. Thus the real power trends for the distribution node with and without the EV load were created. This made it possible to compare the load scenario with and with out the EV load. From the laboratory testing of the EV charger, it was found that th e EV load causes a distortion of the input line current. The General Motors EVl is known to cause very little input voltage di s tortion [2]. VIZ.: The data represented in thi s chapter can be broadly classified into two categories I. When the distribution transfor1r1er is feeding normal domestic loads. II. When the distribution transfortner is feeding no1111al domestic loads that also include Electric Vehicle Load.

PAGE 39

32 Representing the data in this manner has facilitated the comparison of the load scenarios with and without the EV load. 5.2. Data Trends and Characteristics Figure 5.1 through 5.7 are sample data trends created from the actual monitored field data. Each of these figures consist of two parts, viz. part a and part b. Parts a are the data trends without the EV loads, parts b are the data trends when the EV load is a part of the loads on the transfonners. Following is the summary of highlights of each of the figures 5.1 to 5 7. Figure 5.1: This figure represents the Total Harmonic Distortion trends of the line voltage. The voltage THDs for both cases are very close. Maximum voltage THD variation between ''a'' and ' b '' is 2%. The presence of EV loads does not have a significant effect on the input voltage THD. This fact is verified by the laboratory testing of the EV charger Figure 5.2: Unlike the voltage, the curre11t THD trends shown in figure 5.2 are sensitive to EV loads. Comparing 5.2 a and 5.2 b it can be seen that the current TI-ID trend for two cases is not the same. When EV was a part of the domestic load, we see in 5 .2 b the current THD between the hours of 2:15 and 4:40 is very high. The current TI-ID during this time varies between 80% to 40%. For the remaining part of the day, the current THD for both cases is close to each other. Figure 5 3: These figures show the real power trends. The real power trend when EV is a part of the load is not much different from the case, when EV is not a part of the load Tilis could be due to the fact that the customers do not fully discharge the batteries before

PAGE 40

33 recharging them. Thus each time the battery charger load is put on the system, it draws only a fraction of its maximum rated power of 7 kilowatts. Figure 5.4: These show the power factor trends. The EV load does not effect the power factor Thu s we see, in both parts of this figure the power factor trend stays close to .8. Figure 5.5: These figures represent the third ham1onic current trends. As in the case of the overall current THD trends, the third harmonic current trend between the hours of 2: 15 and 4:00, show a high, for the case when EV is part of the load. The third harmonic component touches a high of 1 5 amperes during this period. This again indicates a possible presence of EV loads during these hours. Figure 5 6 and 5.7: These are the fifth and the seventh harmonic current trends. The conclusions drawn from these trends are similar to those drawn from figure 5.5. 5.3. Data Observation Observing the data trends represented in figures 5.1 through 5.7, the following conclusions can be drawn: I. The EV load does not have a distorting effect on the input line voltage Thus the input line voltage THD is not effected by the presence of EV loads. II. The EV loads cause a considerable distortion of the line current The line current THD is particularly high when the amplitude of the fundamental of the line current is small. Thus the current distortion is ma.xi.mum close to the end of the charging cycle when the charger draws a small current.

PAGE 41

34 III The EV users generally would not be expected to fully discharge the batteries before recharging them. Thus the EV loads may not have a significant effect on the expected peak load on the distribution transforrners. IV. The EV battery charger operates at a high power factor. Thus the reactive power drawn by the charger is low. V. The EV charger injects a significant amount of lower odd current harmonics into the power line. 5. 4. Load Characteristics This section contains the computed maximum and the average daily power output of the distribution transfonner. The daily average power is the average of the load over a period of two weeks The load average has been found, first for the period (2weeks) when EV was not a part of the loads and next for the period of 2 weeks, when EV was a part of the load. Also, for these two two-week periods the maximum load profile has been computed. The figure X.8 a, has been derived by computing the load average for a 2-week period when EV was not a part of the loads. The figure 5.9b is the 2-week load average when EV was a part of the loads. Comparing, 5.8a and 5.8b, we see tl1at there is an increase in the average power consumed, when the loads contain the EV load. This increase is particularly visible, during the early morning hours, this reinforces the previous conclusion about the times at which the EV load was connected to the system. Figure 5.9 a and b, show the maximum power recorded during the two two-week periods. Plot in figure 5.9a is the maxima of the daily load curves of the two weeks period when EV was not a part of the loads. Plot 5.9b is the maxima of the daily load curves of the

PAGE 42

35 two-week period when EV was a part of the loads. It must be noted, the EV loads although do not effect the peak loads, nevertl1eless being an additional load it does affect the average daily power consumed. Figures 5 .10 a and 5 .10 b, show the snap shots of the recorded line current waveforms. The line current snap shot of figure 5 .1 Oa was recorded vvhen EV was not expected to be a part of the loads and 5 .11 a was recorded, when EV was expected to be a part of the loads. From the Fast Fourier Transfortns (FFTs) of wavefo1ms in figure 5.1 Oa, shown in figure 5 .1 Ob, we see that the line current does not have a high content of the lower odd harmonics (viz 3 rd 5 th and the 7 th ). This is not the case with figure 5.11 b where these harmonics are significant. Based on the preceding discussion, the load profile of Figure 5.9, qualifies for a generalized load profile of loads on a distribution transforrner, for that particular site.

PAGE 43

00 : 01:41 01 : 10 :38 1-z"j 02 : 19 : 36 I 03 : 28 : 34 00 E; I 04 : 37 : 32 (I) VI I 05 : 46 : 30 ...... < I 06:55 : 28 0 ...... f:I) I 08 : 10 :3 4 00 (I) I 09:19 : 32 0 ,...,. 10 : 28 : 29 f:I) --i ::c ,,--... 3 11 : 37 : 27 '--" CD 0 12 : 46 :25 -+. 0 0 2. 13 : 55 : 23 D.) () '< t::J 15 :0 4 : 21 (/) ,...,. I 16 : 13 : 18 0 I 17 : 22 : 16 0 ::J I 18 : 31 : 14 ..., (t) I 19 : 40 : 12 ::J 0.. en 20 : 49 : 10 21 : 58 : 08 23 :0 7 : 06 P e r c ent THO 0 N W 01 I i I r ... l : I~ ' I L ; ' I I i l ' i ( 1 I ' ' I ~ ,A < 0 ...... Q) cc CD --i I 0 ::::r m < r 0 ru 0. en Oo_ a 0 v. 7 "{9 0 -.,.=, . 7.> ~ 0 'i,T. ~v! 7,;r (5' o,;r . {:> 9'. a a. (5'. 7 0 O.g 9: -i (9. {:> "0.> ,,--... e 3 a v. CD 7J> .s 0 -7 ._::,. -+. 7~~-ou> 0 D.) 'i,T 7 c:?. o 0 '< 76' {s, 'iT :>. 7.>: Se S: 7. ~& .9. 'iT 7 ._::,q-? .S7 ._::, -..;. 7 . 'iT.9 'iT ._::, 7 ~S'Se i,j>_ .S i,j> P e r ce nt THO 0 N . O Ol ~ OlN OlW t i ' I i ' ; 1 I I 1\1 I I < 0 ...... D.) (Q CD --i I 0 --i ""'"I CD :J 0. z 0 m < r 0 D.) 0. en w O'I

PAGE 44

0 I I+-' C (1) (.) I,... (1) 0.. 0 I I+-' C (1) (.) I,... (1) 0.. 37 Current THO T rend : No E V Load 20 15 10 5 0 .... ,I-. ., ' -I Ill 0) r--lO M 0 co (0 M 0 co (0 ......... ......... 0 0 ....0 0 0 0 . . . ' ' . . . 0 0) co r--(0 lO '<:t' M N ..... 0 0) co 0 0 ..N M lO 0 ..... N M M ' . . . ' .. . 0 ..... N M 'V lO (0 co 0) 0 ..... N M 0 0 0 0 0 0 0 0 0 ..... ..... ..... ..C) r--lO lO . (0 lO lO 0 . . (0 ....... lO M ..lO lO lO .. .. M N ,N M . . r--co C) ..... ....0) l{) M l{) l{) . ,(0 lO 'V lO 0 . . 0 ..M N N N Time of Da y (a) Current THO Trend : E V L oad 100 80 60 40 20 0 ._ Time of Day (b) Figure 5.2. Line Current Total Harmonic Distortion Trends

PAGE 45

R eal P o w e r i n KW Oo. I O 0 7. ..... 7. 'if; 0 N O> O> 0 I ; 0 : I r 1 I I I I I <9 0 If~ .I { I .. I I I I I ~d c5' ov. av i i~ I I I I "'7. o. ~d t ~c:> 0 :::0 I 1-Tj CD 7. O> (JQ 0. '-?~ ~d c5' -0 .., (1) 0. 7. 0 Vt ~ v>a l :E v 7. '-?d CD u-) f .., q~ 7 -i 7 if~ -i .., (1) ~c:> CD --.. p.> CT 3 ::J ....._, CD 7\T. ,?6' a. 0 7. 0 7 '-?~ m 'li;'d V < (1) 0 .., 7. w ?IS. -~ r -3 '< 0 v .., 'J>. w (1) 7~~ {9 a. ::s 0.. ISI I (/) ;a. > O> I i I I I I I 7.9_ ?IS I 7. ,J I I I I I .9 -qd;a t'l I I I I I -? 7. 7. 7 0 v t i I I I I I .? 'J>. ~~ Oa l l IT I I I I I ifa c5' W I ' ' --J I I I I I -i 3 CD --.. p.> -.__,; 0 0 w '< I I I I I R eal P owe r i n KW O_.NvJ~<.nO>__, 00 : 00 : 19 01 : 09 : 17 02 : 18 : 15 03 : 27 : 13 04:36 : 10 05 : 45 : 08 1 I I I 06 : 54 : 06 08 : 03 : 13 I 09 : 12 : 10 10 : 21 : 08 11 : 30 : 06 12 : 39 : 04 I 13 : 48 : 01 14 : 56 : 59 16 : 05 : 57 17 : 14 : 55 18 : 23 : 53 19 : 32 : 51 20 : 41 : 49 21 : 56:55 23 : 05 : 53 :::0 CD w -0 0 :E CD .., -i .., CD ::J a. z 0 m < r 0 w a. w 00

PAGE 46

L... 0 ...., u ro LL L... (1.) 0 a. L... 0 u ro LL L... (1.) 0 a. Power Factor Trend : No EV 1 2 ------------------------------. 1 0 8 0 6 0 4 0 2 0 ~..... .. .. ,,.... Time of Day (a) Power Factor Trend : EV Load 1.2 ----------------------------1 lwwrr::-~,_f_J~( :::J 0. 8 ~ 1------!.!. -,0 6 -1-----------------=------------1 0 4 -1--------------'-------------=-t 0 2 --------------------------! 0 ............. ,, ........ ..J ., ,, .. ...._._ ' ..... ... .... .I .. --T im e of Day (b) Figure 5.4. Power Factor Trends 39

PAGE 47

'Tj (JQ C ""1 (I) Vl I.JI --, ;:r ""1 0.. ::r: 0 ::;J n n C (I) ::;J .-+ --, ""1 (I) ::, 0.. en ,,.-... IT Current in Amps 0 -" l'v W 01 00 : 01 : 41 ,, 01 : 10 : 38 02 : 19 : 36 iltL 03 : 28 : 34 04:37 : 32 r l 05 : 46 : 30 06 : 55 : 28 [ I 08 : 10 : 34 09: 19 : 32 1 :, 10 28 : 29 3 (I) 0 11 : 37 : 27 t I O 1 2 : 46 : 25 Q) '< 13 : 55 : 23 15 : 04 : 21 16 : 13:18 17 : 22 : 16 18 : 31 : 14 ... 19 : 40 : 12 20 : 49 : 10 21 : 58 : 08 23 : 07 : 06 ---i ....) -, 0.. I D.> -, 3 0 ::J (') () C -, ro ::J --i -, CD :J 0.. m < r 0 D.> 0.. ,,.-... Current in Amps 00 : 00 : 19 "'w 01 0) 01 : 06 : 17 : L I. j j 02 : 12 : 15 ;:1 03 : 18 : 13 : 1=1 04 : 24 : 11 05 : 30 : 09 : 06 : 36 : 07 07 : 42 : 13 08 : 48 : 11 I ; --i 09 : 54:09 : t l I I I 3 11 :00:07 CD 0 12 : 06 : 05 I 0 13 : 12 : 03 D.> '< 14 : 18 : 01 15 : 23 : 59 16 : 29 : 56 17 : 35 : 54 18 : 41 : 52 19 : 47 : 50 20 : 59 : 57 22 : 05 : 55 23 : 11 : 53 --i ::::r-, 0.. I D.> -, 3 0 :J (') () C -, -, CD :J --i ro :J 0.. . z 0 m < r 0 D.> 0.. .ii. 0

PAGE 48

41 Fifth Harm onic Current Trend : No EV Load (/) a. 2 5 E 2
PAGE 49

Seventh Harm onic Trend : No EV Load Cl) 0. 1. 2 E 1 <( 0 8 C 0.6 --------+-' 0 4 C Q) 0 2 L. -------------L0 :::, ,. h ,.. .. ; .... . u Time of Day (a) ,-------------------------------------~ Cl) 0. E <( C: Seventh Harmonic Current Trend : EV Load 1 2 1 1---,:-r-----------------------1 0.8 0.6 ~ -0.4 ---0 2 0 ---.. ...,. .. ......, .. ~.... --Time of Day (b) Figure 5. 7. Seventh Harmonic Current Trend .. 42

PAGE 50

43 Average Power : N o EV Load 10 8 C 6 ..... 4 Q) 2 0 a.. 0 ...... _,, W Time of Day (a) Average Power : EV Load 14 12 10 C: 8 ..... 6 Q) 4 2 0 a.. 0 '" ... ,,.., u ., ~ . . Tim e of Day (b) Figure 5 8 Average Power Computed on the Basis of Two Week Daily Load

PAGE 51

44 Max Power Demand : No EV Loads 14 12 1 0 8 C ... Q.) 6 3 0 4 a.. 2 0 T l me of Day (a) Max Power Demand : EV Load 25 20 C 15 L.. (1) 10 3: 0 5 Q. 0 Time of Day (b) Figure 5.9. Maximum Power Demand Computed over Two Week Period

PAGE 52

CUrT-entWa\ie rm No EV Loa 30 ~-------,-------------,----------.-----,,.,----,-, v Max:ll3389 20 I Mii: 21.7285 A'-9: 12 JJ72 -l-----------,----,..~,,,.....-------+-------Abs: ll3389 C (1/ 10 I -I {.. {.. I :, 0 u -10 -, I RMS: 13.fQ57 CF : 1.61574 FF : l.l.2 r I ______ __; -30 '---------l------------------------OB C (1/ t.. 06 {.. :, U OA 02 00 120 0 --I I I s 10 Tile s (a) FF I of-ihe i.oad CUrrentWa\rm I I 1111 .11 1,I .1 I 1.11. 11, .11 720 1320 1920 z (b) 15 20 NO EV Load I ~ :ro : 32 6165 1HD : 33.Jl/3 1D : 0 ; 10 : 3.3J773 RMS'l : l.D/885 RMsi: 0 RMSfn: 32.6343 RMS : 32.6343 I A9JM : 39.5032 I 1F : ll).9il6 I T : 4273.32 .1l111 ... ll11l1 I 2520 3120 Figure 5. 1 0. a) Line Cu rrent Waveform Snap Shot with No EV Load b ) Fast Fourier Transform of th e Current Wavefor1n of ''a'' 45

PAGE 53

10 r"\ <{ '-.J 5 C Q) 0 L L ::, u -5 F~ Recorded OJrrentwitl EV L.oads PASS DEturt:m:e Max: 9.52148 Mn: -10.7422 1 A'v<]: 5 33581 \ I ---+'<-'---i------\-/....-----+-------tts: 10.7422 RMS: 6Ja585 CF : 1]5933 ;\1 ,\ ,I i i i Ff : lJ.4432 rt! A J ,, '/ \ i\A l(\ \,1\ \r \ -101 +-------+--------+----'-lr--tl-'i--------~ '--------+------'------------0 5 10 20 150 C 1 Q) L L O ::, u 0 0 25 m -75 so 25 0 00 I 120 Tine qnS ) (a) Harmoni: cai-ettoftle F~ 1:f Recoraea LOad current ,:) : ff) : 15.l:Bl ll-1D : 12.8278 TD : 0 1D : 12.8278 RM91 : 197152 RMsi: 0 RMsfi: 15 .495 RMS : 15A95 AS.JM : 25 .402 TF : 356155 T : 5518.64 I I I I I I 11.111 I,, ,1ll,l1l l 11 l ll1,1,l1,l 1 l1 l 1 I 720 1320 1920 2520 3120 F z (b) Figure 5 11 a .) Field Recorded Current Wavefo11n. b.) FFT of the Field Recorded Wavefor1n in Fi g ure a. 4 6 r l

PAGE 54

CHAPTER6 LABORATORY TESTING OF THE EV CHARGER 6.1. Testing Procedure and Results The EV 1 charger is a 240 volt split single-phase device In the laboratory t esting of the charger, a nearly fully discharged ( 5% of full charge) battery was connected t o the charger. The voltages and currents on both input lines were measured; the equipment set up for the laboratory testing is shown in figure 6.1. Electrical measurements were taken using the power monitor, automatically once every one-minute i nterval for the entire charging cycle. The fully discharged battery was connected to the charger at 10: 17 AM. The battery reached its almost full charge at 2 PM, a charging time of a little less than 3 hours. The EV 1 charger when charging a fully discharged battery draws a current of 30 amps at 240 volts and is also a source of significant harmonics as shown in tl1e following laboratory resu l ts. These results have been presented in the forrn of waveforms and graphs that depict the electrical characteristics of the EV 1 charger. The first graph, Figure 6 2 shows the power consumed by the EV 1 charger over the entire charging cycle. The load varied from the maximum of 71<:vv at the beginning of the cycle to about lkw when the battery is nearly completely charged. During the entire charging cycle the input voltage remained unaffected as far as either the magnitude or the harmonic distortion is concerned (Figure 6.3). 47

PAGE 55

48 The input current however varied significantly in percentage harmonic content. Figures 6.4, 6.5 6 6 and 6.7 show the ctrrrent waveforms recorded at specified instants of t ime. It can be observed that as the amplitude of the current decreases the waveform distortion becomes more and more visible Figure 6.9 shows the RMS values and figure 5.10 the Total Harmonic Distortion (THD) of the input currents for pl1ases a and b over the entire charging cycle. The figure 6 11 shows the harmonic content of the input current. It was seen that the amplitude of various harmonics throughout the charging cycle does not vary significantly Since the amplitude of the fundamental component decreases with time and the THD, which is a percent measure of harmonics to the fundamental shows a significant increase with time. 6.2. Observation of Laboratory Results The EV charger charging into a fully discharged battery draws a line current of 30 Amps at 240 Volts. This corresponds to a power consumption of approximately 7 kW at time 10: 17 as shown in Figure 6.2. The voltage waveform as shown in Figure 6.3., stays distortion free throughout the charging cycle. The input line current distortion at the beginning of the charging cycle is less visible (Figure 6.4). This is because the amplitude of the line current is high and its THD which is expressed as a percent of the fundamental RMS current, is low. Note that in Figure 6.10 the THD at time 10: 18 is just about 2 % for both line currents Ia and l b, As the charge on the battery increases and the amplitt1de of the line current decreases, the current distortion becomes more and more visible (Figure 6.5 to 6.8). This is verified by Figure 6.10 where in at the end of the charging cycle the THD of the line current climbs to almost 30 o/o. Although the THD of the line current increases as the charge on the battery increases, the actual harmonic content of the line current stays fairly constant. Figure 6.11

PAGE 56

49 shows the line current harmonic content 11p to the 50 th harrnnic. The most significant harmonics are the 3 rd 5 th and the 7 th .The even harmonics are almost not pre se nt.

PAGE 57

Line 1 Line 2 Neutral Currents Power Monitor Digital Oscilloscope Dynamic Signal Analyzer Power Harmonics Analyzer 50 EV Charger 240VAC Single Phase Voltages Figure 6. 1 Experimental Setup fo r EV C h arge r Testing

PAGE 58

C ha rg er Power O u tp u t 8000 V, 6000 I I 4000 ro 3: 2000 \ I I 0 I I I I I I "' I I I I I I I I LO N LO N LO N LO 0 0 ( I 1 N N (Y) (Y) I I I I I I I ( I Figure 6.2. Power Curve for the E V I Charger.

PAGE 59

< 0 rt PJ (0 It) < I 150 +-~-~~--~1I I I --.--r--.--,---1---,---,-,~ -,----,---+~ ~,r---.----, l00 --1----1----1---1----1----1----1----1---1---1---1---1-f--+---1---+---1---1----------+---1--+--+---+---1 50 I --I 0 I ,--1 50 _J_ I I -100 I I -150 -1------,-~ ---. ~ -200 0 l --~----~--, I ---.-I _J 5 I --l I I I --I I I I I I I I I 1 7 1 r 1 I I I I ~L-J..._ __ 1 10 15 Time ms Figure 6.3 Typical Charger Output Voltage Waveform I -I I I I .--1 ,I J -1 1 I I I I _j -L I 20 25

PAGE 60

40 I I 20 --1'.--.;_-= ____:_ __:__1-_..;..-_;._ __ _____c __ ~--+---'--'--'---+--_c_-'<-;_-----'---+----'-----'----j C -, -, (D ::, 0 -20 -40 I I I -60 1 0 I r r I I 1--1 I I I I I I I J __ L_l L j __ 1 5 10 15 20 Time (ms) Figure 6.4. Charger Input Current Wavefortn Sampled at 11: 19 I I I I -rT I I I I I I _) _J 25

PAGE 61

40 .--r--.---,--,---,--,--.---,--,---,--,--,---,---,---,---,--,,--,--.---,---,--,--,--,---, 20 -t---,--t--t-t-~ ---l--+--l--t-------1-+--+----------+----t--t--t--l 10 --f--\-t--+--t--+---t---f--t---+---+--t--,4---l---l--+--t--l---t-\--t---l--+--~---------t 0 -+--10 III -l C I I -, -, (D -20 :::J -30 I I I I I I I I I -40 _j_ L I J..__ l l _J_ J_ l 0 5 10 15 20 25 Time (ms) Figt1re 6 .5 Charge r Input C urrent Waveforrn Sampled at 11 :40

PAGE 62

(i C -, -, (l) ::J rt )> I 10 --~--.--1,--.---~--~-.--.--~-1--.----r+ 5 -l---,--,---,--,---l---.---,--,---,----!----.+---.---,----,---1---.--\--,---,---,---j--,----.--r--r--l 0 -1--\-----.-----,----+-----.------'l---11-t---+----t---t---l---+---\----.---l---t----,,--+---,.---1 -5 - I _,__ .... I I -10 I I I I -15 1 _j_-l 0 5 10 15 20 25 Time ms Figure 6.6. Charger Input Current Waveform Sampled at 12:07

PAGE 63

C: 7 7 ro :J rt )> 15 -.--.---,--,.-r--,---,--,------,-_,---,--,--,--,--,--,---,--,--,------.-,-~-,--,------.--, 10 i----..:-_...!.__ s -..._, 0 1 I I I -5 r TI I I I -1 0 I I I I I I -15 _J t_ J_ L 0 5 I I I I I I I I I I I I 10 Time -1I -r I I I _,_ 15 (ms) I I I I I I 20 Figure 6 7 Charger Input Current Wavefortn Sampled at 12:36 I I 25

PAGE 64

C ., ., l.0 -r----,-----,-----,-----,--.---------,--.--,----,--,-------,--.----,-----~--,---,---, O.S -t--.--.--.r---.--+--.-~-1-.---,--1------.------'n---,--H-1---.r---.--.--,--t-,--.----.--n/-i (1) ::l 0 5 ~--+~ -1--l-+---1--1--1----l--l --j--4-l----1-----+--J_..,_-1--1----1---l--j ) +-+--,----, -1 .0 0 I _._ __._ __, __ j__-+ __ 1__--1.. __ 1__..1... _jf---..L-_l_~.l--_J--+--1-~---l 5 10 Time 15 (ms) 20 Figure 6.8. C har ger Input C urrent W avefo r1n Sampled at 13 :53 l 25

PAGE 65

35 30 25 Ill 20 a. E < 15 10 5 0 10:18 Current RMS 11 RMS i2 RMS --\ 10:42 11:02 11:22 11:42 12:02 12:22 12:42 13:02 Figure 6.9. The Charger Current RMS Variation over the Charging Cycle --.... 13:22 13:42 VI 00

PAGE 66

35 30 25 o 20 :::c ... 1 5 1 0 5 0 10 : 1 8 Current THO -l a T H O ---I lb T H O ,. ... I 10: 4 2 11 : 02 11 : 22 11:42 12:02 12 : 22 1 2 : 42 13:02 Fig u re 6 1 0 The C h ar g er Current THD Variat i on Over t h e Char g in g Cyc l e ,, .. 13:22 13:42

PAGE 67

Ia Max Harmonic s 1.2 .------------------1 0.8 Ill a. E 0.6 < 0.4 0 2 1--1 I I CO 0 r-1 1-t N r-1 I I 1 1 1 1 1 , 1 , 1 , 1 , 1 1 1 s:;t lO CO O N s:;t lO CO r-1 r-1 r-1 N N N N N H a rmonic Figure 6.11. The Charger Input Current FFT. I I 1 1 , 1 1 , 0 N s:;t lO CO 0 "1" "1" "1" "1" s:;t LO 0

PAGE 68

CHAPTER 7 COMPUTER SIMULATION This section contains the results of computer simulations of the power distribution system with domestic and Electric Vehicle (EV) loads. The domestic load included in these simulations is arbitrary. The EV model used is based on the laboratory data obtained by testing the actual charging characteristics of an EV charger used by the GMC EV 1. The aim of this simulation is to see the effect of EV s connected to the distribution system. The transformer current with EV s present in the system is simulated. The FFT of the transformer current with and without EV s are shown and different harmonic components are found to study the effect of these harmonics on the distribution transformer in terms of the derating of transformers. The power distribution system with an EV charger load was simt1lated using SABER. This section includes the simulation results and discussion about the simulation model. 7.1. Component Models 7 .1.1. System Model The system model consists of the distribution transforrner, domestic load and the EV charger and battery load. Figure 7 .1 shows the layout of the simulated power distribt1tion system. The distribution transfor1ner has typical domestic loads connected at 61

PAGE 69

' El e l la .. fJ lt t ,I, :a ?) .. E ~= {o,) ,I, 1} I I .. Figure 7.1. Distribut i on System Mode l

PAGE 70

6 3 FU Harmonic Current Source I I J \ I\ t1 SM!A r---. I\ ti J \ r Figure 7.2. Simplified Model of the Simulated EV Charger

PAGE 71

64 120 volts between the center tap and the outer terminals. The outer terminals, which are at 240 volts, are connected to a full bridge rectifier. The bridge rectifier is used to model the inductive charger u sed by the EVI to charge its batteries. 7 .1 2. Electric Vehicle Charger And Battery Models The EV charger is modeled on the basis of laboratory data obtained during the charging of a fu ll y discharged EV battery connected to charger. The full bridge rectifier represents the charger. The inpu t to the rectifier is 240 V from the distribution transformer. The EV charger and battery model is shown in figure 7.2. The battery has been modeled as a series combination of a resistor and a capaci t or. The resistance represents the loss component of the battery. It is presumed that the charging pattern of the capacitor shall closely represent the charg in g pattern of the EV battery. The values of the resistor and the capacitor can thus be chosen to achieve the desired charging current and charging time. In the simulation while the full wave bridge rectifier represents the EV load, a composite current source has been used to represent the har1nonic effect of the EV charger. The composite current source is injecting har1nonics of the same amplitudes and r espective orders as shown in figure 3 .2. The actual waveform of the harmonic current injected by the compos it e current source is shown in Figure 7.3. 7.1.3. Transformer Model The transforrr1er used for the s imul ation is a gap l ess three winding transfor1ner w i th one primary and two secondary windings. The transfor1ner in the simu l ation has the following specifications.

PAGE 72

65 3 I -3 ... ..... ... Figure 7.3. The Harmonic Current Injected by the C harger

PAGE 73

66 These specifications were provided by FPL. KVARATING 37 Turns Ratio 7260, 120 / 240 Material of the Core Grain oriented silicon steel Thickness of Lamination .009 inches Net cross sectional area of the core 28.9 sq. inches Path length 25.74 inches Magnetization Saturation 20.0 Kilogauss 7.1.4. Domestic Load Model The domestic load has been modeled as a combination of some lighting and heating load represented by resistive load. The load has been modeled such as to exhibit the characteristics similar to the actual load characteristics as determined from the field data. Figure 7.4 shows a typical domestic load model used in simulation. 7.2. Mathematical Model Of The EV Charger. The EV charger can be represented by a single-phase bridge rectifier with an R-C load, as shown in figure 7.1 can represent the EV charger. A current source model is used for the injected harmonics when there is very little voltage distortion [9]. Since the EV charger is

PAGE 74

67 . .J I ..I I Figure 7.4. Domestic Loads Model known to cause very l ittle voltage distortion, the hannonics injected by the charger into the line can be modeled by a composite current source. The behavior of the charger can thus be determined by a second order system. 2 d Ve dvc LC 2 +RC ~+vc=vs dt dt (1) where the transformer secondary vo l tage t he capacito r voltage L the eq u iva l ent battery i nd u ctance Ct he equivalent battery capacitance R the equivalent battery internal resistance

PAGE 75

68 The net harmonic current from the cornposite current source can be expressed as: 50 'H = k=2 where I mk sin kmt iH is the net harmonic current k is the order of tl1e har1r1onic (2) 1,nk is the amplitude of the k th harmonic current The SABER software package was used to simulate the system, both transient and steady state operating conditions, including RMS values, THD etc can be obtained using SABER. In the battery model the Lis small and hence can be neglected. Therefore the charging time constant (T) for the battery can be given by: T=RC (3) 7.3. Simulation Results In the computer simulation two types of loads have been considered, viz. the lighting loads and the EV load. The simulation results with two houses and a single EV load are presented. The lighting load current is taken as a pure sinusoidal current. The current harmonics injected into the mains by the EV charger are shovvn in figure 7.3. This is based on tl1e experimental results obtained by testing the EV 1 charger in the laboratory. The harmonic components of the charger input currents (figure 7.5) are supplied by the composite current source. The charger input current in the simulation is

PAGE 76

69 controlled by varying R. We know from the lab test data that the maximum current that the charger draws (for fully discharged battery) is approximately 30 Amps. By modeling the GM EVl charger as a combination of a full wave bridge rectifier with a resistive and inductive load and a composite current source, simulation results very close to the lab rest1lts were obtained. A detailed discussion on the compariso11 of tl1e laboratory results the simulation results and the field-recorded data has been included in Chapter 7. The simulation run s of the distribution sy s tem made it possible to check the various load scenarios in the distribution system. The system design can now be checked for an arbitrary number of EV s connected to it. This data obtained from simulations is important and makes it possible to optimally design/redesign future distribution systems. This includes replacing and/or addir1g transfo11ner units in the existing systems in order to make them able to withstand the new load conditions caused by the EV load and the harn1onic distortion associated with it Based on the field and laboratory data the EV charger has been succe s sfully modeled The simulation results are close to the laboratory results obtained during the actual charging of tl1e EV battery The waveforms on the input side of the transfonner have a low order current harmonics. 7.4. Simulation ofSl1ort Charging Time Chargers The computer simulations were used to check scenarios when the charging time for charging the battery is much shorter. The GM EVI charger takes about three and a half hours to charge a fully di s charged battery (from about 15 % of full charge). The long

PAGE 77

70 charging time is a drawback with the EV and research is being done to reduce this charging time substantially. The fast charging time chargers are not commercially available and their behavior in terms of harmonic effect is not known. The computer si1nulations have been used to provide an estimate of the effect of the fast charging EV on the distribution system. The rest1lts of these simulations are based on the researcher intuition and experience of working with highly non-linear loads. The fast charger study is aimed at providing an idea about how such a technology, when developed can effect the distribution system. These results need to be verified once the fast charging technology becomes available. 7 4 1. Simulation Methodology For Fast Chargers The charger model for the fast charging charger is essentially the same as that for the nom1al charging time charger. The difference however is only in the values of the resistor and the capacitor that represents the battery. Also the amplitude of the ha1nonic currents has been increased proportional to the decrease in the charging time. Thus the fast charging time charger has also been simulated as a combination of a full-wave bridge rectifier and a harmonic current injection source. As mentioned earlier, these simulation results are based on intuition and experience of working with non-linear loads. These results can only be verified once fast charging technology becomes available. Figures 7 .6 through 7.8 show the current waveforms and their respective FFTs for simulated scenarios where the charging time is reduced as indicated in the figures. The distribution transforrr1er derating corresponding to this charging scenario is shown in Section 8.2

PAGE 78

~s. I (A) ; '(a) 2M.,C, ,--------------,, TrlW111f6itiW CU,t anl 100 0 I r I .. 0..0 n ,, -100.0 -200.0 r==== ::;==== =r= = ====;:===~ I I I I J o.o o..w o.o, 0.01 a o.oa t() Figure 7.5 Transfo r 111er Current for the Fast C h arging Charger with C h arging Time of 45 minutes. Meg(A) f(Hz) eo.o. ,---------------, T 1U\lda.-.,.o..,,ent I 40..0. ao.o. 0.QJ~~~--__, ........ ...,.._====;:!===-=::.==~ r ' 0.-at c..i 0.78k 't OIC 1 ZSIIC 1.Slk t(HK) Figure 7.6 FFT of the Transfo r 1ner Current for the Fast Charging Charger with C h arging Time of 45 min u tes. 71

PAGE 79

7 2 {A.) : t() ,------------------Tran I b 11. If' Cu,,aat 100.0 0.0 .0 a.o I -. . . . .. . . . .. . 0.005 . 0..01 t() . . . .... . . 0.015 0.02 Figure 7.7. Transformer C urr e nt fo r th e Fast C har ging C har ger with Charging Time of 1.2 Hours. 30.0 20.0 10.0 0.0 I O.&k I 1.0k f(Hz) : f(H%) I I I I I J J I 1.&k 2.0k F i gu r e 7.8. FFT of the Transformer C urr e nt for th e Fast C h a r ging Charger with Chargi n g Time of 1.2 hours

PAGE 80

I CHAPTER8 COMPARISON OF FIELD, LABO RA TORY AND SIMULATION RESULTS In this section the field data and laboratory test results are compared with those done by computer si1nulation to present the effect of EV loads on the power distribution transformer. The aim of this project is to study the impact of the EV load on the power distribution system. In summary, this goal has been accomplished by recording the field data with and without the EV load, recording the behavior of the EV charger by the laboratory testing of the GM EVl charger charging a fully exhausted battery. The data acquired from the field provides the information abot1t the characteristics of the distribution loads when EV load is a part of the load and when there is no EV load. The laboratory testing of the cl1arger provides the input current characteristics of the charger and the line current harmonic content of the charger. Using the data from the l8boratory test, the computer simulation model of the EV charger has been built. The details of the computer model have been discussed in Section 8 .1. Thus the laboratory testing of the charger forrns the basis on which the computer model is built. The charger computer model has been incorporated into the overall model of the distribution system. The distribution system consists of the transformers, feeders different types of normal domestic loads and the EV loads. First of all, in this section, the accuracy of the computer model of the EV charger will be established. This has been achieved by comparing the laboratory test results of the 73

PAGE 81

74 EV charger with the computer simulation results. Second, the similarity of the field recorded data with the simulation results will be discussed. 8.1. The Field Data The field data was recorded in the form of one-cycle snapshots of current and voltage as well as power factor at the distribution transfor1ners at four different locations in South Florida and then transmitted to a con1puter over the telephone line. The detailed description about the data monitoring and acquiring methodology are mentioned in Section 4.2. 8.2. Laboratory Test For the laboratory test the EV charger is selected as a 240 volt, split single-phase device as described in the previous section. The voltage and current input into the charger were monitored and recorded over the full charging cycle for the EV, which was nearly completely discharged to a level of 5% of full charge. Electrical meast1rements were taken using a power monitor automatically every minute over the entire charging period of about three hours. From the measurement it was found that the current supplied to the EV charger varies throughout the charging cycle. The charger current variation is shown in Figure 5.9. The actual amplitude of each of the harmonics in the charger current does not vary significantly during the charging cycle.

PAGE 82

75 8.3. Simulation Results The power distribution systen1 consisting of the feeder, distribution transformer the EV charger, together with domestic loads (mainly the lighting lamps) wa s simulated using SABER. The model s of various components used for simulation are based on both field data recorded and the laboratory data. To simt1late variation of input current and distortion at different instants over its chargi11g period an adjustable resistor, as well as a composite current source which is based on the FFT analysis of the laboratory test data, is introduced in the loop. By vary in g the resistance, the identical situations to tl1ose in the laboratory test can be established. Therefore current waveforrn snapshots over a whole cycle of 60 Hz are obtained and then compared with those recorded in the laboratory test. 8 4. Compariso11 The field data laboratory test data and simt1lation results are compared in terms of the waveforrn snapshots with the same magnitude of charger input current for all three cases viz. the field data, the laboratory test data and the simulation results. The aim is to establish the r e li ability of the laboratory test data; the simulation model and simulation results by comparing them to each other and the field recorded data. Some of the waveforms that have been included in previot1s chapters have been put in this section for convenience of comparison.

PAGE 83

F-e 1:i Recorded currentwitl cV LOad 30 rta 11\SSD-rul'.E Max: 27.832 ,...... 20 10 C: Q) '0 '::, 10 u Mil: -Il2734 I +------H--------'<....----+------------' A\<]: 15.?fB? .. 1 ..t6: Il2734 I J RMS: 18.1012 ----F-----+------>r-+------~ CF : 1.672'16 FF : 114 \,. 2 0 -1---------l-------_;._------l-----_p...-----' ~. 30 -r-----------------------------W -------+------~,-----------0 5 1 0 1 5 20 150 -...J 125 C 100 Cl) L L 0 7 5 ::J u 050 025 0 00 120 Tile (a) I I ) : V 91:i3 l'HD : 10.ll:8 lD : 0 .. 1D : 10.7".E RMS, : l.94ll2 RMsi: 0 .. RMSfl: 18.])}.) RMS : 18.])}.) A9JM : 7J .77/l TF : 28'.3335 I .. r : 5237 23 1 I I I I I I I I 1 . 111 I., 11.11 .1 1 I I 111. I 1. l.1 I .. I .. 11. I 7 20 1 320 1920 2 520 3120 (b ) Figure 8 I a ) Field Recorded Current b .) Waveform FFT of th e Field Recorded Wavef o rm 7 6

PAGE 84

1 0 6155 r : 5518.64 I I I 1 I .111, I., 111 I I .I 111, 1111, .1. I 1. I , I 1. I i ' 720 1320 1920 2520 3120 F z (b) Fi g ure 8 .2. a ) Field Recorded Cu rr ent Wavefo r rn b.) FFT of t he Fie l d Reco r ded Waveform in Figure 8 .3. 77

PAGE 85

4 0 I 30 .J.... __;___;_ _;__.:..__--1_ _;_ __ ___;_ ___ +-___;_----,----. ,-+------------------,---------1 J 20 .J...._.!__;_ _..___.l_--l_ _!_ __ ___J. _!__+-___J.-,Ll---'---''----4---",,-''----'---'-----'--'------'-----, 10 ..L..)... _______ ..;._ ______ -,--I---I-,---,----,------+---->,,--------,--.,----,,-------, I I O -l--~__:. ___ ....:....._+-__; _;_ _;_ __ 4,:_______ ,----+ _____ _,,_ ____ -,-----,------, ,;' -20 -r-----~-----+---.,---,.:___--+--,-----:----+--------...:.------/--~ -30 -t---;----;----..::::..-;,=z"""--;-~--;--.---.---,-,---.~.----'-=-;:,..::--~'---4 0 '--------''-------~----~-------'------=---=---'-------'------0 s 1 0 15 20 25 Time ( m s) Figure 8.3 Laboratory Test Transfor1ner Current at 11 :40 AM 15 1 0 + -.!_---'----'-.l..--t----'------l--'---i------l--'---~=6..:::;::+--'----'---'-------'--"-----'-~ s +--'---'---'--'---.....J..--'-----'--'------1---->'"----'----+---'<----'-------'-------s I I -10 -~ .. I I 0 5 10 1S 20 2S Ti me (ms) Figure 8.4 Laboratory Test Transformer Current at 12:36 PM la Mas Harm o n ic s I 2 l O 8 ... e o 6 ,c 0 4 0 2 0 "' .. "' J_1_1, __., 1 ___ 1 _. ._. ---L _. --_ __ _a-.i 2 g : 00 0 "' .. "' .. Hll rmonlc Figure 8.5. FFT of the Current Waveform of Figure 8.6. 78

PAGE 86

' 8tmulllad Chargal' CUnwtt )1Jlt11 cocapo."'idlng FF r (COl'l 18po.,dlng 10 f sh Tl 11 Cu. t trrt 12:N) I u.g(A) ; f(Hz) r---------------------. l(ahort.tx,Jn2) 2J). 0.0 =:::==.:: ===.;;;;;;;=: ;;;:::::!;: =:=:::::!::::!=~-=====:!;:::!-~== I I I I l O..Sk 1.ok 1 .&k 2.0k :Z.5k 3.0k f(HI) (A) : qa) r--------~~ -------, l(9hort.tx.Jn2) 0.02 0.025 o..oa O..OS5 Figure 8.6. The Simulated Current Wavefonn and the Corresponding FFT 79

PAGE 87

1.0 0.0 20.0 0.0 40..0 1"........_. eurr.ntW.Tbna and lta PFI (Cot: oepo."dng to 11 :40 Cunn Wwefot ,,, of the Lab T..c) llag(A) :f(Hz) ,---------------~ahort.bt..ln) DA 1.Gk 0.0 1A f(Hz) o.m 1(1) r Ole ?Sic S.Ot, (A) : t() a(1hort.tx_ln) O.GI 0.04 Figure 8.7. The Simulated Current Waveform and the Corresponding FFT 80

PAGE 88

81 The waveforrn of figure 8. lshows the fie l d recorded current for the distribution transformer feeding two homes (FPL_l). This is the distribution transformer current when the load on the transformer contains EV load. Figure 8 2 is the FFT of the wavefor111 in Figure 8.1. The observation after comparing the current waveform of Figure 8.1 and its FFT in figure 8.2 with the laboratory recorded waveform in Figure 8.5 and i ts corresponding FFT in Figure 8.6 are p r esented in a tabular forr 1 1 in Table II. Both the currents (of Figure 8.1 and 8.5) are then compared to the waveform of Figure 8.8. Figure 8.8 shows the simulated transformer current. Also compared are the FFTs of the tlrree waveforms. (Figure 8.2, 8.4 and 8.8). Table. II Laboratory Test Simulation Results Field Data (Figure 8.5 and 8.7) (Figure 8.8) (Figure 8.1 and 8.2) Peak V alu e o f 30 Amps 29 Amps 27 82 Amps C urr e n t RM S V alu e o f 19.8 Amps 20.2 Amps 18.1 Amps C urr e nt C r es t Fa ct o r 1.515 1.435 1.537 THD 14.1 11.2 12.8 Holding the field recorded data as reference the relative percent values of the laboratory and simulation r esults is shown below in Table III.

PAGE 89

Table. III Laboratory Te st Simu l ation Res u lts (Figure 8.5 and 8. 7) (Figure 8.8) Peak Value of Current 7.8% 4.24% RMS Value of Current 9 393/o 11.6 3/o Crest Factor -1.433/o 6.636% THD 10 1 53/o -12.5 % Look in g at the wavefo r ms figu r e 8. 1 through figure 8.8 and Tables II & III, we note that the field recorded and t l 1e sim ul ated cur r ent contain tl1e domestic load current. 82 The l aboratory test waveform h oweve r is only the charger cur r ent. I n th e simulation the domestic l oad is t aken as l inear load (Figure 6.4) while i n reality (Fie l d recorded) the domestic load s may contain some nonlin ear loads. Thus the THD of the fie l d current is higher than the THD of the simulated current. The laboratory test current consists of only the charger current.

PAGE 90

CHAPTER9 DERATING OF DISTRIBUTION TRANSFORMERS 9 1 Method of Calculation 9.1.1 ANSVIEEE C 57.110 Method The method u se d for calculating the derating of the distribution transfo11ner is the ANSI/IEEE C 57.110 (see Section2.3.l) Using the method described in Section 2.3.1, the derating has been calculated as a function of the current drawn by the charger. The derating curve i s s hown in Figure 9.1. In the calculations the derating is ca l culated taking the cl1arger current as tl1e base. So fh in equation 2.1 is calculated dividing each hru1nonic component by tl1e amplitude of the fundamental charger current. Thus, in this method the value of fh for the fundamental component is always 1. This is the procedure employed by ANSI/IEEE C57 110 method From the curve of Figure 9.1. we get the minimum and maximum limits of derating. Table IV. Maximum and Minimum Derating Using ANSI/IEEE C57.110 Maximum Derating Minimum Derating 99.95/o 83

PAGE 91

en C: I ca L. Q) C Derating of Distribution Transformer 1 ...------------------------. 0.8 -------~ ----------0 6 --0 4 -1 ---------------11 0 2 -----------,--, Figure 9 .1. The Transformer Derating Curve for the Complete Charging Cycle

PAGE 92

85 9 .1.2 Modified ANSI/IEEE C 57 110 Method The ANSI/IEEE C57.110 method calculates the derating of the transfor1ner based on the charger fundamental ctrrrent. In reality the transformer in a distribution system is supplying other loads. Also from figure 9 .1 it can be seen that the transformer derating at 13 :53 hours is as low as 203/o. This calculation is done taking the charger fundamental as base, and from Figure 5.9 we see that the charger current at 13 :53 hours is just about 5 Amperes. This would for1n only a small percentage of the total load on a distribution transformer of 37 KVA rating. Thus, although the derating based on C57 11 0 method is very high the effect of the small current on the transformer would be much less significant. In the modified ANSI/IEEE C57 110 method the calculation of transformer derating is done based on the maximum current carrying capacity of tl1e transfo11r1er. Since the transfonner derating should be a concern when the l oad on the transforrner is fi1ll or close to full l oad this method calculates the derating based on full load current of the transformer. The derating of a 3 7 KV A transfonner using the modified method is shown in Figure 9 .2. 9.2. Transfonner Derating with Fast Time Chargers Based on the simulation results of section 6.4, the modified ANSI/IEEE C57.110 meiliod has been used to calculate the derating of distribution transfo1111ers.

PAGE 93

86 The derating has been calculated based on the 3 7 KV A distribution transformer whose specifications are given in Section 6.1.3. Figures 9.3 and 9 .4 show tl1e transformer derating for the two simulated load scenarios. 9.3. Conclusion The derating of the distribution transforrner based on ANSI/IEEE C57. 11 0 method is significant but it does not consider the current carrying capacity of the transforrner when calculating the derating. The modified approach calculates the derating of the transformer based on the full load current of the transforrner. It is obvious that based on this method we see that the GM EVl charger does not significantly derate the transfortner. One thing that may be noted here is that the ANSI/IEEE C57 .110 method is generally applied to fmd the derating of large power transforrr1er and gives highly conservative results for distribution transfortners. It is strongly recommended that tl1e equation for derating of distribution transformers be experimentally determined wiili KV A ratings of concern to FPL. (See Chapter 10). The derating for the fast charging time chargers as seen from figures 9.3 and 9.4 i s much higher. It is expected that the fast chargers will be more prevalent than slow chargers if consumer preference is not modified by the utility industry.

PAGE 94

Derating of Distribution transfonner 1.0001 1 ----------------------------~ C') = .... ; 0.9998 -a --.-t ---=---&--1 --------------1 C 0.9997 ------------------0. I 11 lllllflllUI 111 1111111111 I 1111111 I 111111 Figure 9.2. Transformer Derating using the Modified ANSI/IEEE Method

PAGE 95

4 1 Ofo Derating of Transformer 11 % Derated Value 89% l Derated Value J a. Derating of Transform er l e r ated Value b. Derated Va l ue 5 9 Ofo 88 Figu r e 9 4. Derating Corresponding Short Cha r ging Times a.) 1.2 Hours b.) 45 minutes

PAGE 96

CHAPTERlO CONCLUSIONS 10.1. Summary and ConcltlSions. The objectives of the project outlined in section 1.3 have been acl1ieved. The behavior of four distribution nodes (FPL_l to FPL_ 4) has been monitored. The load behavior with and without the EV load has been recorded. Volume II of the final report contains samples of the recorded data with an executive summary. By performing the laboratory test on the GMC EV 1 charger charging a fully discharged battery, the charging characteristics of the charger has been deter1nined. The charger takes about 3 hours to charge a nearly fully discharged (1 So/o charge) battery to full (1 OOo/o) charge The harmonic spectrum of the transformer current charging the discharged battery has been determined. The snap shots of the current waveforms have been recorded and studied. Some of these waveforms have been included in this report Based on the laboratory data a computer model of the charger and the distribution system has been developed. Comparing the simulation results to the laboratory results and the field data has validated the simulation results. The simulated currents are roughly 4 to 1 Oo/o away from the field recorded ct1rrents, in their magnitude and distortion. Based on the field data, laboratory results and the simulation results the derating of the distribution transformer has been calculated. Based on the field data, laboratory results and the simulation results the derating of the distribution transformer has been calculated. The two methods suggested 89

PAGE 97

90 for the ca lcul ation of derating have been presented and compared. The maximum and minimum limit s of derating have been suggested. Using tl1e frrst method the maximum derated value is 80% and the minimum value is 1 %. Using the modified method the derating remains close to 99 % tluoughout the charging cycle. The transformer derating results corresponding to fast charging (45 minutes and 1.2 hours) chargers has been presented. A Microsoft Exce l based program has been provided to FPL. This program shall enable the user to find the derating of the transformer using the harmonic content inf ormation of the load on the transformer. 10.2. Scope For Further Research The transformer derating has been calculated using the IEEE / ANSI method. The formula is an empirical formula To get a h ighly accurate measure of the derating of transformers, laboratory tests should be performed. These tests can be perfor1ned in the Power Laboratory in the Electrical and Computer Engineering Department at the University of Florida. Using the real time program R4, the computer can generate low voltage cont r olled harmonic signals. These signals will be then amplified using the power amplifier in the laboratory and fed to the transformer under test. Temperature sensors placed selectively on and inside the transfor1ners will sense the rise in the core and winding temperatures. The change in the temperature rise with a change in har1nonic content from pure sinuso id al to very high harmonic content can be gauged. The transfo1111ers can be

PAGE 98

91 tested for different harmonic scenarios and loads. Thus an exact formula for the derating of transformers can be evolved. As the hannonic generation can easily be controlled using the computer, affect of individual harmonics can also be studied. This information can be particularly useful to FPL for the design of harmonic filter s that may need to be installed in place s wl1ere tl1e load is highly sens itive to h ar monic s. Also, at sites where tl1e load is known to generate a significant amount of a particular harmonic these te s t s can provide a strong b as i s for calculating the derating of equipment and thu s the design of appropriate transformers.

PAGE 99

REFERENCES 1. J Pages, ''On the quality of public electricity supply'' Revue Francaise de l 'E lectricite vol.41, no.223 p. 14-21, 1968. 2. C. H. Hilg e r ''Wha t is the quality of electric power?'' Elektroteknikeren, vol.68, no.19, p. 418-22, October 7, 1972. 3. H Heikkila, ''Harmonics in a power system'' Saelikoe, vol.49, no. l, p. 27-30, January 1976 4. G. L. Johnson, A. N. Mortensen, ''A power system digital harmonic analyzer'' IEEE Trans. Instrum. Meas. (USA), IEEE Transactions on Instrumentation and Mea s urement, vol.37, no.4, p. 537-40, 1988. 5. C. K. Duffe y, R. P Stratford, ''Update of harmonic standard IEEE-519 : IEEE Recommended Practices and Requiren1ents for Har111onic Control in Electric Power Systen1s. Industrial Applications Society 35t11 Annual Petroleum and Chemical Industry Conference pp. 249-55, 1988. 6. P. Meynaud, G. Hensman, J. L. Javerzac, L. Lagostena, M J. Lahtinen, F. Lionetto, W. Monbauer P Herager A. Robert, J. Smid, G. Allen, C. Britten, A. Emanuel, K. Murotani, P. Ribeilo, P. Wright, Yu Zhelezko, ''Equ ipment producing harmonics and conditions governing their connection to the mains power supply." Electra, no.123, p. 20-37, March 1989. 7. T. S. Key, R J. Ferraro, R. W. Taylor, ''E lectric power for the 21st Century: changing load s -changing needs.'' Power Qt1ality '89. Official Proceedings of the First International Conference, p. 613, 58-73, 1989. 8. R. J. Gille sk ie ''Power quality and emerging electrotechnologies their relationship impact and suggested precautions." Power Quality '89. Official Proceedings of th e First International Conference, p. 613, 108-23, 1989. 9. W.E. Kazibwe, R. J. Ringlee, G. W. Woodzell, H. M Sendaula, '' Power quality: a review." IEEE Computer Applications in Power vol.3, no. l, p. 39, January 1990 10 M. S. Padbidri, D. G. Pantulu ''Studies on harmonic pollution in power systems." Power System s for the Year 2000 and Beyond Proceeding s of the Sixth National Power System s Conference, p. xvi+704, 31-5, 1990. 92

PAGE 100

93 11. A. Domijan, Jr., E. Embriz-Santander, ''Harrnonic mitigation techniques for the improvement of power quality of adjustable speed drives (ASDs)." APEC '90. Fifth Annual Applied Power Electronics Conference and Exposition. Conference Proceedings. p. 839, 96-105, 1990. 12 D. Ollivier, ''Simulation of harmonic disturbances (in power systems)." Bulletin de la Direction des Etudes et Recherches, Serie B no.I, p. 1-31 1991. 13. J. K. Phipps J. P. Nelson, P. K. Sen, '' Power quality and harmonic distortion on distribution systems.'' 1991 Rtiral Electric Power Conference. Papers Pre se nted at the 35th Annual Conference. p. xi+l64, Dearborn, MI USA; 28 30 April 1991. 14. D A. Van Den Eijkel, D A. Marshall, ''The simulation and measurement of system harmonic impedance ." CIGRE. Proceeding s of the 33rd Session. International Conference on Large High Voltage Electric Systems, p. 2 vol. 2772, 36-302 / 1-4 vol.2 Pari s, France; 26 Aug.-1 Sept. 1990 15 H M Beide s, G. T. Heydt ''Power system harmonics estimation and monitoring .'' Electric Machines and Power Systems, vol.20, no.2 p. 93-102. March-April 1992. 16. B. M. Hughe s, J S Chan, D. 0 Koval, ''Distribution customer power qualit y experience." Conference Rec ord of the 1991 IEEE Industry Applications Society Annual Meetin g. p. 2 vol. 1910 1556-63 vol.2.Dearbom, MI USA; 28 Sept.-4 Oct. 1991 17. R R Verderber 0. C. Morse, W. R. Alling '' Harmonic s from compact fluore scent lamps.,, Conference Record of the 1991 IEEE Indu st ry Applications Society Annual Meeting. p. 2 vo l. 1910, 1853-8 vol.2. Dearborn, MI USA; 28 Sept -4 Oct. 1991 18. R. Christiansen, ' Effects of high levels of harmonic s from lighting equipment and systems.'' Conference Record of the 1991 IEEE Indu s try Applications Society Annual Meeting (Cat. No.91CH3077-5) p. 2 vol. 1910, 1859-62 vol.2. Dearborn, MI, USA; 28 Sept. 4 Oct. 1991 19 K. Johnson, R. Zavadil, ''Assessing the impacts of nonlinear loads on power quality in commercial buildings-An overview.'' Conference Record of the 1991 IEEE Industry Application s Society Annual Meeting (Cat. No.91 CH3077-5), p 2 vol 1910, 1863-9 vol.2.Dearborn, MI USA; 28 Sept.-4 Oct. 1991. 20. R H Kitclun S. Etminan, ''Co mputer si mulation of harmonic waveform distortion, and comparison with mea s ured result s." IEE Colloquium on 'Simulation of Power Systems', p. 62, 4 / 1-7 London UK; 2 Dec 1992. 21. R. S. Wallace ''T h e harrnonic impact of variable s peed air conditioners on residential power di s tribution." APEC '92. Seventh Annual Applied Power Electronics

PAGE 101

94 Conference and Exposition. Conference Proceedings 1992 p. 858, 293-8.Boston, MA, USA; 23-27 Feb. 1992 22. R. R. Verderber, 0. C. Morse, W. R. A l ling, ''Ha r mon i cs from compact fluorescent lamps." IEEE Transactions on Industry Applications, vo l. 29, no.3, p. 670 4. May June 1993. 23 D. D. Shipp, W. S Vilcheck, M. E. Swartz, N. H. Woodley, ''Expert system for analysis of power system harmonics.'' Conference Record of the IEEE Industry Applications Society Annual Meeting p. 2 vol. xxix+ 2071, 1314-21 vol.2. Houston, TX, USA; 4 9 Oct. 1 992. 24. R. D. Henderson, P. J. Rose, ''Harmonics : the effects on power quality and transfor1ners." IEEE 1993 Annua l Text i le, Fibe r and Fi l m Industry Technical Conference p. 96, 4/1 11 Atlanta, GA, USA; 4-6 May 1993. 25 J. K. Phipps, J. P. Nelson, P. K. Sen, ''Power quality and har1nonic distortion on distribution systems." IEEE Transactions on Industry App l ications, vo l .30, no.2, p. 4 76 84. March-April 1994. 26. Kao Chen, ''The i mpact of energy efficient light i ng on power qua l ity." Power Quality '93 USA. Official Proceedings of the Seventh International Povver Quality Conference, p. 875, 68-76 Power Quality '93. Official Proceedings of the Seventh International Power Q u ality Conference. Irvine, CA, USA; 24 29 Oct. 1993. 27. P. P. Barker, J J. Burke, R. T. Mancao, T. A. Short, C. A. Warren, C. W.Burns, J J. Siewierski, ''Power quality monitoring of a distribution system IEEE Transactions on Power Delivery, vol.9, no.2, p. 1 136-42. April 1994. 28. B. Prokuda, ''Power quality site surveys for industrial and commercial bu i ldings." 1994 IEEE Industrial and Commercial Power Systems Teclmical Conference. Conference Record. Papers Presented at the 1994.Annual Meeting p. 292, 97 103. 29. R. D. Henderson, P. J. Rose, ''Har1nonics: the effects on power quality and transformers.'' IEEE Transactions on Industry App l ications, vol.30, no 3, p. 528-32. May-June 1994. 30. C A. Warren, C. W Bums, ''Home power quality the Niagara Mohawk study.'' Proceedings of the 1 994 IEEE Power Engineering Society, Transmission and D i stribution Conference. Proceedings of IEEE/PES Transmission and Distribution Conference. Chicago, IL, USA; 10 15 April 1994. 31. J. Burnett, ''Survey of power quality in high-rise air-conditioned buildings." Fifth International Conference on 'Power Electronics and Variable-Speed Drives' (Conf. Puhl. No.399), p. x i x+715, 163 8 London, UK; 26-28 Oct. 1994.

PAGE 102

95 32. E. W. Gunther, H. Mehta, ''A survey of distribution system power quality preliminary results." IEEE Transactions on Power Delivery, vol.IO, no.l, p. 322-9. Jan. 1995. 33. K. H. Sueker, ''Co mments on "Harmonics: the effects on power quality and transforn1ers" IEEE Transactions on Indu stry Applications, vol.31, no.2, p. 405-6. March-April 1995. 34. D. D. Shipp, W Vilcheck, M. E. Swartz, N. H. Woodley, ''Expert system for analysis of electric power system harmonics.'' IEEE Indu stry Applications Magazine, vol. I no.2, p. 34-9 March-April 1995. 35. P Steciuk, J Burke, ''The impact of railroad electrification on utility system power quality." PCIM/Power Quality '94 with Mass Transit System Compatibility USA. Official Proceedings of the First International Mass Transit System Compatibility Conference, p. vii+412, 61-72. Dallas TX, USA; 17 22 Sept. 1994. 36. F. J. Bourbeau, ''Powe r quality of electric bus battery chargers." PCIM /Powe r Quality '94 with Mass Transit System Compatibility USA. Official Proceedings of the First International Mass Transit System Con1patibility Conference, p. vii+412, 26173.Dallas, TX, USA; 17-22 Sept. 1994. 37. N G Hingor ani, '' Introducing custom power. IEEE Spectrum, vol.32, no.6 p. 41-8. June 1995. 38. C C. Chan, and I<. T. Chau., ''An Overview of Power Electronics in Electric Vehicles." IEEE Trans. on Indu strial Electronics, Vol. 44, NO. 1, pp. 3-12 February 1997. 39. A. Domijan, Jr, and E. E. Santander., ''A Summary and Evaluation of Recent Development s on Harrr1 onic Mitigation Techniques Useful to Adjustable Speed Drives." IEEE Trans. on Energy Conversio11, Vol 7, NO. 1 pp. 6471, March 1992 40. J. K. Phipp s, J. P. Nelson, and P. K Sen., ''Power Quality and Harrnonic Di s tortion on Distribution Systems." IEEE Trans. on Indu stry Applications, Vol. 30, NO. 2, pp 476-484, March/April 1994. 41 G. W. Massey., ''Estimation Methods for Power System Har111onic Effects on Power Distribution Transformers." IEEE Trans. on Industry Applications, Vol. 30, NO. 2, pp. 485-489 March/ April 1994. 42. G. W. Massey., ''Non-Linear Load Harmonic Co11tent Estimation Techniqu es for Use in Specifying Power Di str ibution Transformers." Master Thesis, univer s ity of Missouri Columbia, December 1992.

PAGE 103

96 43. T. Q. Tran, L. E. Comad, and B. K. Stallman., ''Electric Sl1ock and Elevated EMF Levels due to Triplen Harmonics." IEEE Trans. on Power Delivery, Vol. 11, NO. 2, pp. 1 041-1 04 9. 44. T. M. Gruzs., ''A Survey of Neutral Curre11ts in Three-Phase Computer Power Systems." IEEE Trans. on Industry Applications, Vol. 26, NO. 4, pp. 719-725 July / August 1990. 45. A. Mansoor, W. M. Grady, A. H. Chowdhury and M. J. Samotyj., ' An Investigation of Harmonics Attenuation and Diversity Among Distributed SinglePhase Power Electronic loads." IEEE Paper 0-7803-1883-8, 1994. 46. Task Force on Harmonic Modeling and Simulation., ' Modeling and Simulation of the Propagation of harmonics in Electric power Networks." IEEE Trans on Power Delivery, Vol 11, NO. 1, pp. 466-474, January 1996. 4 7. ''IEEE Reconunended Practice for Establishing Transforrner Capability When Supplying Nonsinusoidal Load currents." ANSI / IEEE C57.110-1986. 48. Task Force on Harmonic Modeling and Simulation., '' Modeling and Simulation of the Propagation of harmonics in Electric power Networks." IEEE Trans. on Power Delivery, Vol. 11, NO. 1, pp. 452-465, January 1996. 49. P. T. Staats, W. M. Grady, A. Arapostathis and R. S. Tballam., ''A Statistical Method for Predicting the Net Harmonic Currents Generated by a Concentration of Electric Vehicle Battery Chargers." IEEE Trans. on Power Delivery, Vol. 12, NO. 3 pp 1258-1266, July 1997. 50. G. W. Massey., ' Power Distribution System Desig11 for Operation Under Nonsinusoidal Load Conditions." IEEE Trans. on Industry Applications, Vol. 31 NO. 3, pp. 513-519 May / June 1995. 51. A. Zahedi, J. A. Orr, A. E. Emanuel, D. J. Pillegi and T. j. Gentile., ''Comparative ha1monic flow Studies Using Current Harrnonic Injection Method." Proc Of IEEE ICHPS VI, Bologna, pp. 311-317, September 21-23, 1994. 52. S. Rahman and G. B. Shrestha., ''An Investigation into the Impact of Electric Vehicle Load on the Electric utility Distribution System IEEE Trans. on Power Delivery, Vol. 8, NO. 2, pp. 591-597, April 1993. 53. G. T. Heydt, D. J. Kish, F. Holcomb and Y Hill., ''A Methodology for Assessment of harmonic Impact and Compliance with Standards for Distribution Systems." IEEE Trans. on Power Delivery, Vol 6, NO. 4, pp. 1748-1754, October 1991. 54. K. Rajashekara., ''History of Electric Vehicles in General Motors.'' IEEE Trans. on Industry Applications, Vol. 30, NO. 4, pp. 897-904, July / August 1994.

PAGE 104

97 55 H. Shimizu, J Harada, C. Bland, K Kawakami and L. Chan., ''Advance Concepts in Elect r ic Vehicle Design." IEEE Trans. on Industrial Electronics Vol 44, NO. 1 pp. 14-18, February I 997 56. M. J. Riezenman. ''Electric Vehicle. IEEE Spectrum, pp 18-24 Nove1nber 1992. 57 Y Kim, and H Ha. ''Design of Interface Circuits W i th Electric Battery Model s." IEEE Trans. on Industrial Electronics Vol 44 NO. 1 pp 81-86, February 1997 58. J. Martinon P. Fauquembergue J. Lacl 1a ume and M. Meunier., ''A New Statistical Approach of harmonic Propagation in Transmission Sy s tems ." IEEE Trans. on Power Delivery Vol 11 NO. 2, pp 1032 1040, April 1996. 59 S. M. William s, G. T. Brownfield and J W. Duffu s ., ''Harmonic Propagation of an Electric Di st ribution System: Field Mea s urement s Compared with Computer Simulation." IEEE Tran s. on Power Delivery Vol 8, NO. 2, pp. 547-552, April 1993. 60. S. A Soliman ., '' A Simple and Novel Technique for Linear and Nonlinear Load Modeling in the Time Domain." International Journal of Power and Energy Sy ste ms, vol. 17, NO 2, pp 127-131, 1997 61 R. D. Hender so n and P. J. Ro se., '' Harmonic s: the Effects on Power Quality and Tran s former s." IEEE Trans on Industry Applications, Vol 30 NO. 3, pp 528-532 May / June 1994 62 G. W Ma ssey.,'' Estimation Methodes for Power Sy ste m Harmonic Effects on Power Di st ribution Tran sfo rmers. IEEE Trans on Indu s try Application s, Vol 30, NO. 2, pp. 485 489, March/April 1994. 63. Tan Weinian V. I. John, ''Nonlinear fluorescent systems: their impact on power quality." 1994 Canadian Conference on Electrical and Computer Engineering. Conference Pr ocee d i ngs (Cat. No.94TH8023 ), p. 2 vol. 834 144-7 vol. I Halifax NS, Canada; 25-28 Sept. 1994. 64. L. Cristaldi, A. Ferrero ''A method and related digital instrument for the measurement of the electric power qua l ity." IEEE Tran sac tions on Power Deli very, vol. I 0, no.3 p 1183-9. July 1995 65. A. E. Ema11uel '' On the as sess ment of harmonic polluti o n (of power systems)." I EEE Tran sac tions on Power Delivery vol.IO, n o .3 p 1693-8 July 1995. 66. L. L Grigsby, ''An overview of power quality." Proceedings of the IEEE International Confe rence on Indu s trial Technology p 863, 014 Guangzh o u China; 5 9 Dec. 1994

PAGE 105

98 67. E. P. Li, S. G. Chua, ''Harmonics measurement and analysis of power supply systems and n1odern electronic equipment." IPEC '95. Proceedings of the Int ernational Power Engineering Conference, p. 2 vo l xviii + 750, 517 22 vol.2. Singapore; 27 Feb.-1 March 1995 68. E. P. Li, S. G. Chua ''Har1nonics measurement and analysis of power s uppl y systems and modem electronic eq uipm ent." IPE C '95. Proceedings of tl1e Intern ational Power Engineering Conference, p. 2 vol. xviii+750 517-22 vol.2. Singapore; 27 Feb. -1 March 1 995. 69. D 0. Koval, C. Carter, ''Power q\tality characteristics of computer loads." IAS 95. Confe ren ce Record of the 1995 IEEE Industry Applications Conference. Thirtieth IAS Annual Meeting p. 3 vol. xxxvii+2765, 2265-72 vol.3 Orlando, FL USA; 8-12 Oct. 1995. 70. G. B. Rauch, D D. Sabin, M McGranaghan, A. Sundaram, 'Implementing system wide power quality monitoring plans." IEEE Transmission and Distribution Conference Proceedings. Proceedings of the 1 996 IEEE Power Engineering Society Transmission and Distribution Conference. p. xii + 706, 37-44. Los Angeles, CA USA; 15-20 Sept 1996. 71. G. T. Heydt E. Gunther, ''Post -m easurement processing of elect ri c power quality data." I EEE Transactions on Power D elivery, vol.1 1 no 4, p. 1853-9. Oct. 1996. 72. M. Lotfalian, ''Modeling and sim ulati on of power systems using PSPICE for power quality analysis. Proceedings of the IAST ED Internat ional Conference. Modeling and Simulation MS'94 p. 442, 164-6. Pittsburgh, PA USA; 2-4 May 1994. 73. D. 0. Koval, C. Carter, ''Power quality cl1aracteristics of computer loads." IEEE Transactions on Industry Applications, vol.33, no.3, p. 613-21. May-June 1997. 74. C. Muscas, ''Assessment of electric power quality: indices for identifying disturbing loads." Fourth International Workshop on Power Definitions and Measurements under Non-sinusoidal Condit i ons. 83-7. Milano Italy; 2-4 Sept 1997. 75. G T. Heydt, ' Electric power quality: a tutorial introduction." IEEE Computer Applications in Power vol. I 1 no.I, p. 1 5 -1 9 Jan. 1 998. 76. P. Bowler, ' Nonlinear loads on AC power systems." 32nd Universities Power Engineering Conference. UPEC '97, p. 2 vol. xxxix+ 1129, 838-42 vol.2.Manchester, UK; 10-12 Sept. 1997. 77 R. M. O'Connell, Y. J. Wang, G. Brownfield, ''Modeling and prediction of system voltage distortion l evels caused by nonlinear residential loads.'' Proceedings of the ICHQP. 7th International Conference on Harrnonics and Quality of Power, p x + 790, 6687 4. Las Vegas, NV USA; 16-18 Oct. 1996.

PAGE 106

99 78. R. Caldon, M. Fauri, L. Fellin, ''Harmonic modeling of household and co1nmercial aggregate loads Proceedings of the ICHQP. 7th International Conference on Harrnonics and Quality of Power, p. x+ 790, 681-6. Las Vegas, NV, USA; 16-18 Oct. 1996. 79. A. Domijan, Jr; D. Cheney, P. M. Garcia, ''Power quality characterization of an electric vehicle charger." Proceedings of the ICHQP 7th International Conference on Harmonics and Quality of Power, p. x+ 790, 721 6. Las Vegas, NV, USA; 16-18 Oct. 1996.

PAGE 107

BIOGRAPHICAL SKETCH Tariq Aslam Buchh was born in the beautiful valley of Kashmir. He received his BSSE from the Regional Engineering College, Srinagar India, in 1984. He completed his M.E. in Power Apparatus and Electric Drives from The University of Roorkee, India, in March 1990. For standing first in Power Apparatus and Electric Drives he was awarded the Uruversity Gold Medal by The University of Roorkee, India. Between 1984 and 1988 he worked as a lecturer at Regional Engineering College, Kashmir and later after completing his masters degree l1e rejoined the college as a senior lecturer. In 1992, he became a Research Associate in the Power Quality and Power Electronics Laboratory at the Department of Electrical and Computer Engineering, University of Florida. Gainesville Florida. I-Iis areas of interest include power electronics, microcomputer controlled electric drives, electric vehicles, PV array driven adjustable drives for Lunar Applications, harmonics in power systems. Mr. Buchh has published numerous journal and conference papers. He chaired a session at the International Association of Technology for Development (IASTED), International Conference on High Technology in Power System Industry. Oct 27-30 1997, Orlando Florida USA. Mr. Buchh is a student member IEEE. 100

PAGE 108

I certify that I have read this study and that in my opinion it conforn1s to acceptable standards of scholarly presentation and is fully adequate in scope and quality, a s a dis s ertation for the degree of Doctor of Philosophy. .> Alexander Domijan, r hairman Associate Profe s sor of Electrical and Computer Engineering I certify that I have read this study and that in my opinion it conforms to acceptable standards of scholarly presentation and is fully adequate in scope and quality, as a di s sertation for the degree of Doctor of Philosophy . f Dennis P Carroll Profe s sor of Electrical and Computer Engineering I certify that I have read this study and that in my opinion it confortns to acceptable standards of scholarly presentation and is fully adequate in scope and qt1ality, as a di s sertation for tl1e degree of Doctor of Philosophy. Khai D. T. Ngo Associate Professor of Electrical and Computer Engineering I certify that I have read this study and that in my opinion it conforms to acceptable standa1ds of scholarly presentation and is fully adequate, in scope and quality, as a dis s ertation for the degree of Doctor of sophy. Herman Lam Associate Professor of Electrical and Computer Engineering I certify that I have read this study and that in my opinion it conforms to acceptable standards of scholarly presentation and is fully adequate, in scope and quality, as a dissertation for the degree of Doctor of P osophy. ~t--:Barney L. Professor o Industrial and Sy ems Engineering

PAGE 109

This dissertation was submitted to the Graduate Faculty of the College of Engineering and to the Graduate school and was accepted as partial fulfillment of the requireme11ts for the degree of Doctor of Philo so phy August 1998 L---~-...._..__ ___ Winfred M. Phillip s Dean College of Engineering Karen A. Holbr oo k Dean Graduate School

PAGE 110

LD 1780 199 ..... CJ 19 UNIVERSITY OF FLORIDA II I II IIII I Ill I I II I II I I III II I II I I I I II II Il l lll l l 1 11 1111 1 11 1 3 1262 08554 9490