40 likes | 236 Views
Ch11 장 예제풀이 예제 11.10. 환경공학과 20071464 백하비. [ 예제 11.10 ] 그림 11.12 와 같이 N 개의 수조가 파이프로 연결되어 있을 때 흡입수조에서 연결점까지의 관로에 펌프가 위치하고 있다 . 펌프에는 체크 밸브 (check valve) 가 부착되어 있다 . 이 시스템의 유동에 대한 프로그램을 준비하라 .
E N D
Ch11장 예제풀이예제11.10 환경공학과 20071464 백하비
[예제 11.10]그림 11.12와 같이 N개의 수조가 파이프로 연결되어 있을 때 흡입수조에서 연결점까지의 관로에 펌프가 위치하고 있다. 펌프에는 체크 밸브(check valve)가 부착되어 있다. 이 시스템의 유동에 대한 프로그램을 준비하라. 풀이>> 펌프가 유량이 0에서 작동하고 있다고 가정하자. 만일 연결부로 들어가는 유동의 합이 0보다 크다면, 이 프로그램은 적당하지 않다. 연결부 전체유동의 부호가 변할 때까지 펌프를 통과하는 유동이 계속 증가된다고 가정한다. 이 계산에서 얻어지는 것은 이등분법을 적용할 때 사용되는 최대유동과 최소유동이 된다. 펌프의 방정식은 다음과 같이 주어진다. 프로그램과 출력은 그림 11.14에 주어진다. 그림11.14는 다음페이지에 제시하였다.
10 REM B:EX1110 EXAMPLE 11.10 PUMP-RESERVOIR SYSTEM. SUBSCRIPT O REFERS 20 REM TO SUCTION RESERVOIR AND PUMP LINE TO JUNCTION 30 DEF FNHJ (DQ) = Z(0) + A0 + DQ * (A1 + DQ * (A2 + DQ * A3)) - F1 * DQ ^ 2 / (L0G(C1 + C2 / DQ ^ .9)) ^ 2 40 DEFINT I-J, N: DEF FNQ (D1, D2, D3, DH) = D1 * DH * LOG(D2 + D3 / DH) 'Eq.(5.8.15) 50 READ N, G, KVIS, II, DQQ: DATA 4,9.806,1E-6,20,8. 60 PRINT : PRINT "N, G, KVIS, II, DQQ="; N; G; KVIS; II; DQQ 70 READ A0, A1, A2, A3: DATA 100.,-.2,-.03,-.007 75 'PUMP HEAD COEFFICIENTS 8O PRINT "A0, A1, A2, A3="; A0; A1; A2; A3 90 FOR I = 0 TO N: READ Z(I), L(I), D(I), EP(I) 100 PRINT "I, Z, L, D, EP="; I; Z(I); L(I); D(I); EP(I) 110 C = SQR(G * D(I) / L(I)): E1(I) = -.9650001 * D(I) ^ 2 * C 120 E2(I) = EP(I) / (3.7 * D(I)): E3(I) = 1.784 * KVIS / (D(I) * C): NEXT I 130 F1 = 1.325 * L(0) / (2! * G * .7854 ^ 2 * D(0) ^ 5): C1 = E2(0) 140 C2 = 5.74 * (.7854 * D(0) * KVIS) ^ .9 150 DATA .0,10000.,4.5,.00006,12.,2000.,2.,.00005,18.,2500.,2.5,.00008 160 DATA 23.,1800.,2.2,.00005,25.,2000.,2.3,.00009 170 HJ = Z(0) + A0: S = 0! 'SEEING IF THE CHECK VALVE WILL OPEN 180 SOSUB 320 190 IF S > 0! THEN PRINT "CHECK VALVE STAYS CLOSED.USE PROGRAM PROB. 11.60": STOP 200 IF S <= 0! THEN QMIN = 0!: QMAX = 0!: ' QQ ASSUMED FLOW THRU PUMP 210 FOR J = L TO II: QQ = J * DQQ: S = QQ: HJ = FNHJ(QQ): GOSUB 320 220 IF S < 0! THEN QMIN = QQ ELSE QMAX = QQ: GOTO 240: ' BRACKETS FLOW THRU PUMP 230 NEXT J ㅡ>계 속..
240 PRINT "QQ, HJ, S, QMIN, QMAX="; QQ; HJ; S; QMIN; QMAX 250 FOR J = L TO II: Q(0) = .5 * (QMAX + QMIN): S = Q(0): HJ = FNHJ(Q(0)) 260 GOSUB 320: IF S < 0! THEN QMIN = Q(0) ELSE QMAX = Q(0) 270 NEXT J 280 PRINT : PRINT "SUM OF Q'S INTO J="; S; " HEAD AT JUNCTION=": HJ 290 PRINT : PRINT "Q(I)="; 300 FOR I = 0 TO N: PRINT USING " ###.### "; Q(I); : NEXT I: PRINT 310 END 320 FOR I1 = 1 TO N: HF = Z(I1) - HJ: HFS = SQR(ABS(HF)) 330 Q(I1) = FNQ(E1(I1), E2(I1), E3(I1), HFS) * SGN(HF): S = S + Q(I1): NEXT I1 340 RETURN N,G,KVIS,11,DQQ- 4 9.806001 .000001 20 8 AO,A1,A2,A3= 100 -.2 -.03 -.007 I,Z,L,D,EP= 0 0 10000 4.5 .00006 I,Z,L,D,EP= 1 12 2000 2 .00005 I,Z,L,D,EP= 2 18 2500 2.5 .00008 I,Z,L,D,EP= 3 23 1800 2.2 .00005 I,Z,L,D,EP= 4 25 2000 2.3 .00009 SUM OF Q'S INTO J=1.198769e-03 HEAD AT JUNCTION= 22.97019 Q(I)= 20.325 -14.644 -15.100 0.B76 B.543 그림 11.14 예제 11.10에 대한 프로그램과 출력