Название: Компьютерное моделирование теплофизических процессов - Методические указания (А.В. Чичиндаев)

Жанр: Информатика

Просмотров: 1003


3.2.2. ввод исходных данных

В программе расчета теплоемкости выдается следующая информация о используемых исходных данных:

 

 T0,C  TO,C  N  D0,г  P, кПа  IRT

 -50.0  50.0  100  20.    101.3    1

 

где T0, TO – начальное и конечное значения температур, в интервале между которыми рассчитываются значения теплоемкости [˚C]; N – число шагов разбиения всего интервала; D0 – влагосодержание воздуха, [г/кг с.в.]; P – абсолютное давление воздуха, [кПа]; IRT – вспомогательный ключ для задания варианта расчета теплоемкости (в настоящей работе IRT = 1).

3.2.3. Вывод результатов расчета

Вывод результатов расчета производится подпрограммой PR. В процессе ее работы непосредственно на экран выводится справочная информация, сообщающая о проценте заполнения файла результатов R3U.REZ в виде: FILE R3U.REZ: 100 \%, где значения цифр быстро изменяются от нуля до N. Результаты расчета помещаются в специальный файл результатов R3.REZ, который создается автоматически. После окончания работы с программой в этом файле можно будет просмотреть все просчитанные варианты в хронологическом порядке. Выходные данные представлены в удобной для чтения табличной форме, формат которой задается в подпрограмме SUBROUTINE PR.

3.3. Текст программы (файл R3U.FOR)

C       PROGRAM ТЕПЛОEМКОСТЬ 2 R3U.FOR

          COMMON/DAN/T0,TO,N,D0,P,IRT

          DATA T0,TO,N,D0,P,IRT/-50.,50.,100,20.,101.325,1/

          END

          COMMON/DAN/T0,TO,N,D0,P1,IRT

          COMMON/REZ/B(7,100)

          OPEN (50,FILE='R3U.REZ')

  11    CALL UID(IO)

          IF(IO.EQ.0) GOTO 10

          TSH=(TO-T0)/N

          CALL CX1(TSH)

          CALL PR

          GOTO 11

  10    CONTINUE

           STOP

          END

          SUBROUTINE CX1(TSH)

          COMMON/DAN/T0,TO,N,D0,P1,IRT

          COMMON/REZ/B(7,100)

          T1=T0

          P2=P1

          PRINT 55

  55    FORMAT(1X,' GO R3U : ')

          DO 1 I=1,N

          T2=T1+TSH

          CPV=CPB(T1,P1,D0,IRT)

          CALL FCPV(T1,P1,T2,P2,D0,IRT,CPX)

          Q=CPV*(T2-T1)

          B(1,I)=I*1.

          B(2,I)=T1

          B(3,I)=CPV

          B(4,I)=CPX

          B(5,I)=Q

          T1=T2

          PRINT 56,I

  56    FORMAT(1H+,' GO R3U : ',I4,' \%'$)

  1      CONTINUE

          RETURN

          END

          SUBROUTINE PR

          COMMON/DAN/T0,TO,N,D0,P1,IRT

          COMMON/REZ/B(7,100)

  1      CONTINUE

          WRITE (50,2) P1,D0,IRT

  2      FORMAT(

     *   1X,' ТАБЛИЦА ЗНАЧЕНИЙ ТЕПЛО|МКОСТИ : P=',

     *   F5.0,'КПА, D0=',F4.0,' ГР/КГ.'/

     *'************ IRT = ',I3,' ************'/

     *   1X,66(1H-)/

     *   1X,'I NN I T1, I CPV, I CPX(M),I Q, I'/

     *   1X,'I ПП I "C I КДЖ/КГ I КДЖ/КГ I КВТ I'/

     *   1X,'I----I',4(8(1H-),1HI))

          Z0=0.

           PRINT 6

          DO 4 I=1,N

          PRINT 5,I

          WRITE (50,3) (B(I1,I),I1=1,5)

          Z=Z0+I

          IF(Z.EQ.10.) WRITE(50,7)

          IF(Z.EQ.10.) Z0=(-1.)*I

  3      FORMAT(1X,1HI,F4.0,1HI,4(F7.3,' I'))

  4      CONTINUE

  5      FORMAT(1H+,'FILE R3U.REZ :',I4,' \%'$)

  6      FORMAT(1H ,'FILE R3U.REZ :')

  7      FORMAT(' I----I',4(8(1H-),1HI))

          RETURN

          END

3.4. Варианты заданий

Номер

варианта

ТО,

˚С

Т0,

˚С

N

DO,

г/кг

P,

кПа

 

1

2

3

4

5

6

7

8

9

10

–50

–50

–50

–50

–50

–50

–50

–50

–50

–50

50

50

50

50

50

50

50

50

50

50

100

100

100

100

100

100

100

100

100

100

20

20

20

20

20

20

20

20

20

20

25

50

75

100

125

150

175

200

225

250

 

11

12

13

14

15

16

17

18

19

20

–50

–50

–50

–50

–50

–50

–50

–50

–50

–50

50

50

50

50

50

50

50

50

50

50

100

100

100

100

100

100

100

100

100

100

20

20

20

20

20

20

20

20

20

20

275

300

325

350

375

400

425

450

475

500