summaryrefslogtreecommitdiff
path: root/3872/CH5/EX5.10/Ex5_10.sce
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3872/CH5/EX5.10/Ex5_10.sce
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3872/CH5/EX5.10/Ex5_10.sce')
-rw-r--r--3872/CH5/EX5.10/Ex5_10.sce28
1 files changed, 28 insertions, 0 deletions
diff --git a/3872/CH5/EX5.10/Ex5_10.sce b/3872/CH5/EX5.10/Ex5_10.sce
new file mode 100644
index 000000000..6b9fd554a
--- /dev/null
+++ b/3872/CH5/EX5.10/Ex5_10.sce
@@ -0,0 +1,28 @@
+// Book - Power System: Analysis & Design 5th Edition
+// Authors - J. Duncan Glover, Mulukutla S. Sharma, Thomas J. Overbye
+// Chapter - 5 : Example 5.10
+// Scilab Version 6.0.0 : OS - Windows
+
+clc;
+clear;
+
+Comp = 30/100; // Compensation in percent
+Vs = 765; // Sending end voltage in kV
+Vr = Vs; // Receiving end voltage in kV
+Z = 97.02; // Absolute eqivalent pi circuit value ; Taken from Ex 5.3
+PRmaxun=5738 // Maximum power that can be delivered by uncompensated line(From example 5.5)
+F1 = sind(87.210); // Eqivalent pi circuit angle ; Taken from Ex 5.3
+X1 = Z*F1; // Eqivalent series reactance without compensation in Ohm ; taken from Ex 5.3
+Zcap = -(%i)*(1/2)*Comp*X1; // Impedance of series capacitor in Ohm
+ABCD = [1 Zcap; 0 1]; // From figure 5.4 for series impedance the ABCD matrix
+ABCD2 = [ 0.9313*exp(%i*0.209*%pi/180) 97.0*exp(%i*87.2*%pi/180);
+ 1.37*10^(-3)*exp(%i*90.06*%pi/180) 0.9313*exp(%i*0.209*%pi/180) ]; // The ABCD parameters taken from Ex 5.2
+ABCDeq = ABCD*ABCD2*ABCD; // The eqivalent ABCD matrix of the compensated line
+Aeq = abs(ABCDeq(1,1)); // Absolute value of the line parameter A
+thetaAeq = atand(imag(ABCDeq(1,1))/real(ABCDeq(1,1))); // Angle value of the line parameter A
+Beq = abs(ABCDeq(1,2)); // Absolute value of the line parameter B
+thetaBeq = atand(imag(ABCDeq(1,2))/real(ABCDeq(1,2))); // Angle value of the line parameter B
+PRmax=(Vs^2/Beq)-(Aeq*Vs^2/Beq)*cosd(thetaBeq-thetaAeq); // maximum power that can be delivered
+dif=(PRmax/PRmaxun)*100-100; // Percentage difference in power delivered between compensated and uncompensated line
+printf('The theoretical maximum power that this compensated line can deliver is %d MW',PRmax)
+printf('\nThe power delivered by compensated line is %.2f percent more than that of uncompensated line',dif)