Научный журнал

Фундаментальные исследования

ISSN 1812-7339

"Перечень" ВАК

ИФ РИНЦ = 1,674

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Середа Т.Г. 1

---

1 ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет»

К физическим процессам, протекающим в массиве твердых бытовых отходов (ТБО), относятся тепло и массоперенос, растворение и диффузия. Движущей силой теплопереноса является разность температур, при этом температура оказывает большое влияние на весь спектр процессов, протекающих на полигоне ТБО. Воздействие температуры в значительной степени определяет характер протекания химико-биологических процессов и влияет на состояние системы в целом. В связи с этим разработана задача теплопроводности в массиве ТБО с использованием метода конечных элементов, включающего в себя дифференциальную и вариационную постановки задачи теплопроводности на полигоне ТБО. Решена задача теплопроводности с использованием метода конечных элементов. Для расчета температурных полей разработан пакет программ «work», предназначенный для определения стационарных и нестационарных температурных полей в плоских областях. В пакете реализован стандартный метод конечных элементов (МКЭ), применяемый для решения различных задач теплопроводности.

Статья в формате PDF

513 KB

твёрдые бытовые отходы (ТБО)

АСУТП

полигон ТБО

1. Артемов Н.И. [и др.] Технологии автоматизированного управления полигоном твердых бытовых отходов / Н.И. Артемов, Т.Г. Середа, С.Н. Костарев, О.Б. Низамутдинов // Международный журнал экспериментального образования. – 2010. – № 11. – С. 43.

2. Костарев С.Н., Середа Т.Г., Михайлова М.А. Системный анализ управления отходами. – LAP LAMBERT Academic Publishing GmbH & Co. KG, 2012.

3. Костарев С.Н. Математическая модель управления состоянием полигона твердых бытовых отходов : дис. … канд. техн. наук : 05.13.18. – Пермь, 2003. – 199 с.

4. Костарев С.Н., Середа Т.Г., Михайлова М.А. Программно-аппаратный комплекс управления качеством фильтрационных стоков // Экологические системы и приборы. – 2014. – № 3. – С. 39–46.

5. Костарев С.Н., Середа Т.Г., Михайлова М.А. Разработка автоматизированной системы мониторинга и управления природно-техническими системами утилизации отходов // Фундаментальные исследования. – 2013. – № 6–2. – С. 273–277.

6. Костарев С.Н., Середа Т.Г. Комплексное внедрение автоматизированных систем управления технологическими процессами на санитарных полигонах твердых бытовых отходов // Экологические системы и приборы. – 2014. – № 4. – С. 20–28.

7. Свидетельство о регистрации программы для ЭВМ 2009612494.– М.: ФИПС, 2009.

8. Середа Т.Г. [и др.]. Снижение пожаровзрывоопасности объектов депонирования отходов // Пожарная безопасность. – 2008. – № 3. – С. 84–89.

9. Середа Т.Г. Обоснование технологических режимов функционирования искусственных экосистем хранения отходов: дис. … д-ра техн. наук : 03.00.16; МГУП. – М., 2006.

10. Середа Т.Г., Костарев С.Н. Разработка методов проектирования автоматизированных систем обработки информации и управления искусственными экосистемами хранения отходов // Экологические системы и приборы. – 2006. – № 11. – С. 21–24.

11. Середа Т.Г. Источники загрязнения и методы защиты человека в техносфере: учеб. пособие для студентов вузов. – Пермь, Изд-во «ПНИПУ», 2012.

12. Сегерлинд Л. Применение метода конечных элементов: пер. с англ. – М.: Мир, 1979. – 390 с.

13. Таблицы физических величин: справочник / под ред. И.К. Кикоина. – М.: Атомиздат, 1976. – 1008 с.

14. Kostarev S.N., Sereda T.G. Аutomated process control of sanitary municipal solid waste landfill // World Applied Sciences Journal. – 2013. – Т. 22. – № SPL.ISSUE2. – С. 64–69.

15. Sereda T.G., Kostarev S.N., Elancheva E.N. Study safety environmental protection landfills using models anaerobic digesters // Applied Mechanics and Materials. – 2014. – Vol. 682. – Р. 339–345.

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

Целью представленной работы являлось моделирование процессов теплопереноса на полигонах ТБО.

Материал и методы исследования

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

Результаты исследования
и их обсуждение

Дифференциальная постановка задачи. Процесс теплопереноса на полигоне ТБО задан в области W в R3 c границей S, на части границы Si заданы граничные условия 1-го рода (изолированный экран), на части Sk заданы граничные условия 3-го рода (конвективный теплообмен) (рис. 1) [8, 9].

Уравнение теплопроводности в сплошной среде имеет следующий вид [2]:

(1)

где t – температура; КХХ, КУУ, КZZ – коэффициенты теплопроводности в направлениях х, у, z; Q – внутренний источник тепла, при граничном условии

(2)

где h – коэффициент теплообмена, кВт/м2; t – температура на границе; tос – температура окружающей среды; q – точечный источник тепла.

Рис. 1. Сечение полигона с граничными условиями 1-го и 3-го рода

Для нахождения приближенного решения дифференциальной задачи (1), (2) с помощью таких численных методов, как метод Ритца, Бубнова – Галеркина, конечных элементов (МКЭ), требуется построение вариационных аналогов исходных дифференциальных задач.

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

(3)

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

(4)

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

β = i, j, k. (5)

(6)

Матрица градиентов [В]:

(7)

Матрица свойств материала [D]:

(8)

Теперь объемный интеграл примет следующий вид:

(9)

Предполагая толщину элемента единичной, заменим dV на dA. Подынтегральное выражение в (9) постоянно и может быть вынесено за знак интеграла:

(10)

Вычисляя произведение матриц, имеем

(11)

Второй интеграл должен быть вычислен по поверхности:

Подставляя в матрицу [N] функции формы и выполняя мат­ричное умножение, получаем

(12)

Программная реализация

Для расчета температурных полей разработан пакет программ «work», предназначенный для определения стационарных и нестационарных температурных полей в плоских областях [4, 6, 7]. В пакете реализован стандартный метод конечных элементов (МКЭ), применяемый для решения различных задач теплопроводности. Функционально пакет разделен на две отдельные части: Grid – для разбивки рабочей области на конечные элементы; Vega – для решения плоской задачи теплопроводности. Рассмотрим основные этапы работы с программой «Work».

расчет теплофизических свойств фракций отходов

Для определения коэффициентов теплоемкости и теплопроводности использованы источники [11, 13] (таблица).

Теплофизические свойства фракций отходов

Разбивка рабочей области
на конечные элементы

Для разбивки рабочей области на конечные элементы (КЭ) необходимо предварительно разделить ее на зоны, каждая из которых образована восемью опорными узлами. Данную операцию можно провести либо вручную, либо с использованием специальных утилит для автоматического задания координат и номеров опорных узлов. До разбивки следует выполнить эскиз исследуемой области, выбрать необходимое число зон, опорных узлов, занумеровать зоны и узлы, нумерация сторон в каждой зоне ведется против хода часовой стрелки. Информация по разбивке формируется в файлах «Work» с расширениями in1, in2, in3. Далее запускается модуль «Maingrid». Происходит формирование файла граничных условий. Для формирования файла запускается программа «Bc_trm». Исходные данные для «Bc_trm» собирают в текстовом файле «Work.ggg». Расчет температурных полей создается в файле «Work.reg».

Для построения диаграммы температурных полей рис. 2 использована программа AXUM 5.0.

Рис. 2. Распределение температурных полей по сечению полигона

Применение методов системного анализа и кибернетики в управлении системами ПТО открывает возможность исследования и организации безопасного функционирования объекта депонирования ТБО в условиях, когда информация последовательно накапливается и реализуется в виде алгоритмов для ЭВМ. В этом случае управляющие действия будут направлены на минимизацию жизненного цикла природно-технических систем утилизации отходов, что имеет важное значение [14]. Внедрение автоматизированного рабочего места конструктора-проектировщика полигона ТБО уже на этапе проектирования позволит ускорить выполнение проектной документации и повысить безопасность эксплуатации природно-технической системы депонирования отходов на этапах её эксплуатации и рекультивации [8, 10, 15].

Козлов В.В., д.т.н., профессор кафедры ВМКСиС, Пермский военный институт внутренних войск МВД России, г. Пермь;

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

Работа поступила в редакцию 19.12.2014.

Библиографическая ссылка