blob: eb7f4d99b1bf260d6c9609d0c12345e9de5510cc (
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
//example 8.11
//calculate width of base if no tension is to develop
//check the stability if uplift is neglected
clc;funcprot(0);
//given
c=1;
hw=6; //heigth of water in reservior
Bt=1.5; //width of top of dam
H=6; //heigth of the dam
gamma_m=20; //unit weigth of masonary
gamma_w=9.81; //weigth density of water
f=1800; //compressive strength
miu=0.6; //coefficient of friction
//to develop no tension e=b/6;x=b/3.
//hence on solving the relations we get
P=poly([-19.908 1.5 1],'b','c')
wb=roots(P);
//roots are 3.774 and -5.27
//since negative value cannot be taken
wb=3.77;
mprintf("Neglecting the negative value.\nWidth of base is=3.77 m.");
W1=Bt*gamma_m*H;
W2=gamma_m*H*(wb-Bt)/2;
L1=(wb-Bt)+(Bt/2);
L2=(2*(wb-Bt))/3;
M1=W1*L1,
M2=W2*L2;
W3=gamma_w*H^2/2;
L3=hw/3;
M3=W3*L3;
SumW=W1+W2;
SumM=M1+M2-M3;
pn=2*SumW/wb;
pn=round(pn*10)/10;
mprintf("\nMaximum stress=%f kN/square.m.",pn);
mprintf("\nDam is safe against compression");
FOS=miu*SumW/W3;
FOS=round(FOS*1000)/1000;
mprintf("\nFactor of safety against sliding=%f. > 1",FOS);
mprintf("\nDam is safe against sliding.");
|
