Математическая модель воздействия движителя на уплотняющийся лесной почвогрунт Текст научной статьи по специальности «Строительство и архитектура»

_ТЕХНИЧЕСКИЕ НАУКИ / ENGINEERING_

DOI: https://doi.org/10.23670/IRJ.2019.80.2.004

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВОЗДЕЙСТВИЯ ДВИЖИТЕЛЯ НА УПЛОТНЯЮЩИЙСЯ ЛЕСНОЙ

ПОЧВОГРУНТ

Научная статья

Андронов А.В.1' *, Котенев Е.В.2, Пегов В.Ю.3

1, 2 3 Санкт-Петербургский государственный лесотехнический университет имени С.М. Кирова,

Санкт-Петербург, Россия

* Корреспондирующий автор (andronovalexandr[at]gmail.com)

Аннотация

Предлагаемая статья направлена на разработку математической модели процесса взаимодействия движителя машины с грунтом, свойства которого изменяются непосредственно в процессе взаимодействия. Базой для исследования служат положения теории движения автомобильного транспорта в условиях бездорожья, математический анализ. Основой математической модели является схема взаимодействия штампа-движителя и деформируемого массива грунта ограниченной толщины, размещенного на жестком недеформируемом основании. При разработке модели использована степенная зависимость модуля деформации от плотности грунта, изменяющейся под воздействием штампа-движителя вслед за относительной деформацией сжатия массива грунта. Решение уравнения для осадки штампа получено в виде полиномов Тейлора, это позволяет получить расчетное значение глубины колеи не прибегая к итерациям, что упрощает разработку автоматизированных рабочих мест, предназначенных для оценки показателей взаимодействия техники с грунтами лесосек на стадии планирования технологического процесса заготовки древесины. На примере взаимодействия движителей лесных машин со слабыми лесными грунтами показано, расчетные значения глубины колеи, полученные с использованием предлагаемой математической модели, заметно отличаются от значений, полученных расчетом при постоянном значении модуля деформации. Так, при давлении в диапазоне 0,06-0,1 МПа, характерном для взаимодействия колесных движителей лесных машин с грунтом, прогнозируемые значения глубины колеи ниже на 30-60 %. В заключении отмечаются перспективные направления дальнейших разработок, заключающиеся в исследование изменения физико-механических свойств грунта, характеризующих сопротивление грунтов сдвигу и разработке математических моделей, учитывающих потерю несущей способности грунта, расчетное значение которой также изменяется при уплотнении.

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

MATHEMATICAL MODEL OF THE IMPACT OF THE PROPELLING DEVICE ON THICKENING FOREST

SOIL

Research article

Andronov A.V.1' *, Kotenev E.V.2, Pegov V.Yu.3

1 2 3 St. Petersburg State Forestry University named after S.M. Kirov, St. Petersburg, Russia

* Corresponding author (andronovalexandr[at]gmail.com)

Abstract

The goal of the presented article is to develop the mathematical model of the interaction process between a track mover of a machine and a ground, its properties change directly during the interaction. The provisions of the theory of road transport movement under off-road conditions are the basis for the study, and mathematical analysis. The basis of the mathematical model is the interaction scheme between the track mover and the deformable soil massif of limited thickness placed on a rigid non-deformable base. When developing the model, the authors used a power-law dependence of the deformation modulus on the density of the soil, changing under the influence of the mover stamp following the relative compression of the soil mass. The solution of the equation for stamp precipitation is obtained in the form of Taylor polynomials, this allows to obtain the calculated value of the track depth without resorting to iterations, which simplifies the development of automated workplaces designed to estimate the indicators of the interaction of equipment with the cutting areas at the stage of wood harvesting. On the example of the interaction of thrusters of forest machines with weak forest soils it is shown that the calculated values of the depth of the gauge, obtained with the use of proposed mathematical model, differ markedly from the values obtained by calculation at a constant value of the deformation modulus. So, having a pressure within the range of 0.06-0.1 MPa, which is characteristic of the interaction of wheeled movers of forest machines with the ground, the predicted values of the track depth are lower by 30-60%. In conclusion, promising areas for further development are noted, consisting in the study of changes in the physical and mechanical properties of the soil, which characterize the resistance of soils to shear and the development of mathematical models that take into account the loss of bearing capacity of the soil, the calculated value of which also changes during compaction.

Keywords: forest machines, soil, mathematical analysis, deformation, track depth, soil strengthening, timber harvesting.

Введение

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

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

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

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

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

рисунке 1.

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

Сжатие элементарного слоя массива грунта определяется по формуле Я.С. Агейкина [7], [8]:

dh =-йг

1 -е

(1)

где е - относительная деформация элементарного слоя грунта, г - координата, отсчитываемая от площадки контакта движителя с грунтом.

Суммарная деформация сжатия массива грунта находится по формуле [7], [8]:

Н-h

h =| -йг

о 1 -е

где Н - толщина деформируемого слоя грунта.

Традиционно решение задачи строится с использованием линейного выражения [7], [8]:

(2)

1

е = —а Е

(3)

где Е - модуль деформации грунта, а - нормальное напряжение в массиве грунта, изменяющееся по глубине [7],

[8]:

а =

Зр

1 +

аВ

(4)

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

2

г

Модуль деформации в формуле (3), как правило, считается постоянной величиной, не зависящей от ст, е [4], [5],

[6].

Необходимо учесть изменение модуля деформации Е при уплотнении. В геометрической интерпретации зависимости «напряжение-сжатие» модуль деформации Е представляет собой тангенс угла наклона касательной к графику напряжения ст от деформации е. Тогда при нелинейной функции ст(е) «мгновенный» модуль деформации Е определяется, как:

^ до

Е — — (5)

ае

Относительная деформация е связана с плотностью грунта р формулой [9]:

р-р°

(6)

1 -е

где р0 - начальная плотность грунта до воздействия напряжения ст.

Результаты исследований показывают, что модуль деформации Е выражается через плотность грунта степенной зависимостью [9], [10]:

г- Ь

Е — аЕ р

(7)

где аЕ, ЬЕ - числовые коэффициенты, зависящие от свойств грунта, его типа и состояния.

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

Таблица 1 - Свойства лесных грунтов различных категорий [9]

Категория лесного грунта р, г/см3 Е, МПа Н, м

Слабый (III) 0.75 0.4 0.8

Нормальный (II) 1 1 0.4

Прочный (I) 1.25 3 0.3

По данным таблицы 1 при помощи метода наименьших квадратов установим, что для лесного грунта коэффициенты в формуле (7): аЕ = 1,16, ЬЕ = 3,91.

Переписать уравнение (7) с учетом выражения (6) возможно следующим образом:

/■ \ЬЕ Г Л

Ро 1 ~ ( 1

Е — ас

1 -е

— Еп

1 -е

где Е0 - модуль деформации грунта по формуле (7) при р = р0, то есть до воздействия движителя. На основании формул (5), (8), получили:

де

— ^ (1 -е)Ь

ао Еь

Решение дифференциального уравнения (9) при начальном условии е(0) = 0 имеет вид:

(

е-1 -

Еп

Ль„ -1

■(Ье -1) + Е0

(8)

(9)

(10)

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

(

1 -е

1

(ье - 1)о + Е01 = /

Т П!

То

П—1

(11)

гдеУ - производная функции/порядка п по ст в точке 0, коэффициенты д„ при стп:

1

—

Е

Ье - 2 2Е2

Е

е

_ 2Ь2 - 7ЬЯ + 6

Я3 = зЕ

Я =

0

6ЬЕ - 29ЬЕ + 46ЬЕ - 24

4Е0

Коэффициенты д„ не зависят от г, тогда запишем уравнение (2) следующим образом:

к Н-И

И I (13)

п=1 о

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

, 1 х й^1 (2п-1)(2п - зУ2п - 5)...(2п - 2 ] +1) (2п - 3)| Ь-^йх =-^-£-А-¿V-+ V ) агс1ёх (14)

1(1 + х2)" 2п- 1Й 2^(п- 1)(п-2)...(п-])(1 + х22 (п-1) * ( )

Тогда определенные интегралы от степеней п напряжения а можно определить по формуле (13):

. г „ / Н - К

\ = ЗраВагод-

\ аВ

„ З2 р2а2В2 (Н - к) З2 р2аВ ( Н - h I

¿г = / , - У-^ + —-агсЫ-I (15)

2 2(а2В2 +(Н - к)2) 2 \ аВ )

33 ръаА Б4 (Н - И ) 373 ръа2 В2 (Н - И ) 3/Уа£ г (Н-И! Л = 4 (а2Б2 +(Н - И)2 )2 + 8 (а2Б2 +(Н - И)2) 8 ^ I аБ )

_ З4р4а6В6 (Н - к) 5З4р4а 4В4 (Н - к) 5З4р4а 2В2 (Н - к) ¿ЗраВ ( Н - к I 4 = б(а2В2+(Н - к)2 ) + 24(а2В2+(Н - К)2 )2 + 1б(а2В2+(Н - к)2 ) + 16 ^Ч аВ )

Уравнение для определения сжатия массива грунта (2) с учетом выражений (13), (15) можно представить следующим образом:

к

к (16)

п=1

где Ь - длина пятна контакта движителя с опорной поверхностью.

Расчеты при варьировании В = 0,4...0,8 м, Ь = 0,8...5 м, Е, Н по таблице 1, И = 0...0,8Н, показывают, что в уравнении (16) можно ограничиться к = 2 практически без потери точности, тогда:

И = ЯЛ +Я2Л =

(17)

ЗраВ ( Н - К -агсщ

Еп

Л аВ

+

З 2 р 2 аВГ

4Е2

аВ(Н - К) (Н -, 9 -г + агсЫ-

а2В2+(Н - К)2 \ аВ

Более того, в тех же диапазонах варьирования:

аВ(Н - К) аВ ( Н - К

, , -^-Г « — агсщ-

а2В2+(Н - К)2 Н \ аВ

агС^

ГН-И! 3 (Н-И

аБ ) 4 27

аБ

(18) (19)

С учетом (12), (15), (18), (19), получим из (17):

h = 3 JP(H -h)-1J2P2 (bE -2)(aB + H) (h _h)

4 E0V Мб E02H v '

(20)

Решив уравнение (20) относительно А, получим:

h = H

1 + pC£

(21)

где С - жесткость массива грунта с учетом геометрии пятна контакта, £ - безразмерная функция упрочнения грунта:

C = 3J

4 E0 4=1-XP

где х - коэффициент упрочнения грунта при уплотнении:

(22) (23)

, = 1J (¡Ё+П {ЬЕ - 2)

4 Eo H

(24)

Параметры a, J в выражениях (22), (24) рассчитывается по формулам [7], [8], [10]:

B + H H

0,03B + L 0,43B + 0,6L

a = 0,64:

J = ■

(26)

(25)

Пример результатов расчета и выводы. Сравнение результатов расчета сжатия h по формуле (21) с решением (2) при E = const представлено графиками на рисунке 2 (B = 0,7 м, L = 1,3 м, III категория лесного грунта).

0,2 0,15 0,1 0,05

0,02 0,04 0,06

p, МПа

0,08

0,1

—□—сжатие массива грунта с учетом упрочнения —О—сжатие массива грунта без учета упрочнения Рис. 2 - Деформация сжатия лесного грунта III категории в зависимости от давления движителя

Как показывают графики, расчетные значения глубины колеи И, полученные с использованием предлагаемой математической модели, заметно отличаются от значений, полученных расчетом при постоянном значении модуля деформации. Так, при давлении в диапазоне 0,06-0,1 МПа, характерном для взаимодействия колесных движителей лесных машин с грунтом, прогнозируемые нами значения И ниже на 30-60 %.

0

0

Для наглядности представим результаты расчета изменения Е, р под воздействием движителя - рисунок 3. 0,95

p, МПа

Рис. 3 - Изменение модуля деформации и плотности лесного грунта III категории в зависимости от давления

движителя

Расчеты показывают, что при давлении в диапазоне 0,06-0,1 МПа плотность грунта под воздействием движителя увеличивается ориентировочно на 10-15 %, как следствие - модуль деформации увеличивается на 50-90 %.

По нашему мнению, предлагаемая математическая модель позволяет при расчете учесть упрочнение грунта под воздействием движителя, развивая тем самым базу для теоретических исследований взаимодействия движителей лесных машин с опорными поверхностями. Уравнение (21), в отличие от известных ранее, позволяет получить расчетное значение глубины колеи h не прибегая к итерациям, что упрощает разработку автоматизированных рабочих мест (АРМ), предназначенных для оценки показателей взаимодействия техники с грунтами лесосек на стадии планирования технологического процесса заготовки древесины.

Дальнейшие направления исследований, на наш взгляд, следует сосредоточить в следующих направлениях:

1. Исследование изменения физико-механических свойств грунта, характеризующих сопротивление грунтов сдвигу.

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

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

Конфликт интересов Conflict of Interest

Не указан. None declared.

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

1. Герасимов Ю. Ю, Сюнёв В. С. Экологическая оптимизация технологических машин для лесозаготовок / Ю. Ю. Герасимов, В.С. Сюнёв. -. Йоэнсуу: университет Йоэнсуу, 1998. - 178 с.

2. Анисимов Г. М., Большаков Б.М. Основы минимизации уплотнения почвы трелевочными системами / Г. М. Анисимов, Б.М. Большаков. - СПб. : ЛТА, 2000. - 106 с.

3. Анисимов Г. М. О средствах научного проникновения в исследовании лесных машин и оборудования / Г.М. Анисимов, А. М. Кочнев // Известия ВУЗов. Лесной журнал. - 2010. - №4. - С. 39-46.

4. Божбов В. Е. Повышение эффективности процесса трелевки путем обоснования рейсовой нагрузки форвардеров : Автореферат дисс. ... канд. техн. наук. / Божбов Владимир Евгеньевич. - Архангельск: САФУ., 2015. -20 с.

5. Устинов В. В. Оценка тягово-сцепных свойств колесных движителей лесных машин методами теории движения автотранспорта по бездорожью : Автореферат дисс. . канд. техн. наук. / Устинов Владимир Владимирович. -Архангельск: САФУ., 2016. 20 с.

6. Хахина А. М. Методы прогнозирования и повышения проходимости колесных лесных машин : дис. ... докт. техн. наук : 05.21.01 / Хахина Анна Михайловна. - Архангельск: САФУ., 2018. - 318 с.

7. Агейкин Я. С. Проходимость автомобилей / Я. С. Агейкин - М.: Машиностроение, 1981. - 232 с.

8. Агейкин Я. С. Проходимость автомобиля / Я.С. Агейкин, Н.С. Вольская, И.В. Чичекин. - М.: МГИУ, 2010. -275 с.

9. Хитров Е. Г. Повышение эффективности трелевки обоснованием показателей работы лесных машин при оперативном контроле свойств почвогрунта : Автореферат дисс. ... канд. техн. наук : 05.21.01 / Хитров Егор Германович. - Архангельск: САФУ., 2015. - 20 с.

10. Ларин В. В. Методы прогнозирования и повышения опорной проходимости многоосных колесных машин на местности : дисс. ... докт. техн. наук : 05.05.03 / Ларин Василий Васильевич. - Москва: МГТУ им. Н.Э. Баумана., 2007.

29

- 530 с.

Список литературы на английском языке / References in English

1. Gerasimov Ju. Ju. Jekologicheskaja optimizacija tehnologicheskih mashin dlja lesozagotovok [Ecological optimization of technological machines for logging] / Ju. Ju. Gerasimov, V.S. Sjunjov. -. Jojensuu: universitet Jojensuu [University of Joensuu], 1998. - 178 p. [in Russian]

2. Anisimov G. M. Osnovy minimizacii uplotnenija pochvy trelevochnymi sistemami [Basics of minimizing the soil compaction under skidding systems] / G. M. Anisimov, B. M. Bol'shakov. - SPb. : LTA, 2000. - 106 p. [in Russian]

3. Anisimov G. M. O sredstvah nauchnogo proniknovenija v issledovanii lesnyh mashin i oborudovanija [On the means of scientific methods in the study of forest machines and equipment] / G. M. Anisimov, A. M. Kochnev // Izvestija VUZov. Lesnoj zhurnal [News of Universities. Forest Magazine]. - 2010. - №4. - P. 39-46. [in Russian]

4. Bozhbov V. E. Povyshenie jeffektivnosti processa trelevki putem obosnovanija rejsovoj nagruzki forvarderov [Improving the efficiency of the process of skidding by justifying the forwarders travel load] : Dissertation abstract ... of P hD in engineering. / Bozhbov Vladimir Evgen'evich. - Arhangel'sk: SAFU, 2015. - 20 p. [in Russian]

5. Ustinov V. V. Ocenka tjagovo-scepnyh svojstv kolesnyh dvizhitelej lesnyh mashin metodami teorii dvizhenija avtotransporta po bezdorozh'ju [Assessment of traction characteristics of wheeled rovers of forest machines using the theory of off-the-road vehicle locomotion] : Dissertation abstract ... of PhD in engineering / Ustinov Vladimir Vladimirovich. -Arhangel'sk: SAFU., 2016. 20 p. [in Russian]

6. Hahina A. M. Metody prognozirovanija i povyshenija prohodimosti kolesnyh lesnyh mashin [Methods of forecasting and improving mobility of wheeled forest machines] : Dissertation abstract ... of PhD in engineering : 05.21.01 / Hahina Anna Mihajlovna. - Arhangel'sk: SAFU., 2018. - 318 p. [in Russian]

7. Agejkin Ja. S. Prohodimost' avtomobilej [Passability of Automobiles] / Ja. S. Agejkin - M.: Mashinostroenie, 1981. -232 p. [in Russian]

8. Agejkin Ja. S. Prohodimost' avtomobilja [Passability of Automobile] / Ja. S. Agejkin, N. S. Vol'skaja, I. V. Chichekin. -M.: MGIU, 2010. - 275 p. [in Russian]

9. Hitrov E. G. Povyshenie jeffektivnosti trelevki obosnovaniem pokazatelej raboty lesnyh mashin pri operativnom kontrole svojstv pochvogrunta [Improving the efficiency of skidding by justifying the performance of forest machines in the operational control of the soil properties] : Dissertation abstract ... of PhD in engineering : 05.21.01 / Hitrov Egor Germanovich. - Arhangel'sk: SAFU., 2015. - 20 p. [in Russian]

10. Larin V. V. Metody prognozirovanija i povyshenija opornoj prohodimosti mnogoosnyh kolesnyh mashin na mestnosti [Methods of forecasting and increasing the mobility of multi-axle wheeled vehicles on the ground] : Dissertation abstract ... of PhD in engineering : 05.05.03 / Larin Vasilij Vasil'evich. - Moskva: MGTU im. N.Je. Baumana., 2007. - 530 p. [in Russian]