Дипломы, курсовые, рефераты, контрольные...
Срочная помощь в учёбе

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

ДиссертацияПомощь в написанииУзнать стоимостьмоей работы

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

Содержание

  • Глава 1. ЛИТЕРАТУРНЫЙ ОБЗОР
    • 1. 1. Образование виброкипящего слоя
    • 1. 2. Обзор конструкций вибрационных смесителей
    • 1. 3. Поведение систем под действием вибрации
    • 1. 4. Выводы и формулировка задачи исследования
  • Глава 2. МОДЕЛИРОВАНИЕ ДВИЖЕНИЯ СЫПУЧЕГО МАТЕРИАЛА ПОД ДЕЙСТВИЕМ ВИБРАЦИИ
    • 2. 1. Движение одиночной частицы на вибрирующей поверхности
    • 2. 2. Циркуляция сыпучего материала под действием вертикальных колебаний в аппаратах с вертикальными стенками
    • 2. 3. Циркуляция сыпучего материала под действием вертикальных колебаний в аппаратах с наклонной боковой границей
    • 2. 4. Гидродинамическая модель виброожиженного слоя с позиций механики сплошной среды
  • Глава 3. ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ ДВИЖЕНИЯ СЫПУЧЕГО МАТЕРИАЛА ПОД ДЕЙСТВИЕМ ВИБРАЦИИ
    • 3. 1. Разработка экспериментальной установки и методика проведения экспериментов
    • 3. 2. Экспериментальное исследование движения сыпучего материала под действием вертикальной вибрации
    • 3. 3. Сравнительный анализ параметров движения частиц сыпучего материала по результатам эксперимента и гидродинамической модели
    • 3. 4. Экспериментальное исследование движения сыпучего материала под действием не вертикальной вибрации
  • Глава 4. МОДЕЛИРОВАНИЕ ПРОЦЕССА СМЕШИВАНИЯ СЫПУЧИХ МАТЕРИАЛОВ И РАЗРАБОТКА НОВЫХ КОНСТРУКЦИЙ ВИБРАЦИОННЫХ СМЕСИТЕЛЕЙ
    • 4. 1. Ячеечная модель процесса смешивания
    • 4. 2. Идентификация параметров модели
    • 4. 3. Проверка адекватности модели
    • 4. 4. Анализ эффективности вибрационного смесителя
    • 4. 5. Новые конструкции вибрационных смесителей

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

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

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

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

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

Изучению вибрационного воздействия на различные химико-технологические процессы посвящено значительное количество работ. Из них особое значение имеют работы В. А. Членова, Н. В. Михайлова, В. Д. Варсанофьева, Э.Э. Кольмана-Иванова, И.И. Блехма-на и др. связанные с изучением воздействия вибрации на обработку сыпучих материалов [1−5].

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

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

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

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

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

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

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

Практическая ценность результатов работы состоит в следующем:

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

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

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

4. Новая конструкция вибрационного барабанного смесителя и методика расчета его основных режимных и reo8 метрических параметров используются при разработке специализированного комплексного технологического оборудования для промышленного выпуска многокомпонентных сухих готовых товарных форм специального назначения из природного сырья (приложение 1).

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

ОСНОВНЫЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ РАБОТЫ.

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

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

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

4. Создана оригинальная экспериментальная установка с использованием видеокамеры «PANASONIC NV-R100», позволяющая на основе достоверной информации подтвердить результаты расчета параметров движения частиц сыпучего материала.

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

6. Разработан новый способ определения неоднородности смеси сыпучих компонентов.

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

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

  1. В.А., Михайлов Н. В. Новый принцип создания «кипящего слоя». -Докл. АН, 1964, т. 154, № 3.
  2. В.А., Михайлов Н. В. Виброкипящий слой. М.: Наука, 1972. — 340 с.
  3. В.Д., Кольман-Иванов Э.Э. Вибрационная техника в химической промышленности. М.: Химия, 1985.- 240 с.
  4. И.И., Джанелидзе Г. Ю. Вибрационное перемешивание. М., 1964.
  5. И.И. Что может вибрация? М.: Наука, 1988. — 208 с.
  6. H.H. Вибрационный смеситель для тонкодисперсных материалов. В кн.: Современные машины и аппараты химических производств: Докл. II Всесоюз. науч. конф. -Чимкент, 1980, с. 672−676.
  7. H.H. Моделирование и интенсификация процесса вибрационного смешивания: Автореф. дис. канд. тех. наук.- Иваново, 1980. 17 с.
  8. Ю.И. Аппараты для смешения сыпучих материалов. М.: Машиностроение, 1973. — 215 с.
  9. В.Д. Техника и применение вибрирующего слоя.- Киев: Наук, думка, 1977. 239 с.
  10. A.c. № 655 419 СССР. Вибрационный смеситель. Б.и.№ 13,1979.
  11. A.c. № 1 558 449 СССР. Вибрационный смеситель. Б.и.№ 15,1990.
  12. Патент РФ № 2 035 986. Вибрационный смеситель.
  13. A.c. № 1 499 831 СССР. Вибрационный смеситель.
  14. Патент РФ № 2 137 536. Вибрационный смеситель.
  15. Патент РФ № 2 122 891. Вибрационный смеситель.
  16. И.И., Гортинский В. В., Птушкина Г. Е. Движение частицы в колеблющейся среде при наличии сопротивления типа сухого трения. -Изв. АН СССР. Механика и машиностроение, 1963, № 4, с. 27−30.
  17. A.C., Баскаков А. П., Грачев C.B. Химико-термическая обработка в кипящем слое. М.: Машиностроение, 1985. — 158 с.
  18. Э.М. Кальцинация технического бикарбоната натрия в вибрирующем слое. ЖПХ, 1960, т. ЗЗ, № 6.
  19. В.А., Михайлов Н. В. Некоторые свойства вибро-кипящего слоя. ИФЖ, 1965, т.9, № 2.
  20. Kroll W. Forschung auf dem Gebiete des Ingenieurwesens, 1954, 20,1.
  21. Р.Ф., Украинский Л. Е. О динамике твердых частиц взвешенных в несжимаемой жидкости при вибрационных воздействиях. -Изв. АН СССР. Механика твердого тела, 1975, № 5, с. 31−40.
  22. Р.Ф., Цапенко A.C. О динамике газовых пузырьков в жидкости, подверженной вибрационным воздействиям. -В кн.: Вопросы математической физики и теории колебаний. Ивановский энергетический институт им. В. И. Ленина, 1975, вып. З, с. 5−12.
  23. Р.Ф., Лакиза В. Д., Цапенко A.C. О явлениях вибрационного перемещения и образования периодических структур в условиях, близких к невесомости. Изв. АН СССР. Механика твердого тела, 1977, № 2, с. 56−59.
  24. Вибрации в технике: Справочник. В 6-ти т./Ред. совет: В. Н. Челомей (пред.). М.: Машиностроение, 1981. — Т.1 Колебания линейных систем / Под ред. В. В. Болотина, 1999. — 504 с.
  25. Вибрации в технике: Справочник. В 6-ти т./Ред. совет: В. Н. Челомей (пред.). М.: Машиностроение, 1981. — Т.2 Колебания нелинейных механических систем / Под ред. И. И. Блехмана, 1979. — 351 с.
  26. Действие вибрации на нелинейные механические систе-мы/А.А. Пасько// Сб. науч. тр. ТГТУ. 2000, Вып.б.-с. 8689.
  27. В.Г., Баскаков А. П., Берг Б. В. и др. Псевдоожижение. М.: Химия, 1991. — 400 с.
  28. D. J. Tritton, Physical Fluid Dynamics (Oxford, New York, (1988), p. 63.
  29. J. B. Knight, H. M. Jaeger, and S. Nagel, «Vibration-Induced Size Separation in Granular Media: The Convection Connection», Physical Review Letters, 70, 3728 (1993).
  30. E. E. Ehrichs, H. M. Jaeger, G. S. Karczmar, J. B. Knight, V. Y. Kuperman, and S. R. Nagel, «Granular Convection Observed by Magnetic Resonance Imaging,» Science, 267, 1632 (1995).
  31. J. B. Knight, E. E. Ehrichs, Y. Y. Kuperman, J. K. Flint, H. M. Jaeger, and S. R. Nagel, «Experimental Study of Granular Convection,» Phys. Rev. E, 54, 5726 (1996).
  32. M. Faraday, Philos. Trans. R. Soc. London 52, 299 (1831).
  33. M. Bourzutschky and J. Miller, Phys. Rev. Lett. 74, 2216 (1995).
  34. H. Takahashi, A. Suzuki, and T. Tanaka, Powder Technol. 2,65 (1968/69).
  35. Y. Kuperman, E. E. Ehrichs, H. M. Jaeger, and G. S. Karczmar, «A New Technique for Differentiating Between Diffusion and Flow in Granular Media using Magnetic Resonance Imaging,» Review of Scientific Instruments, 66, 4355 (1995).
  36. J.B. Knight, «External Boundaries and internal shear bands in granular convection», Phys. Rev. E, 55, 6016 (1997).
  37. L. Grossman, «The Effects of Container Geometry on Granular Convection», Phys. Rev. E 56, 3290(1997).
  38. S. Douady, S. Fauve, and C. Laroche, Europhys. Lett. 8, 621 (1989).
  39. JT.И. Механика сплошной среды. Т.1. М.: Наука, 1983. — 528 с.
  40. Л.И. Механика сплошной среды. Т.2. М.: Наука, 1983. — 560 с.
  41. Л.Г. Механика жидкости и газа. М.: Наука, 1973. — 847 с.
  42. Р., Стьюарт В., Лайтфут Е. Явления переноса. М.: Химия, 1974. — 687 с.
  43. А.А. Пасько, В. П. Таров, Т. В. Игнатьева Стенд для исследования механики сыпучих материалов при вибрационном воздействии // 5-я науч. конф. ТГТУ. Тамбов, 2000.
  44. Т. Poschel, D. Rosenkranz Experimental Study of Horisontally Shaken Granular Matter The Swelling Effect.
  45. C. Salueta, T. Poschel: Convection in horizontally shaken granular material. Europ. J. Phys. 2, 1999.
  46. A.A. Исчисление вероятностей. M.: ГИЗ, 1924. -202 с.
  47. Баруча-Рид А. Т. Элементы теории марковских процессов и их приложения. М.: Наука, 1969. — 225 с.
  48. . В. Введение в теорию вероятностей и ее приложения. М.: Мир, 1967. — 620 с.
  49. А.А., Першин В. Ф., Таров В. П., Негров В. Л. Математическое моделирование процесса смешения сыпучих материалов в вибрационном смесителе/ Вестник ТГТУ, № 2, 2000.
  50. В.Ф. Модель процесса смешивания сыпучего материала в поперечном сечении гладкого вращающегося барабана// Теор. основы хим. технологии, 1989, T. XXIII, № 3. с. 370−377.
  51. В.Ф. Методы расчета и новые конструкции машин барабанного типа для переработки сыпучих материалов: Дис. д-ра техн. наук. Тамбов, 1994. — 431 с.
  52. В.Ф. Машины барабанного типа: основы теории, расчета и конструирования. Воронеж: Изд-во ВГУ, 1990. — 168 с.
  53. Determination of mixture inclination to segregation/ V. Per-shin, S. Barishnikova, A. Pasko, Y. Selivanov //Abstracts of Papers World Congress on Particle Technology 3, Brighton, UK, 1998, P.173.
  54. У. Pershin, M. Sviridov, A. Pasko, A. Sherbakov, E. Mandrika Stochastie-determinate and determinate-stochastie mixing// 13th International Congress of Chemical and Process Engineering,
  55. Praha, Czeh Republic, 1998, V.7,P.177.
  56. Ю.И.Гусев, И. Н. Карасев, Э.Э.Кольман-Иванов и др. Конструирование и расчет машин химически производств. -М.: Машиностроение, 1985. 408 с.
  57. А.с. СССР № 1 769 084. Способ определения неоднородности смеси сыпучих компонентов. G 01N15/02, Пасько А. А. и др.
  58. Ю.Т., Казаков М. Т., Першин В. Ф., Ткачев А. Г. Исследование процесса приготовления многокомпонентных смесей. В сб.: Тезисы докладов I научной конференции ТГТУ, Тамбов, 1994, с. 81.
  59. В.Ф., Негров В Л., Селиванов Ю. Т. Исследование процесса смешивания полидисперсных материалов // Роль118мол. конструкторов и исследователей хим. машиностроения. Тезисы докладов V Всесоюзной науч.-техн. конференции. -Северодонецк, 1986. с. 28.
  60. В.Ф., Негров B.JI., Селиванов Ю. Т. Моделирование процесса смешивания полидисперсных материалов. В сб.: Тезисы докладов Всесоюзной конференции «Технология сыпучих материалов ХИМТЕХНИКА-86», Белгород, 1986, с. 49−50.
  61. Об использовании результатов диссертационной работы Пасько A.A. «Разработка новых конструкций вибрационных смесителей барабанного типа для сыпучих материалов и методика их расчёта».
  62. Анализ результатов, полученных соискателем и наши опытно-промышленные данные при использовании рекомендаций Пасько A.A., позволяют сказать, что за счет снижения нежелательного эффекта сегрегации существенно повышается качество готового продукта.
  63. Президент научно-технической 4фирмы «ЛИОНИК» ПЧЕЛОПРОМа Минсельхозпрода России, кандидат технических наук, лауреатпремии Совета Министров СССР, ^^лауреат премии Правительства Е.А.Мандрыка
  64. Российской федерации в областинауки и техники
  65. Концентрации правого контура}
  66. Концентрации правого контура}
  67. Концентрации левого контура}
  68. Концентрации левого контура}single- single — single- single -i: integer- j: integer- DWG: Boolean-f: text-1. Tangens} {*}
  69. Function tan (x:single):extended-}1. Begin {*}tan:=sin (x)/cos (x) — {*}1. End- {*}1. Proc 1 ***}
  70. Procedure Procl (xl, yl, x2, y2:single-var k: single)-1. Begink:=(y2-yl)/(x2-xl) — End-1. End Proc1
  71. Proc2 ***} Procedure Proc2 — var i: integer-beta:single- x, y: single- ul, u2, u3:single- AA., BB, CC, DD: single- Begin NG:=1-ul:=Arctan (kl) — u2:=Pi+Arctan (k2)-u3:=Pi+Arctan (k3) — {
  72. Writeln (ul/pi*18 0:3:3) — Writeln (u2/pi*18 0:3:3) — Writeln (u3/pi*180:3: 3) -readln- }beta:=u2- For i:=l to n do Begin
  73. BB: =2*tan (beta) * (yc-tan (beta) *xc) — CC:=sqr (yc-tan (beta)*xc)-sqr® —
  74. DD:=sqr (BB)-4*AA*CC- x:=(-BB+sqrt (DD))/(2*AA) — y:=tan (beta)*(x-xc)+yc- end-if ((beta > u3) and (beta≤l.5*Pi)) then begin1. AA:=l+sqr (tan (beta)) —
  75. BB:=2*tan (beta)*(yc-tan (beta)*xc) — CC: =sqr (yc-tan (beta) *xc) -sqr ® — DD:=sqr (BB)-4*AA*CC- x:=(-BB-sqrt (DD))/(2*AA) — y:=tan (beta)*(x-xc)+yc- end-if beta > 1.5*Pi then begin1. AA:=l+sqr (tan (beta)) —
  76. BB:=2*tan (beta)*(yc-tan (beta)*xc) — CC:=sqr (yc-tan (beta)*xc)-sqr® — DD:=sqr (BB)-4*AA*CC- x: =(-BB+sqrt (DD))/(2*AA) — y:=tan (beta)*(x-xc)+yc- end-end-
  77. Wi 0. 1.. x: =x- Wi [0] [i]. y: =y-beta:=beta+2*Pi/n-if beta ≥ 2*Pi then beta := beta 2*Pi- End{For)-1. End-
  78. End Proc2 ***} {*** Begin PLATZ
  79. Function Platz (i .-integer) :single-var a, b, c, p: single- Begina:=sqrt (sqr (xc-Wi0.1.x) + sqr (yc-Wi[0][i]. y)) — b:=sqrt (sqr (xc-Wi[0][i-1]. x) + sqr (yc-Wi[0][i-1]. y)) — c:=sqrt (sqr (Wi[0][i]. x-Wi[0][i-1].x) + sqr (Wi[0][i]. y-Wi[0][i-1].y)) — p:=(a+b+c)/2-
  80. PLATZ := sqrt (p*(p-a)*(p-b)*(p-c)) — End-1. END PLATZ)
  81. Proc2NEW ***} Procedure Proc 2 NEW — var i: integer-beta:{single}extended- {!!!} x, y: single- ul, u2, u3:single- AA, BB, CC, DD: single- label 5,10- Begin NG:=1-ul:=Arctan (kl) — u2:=Pi+Arctan (k2) — u3:=Pi+Arctan (k3)-beta:=u2−1. For i:=1 to n do Begin
  82. BB:=2*tan (beta)*(yc-tan (beta)*xc) — CC:=sqr (yc-tan (beta)*xc)-sqr® — DD:=sqr (BB)-4*AA*CC- x:=(-BB+sqrt (DD))/(2*AA) — y:=tan (beta)*(x-xc)+yc- end-if ((beta > u3) and (beta≤l.5*Pi)) then begin1. AA:=l+sqr (tan (beta)) —
  83. BB:=2*tan (beta)*(yc-tan (beta)*xc) —
  84. CC:=sqr (yc-tan (beta)*xc)-sqr®-1. DD:=sqr (BB)-4*AA*CC-x: = (-BB-sqrt (DD))/(2*AA) — ¦ y:=tan (beta)*(x-xc)+yc- end-if beta > 1.5*Pi then begin1. AA:=l+sqr (tan (beta)) —
  85. BB:=2*tan (beta)*(yc-tan (beta)*xc) —
  86. CC:=sqr (yc-tan (beta)*xc)-sqr®-1. DD:=sqr (BB)-4*AA*CC-x:=(-BB+sqrt (DD))/(2*AA)-y:=tan (beta) * (x-xc) +yc- end-end-
  87. S1. :=sqrt (p*(p-a)* (p-b)*(p-c) SS:=SS+Si.- End- End-1. End Geron1. Begin NewCenter ***}
  88. Возвращате ПРАВИЛЬНЫЕ значения центра сечения } Procedure New Center- var i: integer-
  89. SO, Sx, Sy: single- Begin SO:=0- Sx:=0- Sy:=0−1. For i:=l to n do Begin
  90. SO:=SO+S1. — Sx:=Sx+(Si.*Ct[i] .x) — Sy:=Sy+(S[i]*Ct[i] .y) — End- x c:=Sx/SO-yc:=Sy/SO- End-1. End NewCenter ***}
  91. Begin Zona ***) Procedure Zona- var i, j: integer-1: single- Begin For i:=l to n do Beginfor j:=nn-l downto 1 do begin1:= (1 sqrt (j/nn))/sqrt (j/nn) — Winn-j.1.x := (Wi[0][i]. x + 1* xc)/(1+1) — Wi[nn-j][i]. у := (Wi[0][i]. у + 1* yc)/(l+l) — end-1. End- End-
  92. End Vizual {*** Begin Center
  93. Заполняет массив координатами центров зон) Procedure Center (var A: typeiXY- var В: typeCenter) — var i, j: integer- x, у: single-1. Begin For i:=1 to n dofor j:=0 to nn-2 do begin if ion then begin
  94. Bj+l.1.x := (A[j][i]-x + A[j+1][i+1]. x)/2-
  95. Вj +1. [ i ] .у := (A[j] 1. .y + A[j+1] [i + 1] -y)/2-endelsebegin
  96. Вj +1. 1. .x := (A[ j] [i] .x + A[j+1] 1. .x)/2- B[j+1] [i] .y := (A [ j ] [i] .у + A [ j +1] [1] .у)/2- end-
  97. Bnn.1.x :=(2*x+xc)/3 — B[nn][i]. y :=(2*y+yc)/3 — {Circle (getmaxx div 2 +round (k*B[j][i]. x), getmaxy div 2 round (k*Bj.1.y), 2)-} End-1. End-1. End Center1. Begin MeinColor ***}
  98. Function MeinColor (C:single):Byte-var i: byte-1. Begin1. (C<0) or (C>1) then begin sound (440)-Delay (500)-Nosound-Halt-end-for i:=0 to 63 doif (C≥i/64) and (C≤ (i + 1)/64) then MeinColor:=i+64−1. End-
  99. End MeinColor {*** Begin P1W
  100. Function P1W (a:integer):single- Begin
  101. P1W:=P01*(1-WCa.1.) — End- {*** End P1W {*** Begin P2W
  102. Procl (xx, yy, xc, yc, kl) — Procl (0, b, xc, yc, k2) — Procl (0, -R, xc, yc, k3) — Proc2- Geron-1. NewCenter- {*}1. ClearDevice-) {*}
  103. ColZona:=15- {Цвет сетки) SetColor (Col Zona)-1. Proc 1 (xx, У Y, x kl1. Proc 1(0, b, x y k21. Proc .1(0, -R, X y k31. Proc2NEW- Geron-
  104. New Center- Proc2NEW- Geron-1. Zona- {*}for i:=l to n do for j:=nn-l downto 0 do begin
  105. Mij.1.x:=-Wi[j][i]. x- Mi[j] [i]- у:= Wi[j] [i]. y- end-1. Vizual (Wi)-1. Center (Wi, Wcenter)-xc:=-xc-yc: = yc-1. Vizual (Mi)-1. Center (Hi, Mcenter)-if DWG then Close (f) —
  106. Заполним начальными концентрациями} Randomize- For i:=1 to n do for j:=1 to nn do Begin
  107. WCj.1.:=0{random) — MC[j][i]: =0{random} -1. End-
  108. SOUND (440)-DELAY (1000) END-1. HALT-1:step:=step+l- {Покажем на экране} S CW: = 0- S CM: = 0- Xv:=0 j Xw: = 0 Xm:=0
  109. For i:=l to n do for j:=1 to nn Begindo
  110. SetFillStyle (l, MeinColor (WCj. 1.)) —
  111. FloodFill (getmaxx div 2 +round (k*Wcenterj.1.x), getmaxy div 2 round (k*Wcenterj.1.y), ColZona) — SetFillStyle (l, MeinColor (MC[j][i])) —
  112. FloodFill (getmaxx div 2 +round (k*Mcenterj. 1.. x), getmaxy div 2 round (k*Mcenterj.1.y), ColZona) —
  113. Сумма концентраций} SCW:=SCW+WCj.1.- SCM:=SCM+MC[j][i]-
  114. Xv:=Xv+WCj.1.+MC[j][i]- Xw:=Xw+WC[j][i]- Xm:=Xm+MC[j][i]- End-1. SC:=SCW+SCM-
  115. Xv:=Xv/(2*n*nn) — Xw:=Xw/(n*nn) — Xm:=Хш/(n*nn) —
  116. Коэфф. НЕОДНОРОДНОСТИ} Kv:= 0−1. Kw:= 0−1. Km:=0-for i:=1 to n dofor j: =1 to nn do Begin
  117. WCnewj.1.:=WC[j][i]-WC[j][i]*P1W (j+1) + WC[j+1][i]*P2W (j) +MC[j][i]*P2W (j) WC[j][i]*P2M (j)else
  118. WCnewj.1.:=WC[j][i]-WC[j][i]*P1W (j+1) + WC[j+1][i]*P2W (j)-for j:=2 to nn-1 do
  119. WCnew j.1.:=WC[j][i]-WC[j][i]*P1W (j+1)-WC[j][i]*P2W (j-1)+ WC[j+1][i]*P2W (j)+WC[j-1][i]*P1W (j)-j:=nn-
  120. WCnewj.1.:=WC[j][i]+WC[j-l][i]*P1W (j)-WC[j][i]*P2W (j-l) — END----------------------------------------------------------------,
  121. FOR i:=1 to n do BEGIN j :=l-if i≤NG then
  122. MCnewj.1.:=MC[j][i]-MC[j][i]*P1M (j+1)+MC[j+1][i]*P2M (j) +WC[j][i]*P2M (j) MC[j][i]*P2W (j)else
  123. MCnewj.1.:=MC[j][i]-MC[j][i]*P1M (j+1)+MC[j+1][i]*P2M (j)-for j:=2 to nn-1 do
  124. MCnewj.1.:=MC[j][i]-HC[j][i]*P1M (j+1)-MC[j][i]*P2M (j-1)+ MC[j+l][i]*P2M (j)+MC[j-1][i]*PlM (j)-j:=nn-
  125. MCnewj. l]: =MC[j]1.+MC[j-l][i]*P1M (j)-MC[j][i]*P2M (j-l) — END-
  126. WCnewj. [n] :=WC new[j] [n]1. WCj. n-1]*P3W (n) — WCj. n]*P3W (1) — WCj. n]*P4W (n-l) + WC[j]1.*P4W (n)-1. MCnewj. n]: =MCnew[j][n]1. MCj. n-1]*P3M (n) — MCj. n]*P3M (1) — MCj. [n]*P4M (n-1) + MC[j]1.*P4M (n)-1. END-i SetColor (63) —
  127. OutTextXY (10,10, '-------------------') —
  128. SetColor (0) — Str (step, stroka) —
  129. OutTextXY (10,10,'Шаг '+stroka) — SetColor (63) —
  130. OutTextXY (42 0,70,'-----------------')-1. SetColor (0) —
  131. Str (10 0*SCW/SC:5:3,stroka) —
  132. OutTextXY (420,70,'SCP = '+stroka+' %')-1. SetColor (63) —
  133. OutTextXY (250,70,'-----------------')-1. SetColor (0) —
  134. Str (10 0*SCM/SC:5: 3, stroka) —
  135. OutTextXY (250,70,'SCL = '+stroka+' %')-1. SetColor (63) —
  136. OutTextXY (420,90,'-----------------')-1. SetColor (0) —
  137. Str (100*Kw:6: 3, stroka) —
  138. OutTextXY (420,90,'Kp = '+stroka+' %')-1. SetColor (63) —
  139. OutTextXY (250,90,'-----------------')-1. SetColor (0)-1. Str (10 0*Km:6:3,stroka) —
  140. OutTextXY (250,90,'K1 = '+stroka+' %')-1. SetColor (63) —
  141. OutTextXY (335,110,'-----------------')-1. SetColor (0) —
  142. Str (100*Kv: 6: 3, stroka) —
  143. OutTextXY (335,110,'Kv = '+stroka+' %')-1. Kv}
Заполнить форму текущей работой