Исследование гидродинамики течения рабочего тела в однокольцевом самодействующем клапане поршневого компрессора тема диссертации и автореферата по ВАК РФ 05.04.13, кандидат технических наук Левшин, Владимир Петрович

• Левшин, Владимир Петрович
• кандидат технических наук
• 1984

Оглавление диссертации кандидат технических наук Левшин, Владимир Петрович

ВЩЦЕНИЕ.

ГЛАВА. I.АНАЛИЗ СУЩЕСТВУЮЩИХ МЕТОДОВ ОПРЕДЕЛЕНИЯ ПРОПУСКНОЙ СПОСОБНОСТИ САМОДЕЙСТВУЮЩЕГО КЛАНАНА ПОШНЕВОГО КОМПРЕССОРА И СИЛОВОГО ВЗАИМОДЕЙСТВИЯ ЕГО ЗАПОРНОГО ОРГАНА С

ГАЗОЕШ ПОТОКОМ.

1.1.Определение пропускной способности самодействующего клапана поршневого компрессора.

1.2.0пределение силового взаимодействия запорного органа самодействующего клапана поршневого компрессора с газовым потоком

1.3.Методы расчета параметров потока в районе запорного органа самодействующего клапана, основанные на моделировании фиаической картины течения в его проточной части

1.4.Шли и постановка задачи.

ГЛАВА. 2.ПРИМЕНЕНИЕ ТЕОРИИ СТРУЙ ДЛЯ ИССЛЕДОВАНИЯ ТЕЧЕНИЯ РАБОЧЕГО ТЕЛА В РАЙОНЕ ЗАПОРНОГО ОРГАНА САМОДЕЙСТВУЮЩЕГО КЛАПАНА ПОРШНЕЮГО КОМПРЕССОРА.

2.1.Расчетная модель.

2.2.Формирование струи несжимаемой жидкости

2.3.Определение поправки на сжимаемость

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

ГЛАВА. 3. ЭКСИЕРШЕНТАЛЪНОЕ ИССЛЕДОВАНИЕ ТЕЧЕНИЯ РАБОЧЕГО ТЕЛА В ПРОТОЧНОЙ ЧАСТИ МОДЕЛИ ОДНОКОЛЫЦВОГО КЛАПАНА.

3.1. Цель эксперимента и описание экспериментального стенда.

3.2.Методика проведения эксперимента и обработки экспериментальных данных

3.3. Результата экспериментальных исследований:

ГЛАВА 4.СОПОСТАВЛЕНИЕ ЭКСПЕРИМЕНТАЛЬНЫХ: ДАННЫХ

С РАСЧЕТНЫМИ.

4.1.Сопоставление эпюр давления .••.•

4.2.Сопоставление расходовых характеристик.

4.3.Результата сопоставления экспериментальных и расчетных данных

ГЛАВА 5.ПРИБШЕННАЯ ГАЗОДИНАМИЧЕСКАЯ ТЕОРИЯ ТЕЧЕНИЯ РАБОЧЕГО ТЕЛА В РАЙОНЕ ЗАПОРНОГО ОРГАНА ОДНОКОЛЬЦЕВОП) КЛАПАНА ПОШНЕВОГО КОМПРЕССОРА.

5Л.Существование двух режимов течения.

5.2. Соотношения, описывающие процесс формирования струи для несжимаемой жидкости.

5.3.Приближенное исследование характера формирования струи в пространстве между седлом и запорным органом •

5.4.Соотношения, описывающие процесс формирования струи в однокольцевом клапане для сжимаемого рабочего тела.

ГДАВА 6.ИНЖЕНЕРНАЯ МЕТОДИКА РАСЧЕТА ПРОПУСКНОЙ

СПОСОБНОСТИ ОДНОКОЛЬЦЕВОГО САМОДЕЙСТВУЮЩЕГО ШПАНА И СИЛЫ ВЗАИМОДЕЙСТВИЯ ЕГО ЗАПОИЮГО

ОРГАНА С ГАЗОМ ПОТОКОМ.

6«I.Инженерная методика расчета пропускной способности однокольцевого самодействующего клапана поршневого компрессора • •

6.2.Некоторые примеры расчета пропускной способности однокольцевого самодействующего клапана поршневого компрессора .„ б. 3.Расчет газовых сил, действующих на запорный орган однокольцевого самодействующего клапана поршневого компрессора.

ВЫВОДЫ.

Введение диссертации (часть автореферата) на тему «Исследование гидродинамики течения рабочего тела в однокольцевом самодействующем клапане поршневого компрессора»

Директивы Партии и Правительства / I т 5 /, направленные на ускорение научно-технического прогресса и повышения эффективности народного хозяйства, предъявляют серьезные требования к повышению технического уровня и качества продукции компрессорострое-ния.

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

Оамечают / 94 /, что "потери энергии в сопротивлениях газового' тракта слагаются из потерь в клапанах и коммуникациях.Потери энергии в клапанах относятся к наиболее значительным в ПК.Они часто больше, чем сумма всех остальных потерь, и достигают 20 --25$ энергии, затрачиваемой на его привод. Такое положение явилось в значительной мере результатом отсутствия рационального и не трудоемкого метода расчета размеров клапанов и возникающих в них потерь давления.Необходимость такого метода трудно переоценить, так как суммарные потери в клапанах компрессоров, действующих в народном хозяйстве СССР, составляют сотни тысяч киловатт".

Наиболее эффективный путь снижения металлоемкости ПК - увеличение средней скорости поршня и частоты вращения вала.Отмечают / 55 /, что при этом "клапаны цилиндров компрессоров работают в тяжелых условиях: при ударных нагрузках с высокой цикличностью -3 г 50 с"1 (200. - 3000 циклов в минуту) , при высоких давлениях - до 350 мПа и повышенных температурах.Вследствие этого срок службы самодействующих клапанов невелик и является определяющим в уровне надежности работы компрессоров. Самодействующие клапаны, надежно работающие при частоте вращении ваяа 3 г 9 с""1, совершенно неработоспособны и неэкономичны при частотах 16 г 50 с""*, к которым стремятся сегодня"»

По данным Союзпромэнерго / 69 / потери энергии в клапанах достигают 70 % общих потерь в газовом тракте, а 43 г 63 % отказов поршневых горизонтальных компрессоров для сжатия аэотоводо-родной смеси составляют отказы самодействующих клапанов.

В настоящее время развитию исследований работы самодействующих клапанов ПК способствуют развитие измерительной техники быстропротекающих процессов и широкое внедрение метода® математического моделирования ПК с использованием ЭМ / 55; 76 /.Эти исследования включают в себя комплекс вопросов: газодинамики течения рабочего тела, динамики движения запорного органа и его взаимодействия с ограничителем подъема, прочности, надежности и др.

Если деформация запорного органа под действием газовых усилий существенно меньше его геометрических размеров и высоты подъема, газодинамическая задача о течении рабочего тела в проточной части клапана и динамическая задача о движении его запор-ног© органа могут решаться независимо.При решении газодинамической задачи запорный орган можно рассматривать как материальную точку.Это справедливо для кольцевых, дисковых, тарельчатых, пятачковых и цельнолитых клапанов.

Если деформация запорного органа определяется действующими на него газовыми усилиями, то газодинамическая задача и задача о динамике движения запорного органа должны рассматриваться как одна общая. Это относится ко всем прямоточным клапанам с самопружинящим запорным органом (типа.ПИК, ПИК-А и др.) •

Особый случай могут составлять самопружинящие полосовые клапаны и полосовые клапаны с гибким ограничителем подъема.Если профиль запорного органа таких клапанов разбить на ряд ступень-чатах участков, то для каждого участка газодинамическая задача также может решаться независимо от динамической (механической) • Здесь возникают две трудности: а) сшивка участков между собой; б)нахождение профиля запорного органа.

В данной работе для приближенного моделирования физической картины течения рабочего тела в проточной части однокольцевого клапана используется аппарат теории струй, а запорный орган рассматривается как материальная точка. Целью работы является создание приближенного инженерного метода расчета пропускной способности однокольцевого самодействующего клапана Ж и определения силового взаимодействия его запорного органа с газовым потоком. Автор защищает: приближенную газодинамическую теорию течения рабочего тела в одноколытевом самодействующем клапане поршневого компрессора;инженерную методику расчета пропускной способности однокольцевого самодействующего клапана и силового взаимодействия его запорного органа с газовым потоком.

Заключение диссертации по теме «Гидравлические машины и гидропневмоагрегаты», Левшин, Владимир Петрович

ВЫВОДЫ

Г.Проведена теоретическое и экспериментальное исследование гидродинамики течения рабочего тела в однокольцевом самодействующем. клапане поршневого компрессора.

2.Предложена приближенная физическая картина течения рабочего тела в районе запорного органа клапана, позволяющая использовать. аппарат теории струй сжимаемой жидкости для расчета параметров потока.

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

4.В первом режиме расчетная модель удовлетворительно согласуется с экспериментальными данными при постулировании отрыва потока с острой кромки седла.В пространство., расположенное за запорным органом, имеет место ажекция сформировавшейся струи. Реальное "эффективное проходное сечение щели клапана" в этом случае существенно меньше геометрического.

5.Во втором режиме необходимо ввести поправочную функцию для определения асимптотического значения ширины струи.Значенш этой функции найдены сопоставлением экспериментальных эпюр давления на поверхности запорного органа с теоретическими, получаемыми методами теории струй.Давление на срезе запорного органа клапана в этом режиме практически равно давлению в окружающем пространстве, где происходит затопление сформировавшейся струи.Реальное "эффективное проходное сечение щели клапана" остается меньшим геометрического, но используется "рациональнее", чем в первом режиме.

6.Найдено, что критерием, определяющим режим течения, может служить величина "относительной длины канала" £/ .В узком имеет место переходной режим, где течение неустойчиво.

7-В результате проведенных исследований разработана приближенная газодинамическая теория течения рабочего тела в сднакльдевам самодействующем клапане поршневого компрессора, позволяющая: а) исследовать характер формирования струи в районе запорного органа; б) разработать инженерные методики расчета пропускной способности клапана и силового взаимодействия его запорного органа о рабочим телом; в) оценить влияние геометрических параметров клапана и режима работы на его пропускную) способность и силовое взаимодействие запорного органа с рабочим телом.

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

9.Предложены универсальные безразмерные продувочные характеристики клапана для определения его пропускной способности.Эти характеристики могут быть, использованы в различных прикладных задачах: определение абсолютного массового расхода рабочего тела через клапан; расчет газообмена в ступени поршневого компрессора с самодействующими клапанами; оптимизация энергетических показателей ступени поршневого компрессора и т.д.

10 .Полученные результаты: позволяют в первом приближении объяснить наличие минимума у зависимости коэффициента давления от относительной высоты подъема запорного органа Ал .С увеличением величины ^/¿д уменьшается "относительная длина канала" и происходит смена режима течения со второго на первый , в результате чего давление на кромках запорного органа становится меньше, чем давление в окружающем пространстве.Очевидно, газовая сила, действующая на запорный орган и отталкивающая его от поверхности седла при малых высотах подъема с ростом величины: ^/2может изменить направление.

Список литературы диссертационного исследования кандидат технических наук Левшин, Владимир Петрович, 1984 год

1.Материалы ХХУ1 съезда КПСС. - М.: Политиздат, 198Г. - 223 с. 2.Черненко К. У. Наше общее дело и дело каждого коммуниста. -- Коммунист, 1983, № 15, с. 14-26.

2. З.Черненко К.У. Речь на апрельском (1984 г) Пленуме ЦК КПСС.- Правда, II апреля 1984 г.

3. Марчук Г.И» Научно-технический прогресс основа интенсификации общественного, производства. - Коммунист, 1983, $ 4, с. 61-72.

4. Щербаков А.И. Эффективность научной деятельности в СССР.- М.: Экономика, 1982- 223 с.

5. Абрамович Г.Н. Прикладная газовая динамика. М.: Наука, 1976. - 888 с.

6. Абрамович Г.Н. Турбулентное смешение струй. -М.: Наука, 1974. 272 с.

7. Афанасьева И. А. Исследование всасывающих, клапанов для низкотемпературных холодильных поршневых компрессоров.: Автореф. дис. . канд. техн. наук. Л., 1979*- 16 с.

8. Баранов В.Л., Барышников Г.А., Левшин В.П. Расчет процесса газообмена в быстроходном поршневом компрессоре холодильной установки. Тзз. докл. II Всесоюз. науч.-техн. конф. по холодильному машиностроению!. - М.: ЦИНЖШЛНЕФ'ШЙШ, 1978, с. 29—30.

9. Барышников Г.А., Левшин В.П. Математическое моделирование газодинамических: процессов у запорного органа клапана поршневого компрессора. Известия Вузов.Машиностроение, 1982, В II, с. 86-90.

10. Барышников Г.А., Левшин В.П. Учет сжимаемости рабочего тела при моделировании газодинамических процессов в районе запорного органа клапана.поршневого компрессора. Известия Вузов.Машиностроение, 1983, № I, с. 55-59.

11. Барышников Г.А., Левшин В.П. Расчет газовых сил, действующих на запорный орган самодействующего клапана поршневого компрессора. Известия Вузов.Машиностроение, 1984, № 2, с. 99-103.

12. Барышников Г.А., Левшин В.П. » Исаков В.П., Мясников В.Г.

13. Инженерная методика расчета пропускной способности щели колырвого самодействующего клапана поршневого компрессора.-- Известия Вузов.Машиностроение„ 1984, № 4, с.

14. Барышников Г.А., Левшин В.П. Особенности расчета формирования струй в многоэлементном самодействующем клапане поршневого компрессора. Известия Вузов.Машиностроение, 1984,5, с.

15. Белоногов В.И. Исследование работы прямоточного клапана в высокооборотном поршневом компрессоре.: Диссерт. .канд. техн.наук. Л.: ЛПИ им. М.И.Калинина, 19-75. - 162 с.

16. Бергман С. Интегральные операторы в теории линейных уравнений с частными производными. М.: Мир, 1964. - 324 с.

17. Биркгоф Г., Сарантанелло Э. Струи, следы и каверны. -М.: Мир, 1964. 466 с.

18. Борисоглебский А.Н., Кузьмин Р.В. Судовые компрессорные машины и установки. Л.: Судостроение, 1972. - 356 с.

19. Браиловская И.Ю. Расчет обтекания угла потоком вязкого сжимаемого газа. Известия АН СССР, Механика жидкости и газа, 1967, В 3, с. 21-26.

20. Бунимович А.И. Об истечении газа с большими дозвуковыми скоростями. Ученые записки МГУ" им. М.В.Ломоносова, Механика, 1951, т.III, вып. 152, с. 50-61.

21. Ван-Дайк М. Методы возмущений в механике жидкостей. М.: Мир, 1967. - 310 с.

22. Владимиров ф.С. Проникание струи в канал. Прикладная математика и техническая физика, 1967, № 3, с. 103-108.

23. Владимиров ф.С. Струйное обтекание пластинки вблизи стенки.-- Прикладная математика и техническая физика, 1968, № 3,с. 132-136.

24. Владимиров Ф.С. Соударение струй, вытекающих из каналов с параллельными стенками. Прикладная механика и математика, 1969, T.XXXIII, вып. I, с. 21-26.

25. Воронина Л.П., Дмитриев И.Н., Исаков В.П. Клапаны самодействующие кольцевые и дисковые поршневых компрессоров. Отраслевой стандарт СССР: ОСТ 26-12-2030-81, t. 60-112.

26. Гинзбург В.М. и др. Оптическая голография. Практические применения. -М.: Советское радио, 1978. 240 с.

27. Годунов С.К.„ Рябенький B.C. Разностные схемы. М.: Наука, 1977. - 440 с.

28. Губарев Г.В., Френкель М.И., Шварц И.Н. К расчету с помощью' ЭШ Диаграмм движения запорных органов самодействующих клапанов поршневых компрессоров. Труды ЛенНИИхиммаш, 1969, вып. 5, с. 127-132.

29. Гуревич М.И. Тзория струй идеальной жидкости. М.: Наука, 1979. - 536 с.

30. Дейч М.Е. Техническая газодинамика.-М.: Энергия, 1974.- 592 с.

31. Дмитревский В.А., Исаков В.П. Исследование газодинамических характеристик многоканального) дискового клапана.-Труды II Всесоюз. науч.-техн. конф. по компрессоростроенито, М.: ЩНШШНБФТЕМДП, 1970, с. 82-86.

32. Доллежаль H.A. К теории самодействующего клапана поршневого компрессора.-Химическое и нефтяное машиностроение, 1939,1. I, с. 1-8.

33. Доллежаль H.A. Опыт определения сопротивления самодействую,-щих пластинчатых клапанов поршневого компрессора.-В сб.: Котлонадзор: Труды ВИ1М, М., 1940, с. 10-14.

34. Доллежаль H.A. Прикладная теория всасывающего клапана поршневого компрессора.-Общее машиностроение, 1941, № I,с.16-22.

35. Доллежаль H.A. Расчет основных параметров самодействующих пластинчатых клапанов поршневого компрессора.-Общее машиностроение, 1941, J6 9, с.2-5.

36. Исаков В.П. Исследование динамики и прочности самодействующих клапанов.: Диссерт. . канд. техн. наук, Л.: ЛПИ им. М.И.Калинина, 1969.- 172 с.

37. Карапетян Р.Н., Прилуцкий И.К., Хрусталев B.C., фотин B.C. Газодинамические характеристики элементов поршневых компрессоров. Промышленность Армении, 1975, № II, о. 42-44.

38. Карапетян Р.Н., Прилуцкий И.К., Хрустадев Б.С. Выбор клапанов о упругим ограничителем подънма. Промышленность Армении, 1976, $ 8, о. 38-40.

39. Кибель И.Л., Кочин Н.Е., РОзе Н.В. Теоретическая гидромеханика, т. II. М.: ^гаматгиз, 1963, - 727 с.

40. Ковано К., <$гтакова А., Сахо К., Даймон К. Динамика кольцевых нагнетательных клапанов холодильного компрессора поршневого типа. Мищгбиси Дэнки Гихо (япон.), 1972, т.46,3, с. 319-323.

41. Ковеля В.М. Численный метод расчета стационарных: уравнений Навье-Стокса сжимаемого газа. В сб.: Численные методы механики сплошной среды, Новосибирск, 1970, т.1, № 3,с.12-18.

42. Кондратьева Т.Ф. Исследование полосовых самопружинящих клапанов.: Диссерт. . канд. техн. наук, Л.: ЛенНИИхиммаш* 1953, 102 с.

43. Кондратьева Т.Ф. Определение потерь в клапанах поршневого компрессора. Труды НИИхиммаша, 1954, № 18, с;. 88-102.

44. Кондратьева Т.Ф. Об определении потерь энергии в самодействующих' клапанах поршневого, компрессора. Труды НИИхиммаша, 1958, $ 22, с. 113-115.

45. Кондратьева Т.Ф., Исаков В.П. Клапаны поршневых: компрессоров. Л.: Машиностроение, 1983. - 160 с.

46. Лаврентьев М.И., Шабат Б.В. Методы теории функций комплексного переменного. М.: Наука, 1973. - 736 с.

47. Лаврентьев М.И., Шабат Б. В. Проблемы гидродинамики и их математические модели. М.: Наука, 1977. - 408 с.

48. Лаврентьев М.М. О некоторых некорректных задачах математической физики. Изв. Сиб. отд. АН СССР, Новосибирск, 1962, с. 61-68.59• Лешцянский Л.Г. Ламинарный пограничный слой. М.: Физмат-гиз, 1962. - 479 с.

49. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука, 1978. - 736 с.

50. Макеев H.H. 0 соударении газовых струй. Прикладная механика и математика, 1962, т. ХХУ1, вып.2, с. II5-II9.

51. Макеев H.H. Эксцентрическое соударение газовых струй.- В сб.: Работы аспирантов Казанского государственного университета: Математика, механика, физика. Казань.; Казанский гос. ун-т, 1964, с. 9- 14.

52. Мартынов А.К. Экспериментальная аэродинамика. М.: Оборон-гиз, 1958. - 348 с.

53. Махортых Г. В. Истечение жидкости или газа, движущегося вдоль плоскости, через отверстие в этой плоскости. В кн.: Промышленная аэродинамика. -М.: Оборонгиз, 1962, вып.23, а. 203-211.

54. Мельниченко Л.Г., Лантух H.A., Константинов Л.И. Расчет характеристик холодильных компрессоров с применением электронных цифровых вычислительных машин. В сб.: Исследование работы судовых холодильных установок, Калиниград, 1970,с. 17-22.

55. Назаров Г.Н. Давление газовой струи на равнобокий клин.- Прикладная математика и техническая физика, 1962, № I, с. 102-106.

56. Назаров Г.Н. Метод Бершана в теории струй идеального газа.- Докл. науч. конф. по математике и механике. ТЪмск: ЗЬмский гос. ун-т, 1962, с. 31-35.

57. Новиков И.И., Сафин А.Х. Технический уровень и направления развития компрессоростроения для химической, нефтехимической и других отраслей промышленности. Компрессорное машиностроение, М.: ЩНТИХШНЕФТЕШШ, 1978. - 58 с.

58. Овсянников Л.В. Групповые свойства уравнений С, А. Чаплыгина. -- Прикладная математика и техническая физика, I960, & 3, с. 73-80.

59. Островский A.M. Решение уравнений и систем уравнений. -М.: Иностранная литература, 1963. 218 с.

60. Отчет ЛенНИИхиммаш по теме № ~ : Исследование конст7327-8рукшй самодействующих клапанов для компрессоров без смазки цилиндров, этап 7: Газодинамические исследования кольцевых и прямоточных самодействующих клапанов. Л., 1968, инв. Л 4558/1. - 54 с.

61. Пирумов И.Б. Исследование динамики и прочности полосовых самодействутощих клапанов поршневых компрессоров.: Автореф. дис. . канд. техн. наук, Л.: ЛПИ игл. М.И.Калинина, 1965.- 16 с.

62. Пирумов И.Б., Фотин Б.С. и др. Разработка методики расчета клапанов поршневого компрессора по интегральным характеристикам. : Тез, докл. У Бсесоюзн. конф. по компрессоростроению, М., 1978, с. 22.

63. Пластинин П.И., Твалчрелидве А.К. Введение в математическое моделирование поршневых компрессоров. М.: МВТУ им. Н.Э.Баумана, 1976, - с.

64. Пластинин П.И. Расчет и исследование поршневых компрессоров с использованием ЭШ. Насосос троение и компрессорос троение. Холодильное машиностроение, т.2.: Итоги науки и техники. - M.s ВИНИТИ АН СССР, 1981. - 168 с.

65. Поска А.А.Исследование новых конструкций прямоточных и кольцевых клапанов и разработка методов их расчета.: Диссерт. . канд. техн. наук, Укмерге: Завод "Бенибе" (также Л.: ЛПИ им. М.И. Калинина) , 1981. IS8 с.

66. Прохур И.З. 0 расчете дозвуковых течений вязкого газа в приближении обобщенных уравнений Прандтля. В сб.: Динамика сплошной среды, Новосибирск: Институт гидромеханики Сиб. отд. АН СССР, 1977, вып. 30, с. II2-I2I.

67. Рождественский Б.Л., Яненко H.H. Системы квазилинейных уравнений и их приложения в газовой динамике.- М. : Наука, 1968.- 592 с.

68. Седов Л.И. Плоские задачи гидродинамики и аэродинамики.- М.: Наука, 1966. 448 с.

69. СЛезкин H.A. Обтекание плоским прерывным потоком криволинейного препятствия. Докл. АН СССР, 1935, т.П, № 8, с.1230- 1240.

70. Справочник машиностроителя, т.2. -М.: Машиностроение, 1955. 333 с.

71. Зб.Стулов В.П. Пограничные слои в диссипативных средах.- Научные труды НИИ механики МГУ" игл. М.В. Ломоносова, 1977, й 48. 79 о.

72. Трошин В.И. Обтекание пластины струей газа, вытекающей из канала. Прикладная механика и математика, 1959, т.ХХШ, вып. 4, с. 14-19.

73. Трошин В.И. Истечение газа, движущегося в канале, через отверстие в стенке. Прикладная механика и математика, 1959, т. XXIII, вып. 5, с. II8-I23.

74. Трошин В.И. Две задачи о дозвуковых газовых струях. Вестник МГУ игл. М.В. Ломоносова, I960» № 2, с. 19-23.

75. ЭО.Трошин В.И. Удар дозвукового газового потока о пластину, прикрывающую вход в канал с параллельными стенками.- В сб.: Механика машин : Изв. АН СССР , Отд. техн. наук, I960, № 4, с. 18-24.

76. Уиттекер Э.Т., Ватсон Д.Н. Курс современного анализа. : ч.И, Трансцендентные функции. -М.: Физматгиз, 1963.- 430 с.

77. Фалькович C.B. К теории газовых струй. Прикладная механика и математика, 1957, T.XXI, вып.4, с. 3-12.

78. ЭЗ.Фздоренко с.В. Исследование изменения температуры газа в цилиндрах поршневых компрессоров.: Автореф. диссерт. . канд. техн. наук, М.5 МВТУ им. Н.Э.Баумана, 1977. 16 с.

79. Френкель М.И. Поршневые компрессоры. Л.: Машиностроение, 1969, - 744 с.

80. Хрусталев Б.С. Исследование работы группы клапанов порлне-вого компрессора.: Диссерт. . канд. техн. наук, Л.: ЛПИ им. М.И.Калинина, 1974, 172 с.9$.Чаплыгин С.А. 0 газовых струях.: Полное собрание сочинений, Т.Н. Л.: АН СССР, 1933. - 290 с.

81. Чекушин Г.Н. Исследование динамики и прочности пластин кольцевых самодействующих клапанов поршневых компрессоров.: Дис-серт. . канд. техн. наук., JI.j ЛПИ им. М.И. Калинина, 1966» 142 с.

82. Шварц И.Н. Применение ЭВЛ для расчета и оптимизации поршневых компрессоров. -М.: 1ЩЖШНЕШМАШ, 1973. 24 с. 99.1ШШХТИНГ Г. Теория пограничного слоя. -М.: Наука, 1974. -- 712 С;.

83. ЮО.Яненко H.H. Метод дробных шагов решения многомерных задач математической физики. Новосибирск: Сиб. отд. АН СССР, 1967. - 56 с.iol Ве ъyman Оп in iegzeß ojoe га Жогб б/ъеа/г? ft/src/Sor? of co/r?p j^^cftc/'i. ЛГоп&л-еаг1. J Уо-zjfc (¿onafo/?;

84. Внбемаяг? J. Hoc/o&zcLjb/?e/7/7?e-//?oc/e,der, ZA MM,v. л/2, f>. ?39- 7S4.

85. Wee/lame^ fffO, f. J7JfiS4J/D.jss--42o. 4 о С. ß. ^t ze-M^to*?frica, ßz, ¿^u~eпсоюр'геЛ .

86. PuM. UsutT. o/e ве^-га^е, Yir. 2, p. ßr-Yos;p^cz/riesr-f ifypé . —s S?? Y. -■/£>&. Hctnti&oM У. £y/err-StoM û-f tca¿ rr? oc/e of J&zHrt'^tí^

87. Jdcoi С. ¿faz. un ^azettó , >e. ^-i9 36^ tr. 2 o p, H O. Jaco g с. /Уе/Ъ ¿k/r¿ó

88. PoT-/¿/^ /¿'сгъ J Y93*, f.3- /Д

89. H< sJacotf C, cTr? ^ 'ом ¿г? я ^б-е/г/я ^ue^a ¿cu ¿r? f ¿¿¿/¿fe. —rtüxie^- Pœ-^^j p.

90. JS2. faeóie e с/. £tj pitee. ¿Г./f. Ор1ете*а-/еои TÍ:ec/?n(&i{-e4 кг té It -У^виГ1. Wc. sm.-^f.

91. Ю. Mac Ate J7?e ¿ofu-fron of £оиг?с/ягрira вне, J-o^ ¿г ^e/7 e ж £ hoc/op1. G-^cc&ty V/7. — Ptoc. of1. Fh¡¿. ^C, Jl^j^Zjf. 4

92. У/4 NacktC A. g Jfp ¿t'c*/to*4 of- Kt*c/t/?off-feMeft* Т&г^¿z, j-€ot*r 4on*с Ib-efe ОС. Ся/гг^г^гг-е.1. Mié.

93. V3/ ЛахгеъзаЫjf.j&ae&^oàt A ¿S*****/ ca-écwéa-^e'os? of fy/*^ zrafoe^

94. Qaefe ъпс Pc^r/i S4 т2^, «л^p. -/Vi32, Jôcf&et* Л Jr ¿o ¿tupí tr&sn-e/i'é of aa^yb^^óo^ ¿raé'tr-e.1. Ьос/оъ zr/^-e^s^--/33. ¿/jofo£¿/ Л ¿-/ифг offPo%*r ¿t? со/эуЬъе&Зоы

95. Л/екГ ,Socr/A 2J&. £24 fdbrs/of ¿u ttc-t/Dloctntcfiip œ-SJo^кг с//? ¿исто/??гг-/<ъ t&e-cf Va б :of т?А<е Xf? ¿torn gof j&ef 'fe o/7j1. SféZj LT. 3. f><f3S: Mehofez. /Уо^оъ к/, eh £0/7¿re 'z^&sr^1. J)¿7Tez^re/tP Jfcz^rr CcnsCe/rcsr-íS t ^SДДрареъ ^ л/POj

96. В этом приложении приводится программа расчета на ЭШ асимптотических характеристик "У^Г и $ струи несжимаемой жидкости, определяемых системой уравнений (2.18) и (2*19). Алгоритм определения этих величин сводится к минимизации функционала / 71 /

97. Здесь и разность левых и правых частейуравнений (2.18) и (2.19) , соответственно.

98. Программа, составленная на языке Фортран-4, транслирована на ЭШ СМ 3 с операционной системой РАФОС.Результаты вычислений для величины ^ оформлены в виде номограмм и приведены ниже.

99. Задание начального приближения3

100. Определение направления екорейшего спуска1

101. Организация одномерного поиска в направлении скорейшего спускаI1. Определение компонент сопряженных направлений11. Организация одномерного поискав сопряженном1. Нет направлении

102. OGRAM MFMCG Продолжение приложения I

103. ГОЛОВНАЯ ПРОГРАММА ДЛЯ F11CG1. ЗАДАВАЕМЫЕ ПАРАМЕТРЫi

104. ОЗНАЧАЕТ/ЧТО СХОДИМОСТИ НЕ БЫЛО В ПРЕДЕЛАХ ЗАДАННОГО ЗНАЧЕНИЯ 11Н/-1-ОЗНАЧАЕТ ОШИБКУ ПРИ ВЫЧИСЛЕНИИ ГРАДИЕНТА/

105. ОЗНАЧАЕТ/ЧТО МИНИМУМ ПОВИДИМОМУ НЕ СУЩЕСТВУЕТ XI/Х2-ЗНАЧЕНИЯ ПЕРЕМЕННЫХ В МИНИМУМЕ Б

106. COMMON/MCG/X<2)/G(2 )/Н С4)/Н/М/F/EST/EPS/LIМ/1ER/KOUNT/ВА/DA1. EST=8.1. N=21.1=41. TVPE i

107. FORMAT</2Х/' ВВЕДИ ПАРАМЕТРЫ i ВА С F5.3 )') ACCEPT 2/BA F0RMAT<F5.3) TVPE 3

108. FORHATO-2X/-'BA-\> : ACCEPT 2/HA TVPE 4

109. FORMAT <У2Х/'НАЧАЛЬНОЕ ЗНАЧЕНИЕ БЕТА') ACCEPT 2/BET K<2)=BET TVPE 5

110. FORMAT <>'2Х/' НАЧАЛЬНОЕ ЗНАЧЕНИЕ КОРНЯ ИЗ'НО ACCEPT 2/SQH XU>=SQH TVPE ь

111. FORMAT <У2Х/ ' ОШИБКА ВЫЧИСЛЕНИЙ <F5.4)') ACCEPT 7/EPS FORMAT(F5.4> TVPE 8

112. FORMAT <У2Х/' МАКСИМАЛЬНОЕ КОЛИЧЕСТВО ИТЕРАЦИЙ <I5)') ЙССЕРТ Э/LIM FORMAT<15) CALL FMCG1. TVPE 18

113. FORMAT 0-'2Х/ 'КОНЕЦ МИНИМИЗАЦИИ/ЗАДАЙ НОМЕР КАНАЛА ВЫВОДА <11)') ACCEPT 11/НС FORMAT*! II) WRITEСНС/12)

114. FORMAT О' 19Х/ ' КОНЕЦ МИНИМИЗАЦИИ ПО FMCG') I.JR I ТЕ < НС/ 13)1 ER/ LIМ/ EPS/ KOUNT

115. FORMAT С -'ЗХ/' IER=' / 12/ 2Х/ ' LIM=' /15/ 2Х/ ' EPS=' / F7.5/ ЭХ/'КОЛИЧЕСТВО ИТЕРАЦИй='/15) WRITEСНС/14)

116. Продолжение приложения I FORMAT (. ••■'ЗХ,' НАЧАЛЬНЫЕ ЗНАЧЕНИЯ ПЕРЕМЕННЫХ7^ ~~ ::l-JR I ТЕ С НС/ 15 > В A, D А/ BET, SQH

117. FORMAT <.'ЭХ,'В^А=',F5.3,2Х,'ДЕЛЬТА/А=',F5.3,2Х, ЕТА=' ,F5.3, 2Х, КОРЕНЬ ИЗ Н=',F5.3 > URITECHC,16)

118. FORMAT</2Х,'КОНЕЧНЫЕ ЗНАЧЕНИЯ ПЕРЕМЕННЫХ'> Х1=Х<1>1. У ("N1. У RI ТЕ<NC, 17 > Х2, X1

119. FORMAT С /2Х, ' БЕТ А=', F6.4, 2Х,' КОРЕНЬ ИЗ Н=', Fb. 4 >1. STOP1. END1. SUBROUTINE FMCG

120. МЕТОД СОПРЯЖЕННЫХ ГРАДИЕНТОВ ФЛЕТЧЕРА

121. DV=DV+G < I > *H a > " * j^^oJj^'eii pMJiO^e HMH/I1.<DV>17/38/ 281. ( FV-FX > 18.» 28." 281BDA=A!1BDA+ALFA1.A=AHBDA1. < HNRI1* AHBD A-1. E18 > 14/ 14/191.R=21. RETURN1. T=8.1. < AHBD A > 22/ 38/ 22 . T

122. Z=3. * < FX-F V > /AHBDA+DX+DV ¿LF A=AI-1 AX 1 < ABS < Z> / AE'S < DX >/ABS < DY >> DALFA=Z-'ALFA . '

123. DALFA=DALFA*DALFA-DX.'-ALFA*DV.-'ALFA1. < DALFA >23/27/271. DO 24 J=i/N1. H+Jj>=Hao1. CALL FUNCT1.<IER>47/2b/471.R=-11. GOTO 1bJ=ALF ASSORT < D ALF A >1. AIFfi= < BV+2. *!.<f-DX >1. DO 28 1 = 1/N "" Ix <: i >=x a >+(t-aLfa > :*h a >

124. CALL FUNCT . IF < F-FX >23/29/38 IF < F-F V > 38/ 38/ 381. DALFA=8. .

125. DO 31 1=1/N DALFA=DALFA+6 <I>*HCI> IF < D ALFA > 32/ 35/ 35 IF < F-FX >34/33/35 IFCDX-DALFA>34/38/34

126. ВЫЧИСЛЕНИЕ МИНИМИЗИРОВАННОЙ ФУНКЦИИ И ЕЕ ЧАСТНЫХ ПРОИЗВОДНЫХ ДЛЯ FHCG

127. COMMON-'MCG-'X С 2 ),G<2>, Н<4 N, И.- F, EST.» EPS, LI И, IER, KOUNT, В A, DA1. SQH=X<1>1. BET=X<2>1. P1=3.1415331. PI4=PI/41. PIi=l/PI1. Si2P=SQH/PI1. BET2=BET--'21. SIB=SIN<BET>1. COB=COSCBET>

128. TG=S IN С BET2+P14 > У COS < BET2+F' 14) ALTG=AL0G(T6>

129. ARC=ATAN < 2*SQH/ <. 1 -SQH**2 >> FII=2*SIB*ALTG F12= < (1 -SQH**2 > /SQH ) *:ARC FI3=PI*C0B

130. FI =B A- i ■-SQP5K < F11 •-F12-F13 >

131. DFH1=< 2:* < 1 -SQH**2 > > < 1 +SQH**2 У1. DFH2=2*SQH»ARC1. DFH3=2»SIB*ALT61. DFH4=F'I*C0B

132. DFH=PI1ж<DFH1-DFH2-DFH3+DFH4>1. DFB1=2*C0B»ALTG1. DFB2=2^ (SI B/COB )

133. DFB=-SQP* < DFB1+DFB2+PI*SIВ >

134. AL=ALOG С СI+SQH > / < 1 -SQH > > ALC=ALOG С COS ( BET/2 ) /SIN С BET/2 > ) PS 11=С < 1 +SQH**2 > /SQH > ж AL

135. PSI=DA-SQP*< PS11+PI*SIB-2*C0B*ALC>1. DPS1=2*SQH»AL

136. ВР52=2ж С1 +SQH**2 > / С1 -SQH**2 )

137. DPSIH=-P11жС DPS1+DPS2+PI*SIВ-2ЖС0ВЖALC)

138. DPS I B=-SQP* С PI *C0B+2*S I Вж ALC+2:* < COB/SIВ > >1. F=FI»-:*:2+pSI*"*:2

139. G <2 > =2ж<FI«DFB+PSIжПРЗIВ >

140. G С1) =2ж < FI meDFH+PS I ж-DPS IH )1. RETURN ENDо. 4; О А; 0.6; 0.8; /.О40 4.2 4. А 4,$ 2.0

141. Рис.П.1. Номограммы для определения величины ^ ~ ^ ( "а")

142. Продолжение приложения 2 Блок-схема программы "Л/у 41. Ввод параметров |1. С <0.351. Нужно считать1. Стоп.ч по программе "М! Р&2 " /^0.45

143. Выбор из таблида с использованием линейной интерполяции

144. Вычисление промежуточных параметровI^ 0.5291. Печать таблицыТ

145. Запись таблицы в дисковый файл1. Да1. Занесение значений1. Ъг/Ь, оСв таблицу1. Стоп}- ХУ3 т1. Продолжение приложения 21. РкСиЗкйИ il.JP 01расчет таблицы относительного иАссоеого расхода е зависимости от относительного перепАДА давления

146. НА КЛАПАНЕ ДЛЯ РЕЖИМА СЕКЦИИ < 1 А>й.45-К=1.4>

147. РОРИй Т < /2Х' •' 3 А Д Ай ЗНАЧеНИЕ О •••"А' С Р4.2 > > АССеРТ 2' ОАтуре 10

148. Р0РИнТ С-"' 2Х'' ЗАДАЙ ноиер ЗАПИСИ 8 ДИСКОвОИ ФАЙЛЕ Нь, <13 > АссеРт 1 ЬНР РОР:ИАТ<13> туре 12

149. Р ОРИ А Т < .••-2Х' ЗАДАЙ Н АН А ЛьНОб ЗНАЧеНИЕ И8 ' < Р Р. 5 Г > АССеРТ 13'АИ0 РОРИ А Т < Р Р«5 > ту ре 14

150. Р0Р:ИАТ<У2Х'' ЗАДАЙ 8С ЛИЧИНУ ШАГА ПО И8'<РР.5>"> АССеРТ 13'8И0туре 18 й!=1й--"<ей-1 > 1Р(й1. !е.0.з 5>60то б1. < Й1. бе. 0.45 > 80 ТО 8 ТУРЕ 4

151. РОРМй Т < /2Х' •■ Н£':-'С Т0ЙЧИ8ЙЯ 06 ЛАС ТЬ ' .--2Х'' С леДУе Т Производи Ть

152. РАСЧЕТЫ ПО ДАННОЙ ПРОГРйИИе И ПО ПРОГРйИИе И.1Р02''•••••2Х ' -'еслину:+:ны РАСче ты по данной пРогрАине-еееди г> АССеРТ 5'к1. РОЗНАТ< 11 >гр <:к.ео.1>бото 8 туре г?й1

153. РОРИй Т < /2Х' Й1 = ' Р4.2' 2Х'' Ре:+!ИИ ЗА ТОП ПЕНИЯ С ТРУИ' С ЛЁДУ6 Т

154. Производить РАСЧЕТЫ ПП ПРОГРАМИЕ I-1.JP02') ЗТОРсонпнуе

155. С Аи 1 IN ТРР С Т ое' 8' 4' Т О' Г6' О А 8 А' 020 > М28=1/0281. U2 82=^28**2 ——~ 1. ( Ой. LE . 08 > GO ГО 58 Н=<0.168+0.38*Ой>60ТО 51 Н=8.59*0й CONTINUE R 1=03. Р97-Q. 19 Р*8 А >*<1+Н>*3.141 0й2=0йжж2г=г.828*н*а+ойг>

156. PS I1=R1*DA2.'"R2 ЦИКЛ ПО 110- assiш- ш>.1----------1. ASSIGN 32 ГО ISUH J=1сонпние1. ЙИ@=ЙИ0+8И0й1'18г=й1'18*ж2 й!'1=1-н.г*йие 2 АИ12=йИ82.'йИ РИ1=50РГ<АИ 12) РЛ^ЙИйжцИ:**:?. 5

157. Р28=<1-0.РжйИ12* < U202-1>)жйИжжЗ.5 U21=PSI 1*U202*P 20жжС-0.285PI43>иг=и20+йИ1г*иг1lv|22=U2s*!-*2

158. DEFINE FHE 51 (. 238»115»U»IN >ijJR I ТЕ < 51 •' NF > NF»й! * 8й f Ой * TM> TM1» Т J * Т 02 * Т DP > йиРй1. Сйи CLOSE<.'51 >1. STOP1. END

159. SU8R0U ТINE LINTRP <T 08»L D»L8»T 0 * T8»Ой * 8й»F >

160. ПИНЕйИйЯ ИНТЕРПОЛЯЦИЯ Ф'-'НКЦ ИИ A8VX ПЕРЕ MEHHblX» ЗАДАННЫХ ТАБЛИЦЕЙ. DI HENS I ON Т D8 С1 г 1 > * Т 0 < 1 > ■■ Т8 < 1 >»FR ( 4) DO 1 J=1*LD R=TD<J>1.(Dh.L E.R> GO TO 21. CONTINUE1. JU=LD1. GO ТЛ 51. JU=J1.CJiJ.EQ. 1> GO TO 51. JL=JU-11. DO 6 J=1'L8

161. FR с ,J > = T D8 ( JL f J) + С T D8 < JU J > T D8 < JL ' J > > •■■"2.

162. GO TO ? DO 8 J= 11 i8 FR<:.J> = TD8(JU' J> DO 9 J=1»L8 R=T8<J>1РС8Й.L£.P> GO TO 181. CONTINUE1. JU=L81. GO TO 131. JU=J

163. OJiJ. EQ. 1 > GO TO 13 •JL=JU-1

164. F=FR с JL > + < FR <; .JiJ > -FR < JL > > .-'2. GO TO 14 F= FRCJU) RETURN END

165. Продолжение приложения 2 Блок-схема программы "/V Аиг"1. Ввод параметров!0.451. Нужно считатьпо;0.35чЧСТОПзначенияу^ па программе ^^^

166. Вычисление промежуточных параметров1. Печать, таблищ *

167. Запись таблицы в дисковый файл1. Занесение значенийосв таблицу1. PR06РЙМ ¡-1JP 02рйсчег гнблииы относ итепьного массового рйсходйб ЗЙ8ИСИМ0СТИ ОТ ОТНОСИ ТЕЛЬНОГО ПЕРЕПйДй ДЙ6ЛЕНИЯ Нй КЛйПйНЕ ДЛЯ РЕФИМй ЗйТОПЛЕНИЯ СТР VH. (. i-yL<8.35»К= 1.4)

168. ОI MENS I ОН ТИ1 ( 9 )' ТJ 9 >» Т 02 < 9 )» Т DP < 9 ) <■ ТМ ( 9 )' ТР < 8) к' й!Рй ( 9)

169. Ой Тй ТР-"'в. 885' 8.81» 8.82»8.85 * 8. Ь 8.2»8. 3 8.4-''1. DO 25 J=1'91. TMl < J)=0.1. TJ<J)=0.1. TO2CJ)=0.г dp(j)=8.1. Ti'KJ>=8.й'Рй • i)=R.сонпние1. K=8 TVPE 1

170. POR ИйТС-"'2>:!* •' ЗйДйй ЗНйЧЕНИЕ 8-"'A» <Р4.2) •') й ССЕРТ 2'8й P0RMATí'P4.2) TVPE 3

171. FORMA Т (.••'2>:! i •' ЗйДйй ЗНйЧЕНИЕ О-"'А >■ < Р 4.2) ) йССЕРТ 2-D й TVPE 18

172. FORMAТ(2>-:' ■'ЗйДйй НОМЕР ЗйПИСИ 8 ДИСКОВОМ ФАЙЛЕ ИР' (13) •' >

173. ЙССЕР Т 1Ь НР FORMAТ<13) TVPE 12

174. FORMA Т С--'2>:!ЗАДА й НАЧАЛЬНОЕ ЗНчЧЕНИЬ М8» (F?. 5) •') йССЕРТ 13'йМ8 FORMAТ< Р?.5) TVPE 14

175. PORII йТ(/2•' ЗйДйй 8ЕЛИЧИ HV ШйГй ПО М8'<РР.5>-') йССЕР Т 13 <• Si'10 TVPE 18 й1=0й--"(8й-1) IF<AL.LE.О.35)GOТО 8 IP<й!.6E.0.45)60ТО ь TVPE 4

176. FORMAГ(у'2>'' •' HEVCТ0ЙЧИ8А Я ОБЛАСТЬ' <■ .••'2>í<' СЛЕДУЕТ ПРОИЗВОДИТЬ ж РАСЧЕТЫ ПО ДйННОй ПРОГРАММЕ И ПО ПРОГРйММЕ И JPD1 •'' •••"2>-,'' ЕСЛИ ж Н':-':+!НЫ Р'йСЧЕТЫ ПО ДйННОй ПРОГРАММЕ*88 ЕДИ Г') йССЕРТ 5'К FORMATO!) IF<K.EQ.i>GÜTO 8 TVPE P-AL

177. ЦИКП по 1-10 ' ^-г.-,.^-^--^.,-- « 7

178. ASSIGN 38 ГО ISU ASSIGN 32 ТО ISUi J=1

179. CONTINUE AM8=AM0+SM8 А И 8 2=AM0$$2 йИ=1-0 .2жйИ82 AM12=AH02/AH Ri'1 i=SQRT<Ai'1l2>1. Р:.1=й1,10жй1,1жжг. s

180. Р28= < i -8. Р:*йИ 12$ < U282- i > > .*йИ**3. 5

181. U2i-=eSI-lJKU2023i®2Q**s{r:0-« 285£i4JXиг=иг8+йн i2*U2 i U22=U2sís2 RP2= < i -8.2* AM i 2'#. < U22 -1 > > ШЗ. 51. RiP=i-P21. < R DP .LT. TP < J > ) GO TO 24 GOTO ISUI» <30 »31) Tili". 1>=АМй1. Tili <J)=RM1 TJ < J>=RJ

182. R02=P2^'*8. 714285 ? T D2 < J > = 1 < R02*U2* ОЙ > TDP<J>=RDP

183. ALFA < J > =. 836 7* < P. I.*' < 1 / Ой > ■■•'SQR T CR DP > >1. GOTO ISb.il ' <32' 33>1. J=J+11.<J.EQ.9>й3516Н 31 TO ISbJ CONTINUE1.<P2.GT .0.529>60TO 24 ASSIGN 33 TO IЗУ1 GOTO 30 CONTI Ni IE TVP'E 16

184. PORi'1 A T < --'2H ' ' EC ПИ Hv'.tlH A P АСПЕН A TK A ' 8К ЛЮЧИ ж ПРИНТЕР И 88ЕДИ Г"> К =й

185. ACCEPT 5'К IRX.tQ.EOGOrO I? I.1.IR ITt <6' 18>

186. FORMAr</25>í' ПР ОГРАМИА'1 ' 18>í' 'MJPD2' > b.lR I ТЕ < 6 ' 19 > AL ' 6А » DA ' NF FORM A T < •••"SX ' ■' 1 /L= ' ' F4.2 ' 5 >' ' ■' 8A= ' ' F4.2 ' $5X' •' DELTA/A=-' ' F4.2' SX' •'NF=' » I3>h.lRI TE<6' 20 >

187. FOR Мй f < /2% '110 < •' ~ •' > >i.'JRITE<6'2 i>

188. URIГЕ<52' NF>HFp й1» 8йр Ой» ГИ» ГИ1 » t\j» Г02» TDP» й1Рйcall close < 52 >1. STOP1. END

189. Продолжение приложения 3 Блок-схема программы г/71. Нужно считать по программе руг-Организацияшкяа по1. Вычисление значения1. Нет1. Нет1. Занесение в таблицу

190. Запись таблицы в дисковый файл3.0.5291. Р< 0.5291. Стоп- 2Ü21. Продолжение приложения 31. PROGRAM PY.I

191. РАСЧЁТ XlKíPhi ДАВЛЕНИЯ Ни ПОВЕРХНОСТИ ЗАПОРНОГО ОРГйИй ПЁР8ЫЙ 8ЙРИЙН1" s 1 .••''!. )'=@ . 45

192. DIMENSION ТМ<9> * TUi (9) » ТР0<9> * ТР<9> ? TU8<8> * TDAíS) * TUYA<8* 8> ? Г8ЁТ(1 i » 26)? ТV8(9 >? ТD <1 i > *Т8(20 >

193. Ой TA T D--'8. 1*8. i 5 * 8 „ 2 * 8. 5 ? 8 = 4 * 8.5 ? 8 = 6 ? 8 = P * 8.8 * 8. 9 * i./

194. ОйТй те/1.05*i=1? i.15? i.2?1 = 25* Í.3?1.35? i.4* i.45* i.5*1.55?1.6?si.65*i.P* i.P5* 1.8*1.85?1.9?1.95*2./

195. FORMAT! ACCEPT 2* 8й FíjRMAT<F4. 3) TYPE 3r O R i 1A T ( / 2 * 3 A ДАЙ 3 H A H Ë H И E 0 / A * (64.2) ' ) ACCEPT 2*OA AL310 А/ ( 8A~ 1 > IP < A L « 6 Ë . 8.4 5 > G G T 8 4 IP< AL.LE.8.35 > GО T8 5 TYPE 6

196. FORMA T </2Х * ■' HE Y С Г0ЙЧИ8 АЯ OS ПАС Г и * $CЛЕДYET ИЗМЕНИТЬ ВВОДНЫЕ ПАРАМЕТР*!'' ) S TOP1. TYPE P * ALi.56 P * i.564 ? i. t— íüT ! .J.J I . 555 * i.549 ?

197. FORM и Г ( 2 X ? ■' 3 й Д A й НОЧЕР ЗАПИСИ 8 ДИСК080М АйЛЕ 4F?< 13 > } йССЁР Г 9 ? MF FORMA ТЧ13) TVPE i 0

198. FORMA Г < •••"ЭХ p ■' ЗАДАЙ НАЧА ПьИОЕ ЗНАЧЕНИЕ И8 ? ( F 7 , 5 > > АССЕР Г i Í ? AM0 FORMA Г ( F7 = 5 > TVPE 12

199. F0Ri,1Ar<-y2X' •' ЗАДАЙ ВЕЛИЧИНА ШАГА ПО И0р 5> > ACCEPT 11 ' Si'10 TVPE 23

200. FORMA Г < -••'2>:! p •' ЗАДАЙ H A4 A ПьИОЕ 3H A4 ьНИЕ T p ( F P » 5 ) •' > ACCEPT ib Г .1. TVPE 22

201. FORMAT«У2ХЗАДАЙ вЕПИЧИН':-' ШАГА ПО Tp<FF\5>'>1. АССЕР T I i p D Т1. < он.LE.0 =8)60 го 13804=0.168+8а 38«DA1. GOTO 14304=0« 59»; DA

202. CALL L ï H TRP Г8Е T f i 1 » 20 r TD? Г8 ? DA * 8A * 8E Г )sie=sin<eer)1. H=8QH^:»;21. U8Q4= i ••••'SQ41. DA2 = D As; ж 21. C08=C0S<8ED1. S4P=SQ4--'3.141592 78Ы60Р 8EГА 8 ЦИКЛЕ ПО Г1. DO 50 1 = 1 ? 9f=T+DTl'f< г. бе» 1. >60to 194 г=4;*; г

203. SQHT=SQR Г(4 Г) SQT=SQRГ < Г > IF (H. LE. D60T0 16 ARC=A TAH ( 2&SQH T--' < 4- D > GOTO IP TVPE 20

204. FORMAT(••••'2Xp ■' T=1 ■' >5 ГОР

205. ARC=3 „141592 7-А ГАМ < 2&SQ4 Г/ ( T-H ) ') ARC 1 = A TAN ( 2 ж SQH Г--' < 1 -H f ) )

206. Г s; S18+ Г R1 = ALOG < ( 1 "R > •■•"' ( 1+R > ) ARC2=A TAH < 2&SQ Г»;СОв--' < 1 Г ) )

207. V=S4P»; < SQHsARC 1 +USQH&ARC-S 18&R 1 -2®C08*ARC2 -"'вAif<v.lí. rvea;060Tü is

208. UV0=SQR Г ( Г •••'4 > ЦИКЛ ПО М0 ASSIGN 30 ГО IS4assign з? m rsiiüj = 1f?|--|0=Al'10

209. К i '18 =R И @+S i 'Í8 . ^

210. R i '1 ñ P ~ R í i 0 ж ж SйМ= 1 -8. 2&RH02йМ i =RM02-''AM1. РЙИ^ЙИЖ^З* 5

211. Р@=РЙИ1ф; ( i -Q „ ps-йМ i ж ( i.j V8&s;2- I ) )1. R=8. P9F-8. Ï9P»:6h

212. R1 = < í +H > œ 0 A2&3.141592 71. R2=2.828$H& <1 + DH2)1. P S11=R ж R i / R 21. У8Й:~УЛ8Ч

213. FORMAT^** 'ЕСПИ HV:+!HA РАСПЕЧАТКА* 8кПН)ЧИ ПРИНТЕР И 88ЕДИ 1 •') к;=8йССЕР." 35?К1. FORMAГ<I1)1. < К » EQ.8> 60 ГО 36b.lR IТ£ < 6 г 3 7 ) hL ? 8 А » Ой » Г » V » HP

214. FíDRíIh Г </2% f •' 1 •••"i .= •' * F4,2 p 5>:! r ■' 8--'й=•' r F4„ 2 ? ? •' DEL Гй.--'й= * F4 r 2 ? 5>iя ' Г=' V.--8=-' »Ff.S'Sïi» •'NF= •' » 13)

215. HRI TE <ь ? 38) FORMA Г < ■•••'2X * < ' -- ' ) ) UPI TE<ь'39)

216. FORMAT< l.2>-,> ' l'V г lbX» ¡II » IbXf P0- " 7>-, <■ P-' )1. UR I TE < ь <■ 38 )1. DO 41 J=ï » 91. AM»TH<U)1. AM1 = ni I(J)1. P0=TP0< J)p=r TP ( J •;,ilp I TE < 6 ' 42 ) AN * AM 1 ? P0 f V f P P0R|'1AT(4E18. P) HR I TE < ь f 38 ) СОНГÍNUE

217. UiRI ГЕ< ?*^JüL>8üi Ой» Т* ~

218. D£P I Н£ г IL£ 53 < 1 @@« 83 и f I Н )i.i.SR I Те < 53 ■' HP > HF i. AL 8й - Он? ГМ ? ."¡1 i р ГР@ * Г Р * Vcall close <53 >1. HF=NF+Iсонпние1. STOP END

219. SU8R0U т i не seeue < т ? нй »не ¡нй ? те * й е ? п >■ го >1Х0и!дение одной из координат двумерного мйссивй

220. ЗйДАННОЙ 8ТОРОЙ КООРДИНАТЕ И ЗАДАННОЙ':'1 1А Ч £ Н И н) i 1 й С С И 6 А

221. DIMENSION Т<1» 1)" ГА< i) * те<1) IР<НА.£Q . 0> GO TO 1 DO 2 J=1 * NA R=TA<J)1.<A«L£,R>GOTO 3сонпние1. DO 4 1=1»N6 r= у <; , J )1.<TI.L£.R > GОT0 5 СОН TIHUE T0=T<Ü-4p I.) GOTO ь DO 7 ,J= i t N8 R=T8< .J>1.<8.L£.R>GOTO 81. СОН TIHUE1. J8=J

222. DO 9 I = i r N A R= Г(IrJ8)1.<TI.LE.R)GОГО 10сонпние то= r <ifje >1. RETURN END

223. Продолжение приложения 3 Блок-схема программы "руго.45 / ч Нужно считать * ч222Шзначенияпо программе РУ41.<0.35

224. Организация цикла по -Ь при £ «= |1. Вычисление значения )1. Нет

225. Организация цикла по пс ■ ♦ ~1. Вычисление р1. Нет1. Занесение в таблицу

226. Запись таблицы в дисковый файл1. Р £-0.5291. Р <0.529.1. РРОбРЙН Р'т'2рцсчь г эпюры давления на поверхности запорного органа

227. ГОРОЙ 8 АРИАН Г г 1 /1< =8.35

228. ОI НеНЗ I ОН ГН < У > р ГН1< 9) <■ ГР8 >:.' У > <• ГР < 9 > » Г08 (8 > р Г О А (8) р ГО УА (8 р 8 > р з*з ГУ 8 < 9 >

229. РОкНА Г < /2Х р ■' ЗАДАЙ ЗНАЧеНИе О/А р < Р4.2 >)йсссрг 2»ой

230. РОЗНА Г ( /2х р ■' ЗАДАЙ НйЧй ПьНОЕ ЗНЙЧЕНИЕ Г» < Р Р. 5) •') АССЕР Г 11» Г ТУРЕ 22

231. Р03НАГ</2Х» ЗйДйй 8Е ЛИЧИНУ ШАГА по Г»< Р Р. 5) •' > АССЕРТ 11»ОТ 8ЕГ=1.5Р8Р963 ЗОН "У. 544з*з ОА/803 Г (й1.+8. 45) Н ~~ о К! Н ^ -с!изон=1/зон 0Й2— Ойз*$2

232. ЙЗНР=8ОЬР-'3» Г4Г592Р 8ь!60Р 8ЕГА 8 ЦИКПЕ ПО Г DO 58 I = i F9 Í=Í+DÍ1.- с; T.6E» 8.995>60TO i9 НТ=НжТ

233. S Q И Г = 8 G! К Г < H Г ) SQ^SQRK Г> IЕ < H « LE.Г>60 ГО ib ARC ~A ТАН ( 2?iíSQHT--' ( H- Г > ) 60TO IP TVPE 28i?@P>:; p •' f=8.995'' > STOP

234. ARC=3 Л 4 i 592 Р-й TAH ( 2«SQH Г--' < Г-Н > >

235. ARC i =A TAH < 2ж8ЙН Г/< i-H Г >>

236. R i=2жAL06 i +SQ Г > ( i -8Q Г > )у=5НРж ( SQHsARC i +USQH$ARC+R I > •••"'8A1. P(V „LT. TV8<I>>60 TO i 51. UV0=SQR Г < Г.-'И >цикп rio lie1. ASSIGN 38 ГО ISU1. ASS I GH 32 ГО IS I" Л1. J=i

237. Ri-ig=A г-IQ pH0=Ri'1ñ+S¡'18 РН02=РН0жж2 AH= i -8 . 2&RH02 AH.l=RH02/AH PAi'i=Ai'1^s3 г 5

238. Р8=РАИжi -ñ„ РжАИí ж<0жж2-i > )1. R=8.P9P-8*19Рж8А

239. Ri = < 1+Н)ж0й2жЗ. 141592P1. R2=2.828жИж<i + DA2 >1. PS 11 =ржр i /R21. V8A=Vs8A

240. CALL SEEU8 ( TUVA > 8 > 8 я TU8 * Г DA ? 8 » <• DA ? V8А >■ U8 > R=(SQR Г <1 + DA2)>ж<1-U8> R i=SQRT< i + ОА2-1!0жж2 > ж<1+U8 > PSÍ2=R--'RI

241. UV 1 =PS I i ж С i -PS 12 > ж < UУ0жж2 ) жР0жж < -0 „ 285 Р .143 > Р=РЙИж ( i -0, 2ж АН i ж < UV 1 жж2-1 > > жжз . 5 RP=i-P1.(RP.LГ.TP(J >>60 TO i 860Г0 ISHp<38P3Í)1. Ti'1<:.J>=RiÍ01. TU i < J ) = Ail i1. ГР0<J)=P060Г0 ISuU р<32p33)1. J=J+i

242. Р0Р!-1ЙТ(/2Хр •' , Р4s 2р 5::<р •' 8--н=' р Р4» 2? 5>:!Р ■' ОЁ1ТЙ/Й=-' Р4» 2? ;д>:;Р ■■ .-=•■ ,РР„5р5>:!Р •'V-"'8=' '£7«5>5-<' •'НР=-' р 13) URITE (6р38) Р0Р1'1йГ< --'2Чр 76 (•' -•' ) ) URI ТЕ кь р 39)

243. Р0Р!'1йГ< 12>ip 1'Г р 1ь>:!Р •' ill р 1ьН> PQ' р Ftf' ' Р-')1. URI ТЕ < ь р 38)1. DO 41 U=1p91. ЙМ=ГИ<J)1. ЙН1 = Г И i (J >1. P0=1"P0<.J)1. P= ."p (j)

244. UR I ТЕ < 6 p 42) f—ii*1 p ЙН1 p P8 p V p P Р0РИйТ<4Е18«Г) URITE(6p38) СОНТIHUE

245. URI ТЕ < 7p 3 7)hLp 8й p Ой p T»V * HP

246. OEPIHE PILE 53<168p 83p UpIH)

247. URI ТЕ(53' HP)HPp й1р 8йр Ой' T' THp THi P ТР0' TP' Vcall close<53)1. Mp.-f.jp+ j1. CONTIHUE1. STOP1. END