diff options
Diffstat (limited to '2123/CH4')
-rwxr-xr-x | 2123/CH4/EX4.1/Exa_4_1.sce | 23 | ||||
-rwxr-xr-x | 2123/CH4/EX4.10/Exa_4_10.sce | 17 | ||||
-rwxr-xr-x | 2123/CH4/EX4.11/Exa_4_11.sce | 14 | ||||
-rwxr-xr-x | 2123/CH4/EX4.12/Exa_4_12.sce | 17 | ||||
-rwxr-xr-x | 2123/CH4/EX4.13/Exa_4_13.sce | 26 | ||||
-rwxr-xr-x | 2123/CH4/EX4.2/Exa_4_2.sce | 23 | ||||
-rwxr-xr-x | 2123/CH4/EX4.3/Exa_4_3.sce | 29 | ||||
-rwxr-xr-x | 2123/CH4/EX4.4/Exa_4_4.sce | 16 | ||||
-rwxr-xr-x | 2123/CH4/EX4.5/Exa_4_5.sce | 18 | ||||
-rwxr-xr-x | 2123/CH4/EX4.6/Exa_4_6.sce | 20 | ||||
-rwxr-xr-x | 2123/CH4/EX4.7/Exa_4_7.sce | 17 | ||||
-rwxr-xr-x | 2123/CH4/EX4.8/Exa_4_8.sce | 15 | ||||
-rwxr-xr-x | 2123/CH4/EX4.9/Exa_4_9.sce | 18 |
13 files changed, 253 insertions, 0 deletions
diff --git a/2123/CH4/EX4.1/Exa_4_1.sce b/2123/CH4/EX4.1/Exa_4_1.sce new file mode 100755 index 000000000..83e634137 --- /dev/null +++ b/2123/CH4/EX4.1/Exa_4_1.sce @@ -0,0 +1,23 @@ +//Example No. 4.1
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+P=30;//KW
+theta1=30;//degree C
+t1=40;//min
+theta2=45;//degree C
+t2=80;//min(t2=2*t1)
+disp("theta=theta_f*(1-exp(-t/T))");
+//Let exp(-t1/T)=a then exp(-t2/T)=a^2
+//theta1/theta2=(1-a)/(1-a^2)
+//a^2*theta1-a*theta2+theta2-theta1=0
+P=[theta1 -theta2 theta2-theta1];//Polynomial for a
+a=roots(P);
+a=a(2);//discarding value 1 as it cant give value of T
+T=-t1/log(a);//min
+disp(T,"Thermal time constant in min : ");
+theta_f=theta1/(1-exp(-t1/T));//degreeC
+disp(theta_f,"Final temperature rise in degree C : ");
diff --git a/2123/CH4/EX4.10/Exa_4_10.sce b/2123/CH4/EX4.10/Exa_4_10.sce new file mode 100755 index 000000000..9b1fe8eed --- /dev/null +++ b/2123/CH4/EX4.10/Exa_4_10.sce @@ -0,0 +1,17 @@ +//Example No. 4.10
+clc;
+clear;
+close;
+format('v',5);
+
+//Given Data :
+T=90;//min
+t=25;//min
+ShortTimeRating=50;//KW
+Eff=80/100;//Efficiency
+//Let full load rating is P KW and Losses=Pc
+//CuLoss=(P/(P*Eff))^2 & alfa=Pc/CuLoss
+alfa=(Eff)^2;//unitless
+S=sqrt(((1+alfa)/(1-exp(-t/T))-alfa));
+ContinuousRating_fl=ShortTimeRating/S;//KW
+disp(ContinuousRating_fl,"Continuous rating of motor in KW : ");
diff --git a/2123/CH4/EX4.11/Exa_4_11.sce b/2123/CH4/EX4.11/Exa_4_11.sce new file mode 100755 index 000000000..49a27a458 --- /dev/null +++ b/2123/CH4/EX4.11/Exa_4_11.sce @@ -0,0 +1,14 @@ +//Example No. 4.11
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+Rating=25;//KW
+T=90;//min
+ts=30;//min
+S=sqrt(1/(1-exp(-ts/T)));
+HalfHourRating=S*Rating;//KW
+disp(HalfHourRating,"Half hour rating of motor in KW : ");
+//Answer wrong in textbook.
diff --git a/2123/CH4/EX4.12/Exa_4_12.sce b/2123/CH4/EX4.12/Exa_4_12.sce new file mode 100755 index 000000000..19256dcc1 --- /dev/null +++ b/2123/CH4/EX4.12/Exa_4_12.sce @@ -0,0 +1,17 @@ +//Example No. 4.12
+clc;
+clear;
+close;
+format('v',8);
+
+//Given Data :
+T=60;//min
+t=20;//min
+ShortTimeRating=300;//W
+Eff=80/100;//Efficiency
+//Let full load rating is P KW and Losses=Pc
+//CuLoss=(P/(P*Eff))^2 & alfa=Pc/CuLoss
+alfa=(Eff)^2;//unitless
+S=sqrt(((1+alfa)/(1-exp(-t/T))-alfa));
+ContinuousRating_fl=ShortTimeRating/S;//KW
+disp(ContinuousRating_fl,"Continuous rating of motor in W : ");
diff --git a/2123/CH4/EX4.13/Exa_4_13.sce b/2123/CH4/EX4.13/Exa_4_13.sce new file mode 100755 index 000000000..42b7ec603 --- /dev/null +++ b/2123/CH4/EX4.13/Exa_4_13.sce @@ -0,0 +1,26 @@ +//Example No. 4.13
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+P=6;//poles
+f=50;//Hz
+MoI=9.5;//Kg-m^2
+Tr=550;//N-m
+S=5/100;//Slip
+Tmax=720;//N-m
+T_LH=1020;//N-m
+th=12;//sec
+Tmin=220;//N-m
+Snl=3/100;//No load slip
+Ns=120*f/P;//rpm
+Nnl=Ns-Ns*Snl;//rpm
+Nrated=Ns-Ns*S;//rpm
+omega_mo=Nnl*2*%pi/60;//rad/s
+omega_mr=Nrated*2*%pi/60;//rad/s
+J=[Tr/(omega_mo-omega_mr)]*[th/log((T_LH-Tmin)/(T_LH-Tmax))];//Kg-m^2
+MoI_flywheel=J-MoI;//Kg-m^2
+disp(MoI_flywheel,"Moment of inertia of flywheel in Kg-m^2 : ");
+//Answer in the book is wrong.
diff --git a/2123/CH4/EX4.2/Exa_4_2.sce b/2123/CH4/EX4.2/Exa_4_2.sce new file mode 100755 index 000000000..fc4ed0b64 --- /dev/null +++ b/2123/CH4/EX4.2/Exa_4_2.sce @@ -0,0 +1,23 @@ +//Example No. 4.2
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+P=30;//KW
+theta1=20;//degree C
+t1=30;//min
+theta2=30;//degree C
+t2=60;//min(t2=2*t1)
+disp("theta=theta_f*(1-exp(-t/T))");
+//Let exp(-t1/T)=x then exp(-t2/T)=x^2
+//theta1/theta2=(1-x)/(1-x^2)
+//x^2*theta1-x*theta2+theta2-theta1=0
+P=[theta1 -theta2 theta2-theta1];//Polynomial for a
+x=roots(P);
+x=x(2);//discarding value 1 as it cant give value of T
+T=-t1/log(x);//min
+disp(T,"Thermal time constant in min : ");
+theta_f=theta1/(1-exp(-t1/T));//degreeC
+disp(theta_f,"Final temperature rise in degree C : ");
diff --git a/2123/CH4/EX4.3/Exa_4_3.sce b/2123/CH4/EX4.3/Exa_4_3.sce new file mode 100755 index 000000000..e3e3ccc8f --- /dev/null +++ b/2123/CH4/EX4.3/Exa_4_3.sce @@ -0,0 +1,29 @@ +//Example No. 4.3
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+P=30;//KW
+theta1=54-30;//degree C
+t1=1;//hour
+theta2=67-30;//degree C
+t2=2;//hour(t2=2*t1)
+disp("theta=theta_f*(1-exp(-t/T))");
+//Let exp(-t1/T)=a then exp(-t2/T)=a^2
+//theta1/theta2=(1-a)/(1-a^2)
+//a^2*theta1-a*theta2+theta2-theta1=0
+P=[theta1 -theta2 theta2-theta1];//Polynomial for a
+a=roots(P);
+a=a(2);//discarding value 1 as it cant give value of T
+T=-t1/log(a);//hour
+theta_f=theta1/(1-exp(-t1/T));//degreeC
+theta_steady=theta_f+30;//degreeC
+disp(theta_steady,"Final steady state temperature in degree C : ");
+disp(T,"Heating time constant in hour : ");
+theta2=theta_f;//degree C
+t=2.7;//hour
+theta=40-30;//degree C
+Tdash=-t/log(theta/theta2);//hour
+disp(Tdash,"Cooling time constant in hour : ");
diff --git a/2123/CH4/EX4.4/Exa_4_4.sce b/2123/CH4/EX4.4/Exa_4_4.sce new file mode 100755 index 000000000..14884d765 --- /dev/null +++ b/2123/CH4/EX4.4/Exa_4_4.sce @@ -0,0 +1,16 @@ +//Example No. 4.4
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+T=110;//min
+Tdash=150;//min
+t=30;//min
+tdash=45;//min
+theta_f=50;//degree C
+//theta=theta_f-(theta_f-theta1)*exp(-t/T)
+//theta1=theta*exp(-tdash/Tdash);
+theta=(theta_f-theta_f*exp(-t/T))/(1-exp(-tdash/Tdash)*exp(-t/T));//degreeC
+disp(theta,"Maximum temperature rise of the motor in degree C : ");
diff --git a/2123/CH4/EX4.5/Exa_4_5.sce b/2123/CH4/EX4.5/Exa_4_5.sce new file mode 100755 index 000000000..c66d55714 --- /dev/null +++ b/2123/CH4/EX4.5/Exa_4_5.sce @@ -0,0 +1,18 @@ +//Example No. 4.5
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+theta1=20;//degreeC
+theta2=28;//degreeC
+dthetaBYdt1=0.08;//degreeC/min
+dthetaBYdt2=0.06;//degreeC/min
+//theta=theta_f-(theta_f-theta1)*exp(-t/T)
+//dtheta/dt=(theta_f-theta)/T
+//dthetaBYdt1/dthetaBYdt2=(theta_f-theta1)/(theta_f-theta2)
+theta_f=(theta2*dthetaBYdt1-theta1*dthetaBYdt2)/(dthetaBYdt1-dthetaBYdt2)
+disp(theta_f,"Final temperature rise in degree C : ");
+T=(theta_f-theta1)/dthetaBYdt1;//min
+disp(T,"Heating time constant in min : ");
diff --git a/2123/CH4/EX4.6/Exa_4_6.sce b/2123/CH4/EX4.6/Exa_4_6.sce new file mode 100755 index 000000000..11a2127a3 --- /dev/null +++ b/2123/CH4/EX4.6/Exa_4_6.sce @@ -0,0 +1,20 @@ +//Example No. 4.6
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+cycle1=50;//hp
+t1=20;//sec
+cycle2=100;//hp
+t2=20;//sec
+cycle3=150;//hp
+t3=10;//sec
+cycle4=120;//hp
+t4=20;//sec
+cycle5=0;//hp
+t5=15;//sec
+hp_rms=sqrt((cycle1^2*t1+cycle2^2*t2+cycle3^2*t3+cycle4^2*t4+cycle5^2*t5)/(t1+t2+t3+t4+t5));//hp
+disp(hp_rms,"hp(rms) for the motor : ");
+disp("We should choose 100hp motor.")
diff --git a/2123/CH4/EX4.7/Exa_4_7.sce b/2123/CH4/EX4.7/Exa_4_7.sce new file mode 100755 index 000000000..4bb51c350 --- /dev/null +++ b/2123/CH4/EX4.7/Exa_4_7.sce @@ -0,0 +1,17 @@ +//Example No. 4.7
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+t_on=15;//min
+t_off=25;//min
+T=100;//min
+Tdash=140;//min
+theta_f=55;//degree C
+
+//theta=theta_f-(theta_f-theta1)*exp(-t/T)
+//theta1=theta*exp(-tdash/Tdash);
+theta_max=theta_f*[1-exp(-t_on/T)]/(1-exp(-(t_off/Tdash+t_on/T)));//degreeC
+disp(theta_max,"Maximum temperature rise in degree C : ");
diff --git a/2123/CH4/EX4.8/Exa_4_8.sce b/2123/CH4/EX4.8/Exa_4_8.sce new file mode 100755 index 000000000..d1d63919e --- /dev/null +++ b/2123/CH4/EX4.8/Exa_4_8.sce @@ -0,0 +1,15 @@ +//Example No. 4.8
+clc;
+clear;
+close;
+format('v',6);
+
+//Given Data :
+Rating=100;//KW
+alfa=0.9;//unitless
+ts=20;//min
+T=100;//min
+S=sqrt((1+alfa)/(1-exp(-ts/T)));
+ShortTimeRating=S*Rating;//KW
+disp(ShortTimeRating,"Short time rating in KW : ");
+//Answer is wrong in the textbook.
diff --git a/2123/CH4/EX4.9/Exa_4_9.sce b/2123/CH4/EX4.9/Exa_4_9.sce new file mode 100755 index 000000000..1565690fb --- /dev/null +++ b/2123/CH4/EX4.9/Exa_4_9.sce @@ -0,0 +1,18 @@ +//Example No. 4.9
+clc;
+clear;
+close;
+format('v',7);
+
+//Given Data :
+T=80;//min
+Tdash=110;//min
+Rating=50;//KW
+ts=15;//min
+S=sqrt(1/(1-exp(-ts/T)));
+ShortTimeRating=S*Rating;//KW
+disp(ShortTimeRating,"Short time rating of motor in KW : ");
+t_off=20;//min
+S=sqrt((1-exp(-(ts/T+t_off/Tdash)))/(1-exp(-(ts/T))))
+DutyRating=S*Rating;//KW
+disp(DutyRating,"Intermittent periodic duty rating in KW : ");
|