10 files changed, 230 insertions, 0 deletions
diff --git a/1055/CH4/EX4.1/ch4_1.sce b/1055/CH4/EX4.1/ch4_1.sce
new file mode 100755
index 000000000..927c670a3
--- /dev/null
+++ b/1055/CH4/EX4.1/ch4_1.sce
@@ -0,0 +1,18 @@
+//To detremine the the voltage at the generating station and efficiency of transmission
+clear
+clc;
+R=0.496;// resistance
+X=1.536;
+Vr=2000;
+Z=(10*2*2/(11*11)) + %i*30*2*2/(11*11);
+Zt=(.04+(1.3*2*2/(11*11))) + %i*(.125 + (4.5*2*2/(11*11)));//Transformer impedence
+Il=250*1000/2000;// line current(amps.)
+Pl=Il*Il*R;//line loss(kW)
+Po=250*0.8;// output(kW)
+cosr=0.8;// power factor
+sinr=.6;
+%n=200*100/(200+7.7);
+Vs=(Vr*cosr+Il*R)+%i*(Vr*sinr+Il*X);
+V=sqrt((1662^2)+ (1392^2));
+mprintf("efficiency= %.1f percent \n",%n);
+mprintf("Sending end voltage,|Vs|=%.0f volts",V);
diff --git a/1055/CH4/EX4.10/ch4_10.sce b/1055/CH4/EX4.10/ch4_10.sce
new file mode 100755
index 000000000..cf688a33b
--- /dev/null
+++ b/1055/CH4/EX4.10/ch4_10.sce
@@ -0,0 +1,31 @@
+// To determine the sending end voltage and current , power and power factor. Evaluate A, B , C, D parameters.
+clear
+clc;
+R=.1557*160
+GMD=(3.7*6.475*7.4)^(1/3);
+Z1=2*(10^-7)*log(560/.978)*160*1000;
+XL=63.8;
+C=(10^-9)*2*(10^6)*%pi*160*1000/(36*%pi*log(560/.978));
+Z=sqrt((.1557^2) + .39875^2)*(cosd(68.67)+ %i*sind(68.67));
+jwC=%i*314*1.399*(10^-6)/160;
+Zc=sqrt(Z/jwC);
+y=sqrt(Z*jwC);
+yl=y*160;
+A=cosh(yl);
+B=Zc*sinh(yl)
+C=sinh(yl)/Zc;
+Ir=50000/(sqrt(3)*132);
+Vs=(A*76.208) +(B*(10^-3)*Ir*(cosd(-36.87)+%i*sind(-36.87)));
+VS=152.34;
+Is=C*76.208*(10^3) +(A*Ir*(cosd(-36.87)+%i*sind(-36.87)));
+Ps=3*abs(Vs)*abs(Is)*cosd(33.96);
+pf=cosd(33.96);
+Vnl=abs(Vs)/abs(A);
+reg=(Vnl-76.208)*100/76.208;
+n=50000*.8*100/abs(Ps);
+mprintf(" Vs line to line =%.2f kV\n",VS);
+disp(Is,"sending end current Is(A)=");//Answer don't match due to difference in rounding off of digits
+mprintf("sending end power=%.0f kW\n",Ps);
+mprintf("sending end p.f =%.3f\n",pf);
+mprintf("percent regulation=%.1f \n",reg);
+mprintf("percent efficency=%.1f ",n);
diff --git a/1055/CH4/EX4.2/ch4_2.sce b/1055/CH4/EX4.2/ch4_2.sce
new file mode 100755
index 000000000..d66746400
--- /dev/null
+++ b/1055/CH4/EX4.2/ch4_2.sce
@@ -0,0 +1,21 @@
+//To determine power input and output (i) star connected (ii)delta connected
+clear
+clc;
+mprintf("when load is star connected\n");
+Vln=400/sqrt(3);// Line to neutral voltage(V)
+Z=7+ %i*11;//Impedence per phase
+Il=231/Z;// line current(amp.)
+I=abs(231/Z);
+Pi=3*I*I*7;
+Po=3*I*I*6;
+mprintf("power input =%.0f watts\n",Pi);//Answers don't match due to difference in rounding off of digits
+mprintf("power output=%.0f watts\n",Po);//Answers don't match due to difference in rounding off of digits
+mprintf("when load is delta connected\n");
+Ze=2+ %i*3;// equivalent impedence(ohm)
+Zp=3+%i*5;// impedence per phase
+il=231/Zp;//Line current(amps.)
+IL=abs(il);
+pi=3*IL*IL*3;
+po=3*IL*IL*2;
+mprintf("power input=%.1f watts\n",pi);//Answers don't match due to difference in rounding off of digits
+mprintf("power output = %.0f watts \n",po);//Answers don't match due to difference in rounding off of digits
diff --git a/1055/CH4/EX4.3/ch4_3.sce b/1055/CH4/EX4.3/ch4_3.sce
new file mode 100755
index 000000000..66610a557
--- /dev/null
+++ b/1055/CH4/EX4.3/ch4_3.sce
@@ -0,0 +1,19 @@
+// To determine efficiency and regulation of line
+clear
+clc;
+a=100/.5
+Xl=2*(10^-7)*log(100/.5);//inductance(H/meter)
+XL=20*(1000)*Xl;// inductance of 20 km length 
+R=6.65;// resistance(ohm)
+Rc=20*1000/(58*90);// resistance of copper(ohm)
+I=10*1000/(33*.8*sqrt(3));// the current(amps.)
+Pl=3*I*I*Rc/(10^6);//loss (MW)
+n=10/(10+Pl);
+mprintf("efficiency=%.4f percent \n",n);
+Vr=19052;
+cosr=.8;//power factor
+sinr=.6;
+Vs=abs(((Vr*cosr+I*Rc) +%i*(Vr*sinr+ I*R)));
+mprintf("Vs =%.0f volts\n",Vs);//Answer don't match due to difference in rounding off of digits
+Reg=(Vs-Vr)*100/Vr;
+mprintf(" regulation =%.2f percent",Reg)
diff --git a/1055/CH4/EX4.4/ch4_4.sce b/1055/CH4/EX4.4/ch4_4.sce
new file mode 100755
index 000000000..f9e195f67
--- /dev/null
+++ b/1055/CH4/EX4.4/ch4_4.sce
@@ -0,0 +1,14 @@
+//To calculate the voltage across each load impedence and current in the nuetral
+clear
+clc;
+IR=(400)/((sqrt(3)*(6.3+%i*9)));
+IY=231*(cosd(-120) + %i*sind(-120))/8.3;
+IB=231*(cosd(120) + %i*sind(120))/(6.3-%i*8);
+In=abs((IR +IY +IB));//Neutral current
+mprintf("Neutral current =%.2f amps\n",In);
+VR=abs(IR*(6+ %i*9));
+VY=abs(IY*(8));
+VB=abs(IB*(6-%i*8));
+mprintf("Voltage across Phase R =%.1f volts \n",VR);
+mprintf("Voltage across Phase Y =%.2f volts \n",VY);
+mprintf("Voltage across Phase B =%.0f volts \n",VB);
diff --git a/1055/CH4/EX4.5/ch4_5.sce b/1055/CH4/EX4.5/ch4_5.sce
new file mode 100755
index 000000000..d6f54829e
--- /dev/null
+++ b/1055/CH4/EX4.5/ch4_5.sce
@@ -0,0 +1,38 @@
+// To determine efficiency and regulation of 3 phase line
+clear
+clc;
+R=100*.1;//Resistance of line (ohm)
+Xl=2*(10^-7)*100*1000*log(200/.75);//inductance of line
+X2=Xl*314;//inductive reactance
+C=2*(%pi*100)*8.854*(10^-12)*100*1000*(10^6)/(log(200/.75));// capacitance per phase (micro farad)
+mprintf("Using Nominal-T method\n");
+Ir=20*1000/(sqrt(3)*66*.8);
+Vr=66*1000/sqrt(3);
+Vc=(38104*.8+ Ir*5) +%i*(38104*.6+ Ir*17.55);// voltage across condenser
+Ic=%i*314*(Vc)*.9954*(10^-6);
+is=Ir+Ic;
+Is=abs(Ir+Ic);
+Vs=abs(Vc + (is*(5 + %i*17.53)));
+VR=abs(Vs*(-%i*3199)/(5-%i*3181));// no load recieving end voltage
+Reg=(VR-Vr)*100/Vr;
+Pl=3*(Ir*Ir*5 + Is*Is*5)/1000000;
+%n=20*100/(20+Pl);
+mprintf("percent regulation=%.1f  \n",Reg);
+mprintf("percent efficiency=%.1f  \n\n",%n);
+mprintf("Using Nominal-pi method\n");
+Ir1=218.68*(.8-%i*.6);
+Ic1=%i*314*.4977*(10^-6)*Vr;
+Il=Ir1+Ic1;
+vs1=Vr+Il*(10+%i*35.1);
+Vs1=abs(vs1);
+Vr1=Vs1*(-%i*6398)/(10-%i*6363);
+VR1=abs(Vr1);// no load recieving end voltage
+Reg2=(VR1-Vr)*100/Vr;
+IL=abs(Ir1+Ic1);
+Loss=3*IL*IL*10;
+%n=20*100/21.388;
+mprintf("percent regulation=%.2f  \n",Reg2);
+mprintf("percent efficiency=%.1f  \n",%n);
+
+
+
diff --git a/1055/CH4/EX4.6/ch4_6.sce b/1055/CH4/EX4.6/ch4_6.sce
new file mode 100755
index 000000000..2d0d399bd
--- /dev/null
+++ b/1055/CH4/EX4.6/ch4_6.sce
@@ -0,0 +1,25 @@
+//To find the rms value and phase values (i)The incident voltage to neutral at the recieving end (ii)The reflected voltage to neutral at the recieving end (iii)The incident and reflected voltage to neutral at 120 km from the recieving end.
+clear
+clc;
+R=0.2;
+L=1.3;
+C=0.01*(10^-6);
+z=R+%i*L*314*(10^-3);// serie impedence
+y=%i*314*C;// shunt admittance
+Zc=sqrt(z/y);// characterstic impedence
+Y=sqrt(y*z);
+Vr=132*1000/sqrt(3);
+Ir=0;
+Vin=(Vr + Ir*Zc)/2;// incident voltage to neutral at the recieving end
+mprintf("Vr =%.3f volts \n",Vr);//Answer don't match due to difference in rounding off of digits
+mprintf("(i)The incident voltage to neutral at the recieving end =%.3f volts \n",Vin);//Answer don't match due to difference in rounding off of digits
+Vin2=(Vr - Ir*Zc)/2;// The reflected voltage to neutral at the recieving end
+mprintf("(ii)The reflected voltage to neutral at the recieving end=%.3f volts \n",Vin2);//Answer don't match due to difference inrounding off of digits
+Vrp=Vr*exp(.2714*120*(10^-3))*exp(%i*1.169*120*(10^-3))/1000;//Taking Vrp=Vr+
+Vrm=Vr*exp(-0.0325)*exp(-%i*.140)/1000;//Taking Vrm=Vr-
+v1=Vrm/2;// reflected voltage to neutral at 120 km from the recieving end
+phase_v1=atand(imag(v1)/real(v1));//Phase angle of v1
+v2=Vrp/2;//incident voltage to neutral at 120 km from the recieving end
+phase_v2=atand(imag(v2)/real(v2));//Phase angle of v2
+mprintf("(iii) reflected voltage to neutral at 120 km from the recieving end =%.2f at angle of %.2f \n",abs(v1),phase_v1);
+mprintf("incident voltage to neutral at 120 km from the recieving end = %.2f at angle of %.2f\n",abs(v2),phase_v2);
diff --git a/1055/CH4/EX4.7/ch4_7.sce b/1055/CH4/EX4.7/ch4_7.sce
new file mode 100755
index 000000000..8277c6c01
--- /dev/null
+++ b/1055/CH4/EX4.7/ch4_7.sce
@@ -0,0 +1,17 @@
+//To determine of efficiency of line
+clear
+clc;
+Ir=40*1000/(sqrt(3)*132*.8);
+Vr=132*1000/sqrt(3);
+Zc=380*(cosd(-13.06)+ %i*sind(-13.06));
+IR=Ir*(cosd(-36.8)+ %i*sind(-36.8));
+Vsp=(Vr+IR*Zc)*(1.033*(cosd(8.02)+ %i*sind(8.02)))/2;
+Vsm=(Vr-IR*Zc)*(.968*(cosd(-8.02)+ %i*sind(-8.02)))/2;
+vs=Vsp+ Vsm;
+Vs=abs(vs);
+is=(Vsp-Vsm)/Zc;
+Is=abs(is)
+P=3*Vs*Is*cosd(33.72)/10^6;
+n=40*100/P;
+mprintf("efficiency=%.1f",n);//Answer don't match due to difference in rounding off of digits
+
diff --git a/1055/CH4/EX4.8/ch4_8.sce b/1055/CH4/EX4.8/ch4_8.sce
new file mode 100755
index 000000000..5b284b701
--- /dev/null
+++ b/1055/CH4/EX4.8/ch4_8.sce
@@ -0,0 +1,23 @@
+//To determine the ABCD parameters of Line
+clear
+clc;
+yl=(0.2714+ %i*1.169)*120*(10^-3);
+Ir=40*1000/(sqrt(3)*132*.8)
+A=cosh(yl);
+phase_A=atand(imag(A)/real(A));//Phase angle of A
+IR=Ir*(cosd(-36.8)+ %i*sind(-36.8))
+Vr=132*1000/sqrt(3);
+Zc=380*(cosd(-13.06)+ %i*sind(-13.06));
+B=Zc*sinh(yl);
+phase_B=atand(imag(B)/real(B));//Phase angle of B
+Vs=(A*Vr+B*IR);
+f=abs(B);
+d=abs(Vs);
+C=sinh(yl)/Zc;
+phase_C=atand(imag(C)/real(C));//Phase angle of C
+D=cosh(yl);
+phase_D=atand(imag(D)/real(D));//Phase angle of D
+mprintf("A=%.2f at an angle of %.2f \n",abs(A),phase_A)
+mprintf("B=%.1f at an angle of %.0f \n",abs(B),phase_B)
+mprintf("C=%.2f at an angle of %.2f \n",abs(C),phase_C)
+mprintf("D=%.2f at an angle of %.2f \n",abs(D),phase_D)
diff --git a/1055/CH4/EX4.9/ch4_9.sce b/1055/CH4/EX4.9/ch4_9.sce
new file mode 100755
index 000000000..cd55d2a41
--- /dev/null
+++ b/1055/CH4/EX4.9/ch4_9.sce
@@ -0,0 +1,24 @@
+//To determine the sending end voltage and efficiency using Nominal_pi and Nominal-T method
+clear
+clc;
+Ir=218.7*(.8-%i*.6);
+Ic1=%i*314*.6*(10^-6)*76200;
+Il=Ic1+Ir;
+Vs=76200 + Il*(24+ %i*48.38);
+phase_Vs=atand(imag(Vs)/real(Vs));//phase angle of VS
+Pl=3*24*abs(Il)*abs(Il)/1000000;//The Loss(MW)
+n=40*100/(40+Pl);
+mprintf("Using Nominal- pi method\n ");
+mprintf("Vs=%.0f volts at an angle of %.2f \n",abs(Vs),phase_Vs)
+mprintf("efficiency=%.2f percent\n",n);
+mprintf("\nUsing Nominal-T method\n");
+Vc=76200*(.8+%i*.6) + 218.7*(12+ %i*24.49);
+Ic=%i*314*1.2*(10^-6)*(63584+ %i*51076);
+Is=199.46+ %i*23.95;
+Vs=(Vc + Is*(12+ %i*24.49))/1000;
+phase_Vs=atand(imag(Vs)/real(Vs));//Phase angle of Vs
+Pl1=3*12*((200.89^2)+ 218.7^2)/1000000;//The loss(MW)
+n1=40*100/(40+Pl1);
+mprintf("Vs=%.2f  at an angle of %.2f \n",abs(Vs),phase_Vs)
+mprintf("efficiency=%.2f percent\n",n1);
+