5 files changed, 82 insertions, 0 deletions
diff --git a/1730/CH5/EX5.1/Exa5_1.sce b/1730/CH5/EX5.1/Exa5_1.sce
new file mode 100755
index 000000000..f5e3555d4
--- /dev/null
+++ b/1730/CH5/EX5.1/Exa5_1.sce
@@ -0,0 +1,15 @@
+//Exa 5.1
+clc;
+clear;
+close;
+// given data
+Area_hysteresis_curve=9.3;//in cm^2
+Cordinate1_1cm=1000;//in AT/m
+Cordinate2_1cm=0.2;//in T
+//Part (i)
+hysteresis_loss=Area_hysteresis_curve*Cordinate1_1cm*Cordinate2_1cm;//in J/m^3/cycle
+disp(hysteresis_loss,"Hysteresis loss/m^3/cycle in J/m^3/cycle: ");
+//Part (ii)
+f=50;//in Hz
+H_LossPerCubicMeter=hysteresis_loss*f;//in Watts
+disp(H_LossPerCubicMeter*10^-3,"Hysteresis loss Per Cubic Meter in KWatts :");
diff --git a/1730/CH5/EX5.2/Exa5_2.sce b/1730/CH5/EX5.2/Exa5_2.sce
new file mode 100755
index 000000000..242df5a85
--- /dev/null
+++ b/1730/CH5/EX5.2/Exa5_2.sce
@@ -0,0 +1,17 @@
+//Exa 5.2
+clc;
+clear;
+close;
+format('v',11)
+// given data
+Area_hysteresis_loop=93;//in cm^2
+scale1_1cm=0.1;//in Wb/m^2
+scale2_1cm=50;//in AT/m
+
+hysteresis_loss=Area_hysteresis_loop*scale1_1cm*scale2_1cm;//in J/m^3/cycle
+disp(hysteresis_loss,"Hysteresis loss/m^3/cycle in J/m^3/cycle: ");
+
+f=65;//unit less
+V=1500*10^-6;// in m^3
+P_h=hysteresis_loss*f*V;
+disp("Hysteresis loss is : "+string(P_h)+" W");
diff --git a/1730/CH5/EX5.3/Exa5_3.sce b/1730/CH5/EX5.3/Exa5_3.sce
new file mode 100755
index 000000000..65cb90c9d
--- /dev/null
+++ b/1730/CH5/EX5.3/Exa5_3.sce
@@ -0,0 +1,25 @@
+//Exa 5.3
+clc;
+clear;
+close;
+format('v', 11)
+// given data
+nita=628;// in J/m^3
+B_max=1.3;// in Wb/m^2
+f=25;// in Hz
+ironMass=50;// in kg
+densityOfIron=7.8*10^3;// in kg/m^3
+V=ironMass/densityOfIron;
+x=12.5;// in AT/m
+y=0.1;// in T
+// formula Hysteresis loss/second = nita*B_max^1.6*f*V
+H_Loss_per_second = nita*B_max^1.6*f*V ;// in J/s
+H_Loss_per_second=floor(H_Loss_per_second);
+H_Loss_per_hour= H_Loss_per_second*60*60;// in J
+disp("Hysteresis Loss per hour is : "+string(H_Loss_per_hour)+" J");
+// Let Hysteresis Loss per m^3 per cycle = H1
+H1=nita*B_max^1.6;
+// formula  hysteresis loss/m^3/cycle = x*y*area of B-H loop
+Area_of_B_H_loop=H1/(x*y);
+Area_of_B_H_loop=floor(Area_of_B_H_loop);
+disp("Area of B-H loop is : "+string(Area_of_B_H_loop)+" cm^2");
diff --git a/1730/CH5/EX5.4/Exa5_4.sce b/1730/CH5/EX5.4/Exa5_4.sce
new file mode 100755
index 000000000..fe045d5d2
--- /dev/null
+++ b/1730/CH5/EX5.4/Exa5_4.sce
@@ -0,0 +1,12 @@
+//Exa 5.4
+clc;
+clear;
+close;
+// given data
+H_L_per_M_Cube_per_C=380;// in W-S
+f=50;// unit less
+density=7800;// in kg/m^3
+V=1/density;// in m^3
+// formula Hysteresis loss = Hysteresis loss/m^3/cycle * f * V
+P_h=H_L_per_M_Cube_per_C * f * V;
+disp("Hysteresis loss is : "+string(P_h)+" W");
diff --git a/1730/CH5/EX5.5/Exa5_5.sce b/1730/CH5/EX5.5/Exa5_5.sce
new file mode 100755
index 000000000..d17f77547
--- /dev/null
+++ b/1730/CH5/EX5.5/Exa5_5.sce
@@ -0,0 +1,13 @@
+//Exa 5.5
+clc;
+clear;
+close;
+// given data
+P_e1=1600;// in watts
+B_max1=1.2;// in T
+f1=50;// in Hz
+B_max2=1.5;// in T
+f2=60;// in Hz
+// P_e propotional to B_max^2*f^2, so
+P_e2=P_e1*(B_max2/B_max1)^2*(f2/f1)^2
+disp("Eddy current loss is : "+string(P_e2)+" watts");