From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2087/CH12/EX12.5/example12_5.sce | 31 +++++++++++++++++++++++++++++++
 1 file changed, 31 insertions(+)
 create mode 100755 2087/CH12/EX12.5/example12_5.sce

(limited to '2087/CH12/EX12.5')

diff --git a/2087/CH12/EX12.5/example12_5.sce b/2087/CH12/EX12.5/example12_5.sce
new file mode 100755
index 000000000..2cc29b121
--- /dev/null
+++ b/2087/CH12/EX12.5/example12_5.sce
@@ -0,0 +1,31 @@
+
+
+//example 12.5
+//calculate uplift pressure at the junction of inner faces of pile with weir floor using Khosla theory
+clc;funcprot(0);
+//given
+b=16;       //total length of floor
+d=5;        //depth of downstream pile
+D=4;        //depth of upstream pile
+H=2.5;      //head created by weir
+
+//pressure at E
+alpha=b/d;
+lambda=(1+(1+alpha^2)^0.5)/2;
+fie=acos((lambda-2)/lambda)/%pi;
+C=19*(D/b)^0.5*((d+D)/b);
+fie=fie*100-C;
+P=H*fie/100;
+P=round(P*1000)/1000;
+mprintf("Pressure at E=%f m.",P);
+
+//pressure at C1
+alpha=b/D;
+lambda=(1+(1+alpha^2)^0.5)/2;
+fie=acos((lambda-2)/lambda)/%pi;
+fic=1-fie;            //by principle reversibility of flow
+C=19*(d/b)^0.5*((d+D)/b);
+fic=fic*100+C;
+P=fic*H/100;
+P=round(P*1000)/1000;
+mprintf("\n Pressure at C=%f m.",P);
-- 
cgit