Главная  Журналы 

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 [ 87 ] 88 89 90 91 92 93 94 95 96 97 98 99

11.5.4. Листинг подпрограммы ADAPT

с сссссссссссссссссссссссссссссссссссссссссссссссссссссссс

SUBROUTINE ADAPT

с--EXAMPLE 15 - WATER SEEPAGE UNDER A DAM

С----

$INCLUDE:COMMON

DIMENSION P(NI,NJ),U(NI,NJ),V(NI,NJ) EQUIVALENCE (F(1,1,1),P(1,1)),(F(1,1/2),U(1,1)), 1 (F(l,l,3) ,V(1,1) )

ENTRY GRID

HEADER=WATER SEEPAGE UNDER A DAM

PRINTF=PRINT15

PL0TF=PL0T15

CALL INTA4(NZX, 3 , NCVX(1) ,5,NCVX(2) ,2,NCVX(3),5) CALL DATA5(XZONE(1),3.,XZONE(2),.5,XZONE(3),3., 1 POWRX(1),-1.2,POWRX(3),1.2)

CALL INTA3(NZY,2,NCVY(1),6,NCVY(2) ,3) CALL DATA2(YZONE(1),3.,YZONE(2) , 1.) CALL ZGRID RETURN

ENTRY BEGIN TITLE(1)= PRESSURE TITLE(2)= и VELOCITY TITLE(3)= V VELOCITY

CALL INTA8(KSOLVE(1),1,KPRINT(1),1,KPLOT(1),1,KPRINT(2), 1 1, KPL0T(2) , 1,KPRINT (3) , 1, KPLOT (3) , 1, LAST, 3)

DC=1.

DO 100 J=1,M1 DO 100 1=1,Ll P(I,J)=0. 100 CONTINUE

DO 110 1=2,L2

IF(X(I).LT.3.) P(I,M1)=100. 110 CONTINUE RETURN

ENTRY OUTPUT

CALCULATE U AND V VELOCITIES AT THEIR STAGGERED LOCATIONS DO 200 J=3,M2 DO 200 1=2,L2

V(I,J)=-DC*(P(I,J)-P(I,J-1))/(Y(J)-Y(J-l)) 200 CONTINUE

DO 210 1=2,L2

V(I,M1)=-FLUXM1(I,1) 210 CONTINUE

DO 220 J=2,M2



DO 220 1=3,L2

U(I, J)=-DC*(P(I,J)-P(I-1,J) )/(X(I)-X(I-l) ) 220 CONTINUE

DO 230 J=2,M2 DO 230 1=2,L2

IF(X(I).GT.3..AND.X(I).LT.3.5.AND.Y(J).GT.3.) THEN U(I,J)=0. и(I+l,J)=0. V(I,J)=0. V(I, J+1)=0. ENDIF 230 CONTINUE

DO 240 IUNIT=IU1,IU2

IFdTER.EQ. 0) WRITE (lUNIT, 250) 250 FORMATdX, 4TER,2X, P(3,5) , 5X, U(3,5) , 5X,

1 V(3,5) ЧбХ, P(5,5) ,5X, U(5,5) ,5X, •V(5,5) )

WRITE (lUNIT,2 60) ITER,P(3,5),U(3,5),V(3,5), 1 P(5,5),U(5,5),V(5,5) 260 FORMAT(IX,12,1P6E11.2) 240 CONTINUE

IFdTER.EQ.LAST) THEN CALL PRINT

COME HERE TO FILL IBLOCK(I,J) BEFORE CALLING PLOT DO 270 J=2,M2 DO 270 1=2,L2

IF(X(I).GT.3..AND.X(I).LT.3.5.AND.Y(J).GT.3.) 1 IBL0CK(I,J)=1

270 CONTINUE CALL PLOT

ENDIF RETURN

ENTRY PHI DO 300 J=2,M2 DO 300 1=2,L2 GAM(I,J)=DC

IF(X(I).GT.3..AND.X(I).LT.3.5.AND.Y(J).GT.3.) 1 GAM(I,J)=0.

300 CONTINUE

COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 310 1=2,L2 KBCJl(I)=2 310 CONTINUE

DO 320 J=2,M2 KBCIl(J)=2 KBCLl(J)=2 320 CONTINUE RETURN END



11.5.5. Результаты расчетов

RESULTS OF CONDUCT FOR CARTESIAN COORDINATE SYSTEM

WATER SEEPAGE UNDER A DAM

V(3,5) O.OOE+00

ITER P(3,5) U(3,5)

0 0.OOE+OO 0.OOE+OO

1 7.33E+01 2.92E+00 -6.14E+00

2 7.33E+01 2.92E+00 -6.15E+00

3 7.33E+01 2.92E+00 -6.15E+00 1=1 2 3

X = O.OOE+OO 3.52E-01 1.04E+00 1

1=8 9 10

X = 3.38E+00 3.72E+00 4.22E+00 4

P(5,5) U(5,5) V(5,5) O.OOE+OO O.OOE+OO O.OOE+OO 6.37E+01 9.50E+00 -4.94E+00 6.37E+01 9.49E+00 -4.95E+00 6.37E+01 9.49E+00 -4.95E+00

4 5 6 7

69E+00 2.28E+00 2.78E+00 3.13E+00

11 12 13 14

81E+00 5.46E+00 6.15E+00 6.50E+00

J=l 23 4 56 7

Y = O.OOE+OO 2.50E-01 7.50E-01 1.25E+00 1.75E+00 2.25E+00 2.75E+00 J = 8 9 10 11

Y = З.ПЕ+00 3.50E+00 3.83E+00 4.OOE+OO

PRESSURE

******

.OOE+00

.OOE+02

,00E+02

,OOE+02

OOE+02

OOE+02

OOE+00

.78E+01

.78E+01

.7SE+01

.71E+01

65E+01

61E+01

,56E+01

.36E+01

.35E+01

.28E+01

.15E+01

96E+01

80E+01

,92E+01

.95E+01

.93E+01

.82E+01

,60E+01

27E+01

92E+01

81E+01

.46E+01

.44E+01

.29E+01

.97E+01

43E+01

5 9E+01

, 41E+01

.96E+01

.94E+01

.75E+01

.37E+01

78E+01

OOE+01

27E+01

.55E+01

.53E+01

.33E+01

.94E+01

,37E+01

72E+01

,20E+01

.25E+01

.22E+01

7.02E+01

.64E+01

12E+01

57E+01

,15E+01

.04E+01

.02E+01

.82E+01

.46E+01

.98E+01

,49E+01

,13E+01

,94E+01

.92E+01

.72E+01

.37E+01

91E+01

46E+01

.12E+01

.OOE+OO

.90E+01

.71E+01

.36E+01

91E+01

45E+01

,12E+01

.OOE+00

.OOE+OO

.OOE+OO

.OOE+OO

,OOE+OO

,OOE+OO

.OOE+OO

.90E+01

.93E+00

.4 6E+00

.85E+00

,4lE+00

,18E+00

о Z ,

, 15E+00

.50E+01

.20E+01

.04E+01

.51E+00

17E+00

50E+00

,4lE+00

.OOE+01

.08E+01

.73E+01

.40E+01

,18E+01

,07E+01

.05E+01

.59E+01

.41E+01

.57E+01

.03E+01

71E+01

56E+01

,54E+01

.73E+01

.OOE+01

.22E+01

,63E+01

,25E+01

,06E+01

.04E+01

,80E+01

.28E+01

.63Ё+01

.06E+01

,67E+01

,47E+01

,45E+01

.85E+01

.43E+01

.88E+01

.36E+01

,98E+01

,78E+01

.75E+01

.87E+01

.51E+01

.02E+01

.54E+01

Л8Е+01

,98E+01

.96E+01

.88E+01

.54E+01

.09E+01

.63E+01

28E+01

,08E+01

,06E+01

.88E+01

.55E+01

.09E+01

.64E+01

,29E+01

,10E+01

.OOE+OO

****** и VELOCITY ******

,OOE+OO

.OOE+00

,OOE+00

,OOE+00

OOE+00

OOE+OO

,OOE+OO

.OOE+OO

.OOE+OO

.29E-01

.91E-01

,03E+00

,40E-01

O.OOE+OO

,OOE+00

.OOE+OO

.79E-01

,07E+00

, 17E+00

, 16E+00

.OOE+00

.OOE+OO

.OOE+OO

.59E+00

.43E+00

,65E+00

,98E+00

,OOE+OO

,OOE+00

.OOE+00

.25E+00

.97E+00

.03E+00

,69E+01

.43E+01

,00E+00

, OOE+00

,73E+00

,86E+00

98E+00

56E+01

,11E+01

,OOE+OO

.OOE+00

.92E+00

.03E+00

,49E+00

,30E+01

.53E+01





0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 [ 87 ] 88 89 90 91 92 93 94 95 96 97 98 99