Научная статья на тему 'ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ КИНЕМАТИЧЕСКИХ И ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ШНЕКОВОГО РАБОЧЕГО ОРГАНА ЛЕСОПОЖАРНОЙ ГРУНТОМЕТАТЕЛЬНОЙ МАШИНЫ'

ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ КИНЕМАТИЧЕСКИХ И ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ШНЕКОВОГО РАБОЧЕГО ОРГАНА ЛЕСОПОЖАРНОЙ ГРУНТОМЕТАТЕЛЬНОЙ МАШИНЫ Текст научной статьи по специальности «Механика и машиностроение»

CC BY
56
13
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Лесотехнический журнал
ВАК
AGRIS
RSCI
Ключевые слова
МОДЕЛИРОВАНИЕ / НАПОЧВЕННЫЙ ПОКРОВ / ШНЕКОВЫЙ РАБОЧИЙ ОРГАН / ГРУНТОМЕТАТЕЛЬНАЯ МАШИНА

Аннотация научной статьи по механике и машиностроению, автор научной работы — Попиков Петр Иванович, Поздняков Антон Константинович, Усков Владимир Игоревич, Лысыч Михаил Николаевич, Гнусов Максим Александрович

В работе представлены результаты исследований рабочих процессов шнековых рабочих органов в различных отраслях производства. Установлено, что исследований технологических параметров шнековых рабочих органов при удалении напочвенного покрова с горючими растительными остатками проведено недостаточно. Предложена новая конструкция лесопожарной машины, включающая шнековые рабочие органы для удаления напочвенного покрова с горючими растительными остатками, вырезные диски для образования почвенного вала перед ротором-метателем, регулируемый кожух-направитель с возможностью изменять направленный поток грунта. Рабочий процесс гидропривода шнековых рабочих органов лесопожарной грунтометательной машины описан системой дифференциальных уравнений поступательного, вращательного движений шнекового рабочего органа и уравнением расходов рабочей жидкости. Рассмотрена задача Коши на интервале t∈(t_0;t_end ). В силу нелинейности исходная система сводится к системе рекуррентных соотношений, заменяя производные искомых функций их конечными аналогами. Эта система неразрешенная относительно старшего члена. Благодаря свойству фредгольмовости операторного коэффициента перед старшим членом можно применить метод расщепления системы на системы в подпространствах уменьшающихся размерностей и решить систему в этих подпространствах. Аналитическое решение математической модели удаления напочвенного покрова шнековыми рабочими органами лесопожарной грунтометательной машины при прокладке противопожарных полос и тушении низовых лесных пожаров позволяет по конечным формулам определить кинематические и динамические параметры шнековых рабочих органов. Проведено имитационное моделирование рабочего процесса шнека на лабораторном виртуальном стенде. Получены зависимости динамических характеристик шнековых рабочих органов при преодолении препятствий. Установлено, что максимальное среднее значение крутящего момента наблюдается при преодолении препятствия высотой 100 мм и скорости вращения рабочего органа 400 об/мин и составляет 1468,49 Н∙м, при котором не нарушается его работоспособность

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Попиков Петр Иванович, Поздняков Антон Константинович, Усков Владимир Игоревич, Лысыч Михаил Николаевич, Гнусов Максим Александрович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

THEORETICAL STUDY OF THE KINEMATIC AND DYNAMIC CHARACTERISTICS OF THE AUGER WORKING BODY OF A FOREST FIRE GROUND-THROWING MACHINE

The paper presents the results of studies of the working processes of screw working bodies in various industries. Studies of the technological parameters of screw working bodies during the removal of the ground cover with combustible plant residues have been insufficiently carried out. A new design of a forest fire machine is proposed, including screw working bodies for removing the ground cover with combustible plant residues, cut-out discs, for the formation of a soil shaft in front of the rotor-thrower, with an adjustable guide casing with the ability to change the directional flow of soil. Working bodies of the forest fire soil-throwing machine are described by a system of differential equations of the translational, rotational movements of the auger working body and the equation of the flow rates of the working fluid. The Cauchy problem is considered on the interval t∈ (t_0; t_end). Original system is reduced to a system of recurrence relations due to its nonlinearity, replacing the derivatives of the desired functions with their finite analogs. This system is unresolved with respect to the senior member. One can apply the method of splitting the system into systems in subspaces of decreasing dimensions and solve the system in these subspaces due to the Fredholm property of the operator coefficient in front of the leading term. The analytical solution of the mathematical model for removing the ground cover by auger working bodies of a forest fire soil-throwing machine when laying fire strips and extinguishing grassland forest fires makes it possible to determine the kinematic and dynamic parameters of the auger working bodies using the final formulas. Simulation modeling of the screw working process was carried out on a virtual laboratory bench. The dependences of the dynamic characteristics of the auger working bodies when overcoming obstacles have been obtained. It was found that the maximum average value of the torque is observed when overcoming an obstacle with a height of 100 mm and the rotation speed of the working body of 400 rpm is 1468.49 N ∙ m, at which its performance is not impaired.

Текст научной работы на тему «ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ КИНЕМАТИЧЕСКИХ И ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ШНЕКОВОГО РАБОЧЕГО ОРГАНА ЛЕСОПОЖАРНОЙ ГРУНТОМЕТАТЕЛЬНОЙ МАШИНЫ»

DOI: 10.34220^п.2222-7962/2021.3/12 УДК 630*232.318:630*651.72

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

Петр И. Попиков1, popikovpetr@yandex.ru, 0000-0002-6348-8934

Антон К. Поздняков1 И, mcak74@gmail.com, 0000-0003-4590-695Х

Владимир И. Усков1, уитШуаМех.ги, 0000-0002-3542-9662

Михаил Н. Лысыч1, vumli@yandex.ru, 0000-0002-3764-3873

Максим А. Гнусов1, ko407@yandex.ru, 0000-0003-1653-4595

ФГБОУ ВО «Воронежский государственный лесотехнический университет имени Г. Ф. Морозова», ул. Тимирязева, 8, г. Воронеж, 394087, Россия

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

производные искомых функций их конечными аналогами. Эта система неразрешенная относительно старшего члена. Благодаря свойству фредгольмовости операторного коэффициента перед старшим членом можно применить метод расщепления системы на системы в подпространствах уменьшающихся размерностей и решить систему в этих подпространствах. Аналитическое решение математической модели удаления напочвенного покрова шнековыми рабочими органами лесопожарной грунтометательной машины при прокладке противопожарных полос и тушении низовых лесных пожаров позволяет по конечным формулам определить кинематические и динамические параметры шнековых рабочих органов. Проведено имитационное моделирование рабочего процесса шнека на лабораторном виртуальном стенде. Получены зависимости динамических характеристик шнековых рабочих органов при преодолении препятствий. Установлено, что максимальное среднее значение крутящего момента наблюдается при преодолении препятствия высотой 100 мм и скорости вращения рабочего органа 400 об/мин и составляет 1468,49 Н-м, при котором не нарушается его работоспособность.

Ключевые слова: моделирование, напочвенный покров, шнековый рабочий орган, грунтометательная машина

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

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 20-38-90082 "Аспиранты "

Авторы благодарят рецензентов за вклад в экспертную оценку статьи.

Конфликт интересов: авторы заявили об отсутствии конфликта интересов.

Для цитирования: Математическое моделирование удаления напочвенного покрова шнековыми рабочими органами лесопожарной грунтометательной машины / П. И. Попиков, А. К. Поздняков, В. И. Усков, М. Н. Лысыч, М. А. Гнусов // Лесотехнический журнал. - 2021. - Т. 11. - № 3 (43). - С. 140-151. - Библиогр.: с. 147-150 (16 назв.). - DOI: https://doi.Org/10.34220/issn.2222-7962/2021.3/12.

Поступила: 30.06.2021 Принята к публикации: 25.09.2021 Опубликована онлайн: 01.07.2021

THEORETICAL STUDY OF THE KINEMATIC AND DYNAMIC CHARACTERISTICS OF THE AUGER WORKING BODY OF A FOREST FIRE GROUND-THROWING MACHINE

Petr I. Popikov1, popikovpetr@yandex.ru, 0000-0002-6348-8934 Anton K. Pozdnyakov El, mcak74@gmail.com, 0000-0003-4590-695X Vladimir I. Uskov1, vum1@yandex.ru, 0000-0002-3542-9662 Mikhail N. Lysych\ vmnl i@yandex.ru, 0000-0002-3764-3873 Maksim A. Gnusov1, ko407@yandex.ru, 0000-0003-1653-4595

'Voronezh State University of Forestry and Technologies named after G.F. Morozov, Timiryazeva str., 8, Voronezh, 394087, Russian Federation

Abstract

The paper presents the results of studies of the working processes of screw working bodies in various industries. Studies of the technological parameters of screw working bodies during the removal of the ground cover with combustible plant residues have been insufficiently carried out. A new design of a forest fire machine is proposed, including screw working bodies for removing the ground cover with combustible plant residues, cut-out discs, for the formation of a soil shaft in front of the rotor-thrower, with an adjustable guide casing with the ability to change the directional flow of soil. Working bodies of the forest fire soil-throwing machine are described by a system of differential equations of the translational, rotational movements of the auger working body and the equation of the flow rates of the working fluid. The Cauchy problem is considered on the interval te (t_0; t_end). Original system is reduced to a system of recurrence relations due to its nonlinearity, replacing the derivatives of the desired functions with their finite analogs. This system is unresolved with respect to the senior member. One can apply the method of splitting the system into systems in subspaces of decreasing dimensions and solve the system in these subspaces due to the Fredholm property of the operator coefficient in front of the leading term. The analytical solution of the mathematical model for removing the ground cover by auger working bodies of a forest fire soil-throwing machine when laying fire strips and extinguishing grassland forest fires makes it possible to determine the kinematic and dynamic parameters of the auger working bodies using the final formulas. Simulation modeling of the screw working process was carried out on a virtual laboratory bench. The dependences of the dynamic characteristics of the auger working bodies when overcoming obstacles have been obtained. It was found that the maximum average value of the torque is observed when overcoming an obstacle with a height of 100 mm and the rotation speed of the working body of 400 rpm is 1468.49 N • m, at which its performance is not impaired.

Keywords: modeling, ground cover, screw working body, soil throwing machine

Acknowledgments: The authors express their gratitude to colleagues, Candidate of Technical Sciences of the Department of Forestry Mechanization and Machine Design, Head of the Laboratory Gnusov M.A. and the junior researcher Petkov A.F., who advised the work in the process of writing it.

The research was carried out with the financial support of the RFBR in the framework of the scientific project No. 20-38-90082 "Postgraduates"

The authors thank the reviewers for their contribution to the expert evaluation of the article.

Conflict of interest: the author(s) declare(s) no conflict of interest.

For citation: Popikov P. I., Pozdnyakov A. K., Uskov V. I., Lysych M. N., Gnusov M. A. (2021) Mathematical modeling of the removal of ground cover by auger working bodies of a forest fire ground-throwing machine. Lesotekhnicheskii zhurnal [Forest Engineering Journal]. Vol. 11. No. 3 (43). P. 140-151. DOI: https://doi.org/10.34220/issn.2222-7962/202L3/12.

Received: 30.06.2021 Accepted for publication: 25.09.2021 Published online: 01.07.2021

Введение

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

В работе Лобанова Г. Л. [1] математическая модель взаимодействия стеблей растений с ножами шнековых рабочих органов включает уравнения проекций скоростей крайних точек ножа на оси координат в параметрическом виде. Определены численные значения параметров режущего аппарата в уравнения в программе Microsoft Excel и выявлена траектория движения ножа при срезании стеблей растений. Однако, предложенная математическая модель - это система уравнений первого порядка. Она не содержит производных второго порядка, т.е. не учитывает второй закон Ньютона.

В расчетах Хайбуллина Р.Р. [3] на основе анализа математической модели выкопки траншей под фундамент рабочим органом землеройной машины получены зависимости сопротивления грунта, скорости бурения, фрезерования, размыва грунта, а также мощности и энергоемкости. Следует отметить, что данная математическая модель является линейной по искомым функциям х,<р. Линей-

ная модель сильно огрубляет описание реального процесса: для угла <р, близкого к 0, в силу первого замечательного предела допустима замена Бткф на <рк, к = 0,1,2,..., но при больших значениях <р такая замена влечет накопление погрешности в вычислениях.

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

V, но =<ЛУ2, поэтому сводится к линейной, и функция V находится в явном виде. Кроме того шнек движется параллельно своей оси, а в составе лесопожарных грунтометательных машин ось шнеков расположена перпендикулярно направлению движения трактора..

В работе Закирова М.Ф. [2] проводились исследования влияния шага шнека на мощность привода питателя малогабаритного шнекороторного снегоочистителя. Шнековый рабочий орган выполнен в виде барабана диаметром 0,3 м, к которому приварена ленточная лопасть шириной 0,7 м и шагом 0,3 м с углом подъема лезвия шнека около 17,8°. Установлено, что мощность привода шнеко-вого питателя расходуется на вырезание снега из массива и на перемещение снега питателем.

В работах Бартенева И.М. и Nunez-Regueira. L. [4, 8] выполнено математическое моделирование рабочих процессов шнековых рабочих органов лесопожарной машины, с использованием метода конечных элементов, при котором лесной грунт представляется множеством (порядка 6000) шарообразных элементов. Определена оптимальная высота ленты шнекового барабана в пределах 8-10 см. При метании грунта доля напочвенного покрова составила 11-15 %, мощность на привод шнека 11-16 кВт. Однако кинематические и динамические характеристики шнековых рабочих органов авторами не изучались.

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

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

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

Материалы и методы

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

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

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

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

1) (Мш + МГр) ^ = РХ-[СШ + СТ р +Р2]/ +

71+/2 • с<ю(р + агсС^/) •(Рр + РТр) + 71+/2 • соб(ф + агсс

р±Т + РГС05<Р\ (Рр + Ртр)-[^ + Р/С08<р\, ЦМ ЙЮ йр

3) ^"Л+^Р + ^Л'

где Мш, Мгр - масса шнека и напочвенного покрова на шнеке, кг;

Gш, Gгр - силы тяжести шнека и напочвенного покрова на шнеке, Н;

f - коэффициент трения скольжения;

R - радиус шнека, м;

L - длина шнека, м;

Н - глубина удаляемого покрова (0,10...0,15), м;

t - время, с;

1пр - приведенный к валу гидромотора момент инерции, кгм2;

Ь - ширина рабочей грани черенкового ножа, м;

Ь - угол заточки черенкового ножа, град;

а - угол наклона винтовой лопасти от перпендикуляра к оси шнека, град;

ц> - угол поворота вала шнека, рад; Чн„ Чм - рабочие объемы насоса и гидромотора, м3/об;

р - давление рабочей жидкости, Па; ау - коэффициент утечек, м3/(сПа); пн - частота вращения насоса, с-1;

кр - коэффициент податливости упругих элементов гидропривода, м5/(Нс);

К - удельное сопротивление резанию черенкового ножа;

т - сопротивление сдвигу напочвенного покрова, Па;

фтр - угол трения почвы о винтовую поверхность шнека, град;

5 - шаг винтовой ленты, м.

Рисунок 1. Расчетная схема шнекового рабочего органа (собственные разработки). VTp - скорость трактора, Gш, Grp -силы тяжести шнека и грунта на шнеке, Pz - сила пружины предохранителя, Fтp - сила сопротивления перекатыванию шнека, N - реакция поверхности почвы, Мк - крутящий момент гидромотора, Fc - сила сопротивления резанию,

Px - тяговая сила, Ри - сила инерции Figure 1. Design scheme of the screw working body (own developments). Vtr is the tractor speed, Gsh, Ggr is the gravity force of the auger and the soil on the auger, Pz is the spring force of the safety device, Ftr is the resistance to rolling of the auger, N is the reaction of the soil surface, Mk is the torque of the hydraulic motor, Fc is the cutting resistance force, Px is traction force, Pi - force

of inertia

2.1. Постановка задачи

Рассматривается задача Коши на интервале tE (t0;tend):

d2x

= Pr-

(Mm + MTp)dt2 (Gm + GTp + Pz)f + Fc^1 + f2cos{y + arcctgf),

d2v _ Шр- p (R+r\

(.tend ~t0)/n. Выбираем шаг так, чтобы кф Кра}/1. В обозначениях х^ = х(^), ^ = ^(tj), pi = p(ti) имеем

d x 2 2 1 £

dF~ h2 '

fy^ <Pi+i~<Pi d2V Vi+2~2<Pi+i + Vi dt ~

'пр

dt¿

4 2/

+

(1)

h2 ' dt2

dp ^ pi+1

dt ~

h2

Pi

h2

+ СГр + PZ -Fccos<p)Rf, qMdm dp

q»n» = -—t+ayV + Kv-dt,

x(t0) = X0,X (t0) =X1,

P(t0) = <p0,<p'(t0) = (рг,

pito) = p0, (2)

Для решения задачи аппроксимируем производные конечными разностями, для чего отрезок tendí разобьём п узловыми точками t¿ с шагом h: ti+1 = t¿ + h; так, величина шага равна h =

Тогда систему можно переписать, выделив Х[+2 в первом уравнении и во втором и умножив на К2:

(Мш + МТр)х1+2 = 2{МШ + МГр)х1+1 ~(МШ +МТр)хь + +к2Рх-к2(Сш + СТр + Рг)Г +

+ f2cos{^i + arcctgf), InpVí+2 =2hp<Pí+l -InpVí +h2^Pl-h2Fc(^) +

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(3)

к2Н[(вш + Сгр + Р2-Рссозф1), п Чм Чм . ,,

0=—91+1--Фь+КрР1+1

л л

+ (кау-Кр)р1-кЦнП11. Рассмотрим отдельно второе и третье соотношения системы. Их можно записать в векторном виде

АУ1+2 =Ву1+1 + Р1'

(4)

с операторами А, В-.И2 ^И2

и вектор-последовательностями

(6)

в обозначениях:

+к2ЯГ(Сш + Сгр + Р2-Рссо5<р1),

(2) Чм / \

^ + (]гау -Кр)р1 -КЦнПн.

Необходимые сведения

Далее, рассматривается векторное линейное рекуррентное соотношение (далее, ЛРС) второго порядка

Лу1+2=Ву1+1+Р1' 1 = 0,1,2..........(2.1)

с заданными начальными значениями у0, у±, где А, В - операторы, задаваемые квадратными числовыми матрицами, причем det А = 0.

При разрешении ЛРС (2.1) используется факт, что линейный оператор, задаваемый вырожденной квадратной матрицей, фредгольмов [8].

Пусть оператор А обладает одномерным ядром. Для него вводится проектор Q в подпространство Сокег А, полуобратный оператор А~ = - Q), где А = лпйот А; фиксируются элементы еЕ Кег А, еф 0, фЕ Сокег А.

Вводится в Сокег А скалярное произведение <, > так, что <ф,(р >= 1.

Используется следующий результат [9]. Лемма 1 Линейное уравнение

Ар = ж, V Е Е± п dom А, ж Е Е2, равносильно системе

у = А~ш + с • е для любых с ЕС, <Qw,(p >=0.

Пусть выполнено следующее условие.

Теорема 1 Пусть выполнены условия

1. Выражение <QFi,<p > определено при всех

I = 0,1,....

йез

2. й = <QBe,(p >Ф 0.

Тогда ЛРС (3.2.1)равносильно системе У1+2 =Ку1+1 + Ф>, (2.2)

<(}Ву1+1,<р > +< (¿ВРь<р ; (2.3)

в обозначениях К(•) = А~В{•) -а-1 <QBA~B(•),ф > е, Ф; = А~Р1 -й~1(<дВА-р1,ф > +

<QFi+1,<P>)e. Доказательство. В силу леммы 3.2.1 соотношение (3.2.1) равносильно системе

Уь+2 = А~Ву1+1 + А~Р, + С1е, (2.4)

и равенству (3.2.3), где с^ - искомая последовательность. В равенстве (3.2.3) заменим I на ¿ + 1: <QByi+2,ф > +< др1+1,ф >= 0 и вместо у 1+2 подставим выражение (3.2.4):

<QBA~Byi+1,(p >+<QBA-Fi,<p > +С1 <QBe, ср> + +<QFi+1,Ф>=0. (2.5)

Из (3.2.5) выразив С; и подставив в (3.2.4), получим соотношение (2.2). Теорема доказана.

Отметим, что равенство (3.2.3) при 1 = 0 -это условие разрешимости начальной задачи для соотношения (3.2.1):

<дВу1,ср >+<дВР0,^ >= 0. (2.6)

Далее рассматривается скалярное ЛРС второго порядка

и1+2 + ащ+1 + Ьщ = I = 0,1,2,..., (2.7)

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

Л2+аЛ + Ь = 0 (2.8)

п 2 м 1 -а+4о

и = а —4Ь -дискриминант, Я12 = —-— -

корни.

Лемма 2 Пусть 0> 0. Тогда общее решение и^ соотношения (3.2.7)равно

Щ =с1(Л1у+ с2(Л2у 1-2

^ 1=0 = 2,3,..,

где с1, с2 - произвольные постоянные.

Лемма 3 Пусть 0 = 0, Л - корень уравнения (3.2.8). Тогда общее решение щ соотношения (3.2.7) равно

и^ = с1Л1 + с21Л1~1

1-2 1=0

-1)Л1-2~1[1, 1 = 2,3,... Обе леммы доказываются непосредственной подстановкой решения в (3.2.7). Об операторе

А = (1*Р 0). V 0 0/

Результаты и обсуждение Решение начальной задачи для выражений (4), (5), (6)

Для вычисления X; получено соотношение (3), откуда

Х1+2 = 2Х1 +1 ~Х1 ^ Ь '

(Мш +Мгр)~1Рх

2 г л-1 (7)

~к2(Мш +Мгр) •

+h2(Mm + мгру FcVT+я-

cos{(pi + arcctg/). Имеем:

d = (QBe, ф) = КрФ0. Вычисления приводит к системе:

т I т 2 j—1 Чм

ф1+2 = 2cpi+1 -(Pi + Vl^ — Vi

, R+r ~h2Inp-Fc—2--1"

+h2I~v1Rf

+ СГр + PZ — Fccos(pj)

Pi+2 = (T - hKp1ay)pi+1 -

+ hK~1qHnH.

(8)

(9)

Из разложения функций х,ф в ряд Тейлора в окрестности точки £ = £0 при И.^0

х{Ъ0 + Ь) = х{1с0) + х' + 0(И2), ф(р0 + =x(to) + x'(to)ft + ОСЬ2),

и начальных значений (2) вытекают формулы для нахождения значений х1, ф±: х1 =х° +х1И,

(10)

фг = ф° + ф1Ь.

(11)

Далее, из соотношения (10) вытекает формула для вычисления р±:

р1 = (1- КК^1ау)р0 -

К"1 f 0^-0 +

^ П

hK^1qHnH.

(12)

Теорема 2.

1. Решение хь с начальными значениями х0 = х0 и хх, определяемым формулой (10), равно XI = х0 +х1Ш + +й2 Мш +Мгр

Y(i-i-T)[px

1=0

(13)

-(Сш + Сгр + рг)г +

+71 + рРсСОБ (ф1 + аГсс шо],

I = 2,3,...

2. Значения угла поворота шнека ф^ определяются из итерационного процесса (8), где ф0 = Ф0, а ф± из (11).

3. Решение для определения давления рабочей жидкости в гидромоторе Р1 с начальными значениями р0 = р° и р±, определяемым формулой (12), равно

-(1 - нк-^уУ1 +(1 - нк-^уУ1

i-1

i-1-l

Yp--hKviayT (14)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1=0

(-Kp^Vi+i + hK~1qHnHy

Тем самым решение x(t), ф(р), p(t) задачи (1), (2) определяется приближенно в виде значений xi, <Pi, Pi в узловых точках tt

При участии доц. Лысыча М.Н. проведено имитационное моделирование рабочего процесса

шнековых рабочих органов. На рис. 2 представлены зависимости изменения средних значений крутящего момента от скорости вращения шнекового рабочего органа.

С увеличением высоты препятствия крутящий момент заметно возрастает. При этом интенсивность роста увеличивается с увеличением высоты препятствия. Так среднее минимальное значение крутящего момента для высоты препятствия 0 мм и скорости вращения рабочего органа 200 об/мин составляет 880,20 Н-м. Максимальное среднее значение крутящего момента наблюдается для высоты препятствия 100 мм и скорости вращения рабочего органа 400 об/мин, это 1468,49 Н-м.

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

Заключение

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

Рисунок 3. Зависимости изменения средних значений

крутящего момента от высоты препятствия [9] Figure 3. Dependences of the change in the average values on the height of the obstacle [9]

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

1. Лобанов Г. Л., Погоров Т. А. Математическая модель траектории движения ножей шнекового режущего аппарата в зоне резания стеблей растений. Научный журнал Российского НИИ проблем мелиорации.

Рисунок 2. Зависимости изменения средних значений крутящего момента от скорости вращения шнекового

рабочего органа [11] Figure 2. Dependences of changes in the average values on the rotational speed of the auger working body [11]

Анализ полученных зависимостей показал, что крутящий момент гидромотора незначительно возрастает с увеличением скорости вращения шне-кового рабочего органа практически по линейному закону. При возрастании частоты вращения от 200 об/мин до 400 об/мин крутящий момент гидромотора увеличивается от 1400 Нм до 1490 Нм.

Характер изменения средних значений крутящего момента от высоты препятствия, преодолеваемого шнеком, представлен на рис. 3.

2017; 1(25). URL: https://cyberleninka.ru/article/n/matematicheskaya-model-traektorii-dvizheniya-nozhey-shnekovo go-rezhuschego-apparata-v-zone-rezaniya-stebley-rasteniy (дата обращения: 14.05.2021).

2. Закиров М. Ф. Исследование влияния шага шнека на мощность привода питателя малогабаритного шнекороторного снегоочистителя. Интеллектуальные системы в производстве. 2015;2: 56-57. URL: https://www.elibrary.ru/item.asp?id=40999966/ (дата обращения: 27.04.2021).

3. Хайбуллин P.P. Теоретические основы разрушения грунта при фрезеровании и бурении рабочими органами строительных машин: автореф. дис. ... д-ра техн. наук : 05.05.04. Москва, 2008. URL: http://www.scholar.ru/speciality.php?spec_id=276 (дата обращения: 15.05.2021).

4. Bartenev I. М., Popikov P. I., Malyukov S. V. Research and development of the method of soil formation and delivery in the form of a concentrated flow to the edge of moving ground forest fire. IOP Conference Series: Earth and Environmental. 2019;226(1): 012052. DOI: 10.1088/1755-1315/226/1/012052.

5. Согин А. В. Средства гидромеханизиции для очистки водных объектов от донных отложений: дис. ... д-ра техн. наук. Москва, 2011. 360 с. URL: https://www.dissercat.com/content/sredstva-gidromekhanizatsii-dlya-ochistki-vodnykh-obektov-ot-donnykh-otlozhenii (дата обращения: 16.05.2021).

6. Бартенев И. М., Поздняков А. К. Анализ рабочих органов технических средств для тушения лесных пожаров. Актуальные направления научных исследований XXI века: теория и практика. Воронеж, 2020. № 1 (48). С. 119-122. DOI: 10.34220/2308-8877-2020-8-1-119-122.

7. Whitehouse N. J. Forest fires and insects: palaeoentomological research from a subfossil burnt forest. Palaeogeography, Palaeoclimatology, Palaeoecology. 2000;164(1-4): 231-246. URL: https://www.sciencedirect.com/ science/article/pii/S0031018200001887 (дата обращения: 15.05.2021).

8. Nunez-Regueira L. [et al.] Calculation of forest biomass indices as a tool to fight forest fires. Thermochimica Acta. 2001;378(1-2): 9-25. DOI: 10.15372/SJFS20150609.

9. Бартенев И. М. (и др.) Комбинированный лесопожарный грунтомет и рекомендации по его применению. Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. 2012;84. URL: https://cyberleninka.ru/article/n/kombinirovannyy-lesopozharnyy-gruntomet-i-rekomendatsii-po-ego-primeneniyu (дата обращения: 16.05.2021).

10. Патент № 2684940 Российская Федерация, МПК Е02 F 3/18. Пожарный грунтомет-полосопрокладыватель : № 2016104672; заявл. 09.07.2018; опубл. 16.04.2019 / И. М. Бартенев, П. И. Попиков, С. В. Малюков, С. В. Зимарин, Н. А. Шерстюков; заявитель и патентообладатель ФГБОУ ВО «ВГЛТУ им. Г.Ф. Морозова». URL: https://patents.google.com/patent/RU2541987C1/ru (дата обращения: 01.02.2021).

11. Драпалюк, М. В. (и др.) Математическая модель процесса подачи и выброса грунта рабочими органами комбинированной машины для тушения лесных пожаров. Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. 2012;84. URL: https://cyberleninka.ru/article/n/matematicheskaya-model-protsessa-podachi-i-vybrosa-grunta-rabochimi-organami-kombinirovannoy-mashiny-dlya-tusheniya-lesnyh-pozharov (дата обращения: 12.03.2021).

12. Гнусов М. А., Попиков П. И., Малюков С. В., Шерстюков Н. А., Поздняков А. К. Повышение эффективности предупреждения и тушения лесных пожаров с помощью лесопожарной машины. IOP: Материаловедение и инженерия. 2020;919(3): 032025. URL: https://cyberleninka.ru/article/n/vliyanie-rezhimov-raboty-lesopozharnoy-gruntometatelnoy-mashiny-s-gidroprivodom-na-pokazateli-effektivnosti (дата обращения: 22.04.2021).

13. Kruggel-Emden H., Wirtz S., Scherer V. (2008) A study on tangential force laws applicable to the discrete element method (DEM) for materials with viscoelastic or plastic behavior. Chemical Engineering Science, 63(6), 15231541. URL: https://www.edemsimulation.com/papers/a-study-on-tangential-force-laws-applicable-to-the-discrete-element-method-dem-for-materials-with-viscoelastic-or-plastic-behavior/ (дата обращения: 27.05.2021).

14. Tsunazawa Y., Shigeto Y., Tokoro C., Sakai, M. (2015) Numerical simulation of industrial die filling using the discrete element method. Chemical engineering science, 138, 791-809. DOI: 10.14419/ijet.v7i2.23.11876.

15. Guo Y., Curtis J. S. (2015) Discrete element method simulations for complex granular flows. Annual Review of Fluid Mechanics, 47, 21-46. URL: https://www.annualreviews.org/doi/abs/10.1146/annurev-fluid-010814-014644.

16. Орловский С. Н. Определение энергетических и динамических параметров тракторов, режимов резания активных рабочих органов машинно-тракторных агрегатов: моногр. Красноярск: КрасГАУ, 2011. 376 с. URL: http://lesnoizhurnal.ru/issuesarchive/ ELEMENT_ID=294450 (дата обращения: 21.05.2021).

References

1. Lobanov G. L., Pogorov T. A. Mathematical model of the trajectory of the movement of knives of the screw cutting device in the cutting zone of plant stems. Scientific Journal of the Russian Research Institute of Problems of Melioration. 2017;1(25). URL: https://cyberleninka.ru/article/n/matematicheskaya-model-traektorii-dvizheniya-nozhey-shneko- vo go-rezhuschego-apparata-v-zone-rezaniya-stebley-rasteniy (date of access: 14.05.2021) (in Russian).

2. Zakirov M. F. Investigation of the effect of the auger pitch on the power of the feeder drive of a small-sized auger-rotor snowplow. Intelligent systems in production. 2015;2: 56-57. URL: https://www.elibrary.ru/ item.asp?id=40999966/ (date of access: 27.04.2021) (in Russian).

3. Khaibullin R. R. Theoretical foundations of soil destruction during milling and drilling by working bodies of construction machines: abstract. ... doct. tekhn. science: 05.05.04 [D.Sci. thesis abstr.]. Moscow, 2008. URL: http://www.scholar.ru/speciality.php?spec_id=276 (date of access: 11.05.2021) (in Russian).

4. Bartenev I. M., Popikov P. I., Malyukov S. V. Research and development of the method of soil formation and delivery in the form of a concentrated flow to the edge of moving ground forest fire. IOP Conference Series: Earth and Environmental. 2019;226(1): 012052. DOI: 10.1088/1755-1315/226/1/012052.

5. Sogin A. V. Means of hydro-mechanization for cleaning water bodies from bottom sediments: diss. ... Doctor of technical Sciences. Moscow, 2011. 360 p. URL: https://www.dissercat.com/content/sredstva-gidromekhanizatsii-dlya-ochistki-vodnykh-obektov-ot-donnykh-otlozhenii (date of access: 16.05.2021) (in Russian).

6. Bartenev I. M., Pozdnyakov A. K. Analysis of working bodies of technical means for extinguishing forest fires / I. M. Bartenev. Actual directions of scientific research of the XXI century: theory and practice. Voronezh, 2020. №1 (48), pp. 119-122. DOI: 10.34220/2308-8877-2020-8-1-119-122 (in Russian).

7. Whitehouse N. J. Forest fires and insects: palaeoentomological research from a subfossil burnt forest. Palaeogeography, Palaeoclimatology, Palaeoecology. 2000;164(1-4): 231-246. URL: https://www.sciencedirect.com/ science/article/pii/S0031018200001887 (дата обращения: 15.05.2021).

8. Nunez-Regueira L. [et al.] Calculation of forest biomass indices as a tool to fight forest fires. Thermochimica Acta. 2001;378(1-2): 9-25. DOI: 10.15372/SJFS20150609.

9. Bartenev I. M. (et al.) Kombinirovannyy lesopozharnyj gruntomet i rekomendatsii po ego primenenii [Combined forest fire gruntomet and recommendations for its use]. Polythematic network electronic scientific journal of the Kuban State Agrarian University. 2012;84. URL: https://cyberleninka.ru/article/n/kombinirovannyy-lesopozharnyy-gruntomet-i-rekomendatsii-po-ego-primeneniyu (date of access: 16.05.2021) (in Russian).

10. Patent No. 2684940 Russian Federation, IPC E02 F 3/18. Fire grundomat-protopalatial : No. 2016104672; Appl. 09.07.2018; publ. 16.04.2019 / I. M. Bartenev, P. I. Popikov, S. V. Malyukov, S. V. Zimarin, N. A. Sherstyukov; applicant and patentee of the "VGLTU named after. G.F. Morozov". URL: https://patents.google.com/ patent/RU2541987C1/ru (date of the application: 01.02.2021).

11. Drapalyuk M. V. (et al.) Mathematical model of the process of feeding and ejecting soil by working bodies of a combined machine for extinguishing forest fires. Polythematic network electronic scientific journal of the Kuban State Agrarian University. 2012; 84. URL: https://cyberleninka.ru/article/n/matematicheskaya-model-protsessa-podachi-i-vybrosa-grunta-rabochimi-organami-kombinirovannoy-mashiny-dlya-tusheniya-lesnyh-pozharov (date of access: 12.03.2021) (in Russian).

12. Gnusov M. A., Popikov P. I., Malyukov S. V., Sherstyukov N. A., Pozdnyakov A. K. Improving the effectiveness of prevention and extinguishing of forest fires with the help of a forest fire engine. IOP Conference Series: Materials Science and Engineering. 2020;919(3): 032025. URL: https://cyberleninka.ru/article/n/vliyanie-rezhimov-raboty-lesopozharnoy-gruntometatelnoy-mashiny-s-gidroprivodom-na-pokazateli-effektivnosti (date of access: 22.04.2021).

13. Kruggel-Emden H., Wirtz S., Scherer V. (2008) A study on tangential force laws applicable to the discrete element method (DEM) for materials with viscoelastic or plastic behavior. Chemical Engineering Science, 63(6), 15231541. URL: https://www.edemsimulation.com/papers/a-study-on-tangential-force-laws-applicable-to-the-discrete-element-method-dem-for-materials-with-viscoelastic-or-plastic-behavior/ (date of access: 27.05.2021).

14. Tsunazawa Y., Shigeto Y., Tokoro C., Sakai, M. (2015) Numerical simulation of industrial die filling using the discrete element method. Chemical engineering science, 138, 791-809. DOI: 10.14419/ijet.v7i2.23.11876.

15. Guo Y., Curtis J. S. (2015) Discrete element method simulations for complex granular flows. Annual Review of Fluid Mechanics, 47, 21-46. URL: https://www.annualreviews.org/doi/abs/10.1146/annurev-fluid-010814-014644.

16. Orlovsky S. N. Determination of energy and dynamic parameters of tractors, cutting modes of active working bodies of machine-tractor units: monogr. Krasnoyarsk, KrasGAU, 2011. 376 p. URL: http://lesnoizhurnal.ru/issuesarchive/?ELEMENT_ID=294450 (date of access 21.05.2021) (in Russian).

Сведения об авторах

Попиков Петр Иванович - доктор технических наук, профессор, профессор кафедры механизации лесного хозяйства и проектирования машин ФГБОУ ВО «Воронежский государственный лесотехнический университет имени Г.Ф. Морозова», ул. Тимирязева, д. 8, г. Воронеж, Российская Федерация, 394087, ORCID: https://orcid.org/0000-0002-6348-8934, e-mail: popikovpetr@yandex.ru

ЕЕЗ Поздняков Антон Константинович - аспирант кафедры механизации лесного хозяйства и проектирования машин ФГБОУ ВО «Воронежский государственный лесотехнический университет имени Г.Ф. Морозова», ул. Тимирязева, д. 8, г. Воронеж, Российская Федерация, 394087, ORCID: https://orcid.org/0000-0003-4590-695X, e-mail: mcak74@gmail.com

Усков Владимир Игоревич - кандидат физико-математических наук, старший преподаватель ФГБОУ ВО «Воронежский государственный лесотехнический университет имени Г.Ф. Морозова», ул. Тимирязева, д. 8, г. Воронеж, Российская Федерация, 394087, ORCID: https://orcid.org/0000-0002-3542-9662, e-mail: vum1@yandex.ru

Лысыч Михаил Николаевич - кандидат технических наук, доцент кафедры механизации лесного хозяйства и проектирования машин ФГБОУ ВО «Воронежский государственный лесотехнический университет имени Г.Ф. Морозова», ул. Тимирязева, д. 8, г. Воронеж, Российская Федерация, 394087, ORCID: https:/ 0000-0002-3764-3873, e-mail: miklynea@yandex.ru

Гнусов Максим Александрович - кандидат технических наук кафедры механизации лесного хозяйства и проектирования машин, заведующий лабораторией ФГБОУ ВО «Воронежский государственный лесотехнический университет имени Г.Ф. Морозова», ул. Тимирязева, д. 8, г. Воронеж, Российская Федерация, 394087, ORCID: https:/ 0000-0003-1653-4595, e-mail: ko407@yandex.ru

Information about the authors

Popikov Petr Ivanovich - DSc (Engineering), Professor, Professor of the Department of Forestry Mechanization and Machine Design, Voronezh State University of Forestry and Technologies named after G. F. Morozov, Timiryazev str., 8 Voronezh, Russian Federation, 394087, ORCID: https://orcid.org/0000-0002-6348-8934, e-mail: popikovpetr@yandex.ru

ISI Pozdnyakov Anton Konstantinovich - post-graduate student of the Department of Forestry Mechanization and Machine Design, Voronezh State University of Forestry and Technologies named after G. F. Morozov, Timiryazev

str., 8 Voronezh, Russian Federation, 394087, ORCID: https://orcid.org/0000-0003-4590-695X, e-mail: mcak74@gmail.com.

Uskov Vladimir Igorevich - Candidate of Physical and Mathematical Sciences, senior lecturer, Voronezh State University of Forestry and Technologies named after G. F. Morozov, Timiryazev str., 8, Voronezh, Russian Federation, 394087, ORCID: https://orcid.org/0000-0002-3542-9662, e-mail: vum1@yandex.ru.

Lysych Mikhail Nikolaevich - Candidate of Technical Sciences, Associate Professor of the Department of Forestry Mechanization and Machine Design, Voronezh State University of Forestry and Technologies named after G.F. Morozov, Timiryazev str., 8, Voronezh, Russian Federation, 394087, ORCID: https:/ 0000-0002-3764-3873, e-mail: miklynea@yandex.ru.

Gnusov Maksim Aleksandrovich - Candidate of Technical Sciences, Department of Forestry Mechanization and Machine Design, Head of the Laboratory, Voronezh State University of Forestry and Technologies named after G.F. Morozov, Timiryazev str., 8, Voronezh, Russian Federation, 394087, ORCID: https:/ 0000-0003-1653-4595, e-mail: ko407@yandex.ru.

И- Для контактов/Corresponding author

i Надоели баннеры? Вы всегда можете отключить рекламу.