summaryrefslogtreecommitdiff
path: root/3819/CH3/EX3.19/Ex3_19.sce
blob: 1fd7b5dd8ff164901e42056a4a9cc8b86b4952f1 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
// A Textbook of Fluid Mecahnics and Hydraulic Machines - By R K Bansal
// Chapter 3-Hydrostatic Forces on surfaces
// Problem 3.19
//Data given in the Problem
theta=60;
s=sin(theta/180*%pi);
hh=poly(0,"hh");
AC=hh/s;
H=hh/2;
b=1;
d=AC;
IG=b*d^3/12;
//COP=(IG/(A*H)+H)
//COP=(h/sin(theta/180*%pi)^3/12/(h/sin(theta*%pi/180)/(h/2)+h/2
//We know that COP is equal to (h-3),THAT IS ,the depth of centre of pressure
//hence
function f=F(hh)
f=((hh/s)^3/12*s^2/(hh/s*hh/2))+(hh/2)-(hh-3);
endfunction
hh=100;
y=fsolve(hh,F);
mprintf("The height of water for tipping the gate is %f m",y)