Помощь в написании студенческих работ
Антистрессовый сервис

Коэффициенты диффузии в многокомпонентных системах при испарении бинарных азеотропных смесей в инертный газ

Диссертация: Коэффициенты диффузии в многокомпонентных системах при испарении бинарных азеотропных смесей в инертный газ
ДиссертацияПомощь в написанииУзнать стоимостьмоей работы

Эти задачи можно успешно решить, если будет получено более точное описание явлений и внутренних закономерностей, определяющих процессы массопереноса. В значительной мере это касается и описания лимитирующей стадии процесса переноса массы — процесса диффузии. Существующая в настоящее время практика расчета коэффициентов переноса массы в таких массообменных процессах, как ректификация, абсорбция… Читать ещё >

Коэффициенты диффузии в многокомпонентных системах при испарении бинарных азеотропных смесей в инертный газ (реферат, курсовая, диплом, контрольная)

Содержание

  • Глава 1. Коэффициенты диффузии многокомпонентных систем, методы их расчета
    • 1. 1. Термодинамика необратимых процессов
    • 1. 2. Расчет многокомпонентных коэффициентов диффузии по кинетической теории газов
      • 1. 2. 1. Метод Максвелла — Стефана
      • 1. 2. 2. Метод Уилка
      • 1. 2. 3. Метод Тура (метод линеаризации уравнений)
  • Глава 2. Экспериментальные методы исследования многокомпонентных коэффициентов диффузии
    • 2. 1. Метод Стефана
    • 2. 2. Метод Грэма — Максвелла
    • 2. 3. Метод Лошмидта
    • 2. 4. Хроматографический метод
  • Глава 3. Определение равновесных составов в системах бинарная азеотропная жидкость — газ при нормальных условиях
    • 3. 1. Выбор исследуемых веществ и их краткая характеристика
    • 3. 2. Экспериментальная установка для определения равновесных составов многокомпонентных парогазовых систем
    • 3. 3. Исследование изменения состава бинарной азеотропной смеси в процессе испарения и по пути диффузии
  • Глава 4. Экспериментальная установка и методика измерения многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в инертный газ
    • 4. 1. Выбор метода исследования
    • 34. 2. Описание экспериментальной установки
    • 4. 3. Методика проведения эксперимента
    • 4. 4. Расчетная формула и выбор системы отсчета
      • 4. 4. 1. Оценка погрешности опытов
  • Глава 5. Результаты измерений многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в инертный
    • 5. 1. Результаты измерений
    • 5. 2. Сравнение опытных значений полученных многокомпо- 71 нентных коэффициентов диффузии с теоретическими
  • Выводы
  • Список используемой литературы

К точности выполнения расчетов, связанных с процессами тепломассообмена, предъявляются все более высокие требования. Это вызвано рядом задач, решаемых при создании современной тепломассообменной аппаратуры, широким применением в расчетах ЭВМ, а также разработкой надежных методов инженерного расчета высокоэффективных тепломассобменных аппаратов и установок.

Эти задачи можно успешно решить, если будет получено более точное описание явлений и внутренних закономерностей, определяющих процессы массопереноса. В значительной мере это касается и описания лимитирующей стадии процесса переноса массы — процесса диффузии. Существующая в настоящее время практика расчета коэффициентов переноса массы в таких массообменных процессах, как ректификация, абсорбция, экстракция, сушка, [1, 2]базируется на использовании приближенных данных о коэффициентах диффузии.

Картина переноса массы в многокомпонентных системах, где возможны такие предельные явления как диффузионный барьер, осмотическая и реверсивная диффузии, довольно сложна. В таких системах полную картину о механизме переноса массы может дать только знание матрицы многокомпонентных коэффициентов диффузии. (МККД).

Несмотря на то, что в общей теории диффузионных процессов к настоящему времени достигнуты известные успехи [3, 4, 5], теоретические и экспериментальные исследования в области диффузионного массопереноса в многокомпонентных системах далеки от завершенности. Это, с одной стороны, связано со сложностью математического описания процессов переноса в таких системах, а с другой — трудностями экспериментального характера, что обусловлено рядом обстоятельств. Во-первых, методики проведения опытов, 5 на основе которых можно было бы проводить расчет МККД, существуют, в основном, для газовых и жидких систем, во-вторых, для парогазовых систем отсутствуют расчетные уравнения, применимые к конкретному методу исследования.

В настоящее время наиболее полно решена задача динамики частиц с потенциалом взаимодействия «твердых сфер» [5−8], которая, несмотря на упрощенность потенциала, может служить первым приближением для описания свойств реальных газов и жидкостей.

Для описания реальных систем, очевидно, необходимо использование различного рода априорных, феноменологических предположений о динамике частиц, таких, например, как принцип независимости диффузионных потоков [9] или принцип молекулярного хаоса [5, 10, 11]. Подобные предположения позволяют значительно упростить математическое описание взаимодействующих частиц и получить конечные выражения для кинетических коэффициентов через динамические параметры системы.

Как показывает анализ работ, посвященных переносу массы, исследование диффузии в многокомпонентных системах представляет собой сложную задачу [12−24].

На основе строгой кинетической теории разработаны уравнения, описывающие перенос массы в многокомпонентных системах. Эти зависимости проверены только для распределения концентрации компонентов для газов близких по своему состоянию к идеальным [22, 23], в которых роль перекрестных коэффициентов в общем переносе массы мала. Для неидеальных систем и систем под давлением мы располагаем пока линейными соотношениями термодинамики необратимых процессов, для использования которых необходимо знание феноменологических коэффициентов.

В периодической литературе в настоящее время отсутствуют работы, в которых были бы приведены значения МККД для парогазовых систем. 6.

Этим объясняется необходимость глубокого экспериментального и теоретического исследования многокомпонентной диффузии, результаты которого могут быть использованы для удовлетворения запросов практики и уточнения вопросов теории.

Целью данной работы является:

1. Определение равновесных составов в системах бинарная жидкостьпар и оценка распределения концентрации паров бинарной азеотропной смеси по пути диффузии.

2.Разработка методики экспериментального исследования и теоретически обоснованного метода расчета многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в газ.

3.Проведение комплексных экспериментальных исследований для определения многокомпонентных коэффициентов диффузии для интервала температур 295. .ЗЗОК и нормальном давлении.

4.Численное моделирование прямых и перекрестных коэффициентов диффузии при испарении бинарных азеотропных смесей в инертный газ при температуре 295. .ЗЗОК и нормальном давлении.

Научная новизна заключается в следующем:

1.Определены равновесные составы в системе жидкость-пар для бинарных азеотропных смесей: метанол-бензол, этанол-бензол, н-пропанол-бензол.

2.Проведена оценка распределения концентраций паров бинарных азеотропных смесей метанол-бензол, этанол-бензол и н-пропанол-бензол по пути диффузии и установлена нелинейность распределения концентрации компонентов бинарного азеотропа в газ.

3.Предложена методика экспериментального исследования и теоретически обоснован метод расчета многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в газ. 7.

4.Впервые определены прямые (Dn, D22) и перекрестные (Di2, D21) коэффициенты диффузии при испарении бинарных азеотропных смесей в инертный газ.

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

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

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

Личное участие. Все основные результаты работы получены лично автором под научным руководством профессора Я. Д. Золотоносова.

Апробация работы. Основные положения работы были представлены на ежегодных научно-технических конференциях КХТИ им. Кирова в период 1972 ~ 1975 г., 1979., 1981 г., 1983 г., 1986, 1988, 1991; на V Всесоюзной конференции по теплофизическим свойствам веществ, г. Киев, 1974 г.- на Российском национальном симпозиуме по энергетике, РНСЭ, г. Казань, 2001 г., на XXVI Сибирском теплофизическом семинаре, г. Новосибирск, СО РАН, 8 июнь 2002 г.- на Всероссийской школе — семинаре академика В. Е. Алемасова «Проблемы тепломассообмена и гидродинамики в энергомашиностроении», г. Казань, октябрь 2002 г.

Публикации. По теме диссертации имеется 7 публикаций. Объем работы. Диссертация изложена на 103 страницах и состоит из введения, пяти глав, выводов, списка литературы из 100 наименований и приложений. Иллюстрационный материал содержит 17 рисунков и 6 таблиц.

Выводы.

1. На основании анализа литературных данных установлено, что из всего многообразия разработанных на сегодня способов определения коэффициентов диффузии? наиболее информативным и надежным для измерения многокомпонентных коэффициентов диффузии в газовой и парогазовой средах является стационарный метод Стефана.

2. На базе статического метода измерения равновесных концентраций бинарных паров, впервые, для интервала температур 295.330К и нормального давления получены экспериментальные и расчетные значения равновесных составов многокомпонентных парогазовых систем: метиловый спирт — бензол-азот, этиловый спирт — бензол-азот, н-пропиловый спирт — бензол-азот.

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

Установлено, что при испарении первых 10 весовых процентов бинарной азеотропной смеси, состав жидкой и паровой фаз практически не меняется, а концентрация паров бинарной азеотропной смеси по пути диффузии меняется нелинейно.

4. Посредством стационарного метода Стефана разработана методика, создана экспериментальная установка и проведены измерения многокомпонентных коэффициентов диффузии паров бинарных азеотропных смесей в инертный газ.

Показано, что рассматриваемые системы являются постоянно испаряющимися, а предложенная методика может быть использована для измерения МККД.

5. На базе феноменологических уравнений диффузии описан процесс многокомпонентной диффузии паров бинарной азеотропной смеси в инертный газ, определен характер изменения концентрации компонентов по пути диффузии и впервые рассчитаны многокомпонентные коэффициенты диффузии: прямые — Dn, D22 и перекрестные — D12, D21, представленные в виде таблиц.

6. Проведено сравнение прямых коэффициентов диффузии — Дни Q22с бинарными коэффициентами, известными из литературных источников и рассчитанными нами. Показано, что сходимость прямых коэффициентов матрицы МККД с бинарными коэффициентами диффузии не превышает 12%.

Показать весь текст

Список литературы

  1. Рид Р., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей. JL: Химия, ЛО, 1982−592с.
  2. Т., Пигфорд Р., Уилки Ч. Массопередача. М.: Химия, 1982 — 696с.
  3. П., ДеЛенер М. Классическая кинетическая теория жидкостей и газов. М.: Мир, 1980 — 423с.
  4. Р. Равновесная и неравновесная статическая механика. М.: Мир, т.1, 1978−405с.
  5. Дж., Кертисс И., Берд Р. Молекулярная теория газов и жидкостей. М.: ИЛ, 1961 — 930с.
  6. Van Beijeren, Ernst М.Н. The modified Enskog equation for mixtures//Physica. 1973. v.70.-p. 225−242.
  7. Haro M.V., Cohen E.G.D. The Enscog theory for multicomponent mixtures. I Lineor transport theory//J.Chem. Phys. — 1983. v.78, № 5. — p. 2747−2759.
  8. Van Beijeren. Equlibzium distribution of hard spheres systems and revised Enskod theory//Phys. Lett. — 1983. v.51, № 7 — p. 1503−1504.
  9. Дж., Капер Г. Математическая теория процессов переноса в газах. М.: Мир, 1976. — 554с.
  10. С., Каулинг Т. Математическая теория неоднородных газов. М.: ИЛ, 1960−510с.
  11. Л. Лекции по теории газов. М.: ГИТТЛ, 1953 — 556с.
  12. Wilke C.R.//Chem. Eng. Progr. 1950. v.46. — p.95.
  13. Fairbanks D.F., Wilke C.R.//Ind. Eng. Chem. 1950. v.42. — p.471.
  14. Gilliland E.R., Sherwood T.K.//Absorbtion and Exstraction//New York 1937. -p.ll.
  15. Toor H.L.//A.I.Ch.E.Journal. 1957. v.3, № 2. -p.198.
  16. M., Boas A. // Chem. Eng. Progr. 1951. v.47. p.51, p. l 1183
  17. Cichelli M.T., Weathezford, Bowman J.R. 11 Chem. Eng. Progr. 1951. v.47.p.63, p. 123
  18. Han Hsien Wen, Bird R.B. I I A.I.Ch.E. Journal. — 1960. v.6. p.516, p.551
  19. S.A. //A.I.Ch. E. Journal. -1961.v.7.p.l7
  20. E.L., Lightfoot E.M. // A.I.Ch. E. Journal. 1963. v.9. p.703, p.783
  21. W.E., Prober L. // Ind.Eng.Chem.Fund. 1964. v.3. p.224
  22. H.L. // A J. Ch. E. Journal. 1964. v. 10.p.448, p.460
  23. H.L. //Altech.- 1962−1963. v.12. p.21−43,p.51−61
  24. P., Стьюарт В., Лайтфут E. Явления переноса. М. Химия, 1974−688с.
  25. Н.К., Лундин А. Б. Термодинамика необратимых физико-химических процессов. М.: Химия, 1984−334с.
  26. А.Н., Бикбулатов А. Ш. Установка для определения коэффициентов диффузии паров в сжатых газах //Тр. Казан, хим.-технолог. ин-та. 1969. — вып.43. — с. 47−49.
  27. Де Гроот С. Р. Термодинамика необратимых процессов. М.: Гостехиздат, 1956.-370 с.
  28. Де Гроот С. Р., Мазур П. Неравновесная термодинамика. М.: Мир, 1964. -457 с.
  29. Д.Н. Неравновесная статистическая термодинамика. М.: Наука, 1971.-416с.
  30. К.П. Феноменологическая термодинамика необратимых процессов. -М.: Наука, 1978.-128 с.
  31. P.I., Gosting L.J. //Jann. phys. chem. 1959, v.63, p.86
  32. P., Dunlop P., Gosting L. //J. Am. Chem. Soc. 1955, v.77 p.5235.
  33. Fugita, Gosting L. //J. Am. Chem. Soc. 1956, v.78 p.1099.
  34. A.C. Испаряемость топлив для поршневых двигателей и методы ее исследования. -М.: Гостоптехиздат, 1955, 306 с.
  35. К., Sada Е., Nagasaka М. // Kaguku Kogaku. 1971, v. 35, № 9, p. 62.84
  36. Теплопередача /А.Гриб, Б. Уокил // Труды амер. об-ва инж. механиков. 1969. № 2- с. 62.
  37. Н.А., Кувшинский Е. Б. // ИТЖ. 1948. т. XVIII, вып. II. — с. 162.
  38. Ц.М., Померанцев В. В., Франк-Каменецкий Д.А. // ЖТФ. -1942. т. 12, вып. 1. -с.14.
  39. Stefan / J. Sitzungsberg. Akad. Wies. Wien. Abt. 2 1871, Bd. 63. 1873, 68, 385.
  40. Stefan / J. Ann. Phys. 1890. № 41, p.725.
  41. И.Р., Хазанова H.E. // Газовая промышленность. 1962. № 2. -c.29.
  42. G. /с. Scientific Papers, Dover Publications, N.Y. 1952. № 2. p. 57.
  43. J.J., Pigford R.L. / Chem. End. Sci. 1957, v.6, № 3 — 6. p. 215.
  44. Loschmidt / J. Akad. Wiss. Wien. 1870, Bd. 61, Abt. 2, № 3, s. 367 380- Bd. 62, Abt. 2, № 7, s. 468−478.
  45. Obermayer.// Akad. Wiss. Wien. 1880, Bd. 81, Abt. 2, № 5, s. 1102 1127- 1882, Bd. 85, Abt. 2, № 1 — 5, s. 147 — 168- 1883, Bd. 87, Abt. 2, № 1, s. 188 -263.
  46. Lonius / A.Ann. Phys. 1909. Ziff. № 29. p. 664. Ziff. № 13. p. 1414.
  47. J.C. /J. Chem. Phys. 1959. v. 31. p.1462.
  48. J.C., Seager S.L. /J. Chem. Phys. 1961. v. 35. p.2242.
  49. J.C., Seager S.L. /Ind. End. Chem. Fund. 1962. v. 1. p, 277.
  50. Huber. / J.F.K.Metodes Phys. d’Analise. 1968. v.4 № 3. p. 285.
  51. Физика химических растворов / Ким С. А., Сазонов M. JL, Лозгачев В. И., Жуховицкий А.А.// 1972. с.250
  52. В.И., Канчеева О. А. // 1972, ЖФХ, t.XVI. № 5. с. 1242
  53. М.А., Дьяконов С. Г. Расчет коэффициентов бинарной диффузии в плотных газах. // Труды КХТИ им. Кирова С. М. Казань., 1974. Вып.5385
  54. Walker R.E., de Naas N., Westenbergs. // A.A.J.Chem.Phys. 1960. v.32. p.1314
  55. Maxwell G.C. On the dynamical theory of gases // Phys. Trans. Roy. Soc. -1867. vol.157, №= l, p.49−88
  56. Maxwell G.C. On the dynamical theory of gases //Philos.Mag. 1868. vol.35, S = 236, p.185−217
  57. Экспериментальные методы определения коэффициентов взаимной диффузии в газовой фазе: Обзфинформ. / Бережной А. Н., Семенов А.В.- Госстандарт- ВНИЦ MB- Казан. хим-технол. ин-т. М.: Изд-во стндартов, 1986.-39с.
  58. Winkelmann. A. Ann. Phys. und Chem., XXII, 1, 1884, XXII, 154,1884, XXIII, 203, 1884, XXIV, 105, 1885, XXXIII, 485, 1888.
  59. Jl. Таблицы азеотропных смесей. М: ИЛ, 1951−292с.
  60. С.К., Лестева Т. М., Кочан В. Б. Азеотропные смеси, справочник. Л.О., Химия, 1971 .-848с.
  61. В.Б., Фридман В. М., Кафаров В. В. Равновесие между жидкостью и паром (справочник). М.Л.: Наука, 1966.
  62. С.Т. Свойство газов и жидкостей. М.Л.: Химия, 1966. -422 с.
  63. М.С. Работы по теории растворов, Изд. АН СССР, 1953
  64. Д.П. Об упругости пара растворов. Wied. Ann.(3), 14. 35, 1881-с.219.- С.-Петербург, ЖРФХО, ч. химич., I, 16
  65. А.В. Термодинамика гетерогенных систем. Изд. ЛГУ. 1967
  66. Hellund E.J. Phys. Rev., 57, 4, 737, 1940.
  67. Sazan A., Sing Y. Canad J. Chem., 44, 2222, 1966.
  68. Srivastava B.N., Paul RPhysica, 28, 646,1962.
  69. Cunningham R.S., Geankoplis С J. Ind and Eng. Chem. Fund, 7, 3, 429, 1968.
  70. Curyiss C.F.J. Chem. Phys., 49, 7, 2917, 1968.
  71. Onsager L.Phys. Rev., 37, 405, 1931.86
  72. Cussler E.L.J. Chem. Phys., 44, 7, 2829,1966.
  73. By Watts H.Canad.J.Chem., 43, 431, 1965
  74. Stewart W.E., Prober L.Ind.Eng.Chem.Fund, 3, 224, 1964
  75. E.H., Александров И.А. ИФЖ, XXIII, 6, 1012, 1972
  76. Д. Диффузия в многокомпонентных системах, (перевод) АН ЧССР, Прага, 1996
  77. Н.Б. Справочник по теплофизическим свойствам газов и жидкостей. М.: Наука, 1972, — 598 с.
  78. В.В. Азеотропия и полиазеотропия. М.: Химия, 1968. -348 с.
  79. Н.Д., Косов Н. Д. Физика, вып.1, 130 Алма-Ата, 1970
  80. Tarenko В.Т., Bicbulatov A.Ch., Zolotonosov Y.D. Measuring of the bilateral aseothropic mixture during the process of the evaporation // Proceedings of the Russian National Symposium on Power Engineering. Kazan. RNSPE, 2001. V.l. P. 234−236.
  81. Н.Д., Новосад Н. Д., Сб. Вопросы общей и «прикладной физики», Алма-Ата 1969.
  82. В.В. Курс дифференциальных уравнений.- М.: ГТЛ, 1950 254с.
  83. .И., Бикбулатов А. Ш., Золотоносов Я. Д. Измерение равновесных составов трехкомпонентных парогазовых смесей. // Известия вузов. Проблемы энергетики. 2002. № 7 8. С. 36 — 42.
  84. Е. Н. Александров И.А., Гореченков В. Г. Химия и технология топлив и масел № 5, с. 34 36, 1973.
  85. П.В., Шейдлин А. Е. Исследование термодинамических свойств веществ. Гоэнергоиздат, 1963.
  86. П.В. Физика высоких давлений, ОНТИ, М. Л., 193 5
  87. П.В. Новейшие работы в области высоких давлений, И.Л., 1948.
  88. С.Г., Казанцев С. А. / Определение коэффициентов диффузии в многокомпонентных смесях методом молекулярной динамики// Теоретические основы химической технологии. М.: 1981, т. XV. С. 20 — 24.
  89. А.Н. /Экспериментальное определение, обобщение и прогнозирование характеристик молекулярного переноса массы в газах: Автореф. дис.. д-ра техн. Наук. М., 1989. — 30с.
  90. А.В. Термодинамическая химия парообразного состояния. М.: Химия. 1970.-238 с.
  91. В.И., Миначев Х. М. Молекулярная хроматография. М.Л.: Наука. 1964. 138 с.
  92. С.И., Баранаев М. К. К вопросу о скорости испарения жидкостей // ЖФХ 1935.- Том VI, вып. 9. — С. 1180−1191.
  93. М.К. Кинетика испарения // Успехи химии -1938.-Том VII, № 8,. -С. 12 311 259.891. DEFINTI-N
  94. DECLARE FUNCTION SUM (KD, XF!()) DECLARE SUB MINF (NX, KD, XF (), ST!, EP, FMAX!) DECLARE SUB FPC (II, hh, XF ()) DECLARE FUNCTION GODNUP% () DECLARE SUB MENU (IM%, MEN%, TREG$) DECLARE FUNCTION GOYN% ()
  95. DECLARE SUB BBODKOR (NNN, AAA (), TEXTCOL$(), NKOL, NSTR, HAG, ZAGOL$) DECLARE SUB SLU (ISLU%, NSLU%, ASLU (), BSLU (), XSLU ()) CONST MAXNK = 20 DIM SHARED NEXP, DET, SH (MAXNK), SHH (MAXNK)
  96. KEY (15) ON KEY (16) ON KEY (17) ON KEY (18) ON90
  97. FOR I = 1 TO NEXP: P (1,1) = C (1,1) * AM (1) * .001: P (2,1) = C (2,1) * AM (2) * .001:
  98. P (3,1) = C (3,1) * AM (3) * .001: NEXT
  99. BEG: CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25:1. COLOR 15, 1
  100. FOR I = 1 TO NEXP: DAA (1,1) = TC (I): DAA (I, 2) = h (l): DAA (I, 3) = dh (I): DAA (I, 4) = t (l):
  101. DAA (I, 5) = PL (I): DAA (I, 6) = P (l, I): DAA (I, 7) = P (2,1): DAA (I, 8) = P (3,1):1. NEXT1. BEGIN: LNFN = LEN (DDl)
  102. CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25: COLOR 15,1
  103. CATE 2,4: PRINT «Расчёт коэффициентов диффузии в трёхкомпонентной смеси путём аппроксимации» LOCATE3,1: PRINT «„Ввод экспериментальных данных“
  104. CALL MENU (1, IREG1, TREG$)
  105. ON IREG1 GOTO 1000, 2000, 3000, 4000, 5000 1000: CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25: CLS: COLOR 15, 1 LOCATE 2, 24: PRINT „Назначение и краткая теория“ LOCATE 4, 25:
  106. PRINT „D11,D12,D21,D22 для средних температур Т по данным с близкими температурами.“ PRINT „Расчёт коэффициентов диффузии производится путём решения дифференциальных уравнений,“
  107. PRINT „связывающих потоки масс первых двух компонент J1, J2 с градиента -“ PRINT „ми их массовых плотностей С1, С2 (при C0=Cl+C2+C3=const) относительно. D11-D22″
  108. PRINT“ (J 1 * Со-J1 * С1 J2 * С1)/Co=-D 11 *dC 1 /dx-D 12 *dC2/dx»
  109. PRINT" (J2 * Co-J2* C2-J2* C2)/Co=-D21 * dC 1 / dx-D22 * dC2/dx"91
  110. PRINT «где Jl=dh*PL*AN (l)/t, J2=dh*PL*N2/t. При решении используются граничные условия.»
  111. PRINT «С (1,0)С10, С (2,0)=С20, C (l, h)=0, C (2,h)=0, C0=C (3,h)»
  112. PRINT" (Jl*Co-Jl*Cl-J2*Cl)/Co=-Dll*dCl/dx-D12*dC2/dx"
  113. PRINT" (J2 * Co-J2* C2-J2* C2)/Co=-D21 * dC 1 / dx-D22* dC2/dx"
  114. Для окончания ввода или просмотра данных нажмите ."
  115. PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT PRINT N1."
  116. PRINT «При необходимости результаты расчётов можно напечатать на принтере.» COLOR 4, 7: LOCATE 20, 1: PRINT «РАБОТА с программой производится с помощью меню и функциональных клавиш.»
  117. CATE 21, 1: PRINT" ИНФОРМАЦИЯ о возможных ситуациях и действиях сообщается по ходу работы."
  118. CATE 22,1: PRINT «Для аварийного выхода в меню можно нажать Ctrl+S+Enter или Esc+Enter.»
  119. COLOR 15, 1: ON GODNUP% GOTO 3000,1000 3000:CLS 0: CLS 1: SCREEN 9,1: WIDTH 80,25: COLOR 15, 1 LOCATE 2,29: PRINT «РАБОТА С БАЗОЙ ДАННЫХ «LOCATE 3,18: PRINT «Выберите режим работы, используя. и (Enter) «
  120. РАСЧЁТ КОЭФФИЦИЕНТОВ ДИФФУЗИИ»
  121. PRINT: CALL MENU (2, IREG1, TREG$)
  122. ON IREG1 GOTO 3100, 3200, BEGIN 3100: CLS 0: CLS 1: SCREEN 9, 1: WIDTH 80,25: CLS: COLOR 15, 1 LOCATE 1, 24: PRINT «ПЕРВИЧНЫЙ ВВОД ДАННЫХ»
  123. CATE 3, 7: INPUT «Введите имя файла Ваших данных (до пяти символов согласно DOS)" — NFILES
  124. CATE 5, 7: INPUT «Введите число экспериментов" — NEXP NFILES = LEFT$(NFILE$, 5) OPEN «R», #1, NDIR$ + NFILES + «.TAP», LNFN: PUT#1, 1, DD (1) CLOSE #1: KEY (10) ON1. ON KEY (10) GOSUB KONBB31
  125. CALL BBODKOR (l, DAA (), TEXTCOL$(), 8, NEXP, 10, «Ввод результатов экспериментов»)1. GOTO KONB311. KONBB31: KEY (10) OFF1. GOTO KONB31:1. RETURN
  126. KONB31: OPEN «R», #1, NDIR$ + NFILES + «.TAP», LNFN: FOR 1=1 TO NEXP
  127. CATE 5, 2: PRINT «В памяти- NOIPS- имеются следующие файлы типа… ТАР»:1. PRINT1. FILES NDIRS + «…ТАР"3210: LOCATE 12,2: INPUT -» Введите имя файла, который хотите просмотреть (первые символы до точки)" — NFILES NFILES = LEFT$(NFILE$, 5)
  128. OPEN «R», #1, NDIRS + NFILES + «.TAP», LNFN: NEXP = LOF (l) / LNFN FOR 1=1 TO NEXP GET #1,1, DD (I)
  129. DAA (I, 1) = DD (I).TC: DAA (I, 2) = DD (I).h: DAA (I, 3) = DD (I).dh: DAA (I, 4) = DD (I).t: DAA (I, 5) = DD (I).PL: DAA (I, 6) = DD (I).P1: DAA (I, 7) = DD (I).P2: DAA (I, 8) = DD (I).P3: NEXT CLOSE #1: 3220: KEY (10) ON ON KEY (10) GOSUB KONBB33
  130. CALL BBODKOR (l, DAA (), TEXTCOL$(), 8, NEXP, 10,» ВВОД РЕЗУЛЬТАТОВ ЭКСПЕРИМЕНТОВ"93
  131. GOTO KONB33 KONBB33: KEY (IO) OFF GOTO KONB33: RETURN
  132. KONB33: OPEN «R», #1, NDIR$ + NFILE$ + «.TAP», LNFN: FOR 1=1 TO NEXP
  133. CATE 1, 19: PRINT «ПРОВЕДЕНИЕ РАСЧЁТОВ СИСТЕМЫ ПО ДАННЫМ В ФАЙЛАХ» LOCATE 3,2: PRINT «Вы можете проводить проводить расчёты по системам загруженным ранее.»
  134. CATE 5, 2: PRINT «В памяти «- NDIR$- «имеются следующие файлы типа… ТАР»: PRINT1. FILES NDIR$ + «…ТАР»
  135. CATE 12, 2: INPUT — «Введите имя файла, по которому будете считать (первые символы до точки) NFILES NFILE$ = LEFT$(NFILE$, 5)
  136. OPEN «R», #1, NDIRS + NFILE$ + «.TAP», LNFN: NEXP = LOF (l) / LNFN FOR 1= 1 TO NEXP GET #1,1, DD (I)
  137. DAA (I, 1) = DD (I).TC: DAA (I, 2) = DD (I).h: DAA (I, 3) = DD (I).dh: DAA (1,4) = DD (I).t: DAA (I, 5) = DD (I).PL: DAA (I, 6) = DD (I).P1: DAA (I, 7) = DD (I).P2: DAA (I, 8) = DD (I).P3: NEXT CLOSE #1:
  138. FOR I = 1 TO NEXP: TC (I) = DAA (I, 1): h (I) = DAA (I, 2): dh (I) = DAA (I, 3): t (I) = DAA (I, 4):
  139. PL (I) = DAA (I, 5): P (l, I) = DAA (I, 6): P (2,1) = DAA (I, 7): P (3,1) = DAA (I, 8):1. NEXT4001: CLS 0: CLS 1: SCREEN 9,1: WIDTH 80,25: CLS: COLOR 15,1
  140. CATE 1,19: PRINT «ПРОВЕДЕНИЕ РАСЧЁТОВ СИСТЕМЫ ПО ДАННЫМ ФАЙЛА «+ NFILES
  141. PRINT: LOCATE 4, 2: INPUT «ВВЕДИТЕ во сколько раз начальный шаг меньше, чем коэффициент (4, 3,2)
  142. CATE 5, 12: INPUT «ВВЕДИТЕ N1 процент содержания первого компонента «- N1, N2= 100-N1.
  143. FOR I = 1 TO NEXP -1: CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25: CLS: COLOR 15,1
  144. BSLU (l) = (PJ (1,1) * POB (I) (PJ (1,1) + PJ (2,1)) * P (l, I)) / POB (I)
  145. BSLU (2) = (PJ (1,1 + 1) * POB (I + 1) (PJ (1,1 + 1) + PJ (2,1 + 1)) * P (l, I + 1)) / POB (I + 1)
  146. ASLU (1, 1) = -dcx (l, I): ASLU (1,2) = -dcx (2,1):
  147. ASLU (2, 1) = -dcx (l, 1+1): ASLU (2, 2) = -dcx (2,1 + 1)
  148. CALL SLU (ISLU, 2, ASLUQ, BSLUQ, XSLUQ)94
  149. D11(1) = XSLU (l): D12(I) = XSLU (2):
  150. PRINT» Погрешность: (C (l, 0)-C10)/C10, (C (2,0)-C20)/C20 для двух температур «CALL FPC (I, 0, XFO)
  151. OHHl (I) = (FP (1) P (l, I)) / P (l, I): OHH2(I) = (FP (2) — P (2,1)) / P (2,1) CALL FPC (I+ 1,0, XF ())
  152. DH11(1) = XF (1): DH12(I) = XF (2): DH21(I) = XF (3): DH22(I) = XF (4) PRINT «С учётом нелинейности: СКО, ШАГ= «- FMAX- ST PRINT DH11(I) — DH12(I) — DH21(I) — DH22(I)
  153. PRINT» Погрешности (C (l, 0)-C10)/C10, (C2,0 C20)/c20 для двух температур» CALL FPC (I, 0, XF ())
  154. OHl (I) = (FP (1) P (l, I)) / P (l, I): OH2(I) = (FP (2) — P (2,1)) / P (2,1) CALL FPC (I+ 1,0, XF ())
  155. PRINT» Графики для"-: COLOR IC1, 0: PRINT «CI"-: COLOR IC2,0: PRINT «
  156. C2 «-: COLOR IC3, 0: PRINT «C3 «
  157. XM = (h (I) + h (I + 1)) / 2: YM = (P (3,1) + P (3,1 + 1)) * 1.2 / 2:
  158. VIEW (30, 50)-(600, 300), 2,11. WINDOW (0, YM)-(XM, 0)
  159. NE (0, 0)-(0, XM), 0, &HAAAA
  160. NE (0, 0)-(0, YM), 0, &HAAAA RR1. XM/100
  161. CIRCLE (0, P (l, I)), RR, IC1 CIRCLE (0, P (l, I + 1)), RR, IC1 CIRCLE (0, P (2,1)), RR, IC2 CIRCLE (0, P (2,1 + 1)), RR, IC2951. NX= 1000: HX = XM/NX
  162. CIRCLE (NX * HX, P (3,1)), RR, IC3
  163. CIRCLE (NX * HX, P (3,1 + 1)), RR, IC31. FOR KK = 0 TO NX1. X = KK * HX:1. CALL FPC (I, X, XFO)
  164. FP (3) = POB (I) (FP (1) + FP (2))
  165. PSET (X, FP (1)), IC1: PSET (X, FP (2)), IC2: PSET (X, FP (3)), IC3: NEXT KK
  166. CATE 5, 6: PRINT «T=" — TC (I) — TC (I +1) LOCATE 33, 63: PRINT «X^" — h (I) — М (см)» LOCATE 34, 63: PRINT «X=" — h (I + 1) — «(см)»
  167. CLS 0: CLS 1: SCREEN 9,1: WIDTH 80, 25: CLS: COLOR 15, 1
  168. CATE 2, 12: PRINT «Матрица коэффициентов многокомпонентной диффузии (мЛ2/с)» PRINT «ДЛЯ СРЕДНИХ ТЕМПЕРАТУР «
  169. PRINT» Т Dll D12 D21 D22 СКО» FOR 1=1 ТО NEXP -1:
  170. PRINT USING» тм ##.#М#аааа ##.####лллл ##.М##аааа ##.####лллл ##.##ллллм-
  171. TC (I) + TC (I + 1)) / 2- DH11(1) — DH12(I) — DH21(I) — DH22(I)'- SH (I) NEXT
  172. PRINT «Ошибки % «FORI= 1 TO NEXP-1:
  173. PRINT USING «##.## Ш. Ш#АААА ##.####ЛЛАЛ ##М##ЛЛАА ##.Ж#ааАа (TC (I) +
  174. TC (I + 1)) / 2- OHl (I) * 100- OH2(I) * 100- OH3(I) * 100- OH4(I) * 1001. NEXT
  175. CATE 25, 15: PRINT «Вы будете печатать результаты? (Y/N)"-1. ON GOYN% GOTO PR4,4003:
  176. PR4: CLS: LS = 1: LC = 10: NC 8: PRINT1. PRINT «N1=" — AN (1)-«%»
  177. PRINT ««- TEXTCOL$(l) — TEXTCOL$(2) — TEXTCOL$(3) — TEXTCOL$(4) — TEXTCOL$(5) — TEXTCOL$(6) — TEXTCOL$(7) — TEXTCOL$(8) FOR 1ST = 1 TO NEXP:
  178. PRINT USING «##.##AAAA" — DAA (IST, 1) — DAA (IST, 2) — DAA (IST, 3) — DAA (IST, 4) — DAA (IST, 5.- DAA (IST, 6) — DAA (IST, 7) — DAA (IST, 8):1. NEXT
  179. PRINT:PRINT: PRINT» T Dll D12 D21 D22(mA2/c)» FOR 1=1 TO NEXP -1:
  180. PRINT USING «##M ##.#mAAAA №.##№AAAA ##M##AAAA ##M##AAAA «- (TC (I) +
  181. TC (I + 1)) / 2- DH11(1) — DH12(I) — DH21(I) — DH22(I)1. NEXT
  182. RINT:LPRINT:LPRINT «N1=" — AN (1)-«%»
  183. RINT ««- TEXTCOLS (l) — TEXTCOL$(2) — TEXTCOL$(3) — TEXTCOL$(4) — TEXTCOL$(5) — TEXTCOL$(6) — TEXTCOL$(7) — TEXTCOL$(8)961. FOR 1ST = 1 TO NEXP:
  184. RINT USING «##.##AAAA" — DAA (IST, 1) — DAA (IST, 2) — DAA (IST, 3) — DAA (IST, 4) —
  185. DAA (IST, 5) — DAA (IST, 6) — DAA (IST, 7) — DAA (IST, 8):1. NEXT
  186. RINT: LPRINT» T Dll D12 D21 D22(mA2/c)» FOR I = 1 TO NEXP-1:
  187. RINT USING «##.## ##.####ААЛА ##.####ЛАЛА ##.####AAAA ##.####AAAA (TC (I) +
  188. TC (I + 1)) / 2- DH11(1) — DH12(I) — DH21(I) — DH22(I)1. NEXT1. RINT CHR$(12)4003: LOCATE 25, 14: PRINT «Вы будете менять N1? (Y/N) «-
  189. ON GOYN% GOTO 4001, BEGIN:5000: CLS 0: CLS 1: LOCATE 2, 18: PRINT ««
  190. CATE 5, 24: PRINT «Конец работы до свидания «1. CLOSE1. END1. KONES:
  191. KEY (15) OFF: KEY (16) OFF: KEY (17) OFF: KEY (18) OFF: KEY (19) OFF: KEY (20) OFF: CLOSE: GOTO BEG RETURN
  192. SUB BBODKOR (NBB, AAA (), TEXTCOL$(), NCOL, NSTR, HAG, ZAGOL$) LSTR = 1: IF NBB = 3 THEN LSTR = 5 LCOL = HAG
  193. DIM TEXTSTR (4, 20) AS STRING * 24' LSTR 'DIM TEXTCOL$(NBB, NSTR) AS STRING * LSTR FOR I = 1 TO NSTR: TEXTSTR (1,1) = STR$(I): NEXT» SCREEN 0: COLOR 15, 1 ISTR = 1: ICOL = 1
  194. BEGBB: CLS: LOCATE 1, 23: PRINT ZAGOL$
  195. CATE 2, 1: PRINT «используя <стрелки> числа , — Введите или просмотрите данные в произволь-«
  196. CATE 3,1: PRINT «ном порядке. Для окончания работы и записи данных НАЖМИТЕ ."1. FOR ICO = 1 ТО NCOL:
  197. CATE 6, LSTR + (ICO -1) * LCOL: PRINT «TEXTCOLS (ICO): NEXT
  198. FOR ICO = 1 TO NCOL: FOR 1ST =1 TO NSTR:
  199. CATE 7 1 + 1ST, LSTR + (ICO -1) * LCOL: PRINT USING «##.### AAAA" — AAA (IST, ICO): NEXT NEXT
  200. BB4: LOCATE 4, 1, 0, 4, 7: PRINT SPACE$(78)
  201. CATE 4, 20: PRINT «Полное значение текущей переменной" — AAA (ISTR, ICOL) LOCATE 7 -1 + ISTR, LSTR + (ICOL -1) * LCOL, 1, 5, 8: LOCATE, 1 BB6: C$ = INKEY$: IF C$ = «» THEN GOTO BB6 IF LEN (C$) = 2 THEN C$ = RIGHT$(C$, 1) С = ASC (C$)
  202. С = 72 AND ISTR > 1 THEN ISTR = ISTR 1: GOTO BB4
  203. CATE 4, 1: PRINT SPACE$(78):LOCATE 4, 1, 1, 5, 8: LOCATE, 1 INPUT «Введите значение переменной «- AAA (ISTR, ICOL): GOTO BEGBB: END IF GOTOBB6 ENDSUB
  204. SUB FPC (II, XX, XF ()) hh = h (II) * .01 :X = XX * .01
  205. F1=(B2* A12-B1 * A22)/A: F2 = (B1 * A21 -B2* All)/A DET = (A11 + A22) A 2 / 4 + (A12 * A21 A11 * A22) SELECT CASE DET CASE IS < 0
  206. AL1 = (A11 + A22) / 2: AL2 = SQR (-DET)
  207. BSLU (l) = -F1 / EXP (AL1 * hh): BSLU (2) = -F2 * A12 / EXP (AL1 * hh)
  208. ASLU (1, 1) = COS (AL2 * hh): ASLU (1, 2) = SIN (AL2 * hh)
  209. ASLU (2,1) = (AL1 A11) * COS (AL2 * hh) — AL2 * SIN (AL2 * hh)
  210. ASLU (2,2) = AL2 * COS (AL2 * hh) + (AL1 All)* SIN (AL2 * hh)
  211. CALL SLU (ISLU, 2, ASLU (), BSLU (), XSLU ())
  212. CK1 = XSLU (l): CK2 = XSLU (2)
  213. CA1 = (CK1 * (AL1 A11) + CK2 * AL2) / A12:
  214. CA2 = (CK2 * (AL1 All) — CK1 * AL2) / A12:
  215. FP (1) CK1 * EXP (AL1 * X) * COS (AL2 * X) + CK2 * EXP (AL1 * X) * SIN (AL2 FP (2) = CA1 * EXP (AL1 * X) * COS (AL2 * X) + CA2 * EXP (AL1 * X) * SIN (AL2 CASE IS > 0
  216. AL1 = (A11 + A22) / 2 + SQR (DET): AL2 = (A11 + A22) / 2 SQR (DET): CA1 = (AL1 — A11) / A12: CA2 = (AL2 — A11) / A12: BSLU (l) = -F1: BSLU (2) = -F2
  217. ASLU (1,1) = EXP (AL1 * hh): ASLU (1, 2) = EXP (AL2 * hh)
  218. ASLU (2, 1) CA1 * EXP (AL1 * hh): ASLU (2, 2) = CA2 * EXP (AL2 * hh)
  219. CALL SLU (ISLU, 2, ASLUO, BSLU (), XSLU ())
  220. CK1 = XSLU (l): CK2 = XSLU (2)
  221. FP (1) = CK1 * EXP (AL1 * X) + CK2 * EXP (AL2 * X) + F1
  222. FP (2) = CK1 * CA1 * EXP (AL1 * X) + CK2 * CA2 * EXP (AL2 * X) + F2
  223. PRINT «1" — CK1- CA1- AL1- Fl- «kl»
  224. PRINT «2м- CK2- CA2- AL2- F2- «k2"1. CASE 01. AL = (A11 + A22) / 2
  225. BSLU (l) = -Fl / EXP (AL * hh): BSLU (2) = -F2 / EXP (AL * hh) * A12
  226. ASLU (1, 1) = 1: ASLU (1, 2) = hh
  227. ASLU (2, 1) = AL A11: ASLU (2, 2) = 1 + (AL — A11) * hh
  228. CALL SLU (ISLU, 2, ASLU (), BSLU (), XSLU ())
  229. CK1 = XSLU (l): CK2 = XSLU (2)
  230. CA1 = (CK1 * (AL A11) + CK2) / A12: CA2 = CK2 * (AL — A11) / A12:
  231. FP (1) = (CK1 + CK2 * X) * EXP (AL * X) + Fl
  232. FP (2) = (CA1 + CA2 * X) * EXP (AL * X) + F21. END SELECT1. END SUB
  233. FUNCTION GODNUP% LOCATE 25, 17:
  234. PRINT «Листать далыне-PgDn, листать назад -PgUP" — FL: C$ = INKEY$: IF C$ = «» THEN GOTO FL С = ASC (C$)
  235. LEN (C$) = 1 THEN GOTO FL C$ = RIGHT$(C$, 1): С = ASC (C$) IF С = 81 THEN GODNUP% 1: EXIT FUNCTION IF С = 73 THEN GODNUP% = 2: EXIT FUNCTION GOTO FL END FUNCTION1. FUNCTION GOYN%1. STATIC C$ YN: C$ = INKEY$
  236. C$ = «Y» OR C$ = «y» OR C$ = «K» OR C$ ="-» THEN GOYN% = 1
  237. ELSEIF C$ = «N» OR C$ = «n» OR C$ = «'» OR C$ = «в» THEN GOYN% = 2 ELSE GOTO YN END IF END FUNCTION
  238. CATE 7, 13: PRINT «Выберите в меню режим работы, используя CHR$(24)-««- CHR$(25)-» Ё Enter)»: PRINT FOR I = 1 TO N (IM%): LOCATE NSTR + 2 * 1,10: PRINT I-: NEXT FOR I = 1 TO N (IM%): LOCATE NSTR + 2 * 1,14: PRINT TM (IM%, I) 1=1:
  239. LOCATE NSTR + 2 * I, 8, 1: PRINT «#» 6 C$ = INKEY$: IF C$ = «» THEN 6 IF LEN (C$) = 2 THEN C$ = RIGHT$(C$, 1) С = ASC (C$)
  240. С = 13 THEN MEN% = I: TREGS = TM (IM%, I): EXIT SUB IF С = 72 THEN 1С = 1:1 = I 1: GOTO 5 IF С = 80 THEN 1С = I: I = I + 1: GOTO 5 GOTO 699
  241. IF I = N (IM%) + 1 THEN 1=1 IF I = 0 THEN I = N (IM%) LOCATE NSTR + 2* 1С, 8, 1: PRINT SPACE$(1) GOTO 4 STOP END SUB
  242. SUB SLU (ISLU, NSLU, ASLU (), BSLU (), XSLU ()) SELECT CASE ISLU
  243. CASE 1 FOR I = 1 TO NSLU -1 FOR J = I + 1- TO NSLU ASLU (J, I) = -ASLU (J, I) / ASLU (I, I) FOR К = I + 1 TO NSLU:
  244. ASLU (J, K) = ASLU (J, K) + ASLU (J, I) * ASLU (I, K) NEXT К
  245. BSLU (J) = BSLU (J) + ASLU (J, I) * BSLU (I) NEXT J NEXT I
  246. XSLU (NSLU) = BSLU (NSLU) / ASLU (NSLU, NSLU) FOR I = NSLU 1 TO 1 STEP -1: HSLU = BSLU (I) FOR J = I + 1 TO NSLU: HSLU = HSLU — XSLU (J) * ASLU (I, J): NEXT J XSLU (I) = HSLU / ASLU (I, I) NEXT I CASE 2 MS = 0
  247. RS = MS * ASLU (I, 0) LS ASLU (K, 0) ASLU (K, 0) = LS * ASLU (I, 0) + MS * ASLU (K, 0) ASLU (I, 0) = RS1. NEXT К NEXT I
  248. FOR I = NSLU TO 1 STEP -1: MS = 0
  249. FOR К = 0 TO NSLU -1 1: MS = MS + ASLU (0, NSLU — K) * ASLU (I, NSLU — K) NEXT К
  250. ASLU (0,1) = (ASLU (1,0) MS) / ASLU (I, I) NEXT I
  251. FOR К = 1 TO NSLU: XSLU (K) = ASLU (0, K): NEXT CASE 3100
  252. DSLU = ASLU (1, 1) * ASLU (2, 2) ASLU (2,1) * ASLU (1, 2) XSLU (l) = (BSLU (l) * ASLU (2, 2) — BSLU (2) * ASLU (1, 2)) / DSLU XSLU (2) = (BSLU (2) * ASLU (1, 1) — BSLU (l) * ASLU (2, 1)) / DSLU END SELECT END SUB
  253. FUNCTION SUM (II, XFO) SU = 01. CALL FPC (II, 0, XF ())
  254. SU = SU + ((FP (1) P (l, II)) / P (l, II))л 2 + ((FP (2) — P (2, II)) / P (2, II))A 21. CALL FPC (II, h (II), XF ())
  255. FP (3) = POB (II) (FP (1) + FP (2))
  256. SU = SU + ((FP (3) P (3, II)) / P (3, II)) A 21. CALL FPC (II+1,0, XF ())
  257. SU = SU + ((FP (1) P (l, II + 1)) / P (l, II + 1)) A 2 + ((FP (2) — P (2, II + 1)) / P (2, II + 1)) A 2 CALL FPC (II + 1, h (II + 1), XF ()) FP (3) = POB (II + 1) — (FP (1) + FP (2)) SU = SU + ((FP (3) — P (3, II + 1)) / P (3, II + 1))A 2 SUM=SU/6 END FUNCTION
  258. Завлабораторией оптимизации технологичес»»» -----------««кокрасочных производств, к.т.н.
  259. Старший научный сотрудник лаборатории аппаратуттп-фртиплпгического оформления процессов синтеза лаковых смол, к.т.н
  260. Казанское ордена Октябрьской революции производственное объединение &bdquo-Органический синтез"420 051 Казань Телеграф «Синтез» телетайп № 168 с/с 9 225 403
  261. Начальник производства PI^I^S^* ¦ Недбавлюк Б.Е.
Заполнить форму текущей работой

Сессия? Спокойно!