Scribd
Upload
49 views12 pages

Computer Aided Power System Analysis Inst. DR - Weerakorn Ongsakul

This document provides information about power system analysis using the Newton-Raphson method. It includes: 1) An example of constructing the Y-matrix from line data for a 6-bus, 11-line power system. 2) Two assignments - the first asks to modify the Y-matrix for a capacitor added to bus 6 and a line switched off, the second asks to use Newton-Raphson to solve three nonlinear equations. 3) A third assignment similar to the second but using different initial values for the variables to solve the equations. 4) A reminder that the lab project is to create an m-file to solve power flow using Newton-Raphson and is due

Uploaded by

Krish Dah
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
49 views12 pages

Computer Aided Power System Analysis Inst. DR - Weerakorn Ongsakul

This document provides information about power system analysis using the Newton-Raphson method. It includes: 1) An example of constructing the Y-matrix from line data for a 6-bus, 11-line power system. 2) Two assignments - the first asks to modify the Y-matrix for a capacitor added to bus 6 and a line switched off, the second asks to use Newton-Raphson to solve three nonlinear equations. 3) A third assignment similar to the second but using different initial values for the variables to solve the equations. 4) A reminder that the lab project is to create an m-file to solve power flow using Newton-Raphson and is due

Uploaded by

Krish Dah
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 12

ED72.

10
Computer Aided Power System Analysis
Inst. Dr.Weerakorn Ongsakul
Programming in Power System
Songkran Pisanupoj ([email protected])
7 September 2012
Power network representation: Y-matrix
2
E Y I =
3
L L
L
E Y
I

=
|
( ) ( )
2 2
13
1
12
1 13 3 1 12 2 1 1
b
E
b
E y E E y E E I + + + =
| | | |
3 13 2 12 1
13 12
13 12 1
2 2
E y E y E
b b
y y I + +
(

+ + + =
| | | | | |
3 13 2 12 1 11 1
E Y E Y E Y I + + =
(
(
(

(
(
(

=
(
(
(

3
2
1
33 32 31
23 22 21
13 12 11
3
2
1
E
E
E
Y Y Y
Y Y Y
Y Y Y
I
I
I
Constructing Y-matrix
Line
no.
From
bus
To
bus
R
(pu)
X
(pu)
B
(pu)
1 1 2 0.1 0.2 0.02
2 1 3 0.05 0.2 0.02
3 2 3 0.08 0.3 0.03
2
02 0
2
02 0
2 0 05 0
1
2 0 1 0
1
. j . j
. j .
. j .
+ +
+
+
+
2 0 1 0
1
. j . +

2 0 05 0
1
. j . +

2
03 0
2
02 0
3 0 08 0
1
2 0 1 0
1
. j . j
. j .
. j .
+ +
+
+
+
2 0 1 0
1
. j . +

2 0 05 0
1
. j . +

2
03 0
2
02 0
3 0 08 0
1
2 0 05 0
1
. j . j
. j .
. j .
+ +
+
+
+
3 0 08 0
1
. j . +

3 0 08 0
1
. j . +

Y =
3
To solve for x from f(x)=k

f
1
(x
1
,x
2
) = k
1
f
2
(x
1
,x
2
) = k
2


k
1
= f
1
(x
1
,x
2
) = f
1
(x
1
(0)
+x
1
(0)
, x
2
(0)
+x
2
(0)
)

k
2
= f
2
(x
1
,x
2
) = f
2
(x
1
(0)
+x
1
(0)
, x
2
(0)
+x
2
(0)
)


By Taylors expansion, you get



k
1
= f
1
(x
1
,x
2
) = f
1
(x
1
(0)
, x
2
(0)
) + x
1
(0)
[f
1
/ x
1
]
(0)
+ x
2
(0)
[f
1
/ x
2
]
(0)


k
2
= f
2
(x
1
,x
2
) = f
2
(x
1
(0)
, x
2
(0)
) + x
1
(0)
[f
2
/ x
1
]
(0)
+ x
2
(0)
[f
2
/ x
2
]
(0)


Newton-Raphson Method
You are given two functions (f
1
and f
2
)
And ask to find x
1
and x
2
You guess x
1
, then
you need to find x
1
4
To solve for x from f(x)=k
reform in to matrix, you get:










in general,
Newton-Raphson Method
( ) ( )
( )
( ) ( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
(

c c c c
c c c c
=
(

0
2
0
1
0
2 2
0
2 1
0
2 1
0
1 1
0
2
0
1 2 2
0
2
0
1 1 1

x
x
x f x f
x f x f
x , x f k
x , x f k
( )
( )
( )
( )
( )
(

=
(

0
2
0
1 0
0
2
0
1

x
x
k
k
( )
( )
( )
( )
( )
(

=
(


0
2
0
1 0 1
0
2
0
1

k
k
x
x
( ) ( ) ( )
( ) ( ) ( ) 0
2
0
2
1
2
0
1
0
1
1
1

x x x
x x x
+ =
+ =
( ) ( ) ( )
( ) ( ) ( ) n n n
n n n
x x x
x x x
2 2
1
2
1 1
1
1

+ =
+ =
+
+
5
Newton-Raphson Method
Flowchart

6

arbitrarily set x
(0)
calculate f(x
(o)
)
calculate K
(0)
=K-f(x
(0)
)
calculate J
(0)
= [f/x]
(0)
calculate x
(0)
= J
-1(0)
K
(0)
calculate x
(1)
= x
(0)
+ x
(0)
Newton-Raphson Method
To solve for x from f(x) = K

7
calculate f(x
(1)
)

Lab.3 Assignment
1. Create the shortest m-file to construct Y-matrix from line data
given: 6-bus, 11-line system.
8
Line no. From bus To bus R (pu) X (pu) B (pu)
1 1 2 0.10 0.20 0.02
2 1 4 0.05 0.20 0.02
3 1 5 0.08 0.30 0.03
4 2 3 0.05 0.25 0.03
5 2 4 0.05 0.10 0.01
6 2 5 0.10 0.30 0.02
7 2 6 0.07 0.20 0.025
8 3 5 0.12 0.26 0.025
9 3 6 0.02 0.10 0.01
10 4 5 0.20 0.40 0.04
11 5 6 0.10 0.30 0.03
Lab.3 Assignment
2. With using a system given in Assignment(1), how to treat a
capacitor to be placed at bus 6. The capacitor has susceptance
of 0.1 pu. How it affects to the diagram and the Y-matrix?
And what to do on Y-matrix if a line between buses 2 and 6
(line no.7) is switched off?


Please provide your descriptive comments on this inquiry.
Provide of resulting values of Y-matrix or programming source
code is considered complementary.
9
Lab.3 Assignment
3. Write m-file to solve the below problem by using Newton-
Raphson method.
F
1
(x) = 3x
1
+ 4x
2
+ 5x
3
= 49
F
2
(x) = 4x
1
+ 7x
2
+ 3x
3
= 51
F
3
(x) = 5x
1
+ 6x
2
+ 5x
3
= 61
Hint: using any initial X, and applying



Report the value of Jacobian matrix (J),
the final value of K, x and X (x
1
, x
2
, x
3
)
10
( ) ( )
( )
( ) ( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
(

c c c c
c c c c
=
(

0
2
0
1
0
2 2
0
2 1
0
2 1
0
1 1
0
2
0
1 2 2
0
2
0
1 1 1

x
x
x f x f
x f x f
x , x f k
x , x f k
Lab.3 Assignment
4. Write m-file to solve the below problem by using Newton-
Raphson method.
F
1
(x) = 3x
1
+ 4x
1
x
2
+ 5x
3
= 70
F
2
(x) = 4x
1
+ 7x
2
+ 3x
3
= 51
F
3
(x) = 5x
1
+ 6x
2
+ 5x
3
= 61
Hint: using


With using different initial X: (0,0,0), (1,1,1), (5,5,5), (9,9,9);
Report the value of Jacobian matrix (J),
the final value of J and final value of X (x
1
, x
2
, x
3
)
11
( ) ( )
( )
( ) ( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
(

c c c c
c c c c
=
(

0
2
0
1
0
2 2
0
2 1
0
2 1
0
1 1
0
2
0
1 2 2
0
2
0
1 1 1

x
x
x f x f
x f x f
x , x f k
x , x f k
Reminder Lab. Project
Create the shortest m-file(s) to solve for power flow
analysis using Newton-Raphson Method

Submission 9 November 2012

End of Lab.3 Note
Thank you
12

You might also like