d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

Similar documents
Lecture 07: Poles and Zeros

Physics 114 Exam 2 Fall Name:

1 Lyapunov Stability Theory

Newton s Power Flow algorithm

Power Flow S + Buses with either or both Generator Load S G1 S G2 S G3 S D3 S D1 S D4 S D5. S Dk. Injection S G1

Mu Sequences/Series Solutions National Convention 2014

MOLECULAR VIBRATIONS

Econometric Methods. Review of Estimation

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

EECE 301 Signals & Systems

ENGI 4421 Propagation of Error Page 8-01

2/20/2013. Topics. Power Flow Part 1 Text: Power Transmission. Power Transmission. Power Transmission. Power Transmission

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.

Functions of Random Variables

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Multiple Choice Test. Chapter Adequacy of Models for Regression

8.1 Hashing Algorithms

TESTS BASED ON MAXIMUM LIKELIHOOD

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

Chapter 6. Small Signal Stability Analysis

CHAPTER 4 RADICAL EXPRESSIONS

Chapter 9 Jordan Block Matrices

[ L] υ = (3) [ L] n. Q: What are the units of K in Eq. (3)? (Why is units placed in quotations.) What is the relationship to K in Eq. (1)?

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

PRACTICAL CONSIDERATIONS IN HUMAN-INDUCED VIBRATION

Applying the condition for equilibrium to this equilibrium, we get (1) n i i =, r G and 5 i

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

1 Onto functions and bijections Applications to Counting

G S Power Flow Solution

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

Engineering Vibration 1. Introduction

5 Short Proofs of Simplified Stirling s Approximation

Non-degenerate Perturbation Theory

Bounds on the expected entropy and KL-divergence of sampled multinomial distributions. Brandon C. Roy

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

Numerical Analysis Formulae Booklet

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

Part 4b Asymptotic Results for MRR2 using PRESS. Recall that the PRESS statistic is a special type of cross validation procedure (see Allen (1971))

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

Multiple Linear Regression Analysis

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

ESS Line Fitting

We have already referred to a certain reaction, which takes place at high temperature after rich combustion.

Lecture Note to Rice Chapter 8

1. A real number x is represented approximately by , and we are told that the relative error is 0.1 %. What is x? Note: There are two answers.

x y exp λ'. x exp λ 2. x exp 1.

A Method for Damping Estimation Based On Least Square Fit

Department of Agricultural Economics. PhD Qualifier Examination. August 2011

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

1 Solution to Problem 6.40

MA 524 Homework 6 Solutions

Midterm Exam 1, section 1 (Solution) Thursday, February hour, 15 minutes

Chapter 4 Multiple Random Variables

STK4011 and STK9011 Autumn 2016

MEASURES OF DISPERSION

Objectives of Multiple Regression

ε. Therefore, the estimate

MA/CSSE 473 Day 27. Dynamic programming

CHAPTER VI Statistical Analysis of Experimental Data

2SLS Estimates ECON In this case, begin with the assumption that E[ i

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

CH E 374 Computational Methods in Engineering Fall 2007

Periodic Table of Elements. EE105 - Spring 2007 Microelectronic Devices and Circuits. The Diamond Structure. Electronic Properties of Silicon

Can we take the Mysticism Out of the Pearson Coefficient of Linear Correlation?

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

7.0 Equality Contraints: Lagrange Multipliers

Laboratory I.10 It All Adds Up

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Overview of the weighting constants and the points where we evaluate the function for The Gaussian quadrature Project two

L5 Polynomial / Spline Curves

Centroids & Moments of Inertia of Beam Sections

Third handout: On the Gini Index

Beam Warming Second-Order Upwind Method

D KL (P Q) := p i ln p i q i

Introduction to local (nonparametric) density estimation. methods

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits

Correlation and Regression Analysis

Chapter 14 Logistic Regression Models

4 Inner Product Spaces

Summary of the lecture in Biostatistics

For combinatorial problems we might need to generate all permutations, combinations, or subsets of a set.

Simulation Output Analysis

X ε ) = 0, or equivalently, lim

III-16 G. Brief Review of Grand Orthogonality Theorem and impact on Representations (Γ i ) l i = h n = number of irreducible representations.

Simple Linear Regression

ρ < 1 be five real numbers. The

Evaluating Polynomials

Lecture 5: Interpolation. Polynomial interpolation Rational approximation

Ahmed Elgamal. MDOF Systems & Modal Analysis

2. Independence and Bernoulli Trials

THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA

Midterm Exam 1, section 2 (Solution) Thursday, February hour, 15 minutes

The Selection Problem - Variable Size Decrease/Conquer (Practice with algorithm analysis)

STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1

Chapter 5 Properties of a Random Sample

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific

Transcription:

Learzato of the Swg Equato We wll cover sectos.5.-.6 ad begg of Secto 3.3 these otes. 1. Sgle mache-fte bus case Cosder a sgle mache coected to a fte bus, as show Fg. 1 below. E y1 V=1./_ Fg. 1 The admttace matrx s gve by Y Y Y 11 1 1 1 Y 1 Y y1 y (1) 1 Let s assume the mache s modeled by the swg equato wth dampg, gve by H d d m y e m y M s( ) where M E V Y1 Y 1 Y1 1 1 / (eables use of s stead of cos-see p.7 of A&F) Now let the agle δ chage by a small amout. The d d d d, (3) Also recall that by Taylor seres, ds x s x s x x s x x s x (cos x ) x (4) dx x () 1

The we also see that s( ) s( ) s( ) (cos( )) (5) (Eqt. 3.3) Applyg (3) to the left-had-sde of () ad (5) to the rght-hadsde of (), we obta H d d s( m M ) But Therefore, H d [s( ) (cos( )) ] m M m M s( ) M (cos( )) (6) m M s( ) d M (cos( )) Or H d d (cos( )) M (8) Now defe S M cos( ) (9) What s t? To aswer ths questo, observe: e M s( ) (1) de M cos( ) d (11) de d M cos( ) (1) (7)

Therefore, S d d e M cos( ) (13) S s called the sychrozg power coeffcet. I regards to early swg stablty (whch s a olear pheomea), the larger S s, the more stable wll be the geerator for a gve dsturbace. Ths s true because S dcates the slope of the power-agle curve, ad the hgher ths slope, the more deceleratg eergy s avalable to the mache for a gve fault. Ths dea s llustrated Fg.. Fg. But let s see what t meas for small sgal stablty, whch s characterzed by the egevalues (roots) of the system dfferetal equato trasformed to the s-doma through Lalace trasforms. Substtutg (13) to (8) results H d d S (14) Takg the Lalace trasform (assumg all tal coos are ), we obta 3

H s ( s) s ( s) 4 ( s) S (15) Elmatg δ(s), we obta the system s characterstc equato: H s s S (16) Solvg usg the quadratc formula, we get s S (Eqt. 3.7) 1 H (17) ullg ω/h out of the radcal, we have 8 H s S (18) (Eqt. 3.8) We ca make some observatos about (18), as follows: 1. No dampg: If =, the H s 8 S S H, (19a) or s S H (19b) a. Observe (19b) that f S>, (a) ay respose to a small dsturbace wll be oscllatory, ad (b) the oscllatory frequecy becomes lower as H becomes larger. b. Observe (19a) that f S<, the s S H () ad ay respose s ustable.

Fgures 3, 4 llustrate, for both stuatos S>, S<, respectvely, the pole (egevalue) locatos the s-plae ad the operatg pot locato o the power-agle curve. Im{s} X e S> {s} X Im{s} Fg. 3: S> e δ S> δ X X {s} δ δ Fg. 4: S< I Fg. 3, the oscllatory system s characterzed by purely magary poles (left) ad a stable operatg pot (rght). I Fg. 4, the ustable system s characterzed by the RH-pole (left) ad a ustable equlbrum pot (rght). 5

. Wth dampg: If, the 8 H s S (18) Let s look at the most postve root (ad so we wll use + sg before the radcal, ad we esure the cotrbuto from the secod term sde the radcal s postve,.e., S<) ad ask what are the coos uder whch t ca be the rght-half-plae, that s: s 6 8 S H 8 H S 8 H? S 8 H S 8? S H??? The above relato must be true. Because the above relato s depedet of dampg, we coclude that f S<, the system must be ustable, depedet of how much dampg exsts.. Mult-mache case (Secto 3.4) (We wll come back to sectos 3. ad 3.3.1) call that for a geerator coected to a fte bus, we foud that the swg equato s

H d d m e (1) where e M s( ) M E V Y1 Y 1 Y1 1 Lettg ad learzg, we fd that H d d S where d e S M cos( ) d Let s ow cosder the mult-mache system assumg: Classcal models Network reduced to oly teral geerator odes For geerator, we have that the swg equato s H d d m e where e E G EE Y cos( ) 1 E G 1 Y cos( ) () (13) (3) (4) 7

where δ=δ-δ. I (4), all voltages E, E, ad all Y-bus elemets Y are magtudes. Now let s cosder a small chage the agle of mache : δ=δ+ δ. The left-had-sde of (4) s precsely as the case of the sgle geerator vs. fte by s case. But what happeed to the rght-hadsde? Now the rght-had-sde s, by (3), m e m s uaffected by + δ, but e s affected by t. call δ=δ-δ. We cosder a small chage rotor agle at geerator. To be more geeral, we also allow a small chage geerator. However, geerator wll ot chage as a result of the geerator chage; they are are depedet chages ad we could ust as well have oly oe of them. δ=δ+ δ δ=δ+ δ callg that e E G EE Y cos( ) 1 (5) we eed to see what happes to the cos term for the small chage agle. We kow from trgoometry that cos( x y) s xs y cos x cos The cos( ) s s cos cos (6) Applcato of (6) to (5) yelds:. y 8

e E E G G 1 1 Y B s s Y s G cos cos cos (7) (Eqt. 3.1) Now we eed to learze the cosδ ad sδ terms usg δ=δ+ δ. From Taylor seres wth frst order term oly, s s( ) s cos (8) cos cos( ) cos s (9) Substtutg (8) ad (9) to (7), we get e E G EE B s cos 1 Now collect terms δ: G cos s (3) 9

e E G 1 1 B B B s cos cos m 1 (31a) call that the rght-had-sde of the swg equato s m-e. Equato (31a) ca be rewrtte the as (31b) 1 G G G e m 1 therefore the swg equato (3), whch s H d d m e becomes H d d B cos s s cos G s 1 B cos G s (31b) (3) (3) efe everythg sde the expresso wth the summato of (3), except δ, as S, that s

S EE B cos G s (33) The (3) becomes H d d S (34) Gve the mechacal power s costat, the rght-had-sde of (34) gves the egatve of the chage electrc power out of the mache due to the small chages δ, that s e S 1 11 1 (35) (Eqt. 3.3) What s S? We aswer ths questo by observg that the power flowg from geerator teral ode to geerator teral ode s EE B s G cos (36) fferetatg, we get EE B cos G s (37) Evaluatg at δ, we get B cos G s S (38) (Eqt. 3.4) Note that f bus s the fte bus, eglectg resstace, we have: B S cos whch s the same as the sychrozg power coeffcet the fte bus case (we called t S). We wll look at multmache systems, but before we do, we cosder somethg very mportat respose to load chages...

Oe last ssue: what s the dfferece betwee sychrozg power coeffcet, geerato shft factor (GSF) ad power trasfer dstrbuto factor? We aswer ths here. Sychrozg power coeffcet (SC): S EE B cos G s (a1) Observe that the SC gves o chage flow o crcut {,} wth respect to o a chage agular separato across {,}. Geerato shft factor (GSF): { b} t{ b}, (b1) allocato olcy Observe that the GSF gves o chage flow across ay brach b wth respect to o a chage ecto at bus, subect to a reallocato polcy (.e., how the bus chage ecto s compesated). ower trasfer dstrbuto factor (TF) for 1-bus ecto chage: {} b TF { b}, (c) allocato olcy The TF for 1-bus ecto chage s the same as the GSF. ower trasfer dstrbuto factor for -bus ecto chage: TF { },, { } { } (c) allocato olcy allocato olcy The upshot of the above s that relatg SC to GSF s eough to relate SC to TF. We relate SC to GSF as follows: 1

From (a1), we wrte that S S (a) From (b1), we wrte that {} b t{ b}, { b} t{ b}, allocato olcy allocato olcy allocato olcy (b) Now 1. Cosder our power system s experecg coos such that the agular separato betwee buses ad s δ.. Le b s termated by buses ad,.e., b {,}. 3. We make a chage ected power at bus equal to Δ compesated by a reallocato polcy where a equal ad opposte chage, Δ, s made at bus. The (a) ad (b) are equvalet: t allocato S { } { }, allocato b b olcy: olcy: That s: t S { }, allocato b olcy: whch shows us that allocato olcy: t { b}, allocato S olcy: allocato olcy: 13