initial commit / add all books

1 files changed, 85 insertions, 0 deletions

diff --git a/3733/CH32/EX32.34/Ex32_34.sce b/3733/CH32/EX32.34/Ex32_34.sce
new file mode 100644
index 000000000..167637980
--- /dev/null
+++ b/3733/CH32/EX32.34/Ex32_34.sce
@@ -0,0 +1,85 @@
+// Example 32_34

+clc;funcprot(0);

+//Given data

+L_cap=1500;// kW

+// n=0.43*(L)^0.48;(given)

+T=[0 4 8 12 16 20 24];// Time in hours

+L_a=[200 600 1000 400 200 100];// Load in kW

+L_b=[800 400 200 200 600 400];// Load in kW

+L_t=[1000 1000 1200 600 800 500];// Load in kW

+CV=45*10^3;// MJ/kg

+Dc=30;// The cost of diesel in Rs./liter

+SG=0.85;// Specific gravity

+pr=15/100;// The profit required

+oc=30/100;// The other costs

+n_com=92/100;// Combustion efficiency

+

+// Calculation

+t=[0 4 4 8 8 12 12 16 16 20 20 24 24];// Time in hrs for load curve

+L_A=[200 200 600 600 1000 1000 400 400 200 200 100 100 1500];// Load A in kW for load curve

+subplot(3,1,1);

+xlabel('Time in hours');

+ylabel('Load in kW');

+xtitle('Load of consumer-A');

+plot(t,L_A,'b');

+legend('Load curve for (A)');

+t=[0 0 4 4 8 8 12 12 16 16 20 20 24 24];// Time in hrs for load curve

+L_B=[0 800 800 400 400 200 200 200 200 600 600 400 400 1500];// Load B in kW for load curve

+subplot(3,1,2);

+xlabel('Time in hours');

+ylabel('Load in kW');

+xtitle('Load of consumer-B');

+plot(t,L_B,'b');

+legend('Load curve for (B)');

+t=[0 0 4 4 8 8 12 12 16 16 20 20 24 24];// Time in hrs for load curve

+L_AB=[0 1000 1000 1000 1000 1200 1200 600 600 800 800 500 500 1500];// Load A+B in kW for load curve

+subplot(3,1,3);

+xlabel('Time in hours');

+ylabel('Load in kW');

+xtitle('Load of on plant for consumerA and B');

+plot(t,L_AB,'b');

+legend('Load curve for (A+B)');

+//(i)

+E=(L_a(1)*(T(2)-T(1)))+(L_a(2)*(T(3)-T(2)))+(L_a(3)*(T(4)-T(3)))+(L_a(4)*(T(5)-T(4)))+(L_a(5)*(T(6)-T(5)))+(L_a(6)*(T(7)-T(6)));// Total energy consumed a day in kWh

+L_a1=E/24;// kW

+L_max1=1000;// kW

+LF_A=L_a1/L_max1;

+E=(L_b(1)*(T(2)-T(1)))+(L_b(2)*(T(3)-T(2)))+(L_b(3)*(T(4)-T(3)))+(L_b(4)*(T(5)-T(4)))+(L_b(5)*(T(6)-T(5)))+(L_b(6)*(T(7)-T(6)));// Total energy consumed a day in kWh

+L_b1=E/24;// kW

+L_max2=800;// kW

+LF_B=L_b1/L_max2;

+E=(L_t(1)*(T(2)-T(1)))+(L_t(2)*(T(3)-T(2)))+(L_t(3)*(T(4)-T(3)))+(L_t(4)*(T(5)-T(4)))+(L_t(5)*(T(6)-T(5)))+(L_t(6)*(T(7)-T(6)));// Total energy consumed a day in kWh

+L_ab=E/24;// kW

+L_max=1200;

+LF_AB=L_ab/L_max;

+PL=((LF_AB-LF_A)/LF_A)*100;// Maximum percentage increase in load factor

+DF=(L_max1+L_max2)/L_max;// Diversity factor

+O_1=(L_t(1)*(T(3)-T(1)));// kWh

+n_1=0.43*(L_t(1)/L_cap)^0.48;

+I_1=O_1/n_1;// kWh

+O_2=(L_t(3)*(T(4)-T(3)));// kWh

+n_2=0.43*(L_t(3)/L_cap)^0.48;

+I_2=O_2/n_2;// kWh

+O_3=(L_t(4)*(T(5)-T(4)));// kWh

+n_3=0.43*(L_t(4)/L_cap)^0.48;

+I_3=O_3/n_3;// kWh

+O_4=(L_t(5)*(T(6)-T(5)));// kWh

+n_4=0.43*(L_t(5)/L_cap)^0.48;

+I_4=O_4/n_4;// kWh

+O_5=(L_t(6)*(T(7)-T(6)));// kWh

+n_5=0.43*(L_t(6)/L_cap)^0.48;

+I_5=O_5/n_5;// kWh

+I_t=(I_1+I_2+I_3+I_4+I_5)*3600;// Total input in kJ

+m_f=I_t/(CV*n_com*24);// kg/hr

+V_f=m_f/0.85;// liters/hr

+V_f=V_f*24;// liters

+C_f=V_f*Dc;// Cost of fuel in Rs./day

+Oc=C_f*oc;// The other cost running the plants in Rs./day

+Tc=C_f+Oc;// The total cost running the plants in Rs./day

+Pr=Tc*pr;// The profit required in Rs./day

+Tcs=Tc+Pr;// Total cost of saling the energy generated/day in rupees

+O_t=O_1+O_2+O_3+O_4+O_5;// Total energy generated in kWh

+Cs=Tcs/O_t;// The cost of sailing the energy in Rs./kWh

+printf('\n(i)The individual load factor of consumer A=%0.3f \n   The individual load factor of consumer B=%0.3f \n(ii)Load factor of the system=%0.3f \n    Diversity factor of the system=%0.1f \n(iii)The cost of selling the power=Rs.%0.2f/kWh',LF_A,LF_B,LF_AB,DF,Cs);

+// The answer vary due to round off error

+