1 / 3

Ch11 장 예제풀이 예제 11.9

Ch11 장 예제풀이 예제 11.9. 환경공학과 20071464 백하비. [ 예제 11.9 ] N 개의 수조가 한 개의 공통 연결점에 연결되어 있는 시스템내의 유동을 구하는 프로그램을 준비하라 . 예제 11.8 의 데이터를 이용해서 프로그램을 검토하라 .

leoma
Download Presentation

Ch11 장 예제풀이 예제 11.9

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Ch11장 예제풀이예제11.9 환경공학과 20071464 백하비

  2. [예제 11.9] N개의 수조가 한 개의 공통 연결점에 연결되어 있는 시스템내의 유동을 구하는 프로그램을 준비하라. 예제 11.8의 데이터를 이용해서 프로그램을 검토하라. 풀이>> 가장 높은 수조의 높이를 최대연결부 수두로, 가장 낮은 수조의 높이를 최소연결부 수두로 한 두 수두에 이등분법을 사용한다. 식 (5. 8. 15)는 의 각각의 시행값을 구하는 데 이용되며 연속방정식은 이분법의 다음 단계을 결정하는 데 이용된다(그림 11.13). 그림 11.13은 다음페이지에 제시하였다.

  3. 10 REM B: EX11.9 EXAMPLE 11.9 N reservoirs with a common junction 20 DEFINT I, N: DEF FNQ (D1, D2, D3, DH) = D1 * DH * LOG(D2 + D3 / DH)'EQ.(5.8.15) 30 READ N, G, KVIS, II: DATA 3,9.806,1.007E-6,15 40 PRINT : PRINT "N,G,KVIS,II="; N; G; KVIS; II 50 ZMA = -1!: ZMI = 1! 60 FOR I = 1 TO N: READ L(I), D(I), EP(I), Z(I) 70 PRINT "L,D,EP,Z="; L(I); D(I); EP(I); Z(I) 80 NEXT I 90 DATA 3000.,1.,.0002,30.,600.,.45,.0009,18.,1000.,.6,.0006,9. 100 FOR I = 1 TO N: C = SQR(G * D(I) / L(I)): El(I) = -.9650001 * D(I) ^ 2 * C 110 E2(I) = EP(I) / (3.7 * D(I)): E3(I) = 1.784 * KVIS / (D(I) * C)' FOR EQ. <5. 8. 15) 120 IF Z(I) > ZMA THEN ZMA = Z(I) ' ELEV. OF HIGHEST RESERVOIR 130 IF Z(I) < ZMI THEN ZMI = Z(I) ' ELEV. OF LOWEST RESERVOIR 140 NEXT I 150 FOR I1 = 1 TO 15: HJUN = .5 * (ZMA + ZMI): S = 0! ' START BISECTION METHOD 160 FOR I = 1 TO N: HF = Z(I) - HJUN: HFS = SQR(ABS(HF)) 170 Q(I) = FNQ(E1(I), E2(I), E3(I), HFS) * SGN(HF) 180 S = S + Q(I)': PRINT "HF,Q,S="; HF; Q(I); S 190 NEXT I 200 IF S > 0! THEN ZMI = HJUN ELSE ZMA = HJUN 210 NEXT I1 220 PRINT "ELEVATION OF JUNCTION IS "; HJUN 230 PRINT "PIPE DISCHARGES ARE (POSITIVE INTO THE JUNCTION)" 240 FOR I = 1 TO N: PRINT "Q("; I; ")="; : PRINT USING " ###.### "; Q(I); 250 NEXT I: PRINT N,G,KVIS,II= 3 9.806001 1.007E-06 15 L,D,EP,Z= 3000 1 .0002 30 L,D,EP,Z= 600 .45 .0009 18 L,D,EP,Z= 1000 .6 .0006 9 ELEVATION OF JUNCTION IS 24.8801 PIPE DISCHARGES ARE (POSITIVE INTO THE JUNCTION) Q( 1 )= 1.198 Q( 2 )= -0.329 Q( 3 )= -0.869 그림 11.13 예제 11.9에 대한 프로그램과 출력

More Related