summaryrefslogtreecommitdiff
path: root/3793/CH12
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3793/CH12
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 '3793/CH12')
-rw-r--r--3793/CH12/EX12.1/exp_12_1.sce22
-rw-r--r--3793/CH12/EX12.2/exp_12_2.sce27
-rw-r--r--3793/CH12/EX12.3/exp_12_3.sce18
3 files changed, 67 insertions, 0 deletions
diff --git a/3793/CH12/EX12.1/exp_12_1.sce b/3793/CH12/EX12.1/exp_12_1.sce
new file mode 100644
index 000000000..7aca342f3
--- /dev/null
+++ b/3793/CH12/EX12.1/exp_12_1.sce
@@ -0,0 +1,22 @@
+clear;
+clc;
+Vs=1.05;//in pu
+Vr=1.0;
+//for del 10 degree
+X=1;
+del=10;
+Qs=((Vs^2)/X)-(Vs*Vr*cosd(del)/X);
+Qr=(-(Vr^2)/X)+(Vs*Vr*cosd(del)/X);
+mprintf("for 10 degree Qs = %.3f and Qr = %.3f\n",Qs,Qr);
+del=20;
+Qs=((Vs^2)/X)-(Vs*Vr*cosd(del)/X);
+Qr=(-(Vr^2)/X)+(Vs*Vr*cosd(del)/X);
+mprintf(" for 20 degree Qs = %.3f and Qr = %.3f\n",Qs,Qr);
+del=30;
+Qs=((Vs^2)/X)-(Vs*Vr*cosd(del)/X);
+Qr=(-(Vr^2)/X)+(Vs*Vr*cosd(del)/X);
+mprintf(" for 30 degree Qs = %.3f and Qr = %.3f\n",Qs,Qr);
+del=40;
+Qs=((Vs^2)/X)-(Vs*Vr*cosd(del)/X);
+Qr=(-(Vr^2)/X)+(Vs*Vr*cosd(del)/X);
+mprintf(" for 40 degree Qs = %.3f and Qr = %.3f",Qs,Qr);
diff --git a/3793/CH12/EX12.2/exp_12_2.sce b/3793/CH12/EX12.2/exp_12_2.sce
new file mode 100644
index 000000000..faa4876ab
--- /dev/null
+++ b/3793/CH12/EX12.2/exp_12_2.sce
@@ -0,0 +1,27 @@
+clear;
+clc;
+
+phi=38.74;//in degree
+k=0.8;
+vs=220;
+
+a=gca();
+a.auto_scale="off";
+a.data_bounds=[0,0;.3,1];
+xlabel("Voltage");
+ylabel("Power");
+
+p=[0:.1:1];
+v1=( 0.5-(p*0.8)+((.25-(p*0.8)-(p^2))^(1/2)))^(1/2);
+v2=( 0.5-(p*0.8)-((.25-(p*0.8)-(p^2))^(1/2)))^(1/2);
+//change the axis boundary of x axis to 0.3
+plot(p,v1,p,v2);
+cp=.24;//from graph
+cv=.55;//from graph
+x=60;
+Pr=(cp*(vs^2))/x;
+Vr=cv*vs;
+mprintf("Critical power = %.3f MW and Critical voltage = %.3f KV",Pr,Vr);
+
+
+
diff --git a/3793/CH12/EX12.3/exp_12_3.sce b/3793/CH12/EX12.3/exp_12_3.sce
new file mode 100644
index 000000000..72f0e1563
--- /dev/null
+++ b/3793/CH12/EX12.3/exp_12_3.sce
@@ -0,0 +1,18 @@
+clear;
+clc;
+Pr=.24;
+Vr=.55;
+del=asind(Pr/Vr);
+a=gca();
+a.auto_scale="off";
+a.data_bounds=[0,0.3;1.5,1];
+xlabel("Voltage");
+ylabel(" ReactivePower");
+V=[0:.01:5];
+Q=((V)^2)+(V*cosd(del));
+plot(V,Q);
+q=.38;
+v=.95;
+Qr=((220^2)*q)/60;
+Vrr=v*220;
+mprintf("Critical power = %.3f MVAR and Critical voltage = %.3f KV",Qr,Vrr);