Part A: From (14.19) in your text (see notes Stability1 ), we have

Size: px
Start display at page:

Download "Part A: From (14.19) in your text (see notes Stability1 ), we have"

Transcription

1 Prob Consider a round-rotor generator delivering a steady-state power P G =.5 to an infinite bus through a transmission line with reactance X L =.4. Assume that E a =1.8, V =1/_, H=5sec, and X d =X q =1.. Neglect resistance. A.Find the two possible steady-state (equilibrium values of power angle δ=/_e a -/_V that lie in the interval [,π]. B. Which of these two equilibrium values are we likely to observe in practice? Hint: Linearize around the equilibrium value and consider the natural frequencies found from the linear differential equation in Δδ. Part A: From (14.19 in your text (see notes Stability1, we have M ( t D ( t P G ( PM The problem implies that P M =.5, therefore we have M D PG (.5 Here: Ea V 1.8*1. PG ( sin( sin( 1.86 sin( X X 1..4 d L 1.86 sin( sin

2 Part B: What does it mean to linearize around the equilibrium value? It means to find a linear approximation to the function at that value. If the function is a nonlinear function, then the behavior will only be valid for small distances away from the value. The figure below illustrates: If our function is in terms of δ, then our linearized function is in terms of Δδ. Wherever there is a linear terms in δ of our original function, we may replace it directly with Δδ. But if there is a nonlinear term in our original function, then we must use a Taylor series expansion. Our function is M ( t D ( t P G ( PM M D 1.86 sin. 5 (1

3 Now let the angle δ change by a small amount about an equilibrium δ. Then d d dt dt, d d dt dt 3 ( Recall the Taylor series expansion of a function is: f ( x ( x ( f ( x ( f' ( x ( x ( 1 f'' ( x ( ( x and recognizing that for small Δx, the higher-order terms are extremely small and may be neglected, we obtain, ( ( ( ( d sin sin sin d ( f ( x x f ( x f'( x x And so the sin function can be expanded to (3 sin sin (cos Substituting (3 and ( into (1, we obtain: M D 1.86(sin (cos Distributing through the 1.86, we get ( M D 1.86 sin 1.86 cos....5 And bringing the 1.86 sinδ to the right-hand-side, we obtain: M D 1.86 cos sin (4 Now, consider that δ is either.89 or Either way, 1.86 sin Therefore,.5 M D 1.86 cos.5

4 This is the result of Linearization of swing equation given on pg The LaPlace Transform of this equation is: M ( s s D ( s s 1.86 cos ( s Factoring out the Δδ(s, we have ( s Ms Ds 1.86 cos (5 Now what if the right-hand-side of (4 is not zero but rather, what if we perturb this system ever so slightly so that 1.86 sin.5 In this case, u( t Ms Ds 1.86 cos U( s ( s and we have ( s Ms Ds U( s cos We observe that the denominator of the previous equation is the characteristic polynomial of this system. When the system is in equilibrium, we may divide both sides of (5 by Δδ(s to obtain the characteristic equation of the system:

5 o Ms Ds 1.86 cos( H 5 M.65 f o 6 THEN.65s Ds 1.86 cos( o If we plug in δ o =.89 degrees into the above equation, we obtain:.65s.65s This has roots of D D D Ds 1.86 cos(.89 Ds.5 4(.65 (.5 (.65 D D.53 D Consider D is very small, e.g., D=.1. Then the above is j3.915 In this case, both roots are in the Left-Hand-Plane. It should be clear that any smaller value of D will not change this situation at all (since the first term will still be real and negative and the second term will still be imaginary. Now consider that D is large, e.g., D=1. Then the above is 5

6 , In this case, both roots are in the Left-Hand-Plane. It should be clear that any larger value of D will not change this situation at all (since the first term will still be large and negative, and the second term will never be greater than the first term in absolute value. A system having all of the roots of the characteristic equation in the Left-Half-Plane is a stable system. This corresponds to our understanding of the stable equilibrium point that we discussed in class. If on the other hand we make δ o = degrees one coefficient is negative, meaning that the natural frequencies are in the right hand plane. If we plug in δ o =.89 degrees into the above equation, we obtain:.65s.65s This has roots of D D D Ds 1.86 cos( Ds (.65 ( (.65 D D.53 D It should be clear that the second term will be real and it will be larger in magnitude than the first term, 6

7 independent of what D is. For example, if D=1, then we have , Thus, there will be one root of the characteristic equation in the Right-Half-Plane. A system having one or more roots in the Right-Half-Plane in an unstable system. This corresponds to our understanding of the unstable equilibrium point that we discussed in class. To answer the question, we will never observe the unstable equilibrium in practice because small perturbations are always happening in the power system and would perturb the system off of this point, and once perturbed, it would either move to the stable equilibrium or it would go unstable. The stable equilibrium is the only one we will ever observe in practice. 7

SCHOOL OF ELECTRICAL, MECHANICAL AND MECHATRONIC SYSTEMS. Transient Stability LECTURE NOTES SPRING SEMESTER, 2008

SCHOOL OF ELECTRICAL, MECHANICAL AND MECHATRONIC SYSTEMS. Transient Stability LECTURE NOTES SPRING SEMESTER, 2008 SCHOOL OF ELECTRICAL, MECHANICAL AND MECHATRONIC SYSTEMS LECTURE NOTES Transient Stability SPRING SEMESTER, 008 October 7, 008 Transient Stability Transient stability refers to the ability of a synchronous

More information

( ) ( = ) = ( ) ( ) ( )

( ) ( = ) = ( ) ( ) ( ) ( ) Vρ C st s T t 0 wc Ti s T s Q s (8) K T ( s) Q ( s) + Ti ( s) (0) τs+ τs+ V ρ K and τ wc w T (s)g (s)q (s) + G (s)t(s) i G and G are transfer functions and independent of the inputs, Q and T i. Note

More information

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism

More information

Power Angle & Time Domain Plots 1.0 The power-angle curves

Power Angle & Time Domain Plots 1.0 The power-angle curves Power Angle & Time Domain Plots 1.0 The power-angle curves Consider a generator connected to an infinite bus through a network, as shown in Fig. 1 below (all in pu). E E NETWORK X V V 0 Fig. 1 We assume

More information

Now multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:

Now multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get: Equal Area Criterion.0 Developent of equal area criterion As in previous notes, all powers are in per-unit. I want to show you the equal area criterion a little differently than the book does it. Let s

More information

Laplace Transforms Chapter 3

Laplace Transforms Chapter 3 Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important

More information

EE 451 Power System Stability

EE 451 Power System Stability EE 451 Power System Stability Power system operates in synchronous mode Power system is subjected to a wide range of disturbances (small and large) - Loads and generation changes - Network changes - Faults

More information

Transient Stability Analysis of Single Machine Infinite Bus System by Numerical Methods

Transient Stability Analysis of Single Machine Infinite Bus System by Numerical Methods International Journal of Electrical and Electronics Research ISSN 348-6988 (online) Vol., Issue 3, pp: (58-66), Month: July - September 04, Available at: www.researchpublish.com Transient Stability Analysis

More information

Section 7.4: Inverse Laplace Transform

Section 7.4: Inverse Laplace Transform Section 74: Inverse Laplace Transform A natural question to ask about any function is whether it has an inverse function We now ask this question about the Laplace transform: given a function F (s), will

More information

There are two main properties that we use when solving linear equations. Property #1: Additive Property of Equality

There are two main properties that we use when solving linear equations. Property #1: Additive Property of Equality Chapter 1.1: Solving Linear and Literal Equations Linear Equations Linear equations are equations of the form ax + b = c, where a, b and c are constants, and a zero. A hint that an equation is linear is

More information

APPM 2360: Midterm 3 July 12, 2013.

APPM 2360: Midterm 3 July 12, 2013. APPM 2360: Midterm 3 July 12, 2013. ON THE FRONT OF YOUR BLUEBOOK write: (1) your name, (2) your instructor s name, (3) your recitation section number and (4) a grading table. Text books, class notes,

More information

Automatic Control Systems theory overview (discrete time systems)

Automatic Control Systems theory overview (discrete time systems) Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations

More information

LIMITS AT INFINITY MR. VELAZQUEZ AP CALCULUS

LIMITS AT INFINITY MR. VELAZQUEZ AP CALCULUS LIMITS AT INFINITY MR. VELAZQUEZ AP CALCULUS RECALL: VERTICAL ASYMPTOTES Remember that for a rational function, vertical asymptotes occur at values of x = a which have infinite its (either positive or

More information

The Operation of a Generator on Infinite Busbars

The Operation of a Generator on Infinite Busbars The Operation of a Generator on Infinite Busbars In order to simplify the ideas as much as possible the resistance of the generator will be neglected; in practice this assumption is usually reasonable.

More information

Autonomous Systems and Stability

Autonomous Systems and Stability LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an

More information

ζ (s) = s 1 s {u} [u] ζ (s) = s 0 u 1+sdu, {u} Note how the integral runs from 0 and not 1.

ζ (s) = s 1 s {u} [u] ζ (s) = s 0 u 1+sdu, {u} Note how the integral runs from 0 and not 1. Problem Sheet 3. From Theorem 3. we have ζ (s) = + s s {u} u+sdu, (45) valid for Res > 0. i) Deduce that for Res >. [u] ζ (s) = s u +sdu ote the integral contains [u] in place of {u}. ii) Deduce that for

More information

Learn how to use Desmos

Learn how to use Desmos Learn how to use Desmos Maclaurin and Taylor series 1 Go to www.desmos.com. Create an account (click on bottom near top right of screen) Change the grid settings (click on the spanner) to 1 x 3, 1 y 12

More information

CHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2

CHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2 CHEE 39 Tutorial 3 Solutions. Using partial fraction expansions, find the causal function f whose Laplace transform is given by: F (s) 0 f(t)e st dt (.) F (s) = s(s+) ; Solution: Note that the polynomial

More information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi System Stability - 26 March, 2014 Prof. Dr. Eleni Chatzi System Stability - 26 March, 24 Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can

More information

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace

Unit 2: Modeling in the Frequency Domain Part 2: The Laplace Transform. The Laplace Transform. The need for Laplace Unit : Modeling in the Frequency Domain Part : Engineering 81: Control Systems I Faculty of Engineering & Applied Science Memorial University of Newfoundland January 1, 010 1 Pair Table Unit, Part : Unit,

More information

Basic Procedures for Common Problems

Basic Procedures for Common Problems Basic Procedures for Common Problems ECHE 550, Fall 2002 Steady State Multivariable Modeling and Control 1 Determine what variables are available to manipulate (inputs, u) and what variables are available

More information

QUESTION BANK ENGINEERS ACADEMY. Power Systems Power System Stability 1

QUESTION BANK ENGINEERS ACADEMY. Power Systems Power System Stability 1 ower ystems ower ystem tability QUETION BANK. A cylindrical rotor generator delivers 0.5 pu power in the steady-state to an infinite bus through a transmission line of reactance 0.5 pu. The generator no-load

More information

6.1 Polynomial Functions

6.1 Polynomial Functions 6.1 Polynomial Functions Definition. A polynomial function is any function p(x) of the form p(x) = p n x n + p n 1 x n 1 + + p 2 x 2 + p 1 x + p 0 where all of the exponents are non-negative integers and

More information

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27

ECEN 605 LINEAR SYSTEMS. Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability 1/27 1/27 ECEN 605 LINEAR SYSTEMS Lecture 20 Characteristics of Feedback Control Systems II Feedback and Stability Feedback System Consider the feedback system u + G ol (s) y Figure 1: A unity feedback system

More information

Math 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt

Math 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning

More information

Weierstrass Products and More on Analytic Continuation

Weierstrass Products and More on Analytic Continuation Weierstrass Products and More on Analytic Continuation Thursday, December 05, 2013 1:58 PM Hard deadline for Homework 4: Wednesday, December 11 at 5 PM. rapid grading on or after Friday, December 6 Final

More information

1. There are 8 questions spanning 9 pages total (including this cover page). Please make sure that you have all 9 pages before starting.

1. There are 8 questions spanning 9 pages total (including this cover page). Please make sure that you have all 9 pages before starting. Instructor: K. Rotz Name: Solution PUID: 00000-00000 Instructions and tips: 1. There are 8 questions spanning 9 pages total (including this cover page). Please make sure that you have all 9 pages before

More information

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks. Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the

More information

Math 216 Second Midterm 16 November, 2017

Math 216 Second Midterm 16 November, 2017 Math 216 Second Midterm 16 November, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material

More information

The synchronous machine (SM) in the power system (2) (Where does the electricity come from)?

The synchronous machine (SM) in the power system (2) (Where does the electricity come from)? 1 The synchronous machine (SM) in the power system (2) (Where does the electricity come from)? 2 Lecture overview Synchronous machines with more than 2 magnetic poles The relation between the number of

More information

1. Introduction. Keywords Transient Stability Analysis, Power System, Swing Equation, Three-Phase Fault, Fault Clearing Time

1. Introduction. Keywords Transient Stability Analysis, Power System, Swing Equation, Three-Phase Fault, Fault Clearing Time Energy and Power 17, 7(1): -36 DOI: 1.593/j.ep.1771.3 Numerical Simulations for Transient Stability Analysis of Two-Machine Power System Considering Three-Phase Fault under Different Fault Clearing Times

More information

Our goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always

Our goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.

More information

Problem Set 9 Solutions

Problem Set 9 Solutions 8.4 Problem Set 9 Solutions Total: 4 points Problem : Integrate (a) (b) d. ( 4 + 4)( 4 + 5) d 4. Solution (4 points) (a) We use the method of partial fractions to write A B (C + D) = + +. ( ) ( 4 + 5)

More information

POWER SYSTEM STABILITY

POWER SYSTEM STABILITY LESSON SUMMARY-1:- POWER SYSTEM STABILITY 1. Introduction 2. Classification of Power System Stability 3. Dynamic Equation of Synchronous Machine Power system stability involves the study of the dynamics

More information

Chapter 1. Introduction to Nonlinear Space Plasma Physics

Chapter 1. Introduction to Nonlinear Space Plasma Physics Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large

More information

Engineering Tripos Part IIB Nonlinear Systems and Control. Handout 3: Describing Functions

Engineering Tripos Part IIB Nonlinear Systems and Control. Handout 3: Describing Functions Engineering Tripos Part IIB Module 4F2 Nonlinear Systems and Control Handout 3: Describing Functions 1 Harmonic balance The describing function method (also called the method of harmonic balance) uses

More information

Damped Harmonic Oscillator

Damped Harmonic Oscillator Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or

More information

Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs

Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs 3 Analytical Solutions of Linear ODE-IVPs Before developing numerical

More information

MATH 31B: MIDTERM 2 REVIEW. sin 2 x = 1 cos(2x) dx = x 2 sin(2x) 4. + C = x 2. dx = x sin(2x) + C = x sin x cos x

MATH 31B: MIDTERM 2 REVIEW. sin 2 x = 1 cos(2x) dx = x 2 sin(2x) 4. + C = x 2. dx = x sin(2x) + C = x sin x cos x MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate sin x and cos x. Solution: Recall the identities cos x = + cos(x) Using these formulas gives cos(x) sin x =. Trigonometric Integrals = x sin(x) sin x = cos(x)

More information

Comparative Study of Synchronous Machine, Model 1.0 and Model 1.1 in Transient Stability Studies with and without PSS

Comparative Study of Synchronous Machine, Model 1.0 and Model 1.1 in Transient Stability Studies with and without PSS Comparative Study of Synchronous Machine, Model 1.0 and Model 1.1 in Transient Stability Studies with and without PSS Abhijit N Morab, Abhishek P Jinde, Jayakrishna Narra, Omkar Kokane Guide: Kiran R Patil

More information

Math 4263 Homework Set 1

Math 4263 Homework Set 1 Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that

More information

Time Response of Systems

Time Response of Systems Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =

More information

Chemical Engineering 436 Laplace Transforms (1)

Chemical Engineering 436 Laplace Transforms (1) Chemical Engineering 436 Laplace Transforms () Why Laplace Transforms?? ) Converts differential equations to algebraic equations- facilitates combination of multiple components in a system to get the total

More information

Laplace Transform Part 1: Introduction (I&N Chap 13)

Laplace Transform Part 1: Introduction (I&N Chap 13) Laplace Transform Part 1: Introduction (I&N Chap 13) Definition of the L.T. L.T. of Singularity Functions L.T. Pairs Properties of the L.T. Inverse L.T. Convolution IVT(initial value theorem) & FVT (final

More information

JUST THE MATHS UNIT NUMBER LAPLACE TRANSFORMS 3 (Differential equations) A.J.Hobson

JUST THE MATHS UNIT NUMBER LAPLACE TRANSFORMS 3 (Differential equations) A.J.Hobson JUST THE MATHS UNIT NUMBER 16.3 LAPLACE TRANSFORMS 3 (Differential equations) by A.J.Hobson 16.3.1 Examples of solving differential equations 16.3.2 The general solution of a differential equation 16.3.3

More information

Antiderivatives. Definition A function, F, is said to be an antiderivative of a function, f, on an interval, I, if. F x f x for all x I.

Antiderivatives. Definition A function, F, is said to be an antiderivative of a function, f, on an interval, I, if. F x f x for all x I. Antiderivatives Definition A function, F, is said to be an antiderivative of a function, f, on an interval, I, if F x f x for all x I. Theorem If F is an antiderivative of f on I, then every function of

More information

Transient Response of a Second-Order System

Transient Response of a Second-Order System Transient Response of a Second-Order System ECEN 830 Spring 01 1. Introduction In connection with this experiment, you are selecting the gains in your feedback loop to obtain a well-behaved closed-loop

More information

Symmetrical Components Fall 2007

Symmetrical Components Fall 2007 0.1 Variables STEADYSTATE ANALYSIS OF SALIENT-POLESYNCHRONOUS GENERATORS This paper is intended to provide a procedure for calculating the internal voltage of a salientpole synchronous generator given

More information

Control Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich

Control Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017

More information

Handout 2: Invariant Sets and Stability

Handout 2: Invariant Sets and Stability Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state

More information

MATH 18.01, FALL PROBLEM SET #5 SOLUTIONS (PART II)

MATH 18.01, FALL PROBLEM SET #5 SOLUTIONS (PART II) MATH 8, FALL 7 - PROBLEM SET #5 SOLUTIONS (PART II (Oct ; Antiderivatives; + + 3 7 points Recall that in pset 3A, you showed that (d/dx tanh x x Here, tanh (x denotes the inverse to the hyperbolic tangent

More information

Lecture 5: Function Approximation: Taylor Series

Lecture 5: Function Approximation: Taylor Series 1 / 10 Lecture 5: Function Approximation: Taylor Series MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Better

More information

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Laplace Transforms Dr. M. A. A. Shoukat Choudhury 1 Laplace Transforms Important analytical

More information

Power System Stability and Control. Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India

Power System Stability and Control. Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India Power System Stability and Control Dr. B. Kalyan Kumar, Department of Electrical Engineering, Indian Institute of Technology Madras, Chennai, India Contents Chapter 1 Introduction to Power System Stability

More information

POWER SERIES REVIEW SOLUTIONS

POWER SERIES REVIEW SOLUTIONS POWER SERIES REVIEW SOLUTIONS 1. Convergence of power series: For the following, find the radius of convergence: a) (m + 1)mx m In CME 10, we only teach you the ratio test, so that is the only test you

More information

FINAL EXAM SOLUTIONS, MATH 123

FINAL EXAM SOLUTIONS, MATH 123 FINAL EXAM SOLUTIONS, MATH 23. Find the eigenvalues of the matrix ( 9 4 3 ) So λ = or 6. = λ 9 4 3 λ = ( λ)( 3 λ) + 36 = λ 2 7λ + 6 = (λ 6)(λ ) 2. Compute the matrix inverse: ( ) 3 3 = 3 4 ( 4/3 ) 3. Let

More information

Chapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x

Chapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x Use of Tools from Interactive Differential Equations with the texts Fundamentals of Differential Equations, 5th edition and Fundamentals of Differential Equations and Boundary Value Problems, 3rd edition

More information

Completion Date: Monday February 11, 2008

Completion Date: Monday February 11, 2008 MATH 4 (R) Winter 8 Intermediate Calculus I Solutions to Problem Set #4 Completion Date: Monday February, 8 Department of Mathematical and Statistical Sciences University of Alberta Question. [Sec..9,

More information

DESIGNING POWER SYSTEM STABILIZER WITH PID CONTROLLER

DESIGNING POWER SYSTEM STABILIZER WITH PID CONTROLLER International Journal on Technical and Physical Problems of Engineering (IJTPE) Published by International Organization on TPE (IOTPE) ISSN 2077-3528 IJTPE Journal www.iotpe.com ijtpe@iotpe.com June 2010

More information

July 21 Math 2254 sec 001 Summer 2015

July 21 Math 2254 sec 001 Summer 2015 July 21 Math 2254 sec 001 Summer 2015 Section 8.8: Power Series Theorem: Let a n (x c) n have positive radius of convergence R, and let the function f be defined by this power series f (x) = a n (x c)

More information

Aldi Mucka Tonin Dodani Marjela Qemali Rajmonda Bualoti Polytechnic University of Tirana, Faculty of Electric Engineering, Republic of Albania

Aldi Mucka Tonin Dodani Marjela Qemali Rajmonda Bualoti Polytechnic University of Tirana, Faculty of Electric Engineering, Republic of Albania 8. СОВЕТУВАЊЕ Охрид, 22 24 септември Aldi Mucka Tonin Dodani Marjela Qemali Rajmonda Bualoti Polytechnic University of Tirana, Faculty of Electric Engineering, Republic of Albania REHABILITATIONS OF EXCITATION

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems

Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now

More information

Due Wednesday, February 6th EE/MFS 599 HW #5

Due Wednesday, February 6th EE/MFS 599 HW #5 Due Wednesday, February 6th EE/MFS 599 HW #5 You may use Matlab/Simulink wherever applicable. Consider the standard, unity-feedback closed loop control system shown below where G(s) = /[s q (s+)(s+9)]

More information

Lecture 38. Almost Linear Systems

Lecture 38. Almost Linear Systems Math 245 - Mathematics of Physics and Engineering I Lecture 38. Almost Linear Systems April 20, 2012 Konstantin Zuev (USC) Math 245, Lecture 38 April 20, 2012 1 / 11 Agenda Stability Properties of Linear

More information

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response

Raktim Bhattacharya. . AERO 422: Active Controls for Aerospace Vehicles. Dynamic Response .. AERO 422: Active Controls for Aerospace Vehicles Dynamic Response Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. . Previous Class...........

More information

Nonlinear System Analysis

Nonlinear System Analysis Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control

More information

ECE 585 Power System Stability

ECE 585 Power System Stability Homework 1, Due on January 29 ECE 585 Power System Stability Consider the power system below. The network frequency is 60 Hz. At the pre-fault steady state (a) the power generated by the machine is 400

More information

Hermite Interpolation

Hermite Interpolation Jim Lambers MAT 77 Fall Semester 010-11 Lecture Notes These notes correspond to Sections 4 and 5 in the text Hermite Interpolation Suppose that the interpolation points are perturbed so that two neighboring

More information

Simple Harmonic Motion Concept Questions

Simple Harmonic Motion Concept Questions Simple Harmonic Motion Concept Questions Question 1 Which of the following functions x(t) has a second derivative which is proportional to the negative of the function d x! " x? dt 1 1. x( t ) = at. x(

More information

Math 273 (51) - Final

Math 273 (51) - Final Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all

More information

STABILITY ANALYSIS TECHNIQUES

STABILITY ANALYSIS TECHNIQUES ECE4540/5540: Digital Control Systems 4 1 STABILITY ANALYSIS TECHNIQUES 41: Bilinear transformation Three main aspects to control-system design: 1 Stability, 2 Steady-state response, 3 Transient response

More information

Identification Methods for Structural Systems

Identification Methods for Structural Systems Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from

More information

APPM 2360 Section Exam 3 Wednesday November 19, 7:00pm 8:30pm, 2014

APPM 2360 Section Exam 3 Wednesday November 19, 7:00pm 8:30pm, 2014 APPM 2360 Section Exam 3 Wednesday November 9, 7:00pm 8:30pm, 204 ON THE FRONT OF YOUR BLUEBOOK write: () your name, (2) your student ID number, (3) lecture section, (4) your instructor s name, and (5)

More information

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems. Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics

More information

Midterm Examination #2

Midterm Examination #2 Anthony Austin MATH 47 April, 9 AMDG Midterm Examination #. The integrand has poles of order whenever 4, which occurs when is equal to i, i,, or. Since a R, a >, the only one of these poles that lies inside

More information

MATH 3150: PDE FOR ENGINEERS FINAL EXAM (VERSION D) 1. Consider the heat equation in a wire whose diffusivity varies over time: u k(t) 2 x 2

MATH 3150: PDE FOR ENGINEERS FINAL EXAM (VERSION D) 1. Consider the heat equation in a wire whose diffusivity varies over time: u k(t) 2 x 2 MATH 35: PDE FOR ENGINEERS FINAL EXAM (VERSION D). Consider the heat equation in a wire whose diffusivity varies over time: u t = u k(t) x where k(t) is some positive function of time. Assume the wire

More information

The Harmonic Oscillator

The Harmonic Oscillator The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can

More information

max clear Figure 1 below illustrates a stable case. P post (f) (g) (e) A 2 (d) (c)

max clear Figure 1 below illustrates a stable case. P post (f) (g) (e) A 2 (d) (c) Stability 4. Introduction In the previous notes (Stability 3, we developed the equal area iterion, which says that For stability, A =A, which means the decelerating energy (A must equal the accelerating

More information

Examples of Applications of Potential Functions in Problem Solving (Web Appendix to the Paper)

Examples of Applications of Potential Functions in Problem Solving (Web Appendix to the Paper) Examples of Applications of otential Functions in roblem Solving (Web Appendix to the aper) Ali Mehrizi-Sani and Reza Iravani May 5, 2010 1 Introduction otential functions may be exploited to formulate

More information

θ α W Description of aero.m

θ α W Description of aero.m Description of aero.m Determination of the aerodynamic forces, moments and power by means of the blade element method; for known mean wind speed, induction factor etc. Simplifications: uniform flow (i.e.

More information

AP Calculus Chapter 9: Infinite Series

AP Calculus Chapter 9: Infinite Series AP Calculus Chapter 9: Infinite Series 9. Sequences a, a 2, a 3, a 4, a 5,... Sequence: A function whose domain is the set of positive integers n = 2 3 4 a n = a a 2 a 3 a 4 terms of the sequence Begin

More information

Control Systems. Laplace domain analysis

Control Systems. Laplace domain analysis Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output

More information

Stability of Feedback Control Systems: Absolute and Relative

Stability of Feedback Control Systems: Absolute and Relative Stability of Feedback Control Systems: Absolute and Relative Dr. Kevin Craig Greenheck Chair in Engineering Design & Professor of Mechanical Engineering Marquette University Stability: Absolute and Relative

More information

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions

MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions MAE143a: Signals & Systems (& Control) Final Exam (2011) solutions Question 1. SIGNALS: Design of a noise-cancelling headphone system. 1a. Based on the low-pass filter given, design a high-pass filter,

More information

The Effects of Machine Components on Bifurcation and Chaos as Applied to Multimachine System

The Effects of Machine Components on Bifurcation and Chaos as Applied to Multimachine System 1 The Effects of Machine Components on Bifurcation and Chaos as Applied to Multimachine System M. M. Alomari and B. S. Rodanski University of Technology, Sydney (UTS) P.O. Box 123, Broadway NSW 2007, Australia

More information

STABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse

STABILITY. Have looked at modeling dynamic systems using differential equations. and used the Laplace transform to help find step and impulse SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 4. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.sigmedia.tv STABILITY Have looked at modeling dynamic systems using differential

More information

Modeling and Analysis of Dynamic Systems

Modeling and Analysis of Dynamic Systems Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability

More information

Cylinder with Conical Cap

Cylinder with Conical Cap Cylinder with Conical Cap MUS420/EE367A Supplement to Lecture 11 Stability Proof for a Cylindrical Bore with Conical Cap Julius O. Smith III (jos@ccrma.stanford.edu) Center for Computer Research in Music

More information

ECE317 : Feedback and Control

ECE317 : Feedback and Control ECE317 : Feedback and Control Lecture : Stability Routh-Hurwitz stability criterion Dr. Richard Tymerski Dept. of Electrical and Computer Engineering Portland State University 1 Course roadmap Modeling

More information

Math Real Analysis II

Math Real Analysis II Math 432 - Real Analysis II Solutions to Homework due February 3 In class, we learned that the n-th remainder for a smooth function f(x) defined on some open interval containing is given by f (k) () R

More information

Solutions to Math 41 Exam 2 November 10, 2011

Solutions to Math 41 Exam 2 November 10, 2011 Solutions to Math 41 Eam November 10, 011 1. (1 points) Find each of the following its, with justification. If the it does not eist, eplain why. If there is an infinite it, then eplain whether it is or.

More information

Name: AK-Nummer: Ergänzungsprüfung January 29, 2016

Name: AK-Nummer: Ergänzungsprüfung January 29, 2016 INSTRUCTIONS: The test has a total of 32 pages including this title page and 9 questions which are marked out of 10 points; ensure that you do not omit a page by mistake. Please write your name and AK-Nummer

More information

Math 128A Spring 2003 Week 12 Solutions

Math 128A Spring 2003 Week 12 Solutions Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems

More information

APPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.

APPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015. APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.

More information

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang

CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION. Professor Dae Ryook Yang CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 5-1 Road Map of the Lecture V Laplace Transform and Transfer

More information

= 2e t e 2t + ( e 2t )e 3t = 2e t e t = e t. Math 20D Final Review

= 2e t e 2t + ( e 2t )e 3t = 2e t e t = e t. Math 20D Final Review Math D Final Review. Solve the differential equation in two ways, first using variation of parameters and then using undetermined coefficients: Corresponding homogenous equation: with characteristic equation

More information

Understanding the Inductances

Understanding the Inductances Understanding the Inductances We have identified six different inductances (or reactances) for characterizing machine dynamics. These are: d, q (synchronous), ' d, ' q (transient), '' d,'' q (subtransient)

More information

εx 2 + x 1 = 0. (2) Suppose we try a regular perturbation expansion on it. Setting ε = 0 gives x 1 = 0,

εx 2 + x 1 = 0. (2) Suppose we try a regular perturbation expansion on it. Setting ε = 0 gives x 1 = 0, 4 Rescaling In this section we ll look at one of the reasons that our ε = 0 system might not have enough solutions, and introduce a tool that is fundamental to all perturbation systems. We ll start with

More information

Equivalent Circuits with Multiple Damper Windings (e.g. Round rotor Machines)

Equivalent Circuits with Multiple Damper Windings (e.g. Round rotor Machines) Equivalent Circuits with Multiple Damper Windings (e.g. Round rotor Machines) d axis: L fd L F - M R fd F L 1d L D - M R 1d D R fd R F e fd e F R 1d R D Subscript Notations: ( ) fd ~ field winding quantities

More information

Exam Basics. midterm. 1 There will be 9 questions. 2 The first 3 are on pre-midterm material. 3 The next 1 is a mix of old and new material.

Exam Basics. midterm. 1 There will be 9 questions. 2 The first 3 are on pre-midterm material. 3 The next 1 is a mix of old and new material. Exam Basics 1 There will be 9 questions. 2 The first 3 are on pre-midterm material. 3 The next 1 is a mix of old and new material. 4 The last 5 questions will be on new material since the midterm. 5 60

More information