summaryrefslogtreecommitdiff
path: root/3640/CH2/EX2.10/Ex2_10.sce
diff options
context:
space:
mode:
Diffstat (limited to '3640/CH2/EX2.10/Ex2_10.sce')
-rw-r--r--3640/CH2/EX2.10/Ex2_10.sce24
1 files changed, 24 insertions, 0 deletions
diff --git a/3640/CH2/EX2.10/Ex2_10.sce b/3640/CH2/EX2.10/Ex2_10.sce
new file mode 100644
index 000000000..4ecab79e6
--- /dev/null
+++ b/3640/CH2/EX2.10/Ex2_10.sce
@@ -0,0 +1,24 @@
+clc
+xd=1 //in ohms per unit
+xd1=0.3 //in ohms per unit
+xd2=0.2 //in ohms per unit
+Td2=0.03 //time in seconds
+Td1=1 //time in seconds
+MVA=100 //rating in mega volt ampere
+V=16000 //voltage in volts
+I2pu=1/xd2
+mprintf("Ipú́=%dper unit\n",I2pu)
+Ib=(MVA*(10^6))/(sqrt(3)*V)
+mprintf("Ib=%fA\n",Ib)//ans may vary due to roundoff error
+mprintf("Í́=%fA\n",I2pu*Ib)//ans in textbook is wrong
+I1=1/xd1 //current in per unit
+mprintf("Í=Efo/xd́=%fper unit\n",I1)//ans may vary due to roundoff error
+Iss=1/xd//current in per unit
+mprintf("Iss=Efo/xd=1 per unit\n")
+t=2/60 //time in seconds
+mprintf("I=%fper unit\n",(I2pu-I1)*exp(-t/Td2)+(I1-Iss)*exp(-t/Td1)+1)//ans may vary due to roundoff error
+t=10 //time in seconds
+mprintf("I=%fper unit\n",(I2pu-I1)*exp(-t/Td2)+(I1-Iss)*exp(-t/Td1)+1)//ans may vary due to roundoff error
+
+
+