diff options
Diffstat (limited to '2087/CH14/EX14.35/example14_35.sce')
-rwxr-xr-x | 2087/CH14/EX14.35/example14_35.sce | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/2087/CH14/EX14.35/example14_35.sce b/2087/CH14/EX14.35/example14_35.sce new file mode 100755 index 000000000..56137be93 --- /dev/null +++ b/2087/CH14/EX14.35/example14_35.sce @@ -0,0 +1,28 @@ +
+
+//example 14.35
+//design the distributory using Laecey theory
+clc;funcprot(0);
+//given
+f=0.85; //silt factor
+AR=3600; //area for rabi
+AK=1400; //area for kharif
+delta_r=0.135; //kor depth for rabi
+delta_k=0.19; //kor depth for kharif
+tr=4; //kor period for rabi
+tk=2.5; //kor period for kharif
+Du_r=8.64*tr*7/delta_r; //duty for rabi
+Du_k=8.64*tk*7/delta_k; //duty for kharif
+q_r=AR/Du_r; //discharge for rabi
+q_k=AK/Du_k; //discharge for kharif
+Q=q_r; //since q_r>q_k
+V=(Q*f^2/144)^(1/6);
+A=Q/V;
+P=4.75*(Q)^0.5;
+D=(P-(P^2-6.944*A)^0.5)/3.472;
+S=f^(5/3)/(3340*Q^(1/6));
+P=round(P*100)/100;
+D=round(D*100)/100;
+mprintf("\nBed slope=%f.",S);
+mprintf("\nPerimeter of channel section=%f m.",P);
+mprintf("\nDepth of channel section=%f m.",D);
|