diff options
Diffstat (limited to '2087/CH15/EX15.3/example15_3.sce')
| -rwxr-xr-x | 2087/CH15/EX15.3/example15_3.sce | 83 |
1 files changed, 83 insertions, 0 deletions
diff --git a/2087/CH15/EX15.3/example15_3.sce b/2087/CH15/EX15.3/example15_3.sce new file mode 100755 index 000000000..05d5789f5 --- /dev/null +++ b/2087/CH15/EX15.3/example15_3.sce @@ -0,0 +1,83 @@ +
+
+//example 15.3
+//design and prepare the longitudnal section;schedule of area statistics and channel dimension of irrigation channel
+clc;funcprot(0);
+//given
+dl=157.7; //datum level
+fsl=157; //full supply level of parent channel
+bl=156; //bed level of parent channel
+kor_r=4; //kor period of rabi
+kor_k=2.5; //kor period of kharif
+kord_r=13.4; //kor depth of rabi
+kord_k=19; //kor depth of kharif
+s=0.5; //side slope
+m=1; //critical velocity ratio
+N=0.0225; //Kutter n
+qo_r=8.64*7*kor_r*100/kord_r; //outlet discharge for rabi(calculation is wrong in book)
+qo_k=8.64*7*kor_k*100/kord_k; //outlet discharge for kharif(calculation is wrong in book)
+ca=16000; //culturable commanded area
+Ir=0.3; //intensity of irrigation in rabi
+Ik=0.125; //intensity of irrigation in rabi
+Ar=Ir*ca; //area under rabi
+Ak=ca*Ik; //area under kharif
+q_r=Ar/qo_r;
+q_k=Ak/qo_k;
+q_r=round(q_r*100)/100;
+q_k=round(q_k*100)/100;
+mprintf("discharge neede for rabi crop=%f cumecs.",q_r);
+mprintf("\ndischarge neede for kharif crop=%f cumecs.",q_k);
+mprintf("\noutlet discharge factor adopted=%i hectares per cumecs.",qo_r);
+//at km 5
+ca=8000; //culturable area
+Ar=Ir*ca; //area under rabi
+q_r=Ar/qo_r;
+l=0.5 //total loss after 5 km
+q=q_r+l; //total discharge
+dq=1.1*q; //desigm discharge
+S=1/4000; //slope
+B=[5.5 4.9 4.55]; //Bed width
+D=[0.73 0.79 0.84]; //water depth
+Vo=[0.448 0.472 0.488]; //critical velocity
+mprintf("\n\nBed width water depth area velocity critical velocity C.V.R");
+for i=1:3
+ A(i)=B(i)*D(i)+D(i)^2/2;
+ V(i)=dq/A(i);
+ m(i)=V(i)/Vo(i);
+ A(i)=round(A(i)*100)/100;
+ V(i)=round(V(i)*1000)/1000;
+ m(i)=round(m(i)*100)/100;
+ mprintf("\n%f %f %f %f %f %f",B(i),D(i),A(i),V(i),Vo(i),m(i));
+end
+B=4.55;D=0.84;
+mprintf("\nhence take B=%f .; D=%f m.",B,D);
+//at km 4
+q=round(q*100)/100;
+mprintf("\ndischarge at 5 km=%f cumecs.",q);
+ca=10000; //culturable area
+Ar=Ir*ca; //area under rabi
+q_r=Ar/qo_r;
+l=0.5 //total loss below 5 km
+P=B+D*5^0.5; //wetted perimeter
+l1=P*1000*2/1000000; //loss between 5 km and 4km
+l2=l1+l;
+q=q_r+l2;
+dq=1.1*q;
+q=round(q*1000)/1000;
+mprintf("\ndischarge at 4 km =%f cumecs",q);
+mprintf("\nother discharge are calculated and are tabulated as:");
+x=[0:1:5];
+A1=[4800 4200 3600 3300 3000 2400];
+A2=[2000 1750 1500 1375 1250 1000];
+S=[22.5 22.5 22.5 24 24 25];
+B=[5.5 5.2 4.85 4.7 4.55 4.55];
+D=[1.04 1.007 0.975 0.945 0.915 0.840];
+dq=[3.56 3.17 2.8 2.6 2.4 2.02];
+V=[0.570 0.555 0.538 0.530 0.521 0.484];
+m=[1.015 1 1 1 1 0.992];
+mprintf("\n\nBelow km area to irrigate rabi area to irrigate kharif bed slope bed width water depth design discharge velocity C.V.R");
+for i=1:6
+ mprintf("\n%i %i %i %f %f %f %f %f %f",x(i),A1(i),A2(i),S(i),B(i),D(i),dq(i),V(i),m(i));
+end
+
+
|