Scribd
Upload
42 views4 pages

Machine Problem No. 1

The document describes calculating the bus admittance matrix (YBUS) for a 5-bus system with various coupling configurations: 1) Originally with no coupling between branches 2) With branch 3 coupled to branch 1 3) With branch 2 coupled to branch 5 It provides the circuit diagram and shows the code to calculate the original YBUS and modified YBUS matrices for each coupling configuration.

Uploaded by

Neil Doc
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
42 views4 pages

Machine Problem No. 1

The document describes calculating the bus admittance matrix (YBUS) for a 5-bus system with various coupling configurations: 1) Originally with no coupling between branches 2) With branch 3 coupled to branch 1 3) With branch 2 coupled to branch 5 It provides the circuit diagram and shows the code to calculate the original YBUS and modified YBUS matrices for each coupling configuration.

Uploaded by

Neil Doc
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

MACHINE PROBLEM NO.

1
PROBLEM: Using building-block procedure, determine YBUS for the circuit shown below.

a. Assume there is no mutual coupling between any branches.


b. If the branch 3 and 1 is mutually coupled to branch 3 and 2 through Zm =j0.1 p.u.
c. If the branch 2 and 3 is mutually coupled to branch 2 and 5 through Zm =j0.1 p.u.

A1=(-j*1)*[1 0 0 0 0;0 0 0 0 0;0 0 0 0 0;0 0 0 0 0;0 0 0 0 0];

A2=(-j*2.5)*[1 -1 0 0 0;-1 1 0 0 0;0 0 0 0 0;0 0 0 0 0;0 0 0 0 0];

A3=(-j*2)*[1 0 -1 0 0;0 0 0 0 0;-1 0 1 0 0;0 0 0 0 0;0 0 0 0 0];

A4=(-j*4)*[0 0 0 0 0;0 1 -1 0 0 ;0 -1 1 0 0;0 0 0 0 0;0 0 0 0 0];

A5=(-j*5)*[0 0 0 0 0;0 1 0 0 -1;0 0 0 0 0;0 0 0 0 0;0 -1 0 0 1];

A6=(-j*8)*[0 0 0 0 0;0 0 0 0 0;0 0 1 -1 0;0 0 -1 1 0;0 0 0 0 0];

A7=(-j*2)*[0 0 0 0 0;0 0 0 0 0;0 0 0 0 0;0 0 0 1 -1;0 0 0 -1 1];

A8=(-j*0.8)*[0 0 0 0 0;0 0 0 0 0;0 0 0 0 0;0 0 0 0 0;0 0 0 0 1];

Y=A1+A2+A3+A4+A5+A6+A7+A8;

Ybusorig=Y;

Z=j*[0.5 0.1;0.1 0.25];

Ym=inv(Z);
B=[Ym(1,1) -Ym(1,1) Ym(1,2) -Ym(1,2);-Ym(1,1) Ym(1,1) -Ym(1,2) Ym(1,2);Ym(2,1) -Ym(2,1) Ym(2,2) -
Ym(2,2);-Ym(2,1) Ym(2,1) -Ym(2,2) Ym(2,2)];

Y1(1,1)=B(2,2);

Y1(2,1)=B(4,2);

Y1(3,1)=B(1,2)+B(3,2);

Y1(1,2)=Y1(2,1);

Y1(2,2)=B(4,4);

Y1(3,2)=B(1,4)+B(3,4);

Y1(1,3)=Y1(3,1);

Y1(2,3)=Y1(3,2);

Y1(3,3)=B(1,1)+B(3,1)+B(1,3)+B(3,3);

Y1(1,4)=0;

Y1(4,1)=Y1(1,4);

Y1(5,1)=0;

Y1(1,5)=Y1(5,1);

Ybus=Y+Y1;

Y2=(-j*4)*[1 -1 0 0 0;-1 1 0 0 0;0 0 0 0 0;0 0 0 0 0;0 0 0 0 0];

Y3=(-j*2)*[1 0 -1 0 0;0 0 0 0 0;-1 0 1 0 0;0 0 0 0 0;0 0 0 0 0];

Ybus2=Ybus-Y2-Y3;

Ybm32and31=Ybus2;

Z1=j*[0.25 0.1;0.1 0.2];

Ym1=inv(Z1);

B1=[Ym1(1,1) -Ym1(1,1) Ym1(1,2) -Ym1(1,2);-Ym1(1,1) Ym1(1,1) -Ym1(1,2) Ym1(1,2);Ym1(2,1) -Ym1(2,1)


Ym1(2,2) -Ym1(2,2);-Ym1(2,1) Ym1(2,1) -Ym1(2,2) Ym1(2,2)];

Y4(2,2)=B1(1,1)+B1(1,3)+B1(3,1)+B1(3,3);

Y4(2,3)=B1(1,2)+B1(3,2);

Y4(3,2)=Y4(2,3);

Y4(2,5)=B1(1,4)+B1(3,4);

Y4(5,2)=Y4(2,5);

Y4(3,5)=B1(2,4);

Y4(5,3)=Y4(3,5);
Y4(3,3)=B1(2,2);

Y4(5,5)=B1(4,4);

Ybus3=Ybus2+Y4;

Y5=(-j*4)*[0 0 0 0 0;0 1 -1 0 0;0 -1 1 0 0;0 0 0 0 0;0 0 0 0 0];

Y6=(-j*5)*[0 0 0 0 0;0 1 0 0 -1;0 0 0 0 0;0 0 0 0 0;0 -1 0 0 1];

Ybus4=Ybus3-Y5-Y6;

Ybm23and25=Ybus4;

Ybusorig % Original Ybus of the system %

Ybm32and31 % Ybus when branch 3 and 1 is mutually coupled to branch 3 and 2 %

Ybm23and25 % Ybus when branch 2 and 3 is mutually coupled to branch 2 and 5 %

%RUN….%
UNIVERSITY OF THE EAST
COLLEGE OF ENGINEERING
MANILA CAMPUS
ELECTRICAL ENGINEERING DEPARTMENT

NEE – 524

MACHINE PROBLEM NO. 1

SUBMITTED BY:
BONDOC, NEIL EDGARD C.
20140151275

ENGR. DARWIN ALPIS


INSTRUCTOR

You might also like