Разработка новых конструкций вибрационных смесителей барабанного типа для сыпучих материалов и методики их расчета тема диссертации и автореферата по ВАК РФ 05.04.09, кандидат технических наук Пасько, Александр Анатольевич

• Пасько, Александр Анатольевич
• кандидат технических наук
• 2000

Оглавление диссертации кандидат технических наук Пасько, Александр Анатольевич

ВВЕДЕНИЕ

Глава 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. Разработан новый способ определения неоднородности смеси сыпучих компонентов.

Список литературы диссертационного исследования кандидат технических наук Пасько, Александр Анатольевич, 2000 год

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. №655419 СССР. Вибрационный смеситель. Б.и.№13,1979.

11. A.c. №1558449 СССР. Вибрационный смеситель. Б.и.№15,1990.

12. Патент РФ №2035986. Вибрационный смеситель.

13. A.c. №1499831 СССР. Вибрационный смеситель.

14. Патент РФ №2137536. Вибрационный смеситель.

15. Патент РФ №2122891. Вибрационный смеситель.

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. Каляпин Д.К., Пасько А.А., Негров B.JI., Першин В.Ф. Моделирование механических процессов в виброожижен-ном слое// Международная научная конференция «Жидкофазные системы и нелинейные процессы в химии и химической технологии»: Тез. Докл, Иваново, 1999.

54. 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.

55. У. Pershin, M. Sviridov, A. Pasko, A. Sherbakov, E. Mandrika Stochastie-determinate and determinate-stochastie mixing// 13th International Congress of Chemical and Process Engineering,

56. Praha, Czeh Republic, 1998,V.7,P.177.

57. Ю.И.Гусев, И.Н.Карасев, Э.Э.Кольман-Иванов и др. Конструирование и расчет машин химически производств. -М.: Машиностроение, 1985. 408 с.

58. А.с. СССР №1769084. Способ определения неоднородности смеси сыпучих компонентов. G 01N15/02, Пасько А.А. и др.

59. Селиванов Ю.Т., Казаков М.Т., Першин В.Ф., Ткачев А.Г. Исследование процесса приготовления многокомпонентных смесей. В сб.: Тезисы докладов I научной конференции ТГТУ, Тамбов, 1994, с. 81.

60. Першин В.Ф., Негров В Л., Селиванов Ю.Т. Исследование процесса смешивания полидисперсных материалов // Роль118мол. конструкторов и исследователей хим. машиностроения. Тезисы докладов V Всесоюзной науч.-техн. конференции. -Северодонецк, 1986. с. 28.

61. Першин В.Ф., Негров B.JI., Селиванов Ю.Т. Моделирование процесса смешивания полидисперсных материалов. В сб.: Тезисы докладов Всесоюзной конференции «Технология сыпучих материалов ХИМТЕХНИКА-86», Белгород, 1986, с. 49-50.

62. Об использовании результатов диссертационной работы Пасько A.A. «Разработка новых конструкций вибрационных смесителей барабанного типа для сыпучих материалов и методика их расчёта».

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

64. Президент научно-технической 4фирмы «ЛИОНИК» ПЧЕЛОПРОМа Минсельхозпрода России,кандидат технических наук, лауреатпремии Совета Министров СССР, ^^лауреат премии Правительства Е.А.Мандрыка

65. Российской федерации в областинауки и техники

66. Концентрации правого контура}

67. Концентрации правого контура}

68. Концентрации левого контура}

69. Концентрации левого контура}single; single ; single; single ;i : integer; j: integer; DWG:Boolean;f: text;1. Tangens} {*}

70. Function tan( x:single):extended;}1. Begin {*}tan:=sin(x)/cos(x); {*}1. End; {*}1. Proc 1 ***}

71. Procedure Procl ( xl,yl,x2,y2:single;var k:single);1. Begink:=(y2-yl)/(x2-xl) ; End;1. End Proc1

72. 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); {

73. 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

74. BB: =2*tan (beta) * (yc-tan (beta) *xc) ; CC:=sqr(yc-tan(beta)*xc)-sqr(R);

75. 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) ) ;

76. BB:=2*tan(beta)*(yc-tan(beta)*xc); CC : =sqr (yc-tan (beta) *xc) -sqr (R) ; 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));

77. BB:=2*tan(beta)*(yc-tan(beta)*xc); CC:=sqr(yc-tan(beta)*xc)-sqr(R); DD:=sqr(BB)-4*AA*CC; x: =(-BB+sqrt(DD))/(2*AA); y:=tan(beta)*(x-xc)+yc; end;end;

78. 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;

79. End Proc2 ***} {*** Begin PLATZ

80. 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;

81. PLATZ := sqrt( p*(p-a)*(p-b)*(p-c) ); End;1. END PLATZ)

82. 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

83. BB:=2*tan(beta)*(yc-tan(beta)*xc); CC:=sqr(yc-tan(beta)*xc)-sqr(R); 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) ) ;

84. BB:=2*tan(beta)*(yc-tan(beta)*xc);

85. CC:=sqr(yc-tan(beta)*xc)-sqr(R);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) ) ;

86. BB:=2*tan(beta)*(yc-tan(beta)*xc);

87. CC:=sqr(yc-tan(beta)*xc)-sqr(R);1. DD:=sqr(BB)-4*AA*CC;x:=(-BB+sqrt(DD))/(2*AA);y:=tan (beta) * (x-xc) +yc; end;end;

88. S1. :=sqrt( p*(p-a)* (p-b)*(p-c) SS:=SS+Si.; End; End;1. End Geron1. Begin NewCenter ***}

89. Возвращате ПРАВИЛЬНЫЕ значения центра сечения } Procedure New Center; var i:integer;

90. SO,Sx,Sy:single; Begin SO:=0; Sx:=0; Sy:=0;1. For i:=l to n do Begin

91. 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 ***}

92. 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;

93. End Vizual {*** Begin Center

94. Заполняет массив координатами центров зон) 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

95. Bj+l.1.x := (A[j][i]-x + A[j+1][i+1].x)/2;

96. Вj +1. [ i ] .у := (A[j] 1. .y + A[j+1] [i + 1] -y)/2;endelsebegin

97. В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;

98. 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 ***}

99. 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;

100. End MeinColor {*** Begin P1W

101. Function P1W( a:integer):single; Begin

102. P1W:=P01*(1-WCa.1.); End; {*** End P1W {*** Begin P2W

103. Procl(xx, yy, xc, yc, kl) ; Procl(0 , b, xc, yc, k2); Procl(0 , -R, xc, yc, k3) ; Proc2; Geron;1. NewCenter; {*}1. ClearDevice;) {*}

104. 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;

105. New Center; Proc2NEW; Geron;1. Zona; {*}for i:=l to n do for j:=nn-l downto 0 do begin

106. 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);

107. Заполним начальными концентрациями} Randomize; For i:=1 to n do for j:=1 to nn do Begin

108. WCj.1.:=0{random) ; MC[j][i]:=0{random} ;1. End;

109. SOUND(440);DELAY(1000) END;1. HALT;1:step:=step+l; {Покажем на экране} S CW: = 0; S CM: = 0; Xv:=0 j Xw: = 0 Xm:=0

110. For i:=l to n do for j:=1 to nn Begindo

111. SetFillStyle(l,MeinColor(WCj. 1. ) ) ;

112. FloodFill(getmaxx div 2 +round(k*Wcenterj.1.x),getmaxy div 2 round(k*Wcenterj.1.y),ColZona); SetFillStyle(l,MeinColor(MC[j][i]));

113. FloodFill(getmaxx div 2 +round(k*Mcenterj. 1. . x) ,getmaxy div 2 round(k*Mcenterj.1.y),ColZona);

114. Сумма концентраций} SCW:=SCW+WCj.1.; SCM:=SCM+MC[j][i];

115. Xv:=Xv+WCj.1.+MC[j][i]; Xw:=Xw+WC[j][i]; Xm:=Xm+MC[j][i]; End;1. SC:=SCW+SCM;

116. Xv:=Xv/(2*n*nn); Xw:=Xw/(n*nn); Xm:=Хш/(n*nn);

117. Коэфф. НЕОДНОРОДНОСТИ} Kv:= 0;1. Kw:= 0;1. Km:=0;for i:=1 to n dofor j : =1 to nn do Begin

118. 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

119. WCnewj.1.:=WC[j][i]-WC[j][i]*P1W(j+1) + WC[j+1][i]*P2W(j);for j:=2 to nn-1 do

120. 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;

121. WCnewj.1.:=WC[j][i]+WC[j-l][i]*P1W(j)-WC[j][i]*P2W(j-l); END;---------------------------------------------------------------,

122. FOR i:=1 to n do BEGIN j :=l;if i<=NG then

123. 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

124. MCnewj.1.:=MC[j][i]-MC[j][i]*P1M(j+1)+MC[j+1][i]*P2M(j);for j:=2 to nn-1 do

125. 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;

126. MCnewj.[l]:=MC[j]1.+MC[j-l][i]*P1M(j)-MC[j][i]*P2M(j-l); END;

127. 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);

128. OutTextXY(10,10, '-------------------') ;

129. SetColor (0); Str(step,stroka);

130. OutTextXY(10,10,'Шаг '+stroka); SetColor (63);

131. OutTextXY(42 0,70,'-----------------');1. SetColor (0);

132. Str(10 0*SCW/SC:5:3,stroka);

133. OutTextXY(420,70,'SCP = '+stroka+' %');1. SetColor (63) ;

134. OutTextXY(250,70,'-----------------');1. SetColor(0);

135. Str(10 0*SCM/SC:5 : 3,stroka) ;

136. OutTextXY(250,70,'SCL = '+stroka+' %');1. SetColor(63);

137. OutTextXY(420,90,'-----------------');1. SetColor (0) ;

138. Str(100*Kw:6 : 3,stroka) ;

139. OutTextXY(420,90,'Kp = '+stroka+' %');1. SetColor(63);

140. OutTextXY(250,90,'-----------------');1. SetColor(0);1. Str(10 0*Km:6:3,stroka);

141. OutTextXY(250,90,'K1 = '+stroka+' %');1. SetColor(63);

142. OutTextXY(335,110,'-----------------');1. SetColor(0);

143. Str(100*Kv: 6 : 3,stroka) ;

144. OutTextXY(335,110,'Kv = '+stroka+' %');1. Kv}