Файл: Фаддеев И.П. Эрозия влажнопаровых турбин.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

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

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

04 с'и {

05 с* I

06 с'г )

07 Re,п р

08 2

09 V 1 0 -6

12 Asu

г

^ ^"нач 1

14 а

32 В г

составляющие скорости капли в начальной точке расчета.

диаметр капли, мкм.

коэффициент кинематической вязкости пара, шаг расчета.

--начальная и текущие координаты по оси г.

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

число циклов.

- путь по оси и.

з.ззр' 1 у -|хрS 3 _

Командная память

00 Чт 02

01 Выч1 05

02 Зп 20

03 Умн 20

04 Зп 21

 

 

П е р ф о к а р т а 1

 

 

 

 

 

с2

 

05

Чт

01

с и

 

 

 

сг - <

06

Выч1

04

Си Си

 

Cz ~

C'z

07

Зп

22

Си

Са

 

1

<N

08

Умн

22

и

С

и

 

 

 

 

 

V

 

(Сг

Сг)2

09

Зп

23

и

с

и

 

 

 

\

 

 

Пе р ф о к а р т а 2

 

 

П е р ф о к а р т а 3

10 Чт.

03 сг

 

 

 

20

Чт

25

v

 

 

11 Выч1

06 сг с'г

 

 

21

Умн

08

о21

 

 

 

22 Дел

09

o2|/v = Re

 

 

 

 

 

 

12 Зп

24 сг -

сг

 

 

23 Выч11 07 Renp — Re

 

 

24 Уп1

26

если Renp—Re < 0, то 26

 

 

 

 

 

 

13 Умн

24

(cr -

e r f

 

 

25 Ост

00

остановка счета при Re <

14 Сл

21

(cr ~

‘г ? +

(cz -

czf

 

 

^

Renp

 

 

 

 

 

 

15 Сл

23

Y r ~ c r)2+ { с г -

 

26

Чт

 

 

 

 

с

, у +

 

2 2

с и

-

с 'и

 

 

 

- C'u f

=

y2

27 Умн

2 6 ( С и -

с и ) V ~ v

 

 

 

 

 

16 у

 

V V 2 =

 

 

 

 

 

 

 

V

 

 

28 Дел

0 4 { c u

 

 

17 Зп

25 v

 

 

 

 

 

 

- c ’u ) y v / c u

18 у

 

Y v

 

 

 

29 Умн

 

 

 

 

19 Зп

26 Y v

 

 

 

1 2 A s u

( с и - с и ) V ~ v l c u

 

 

 

 

 

194


 

 

П е р ф о к а р т а 4

 

 

П е р ф о к а р т а 7

30

Умн 32

В г ^

и(си- с 'и) У у с = А с и’

60

Зп

13

Чг+1

 

31

Фр

22 формирование знака Дсц

61

Умн

14

rk+iau

 

62

Сл

15

rk+iau “ Ь

bt

 

 

 

 

 

 

32 Сл

04 с

Дс

~ с и

63

Зп

01

Cuk+l

 

 

 

 

uk

uk+i

64

Чт

13

r k+l

 

33 Зп

30

с

 

 

 

 

 

65

Умн

16

^+1«г

 

 

 

 

uk+i

 

 

66

Сл

17

rk+laz +

Ь-

34

Чт

12 Asu

 

 

 

 

67

Зп

02

 

 

35 Дел

30

As../с

= Дт

СЧ+1

 

 

 

 

uk+i

68

Чт

30

сk+i

 

36

Фр

86

формирование знака Дт

 

37

Зп

33

Дт

 

 

69

Зп

04

с„

 

38

Чт

20

сг — с'г

 

 

 

 

 

 

“ft+l

 

39 Умн 26

(сг — с')

У о

 

 

 

 

 

 

 

П е р ф о к а р т а 5

 

40

Умн

33

Дт (сг — с2) У v

 

41

Умн

32

Sj Дт (с2 — с2) У v =

Асг

42

Фр

20 формирование знака

Дcz

43 Сл

05 с,

+

Ас,

=

с

 

 

 

zk

 

г

 

zk + i

 

44 Зп

34 сzk+i

Дт =

As,

 

45 Умн 33 с

 

 

 

 

zk+i

 

 

г

 

46 Сл

40

s

+

Asг =

szfe+i

47 Зп

40

sZk+1

 

 

 

 

48 Чт

33

Дт

 

Дт =

хк+1

 

49 Сл

46 х/г +

 

 

П е р ф о к а р т а 6

50 Зп

46 Т£+1

 

 

 

 

51 Чт

24 сг сг

 

 

 

52 Умн 32 Bj (сг — с')

 

 

53 Умн 33 AtBj

(сг сг) =

Дc'f

54 Сл

06 с

+

Дс

-

с

 

 

 

rk

 

r

 

V i

 

55 Зп

35 с,

 

 

 

 

 

56 Сл

 

*+1

 

 

 

 

06 с

4- сг

 

 

 

 

 

rk

 

rk+i

 

 

57 Умн 88

0,5

{У У с'

\

= с,

 

 

 

\

rk

rk+i J

rcp

58 Умн 33 cr

Дт =

Asr

 

59 Сл

13 tfr у

Asp = r

 

 

13*

 

 

 

 

 

 

 

 

 

П е р ф о к а р

70

Чт

34 с,

 

 

 

 

zk+i

 

71

Зп

 

05 с,

 

 

 

 

zk+i

 

72

Чт

35 с'

 

 

 

 

rk+i

 

73

Зп

06 с,

 

 

 

 

Д+1

 

74

Чт

86 1

1 = NH l

75

Сл

55 Nk +

76

Зп

55 Nk+1

N

77

Выч11 18 Nk+1 -

78

УпН

80 если Nk+i N = 0, то 80

79

БП

00 переключатель цикла

 

 

П e p ф о к а р т а 9

80

Чт

18 N

 

81

Умн

12 N Asu

 

8 2

Сл

19 su у N Asu

83

Зп

19

% s u

 

84 Ост

01

печать пути

85 Чт

04 си

 

86 Ост

01 печать с'а

87 Чт

05

 

 

88 Ост 01 печать с'г

89 Чт 40 s*k+i

195


 

 

 

 

 

П е р ф о к а р т а

10

 

90 Сл

41

$zk +

SZ£+1

95

Чт

00 0

91

Зп

41

«2

 

96

Зп

55

О

 

97

Зп

40

О

92

Ост

01

печать $2

98

БП

00 переход к следующей груп­

93

Чт

13

г

Г

99

 

 

пе циклов

94

Ост

01

печать

Остановка машины с пульта

Программа расчета траекторий движения влаги по выпуклой поверхности рабочих лопаток

Для расчета систем влагоулавливания и для определения на­ чальных условий расчета движения влаги после схода ее с рабо­ чих лопаток необходимо знание траекторий и скоростей движения влаги по профильным поверхностям рабочих лопаток. Для рас­ чета скоростей и траекторий использованы уравнения (II 1.32) и (III.33), характеризующие движение материальной точки с тре­ нием по поверхности вращающейся пластины. При расчете рас­ сматриваемая профильная поверхность рабочей лопатки заме­ няется вписанным многогранником. Движение частицы опреде­ ляется последовательно по граням, причем движение на каждой грани — это указанное движение на плоской вращающейся по­

верхности в координатных осях х, г (см. гл. III).

При составлении программы командные ячейки машины заняты под решение системы уравнений (III.33) относительно г, х, wx,

wr за интервал времени Дт. В оперативную память машины вве­ дены начальные условия и геометрические характеристики впи­ санных плоскостей.

Под действием кориолисовых сил при определенных началь­ ных условиях поступления капли на вращающуюся пластину влага может начать двигаться по пластине в сторону лопаток НА ступени. В этом случае в программе предусмотрено сравнение

величин х (т -f- Дт) и х (т). Если приращение координаты вдоль оси окажется отрицательным, то для дальнейшего расчета в ко­ мандной памяти следует сменить перфокарты 5 и 6 на 5* и 6*, а также изменить начальные условия для расчета обратного дви­ жения.

Переход с одной вписанной плоскости на другую осущест­ вляется автоматически сравнением текущей координаты с уравне­

нием границы вписанной плоскости хтрп = аш + а1пг. При

х = хгрп расчет траектории движения влаги проводится на сосед­ ней вписанной плоскости. В ячейках оперативной памяти следует

заменить коэффициенты аш и а 1п на

а0п+1 и а1п+1 для уравнения

границы хгр п+1 = а0п+1 -ф а 1п+1г. В

программе на печать через

определенное число циклов выведены координаты х и г. В конце расчетов, кроме координат, печатаются безразмерные проекции скоростей wx и wr-

196


 

 

 

 

Оперативная память

 

 

 

А д р е с

 

С о д е р ж а н и е я ч е й к и

 

01

w'x W (а^,

— в

начале

расчета).

 

 

02

wr

(т) (wr0 — в

начале

расчета).

 

 

соО

Д1 <“11

о

— в

начале

расчета).

 

 

04

х (т) (х0

— в

начале

расчета).

 

 

05

Лр

 

— количество вписанных плоскостей в профильную

 

 

 

часть

лопатки.

 

 

 

06

Ар

 

— номер

ячейки,

содержащей

косинус угла

на­

 

 

 

клона

площадки.

 

 

07

Аа,

 

— номер

ячейки,

содержащей

коэффициент

а г

 

 

 

правой границы плоскости при движении влаги

 

 

 

вправо и левой границы при движении влево.

08

Аао

 

— номер

ячейки,

содержащей

коэффициент

а0

 

 

 

правой границы плоскости при движении влаги

 

 

 

вправо и левой границы при движении влево.

09

N

 

— число циклов.

 

 

 

10

k !

 

— коэффициент трения.

 

 

11

Ат

 

— интервал

времени.

 

 

12—21

 

— значения

косинусов углов (3 от первой до по­

22—31

 

следней плоскости.

 

 

 

— значения

коэффициентов а х уравнений правых

 

 

 

граничных линий плоскостей при движении

32—41

 

вправо и левых при движении влево.

 

 

— значения

коэффициентов а0 уравнений правых

 

 

 

граничных линий плоскостей при движении

 

 

 

влаги

вправо и левых при движении влево.

Командная память

 

 

П е р ф о к а р т а 1

00

Чт11

06

cos Р

 

01

Зп

47

cos Р

 

02

Умн

87

2 cos Р

03

Умн

02

2wf

(т) cos р

04

Зп

48

2wr

(т) cos р

05

Чт

04

х(т )

 

06

Зп

53 х (т)

 

07

Умн

47

X (т) COS Р

08

Умн

47

X (т) COS2 Р

09Выч1 48;X (т) COS2P—-2wz (т) cos Р =

=a t

 

 

П е р ф о к а р т а 2

10 Зп

48 а х

11

Чт

10

К

12

Умн

01

k'w x (т)

13

Выч2

48 а -L— k’тх (7) = а2

14

Умн

11

а 2 Дт

15

Сл

01 а 2 Дт+ау ' (т) = ш' (т+Д т)

16

Зп

42

wx (т + Дт)

17

Чт

47

cos Р

18

Умн

01

wx (т) cos Р

19 Умн

87

2w'x (т) COS Р

197


П е р ф о к а р т а 3

П е р ф о к а р т а 4

‘20

Сл

03

2wx (т) cos Р +

т

:21

Зп

48

2wx (т) cos Р +

г

.22

Чт

10

k'

 

 

23

Умн

02

k wr (т)

 

 

:24

Выч2

48

2wx (т) cos (5 +

г — k'w'rX

 

 

 

X (т) =

а 3

 

:25

Умн

11

а 3 Дт

 

 

:26 Сл

02 а 3 Дт +

w,, (т) = wr (т ф-

 

 

 

+

 

 

:27

Зп

43

wr (т -f- Дт)

 

:28

Чт .

01

ш'х (т)

 

 

:29

Умн

11

wx (т) Дт

 

П е р ф о к а р т а 5

-40 Чт

86 1

41 Выч1 03 1 — 7

42 Уп1 79 если (1 — г) <Г 0

43 Чт 04 х ( т + Дт)

44 Выч1 53 т (т + Дт) — х (т) = Дт

45 Уп1 79 если —Дт, то 79

46 Чт2 07 аг

47 Умн 03 агг (т)

48 Зп

50 ахг (т)

49 Чт2 08 а0

30 Сл

04 wx (т) Дт + х (т) = т (т

 

 

 

+ Дт)

 

31

Зп

04 т (т -|- Дт)

 

32

Чт

02 wr (т)

 

33

Умн

11 wr (т) Дт

 

34 Сл

03 wr (т) Дт + г (т) = г ( т

 

 

 

4- Дт)

 

35

Зп

03 7 (7 -I- Дт”)

 

36

Чт

42

w'x ( / -f Дт)

 

37

Зп

01 шх (т -(- Дт)

 

38

Чт

43 wr (т -f- Дт)

 

39

Зп

02

w’r (т + Дт)

 

 

 

П е р ф о к а р т а 6

50 Сл

50 а0 +

сцг (т) = тгр

51

Выч1

04

тгр — т (т) =

Дт

52 Уп1

65 если

—Дт

 

53

Чт

86 1

 

 

54 Сл

46

1 +

N k = Nk+1

55 Зп

46

N k+1

 

—ДУ

56 Выч1

09

N k+1 N =

57 Уп1

00

если

Nb+1 =

N, то 58

58

Чт

03 7 (т)

 

 

59

Ост

01

печать г (т)

 

 

 

П е р ф о к а р т а 7

 

П е f> ф о к а р т а 8

60

Чт

04

т (т)

 

70

Зп

08

033

 

 

 

 

71

Чт

06

012

+

001

 

61

Ост

01

печать т (т)

72

СлФ

90

012

 

73

Зп

06

013

 

 

 

62

Чт

00

0

 

 

 

 

 

 

 

74

Чт

86

1

 

 

 

63

Зп

46

0

 

 

75

Сл

52

1 +

nk = nk+1

64

БП

00

022

 

 

76

Зп

52

n k+l

 

—Дир

65

Чт

07

+

001

77

Выч1

05

nk+i — «з =

66

СлФ

90

022

78

Уп1

00

если

«£+1 =

яр, то

67

Зп

07

023

 

 

79

 

 

 

 

 

 

68

ЧТ

08

032

+

001

Чт

03

г W

 

 

 

69

СлФ

90

032

 

 

 

 

 

 

 

198