blob: 4cea618b9525ff7bdededea44417d0a3798b0503 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
//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)
|