summaryrefslogtreecommitdiff
path: root/2078/CH5
diff options
context:
space:
mode:
Diffstat (limited to '2078/CH5')
-rwxr-xr-x2078/CH5/EX5.1/Example5_1.sce20
-rwxr-xr-x2078/CH5/EX5.10/Example5_10.sce33
-rwxr-xr-x2078/CH5/EX5.11/Example5_11.sce30
-rwxr-xr-x2078/CH5/EX5.12/Example5_12.sce24
-rwxr-xr-x2078/CH5/EX5.13/Example5_13.sce27
-rwxr-xr-x2078/CH5/EX5.14/Example5_14.sce22
-rwxr-xr-x2078/CH5/EX5.15/Example5_15.sce31
-rwxr-xr-x2078/CH5/EX5.16/Example5_16.sce33
-rwxr-xr-x2078/CH5/EX5.17/Example5_17.sce24
-rwxr-xr-x2078/CH5/EX5.18/Example5_18.sce10
-rwxr-xr-x2078/CH5/EX5.19/Example5_19.sce18
-rwxr-xr-x2078/CH5/EX5.2/Example5_2.sce21
-rwxr-xr-x2078/CH5/EX5.20/Example5_20.sce36
-rwxr-xr-x2078/CH5/EX5.21/Example5_21.sce29
-rwxr-xr-x2078/CH5/EX5.22/Example5_22.sce44
-rwxr-xr-x2078/CH5/EX5.23/Example5_23.sce29
-rwxr-xr-x2078/CH5/EX5.24/Example5_24.sce31
-rwxr-xr-x2078/CH5/EX5.25/Example5_25.sce36
-rwxr-xr-x2078/CH5/EX5.3/Example5_3.sce23
-rwxr-xr-x2078/CH5/EX5.4/Example5_4.sce31
-rwxr-xr-x2078/CH5/EX5.5/Example5_5.sce22
-rwxr-xr-x2078/CH5/EX5.6/Example5_6.sce25
-rwxr-xr-x2078/CH5/EX5.7/Example5_7.sce19
-rwxr-xr-x2078/CH5/EX5.8/Example5_8.sce30
-rwxr-xr-x2078/CH5/EX5.9/Example5_9.sce34
25 files changed, 682 insertions, 0 deletions
diff --git a/2078/CH5/EX5.1/Example5_1.sce b/2078/CH5/EX5.1/Example5_1.sce
new file mode 100755
index 000000000..2747eb4ef
--- /dev/null
+++ b/2078/CH5/EX5.1/Example5_1.sce
@@ -0,0 +1,20 @@
+//Exa 5.1
+clc;
+clear;
+close;
+//Given data :
+P=1100;//kW
+VR=11*1000;//V
+pf=0.8;//power factor
+R=2;//ohm
+X=3;//ohm
+I=P*1000/VR/pf;//A
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+VS=sqrt((VR*cos_fi_r+I*R)^2+(VR*sin_fi_r+I*X)^2);//V
+disp(VS,"Voltage at sending end(V)");
+Reg=(VS-VR)/VR*100;//%
+disp(Reg,"% Regulation");
+LineLoss=I^2*R/1000;//kW
+Eta_T=P*100/(P+LineLoss);//%
+disp(Eta_T,"Transmission Efficiency(%)");
diff --git a/2078/CH5/EX5.10/Example5_10.sce b/2078/CH5/EX5.10/Example5_10.sce
new file mode 100755
index 000000000..f32f97e6f
--- /dev/null
+++ b/2078/CH5/EX5.10/Example5_10.sce
@@ -0,0 +1,33 @@
+//Exa 5.10
+clc;
+clear;
+close;
+//Given data :
+l=200;//km
+P=50;//MVA
+VRL=132*10^3;//Volt
+f=50;//Hz
+R=l*0.15;//ohm
+X=l*0.50;//ohm
+Y=l*2*10^-6;//mho
+pf=0.85;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+VR=VRL/sqrt(3);//Volt
+IR=P*10^6/(sqrt(3)*VRL);//A
+Z=R+%i*X;//ohm
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+Vdash=VR+1/2*IR*Z;//Volt
+IC=Vdash*%i*Y;//A
+IS=IR+IC;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+VS=Vdash+1/2*IS*Z;//Volt
+VSL=abs(VS)*sqrt(3);//Volt
+disp(VSL/1000,"Sending end line voltage(kV) :");
+Reg=(VSL-VRL)/VRL*100;//%
+disp(Reg,"Regulation(%) : ");
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fi_s=cosd(fi_s);//sending end pf
+Eta_T=sqrt(3)*VRL*abs(IR)*cos_fi_r/(sqrt(3)*VSL*abs(IS)*cos_fi_s)*100;//%
+disp(Eta_T,"Transmission Efficiency(%) : ");
+//Ans is wrong in the book.Angle of VS is calculated wrong leads to wrong answers.
diff --git a/2078/CH5/EX5.11/Example5_11.sce b/2078/CH5/EX5.11/Example5_11.sce
new file mode 100755
index 000000000..e64922209
--- /dev/null
+++ b/2078/CH5/EX5.11/Example5_11.sce
@@ -0,0 +1,30 @@
+//Exa 5.11
+clc;
+clear;
+close;
+//Given data :
+S=1*10^3;//kVA
+pf=0.71;//power factor
+VRL=22*10^3;//Volt
+f=50;//Hz
+R=15;//ohm
+L=0.2;//H
+C=0.5*10^-6;//F
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+IR=S*10^3/VRL;//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+X=2*%pi*f*L;//ohm
+//Z=sqrt(R^2+X^2);//ohm
+Z=R+%i*X;//ohm
+Y=2*%pi*f*C;//S
+ICR=1/2*%i*Y*VRL;//A
+IL=IR+ICR;//A
+VS=VRL+IL*Z;//Volt
+disp("Sending end voltage(Volt), magnitude is "+string(abs(VS))+" and angle in degree is "+string(atand(imag(VS),real(VS))));
+ICS=1/2*%i*Y*VS;//A
+IS=IL+ICS;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fi_s=cosd(fi_s);//sending end pf
+disp(cos_fi_s,"Sending end power factor(lag) : ");
diff --git a/2078/CH5/EX5.12/Example5_12.sce b/2078/CH5/EX5.12/Example5_12.sce
new file mode 100755
index 000000000..d3cec0408
--- /dev/null
+++ b/2078/CH5/EX5.12/Example5_12.sce
@@ -0,0 +1,24 @@
+//Exa 5.12
+clc;
+clear;
+close;
+//Given data :
+P=50*10^6;//W
+f=50;//Hz
+l=150;//km
+pf=0.8;//power factor
+VRL=110*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+R=0.1*l;//ohm
+XL=0.5*l;//ohm
+Z=R+%i*XL;//ohm
+IR=P/(sqrt(3)*VRL*pf);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+Y=3*10^-6*l;//S
+ICR=1/2*%i*Y*VR;//A
+IL=IR+ICR;//A
+VS=VR+IL*Z;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end line to line voltage(kV) :");
diff --git a/2078/CH5/EX5.13/Example5_13.sce b/2078/CH5/EX5.13/Example5_13.sce
new file mode 100755
index 000000000..2626f11c1
--- /dev/null
+++ b/2078/CH5/EX5.13/Example5_13.sce
@@ -0,0 +1,27 @@
+//Exa 5.13
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+l=30;//km
+Z=40+%i*125;//ohm
+Y=10^-3;//mho
+P=50*10^6;//W
+VRL=220*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+IR=P/(sqrt(3)*VRL*pf);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+ICR=1/2*%i*Y*VR;//A
+IL=IR+ICR;//A
+VS=VR+IL*Z;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end line to line voltage(kV) :");
+IS=IL+1/2*%i*Y*VS;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fis=cosd(fi_s);//sending end pf
+disp(cos_fis,"Sending end power factor(lag) : ");
diff --git a/2078/CH5/EX5.14/Example5_14.sce b/2078/CH5/EX5.14/Example5_14.sce
new file mode 100755
index 000000000..b58e35487
--- /dev/null
+++ b/2078/CH5/EX5.14/Example5_14.sce
@@ -0,0 +1,22 @@
+//Exa 5.14
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+l=30;//km
+Z=40+%i*125;//ohm
+Y=10^-3;//mho
+P=50*10^6;//W
+VRL=220*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+IR=P/(sqrt(3)*VRL*pf);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+ICR=1/2*%i*Y*VR;//A
+IL=IR+ICR;//A
+VS=VR+IL*Z;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end line to line voltage(kV) :");
diff --git a/2078/CH5/EX5.15/Example5_15.sce b/2078/CH5/EX5.15/Example5_15.sce
new file mode 100755
index 000000000..09a620c88
--- /dev/null
+++ b/2078/CH5/EX5.15/Example5_15.sce
@@ -0,0 +1,31 @@
+//Exa 5.15
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+l=100;//km
+P=50*10^6;//W
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+VRL=132*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+R=0.1*l;//ohm
+XL=0.3*l;//ohm
+Z=R+%i*XL;//ohm
+Y=3*10^-6*l;//S
+IR=P/(sqrt(3)*VRL*pf);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+ICR=1/2*%i*Y*VR;//A
+IL=IR+ICR;//A
+VS=VR+IL*Z;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end line voltage(kV) :");
+ICS=1/2*%i*Y*VS;//A
+IS=IL+ICS;//A
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fi_s=cosd(fi_s);//sending end pf
+disp(cos_fi_s,"Sending end power factor(lag) : ");
+Eta_T=sqrt(3)*VRL*abs(IR)*cos_fi_r/(sqrt(3)*VSL*abs(IS)*cos_fi_s)*100;//%
+disp(Eta_T,"Transmission Efficiency(%) : ");
diff --git a/2078/CH5/EX5.16/Example5_16.sce b/2078/CH5/EX5.16/Example5_16.sce
new file mode 100755
index 000000000..950f48d97
--- /dev/null
+++ b/2078/CH5/EX5.16/Example5_16.sce
@@ -0,0 +1,33 @@
+//Exa 5.16
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+l=10;//km
+S1=5000*10^3;//VA
+S2=10000*10^3;//VA
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+pf2=0.7071;//power factor
+cos_fi_r2=pf2;
+sin_fi_r2=sqrt(1-cos_fi_r2^2);
+R=0.6*l;//ohm
+XL=1.5*l;//ohm
+VRL=33*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+I1=S1/(sqrt(3)*VRL);//A
+I1=I1*(cos_fi_r-%i*sin_fi_r);//A
+Z1=R+%i*XL;//ohm
+VB=VR+I1*Z1;//Volt
+VBL=sqrt(3)*abs(VB);//Volt
+disp(VBL/1000,"Line voltage at mid point(kV) : ");
+I2=S2/(sqrt(3)*VBL);//A
+I2=I2*(cos_fi_r2-%i*sin_fi_r2);//A
+I=I1+I2;//A
+disp("Total current(A), magnitude is "+string(abs(I))+" and angle in degree is "+string(atand(imag(I),real(I))));
+Z2=R+%i*XL;//ohm
+VS=VB+I*Z2;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end line voltage(kV) :");
diff --git a/2078/CH5/EX5.17/Example5_17.sce b/2078/CH5/EX5.17/Example5_17.sce
new file mode 100755
index 000000000..d764de559
--- /dev/null
+++ b/2078/CH5/EX5.17/Example5_17.sce
@@ -0,0 +1,24 @@
+//Exa 5.17
+clc;
+clear;
+close;
+//Given data :
+P=10;//MWatt
+pf=0.8;//power factor
+VRL=30*10^3;//Volt
+R1=5.5;//ohm
+XL1=13.5;//ohm
+R2=6;//ohm
+XL2=11;//ohm
+ZA=R1+%i*XL1;//ohm
+ZB=R2+%i*XL2;//ohm
+S=P*10^3/pf*expm(%i*%pi/180*(-36.52));//kVA
+SA=S*ZB/(ZA+ZB);//kVA
+disp("Load supply by line A(kVA), magnitude is "+string(abs(SA))+" at pf "+string(cosd(atand(imag(SA),real(SA)))));
+SB=S*ZA/(ZA+ZB);//kVA
+disp("Load supply by line B(kVA), magnitude is "+string(abs(SB))+" and angle in degree is "+string(cosd(atand(imag(SB),real(SB)))));
+PA=abs(SA)*(cosd(atand(imag(SA),real(SA))));//kW
+disp(PA,"Power supplied by line A(kW) : ");
+PB=abs(SB)*(cosd(atand(imag(SB),real(SB))));//kW
+disp(PB,"Power supplied by line B(kW) : ");
+//Answer is not accurate in the book.
diff --git a/2078/CH5/EX5.18/Example5_18.sce b/2078/CH5/EX5.18/Example5_18.sce
new file mode 100755
index 000000000..8ba807161
--- /dev/null
+++ b/2078/CH5/EX5.18/Example5_18.sce
@@ -0,0 +1,10 @@
+//Exa 5.18
+clc;
+clear;
+close;
+//Given data :
+L=200;//km
+f=50;//Hz
+omega=2*%pi*f;//rad/s
+Rise=omega^2*L^2*10^-8/18;//%
+disp(Rise,"Percentage rise in voltage : ");
diff --git a/2078/CH5/EX5.19/Example5_19.sce b/2078/CH5/EX5.19/Example5_19.sce
new file mode 100755
index 000000000..188b22591
--- /dev/null
+++ b/2078/CH5/EX5.19/Example5_19.sce
@@ -0,0 +1,18 @@
+//Exa 5.19
+clc;
+clear;
+close;
+//Given data :
+L=80;//km
+f=50;//Hz
+Z=(0.15+%i*0.78)*L;//ohm
+Y=(%i*5*10^-6)*L;//mho
+A=1+1/2*Y*Z;//parameter of 3-phase line
+D=A;//parameter of 3-phase line
+B=Z*(1+1/4*Y*Z);//parameter of 3-phase line
+C=Y;//parameter of 3-phase line
+disp(A,"Parameter A : ");
+disp(B,"Parameter B : ");
+disp(C,"Parameter C : ");
+disp(D,"Parameter D : ");
+//Answer of B is wrong in the book.
diff --git a/2078/CH5/EX5.2/Example5_2.sce b/2078/CH5/EX5.2/Example5_2.sce
new file mode 100755
index 000000000..c16828eab
--- /dev/null
+++ b/2078/CH5/EX5.2/Example5_2.sce
@@ -0,0 +1,21 @@
+//Exa 5.2
+clc;
+clear;
+close;
+//Given data :
+R=0.4;//ohm
+X=0.4;//ohm
+P=2000;//kVA
+pf=0.8;//power factor
+VL=3000;//V
+VR=VL/sqrt(3);//V
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+I=P*1000/3/VR;//A
+VS=VR+I*(R*cos_fi_r+X*sin_fi_r);//V
+Reg=(VS-VR)/VR*100;//%
+disp(Reg,"% Regulation");
+LineLoss=3*I^2*R/1000;//kW
+Pout=P*cos_fi_r;//kW
+Eta_T=Pout*100/(Pout+LineLoss);//%
+disp(Eta_T,"Transmission Efficiency(%)");
diff --git a/2078/CH5/EX5.20/Example5_20.sce b/2078/CH5/EX5.20/Example5_20.sce
new file mode 100755
index 000000000..74146337e
--- /dev/null
+++ b/2078/CH5/EX5.20/Example5_20.sce
@@ -0,0 +1,36 @@
+//Exa 5.20
+clc;
+clear;
+close;
+//Given data :
+Z=200*expm(%i*%pi/180*80);//ohm
+Y=0.0013*expm(%i*%pi/180*90);//mho/phase
+P=80*10^6;//W
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+VRL=220*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+f=50;//Hz
+IR=P/(sqrt(3)*VRL*pf);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+A=1+1/2*Y*Z;//parameter of 3-phase line
+D=A;//parameter of 3-phase line
+B=Z*(1+1/4*Y*Z);//parameter of 3-phase line
+C=Y;//parameter of 3-phase line
+disp("Parameter A, magnitude is "+string(abs(A))+" and angle in degree is "+string(atand(imag(A),real(A))));
+disp("Parameter B, magnitude is "+string(abs(B))+" and angle in degree is "+string(atand(imag(B),real(B))));
+disp("Parameter C, magnitude is "+string(abs(C))+" and angle in degree is "+string(atand(imag(C),real(C))));
+disp("Parameter D, magnitude is "+string(abs(D))+" and angle in degree is "+string(atand(imag(D),real(D))));
+VS=A*VR+B*IR;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end Line voltage(kV) : ");
+IS=C*VR+D*IR;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fis=cosd(fi_s);//sending end pf
+disp(cos_fis,"Sending end power factor(lag) : ");
+Pin=sqrt(3)*VSL*abs(IS)*cos_fis*10^-6;//MW
+disp(Pin,"Power Input(MW) : ");
+Eta=P/(Pin*10^6)*100;//%
+disp(Eta,"Transmission Efficiency(%) : ");
diff --git a/2078/CH5/EX5.21/Example5_21.sce b/2078/CH5/EX5.21/Example5_21.sce
new file mode 100755
index 000000000..766f99576
--- /dev/null
+++ b/2078/CH5/EX5.21/Example5_21.sce
@@ -0,0 +1,29 @@
+//Exa 5.21
+clc;
+clear;
+close;
+//Given data :
+P=50*10^6;//VA
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+A=0.98*expm(%i*%pi/180*3);//parameter of 3-phase line
+D=0.98*expm(%i*%pi/180*3);//parameter of 3-phase line
+B=110*expm(%i*%pi/180*75);//parameter of 3-phase line
+C=0.0005*expm(%i*%pi/180*80);//parameter of 3-phase line
+VRL=110*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+IR=P/(sqrt(3)*VRL);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+VS=A*VR+B*IR;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end Line voltage(kV) : ");
+IS=C*VR+D*IR;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fis=cosd(fi_s);//sending end pf
+disp(cos_fis,"Sending end power factor(lag) : ");
+Pin=sqrt(3)*VSL*abs(IS)*cos_fis*10^-6;//MW
+disp(Pin,"Power Input(MW) : ");
+Eta=P*pf/(Pin*10^6)*100;//%
+disp(Eta,"Transmission Efficiency(%) : ");
diff --git a/2078/CH5/EX5.22/Example5_22.sce b/2078/CH5/EX5.22/Example5_22.sce
new file mode 100755
index 000000000..e095b4369
--- /dev/null
+++ b/2078/CH5/EX5.22/Example5_22.sce
@@ -0,0 +1,44 @@
+//Exa 5.22
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+L=300;//km
+r=0.15;//ohm/km
+x=0.5;//ohm/km
+y=3*10^-6;//mho/km
+VRL=220*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+P=200*10^6;//W
+pf=0.85;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+R=r*L;//ohm
+X=x*L;//ohm
+Y=y*L;//mho
+Z=R+%i*X;//ohm
+//part (i)
+A=1+1/2*%i*Y*Z;//parameter of 3-phase line
+D=A;//parameter of 3-phase line
+B=Z;//parameter of 3-phase line
+C=%i*Y*(1+1/4*%i*Y*Z);//parameter of 3-phase line
+disp("Parameter A, magnitude is "+string(abs(A))+" and angle in degree is "+string(atand(imag(A),real(A))));
+disp("Parameter B, magnitude is "+string(abs(B))+" and angle in degree is "+string(atand(imag(B),real(B))));
+disp("Parameter C, magnitude is "+string(abs(C))+" and angle in degree is "+string(atand(imag(C),real(C))));
+disp("Parameter D, magnitude is "+string(abs(D))+" and angle in degree is "+string(atand(imag(D),real(D))));
+//part (ii)
+IR=poly(0,'IR');
+p=0.024525*IR^2+11.427*IR-2102;//from VS=A*VR+B*IR;//Volt
+IR=roots(p);
+IR=IR(2);//taking +ve value
+P=sqrt(3)*VRL*IR*10^-6;//MW
+disp(P,"Power received in MW : ");
+///part (iii)
+P=200*10^6;//W
+IR=P/sqrt(3)/VRL/pf;//A
+fi=acosd(pf);//degree
+IR=IR*expm(%i*-fi*%pi/180);
+VS=A*VR+B*IR;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end Line voltage(kV) : ");
diff --git a/2078/CH5/EX5.23/Example5_23.sce b/2078/CH5/EX5.23/Example5_23.sce
new file mode 100755
index 000000000..4a48c39d5
--- /dev/null
+++ b/2078/CH5/EX5.23/Example5_23.sce
@@ -0,0 +1,29 @@
+//Exa 5.23
+clc;
+clear;
+close;
+//Given data :
+A=0.936+%i*0.016;//parameter of 3-phase line
+D=A;//parameter of 3-phase line
+B=33.5+%i*138;//parameter of 3-phase line
+C=(-0.9280+%i*901.223)*10^-6;//parameter of 3-phase line
+VRL=200*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+P=40*10^6;//W
+pf=0.86;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+IR=P/sqrt(3)/VRL/pf;//A
+fi=acosd(pf);//degree
+IR=IR*expm(%i*-fi*%pi/180);
+VS=A*VR+B*IR;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end Line voltage(kV) : ");
+IS=C*VR+D*IR;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+fi_s=atand(imag(IS),real(IS))-atand(imag(VS),real(VS));//degree
+disp(cosd(fi_s),fi_s,"Sending end phase angle(degree) & power factor(leading): ");
+Ps=sqrt(3)*abs(VSL)*abs(IS)*cosd(fi_s)*10^-6;//MW
+disp(Ps,"Sending end power(MW) : ");
+Vreg=(VSL-VRL)*100/VRL;//%
+disp(Vreg,"Voltage regulation in % : ");
diff --git a/2078/CH5/EX5.24/Example5_24.sce b/2078/CH5/EX5.24/Example5_24.sce
new file mode 100755
index 000000000..2c01ff723
--- /dev/null
+++ b/2078/CH5/EX5.24/Example5_24.sce
@@ -0,0 +1,31 @@
+//Exa 5.24
+clc;
+clear;
+close;
+//Given data :
+A1=0.98*expm(%i*2*%pi/180);//parameter of 3-phase line
+D1=A1;//parameter of 3-phase line
+B1=28*expm(%i*69*%pi/180);//parameter of 3-phase line
+C1=0.0002*expm(%i*88*%pi/180);//parameter of 3-phase line
+A2=0.95*expm(%i*3*%pi/180);//parameter of 3-phase line
+D2=A2;//parameter of 3-phase line
+B2=40*expm(%i*85*%pi/180);//parameter of 3-phase line
+C2=0.0004*expm(%i*90*%pi/180);//parameter of 3-phase line
+VRL=110*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+IR=200;//A
+pf=0.95;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+fi=acosd(pf);//degree
+A=A1*A2+B1*C2;//generalized parameter of 2 line
+B=A1*B2+B1*D2;//generalized parameter of 2 line
+C=C1*A2+D1*C2;//generalized parameter of 2 line
+D=C1*B2+D1*D2;//generalized parameter of 2 line
+IR=IR*expm(%i*-fi*%pi/180);
+VS=A*VR+B*IR;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+disp(VSL/1000,"Sending end Line voltage(kV) : ");
+IS=C*VR+D*IR;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+//Answer for VSL is wrong in the book.
diff --git a/2078/CH5/EX5.25/Example5_25.sce b/2078/CH5/EX5.25/Example5_25.sce
new file mode 100755
index 000000000..e471219af
--- /dev/null
+++ b/2078/CH5/EX5.25/Example5_25.sce
@@ -0,0 +1,36 @@
+//Exa 5.25
+clc;
+clear;
+close;
+//Given data :
+A1=0.98*expm(%i*1*%pi/180);//parameter of 3-phase line
+D1=A1;//parameter of 3-phase line
+B1=100*expm(%i*75*%pi/180);//parameter of 3-phase line
+C1=0.0005*expm(%i*90*%pi/180);//parameter of 3-phase line
+A2=0.98*expm(%i*1*%pi/180);//parameter of 3-phase line
+D2=A2;//parameter of 3-phase line
+B2=100*expm(%i*75*%pi/180);//parameter of 3-phase line
+C2=0.0005*expm(%i*90*%pi/180);//parameter of 3-phase line
+P=100*10^6;//W
+VRL=132*10^3;//Volt
+VR=VRL/sqrt(3);//Volt
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+fi=acosd(pf);//degree
+A=(A1*B2+A2*B1)/(B1+B2);//generalized parameter of 2 line
+B=B1*B2/(B1+B2);//generalized parameter of 2 line
+C=C1+C2-(A1-A2)*(D1-D2)/(B1+B2);//generalized parameter of 2 line
+D=(B1*D2+B2*D1)/(B1+B2);//generalized parameter of 2 line
+disp("Generalised constants ot two lines combined are : ");
+disp("Parameter A, magnitude is "+string(abs(A))+" and angle in degree is "+string(atand(imag(A),real(A))));
+disp("Parameter B, magnitude is "+string(abs(B))+" and angle in degree is "+string(atand(imag(B),real(B))));
+disp("Parameter C, magnitude is "+string(abs(C))+" and angle in degree is "+string(atand(imag(C),real(C))));
+disp("Parameter D, magnitude is "+string(abs(D))+" and angle in degree is "+string(atand(imag(D),real(D))));
+IR=P/sqrt(3)/VRL/pf;//A
+IR=IR*expm(%i*-fi*%pi/180);
+VS=A*VR+B*IR;//Volt
+VSL=sqrt(3)*abs(VS);//Volt
+IS=C*VR+D*IR;//A
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));
+disp(cosd(fi_s),"Sending end power factor(lagging) : ");
diff --git a/2078/CH5/EX5.3/Example5_3.sce b/2078/CH5/EX5.3/Example5_3.sce
new file mode 100755
index 000000000..ca401c0a1
--- /dev/null
+++ b/2078/CH5/EX5.3/Example5_3.sce
@@ -0,0 +1,23 @@
+//Exa 5.3
+clc;
+clear;
+close;
+//Given data :
+l=15;//km
+P=5;//MW
+V=11;//kV
+f=50;//Hz
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+L=1.1;//mH/Km
+VR=V*1000/sqrt(3);//V
+I=P*1000/sqrt(3)/V/cos_fi_r;//A
+LineLoss=12/100*P*10^6;//W
+R=LineLoss/3/I^2;//ohm
+X=2*%pi*f*L*10^-3*l;//ohm/phase
+VS=VR+I*(R*cos_fi_r+X*sin_fi_r);//V
+VSL=sqrt(3)*VS/1000;//KV
+disp(VSL,"Line voltage at sending end(kV)");
+Reg=(VSL-V)/V*100;//%
+disp(Reg,"% Regulation");
diff --git a/2078/CH5/EX5.4/Example5_4.sce b/2078/CH5/EX5.4/Example5_4.sce
new file mode 100755
index 000000000..3550f6f15
--- /dev/null
+++ b/2078/CH5/EX5.4/Example5_4.sce
@@ -0,0 +1,31 @@
+//Exa 5.4
+clc;
+clear;
+close;
+//Given data :
+l=50;//km
+S=10000;//kVA
+pf=0.8;//power factor
+d=1.2*100;//cm
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+V=33000;//Volts
+VR=V/sqrt(3);//V
+f=50;//Hz
+I=S*1000/sqrt(3)/V;//A
+LineLoss=10/100*S*10^3*pf;//W
+R=LineLoss/3/I^2;//ohm
+rho=1.73*10^-6;//kg/m^3
+a=rho*l*1000*100/R;//cm^2
+r=sqrt(a/%pi);//cm
+L=0.2*log(d/r/0.7788)*l;//mH
+X=2*%pi*f*L*10^-3;//ohm
+VS=VR+I*(R*cos_fi_r+X*sin_fi_r);//V
+VSL=sqrt(3)*VS/1000;//kV
+disp(VSL,"Line voltage at sending end(kV)");
+pf_s=(VR*cos_fi_r+I*R)/VS;//lagging(sendinf end pf)
+disp(pf_s,"Sending end pf(lagging) ");
+Eta_T=S*pf/(S*pf+LineLoss/1000)*100;
+disp(Eta_T,"Transmission Efficiency(%)");
+Reg=(VSL-V/1000)/(V/1000)*100;//%
+disp(Reg,"% Regulation");
diff --git a/2078/CH5/EX5.5/Example5_5.sce b/2078/CH5/EX5.5/Example5_5.sce
new file mode 100755
index 000000000..fdbb263f4
--- /dev/null
+++ b/2078/CH5/EX5.5/Example5_5.sce
@@ -0,0 +1,22 @@
+//Exa 5.5
+clc;
+clear;
+close;
+//Given data :
+VRL=30000;//Volts
+VSL=33000;//Volts
+f=50;//Hz
+P=10*10^6;//W
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+VR=VRL/sqrt(3);//V
+I=P/sqrt(3)/VRL/pf;//A
+Eta_T=0.96;//Efficiency
+LineLoss=P*(1/Eta_T-1);//W
+R=LineLoss/3/I^2;//ohm/phase
+disp(R,"Resistance per phase(ohm/phase)");
+VS=VSL/sqrt(3);//V
+X=(VS-VR-I*R*cos_fi_r)/I/sin_fi_r;//V
+L=X/2/%pi/f;//H/phase
+disp(L*1000,"Inductance per phase(mH/phase)");
diff --git a/2078/CH5/EX5.6/Example5_6.sce b/2078/CH5/EX5.6/Example5_6.sce
new file mode 100755
index 000000000..4a116f3c9
--- /dev/null
+++ b/2078/CH5/EX5.6/Example5_6.sce
@@ -0,0 +1,25 @@
+//Exa 5.6
+clc;
+clear;
+close;
+//Given data :
+l=3;//km
+P=3000;//KW
+VSL=11*10^3;//volt
+R=l*0.4;//ohm
+X=l*0.8;//ohm
+VS=VSL/sqrt(3);//Volts
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+//VS=VR+I*(R*cos_fi_r+X*sin_fi_r);//V
+I_into_VR=P*1000/3/cos_fi_r;//VA
+//VR^2-VS*VR+I_into_VR*(R*cos_fi_r+X*sin_fi_r);
+p=[1 -VS I_into_VR*(R*cos_fi_r+X*sin_fi_r)];
+VR=roots(p);
+VR=VR(1);//taking greater value
+I=I_into_VR/VR;//A
+VRL=sqrt(3)*VR;//volt
+disp(VRL,"Line voltage at load end(volt) : ");
+Eta_T=P*1000/(P*1000+3*I^2*R)*100;//%
+disp(Eta_T,"Transmission Efficiency(%) : ");
diff --git a/2078/CH5/EX5.7/Example5_7.sce b/2078/CH5/EX5.7/Example5_7.sce
new file mode 100755
index 000000000..b1ec8d64c
--- /dev/null
+++ b/2078/CH5/EX5.7/Example5_7.sce
@@ -0,0 +1,19 @@
+//Exa 5.7
+clc;
+clear;
+close;
+//Given data :
+R=5;//ohm/phase
+X=20;//ohm/phase
+VSL=46.85;//kV
+VRL=33;//kV
+VRL=VRL*1000;//v
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+VR=VRL/sqrt(3);//V
+I=(VSL*1000/sqrt(3)-VR)/(R*cos_fi_r+X*sin_fi_r);//A
+Pout=sqrt(3)*VRL*I*pf/1000;//kW
+disp(Pout,"Power output(kW)");
+cosfi_s=(VR*pf+I*R)/(VSL*1000/sqrt(3));//power factor
+disp(cosfi_s,"Power factor at sending end(lagging)");
diff --git a/2078/CH5/EX5.8/Example5_8.sce b/2078/CH5/EX5.8/Example5_8.sce
new file mode 100755
index 000000000..f36b431c7
--- /dev/null
+++ b/2078/CH5/EX5.8/Example5_8.sce
@@ -0,0 +1,30 @@
+//Exa 5.8
+clc;
+clear;
+close;
+//Given data :
+l=80;//km
+P=15;//MW
+VR=66*10^3;//Volt
+R=l*0.3125;//ohm
+X=l*1;//ohm
+Y=l*17.5*10^-6;//S
+pf=0.8;//power factor
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+IR=P*10^6/(VR*pf);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+IC=%i*Y*VR;//A
+IS=IR+IC;//A
+disp("Sending end current(A), magnitude is "+string(abs(IS))+" and angle in degree is "+string(atand(imag(IS),real(IS))));
+VS=VR+IS*(R+%i*X);//volt
+disp("Sending end voltage(V), magnitude is "+string(abs(VS))+" and angle in degree is "+string(atand(imag(VS),real(VS))));
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fis=cosd(fi_s);//sending end pf
+disp(cos_fis,"Sending end power factor(lag) : ");
+Reg=(abs(VS)-VR)/VR*100;//%
+disp(Reg,"Regulation(%) : ");
+LineLoss=abs(IS)^2*R/1000;//kW
+disp(LineLoss,"Line Losses in kW : ");
+Eta_T=P*1000/(P*1000+LineLoss)*100;//%
+disp(Eta_T,"Transmission Efficiency(%) : ");
diff --git a/2078/CH5/EX5.9/Example5_9.sce b/2078/CH5/EX5.9/Example5_9.sce
new file mode 100755
index 000000000..23738c94c
--- /dev/null
+++ b/2078/CH5/EX5.9/Example5_9.sce
@@ -0,0 +1,34 @@
+//Exa 5.9
+clc;
+clear;
+close;
+//Given data :
+l=100;//km
+P=20;//MW
+VRL=66*10^3;//volt
+f=50;//Hz
+R=10;//ohm
+L=111.7*10^-3;//H
+C=0.9954*10^-6;//F
+pf=0.8;//power factor
+X=2*%pi*f*L;//ohm
+Y=2*%pi*f*C;//S
+cos_fi_r=pf;
+sin_fi_r=sqrt(1-cos_fi_r^2);
+VR=VRL/sqrt(3);//volt
+IR=P*10^6/(sqrt(3)*VRL*pf);//A
+IR=IR*(cos_fi_r-%i*sin_fi_r);//A
+Z=R+%i*X;//ohm
+Vdash=VR+1/2*IR*Z;//Volt
+IC=Vdash*%i*Y;//A
+IS=IR+IC;//A
+VS=Vdash+1/2*IS*Z;//Volt
+VSL=abs(VS)*sqrt(3);//Volt
+disp(VSL,"Sending end line voltage(Volt) :");
+Reg=(VSL-VRL)/VRL*100;//%
+disp(Reg,"Regulation(%) : ");
+fi_s=atand(imag(VS),real(VS))-atand(imag(IS),real(IS));//
+cos_fi_s=cosd(fi_s);//sending end pf
+Eta_T=sqrt(3)*VRL*abs(IR)*cos_fi_r/(sqrt(3)*VSL*abs(IS)*cos_fi_s)*100;//%
+disp(Eta_T,"Transmission Efficiency(%) : ");
+//Ans is not accurate in the book.
generated by cgit (git 2.34.1) at 2025-01-24 15:50:40 +0000