From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:27:19 +0530 Subject: initial commit / add all books --- 3793/CH10/EX10.10/exp_10_10.sce | 62 +++++++++++++++++++++++++++++++++++++++++ 1 file changed, 62 insertions(+) create mode 100644 3793/CH10/EX10.10/exp_10_10.sce (limited to '3793/CH10/EX10.10') diff --git a/3793/CH10/EX10.10/exp_10_10.sce b/3793/CH10/EX10.10/exp_10_10.sce new file mode 100644 index 000000000..07622ccf9 --- /dev/null +++ b/3793/CH10/EX10.10/exp_10_10.sce @@ -0,0 +1,62 @@ +clear; +clc; +Sb=30; +vb=11; +sg=20; +p=10; +R=6.6; +//generator +X1=complex(0,.1); +X2=complex(0,.1); +X0=complex(0,.15); +x1=X1*(Sb/sg); +x2=X2*(Sb/sg); +x0=X0*(Sb/sg); +//transformer12 +xt1=complex(0,.12); +xt2=complex(0,.12); +xt0=complex(0,.12); +//transmission line +vtr=22; +Ztr=vtr^2/Sb; +Z=complex(1,5); +Zpu=Z/Ztr; +//transformer34 +Xt1=complex(0,.05); +Xt2=complex(0,.05); +Xt0=complex(0,.05); +xtt1=Xt1*(Sb/sg); +xtt2=Xt2*(Sb/sg); +xtt0=Xt0*(Sb/sg); +Vf3=1; +Rpu=((Vf3^2)*Sb)/p; +Rf=(R*Sb)/vtr^2; +Il=p/Sb; +Vf4=Vf3+(Il*xtt1); +Zfp=((x1+xt1+Zpu)*(Rpu+xtt1))/(x1+xt1++Zpu+Rpu+xtt1); +Zfn=Zfp; +zf0=Zpu+xt1; +Ia1=Vf3/complex(1.611,1.5); //from figure +Ia0=Ia1; +Ia2=Ia1; +Ia=3*Ia1; +Ip23=(Rpu+xtt1)*Ia1/((x1+xt1+Zpu)+Rpu+xtt1); +Ip34=Ia1-Ip23; +In23=Ip23; +In34=Ip34; +Iz23=Ia1; +Iz34=0; +a=complex(-.5,.866); +A=[1 1 1;1 a^2 a;1 a a^2]; +I23=[Iz23;Ip23;In23]; +II=A*I23; +mprintf("Current for node 2-3 is Ia=%f%f A Ib=%f+%f A Ic=%f+%f A\n",real(II(1,1)),imag(II(1,1)),real(II(2,1)),imag(II(2,1)),real(II(3,1)),imag(II(3,1))); +I34=[Iz34;Ip34;In34]; +III=A*I34; +mprintf("Current for node 3-4 is Ia=%f+%f A Ib=%f%f A Ic=%f%f A\n",real(III(1,1)),imag(III(1,1)),real(III(2,1)),imag(III(2,1)),real(III(3,1)),imag(III(3,1))); +Vp2=Vf3+(Zpu*Ip23); +Vn2=Zpu*Ip23; +Vz2=Zpu*Ia1; +Vz=[Vz2;Vp2;Vn2]; +V=A*Vz; +mprintf("Phase voltages at node 2 are Va=%f+%f Kv Vb=%f%f Kv Vc=%f+%f Kv",real(V(1,1)),imag(V(1,1)),real(V(2,1)),imag(V(2,1)),real(V(3,1)),imag(V(3,1))); -- cgit