Fl! Hp Prime - Backup Prog. Diseño de Canales 17.07.15 - 10.38 Pm

download Fl! Hp Prime - Backup Prog. Diseño de Canales 17.07.15 - 10.38 Pm

If you can't read please download the document

description

Hp Prime

Transcript of Fl! Hp Prime - Backup Prog. Diseño de Canales 17.07.15 - 10.38 Pm

EXPORT Hidrulica_de_Canales()BEGINPRINT;// VARIABLES A DECLARARLOCAL b1,b2;LOCAL y,y33;LOCAL g,,Dd;LOCAL u,u2;LOCAL uu,uu2;LOCAL y2,y3,y4,y5;LOCAL y22,yc1,yc2,yc3;LOCAL A2,A3,A4,A5,Aq2;LOCAL A22,Ac1,Ac2,Ac3;LOCAL P2,P3,P4,P5,Pq2,Pq22;LOCAL P22,Pc1,Pc2,Pc22,Pc3;LOCAL R2,R3,R4,R5,Rc3,Rq22;LOCAL R22,Rc1,Rc2,Rc22,Rq2;LOCAL T1,T3,Ba,T4,T5,Tc3;LOCAL k2,Tc1,kc2,kq2;LOCAL k22,Sc1,Sc2,Sc22,Sc3;LOCAL V2,V3,V4,V5,Vc3,Vq22;LOCAL V22,Vc1,Vc2,Vq2,Fq2;LOCAL E2,E3,E4,E5,Em3;LOCAL E22,Em1,Em2,Eq2,Eq22;LOCAL F2,F3,F4,F5,Fq22,Eq3;LOCAL F22,Fc1,Fc2,Fc3,Fq3;LOCAL Q1,Aq1,Pq1,Rq1,Vq1,Tq1;LOCAL Eq1,Fq1,Q2,Q22,uq;LOCAL Q3,Aq3,Pq3,Rq3,Tq3,Vq3;LOCAL b1,Ab,Pb,Rb,Tb,Vb,Eb;LOCAL Fb,Ss,As,Ps,Rs,Ts,Vs;LOCAL Es,Fs,Nn,An,Pn,Rn,Tn;LOCAL Vn,En,Fn,yd,Ad,Pd,Rd;LOCAL Td,Vd,Ed,Fd,Ac,Pc,Rc;LOCAL Sc,Tc,Vc,Ec,Fc,Ass,Pss;LOCAL Rss,Nss,Tss,Vss,Ess,Fss;// VARIABLES A DECLARAR#########INPUT( {{W, {" TIRANTE NORMAL"," TIRANTE CRITICO"," CAUDALES"," OTRAS INTERROGANTES"} }},"QUE DESEA REALIZAR?",{" Hallar "},{"Emp c mos a Dis ar :3"});// _____________________________________________________________________________________________________IF W==1 THENINPUT( {{C, {" Trap, R ct. o Trian."," Parablica"," Circular"," S ccin d Mxima Efici ncia"," S ccin d Mnima Infiltracin"} }},"TIRANTE NORMAL",{" S ccin "},{"Con qu D s a Trabajar? *.*"});IF C==1 THENINPUT({Q,B,Z,N,S},"DATOS P/SECC TRAP, RECT O TRIAN",{"Caudal 'Q' ","Sol ra 'b' ","Talud 'Z' ","Rugosidad 'n' ","P ndi nt 'S' "},{"Caudal n m tros cbicos/s gundo","Ancho d Sol ra n m tros","Talud d Dis o","Co fici nt d Rugosidad","P ndi nt n m tros/m tros"});END;IF C==2 THENINPUT({Q,T,N,S},"DATOS P/SECC PARABOLICA",{"Caudal 'Q' ","Esp jo 'T' ","Rugosidad 'n' ","P ndi nt 'S' "},{"Caudal n m tros cbicos/s gundo","Esp jo d Agua n m tros","Co fici ntosidad","P ndi nt n m tros/m tros"});d RugEND;IF C==3 THENINPUT({Q,D,N,S},"DATOS P/SECC CIRCULAR",{"Caudal 'Q' ","Dim tro 'd' ","Rugosidad 'n' ","P ndi nt 'S' "},{"Caudal n m tros cbicos/s gundo","Dim tro n m tros","Co fici nt d Rugosidad","P ndi nt n m tros/m tros"});END;IF C==4 THENINPUT({Q,Z,N,S},"DATOS P/SECC MAX EFICIENCIA",{"Caudal 'Q' ","Talud 'Z' ","Rugosidad 'n' ","P ndi nt 'S' "},{"Caudal n m tros cbicos/s gundo","Talud d Dis o | TRAP=3/3, TRIAN=1","Co fici ntd Rugosidad","P ndi nt n m tros/m tros"});END;IF C==5 THENINPUT({Q,Z,N,S},"DATOS P/SECC MIN INFILTRACION",{"Caudal 'Q' ","Talud 'Z' ","Rugosidad 'n' ","P ndi nt 'S' "},{"Caudal n m tros cbicos/s gundo","Talud d Dis o","Co fici nt d Rugosidad","P ndi nt n m tros/m tros"});END;// ECUACION PARA CALCULO DE TIRANTES EN SECCION TRAPEZOIDAL, RECTANGULAR Y TRIANGULARIF C==1 THENy:=SOLVE(Q=((B*Y+Z*(Y^2))*(((B*Y+Z*(Y^2))/(B+(2*Y*((1+(Z^2))^(1/2)))))^(2/3))*(S^(1/2)))/N,Y); //Tirant P/S cc Trap, R ct y Triang:=9.81;A:=B*y+Z*(y^2); //Ar a P/S cc Trap, R ct y TrianP:=B+2*y*((1+(Z^2))^(1/2)); //P rm tro P/S cc Trap, R ct y TrianR:=A/P; //Radio hidrulico P/S cc Trap, R ct y TrianT1:=B+2*Z*y; //Esp jo d Agua P/S cc Trap, R ct y TrianV:=Q/A; //V locidad P/S cc Trap, R ct y TrianE:=y+((V^2)/(2*g)); //En rga P/S cc Trap, R ct y TrianF:=V/((g*(A/T1))^(1/2)); //# d Froud P/S cc Trap, R ct y Trian// R sulados para S ccion s TRAPEZOIDALES, RECTANGULARES O TRIANGULARESPRINT("RESULTADOS");PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("********");Tirant Normal 'y' = "+ROUND(y,2)+" m");Ar a Hidrulica 'A' = "+ROUND(A,2)+" m2");P rm tro 'P' = "+ROUND(P,2)+" m");Radio Hidrulico 'R' = "+ROUND(R,2)+" m");Esp jo d Agua 'T' = "+ROUND(T1,2)+" m");V locidad 'V' = "+ROUND(V,2)+" m/s");En rga Esp cfica 'E' = "+ROUND(E,3)+" m-Kg/Kg");# d Froud 'F' = "+ROUND(F,3));CASEIF F < 1 THENPRINT(" * Tipo d Flujo = Subcrtico");END;IF F > 1 THENPRINT(" * Tipo d Flujo = Sup rcrtico");END;DEFAULTPRINT(" * Tipo d Flujo = Crtico");END;END;// R sulados para S ccion s TRAPEZOIDALES, RECTANGULARES O TRIANGULARES###############// ECUACION PARA CALCULO DE TIRANTES EN SECCION PARABOLICAIF C==2 THENy2:=SOLVE(Q=(((2*T*Y)/3)*(((2*(T^2)*Y)/(3*(T^2)+8*(Y^2)))^(2/3))*(S^(1/2)))/N,Y);y22:=SOLVE(Q=((((2*T*Y)/3))*(((((2*T*Y)/3))/(((T/2)*(((1+(4*Y/T)^2)^(1/2))+(1/(4*Y/T))*LN((4*Y/T)+((1+(4*Y/T)^2)^(1/2)))))))^(2/3))*(S^(1/2)))/N,Y);//Tirant P/S cc Parablicag:=9.81;u:=(4*y2)/T;uu:=(4*y22)/T;A2:=(2*T*y2)/3; //Ar a P/S cc ParablicaA22:=(2*T*y22)/3; //Ar a P/S cc ParablicaP2:=T+((8*(y2^2))/(3*T)); //P rm tro P/S cc ParablicaP22:=(T/2)*(((1+(4*y22/T)^2)^(1/2))+(1/(4*y22/T))*LN((4*y22/T)+((1+(4*y22/T)^2)^(1/2)))); //P rm tro P/S cc ParablicaR2:=A2/P2; //Radio hidrulico P/S cc ParablicaR22:=A22/P22; //Radio hidrulico P/S cc Parablicak2:=(T^2)/(8*y2); //Foco d la Parbola P/S cc Parablicak22:=(T^2)/(8*y22); //Foco d la Parbola P/S cc ParablicaV2:=Q/A2; //V locidad P/S cc ParablicaV22:=Q/A22; //V locidad P/S cc ParablicaE2:=y2+((V2^2)/(2*g)); //En rga P/S cc ParablicaE22:=y22+((V22^2)/(2*g)); //En rga P/S cc ParablicaF2:=V2/((g*A2/T)^(1/2)); //# d Froud P/S cc ParablicaF22:=(V22)/((g*((A22)/T))^(1/2)); //# d Froud P/S cc Parablica// R sulados para S ccion s PARABOLICASIF 0 < u < 1 THENPRINT("RESULTADOS");PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("CASEIF F2 1PRINT(" *END;DEFAULTPRINT(" *END;END;RESULTADOS");");Tirant Normal 'y' = "+ROUND(y22,2)+" m");Ar a Hidrulica 'A' = "+ROUND(A22,2)+" m2");P rm tro 'P' = "+ROUND(P22,2)+" m");Radio Hidrulico 'R' = "+ROUND(R22,2)+" m");Foco d la Parbola 'k' = "+ROUND(k22,3)+" m");V locidad 'V' = "+ROUND(V22,2)+" m/s");En rga Esp cfica 'E' = "+ROUND(E22,3)+" m-Kg/Kg");# d Froud 'F' = "+ROUND(F22,3));THENTipo dFlujo = Subcrtico");THENTipo d Flujo = Sup rcrtico");Tipo d Flujo = Crtico");END;// R sulados para S ccion s PARABOLICAS#################// ECUACION PARA CALCULO DE TIRANTES EN SECCION CIRCULARIF C==3 THEN:=(Q*N)/((D^(8/3))*(S^(1/2)));y33:=0.9965*D*(1.56*(^0.4666)*(1-0.565*((0.3353-)^0.4971))); //Tirant P/S cc Circulary3:=ROUND(0.995*D*(1.56*(^0.4666)*(1-0.565*((0.3353-)^0.4971))),3); //Tirantcc CircularP/Sg:=9.81;A3:=0.5*((D/2)^2)*((2*(/180)*ACOS((1-(y3/(D/2)))))-SIN(((180/)*2*(/180)*ACOS((1-(y3/(D/2))))))); //Area P/Secc CircularP3:=D*((2*(/180)*ACOS((1-(y3/(D/2))))))/2; //Permetro P/Secc CircularR3:=A3/P3; //Radio hidrulico P/Secc CircularT3:=2*((y3*(D-y3))^(1/2)); //Esejo de Agua P/Secc CircularV3:=Q/A3; //Velocidad P/Secc CircularE3:=y3+((V3^2)/(2*g)); //Energa P/Secc CircularF3:=(V3)/((g*((A3)/T3))^(1/2)); //# de Froude P/Secc Circular// Resulados ara Secciones CIRCULARESPRINT("RESULTADOS");PRINT("");PRINT(" * Tirante Normal 'y' = "+ROUND(y33,2)+" m");PRINT(" * Area Hidrulica 'A' = "+ROUND(A3,2)+" m2");PRINT(" * Permetro 'P' = "+ROUND(P3,2)+" m");PRINT(" * Radio Hidrulico 'R' = "+ROUND(R3,2)+" m");PRINT(" * Esejo de Agua 'T' = "+ROUND(T3,2)+" m");PRINT(" * Velocidad 'V' = "+ROUND(V3,2)+" m/s");PRINT(" * Energa Esecfica 'E' = "+ROUND(E3,3)+" m-Kg/Kg");PRINT(" * # de Froude 'F' = "+ROUND(F3,3));CASEIF F3 PRINT("END;DEFAULTPRINT("END;*********RESULTADOS");");Tirante Normal 'y' = "+ROUND(y4,2)+" m");Ancho de Solera 'b' = "+ROUND(b1,2)+" m");Area Hidrulica 'A' = "+ROUND(A4,2)+" m2");Permetro 'P' = "+ROUND(P4,2)+" m");Radio Hidrulico 'R' = "+ROUND(R4,2)+" m");Esejo de Agua 'T' = "+ROUND(T4,2)+" m");Velocidad 'V' = "+ROUND(V4,2)+" m/s");Energa Esecfica 'E' = "+ROUND(E4,3)+" m-Kg/Kg");# de Froude 'F' = "+ROUND(F4,3));1 THEN* Tio de Flujo = Subcrtico");1 THEN* Tio de Flujo = Suercrtico");* Tio de Flujo = Crtico");END;// Resulados ara Secciones EN MAXIMA EFICIENCIA HIDRAULICA###############// ECUACION PARA CALCULO DE TIRANTES EN MINIMA INFILTRACIONIF C==5 THENy5:=SOLVE(Q=(((4*Y*(((1+(Z^2))^(1/2))-Z))*Y+Z*(Y^2))*((((4*Y*(((1+(Z^2))^(1/2))-Z))*Y+Z*(Y^2))/((4*Y*(((1+(Z^2))^(1/2))-Z))+(2*Y*((1+(Z^2))^(1/2)))))^(2/3))*(S^(1/2)))/N,Y);//Tirante P/Secc Tra, Rect y Triang:=9.81;b2:=4*y5*(((1+(Z^2))^(1/2))-Z); //Ancho de Solera P/Secc Tra, Rect y TrianA5:=b2*y5+Z*(y5^2); //Area P/Secc Tra, Rect y TrianP5:=b2+2*y5*((1+(Z^2))^(1/2)); //Permetro P/Secc Tra, Rect y TrianR5:=A5/P5; //Radio hidrulico P/Secc Tra, Rect y TrianT5:=b2+2*Z*y5; //Esejo de Agua P/Secc Tra, Rect y TrianV5:=Q/A5; //Velocidad P/Secc Tra, Rect y TrianE5:=y5+((V5^2)/(2*g)); //Energa P/Secc Tra, Rect y TrianF5:=V5/((g*(A5/T5))^(1/2)); //# de Froude P/Secc Tra, Rect y Trian// Resulados ara Secciones EN MINIMA INFILTRACIONPRINT("RESULTADOS");PRINT("");PRINT(" * Tirante Normal 'y' = "+ROUND(y5,2)+" m");PRINT(" * Ancho de Solera 'b' = "+ROUND(b2,2)+" m");PRINT(" * Area Hidrulica 'A' = "+ROUND(A5,2)+" m2");PRINT(" * Permetro 'P' = "+ROUND(P5,2)+" m");PRINT(" * Radio Hidrulico 'R' = "+ROUND(R5,2)+" m");PRINT(" * Esejo de Agua 'T' = "+ROUND(T5,2)+" m");PRINT(" * Velocidad 'V' = "+ROUND(V5,2)+" m/s");PRINT(" * Energa Esecfica 'E' = "+ROUND(E5,3)+" m-Kg/Kg");PRINT(" * # de Froude 'F' = "+ROUND(F5,3));CASEIF F5 1 THENPRINT(" * Tio de Flujo = Suercrtico");END;DEFAULTPRINT(" * Tio de Flujo = Crtico");END;END;// %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%IF u2 > 1 THENPRINT("RESULTADOS");PRINT("");PRINT(" * Tirante Crtico 'Yc' = "+ROUND(yc2,2)+" m");PRINT(" * Area Crtica 'Ac' = "+ROUND(Ac2,2)+" m2");PRINT(" * Permetro Crtico 'Pc' = "+ROUND(Pc22,2)+" m");PRINT(" * Radio Hidrulico Crtico 'Rc' = "+ROUND(Rc22,2)+" m");CASEIF N = 0 THENPRINT(" * Pendiente Crtica 'Sc' = 'n' no esecificado");END;DEFAULTPRINT(" * Pendiente Crtica 'Sc' = "+ROUND(Sc22,6)+" m/m");END;PRINT(" *PRINT(" *PRINT(" *PRINT(" *CASEIF Fc2 < 1PRINT(" *END;IF Fc2 > 1PRINT(" *END;DEFAULTPRINT(" *END;END;Foco de la Parbola 'k' = "+ROUND(kc2,3)+" m");Velocidad Crtica 'Vc' = "+ROUND(Vc2,2)+" m/s");Energa Esec Mnima 'Em' = "+ROUND(Em2,3)+" m-Kg/Kg");# de Froude 'F' = "+ROUND(Fc2,3));THENTio de Flujo = Subcrtico");THENTio de Flujo = Suercrtico");Tio de Flujo = Crtico");END;// Resulados ara Secciones PARABOLICAS###############// ECUACION PARA CALCULO DE TIRANTES CRITICOS EN SECCION CIRCULARIF F==3 THENg:=9.81;yc3:=D*((1+13.6*((Q^2/(g*(D^5)))^-2.1135)-13*((Q^2/(g*(D^5)))^-2.1))^-0.1156);//Tirante Crtica P/Secc CircularAc3:=0.5*((D/2)^2)*((2*(/180)*ACOS((1-(yc3/(D/2)))))-SIN(((180/)*2*(/180)*ACOS((1(yc3/(D/2))))))); //Area Crtica P/Secc CircularPc3:=D*((2*(/180)*ACOS((1-(yc3/(D/2))))))/2; //Permetro CircularRc3:=Ac3/Pc3; //Radio hidrulico Crtico P/Secc CircularTc3:=2*((yc3*(D-yc3))^(1/2)); //Esejo de Agua Crtico P/Secc CircularSc3:=(((Q*N)/(Ac3*(Rc3^(2/3))))^2); //Pendiente Crtica P/Secc CircularVc3:=Q/Ac3; //Velocidad Crtica P/Secc CircularEm3:=yc3+((Vc3^2)/(2*g)); //Energa Mnima P/Secc CircularFc3:=Vc3/((g*(Ac3/Tc3))^(1/2)); //# de Froude P/Secc Circular// Resulados ara Secciones CIRCULARPRINT("RESULTADOS");PRINT("");PRINT(" * Tirante Crtico 'Yc' = "+ROUND(yc3,2)+" m");PRINT(" * Area Crtica 'Ac' = "+ROUND(Ac3,2)+" m2");PRINT(" * Permetro Crtico 'Pc' = "+ROUND(Pc3,2)+" m");PRINT(" * Radio Hidrulico Crtico 'Rc' = "+ROUND(Rc3,2)+" m");PRINT(" * Esejo de Agua Crtico 'Tc' = "+ROUND(Tc3,2)+" m");CASEIF N = 0 THENPRINT(" * Pendiente Crtica 'Sc' = 'n' no esecificado");END;DEFAULTPRINT(" * Pendiente Crtica 'Sc' = "+ROUND(Sc3,6)+" m/m");END;PRINT(" * Velocidad Crtica 'Vc' = "+ROUND(Vc3,2)+" m/s");PRINT(" * Energa Esec Mnima 'Em' = "+ROUND(Em3,3)+" m-Kg/Kg");PRINT(" * # de Froude 'F' = "+ROUND(Fc3,3));CASEIF Fc3 < 1 THENPRINT(" * Tio de Flujo = Subcrtico");END;IF Fc3 > 1 THENPRINT(" * Tio de Flujo = Suercrtico");END;DEFAULTPRINT(" * Tio de Flujo = Crtico");END;END;// Resulados ara Secciones CIRCULAR###############END;// _________________________________________________________________________________________________________________________________IF W==3 THENINPUT( {{L, {" Tra, Rect. o Trian."," Parablica"," Circular"} }},"CAUDAL",{" Seccin "},{"Con qu Desea Trabajar? *.*"});IF L==1 THENINPUT({Y,B,Z,N,S},"DATOS P/SECC TRAP, RECT O TRIAN",{"Tirante 'y' ","Solera 'b' ","Talud 'Z' ","Rugosidad 'n' ","Pendiente 'S' "},{"Tirante en metros","Ancho de Solera en metros","Talud de Diseo","Coeficiente deRugosidad","Pendiente en metros/metros"});END;IF L==2 THENINPUT({Y,T,N,S},"DATOS P/SECC PARABOLICA",{"Tirante 'y' ","Esejo 'T' ","Rugosidad 'n' ","Pendiente 'S' "},{"Tirante en metros","Esejo de Agua en metros","Coeficiente de Rugosidad","Pendiente en metros/metros"});END;IF L==3 THENINPUT({Y,D,N,S},"DATOS P/SECC CIRCULAR",{"Tirante 'y' ","Dimetro 'd' ","Rugosidad 'n' ","Pendiente 'S' "},{"Tirante en metros","Dimetro en metros","Coeficiente de Rugosidad","Pendiente en metros/metros"});END;// ECUACION PARA CALCULO DE CAUDALES EN SECCION TRAPEZOIDAL, RECTANGULAR Y TRIANGULARIF L==1 THENQ1:=((B*Y+Z*(Y^2))*(((B*Y+Z*(Y^2))/(B+(2*Y*((1+(Z^2))^(1/2)))))^(2/3))*(S^(1/2)))/N;//Caudal P/Secc Tra, Rect y Triang:=9.81;Aq1:=B*Y+Z*(Y^2); //Area P/Secc Tra, Rect y TrianPq1:=B+2*Y*((1+(Z^2))^(1/2)); //Permetro P/Secc Tra, Rect y TrianRq1:=Aq1/Pq1; //Radio hidrulico P/Secc Tra, Rect y TrianTq1:=B+2*Z*Y; //Esejo de Agua P/Secc Tra, Rect y TrianVq1:=Q1/Aq1; //Velocidad P/Secc Tra, Rect y TrianEq1:=Y+((Vq1^2)/(2*g)); //Energa P/Secc Tra, Rect y TrianFq1:=Vq1/((g*(Aq1/Tq1))^(1/2)); //# de Froude P/Secc Tra, Rect y Trian// Resulados ara Secciones TRAPEZOIDALES, RECTANGULARES O TRIANGULARESPRINT("RESULTADOS");PRINT("");PRINT(" * Caudal 'Q' = "+ROUND(Q1,2)+" m3/s"+" = "+ROUND(Q1*1000,2)+" lts/s");PRINT(" * Area Hidrulica 'A' = "+ROUND(Aq1,2)+" m2");PRINT(" * Permetro 'P' = "+ROUND(Pq1,2)+" m");PRINT(" * Radio Hidrulico 'R' = "+ROUND(Rq1,2)+" m");PRINT(" * Esejo de Agua 'T' = "+ROUND(Tq1,2)+" m");PRINT(" * Velocidad 'V' = "+ROUND(Vq1,2)+" m/s");PRINT(" * Energa Esecfica 'E' = "+ROUND(Eq1,3)+" m-Kg/Kg");PRINT(" * # de Froude 'F' = "+ROUND(Fq1,3));CASEIF Fq1 < 1PRINT(" *END;IF Fq1 > 1PRINT(" *END;DEFAULTPRINT(" *END;THENTio de Flujo = Subcrtico");THENTio de Flujo = Suercrtico");Tio de Flujo = Crtico");END;// Resulados ara Secciones TRAPEZOIDALES, RECTANGULARES O TRIANGULARES###############// ECUACION PARA CALCULO DE CAUDALES EN SECCION PARABOLICAIF L==2 THENuq:=(4*Y)/T;Q2:=(((2*T*Y)/3)*(((2*(T^2)*Y)/(3*(T^2)+8*(Y^2)))^(2/3))*(S^(1/2)))/N;Q22:=((((2*T*Y)/3))*(((((2*T*Y)/3))/(((T/2)*(((1+(4*Y/T)^2)^(1/2))+(1/(4*Y/T))*LN((4*Y/T)+((1+(4*Y/T)^2)^(1/2)))))))^(2/3))*(S^(1/2)))/N;//Caudal P/Secc Parablicag:=9.81;Aq2:=(2*T*Y)/3; //Area P/Secc ParablicaPq2:=T+((8*(Y^2))/(3*T)); //Permetro P/Secc ParablicaPq22:=(T/2)*(((1+(4*Y/T)^2)^(1/2))+(1/(4*Y/T))*LN((4*Y/T)+((1+(4*Y/T)^2)^(1/2)))); //Permetro P/Secc ParablicaRq2:=Aq2/Pq2; //Radio hidrulico P/Secc ParablicaRq22:=Aq2/Pq22; //Radio hidrulico P/Secc Parablicakq2:=(T^2)/(8*Y); //Foco de la Parbola P/Secc ParablicaVq2:=Q2/Aq2; //Velocidad P/Secc ParablicaVq22:=Q22/Aq2; //Velocidad P/Secc ParablicaEq2:=Y+((Vq2^2)/(2*g)); //Energa P/Secc ParablicaEq22:=Y+((Vq22^2)/(2*g)); //Energa P/Secc ParablicaFq2:=Vq2/((g*Aq2/T)^(1/2)); //# de Froude P/Secc ParablicaFq22:=Vq22/((g*Aq2/T)^(1/2)); //# de Froude P/Secc Parablica// Resulados ara Secciones PARABOLICASIF 0 < uq < 1 THENPRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("********CASEIF Fq2 < 1PRINT(" *END;IF Fq2 > 1PRINT(" *END;DEFAULTPRINT(" *END;END;RESULTADOS");");Caudal 'Q' = "+ROUND(Q2,2)+" m3/s"+" = "+ROUND(Q2*1000,2)+" lts/s");Area Hidrulica 'A' = "+ROUND(Aq2,2)+" m2");Permetro 'P' = "+ROUND(Pq2,2)+" m");Radio Hidrulico 'R' = "+ROUND(Rq2,2)+" m");Foco de la Parbola 'k' = "+ROUND(kq2,2)+" m");Velocidad 'V' = "+ROUND(Vq2,2)+" m/s");Energa Esecfica 'E' = "+ROUND(Eq2,3)+" m-Kg/Kg");# de Froude 'F' = "+ROUND(Fq2,3));THENTio de Flujo = Subcrtico");THENTio de Flujo = Suercrtico");Tio de Flujo = Crtico");// %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%IF uq > 1 THENPRINT("PRINT("PRINT(";PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("PRINT("CASEIF Fq22PRINT("END;IF Fq22PRINT("END;DEFAULTPRINT("END;END;RESULTADOS");");* Caudal 'Q' = "+ROUND(Q22,2)+" m3/s"+" = "+ROUND(Q22*1000,2)+" lts/s")*******Area Hidrulica 'A' = "+ROUND(Aq2,2)+" m2");Permetro 'P' = "+ROUND(Pq22,2)+" m");Radio Hidrulico 'R' = "+ROUND(Rq22,2)+" m");Foco de la Parbola 'k' = "+ROUND(kq2,2)+" m");Velocidad 'V' = "+ROUND(Vq22,2)+" m/s");Energa Esecfica 'E' = "+ROUND(Eq22,3)+" m-Kg/Kg");# de Froude 'F' = "+ROUND(Fq22,3));< 1 THEN* Tio de Flujo = Subcrtico");> 1 THEN* Tio de Flujo = Suercrtico");* Tio de Flujo = Crtico");END;// Resulados ara Secciones PARABOLICAS#################// ECUACION PARA CALCULO DE TIRANTES EN SECCION CIRCULARIF L==3 THENQ3:=((((D^2)/8)*(((2*(/180)*ACOS(1-2*(Y/D))))-SIN(((2*ACOS(1-2*(Y/D)))))))*(((((D^2)/8)*(((2*(/180)*ACOS(1-2*(Y/D))))-SIN(((2*ACOS(1-2*(Y/D)))))))/((D/2)*(2*(/180)*ACOS(1-2*(Y/D)))))^(2/3))*(S^(1/2)))/N;//Caudal P/Secc Circularg:=9.81;Aq3:=(((D^2)/8)*(((2*(/180)*ACOS(1-2*(Y/D))))-SIN(((2*ACOS(1-2*(Y/D))))))); //Area P/Secc CircularPq3:=((D/2)*(2*(/180)*ACOS(1-2*(Y/D)))); //Permetro P/Secc CircularRq3:=Aq3/Pq3; //Radio hidrulico P/Secc CircularTq3:=2*((Y*(D-Y))^(1/2)); //Esejo de Agua P/Secc CircularVq3:=Q3/Aq3; //Velocidad P/Secc CircularEq3:=Y+((Vq3^2)/(2*g)); //Energa P/Secc CircularFq3:=(Vq3)/((g*((Aq3)/Tq3))^(1/2)); //# de Froude P/Secc Circular// Resulados ara Secciones CIRCULARESPRINT("RESULTADOS");PRINT("");PRINT(" * Caudal 'Q' = "+ROUND(Q3,2)+" m3/s"+" = "+ROUND(Q3*1000,2)+" lts/s");PRINT(" * Area Hidrulica 'A' = "+ROUND(Aq3,2)+" m2");PRINT(" * Permetro 'P' = "+ROUND(Pq3,2)+" m");PRINT(" * Radio Hidrulico 'R' = "+ROUND(Rq3,2)+" m");PRINT(" * Esejo de Agua 'T' = "+ROUND(Tq3,2)+" m");PRINT(" * Velocidad 'V' = "+ROUND(Vq3,2)+" m/s");PRINT(" * Energa Esecfica 'E' = "+ROUND(Eq3,3)+" m-Kg/Kg");PRINT(" * # de Froude 'F' = "+ROUND(Fq3,3));CASEIF Fq3 < 1PRINT(" *END;IF Fq3 > 1PRINT(" *END;DEFAULTPRINT(" *END;THENTio de Flujo = Subcrtico");THENTio de Flujo = Suercrtico");Tio de Flujo = Crtico");END;// Resulados ara Secciones CIRCULARES###############END;// _________________________________________________________________________________________________________________________________// _____________________________________________________________________________________________________IF W==4 THENINPUT( {{X, {" Tra, Rect. o Trian. 'b, S, n'"," Circular 'relacin y/d'"," Circular 'S, n'"} }},"OTRAS INTERROGANTES",{" Seccin "},{"Con qu Desea Trabajar? *.*"});IF X==1 THENINPUT( {{K, {" Ancho de Solera 'b'"," Pendiente 'S'"," Rugosidad 'n'"} }},"SECC TRAP, RECT O TRIAN",{" Calcular "},{"A calcular :3"});IF K==1 THENINPUT({Q,Y,Z,N,S},"CALCULO DEL ANCHO DE SOLERA",{"Caudal 'Q' ","Tirante 'y' ","Talud 'Z' ","Rugosidad 'n' ","Pendiente 'S' "},{"Caudal en metros cbicos/segundo","Tirante en metros","Talud de Diseo","Coeficiente de Rugosidad","Pendiente en metros/metros"});END;IF K==2 THENINPUT({Q,Y,B,Z,N},"CALCULO DE LA PENDIENTE",{"Caudal 'Q' ","Tirante 'y' ","Solera 'b' ","Talud 'Z' ","Rugosidad 'n' "},{"Caudal en metros cbicos/segundo","Tirante en metros","Ancho de Solera en metros","Talud de Diseo","Coeficiente de Rugosidad"});END;IF K==3 THENINPUT({Q,Y,B,Z,S},"CALCULO DE LA RUGOSIDAD",{"Caudal 'Q' ","Tirante 'y' ","Solera 'b' ","Talud 'Z' ","Pendiente 'S' "},{"Caudal en metros cbicos/segundo","Tirante en metros","Ancho de Solera en metros","Talud de Diseo","Pendiente en metros/metros"});END;END;IF X==2 THENINPUT({Q,I,N,S},"DATOS P/SECC CIRCULAR",{"Caudal 'Q' ","Relacin 'y/d' ","Rugosidad 'n' ","Pendiente 'S' "},{"Caudal en metros cbicos/segundo","Relacin Tirante/Dimetro","Coeficiente de Rugosidad","Pendiente en metros/metros"});END;IF X==3 THENINPUT( {{H, {" Pendiente 'S'"," Rugosidad 'n'"} }},"SECCIONES CIRCULARES",{" Calcular "},{"r :3"});A calculaIF H==1 THENINPUT({Q,Y,D,N},"CALCULO DE LA PENDIENTE",{"Caudal 'Q' ","Tirante 'y' ","Dimetro 'd' ","Rugosidad 'n' "},{"Caudal en metros cbicos/segundo","Tirante en metros","Dimetro en metros","Coeficiente de Rugosidad"});END;IF H==2 THENINPUT({Q,Y,D,S},"CALCULO DE LA RUGOSIDAD",{"Caudal 'Q' ","Tirante 'y' ","Dimetro 'd' ","Pendiente 'S' "},{"Caudal en metros cbicos/segundo","Tirante en metros","Dimetro en metros","Pendiente en metros/metros"});END;END;//*******************************IF X==1 THEN// ECUACION PARA CALCULO DE ANCHO DE SOLERA EN SECCION TRAPEZOIDAL, RECTANGULARY TRIANGULARIF K==1 THENb1:=SOLVE(Q=((B*Y+Z*(Y^2))*(((B*Y+Z*(Y^2))/(B+(2*Y*((1+(Z^2))^(1/2)))))^(2/3))*(S^(1/2)))/N,B); //Ancho de Solera P/Secc Tra, Rect y Triang:=9.81;Ab:=b1*Y+Z*(Y^2); //Area P/Secc Tra, Rect y TrianPb:=b1+2*Y*((1+(Z^2))^(1/2)); //Permetro P/Secc Tra, Rect y TrianRb:=Ab/Pb; //Radio hidrulico P/Secc Tra, Rect y TrianTb:=b1+2*Z*Y; //Esejo de Agua P/Secc Tra, Rect y TrianVb:=Q/Ab; //Velocidad P/Secc Tra, Rect y TrianEb:=Y+((Vb^2)/(2*g)); //Energa P/Secc Tra, Rect y TrianFb:=Vb/((g*(Ab/Tb))^(1/2)); //# de Froude P/Secc Tra, Rect y Trian// Resulados ara Secciones TRAPEZOIDALES, RECTANGULARES O TRIANGULARESPRINT("RESULTADOS");PRINT("");PRINT(" * Ancho de Solera 'b' = "+ROUND(b1,2)+" m");PRINT(" * Area Hidrulica 'A' = "+ROUND(Ab,2)+" m2");PRINT(" * Permetro 'P' = "+ROUND(Pb,2)+" m");PRINT(" * Radio Hidrulico 'R' = "+ROUND(Rb,2)+" m");PRINT(" * Esejo de Agua 'T' = "+ROUND(Tb,2)+" m");PRINT(" * Velocidad 'V' = "+ROUND(Vb,2)+" m/s");PRINT(" * Energa Esecfica 'E' = "+ROUND(Eb,3)+" m-Kg/Kg");PRINT(" * # de Froude 'F' = "+ROUND(Fb,3));CASEIF Fb PRINT("END;DEFAULTPRINT("END;*********RESULTADOS");");Dimetro 'd' = "+ROUND(Dd,2)+" m");Tirante 'y' = "+ROUND(yd,2)+" m");Area Hidrulica 'A' = "+ROUND(Ad,2)+" m2");Permetro 'P' = "+ROUND(Pd,2)+" m");Radio Hidrulico 'R' = "+ROUND(Rd,2)+" m");Esejo de Agua 'T' = "+ROUND(Td,2)+" m");Velocidad 'V' = "+ROUND(Vd,2)+" m/s");Energa Esecfica 'E' = "+ROUND(Ed,3)+" m-Kg/Kg");# de Froude 'F' = "+ROUND(Fd,3));1 THEN* Tio de Flujo = Subcrtico");1 THEN* Tio de Flujo = Suercrtico");* Tio de Flujo = Crtico");END;// Resulados ara Secciones CIRCULARES###############//*******************************//*******************************IF X==3 THEN// ECUACION PARA CALCULO DE PENDIENTE EN SECCION CIRCULARIF H==1 THENg:=9.81;Ac:=(((D^2)/8)*(((2*(/180)*ACOS(1-2*(Y/D))))-SIN(((2*ACOS(1-2*(Y/D))))))); //Area P/Secc CircularPc:=((D/2)*(2*(/180)*ACOS(1-2*(Y/D)))); //Permetro P/Secc CircularRc:=Ac/Pc; //Radio hidrulico P/Secc CircularSc:=((Q*N)/(Ac*(Rc^(2/3))))^2; //Pendiente P/Secc CircularTc:=2*((Y*(D-Y))^(1/2)); //Esejo de Agua P/Secc CircularVc:=Q/Ac; //Velocidad P/Secc CircularEc:=Y+((Vc^2)/(2*g)); //Energa P/Secc CircularFc:=(Vc)/((g*((Ac)/Tc))^(1/2)); //# de Froude P/Secc Circular// Resulados ara Secciones CIRCULARESPRINT("RESULTADOS");PRINT("");PRINT(" * Pendiente 'S' = "+ROUND(Sc,6)+" m/m");PRINT(" * Area Hidrulica 'A' = "+ROUND(Ac,2)+" m2");PRINT(" * Permetro 'P' = "+ROUND(Pc,2)+" m");PRINT(" * Radio Hidrulico 'R' = "+ROUND(Rc,2)+" m");PRINT(" * Esejo de Agua 'T' = "+ROUND(Tc,2)+" m");PRINT(" * Velocidad 'V' = "+ROUND(Vc,2)+" m/s");PRINT(" * Energa Esecfica 'E' = "+ROUND(Ec,3)+" m-Kg/Kg");PRINT(" * # de Froude 'F' = "+ROUND(Fc,3));CASEIF Fc