автореферат диссертации по энергетическому, металлургическому и химическому машиностроению, 05.04.13, диссертация на тему:Исследование гидродинамики течения рабочего тела в однокольцевом самодействующем клапане поршневого компрессора

ВЩЦЕНИЕ.

ГЛАВА. 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может изменить направление.

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