Расчет воздухоразделительной установки средней производительности КжКАж-0,25
Курсовой проект - Разное
Другие курсовые по предмету Разное
(DZITA2*H*RO2*W2**2.)/(TZ*10**3.) IF (NV.EQ.5) DP2= F2*W2**2.*F1*RO2/ (2.*FPY*10**3.) AA= EL-12. IF (AA.GT.0) MPS= 1.1 IF (AA.LT.0) MPS= .00866*EL**2.-.1952*EL+2.3234 DP1= (DZITA1*MPS*RO1*LO*W1**2.)/ (2*D1*10**3.) WRITE (20, 40) ALFA1, ALFA2, KQ, F1, F2 40 FORMAT (5x, ALFA1=, F8.2, 2x, ALFA2=, F8.2, 2x, KQ=, F8.2, 2x, F1=, F7.4, 2x, DZAR2=, F10.2) WRITE (20, 43) DP1, DP2, H, LO, NT 43 FORMAT (5x, DP1=, F10.3, 2x, DP2=, F10.3, 2x, H=, F7.4, 2x, LO=, F7.4, 2x, NT=, I4) WRITE (20, 55) 55 FORMAT (15x, TABLICA NAVIVKI) WRITE (20, 56) WRITE (20, 57) WRITE (20, 58) WRITE (20, 59) 56 FORMAT (5x, ___________________________________________________________________________,/) 57 FORMAT (5x, NOMER DIAMETR KOLIC KOLIC DLINA) 58 FORMAT (5x, SLOIA SLOIA VITKOV ZAHODOV TRUBKI) NP=0 DO 70 I=1, NKK D(I)= DC+TR*(2*I-1) HH=2*NT*D(I)/ (DHAP+DC)/ NKK NS= HH+.5 IF (I.EQ. NKK) NS= NT-NP NP= NP+NS JZN (I)= NS BS (I)= H/ (TZ*JZN(I)) SL (I)= 3.14*BS (I)*D(I) WRITE (20, 60) I, D (I), BS (I), JZN (I), SL (I) 60 FORMAT (5x, I3, 2x, F7.3, 2x, F7.1, 3x, I3, 5x, F6.1) 70 CONTINUE RETURN END SUBROUTINE ALFAT (W1, D1, DC, MU1, RO1, CP1, LAMDA1, ALFA1, DZITA1, RE1) REAL MU1, LAMDA1, NU1 PR1= CP1*MU1/ LAMDA1 RE1= W1*D1*RO1/ MU1 RCR= DC/ 2. REKR= 2300. +10500.*(D1/RCR) **.3 IF (RE1.LE. 2300) GO TO 20 IF (RE1.GT. 2300. AND. RE1.LE. REKR) GO TO 21 IF (RE1.GT. REKR) GO TO 22 22 NU1= .23*(1. +1.77*D1/RCR)*RE1**.8*PR1**.47) DZITA1= .3164/ RE1**.25 GO TO 23 21 NU1= (RE1/300.)*PR1**.37*(1. +1.77*D1/ RCR) DZITA1= .3164/ RE1**.25 GO TO 23 20 NU1= 3.66 DZITA1= 64./ RE1 23 ALFA1= NU1*LAMDA1/ D1 WRITE (20, 72) D1, W1, RE1, NU1, PR1 FORMAT (5x, D1=, F6.3, 2x, W1=, F6.3, 2x, RE1=, F8.1, 2x, NU1=, F6.2, 5x, PR1=, F6.3) RETURN END SUBROUTINE ALFAM (NV, W2, TR, TZ, D2, MU2, LAMDA2, RO2, ALFA2, DZITA2,RE2) REAL MU2, LAMDA2, NU2 RE2= W2*D2*RO2/ M2 TRO= TR/ D2 TZO=TZ/ D2 IF (NV.EQ.6.AND.RE2.LE.400.OR.RE2.GT.5000) WRITE (6, 257) IF (NV.EQ.6) GO TO 300 IF (NV.EQ.1) GO TO 30 IF (NV.EQ.2) GO TO 31 IF (NV.EQ.3) GO TO 41 30 A= .185 AM= .95 B= 8.1 BM= .2 IF (RE2.LE. 2000. OR. RE2.GT. 10000) WRITE (6, 257) GO TO 300 31 A= .083 AM= .85 B= 13. BM= .21 IF (RE2.LE. 1000. OR. RE2. GT.7000) WRITE (6, 257) NTR= 100.*TRO NTZ= 100.*TZO IF (NTR.EQ. 110. AND. NTZ. EQ.120) GO TO 431 IF (NTR.EQ. 120. AND. NTZ. EQ.120) GO TO 432 IF (NTR.EQ. 115. AND. NTZ. EQ.130) GO TO 433 431 IF (RE2.LE. 1000. OR. RE2. GT.8000) WRITE (6, 257) B= 33.8 GO TO 300 433 IF (RE2.LE. 1500. OR. RE2. GT.4000) WRITE (6, 257) B= 15.9 GO TO 300 432 IF (RE2.LE. 1000. OR. RE2. GT.26000) WRITE (6, 257) B= 13. GO TO 300 41 A= .1 AM= .88 B= 19.2 BM= .1 IF (TZO.LE. 1.8. AND. TZO. GT.1.6) GO TO 334 IF (TZO.LE. 1.6. AND. TZO. GT.1.4) GO TO 333 IF (TZO.LE. 1.8. AND. TZO. GT.1.6) GO TO 334 IF (TZO.LE. 1.2) GO TO 331 334 IF (RE2.LE. 1000. OR. RE2. GT.7000) WRITE (6, 257) A= .195 AM= .8 B= 13.7 BM= .1 GO TO 300 333 IF (RE2.LE. 1000. OR. RE2. GT.7000) WRITE (6, 257) A= .1 AM= .83 B= 17.1 BM= .1 GO TO 300 331 IF (RE2.LE. 800. OR. RE2. GT.44000) WRITE (6, 257) A= .009 AM= 1.1 B= 48.5 BM= .21 GO TO 300 332 IF (RE2.LE.10. OR. RE2. GT.7000) WRITE (6, 257) A= .1 AM= .88 B= 19.2 BM= .1 GO TO 300 300 NU2= A*RE2**AM ALFA2= NU2*LAMDA2/ D2 DZITA2= B*RE2**(-BM) 257 FORMAT (11x, VNIMANIE RE2 NE SOOTVETSTVUET B) WRITE (20, 72) D2, W2, RE2, NU2, A, B 72 FORMAT (5x, D2=, F6.3, 2x, W2=, F6.3, 2x, RE2=, F8.1, 2x, NU2=, F6.2, 2x, A=, F6.3, 2x, B=, F7.4) RETURN END SUBROUTINE ALFAR (W2, DE, MU2, LAMDA2, RO2, CP2, ALFA2, F2,RE2) REAL MU2, LANDA2 RE2=W2*RO2*DE/ MU2 PR2= CP2*MU2/ LAMDA2 ST2= .168*RE2**(-.3)*PR2**(-.666) ALFA2= ST2*W2*CP2*RO2 IF (RE2.GE.30. AND.RE2. LT.100) F2= 12.8/RE**(.64) IF (RE2.GE.100. AND.RE2. LT.5000) F2= 2.65/RE**(.3) WRITE (20, 72) DE, W2, RE2, ST2, PR2, F2 72 FORMAT (5x, DE=, F6.3, 2x, W2=, F6.3, 2x, RE2=, F8.1, 2x, ST2=, F6.2, 2x, PR2=, F6.3, 2x, DZAR2=, F10.8) RETURN END
Приложение 2
Текст подпрограммы конструкторского расчета конденсатора-испарителя из напыленно-оребренных труб RKINOT.
SUBROUTINE RKINOT(trez)/dan/NK1,NK2,Q,p1,d1,d2,st,dTv,HH_,Hok_,cm_,Htr,Hsg_,ntr_/danr/hr,cr,br,dr,Nr/in2/C1,dTsr,GCFsdT
character(80) trez !переменная для сообщений об ошибках
character(10) prodcase(NK1)(1); prod="Aer"(2); prod="N2"(3); prod="Ar"(4); prod="O2"(5); prod="He1"(6); prod="He2"agent(prod) !вычисление свойств криопродукта =Ts(p1) !Tисп в верхнем сечении
R1=R(T1)*1000. !теплота парообразования=rov(p1,T1) !плотность пара=rol(p1,T1) !плотность жидкости
select case(NK2)(1); prod="Aer"(2); prod="N2"(3); prod="Ar"(4); prod="O2"(5); prod="He1"(6); prod="He2"agent(prod)2=T1+dTv/2. !Тконд среднее в верхнем сечении=ps(T2) !рконд в верхнем сечении=R(T2)*1000. !теплота парообразования=rov(p2,T2) !плотность пара=rol(p2,T2) !плотность жидкости=fmu(rol2,T2) !динамическая вязкость=fla(rol2,T2)/1000. !теплопроводность=7000. !можно задать в исходных данных или оценить как в RKIMT
wp1=q01*d1*Htr/(R1*rov1*(995.*st**2.-900.*d1**2.)) !привед. скорость (3)=(SQRT(1.+5.4*wp1*rol1)-1.)/(2.7*wp1*rol1) !относительный уровень (2) =T1*Htr*Ho*rol1/(.102*R1*rov1)!температурная депрессия (1)=dTv-0.5*dlT !средний температурный напор (4)
!комплекс E в формуле (20) Z=E*dTконд
E=4.*fla2*fmu2/(R2*9.8*rol2**2)*(89*.001/(fla2*hr**2))**4*Htr=.6171*9.8*R2*rol2**2*hr/fmu2*(fla2*hr**2/(89.*.001))**3*E**(7./8.) !Aр в (21)
!Lambda ребра принято const=89, полутолщина ребра (t) - const=0.001=8866.8 !коэффициент в (16)=3.14*d1*Htr !поверхность кипения 1 трубки
!пов. конденсации -- ребро + межреберное пр-во на ед. высоты трубки
Fr1=dr+br+2.*sqrt(((cr-br)/2.)**2+hr**2)
F2=Nr*Fr1*Htr !поверхность конденсации 1 трубки=Nr*Ar/(C*F1) !первые две дроби в (26)=.5*d1/89.*log(d1/d2)*c*d1*2./(d1+d2) !множитель в (25) =1.5; dTa=.1; eps=.01 !диапазон поиска dTкип и заданная погрешность=bisection(dTi,dTa,sdT,eps) !расчет dTкип половинным делением=dTsr-dT1-GCF*dT1**1.4655 !dTконд из баланса температур =C*F1*dT1**1.4655 !Qкип 1 трубки=Nr*Ar*dT2**0.875!Qконд 1 трубки=Q*1000./Q1 !количество трубок=Q1/(F1*dT1) !коэффициент теплоотдачи при кипении=Q2/(F2*dT2) !коэффициент теплоотдачи при конденсации(1,*)(1,*) РЕЗУЛЬТАТЫ КОНСТРУКТОРСКОГО РАСЧЕТА(1,*)(1,*)Температура в верхней части испарителя, К ,T1
!write(1,*)Теплота парообразования там же, ,R1
!write(1,*)Плотность пара там же, кг/куб.м ,rov1
!write(1,*)Плотность жидкой фазы там же, кг/куб.м ,rol1
!write(1,*)(1,*)Температура в верхней части конденсатора, К ,T2(1,*)Давление в конденсаторе, МПа ,p2
!write(1,*)Теплота парообразования там же, ,R2
!write(1,*)Плотность пара там же ,rov2
!write(1,*)Плотность жидкой фазы там же, кг/куб.м ,rol2
!write(1,*)Кинематическая вязкость там же, ,fmu2
!write(1,*)Теплопроводность там же, Вт/(м*К) ,fla2
!write(1,*)
!write(1,*)Эквивалентный диаметр, м ,de
!write(1,*)Плотность теплового потока кипения, Вт/кв.м ,q01(1,*)Тепловой поток кипения одной трубки ,Q1(1,*)Тепловой поток конденсации одной трубки ,Q2(1,*)Поверхность кипения одной трубки ,F1(1,*)Поверхность конденсации одной трубки ,F2
!write(1,*)
!write(1,*)Высота трубки, м ,Htr
!write(1,*)Высота столба светлой жидкости, м ,Hsg(1,*)Скорость пара в межтрубном пространстве, м/с ,wp1(1,*)Относительный уровень межтрубного кипения, м ,Ho
!write(1,*)
!write(1,*)Плотн. теплового потока конденсации, Вт/кв.м ,q02(1,*)Температурная депрессия, К ,dlT(1,*)Средний температурный напор, К ,dTsr
!write(1,*)(1,*)Температурный напор кипения, К ,dT1(1,*)Температурный напор конденсации, К ,dT2(1,*)Коэффициент теплоотдачи кипения, К ,ALFA1(1,*)Коэффициент теплоотдачи конденсации, К ,ALFA2(1,*)Поверхность кипения одной трубки, кв.м ,F1(1,*)Кол