Математическая модель процесса подачи и выброса грунта рабочими органами комбинированной машины для тушения лесных пожаров

Размещено на http://www.allbest.ru/

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА ПОДАЧИ И ВЫБРОСА ГРУНТА РАБОЧИМИ ОРГАНАМИ КОМБИНИРОВАННОЙ МАШИНЫ ДЛЯ ТУШЕНИЯ ЛЕСНЫХ ПОЖАРОВ

Драпалюк Михаил Валентинович, д.т.н., профессор

Бартенев Иван Михайлович, д.т.н., профессор

Гнусов Максим Александрович, аспирант

Дручинин Денис Юрьевич, ассистент

Воронежская государственная лесотехническая академия, Воронеж, Россия

Марков Олег Борисович, к.т.н., доцент

Петрозаводский государственный университет, Петрозаводск, Россия

Клубничкин Евгений Евгеньевич, к.т.н, доцент

Московский государственный университет леса, Мытищи, Россия

Введение

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

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

Постановка и решение задачи

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

Лесная почва является сложным для моделирования объектом, ввиду своего типологического разнообразия и большого количества физико-математических параметров (тип почвы, ее твердость, влажность, плотность, рельеф поверхности) [1, 2]. В то же время движение потоков почвы лежит в основе принципа работы грунтомета, поэтому модель почвы должна быть максимально близкой к реальности. Для моделирования почвы используется метод конечных элементов [3, 4]. В рамках метода почва представляется совокупностью большого количества шарообразных элементов малого размера, способных взаимодействовать как между собой, так и с рабочими поверхностями машины (с роторами и дисками). С геометрической точки зрения взаимодействие рабочей поверхности с почвой представляет собой задачу о нахождении расстояния r_в от некоторой поверхности до сферической поверхности произвольного шарового элемента почвы. При этом вязкоупругая сила , действующая на элемент почвы, определяется по формуле

, (1)

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

Наибольшую сложность в данной модели представляет расчет сил, действующих со стороны рабочей поверхности на элементы почвы и приводящих к эволюции почвенной системы. Кроме того производится обратный расчет: суммарной силы со стороны почвенных элементов на рабочие поверхности, оказывающих сопротивление их движению. Состояние каждого элемента почвы i определяется шестью переменными: координатами его центра (x_i, y_i, z_i) и тремя составляющими скорости (v_xi, v_yi, v_zi). Взаимодействие шаров между собой носит вязкоупругий характер, что позволяет адекватно учитывать упругость почвы и, при рассмотрении ее как среды для движения рабочей поверхности, высокую ее вязкость [2]. Кроме того, в модели шары почвы первоначально связаны между собой, и при удалении двух шаров друг от друга рабочими поверхностями происходит отрыв - вязкоупругое взаимодействие отключается. Это решение позволяет адекватно учесть в модели первоначальное связанное состояние почвы и способность ее к разделению на фрагменты, рассыпчатость почвы. Также предусмотрена возможность имитации "слипания" почвы - отделенные шары притягиваются друг к другу на малых расстояниях, хотя это притяжение значительно более слабое, чем было ранее, когда почва не была разделена на фрагменты. Таким образом, можно производить моделирование почв с высоким содержанием влаги.

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

, (2)

где F^У_ij и F^В_ij - силы упругого и вязкого взаимодействия шаров i и j (жирным шрифтом здесь и ниже выделены векторные величины); N_Ш - общее количество шаров в модели почвы.

При расчете сил для каждой пары шаров предварительно вычисляется расстояние r_ij между их центрами S_i(x_i, y_i, z_i) и S_j(x_j, y_j, z_j) (рисунок 2.1):

. (3)

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

Шары, находящиеся на расстояниях r_ij > d_Ш + d_М (где d_Ш - диаметр шара, d_М - расстояние притяжения), считаются не взаимодействующими (F^У_ij = 0, F^В_ij = 0) (рисунок 2). При меньших расстояниях упругие силы F^У_ij вычисляются следующим образом

1) Если r_ij < d_Ш, то

;

; (4)

.

2) Если d_Ш < r_ij < d_Ш + d_М/2, то

;

; (5)

.

3) Если d_Ш + d_М/2 < r_ij < d_Ш + d_М, то

;

; (6)

.

Здесь F^У_xij, F^У_yij, F^У_zij - декартовы составляющие силы F^У_ij; c - жесткость упругого взаимодействия шаров, которая может принимать два значения: с = c_c, если шары i и j связаны, с = c_р - если шары i и j не связаны.

Рисунок 2 - Зависимость силы взаимодействия двух шаров i и j от расстояния между ними

Функция F_ij(r) не должна иметь разрывов во всей области изменения r_ij, чтобы облегчить интегрирование уравнений движения. Поэтому при выборе вида зависимости F_ij(r) в интервале d_Ш + d_М / 2 < r_ij < d_Ш + d_М использовали зеркально отраженную относительно вертикали ветвь из интервала d_Ш < r_ij < d_Ш + d_М / 2.

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

;

; (7)

,

где v_xi, v_yi, v_zi и v_xj, v_yj, v_zj - декартовы составляющие скоростей i-го и j-го шара; k_В - коэффициент вязкости.

Для F^В_ij выбрана общепринятая прямо-пропорциональная зависимость вязкой силы от скорости движущегося в среде тела. При этом введен дополнительный коэффициент (r_ij - (d_Ш + d_М)), характеризующий взаимное проникновение шаров друг в друга, так как торможение шаров друг от друга зависит от степени их контакта.

При взаимодействии двух шаров, согласно третьему закону Ньютона, силы возникают парами, то есть

F_ji = - F_ij.

В соответствии со вторым законом Ньютона запишем уравнения движения i-го шара

;

; (8)

,

где m_Ш - масса шара; t - время; g - ускорение свободного падения.

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

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

В ходе компьютерного эксперимента шары могут двигаться внутри прямоугольного параллелепипеда размерами L_X L_Y L_Z (рисунок 3). В начале эксперимента шары, в несколько слоев, располагаются в нижней части параллелепипеда. При этом шары образуют случайную плотную упаковку (подготовленную заранее по алгоритму встряски). Движущаяся вдоль направления OY машина, взаимодействуя дисками и роторами с шарами почвы, удаляет их из слоя почвы и заставляет двигаться по баллистической траектории. При этом за машиной образуются характерная борозда, близкая к реальной по своей конфигурации.

Рисунок 3 - Объем L_X L_Y L_Z, в котором проводится моделирование

Если в ходе компьютерного эксперимента какой-либо шар выходит за пределы параллелепипеда, то есть выполняется одно из условий: x_i < 0, x_i > L_X, y_i < 0, y_i > L_Y, z_i < 0, z_i > L_Z, то этом случае шар возвращается в параллелепипед путем зеркального отражения от соответствующей стенки. Так, например, если выполнилось условие x_i > L_X, координате x_i шара присваивается значение x_i = L_x - (x_i - L_x), а скорость шара в направлении OX меняет знак: v_xi = - v_xi. При отражении шара от стенки также производится незначительное гашение его энергии по адаптированным формулам (7).

Диаметр шаров d_Ш влияет на адекватность модели почвы. Для более точного представления почвы d_Ш должен быть как можно меньше. В то же время количество шаров N_Ш, необходимое для представления заданного объема почвы V_П, резко возрастает по кубическому закону с уменьшением d_Ш. При увеличении же количества шаров увеличивается количество уравнений, описывающих почву и, соответственно, время расчета на ЭВМ. Даже уменьшение диаметра шара вдвое приводит к увеличению длительности расчета минимум в восемь раз. Поэтому, если ориентироваться на вычислительные возможности современных ЭВМ, оптимальный диаметр шаров, позволяющий с достаточной точностью описать почву, но одновременно не приводящий к сверхдлительным компьютерным расчетам, составляет порядка 7 см. При этом общее количество шаров в модели составляет 7000.

Перед проведением компьютерного эксперимента производится специальная процедура подготовки плотноупакованного объема почвы. Сначала необходимое количество шаров распределяется случайным образом по всему объему параллелепипеда. Затем производится интегрирование уравнений движения (8) в течении промежутка времени 15 с. При этом под действием сил тяжести шары начинают оседать и формировать в нижней части параллелепипеда слой почвы. Для создания более плотной упаковки в промежуток времени от 5 до 10 с производится встряска почвы: координатам (x_i, y_i, z_i) всех шаров на каждом шаге интегрирования сообщают незначительные случайные смещения. В процессе процедуры подготовки почвы скорости шаров, первоначально возникшие вследствие свободного падения, постепенно гасят, чтобы к концу процедуры шары образовали неподвижный плотноупакованный слой почвы. Впоследствии, в ходе компьютерного эксперимента, шары приобретают скорости за счет взаимодействия с рабочими органами грунтомета.

В конце процедуры подготовки почвы производится связывание шаров. Для этого заполняется двумерный массив K(i, j), каждая ячейка которого принимает значение "истинно" в случае, если шары i и j находятся на расстоянии

r_ij < d_Ш + d_М, (9)

и значение "ложно" в противоположном случае. Во время проведения эксперимента, на каждом шаге интегрирования для каждой связанной пары i-j проверяется, не разделилась ли она, по тому же условию (9), и в случае невыполнения (9) в массив K(i, j) записывается значение "ложно".

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

В модели необходимо корректно воспроизвести сложную геометрическую форму рабочих поверхностей, а также взаимодействие поверхностей с элементами-шарами почвы. В рамках конечно-элементного подхода поверхности сложной формы обычно заменяются большим количеством плоских фигур [3]. Как правило, в качестве элементарной фигуры используют треугольники, так как они легко стыкуются между собой и позволяют воспроизводить поверхность любой формы. Поэтому в рамках разрабатываемой модели рабочие поверхности представляются в виде совокупности элементарных треугольников. Так как в данной работе необходимо оптимизировать форму и геометрическое расположение поверхностей, то представление поверхностей в виде совокупности треугольников очень удобно тем, что можно легко изменить, например, размер ротора, изменяя координаты лишь нескольких вершин треугольников.

Треугольник в трехмерном пространстве задается координатами трех его вершин T_i1(x_i1, y_i1, z_i1), T_i2(x_i2, y_i2, z_i2), T_i3(x_i3, y_i3, z_i3), где T - обозначение точки-вершины треугольника; i - номер треугольника; индексы 1, 2 или 3 - означают номер вершины для i-го треугольника. Для образования сложных поверхностей треугольники стыкуются между собой по какому либо ребру: для соседних треугольников совпадают по две вершины (например T₇₁ = T₈₁, T₇₂ = T₈₂). Так как поверхности не являются замкнутыми, после стыковки всех треугольников остаются свободные ребра, которые выступают, например, в качестве режущих кромок роторов и дисков.

Рабочие поверхности в модели с некоторой степенью загрубления были разбиты на 64 элементарных треугольника (рисунок 4). При этом поверхность каждого ротора состояла из восьми треугольников, а поверхность каждого диска из 24 треугольников.

Рисунок 4 - Представление в модели рабочих поверхностей грунтомета как совокупности элементарных треугольников (показаны три проекции)

С течением времени рабочие поверхности совершают сложное поступательно-вращательное движение по отношению к почве. При этом необходимо непрерывно корректировать координаты точек-вершин треугольников T_ij (где i - номер треугольника, j - номер вершины). В рамках приближения о равномерном прямолинейном движении рабочих органов грунтомета вдоль оси OY необходимо изменять все игрековые координаты треугольников:

y_ij = y_0ij + v_с • t, (10)

где y_0ij - координаты вершин треугольников в момент начала моделирования; v_с -скорость поступательного движения машины; t - время, прошедшее с начала моделирования.

В модели также необходимо учесть постепенное заглубление машины в первые секунды движения. В момент времени t = 0 диски и роторы находятся над почвой (нижние точки касаются поверхности почвы), и за промежуток времени t_заг они должны выйти на заданную глубину a. Поэтому до момента времени t = t_заг координаты z_ij всех вершин треугольников изменяются по следующему линейному закону:

z_ij = z_0ij - a • t / t_заг, (11)

где z_0ij - координаты вершин треугольников в момент начала моделирования.

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

Проверка контакта шара j с треугольником i производится в два этапа: сначала находится проекция центра шара (x_j, y_j, z_j) на плоскость, содержащую треугольник, затем проверяется, попадает ли точка проекции внутрь треугольника. Поэтому сначала составляется уравнение плоскости, проходящей через три точки - вершины T_i1, T_i2, T_i3 треугольника по известной в аналитической геометрии формуле

(12)

После раскрытия определителя последнее уравнение приводится к виду

(13)

где A, B, C, D - коэффициенты уравнения плоскости.

Затем уравнение плоскости нормируется путем деления всех коэффициентов уравнения на длину нормального вектора :

(14)

или, введя новые обозначения коэффициентов (с индексом "н"), нормированное уравнение плоскости можно записать в следующем виде

. (15)

Тогда отклонение r_вн центра шара (x_j. y_j. z_j) от плоскости рассчитывается по формуле

(16)

Так как большинство элементов почвы не контактируют с рабочими органами грунтомета, на данном этапе расчета производится отсечение всех не контактирующих элементов почвы по условию r_вн > d_Ш/2.

Затем необходимо найти проекцию центра шара на плоскость данного треугольника. Так как известны координаты нормального вектора A_н, B_н и C_н и отклонение r_вн от точки до плоскости, проекция P(x_П, y_П, z_П) находится смещением по нормальному вектору

(17)

Для того, чтобы проверить, попадает ли точка P внутрь треугольника T_i, используется метод проверки по площадям: сумма площадей трех треугольников с вершиной P должна равняться площади всего треугольника T_i, то есть,

(18)

Если суммарная площадь треугольников , и оказывается больше площади элементарного треугольника , это означает, что точка P не лежит внутри треугольника T_i. и шар не контактирует с данным треугольником.

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

Разработанная модель в целом представляет собой систему из большого количества дифференциальных и алгебраических уравнений, а также условий включения тех или иных сил. Решение системы дифференциальных уравнений производится численно. Используется модифицированный метод Эйлера-Коши [5], который особенно эффективен при решении дифференциальных уравнений второго порядка. При этом координаты и скорости рассчитываются по формулам вида

(19)

где x, v, a - координата, скорость и ускорение некоторого тела (элемента почвы или сеянца); i - номер шага интегрирования (i - текущий шаг, i + 1 - последующий шаг); Дt - шаг интегрирования.

Шаг численного интегрирования Дt системы дифференциальных уравнений определяется путем многократного проведения экспериментов с последовательно уменьшающимся в 2 раза шагом. Останавливаются на том шаге, после которого результаты моделирования практически не изменяются (изменение составляет не более 1-2 %). Определенный таким образом шаг составил Дt = 5·10-⁴ с и был использован во всех расчетах при проведении теоретических исследований.

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

Рисунок 5 - Схема алгоритма моделирования работы грунтомета

Заключение

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

Работа выполнена при финансовой поддержке Министерства образования и науки Российской Федерации по государственному соглашению № 14.В37.21.2095 в рамках ФЦП «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы»

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

1. Советов Б.Я., Яковлев С.А. Моделирование систем: учебное пособие. М.: Высшая школа, 1998. 319 с.

2. Малиновский Е.Ю. Расчет и проектирование строительных и дорожных машин на ЭВМ. М.: Машиностроение, 1980. 216 с.

3. Синеоков Г.Н., Панов И.М. Теория и расчет почвообрабатывающих машин. М.: Машиностроение, 1977. 328 с.

4. Лурье А.Б. Моделирование сельскохозяйственных агрегатов и их систем управления: Учеб. для вузов. Л.: Колос. Ленингр. отд-ние, 1979. 312 с.

5. Троицкий В.А. Инженерные расчеты на ЭВМ: Справочное пособие. Л.: Машиностроение, 1979. 288 с.

Аннотации

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА ПОДАЧИ И ВЫБРОСА ГРУНТА РАБОЧИМИ ОРГАНАМИ КОМБИНИРОВАННОЙ МАШИНЫ ДЛЯ ТУШЕНИЯ ЛЕСНЫХ ПОЖАРОВ

Драпалюк Михаил Валентинович, д.т.н., профессор

Бартенев Иван Михайлович, д.т.н., профессор

Гнусов Максим Александрович, аспирант

Дручинин Денис Юрьевич, ассистент

Воронежская государственная лесотехническая академия, Воронеж, Россия

Марков Олег Борисович,, к.т.н., доцент

Петрозаводский государственный университет, Петрозаводск, Россия

Клубничкин Евгений Евгеньевич, к.т.н, доцент

Московский государственный университет леса, Мытищи, Россия

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

Ключевые слова: ЛЕСНОЙ ПОЖАР, ГРУНТОМЕТ, МАТЕМАТИЧЕСКАЯ МОДЕЛЬ, ГРУНТ, РАБОЧИЕ ОРГАНЫ

MATHEMATICAL MODEL OF GIVING AND SOIL EMISSION PROCESS WITH THE WORKING TOOLS OF THE COMBINED MACHINE FOR FOREST FIRES SUPPRESSION

Drapalyuk Mikhail Valentinovich, Dr.Sci.Tech., professor

Bartenev Ivan Mikhaylovich, Dr.Sci.Tech., professor

Gnusov Maksim Aleksandrovich, postgraduate student

Druchinin Denis Yurievich, assistant

Voronezh State Forestry Academy, Voronezh, Russia

Markov Oleg Borisovich, Cand.Tech.Sci., associate professor

Petrozavodsk State University, Petrozavodsk, Russia

Klubnichkin Evgeniy Evgeniyevich, Cand.Tech.Sci., associate professor

Moscow State Forest University, Mytischi, Russia

The mathematical model of the process of giving and soil emission with the working tools of a combined soil-thrower is developed for optimization of various parameters of the developed assembly

Keywords: FOREST FIRE, SOIL-THROWER, MATHEMATICAL MODEL, SOIL, WORKING TOOL

Размещено на Allbest.ru