summaryrefslogtreecommitdiff
path: root/2132/CH3/EX3.11/Example3_11.sce
blob: efce44a8d659159345affd6d3cb1ef0ec605745e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
//Example 3.11
clc;
clear;
close;
format('v',6);
//Given data :
g=9.81;//gravity
D=4;//meter
xbar=(10+7)/2;//meter
A=%pi*D^2/4;//meter^2
w=g*1000;//in N/m^3
p=w*A*xbar/10^6;//in MN
disp(p,"Total pressure in MN : ");
BC=3;//meter
AB=4;//mete
sin_theta=BC/AB;
IG=%pi*D^4/64;//in m^4
h_bar=IG/A/xbar*sin_theta^2+xbar;//in meter
disp("Position of centre of pressure is "+string(h_bar)+" meter.");