diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /2087/CH14/EX14.29/example14_29.sce | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '2087/CH14/EX14.29/example14_29.sce')
-rwxr-xr-x | 2087/CH14/EX14.29/example14_29.sce | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/2087/CH14/EX14.29/example14_29.sce b/2087/CH14/EX14.29/example14_29.sce new file mode 100755 index 000000000..a2f1439c6 --- /dev/null +++ b/2087/CH14/EX14.29/example14_29.sce @@ -0,0 +1,27 @@ +
+
+//example 14.29
+//design an irrigation canal for given data
+clc;funcprot(0);
+//given
+Q=15; //discharge
+m=1; //critical velocity ratio
+r=5.7; //B/D
+
+D=(Q/(0.55*6.2))^(1/2.64);
+B=D*r;
+R=(B*D+D^2/2)/(B+D*5^0.5);
+Vo=0.55*m*D^0.64;
+//applying kutters formula; V=C(RS)^0.5
+//where C=(23+1/N+0.00155/S)*(R*S)^0.5/(1+(23+0.00155/S)*N/R^0.5);
+//assuming S^0.5=y
+y=poly([-2D-5,1.55D-3,-0.968,67.5],'x','c');
+roots(y);
+//taking real values of y
+S=0.0141937^2;
+B=round(B*100)/100;
+D=round(D*100)/100;
+mprintf("Width of channel section=%f m.",B);
+mprintf("\nDepth of channel section=%f m.",D);
+mprintf("\nBed slope=%f.",S);
+
|