Программа № 2

PROGRAM PUMP 2

REAL K0,KI,K2,K3,KPQ,KH,KA

BYTE ESC,CL,HOM

DATA ESC,CL,HOM/27,69,72/

DIMENSION AS(100),HS(100),QS(100), ETAS(I00), DSS(100),

*DKSS(100),P04S(100),P01S(100), Q0S(100)

61 TYPE 54.ESC.CL

54 FORMAT (2X,2A)

OPEN (UNIT= 1,NAME='PUMP 2-DAT;l)

READ (5,20) KPQ,KH,KA,P4i,Q!,RO,E,P04X,P04N,Q0X,Q0N

20 FORMAT (F 16.5)

TYPE 54,ESC,CL

K0=0.03

Kl=0.06

PI=3.1415926

B=l.

43 I=1

WRITE (5,44)

WRITE (1,44)

44 FORMAT (//3X,'A',6X,'Q',6X,'H',6X,’ETA',6X,'DS',6X,'DKS',

*6X,’P04’,5X,’P01’,4X,’Q0’)

P04=P04X

Q0=Q0N

IF(P04N-P04X) 21,21,6

21 IF(Q0N-Q0X) 22,22,6

22 DQ=KPQ*(Q0X-Q0N)

DP=KPQ*(P04X-P04N)

GOTO 1

2 Q0=Q0+DQ

IF(Q0-Q0X) 23,23,6

23 P04=P04-DP

IF(P04-P04N) 6,1,1

1 H=P41/(P41+P04)

IF (H-l.) 24,6,6

24 Q=Ql/Q0

IF (Q*H/(I.-H)-0.4) 26,25,25

25 IF (Q0-Q0X) 2,6,6

26 AN=Q+1.

АХ=(1.+(1.-Н)**0.5/Н

A=AN

IF (AX-AN) 2,2,5

5 A=A+KA*(AX-AN)*B

IF (A-AX) 27,27,2

27 IF (A-1.4) 28,29,29

28 K3=0.22

K2=0.08

GOTO 30

29 K3=0.125+Q/(5.6*(A-l.)-!.14)

K2=0.005+Q/(13.48*(A-l.)-3.41)

IF(A-24.) 30,30,31

31 K2=0.02

K3=0.16

30 Y1=(A-I.-Q)*(A-1.+2.*(A-1.)**2.+Q)

Y2=(K2+K3)*((Q+1.)*(A-1.))**2.

Y3=Kl*(Q*A)**2.

Y4=A**2.*((l.+K0)*(A-l.)**2.-(l.+Kl)*Q**2.)

HA=(Y1-Y2-Y3)/Y4

IF((H-HA)**2.-(0.2*H)**2.) 51,51.50

50 B=l.

GOTO 5

51 B=0.1

IF((H-HA)*(H-HA)/H/H-KH*KH) 33,33,5

33 AS(I)=A

HS(I)=HA

QS(I)=Q

Q0S(I)=Q0

ETAS(I)=QS(I)*HS(I)/(1.-HS(I))

P01=P04+P4l

Y5=Q0**2.*8.*RO/(E*PI)**2./P01

Y6=I.+K.0-(Q/(A-1.))**2.*(1.+K1)

DSS(I)=(Y5*Y6)**0.25

DKSS(I)=DSS(I)*(AS(I)*E)**0,5

P04S(I)=P04/10.**6.

WRITE (5.34) AS(I),QS(I),HS(I),ETAS(I),DSS(I),DKSS(I), /// *P04S(I),P01S(I),Q0S(I)

WRITE (1,34) AS(I),QS(I),HS(I),ETAS(I),DSS(I),DKSS(I),

*P04S(I).P01S(I).Q0S.(I)

I=I+1

KA=-KA

A=AX

GOTO 5

6 WRITE (5.35) P41,Q1,RO,E,P04X,P04N,Q0X,Q0N

WRITE(1,35) P41,Q1,RO,E,P04X,P04N,Q0X,Q0N

34 FORMAT (4F7.3,1X,2F8.5,2F8.4,F7.4)

35 FORMAT (6X,F16.5)

D=0.

TYPE

55 FORMAT (1Х, ۬ ВЫВЕСТИ ДАННЫЕ НА ПЕЧАТЬ (1-ДА; 0-НЕТ)? ۬. *) ACCEPT*,D

IF(D.NE.1)GOTO41

CLOSE (UNIT=1,DISPOSE='PRINT')

GOTO 42

41 CLOSE (UNIT= 1,DISPOSE='DELETE')

42 TYPE 60,ESC,CL

60 FORMAT (2Х,2А,'ПРОДОЛЖАТЬ РАБОТУ (1-ДА; 0-HET)?', °)

ACCEPT*,D

IF(D.EQ.1)GOT0 61

STOP

END

В перечень исходных данных для расчета на ЭВМ включаются параметры струйного насоса, в соответствии с заданным режимом работы, и параметры силового насоса, строго соответствующие его характеристике.

Так же, как и в программе №1, результаты расчета на ЭВМ по программе №2 выводятся в виде таблицы, под которой печатаются исходные данные.

Каждому режиму работы силового насоса в пределах его характеристики будет соответствовать определенная конструкция струйного насоса, обеспечивающая заданные выходные параметры.

Приемы применения программы проектирования струйных насосов будут более подробно рассмотрены в последующих примерах раздела 5.