blob: ee0f0d9ef10f5491b01a6c67a5228246fbad2866 (
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
|
//example 8.13
//calculate hydrodynamic earthquake pressure
//moment at 50m below water surface
clc;funcprot(0);
//given
c=1;
H=100; //heigth of dam
hw=100; //heigth of water in reservior
FB=1; //free board
s=0.15; //slope of upstream face
gamma_w=9.81; //unit weigth of water
alphah=0.1;
theta=atan(s);
y=50;
Cm=0.735*(1-(theta*2/%pi));
Cy=(Cm/2)*((y*(2-y/hw)/hw)+(y*(2-y/hw)/hw)^0.5);
pe=Cy*alphah*gamma_w*hw;
F=0.726*pe*y;
M=0.299*pe*y^2;
pe=round(pe*1000)/1000;
F=round(F*10)/10;
M=round(M*10)/10;
mprintf("hydrodynamic earthquake pressure=%f kN/square.m\nshear=%f kN/m.\nMoment=%f kN-m/m.",pe,F,M);
|
