From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1187/CH2/EX2.1/1.sce | 21 +++++++++++++++++++++
 1187/CH2/EX2.2/2.sce | 39 +++++++++++++++++++++++++++++++++++++++
 1187/CH2/EX2.3/3.sce | 31 +++++++++++++++++++++++++++++++
 1187/CH2/EX2.4/4.sce | 26 ++++++++++++++++++++++++++
 1187/CH2/EX2.5/5.sce | 22 ++++++++++++++++++++++
 1187/CH2/EX2.6/6.sce | 27 +++++++++++++++++++++++++++
 6 files changed, 166 insertions(+)
 create mode 100755 1187/CH2/EX2.1/1.sce
 create mode 100755 1187/CH2/EX2.2/2.sce
 create mode 100755 1187/CH2/EX2.3/3.sce
 create mode 100755 1187/CH2/EX2.4/4.sce
 create mode 100755 1187/CH2/EX2.5/5.sce
 create mode 100755 1187/CH2/EX2.6/6.sce

(limited to '1187/CH2')

diff --git a/1187/CH2/EX2.1/1.sce b/1187/CH2/EX2.1/1.sce
new file mode 100755
index 000000000..86f554e5d
--- /dev/null
+++ b/1187/CH2/EX2.1/1.sce
@@ -0,0 +1,21 @@
+clc
+d=1.5; // m
+m=1.2; // kg
+rate=0.0065; // K/m
+R=287; // J/(kg.K)
+T_0=288.15; // K
+p_0=101*10^3; // Pa
+g=9.81; // m/s^2
+
+rho=m/(%pi*d^3/6);
+rho_0=p_0/R/T_0;
+
+// log(rho/rho_0)=(g/R*rate - 1)*log((T_0-rate*z)/T_0)
+
+z=1/rate*(T_0-T_0*exp(log(rho/rho_0)/(g/R/rate-1)));
+
+disp("The height above sea level to which the ballon will rise = ")
+disp(z)
+disp("m")
+
+printf("The height above sea level to which the ballon will rise = %f m", z)
\ No newline at end of file
diff --git a/1187/CH2/EX2.2/2.sce b/1187/CH2/EX2.2/2.sce
new file mode 100755
index 000000000..b57396015
--- /dev/null
+++ b/1187/CH2/EX2.2/2.sce
@@ -0,0 +1,39 @@
+clc
+d=2; // m
+a=1; // radius in m
+rho=880; // density of oil in kg/m^3
+g=9.81; // m/s^2
+rho_w=1000; // density of water in kg/m^3
+
+C_0=4*a/3/%pi; // centroid of the upper semicircle
+h1=a-C_0; // distance of the centroid from the top
+
+P1=rho*g*h1; // Pressure of the oil at this point
+F1=P1*%pi*a^2/2; // Force exerted by the oil on the upper half of the wall
+
+cp1=a^4*(%pi/8-8/(9*%pi)); // (AK^2)_C
+
+cp2=cp1/(%pi*a^2/2*h1); // Centre of Pressure below the centroid
+
+cp0=cp2+h1; // Centre of Pressure below the top
+
+P_w=(rho*g*a)+(rho_w*g*C_0);
+F_w=P_w*%pi*a^2/2;
+
+h2=C_0+rho/rho_w;
+cp2_w=cp1/(%pi*a^2/2*h2);
+cp0_w=a+C_0+cp2_w; // below the top of cylinder
+
+F_total=F1+F_w;
+
+// F1*cp0 + F_w*cp0_w = F_total*x
+
+x=(F1*cp0 + F_w*cp0_w)/F_total;
+
+disp("Total force =")
+disp(F_total)
+disp("N")
+
+disp("Distance of line of action of total force from top of cylinder =")
+disp(x)
+disp("m")
\ No newline at end of file
diff --git a/1187/CH2/EX2.3/3.sce b/1187/CH2/EX2.3/3.sce
new file mode 100755
index 000000000..09d889c0b
--- /dev/null
+++ b/1187/CH2/EX2.3/3.sce
@@ -0,0 +1,31 @@
+clc
+
+rho=1000; // kg/m^3
+g=9.81; // m/s^2
+r=4; // m
+h=2; // m
+l=5; // m
+theta=%pi/6;
+
+A=h*l;
+
+F_h=rho*g*h*A; // Horizontal force
+
+C0=(2^2/(12*2))+2; // distance of line of action below the free surface
+
+AB=4-4*cos(theta);
+
+F_v=rho*g*l*(AB*1+%pi*r^2*theta/(2*%pi)-1/2*h*r*cos(theta));
+BC=0.237; // m
+
+F_net=sqrt(F_h^2+F_v^2);
+
+phi=atand(F_v/F_h);
+
+disp("Net force =")
+disp(F_net)
+disp("N")
+
+disp("Angle between net force and horizontal =")
+disp(phi)
+disp("degrees")
\ No newline at end of file
diff --git a/1187/CH2/EX2.4/4.sce b/1187/CH2/EX2.4/4.sce
new file mode 100755
index 000000000..20be00e1e
--- /dev/null
+++ b/1187/CH2/EX2.4/4.sce
@@ -0,0 +1,26 @@
+clc
+
+m=10; // kg
+M=80; // kg
+H=1.5; // m
+rho=1026; // kg/m^3
+g=9.81; // m/s^2
+d=1; // m
+
+// m*H + M*H/2 =(M+m)(OG)
+
+OG=(m*H + M*H/2)/(M+m);
+
+// For vertical equilibrium, buoyancy = weight
+h=(M+m)/(rho*%pi/4*d^2);
+
+BM=(%pi*d^4/64)/(%pi*d^2*h/4);
+OB=h/2;
+
+GM=OB+BM-OG;
+
+disp("GM =")
+disp(GM)
+disp("m")
+
+disp("Since this is negative (i.e. M is below G) buoy is unstable.")
\ No newline at end of file
diff --git a/1187/CH2/EX2.5/5.sce b/1187/CH2/EX2.5/5.sce
new file mode 100755
index 000000000..95c96b8bd
--- /dev/null
+++ b/1187/CH2/EX2.5/5.sce
@@ -0,0 +1,22 @@
+clc
+m=10; // kg
+M=80; // kg
+OG=0.8333; // m
+rho=1026; // kg/m^3
+g=9.81; // m/s^2
+d=1; // m
+W=(m+M)*g;
+
+// W(OG) = (W + F)(OB + BM) = rho*g*%pi/4*d^2*h1*(h1/2+d^2/(16*h1))
+
+h1=sqrt(2*(W*OG/(rho*g*%pi/4*d^2) - d^2/16));
+
+F=rho*g*%pi/4*d^2*h1 - W;
+
+disp("Least vertical downward force =")
+disp(F)
+disp("N")
+
+disp("Depth of immersion =")
+disp(h1)
+disp("m")
diff --git a/1187/CH2/EX2.6/6.sce b/1187/CH2/EX2.6/6.sce
new file mode 100755
index 000000000..7156c4856
--- /dev/null
+++ b/1187/CH2/EX2.6/6.sce
@@ -0,0 +1,27 @@
+clc
+
+a=5; // m/s^2
+s=0.5; // m
+phi=atand(1/4); // degrees
+g=9.81; // m/s^2
+rho=880; // kg/m^3
+
+a_x=a*cosd(phi); // Horizontal component of acceleration
+a_z=a*sind(phi); // Vertical component of acceleration
+
+theta=atand(a_x/(a_z+g)); // b=tan(theta)
+
+d=(tand(phi)+tand(theta))/(1-tand(phi)*tand(theta));
+
+c=s*d;
+
+V=s*(s^2-s*c/2);
+
+disp("(a) Volume left in the tank =")
+disp(V*1000)
+disp("L")
+
+P=rho*g*s*cosd(phi);
+disp("(b)Pressure at the lowest corners of the tank =")
+disp(P)
+disp("Pa")
\ No newline at end of file
-- 
cgit