From f35ea80659b6a49d1bb2ce1d7d002583f3f40947 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:38:01 +0530 Subject: updated the code --- 746/CH3/EX3.01/3_01.jpg | Bin 23737 -> 0 bytes 746/CH3/EX3.01/3_01.sce | 21 --------------------- 746/CH3/EX3.03/3_03.sce | 8 -------- 746/CH3/EX3.04/3_04.sce | 40 ---------------------------------------- 746/CH3/EX3.05/3_05.sce | 25 ------------------------- 746/CH3/EX3.06/3_06.sce | 9 --------- 746/CH3/EX3.07/3_07.sce | 18 ------------------ 7 files changed, 121 deletions(-) delete mode 100755 746/CH3/EX3.01/3_01.jpg delete mode 100755 746/CH3/EX3.01/3_01.sce delete mode 100755 746/CH3/EX3.03/3_03.sce delete mode 100755 746/CH3/EX3.04/3_04.sce delete mode 100755 746/CH3/EX3.05/3_05.sce delete mode 100755 746/CH3/EX3.06/3_06.sce delete mode 100755 746/CH3/EX3.07/3_07.sce (limited to '746/CH3') diff --git a/746/CH3/EX3.01/3_01.jpg b/746/CH3/EX3.01/3_01.jpg deleted file mode 100755 index 74d8becc0..000000000 Binary files a/746/CH3/EX3.01/3_01.jpg and /dev/null differ diff --git a/746/CH3/EX3.01/3_01.sce b/746/CH3/EX3.01/3_01.sce deleted file mode 100755 index 089d4a614..000000000 --- a/746/CH3/EX3.01/3_01.sce +++ /dev/null @@ -1,21 +0,0 @@ -//liquid level// -pathname=get_absolute_file_path('3.01.sce') -filename=pathname+filesep()+'3.01-data.sci' -exec(filename) -//Tube diameter(in mm): -D=1:25; -D1=D/1000 -[m n]=size(D1) -for i=1:n -//Change in liquid level for water(in mm): -dhw(i)=4*STw*cosd(thetaw)/dw/g/D1(i); -//Change in liquid level for mercury(in mm): -dhm(i)=4*STm*cosd(thetam)/dm/g/D1(i); -end; - -//Plotting tube daimeter and water level: -plot(D1*1000,dhw,'-o') -//Plotting tube daimeter and mercury level: -plot(D1*1000,dhm,'-*') -legend(['Water';'Mercury']); -xtitle('Liquid level vs Tube diameter','Liquid level(in mm)','Tube diameter(in mm)') diff --git a/746/CH3/EX3.03/3_03.sce b/746/CH3/EX3.03/3_03.sce deleted file mode 100755 index 92dc4e6bd..000000000 --- a/746/CH3/EX3.03/3_03.sce +++ /dev/null @@ -1,8 +0,0 @@ -//pressure difference// -pathname=get_absolute_file_path('3.03.sce') -filename=pathname+filesep()+'3.03-data.sci' -exec(filename) -//Pressure difference(in lbf/in^2): -dp=g*d*(-d1+SGm*d2-SGo*d3+SGm*d4+d5)/12/144 -printf("\n\nRESULTS\n\n") -printf("\n\nPressure difference between A and B: %f lbf/in^2\n\n",dp) diff --git a/746/CH3/EX3.04/3_04.sce b/746/CH3/EX3.04/3_04.sce deleted file mode 100755 index 4030fcd51..000000000 --- a/746/CH3/EX3.04/3_04.sce +++ /dev/null @@ -1,40 +0,0 @@ -//temperature and pressure// -pathname=get_absolute_file_path('3.04.sce') -filename=pathname+filesep()+'3.04-data.sci' -exec(filename) -//Assuming temperature varies linearly with altitude: -//Temperature gradient(in F/ft): -m=(T1-T2)/(z2-z1) -//Value of g/(m*R): -v=g/m/R/32.2 -//Pressure at Vail Pass(in inches of Hg): -p12=p1*((T2+460)/(T1+460))^v -//Percentage change in density: -pc1=(p12/p1*(T1+460)/(T2+460)-1)*100 -//Assuming density is constant: -//Pressure at Vail Pass(in inches of Hg): -p22=p1*(1-(g*(z2-z1)/(R*32.2)/(T1+460))) -//Percentage change in density: -pc2=0; -//Assuming temperature is constant: -//Pressure at Vail Pass(in inches of Hg): -p32=p1*%e^(-g*(z2-z1)/(R*32.2)/(T2+460)) -//Percentage change in density: -pc3=(p32/p1*(T1+460)/(T1+460)-1)*100 -//For an adiabatic atmosphere: -p42=p1*((62+460)/(80+460))^(k/(k-1)) -//Percentage change in density: -pc4=(p42/p1*(T1+460)/(T2+460)-1)*100 -printf("\n\nRESULTS\n\n") -printf("\n\n1) If temperature varies linearly with altitude\n\n") -printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p12) -printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc1) -printf("\n\n2) If density is constant\n\n") -printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p22) -printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc2) -printf("\n\n3) If temperature is constant\n\n") -printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p32) -printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc3) -printf("\n\n4) For an adiabatic atmosphere\n\n") -printf("\n\n\tAtmospheric pressure at Vail Pass: %f inches of Hg\n\n",p42) -printf("\n\n\tPercentage change in density wrt Denver: %f percent\n\n",pc4) diff --git a/746/CH3/EX3.05/3_05.sce b/746/CH3/EX3.05/3_05.sce deleted file mode 100755 index 4cea618b9..000000000 --- a/746/CH3/EX3.05/3_05.sce +++ /dev/null @@ -1,25 +0,0 @@ -//force and pressure// -pathname=get_absolute_file_path('3.05.sce') -filename=pathname+filesep()+'3.05-data.sci' -exec(filename) -//Net force on the gate(in kN): -Fr=d*g*w*(D*L+L^2/2*sind(theta)) -//Centre of pressure: -//Calculation for y coordinate: - yc=D/sind(theta)+L/2 - //Area(in m^2): - A=L*w - //Moment of inertia of rectangular gate(in m^4): - Ixx=w*L^3/12 - //y coordinate(in m): - y=yc+Ixx/A/yc -//Calculation for x coordinate: - Ixy=0 - xc=w/2 - //x coordinate(in m): - x=xc+Ixy/A/xc -printf("\n\nRESULTS\n\n") -printf("\n\nNet force on the gate: %f kN\n\n",Fr/1000) -printf("\n\nCoordinate of centre of pressure:(%0.1f,%0.1f)\n\n",x,y) - - diff --git a/746/CH3/EX3.06/3_06.sce b/746/CH3/EX3.06/3_06.sce deleted file mode 100755 index 89f2daf73..000000000 --- a/746/CH3/EX3.06/3_06.sce +++ /dev/null @@ -1,9 +0,0 @@ -//force// -pathname=get_absolute_file_path('3.06.sce') -filename=pathname+filesep()+'3.06-data.sci' -exec(filename) -//Force required to keep the door shut(in lbf): -function y=f(z),y=b/L*p0*z+d*b/L*(L*z-z^2),endfunction -Ft=intg(0,L,f) -printf("\n\nRESULTS\n\n") -printf("\n\nForce requiredto kep the door shut: %.1f lbf\n\n",Ft) diff --git a/746/CH3/EX3.07/3_07.sce b/746/CH3/EX3.07/3_07.sce deleted file mode 100755 index 8624fca4a..000000000 --- a/746/CH3/EX3.07/3_07.sce +++ /dev/null @@ -1,18 +0,0 @@ -//force at equilibrium// -pathname=get_absolute_file_path('3.07.sce') -filename=pathname+filesep()+'3.07-data.sci' -exec(filename) -//Horizontal component of resultant force(in kN): -Frh=0.5*d*g*w*D^2 -//Line of action of Frh(in m): -y1=0.5*D+w*D^3/12/(0.5*D)/(w*D) -//Vertical component of resultant force(in kN): -function y=q(x), y=d*g*w*(D-sqrt(a*x)),endfunction -Frv=intg(0,D^2/a,q) -//Line of acion of Frv(in m): -function k=f(x), k=d*g*w/Frv*x*(D-sqrt(a*x)),endfunction -xa=intg(0,D^2/a,f) -//Force required to keep the gate in equilibrium(in kN): -Fa=1/l*(xa*Frv+(D-y1)*Frh) -printf("\n\nRESULTS\n\n") -printf("\n\nForce required to keep the gate at equilibrium: %f kN\n\n",Fa/1000) -- cgit