Оптимальное управление процессом обжига шихты во вращающихся трубчатых печах в производстве глинозема и попутной продукции тема диссертации и автореферата по ВАК РФ 05.13.06, кандидат технических наук Титов, Олег Васильевич

• Титов, Олег Васильевич
• кандидат технических наук
• 2011

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

Введение.

Глава 1. СОВРЕМЕННОЕ СОСТОЯНИЕ И ПЕРСПЕКТИВЫ РАЗВИТИЯ ОБЖИГА ШИХТЫ ВО ВРАЩАЮЩИХСЯ ПЕЧАХ.

1.1 Основные характеристики вращающихся печей.

1.2 Температурный режим работы вращающихся печей.

1.3 Движение материала.

1.4 Сжигание газа в турбулентном потоке^.

1.5 Учет пылеуноса и пылевозврата.

1.6 Влияние цепных завес на процессы теплообмена, расчет коэффициента теплообмена с учетом конструкции и расположения цепных завес.

Глава 2. КИНЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ ПРОЦЕССА ОБЖИГА И ОПРЕДЕЛЕНИЕ КИНЕТИЧЕСКИХ КОНСТАНТ ДЛЯ ШИХТ РАЗЛИЧНОГО

СОСТАВА.

Глава 3. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ОБЖИГА

ВО ВРАЩАЮЩИХСЯ ТРУБЧАТЫХ ПЕЧАХ.

Глава 4. ИССЛЕДОВАНИЕ ПРОСТРАНСТВЕННОГО РАСПРЕДЕЛЕНИЯ КОНЦЕНТРАЦИЙ И ТЕМПЕРАТУР В ПЕЧИ С ИСПОЛЬЗОВАНИЕМ

ДВУХМЕРНОЙ МАТЕМАТИЧЕСКОЙ МОД ЕЛИ:.

Глава 5. РАЗРАБОТКА СИСТЕМЫ УПРАВЛЕНИЯ ПРОЦЕССАМИ ХИМИКО-МЕТАЛЛУРГИЧЕСКИХ ПРЕВРАЩЕНИЙ В ТРУБЧАТЫХ ВРАЩАЮЩИХСЯ ПЕЧАХ ПРИ ОБЖИГЕ ШИХТЫ.

5.1 Формирование критерия оптимизации.

5.2 Постановка задач управления.

5.3 Синтез системы автоматического управления температурным режимом трубчатой'вращающейся печи.

5.3.1' Структура системы управления.

5.3.2 Принцип работы регулятора с прогнозирующей моделью.

5.3.3 Прогнозирующая’Модель.

5.3.4 Оценка состояния.

5.3.5 Модель шума измерений.

5.3.6 Модель выходных возмущающих воздействий.

5.3.7 Состояние наблюдателя.

5.3.8 Задача оптимизации.

5.3.9 Квадратичное программирование.

5.3.10 Ограничения.

5.3.11 Решение алгоритма управления по прогнозирующей модели.

5.3.12 Алгоритм управления по прогнозирующей модели без ограничений.

5.3.13 Алгоритм управления по прогнозирующей модели с ограничениями

5.3.14 Применение разработанного алгоритма.

5.3.15 Построение системы управления в MATLAB.

5.3.16 Сравнение качества управления.

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

В мировой алюминиевой промышленности основным сырьем для производства глинозема служат высококачественные бокситы, перерабатываемые по способу Байера. Развитие глиноземного производства в России ориентировано в основном на использование собственной сырьевой базы. Вследствие ограниченных запасов байеровских бокситов на глиноземных комбинатах России широко используется« небокситовое высококремнистое сырье — нефелины. Запасы нефелинов» в России практически неисчерпаемы. Усилиями отечественных ученых были достигнуты огромные результаты по созданию эффективных схем их комплексной переработки, получившие мировое признание. В соответствии^ с разработанной технологией все компоненты сырья без остатка используются для получения глинозема, портландцемента и содопоташной продукции. Благодаря этому, нефелиновое сырье остаётся конкурентноспособным по отношению к бокситовому сырью.

Дальнейшее повышение эффективности комплексной переработки 1 нефелинового сырья*, энергоэффективность процесса в целом, тесно связаны с совершенствованием систем управления процессом переработки, таю как все остальные технологические ресурсы почти исчерпаны. Вклад в развитие и совершенствование алгоритмов управления объектами алюминиевой промышленности внесли специализированные научно-исследовательские и проектные организации, среди них «ВАМИ», «СибВАМИ», «Цветметавтоматика», «Союзцветметавтоматика» и др. Однако до настоящего времени остаются до конца нерешенными и требуют дальнейшего развития вопросы повышения качества управления в условиях отсутствия достаточной информации о параметрах процесса спекания^нефелиновой шихты в трубчатых вращающихся печах. В связи с этим возникает необходимость в разработке высокоэффективной автоматизированной системы управления процессом спекания нефелиновой шихты в трубчатой вращающейся печи, позволяющей вести этот процесс наиболее эффективно.

Работа выполнена в соответствии с госбюджетной тематикой НИР СПГУ по теме 6.30.020. «Разработка систем управления сложными техническими, объектами с использованием математических моделей в контуре управления» (I кв. 2008 - IV кв. 2010 гг.);

Цель работы:

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

Основные положения, выносимые на защиту •.

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

2. Для обоснования допустимости применения одномерной стационарной; математической модели трубчатой* вращающейся печи следует проводить, сопоставление результатов моделирования« в одномерной постановке с результатами моделирования в двумерной постановке с применением пакета программ типа Сотзо1 МиШрЫэюз и суперкомпьютера высокой производительности. Было определено оптимальное распределение температуры по длине печи путем решения задачи, оптимизации с применением технологического критерия оптимизации, (степени завершенности процесса) и уравнений связи варьируемых параметров с критерием оптимизации в виде: раз;работанной детальной математической модели процесса спекания.*

3 ^Введение в структуру автоматизированной системы управления процессом спекания нефелиновой шихты в трубчатой вращающейся* печи регулятора, реализующего алгоритмы управления по прогнозирующей модели, обеспечивает повышение качества управления по величине перерегулирования, времени переходного процесса и среднеквадратическому отклонению ошибки регулирования по сравнению с обычным ПИД - регулированием.

Методы исследований.

В процессе работы проводились экспериментальные исследования кинетики химических превращений при термической обработке шихты при получении глинозема и цементного клинкера. Эксперименты проводились с использованием'установки синхронного термического анализа* (TG DSC)STA-429 фирмы NETZSCH с получением двух экспериментальных откликов: изменения массы образца и скорости поглощения (выделения тепла) в процессе преобразования шихты различного состава < в алюминат натрия и компоненты клинкера при нагреве образца в ячейке с заданнной скоростью. Проводился также рентгенофазовый анализ образцов шихты с использованием рентгеновского дифрактрометра XRD -6000 (Shimadzu).

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

Научная новизна

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

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

Практическая значимость

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

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

3. Научные результаты работы используются в учебном процессе химикометаллургического факультета СПГУ для студентов специальности 220301 «Автоматизация технологических процессов и производств».

Достоверность результатов работы.

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

Апробация работы.

Содержание и основные положения диссертации докладывались, и обсуждались: на международной научной специализированной конференции в Германии в г. Фрайберг в 2009 г.; на научной конференции студентов и молодых ученых в СПГГИ (ТУ) им. Г. В. Плеханова в 2008г.; научных семинарах кафедры А11111 СПГГИ (ТУ) им. Г. В. Плеханова 2008-2010 г.г.; на научно-технических советах по работе с аспирантами СПГГИ (ТУ) им. Г. В. Плеханова 2008-2010 г.г.

Публикации. Основные положения диссертации опубликованы, в 4 статьях.

Структура и объем работы

Диссертация состоит из введения, пяти глав, заключения,

Заключение диссертации по теме «Автоматизация и управление технологическими процессами и производствами (по отраслям)», Титов, Олег Васильевич

Основные результаты выполненных исследований сводятся к следующему:

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

2. Созданные алгоритмы управления трубчатой вращающейся печью позволяют повысить эффективность управления на 40-50% по величине перерегулирования и на 10-15% по времени стабилизации температуры, с помощью прогнозирующей математической модели, включенной в контур управления процессом спекания нефелиновой шихты.

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

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

Заключение

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

Список литературы диссертационного исследования кандидат технических наук Титов, Олег Васильевич, 2011 год

1. Steffen; W. Hydrogen energy from coupled waste gasification and cementproduction a thermochemical concept study / W. Steffen, S. Hamel and W. Krumm // International Journal of Hydrogen Energy. — 2006. — Vol. 31. — Issue12. P. 1674-1689. .

2. Brimacombe J.K, Watbmison Л.Р. Heat transfer ini Direct-Fired Rotary kiln: III

3. Heat Flow Results and their Interpretation. Metallurgical Transaction. B, 1978 Vol, № 66 p. 209-218 ,

4. U. Kaantee, R. Zevenhoven, R. Backman, M. Hupa. Cement manufacturing using alternative fuels and the advantages of process modeling// Fuel Processing Technology №85. 2004 r. P. 293-301

5. Сизяков В.М., Шморгуненко Н.С. Научно-технический прогресс в производстве глинозема//Цветные металлы. М. 1981. №8, С.43-45. .

6. Мальц Н.С., Сизяков В.М., Шморгуненко Н.С. Новое о взаимодействиидвухкальциевого силиката с алюминатными растворами // Travo. 1974. №11.1. С.79-88. , .

7. Сизяков В.М., Яшунин П.В., Алексеев А.И. Повышение качества белитового шлама при комплексной переработке нефелинов // Цветная металлургия. 1980. №13. С.24-26.

8. Сизяков В.М. Диссертация доктора техн. Наук // Ленинград : ЛГИ им. Г.В.

9. Плеханова. 1983. . . ' . .

10. Сизяков В:М., Насыров Г.З. Эффективные способы комплексной-переработки небокситового алюминиевого сырья на глиноземные и попутные продукты// Цветные метал л ы. 2001. №12: С.63-68.

11. Фолио А. Теплопередача во вращающейся печи. Л.: Гипроцемент, 1957.

12. Ю.Ходоров Е.И. Печи цементной промышленности. Л.,.1968.

13. П.Екимов В.А., Ходоров Е.И. / Исследование степени равномерности температуры материала на открытой, поверхности слоя во- вращающейся;печи. // Производство глинозема. Труды ВАМИ. Л.:

14. ВАМИ, 1974. Вып. 88. с. 58-71.

15. Кауфман Б.Н. Теплопроводность строительных материалов. М.: Госстройиздат, 1965.

16. Лыков A.B.//Инженерно-физический журнал, 1960, № 10. с. 88-93.

17. Екимов В.А., Абрамов В.Я., Сидорова Э.В. // Вращающиеся; печи для; спекания нефелиновых шихт. ЦНИИНцветмет, 1963. с. 35-44.15 .Чудновский А.Ф. Теплофизические характеристики дисперсных материалов. Физматгиз; 1962.

18. Срибнер Н.Г. Анализ движения сыпучих материалов во вращающихся барабанах / Химическая промышленность, №4, 1979. с: 232-235.

19. Ходоров Е. И., Екимов В. А., Лимонов Б; С. / Исследование длительностипребывания материала: на открытой поверхности слоя во вращающейся печи . // Производство глинозема. Труды ВАМИ. Л.: ВАМИ, 1974. Вып. 88. с. 4857. . • ; ■. .

20. Ходоров Е.И. Печи цементной промышленности. Л., 1968.

21. Ходоров Е.И. Движение материала во вращающихся печах. Промсгройиздат,1957. ' • :

22. Клейн F.K. Строительная механика сыпычих тел. М.: Госстройиздат, 1956.

23. Барановский В.В., Екимов В.А. / Изучение процесса движения материала во■, (вращающейся; печи спекания глиноземного: производства // Цветные металлы, 1962, № 6.

24. Линчевский В.П. Топливо и его сжигание. М.: Металлургия; 1959. — 400 с. 231Вулис Л:А., Ярин Л:П; Аэродинамика факела. Л.Энергия, 19781- 216 с.

25. Абрамович Г.Н. Теория турбулентных струй. М.: Физматгиз, 1960. 715 с.

26. Вулис Л.А., Кашкаров В.П. Теория струй вязкой жидкости; М.: Наука, 1965. -429 с.

27. Гиневский A.C. Теория турбулентных струй и следов. М:: Машиностроение; 1969.-398 с.

28. Госмен А.Д., Пан В;М., Ранчел А.К. Численные методыисследования течения вязкой жидкости. М:: Мир; 1972. — 320 с.

29. Гаусорн В., Уиддел Д., Хотел Г. Смешение и: горение в турбулентныхразовых струях. — В< кн.: Вопросы горения. Т. 1. М.: Изд-во иностр- лит., 1953. с. 146-193. !

30. Вулис Л:А. Об эффективном; управлении распространением свободной .турбулентной? струи / Л.А. Вулис, Ю;И; Михасенко; В:И1Хитриков // Механикажидкостшшгаз^.1966. Жб^С. Ь73-Г781 .

31. Кузнецов О.А., Ярин Л.П. Исследование влияния низкочастотных пульсаций на аэродинамику газового/факела. В кн.: Теория и практика сжигания газа. Вып. 5. Л.: Недра, 1972. с. 53-56.

32. Данциг С.Я. Пивоваров. В.В; . Нефелиновые . породы комплексное; алюминиевое сырье : монография / Данциг С.Я., Андреева Е.Д., Пивоваров В.В. .и др//М. : Недра, 1988: - 192 с.,

33. Сизяков В.М., Шморгуненко Н.С. Научно-технический прогресс в производстве глинозема // Цветные металлы. М: 1981. №8, 0.43-45.

34. Лайнер А.И;, Еремин Н.И:,, Лайнер Ю.А., Певзнер; И З. Производство:глинозема//М., Металлургия, 1978г., 344с. '

35. Влодавец, Н.И. Сборник “Нерудные материалы” //

36. Ленинград : Академия Наук СССР, 1929, Т. IV, стр. 56-66.

37. Влодавец, Н:И. Сборник “Хибинские апатиты” // Ленинград' : Госниитехиздат, 1931 г., Т. II, стр. 69-79.

38. Влодавец Н.И. // Доклады Акад.Наук СССР, 1931, стр. 127-130.

39. Сборник “Хибинские апатиты”. 1932 г., Т. IV, стр. 47-97.

40. Строков Ф.Н., и др. Труды ГИПХ. 29. 1936, стр. 13-26.

41. Мазель, В.А. К вопросу о рациональном составе алюминатного' спека« // Труды научно-исследовательского института легких металлов "НИИСАЛЮМИНИЙ", 1932'г., стр. 37.

42. Мазель, В.А. Предварительные опыты по совместному спеканию бокситов с нефелинами // Труды научно-исследовательского института легких металлов "НИИСАЛЮМИНИЙ", 1933 г., Т. 3, стр. 51-59.

43. Строков Ф.Н., и др.' Труды ГИПХ. 29. 1936, стр. 183-185.

44. Друцкая, Г.В., и'др. Отчет ВАМИ по теме 34Я 1-457 // Лениград : ВАМИ. Отчет ВАМИ. 1949.

45. Друцкая, Г.В. и др. Отчет ВАМИ по теме 272/1, часть 5 // Лениград : ВАМИ. Отчет ВАМИ. стр. 1-40.

46. Мазель, В.А. и др. Краткий отчет по технологическому обследованию печей спекания ВАЗ. // Лениград : ВАМИ. Отчет ВАМИ. 1950. стр. 1-28.

47. Отчет ВАМИ по теме 25 // Лениград : ВАМИ: 1951. стр. 1-18.

48. Отчет ВАМИ по теме 25 // Лениград : ВАМИ. 1952. стр. Г-47.50:Канцеленбоген П.Д., Крочесвский-Н.А., Смирнов М.Н. // Цветные металлы, 1957, №4, стр. 37-43.

49. Талмуд И.Л. //Биллютень ЦНИИЦМ, 1957, №8, стр. 30-36.

50. Мазель, В.А. // Труды ВАМИ. Ленинград. 1953. №34. С. 14-23.

51. Теория цемента / Под ред. А. А. Пащенко. — Киев: Будівельник, 1991. — 168 с

52. Н:И.Ерёмин, А.Н.Наумчик, В.Г.Казаков ’’Процессы и аппараты глинозёмного производства М.: Металлургия, 1980’

53. Шариков Ю.В. Моделирование процессов в металлургии/

54. Ю.В.Шариков, Н.В. Данилова, B.C. Зуев.СПб,СПГГИ,2007,81 с

55. Власов К. П., Теория автоматического управления, учебное пособие: -Харьков: Гуманитарный центр, 2006 531 с.

56. Страшинин Е. Э., Основы теории автоматического управления, 4.1линейные непрерывные системы управления, учебное пособие:

57. Екатеринбург: УГТУ, 2000 217 с.

58. Лазарева Т. Я., Мартемьянов Ю. Ф., Основы^ теории* автоматического управления, учебное пособие, 2-е изд., перераб. и доп.: Тамбов: ТГТУ, 2004 -352 с.

59. Александров А. Г., Артемьев В. М., Афанасьев В. Н-., Ашимов.А. А. и др., Справочник по теории автоматического управления, под ред. Красовского А. А.: М.: Наука, 1987 - 712 с.

60. Юревич Е. И.', Теория автоматического управления, учебное пособие: -СПб.: БХВ, 2007-326 с.

61. Пупков К. А., Егупов Н*. Д., Баркин А. И., Воронов Е. М., Методы классической и современной теории автоматического управления; учебник, Т.4, 2-е изд., перераб. и доп.: М.: МГТУ, 2004 - 744 с.

62. Васильев Ф. П., Методы оптимизации: М.: Факториал Пресс, 2002 - 824 с.

63. Гилл Ф., Мюррей У., Райт М., Практическая оптимизация, пер. с англ.: -М.: Мир, 1985 — 509 с.

64. Громов Ю: Ю., Земской Н. А., Лагутин*А. В., Иванова О. Г., Тютюнник В. М., Специальные разделы теории управления. Оптимальное управление динамическими системами, учебное пособие: Тамбов: ТГТУ, 2004 - 108 с.

65. Ким Д. П., Теория автоматического управления, Т.1, Линейные системы: -М.: Физматлит, 2003 288 с.

66. Ким Д. П., Теория автоматического управления, Т. 2. Многомерные, нелинейные, оптимальные и адаптивные системы, учебное пособие: М.: ФИЗМАТЛИТ, 2004 - 463 с.

67. Андрющенко В. A., Теория систем автоматическогоуправления, учебное пособие: Л:: ЛГУ, 1990 - 256 с.

68. Jay С. Hsu, Andrew U. Meyer, Modem control’principles and applications: -New-York: McGraw-hill book company, 1972 544 c.

69. Бесекерский В. A., Попов Е. П., Теория систем автоматического управления: СПб: Профессия, 2003 - 752 с.70; Доронин С. В.,. Теория автоматического управления и регулирования, учебное пособие: Хабаровск: ДВГУПС, 2005 - 127 с.

70. Солодовиков В. В., Теория автоматического управления техническими системами, учебное пособие: М.: МГТУ, 1993 - 492 с.

71. Цыпкин Я: 3., Основы теории автоматических систем: М;, Наука, 1977 -560 с.

72. Рощин. А. В., Основы теории автоматического управления, учебное пособие, 2-е ИЗД;, стер- М:: МГУПИ, 2007 - 100 с.

73. Лукас В. А1., Теория управления техническими системами, учебный' курс для вузов, 3-є изд., перераб. и доп.: Екатеринбург: УГГГА, 2002 - 675 с.

74. Alberto Bemporad, Manfred Morari, N. Lawrence Ricker, Model predictive control; Math Works, 2010 205 с.

75. Брюханов В. H., Косов М. Г., Протопопов С. П., Соломенцев Ю. М.,

76. Теория автоматического управления, 3-є изд., стер.: М.: Высшая школа, 2000 - 270 с. .

77. Филлипс Ч., Харбор Р., Системы управления с обратной связью: М.:: Лаборатория ¡Базовых Знаний, 200Ъ - 616 с.

78. Заболотнов Ю. М., Оптимальное управление непрерывными динамическими системами: Самара: СГАУ, 2005 -129 с.79: Черноруцкий И: Г., Методы оптимизации в теории управления, учебное пособие::- СПб.: Питер, 2004 256 с.

79. Reklaitis G: V., Ravindan A., Ragsdell K. M., Engineering Optimization. Method and Applications;, P.l: NewYork: Wiley-Interscience, 1983 -350 c.

80. Reklaitis G. V., Ravindan A., Ragsdell К. М., Engineering

81. Optimization. Method and Applications, P.2: NewYork: Wiley-Interscience, 1983 - 320 c.

82. А. А. Гальнбек, JI.M. Шалыгин, Ю.Б. Шмонин „Расчётыпирометаллургических процессов и аппаратуры цветной металлургии”: Челябинск: Металлургия, Челябинское отделение, 1990'

83. Классен, В:К. Совершенствование методов испытаний и оптимизации режима работы вращающихся* печей / В.К. Классен, И.Н. Борисов, А.Н. Классен // Тр. НИИЦемента. 1985. — № 88. - С. 97—118.

84. Борисов, И.Н. Некоторые особенности поведения шламов на цепях при высушивании / И.Н. Борисов, Л.С. Дурнева // Цемент и его применение. — 2002.-№5. С. 26-28.

85. Борисов, И.Н. Пути повышения эффективности работы цепных завес вращающихся печей / И.Н. Борисов // Цемент и его применение. 2003. — №1. З.-С. 17-20

86. Ходоров Е.И., Вольперт М.Е., Тимофеев Г.А. / Исследование процесса пылеуноса на модели вращающейся печи // Производство глинозема. Труды ВАМИ. Л.: ВАМИ, 1974. Выт 88. с. 72-82'.

87. Клейн Г.К. Строительная механика сыпычих тел. М.: Госстройиздат, 1956.

88. Срибнер Н.Г. Унос сыпучих материалов из вращающихся барабанов / Химическая промышленность, №11, 1978. с. 861-864.

89. Ходоров Е.И., Вольперт М.Е., Тимофеев Г.А. Труды. ВАМИ, 1974, №88; с. 72.

90. Ходоров Е.И., Стрелкова Е.Н. Вращающиеся печи- спекания глиноземного производства. ЦИИНцветмет, 1962. вып. 1, с: 4-16.

91. Красовский А. А., Статистическая, теория переходных процессов в системах управления: М., Наука, 1968 - 240 с.

92. Симановский А. Ю., Типы регуляторов, методика настройки регуляторов, инструкция: Ивано-Франковск: html-формат, 2004 - 63 с.118холодильника)function sys,xO,str,ts. = sfunmy2(t,x,u,flag,Nt,TbO)

93. SFUNDSC2 Example unit delay М-File S-fimction % The M-file S-fimction is an exaple of how to implement a unit % delay.

94. Handle discrete state updates, sample time hits, and major time step % requirements.end mdlUpdate %0/0===========mdlOutputs

95. Return the output vector for the S-function0/0=============function sys = mdlOutputs(t,x,u,Nt,TbO)в нуле инициализация if t = 0; dp = u(4) rokl = 1650;

96. Hb = u(2)./(rokl.*4.*u(3)); for i = 1 : Nt;1. H(i) = Hb;1. Tb(i)= TbO;rw(i)=row(dp,Tb(i),H(i));end;assignin('base','H',H); assignin('base','Tb',Tb); assignin('baseVrw',rw); for ix = 1 :Nt;1. Tb = evalin('base','Tb');

97. H = evalin('baseVH'); rw = evalin('base','rw');dp u(4);for i=l :Nt;rw(i) =ro w(dp,Tb(i),H(i)); end;tbout = tbix(Tb,rw,H,Nt,TbO); proiz = dtc(tbout,Tb,rw,H,Nt,TbO); SampleTime = evalin('baseVSampleTime');

98. Tb = Tb SampleTime.*proiz;assignin('baseVtbout',tbout); tccp=tvtor(tbout,rw,Nt);zrasx=rasx(rw,Nt);сдвиг вправо for i = Nt: -1 : 2;1. Tb(i) = Tb(i-l);1. H(i) = H(i-l);end;

99. Tb(l) = u(l); rokl = 1650;

100. Hb = u(2) ./(rokl. * 4. * u(3)); H(l) = Hb;assignin('base',Tb',Tb);assignin('baseVH',H);assignin('base','rw',rw);

101. SampleTime = evalin('baseVSampleTime');x(l) = Tb(Nt);% ткл вых x(2) = H(Nt);x(3) = tccp; % вторич воздуха x(4) = zrasx; %расход воздуха sys x;end;end;сдвиг вправо1. Tb = evalin('base','Tb');

102. H = evalin('baseVH'); rw = evalin('baseVrw');dp = u(4); for i=l :Nt;rw(i) =row(dp,Tb(i),H(i)); end;tbout = tbix(Tb,rw,H,Nt,TbO); proiz = dtc(tbout,Tb,rw,H,Nt,TbO); SampleTime = evalin('base','SampleTime');

103. Tb = Tb SampleTime.*proiz;assignin('baseVtbout',tbout); tccp=tvtor(tbout,rw,Nt); zrasx=rasx(rw,Nt);сдвиг вправо for i = Nt: -1 : 2;1. Tb(i) = Tb(i-l);1. H(i) = H(i-l);end;

104. Tb(l) = u(l); rokl = 1650;

105. Hb = u(2)./(rokl.*4.*u(3)); H(l) = Hb;assignin('base','Tb',Tb); assignin('base','H',H); assignin('base' ,'rw',rw);

106. SampleTime = evalin('base','SampleTime');x(l) = Tb(Nt);% ткл вых x(2) = H(Nt);x(3) = tccp; % вторич воздуха x(4) = zrasx; %расход воздухаsys = x; ,end mdlOutputs %function sys=mdlGetT imeOfNext V arHit(t,x,u,Nt)sampletime = sampleTime;

107. SampleTime = evalin('baseVSampleTime'); dx = evalin('baseVdx');

108. SampleTime = dx./(u(3)) ; % Example, set the next hit to be one second later.assignin('baseVSampleTime',SampleTime); sys = t + 9* SampleTime;

109. Вычисление температуры газа на выходе function mean=tbix(tca,rwa,Hh,Nt,TbO);

110. TbO = TbO + 273; kil = 22.434; ki2 = 30.358e-3; kml = 2.394e-8; km2 = 8.966e-ll; kmu = 1.497e-6;кс = 4677; klam = 1.199е-3; d = 1.6е-3; dbig = 0.0031181; tea = tea + 273;

111. A = 4*0.13*klam/(0.2*(d.A(5/3))*(kmu.A(l/3))); for i= 1 :Nt;1. A1 = A / (rwa(i).A(2/3));

112. В = kil *kc/(tca(i)) + kil + ki2*kc + tca(i)*ki2; ■tbi(i) = tca(i) (tca(i) - TbO)*exp(-Al*Hh(i)/B); end;tbi = tbi 273; mean = tbi;

113. Вычисление температуры через время dt для клинкера function mean=dtc(tboz, tea, rwa,Hh,Nt,TbO); d = 1.6e-3;kc = 4677; cc = 0.96*1000; dbig = 0.0031181; rokl = 1650;

114. TbO = TbO + 273; tea = tea+ 273; tboz = tboz + 273; for i = 1 :Nt;proiz(i) = d*d*rwa(i)*0.2*((tboz(i)+kc)*tboz(i)-(Tb0 + kc)*TbO)/(rokl*cc*(dbig*dbig-d*d)*Hh(i)); end;mean = proiz;

115. Nd = 4*4*1.77/(pi*dbig*dbig); mean = 0; for i = 1 :Nt;mean = mean + (0.25*pi*d*d)*rw(i)*dx*4/(0.25*pi*dbig*dbig* 1.293); end;