CC BY
50
УЧЕНЫЕ ЗАПИСКИ ПЕТРОЗАВОДСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Март, № 2 Технические науки 2014
УДК 630*372/375
АНТОНИНА ИВАНОВНА НИКИФОРОВА
кандидат технических наук, доцент кафедры технологии лесозаготовительных производств, Санкт-Петербургский государственный лесотехнический университет им. С. М. Кирова (Санкт-Петербург, Российская Федерация) tlzp@inbox.ru
ЕГОР ГЕРМАНОВИЧ ХИТРОВ
аспирант кафедры технологии лесозаготовительных производств, Санкт-Петербургский государственный лесотехнический университет им. С. М. Кирова (Санкт-Петербург, Российская Федерация) yegorkhitrov@gmail. com
АЛЕКСАНДР АНДРЕЕВИЧ ПЕЛЫМСКИЙ аспирант кафедры технологии лесозаготовительных производств, Санкт-Петербургский государственный лесотехнический университет им. С. М. Кирова (Санкт-Петербург, Российская Федерация) tlzp@inbox.ru
ОЛЬГА ИВАНОВНА ГРИГОРЬЕВА
доцент кафедры лесоводства, Санкт-Петербургский государственный лесотехнический университет им. С. М. Кирова (Санкт-Петербург, Российская Федерация) tlzp@inbox.ru
ОПРЕДЕЛЕНИЕ ОСАДКИ ПРИ ДВИЖЕНИИ ЛЕСОЗАГОТОВИТЕЛЬНОЙ МАШИНЫ ПО ДВУСЛОЙНОМУ ОСНОВАНИЮ
Известен способ укрепления трелевочных волоков порубочными остатками с целью снижения экологического ущерба почвогрунту лесосеки, вызванного движением лесозаготовительной техники. Однако известные модели взаимодействия машины и лесного почвогрунта не достаточно полно раскрывают влияние различия физико-механических свойств слоя порубочных остатков и слоя грунта лесосеки, а также их изменение, вызванное уплотнением после прохода машины. В статье рассмотрено взаимодействие лесозаготовительной машины с основанием, состоящим из двух слоев: слоя измельченных порубочных остатков (верхний слой) и слоя грунта лесосеки (нижний слой). Предложена модель колееобразования, в явном виде учитывающая влияние различия физико-механических характеристик слоев основания, а также изменение свойств грунта после прохода машины. Рассмотрен алгоритм расчета глубины колеи при многократном прохождении техники по участку трассы и пример расчета глубины колеи после 25 проходов машины по трассе при варьировании удельного давления со стороны машины на основание и толщины слоя уложенных на него порубочных остатков.
Ключевые слова: грунт, колееобразование, циклическая нагрузка, многослойное основание, порубочные остатки
ВВЕДЕНИЕ
В последние годы в России все большее внимание уделяется экологическому аспекту проведения лесозаготовок [2]. Негативное влияние лесозаготовительной техники на лесной почвогрунт широко известно: переуплотнение, срез слоев почвогрунта лесосеки вызывают заметное снижение его качеств с точки зрения последующего лесовосстановления. Сохранение плодородия лесных почв является одним из магистральных направлений научно-технического прогресса в лесозаготовительном производстве [2].
Одним из путей решения задачи снижения экологического ущерба лесу, вызванного движением лесозаготовительной техники при осуществлении заготовки древесины, является укрепление трелевочных волоков лесосечными
отходами. Этот способ представляется более целесообразным в свете необходимости утилизации отходов лесозаготовительного производства, что также является одной из важнейших практических задач отрасли [3].
Научно обоснованное планирование организационно-технических мероприятий по укреплению трелевочных волоков лесосек невозможно без четкого описания процесса взаимодействия лесозаготовительных машин с лесным почвогрунтом, при этом одной из главных характеристик процесса является глубина образующейся при проходе машины колеи [2].
Как показал анализ литературных источников обозначенной тематики, вопрос колееобразования при многократном проходе машины по участку трассы изучен недостаточно [3], [4].
© Никифорова А. И., Хитров Е. Г., Пелымский А. А., Григорьева О. И., 2014
88
А. И. Никифорова, Е. Г. Хитров, А. А. Пелымский, О. И. Григорьева
В частности, отсутствуют математические модели, в явном виде учитывающие неоднородность основания, по которому движется машина, что имеет место при работе на укрепленном порубочными остатками участке волока.
В связи с изложенным дальнейшие исследования процесса взаимодействия лесозаготовительной техники с многослойным основанием представляются актуальными. В статье рассматривается процесс колееобразования при циклическом воздействии лесозаготовительной техники на грунт лесосеки, укрепленный измельченными порубочными остатками.
ПОСТРОЕНИЕ МОДЕЛИ И ЕЕ РЕАЛИЗАЦИЯ
Для определения глубины колеи, образующейся при проходе лесозаготовительной техники по волоку, укрепленному порубочными остатками, рассмотрим уплотнение двухслойного основания под действием сжимающей нагрузки.
Для построения модели вначале зададимся законом распределения вертикальных напряжений сжатия о по глубине в следующем виде [1]:
а(z ) = p — 1 +
(1)
где параметры J, a, D' определяются известными
выражениями:
= +---,
0,03 +1 b
0,6 + 0,43 1 b
(2)
где F - площадь штампа, b - ширина штампа, l - длина штампа, H - толщина деформируемого слоя, p - давление со стороны машины на грунт.
Известно [1], что в общем виде осадка слоя грунта высотой z2 - z1 (z - вертикальная координата сечения) от действия нормального давления о (z) может быть найдена из следующего уравнения:
Z1
Тогда, после подстановки (1) в (3) и интегрирования, для осадки можем записать:
h = h
Ps
P-Ps
(5)
При этом несущая способность основания определяется следующим образом:
1 2
Ps =-Ps 0П arctan
П(Z2 - Z1 )
2b
(6)
где p 0 находится из формулы [4]:
Ps
J 1Xу + J2 X2 + Xз hr.
(7)
Значения параметров J1, J2 для выражения (7) определяются так:
J =
l
l + 0,4b
, J 2
l + b l + 0,5b
(8)
Для расчета параметров X, X2, X3 используют следующие формулы [1]:
X Л(-Ш4) X _2C»+Ш2)
XI ^.5 , X 2
Xз =
Ш5
Y
Ш
Ш2
Ш = tan
n
(9)
где у - объемный вес грунта, C0 - внутреннее сцепление грунта, ф0 - угол внутреннего трения грунта.
Поскольку рассматривается двуслойное основание (схема на рис. 1), в выражении (2) примем
н=H + н2.
a = 1 +
b
H + H 2
(10)
где H Н2 - толщина соответственно первого (верхнего) и второго слоев основания.
Рис. 1. Схема к определению осадки двуслойного основания: 1 - слой измельченных порубочных остатков,
2 - слой грунта лесосеки
кл = p
гл *
JaD’
aD
arctan \ I - arctan I Zl
aD'
(4)
Необходимо учесть, что по мере приближения действующих сжимающих напряжений к несущей способности грунта ps вертикальное перемещение штампа-деформатора начинает интенсивно возрастать (потеря несущей способности) [1], поэтому для определения величины истинной деформации hr используем формулу:
В выражении (4) для пределов интегрирования нужно положить для первого слоя z1 = 0, z2 = H1 - hra1; для второго слоя z1 = H1 - hra1,
z2 = H2 - h 2.
2 2 гл2
Тогда для линейной деформации первого слоя основания h можем записать:
JaD’
= p-----arctan
E
H
aD'
(11)
Определение осадки при движении лесозаготовительной машины по двуслойному основанию
89
Для линейной деформации второго слоя hra2 также получим:
Ко. = Р
JaD’
~ёА
arctan
H - hrj]2
aD'
- arctan
h i - aml
aD1
.(12)
Выражения для определения несущей способности слоев примут вид: для первого слоя:
1 2 . П(Я1 - hrl) (13)
Psl = ^ Ps0in arCtan 2 b
для второго слоя:
При многократном проходе лесозаготовительной техники по одному и тому же участку волока необходимо учесть изменение физикомеханических свойств основания в результате предыдущего прохода машины [3].
В работе [5] предлагается определять модуль деформации E, внутреннее сцепление С0 и угол внутреннего трения ф0 грунта в зависимости от его коэффициента пористости е.
Для учета изменения пористости грунта после уплотнения примем зависимость физико-механических характеристик грунта в следующем виде:
ps2 = 4ps02п arctan
-H2 + hri-hгг)
2b
(14)
Учтем также, что величина ps0 в уравнениях (13), (14) для слоев грунта определяется так: для первого слоя:
Ps01 = J1 X1(1) + J2X2(1) + ^(Г^гЬ (15)
для второго слоя:
Ps 02 J1 X1(2) + J2 X2(2) + X3(2)hr2. (16)
Параметры Xl(l), X3(1), X1(2)P X2(2), X3(2) опре-
деляются по формулам (9) после подстановки значений физико-механических характеристик слоев основания (нижний индекс в скобках соответствует номеру слоя).
Выразим h через hг из (5) отдельно для каждого слоя:
= к
Р - Psl Psl
, к,
= к
Р - Ps2
Ps 2
(17)
В итоге получим следующие уравнения:
, p -psl JaD'
hrl-----— = p-------arctan
Psl El
H - hri
P - Psl
Psl
aD'
, (18)
( ~ \ Н h В Bs2
Р-Р. 2 „JaD' arctan 2 Пт2 Psl
Р 17 Psl Е1 aD'
\ )
-arctan
_________Psl
aD'
h = h, + k,
г г1 г2-
(19)
(20)
Нетрудно увидеть, что трансцендентные уравнения (18), (19) содержат искомые величины h h как в левой, так и в правой своей части, получить их решение в замкнутом виде не представляется возможным.
E(е) = (Л„ * BJl )[(1 -s)(е,„ *1) -1](Л“*‘“‘11, (21)
C(е) = (Л,, * B.:h)[(1 s)(е„ +1)- 1fC*“1, (22)
Ф, (е) = (Л., * В,L)[(1 - г) (е,„ *1) - 1](Л**V■ >, (23)
где £ - относительная деформация слоя грунта (2-й слой основания); е - коэффициент пористости, соответствующий максимально рыхлому состоянию грунта; А^ AbE, BbE, AaC, BaC,
Ас Bbc, A^ Baф, AbФ, ВьФ - некоторые постоянные для определенного типа грунта коэффициенты,
определяемые экспериментально, IL - показатель, характеризующий консистенцию грунта, определяемый, в свою очередь, так [5]:
I
L
W - Wp
Wl - Wp ’
(24)
где W - абсолютная влажность грунта, W -влажность грунта на границе пластичности, WL - влажность грунта на границе текучести.
После определения всех необходимых зависимостей можем перейти к расчетам глубины колеи, образующейся при многократном проходе лесозаготовительной техники по волоку, укрепленному слоем измельченных порубочных остатков.
Алгоритм расчетов следующий:
• задаемся типом грунта, его механическими свойствами (плотностью, пористостью, влажностью, влажностью на границе пластичности и раскатывания, толщиной деформируемого слоя) по [5]. По уравнениям (21)-(23) определяем значение модуля деформации, угла внутреннего трения и величину внутреннего сцепления;
• задаемся свойствами слоя порубочных остатков, укладываемых на волок (толщиной слоя, углом внутреннего трения, сцеплением и модулем деформации);
• принимаем удельное давление на грунт, параметры деформатора - ширину и длину, площадь контакта;
• по формулам (8), (10), (2), (13)-(16), (9) рассчитываем величины, входящие в уравнения (18)-(19);
90
А. И. Никифорова, Е. Г. Хитров, А. А. Пелымский, О. И. Григорьева
• далее производится численное решение уравнений (18)—(20) и c использованием найденных значений h ,, h h вычисляется относительная деформация второго слоя основания (грунта) после первого прохода машины;
• с использованием выражений (21)-(23) определяются значения физико-механических параметров слоя грунта после первого прохода машины, при этом относительная деформация определяется с использованием рассчитанной величины h ; полученные таким образом значения модуля деформации, угла внутреннего трения и внутреннего сцепления грунта используются при определении деформации слоя грунта после следующего прохода машины;
• свойства слоя порубочных остатков считаем постоянными, изменяется только его толщина (уменьшается после каждого прохода на величину hrf).
Для всех последующих проходов ход расчетов аналогичен описанному.
Как было отмечено ранее, уравнения (18)-(19) не имеют замкнутого решения, для их решения используются численные методы. Объем вычислений по описанному выше алгоритму значителен, на практике осуществление таких расчетов связано с существенными затратами времени, поэтому для проведения расчетов была составлена специальная программа в пакете смешанных вычислений Maple.
При проведении расчетов использовались следующие исходные данные:
F = 0,04 м2, b = 0,08 м, l = 0,5 м, Н1 = 0,2 ^ 0,4 м, Н2 = 0,75 м, р2 (0) = 1200 кг/м3, WL = 50 %, W = 20 %, W = 40 %, E1 = 6000 кПа, ф10 = 28°, С10 = 0, p = 100 ^ 250 кПа, грунт - суглинок.
Определим, как изменяется осадка первого и второго слоев основания в зависимости от давления на грунт. Решая численно уравнения
(18)-(20) при различных значениях давления р, получим графики, представленные на рис. 2.
Очевидно нелинейное влияние давления p на осадку первого и второго слоев основания. Осадка, начиная с давления р = 100 кПа, интенсивно увеличивается, что качественно согласу-
Рис. 2. Величина осадки слоев основания после однократного прохода машины (Н = 0,2 м)
ется с физической картиной деформации грунта
[1]: по мере приближения давления к несущей способности основания заметно влияние дополнительных касательных напряжений, происходит срез слоев грунта и осадка штампа растет непропорционально давлению.
Расчеты произведены при числе проходов N, равном 25, по трассе, основание которой представлено двумя слоями, в диапазоне давлений на грунт от 100 до 250 кН, результаты расчетов (глубина колеи, равная истинной деформации второго слоя основания h) представлены на графиках (рис. 3-5).
Рассмотренный в статье подход, на наш взгляд, позволит достаточно просто определять
Рис. 3. Глубина колеи в зависимости от давления на грунт, см (толщина слоя порубочных остатков 0,2 м): 1 - 0-5 см, 2 - 5-10 см, 3 - 10-15 см, 4 - 15-20 см, 5 - 20-25 см
Рис. 4. Глубина колеи в зависимости от давления на грунт (толщина слоя порубочных остатков 0,3 м): 1 - 0-5 см,
2 - 5-10 см, 3 - 10-15 см, 4 - 15-20 см
Рис. 5. Глубина колеи в зависимости от давления на грунт (толщина слоя порубочных остатков 0,4 м): 1 - 0-5 см,
2 - 5-10 см, 3 - 10-15 см, 4 - 15-20 см
Определение осадки при движении лесозаготовительной машины по двуслойному основанию
91
глубину образующейся колеи при многократном прохождении техникой участка волока. При этом будет в явном виде учтено:
• влияние на колееобразование толщины слоя уложенных на волок измельченных порубочных остатков;
• влияние на глубину колеи изменения физикомеханических свойств грунта после прохода машины.
Это, в свою очередь, позволит рассчитать достаточное для надежного укрепления волока количество укладываемых на него порубочных остатков.
СПИСОК ЛИТЕРАТУРЫ
1. Агейкин Я. С. Вездеходные колесные и комбинированные движители: теория и расчет. М.: Машиностроение, 1972. 184 с.
2. Григорьев И. В. Снижение отрицательного воздействия на почву колесных трелевочных тракторов обоснованием режимов их движения и технологического оборудования. СПб.: ЛТА, 2006. 236 с.
3. Григорьев И. В., Жукова А. И., Григорьева О. И., Иванов А. В. Средощадящие технологии разработки лесосек в условиях Северо-Западного региона Российской Федерации. СПб.: ЛТА, 2008. 176 с.
4. Григорьев И. В., Шапиро В. Я., Ру до в С. Е., Жукова А. И. Модель процесса циклического уплотнения грунта в полосах, прилегающих к трелевочному волоку // Вестник КрасГАУ № 2. Красноярск, 2010. С. 8-14.
5. Ларин В. В. Методы прогнозирования и повышения опорной проходимости многоосных колесных машин на местности: Дисс. ... д-ра техн. наук. М.: МГТУ им. Баумана, 2007. 419 с.
Nikiforova A. I., Saint-Petersburg State Forest Technical University under name of S. M. Kirov
(Saint-Petersburg, Russian Federation) Khitrov E. G., Saint-Petersburg State Forest Technical University under name of S. M. Kirov
(Saint-Petersburg, Russian Federation) Pelymskiy A. A., Saint-Petersburg State Forest Technical University under name of S. M. Kirov
(Saint-Petersburg, Russian Federation) Grigor’eva O. I., Saint-Petersburg State Forest Technical University under name of S. M. Kirov
(Saint-Petersburg, Russian Federation)
TWO-LAYER GROUND DRAUGHT DEFINITION UNDER LOGGING MACHINE PASSAGE
A method to strengthen skid trails with logging waste to reduce environmental damage to soil and ground caused by the logging machine passage is well known. However, exciting models of machines and forest soils’ interaction do not adequately reveal the impact of differences in physical and mechanical properties on the logging waste layer and cutting area of the soil layer. The depth of its change caused by compression after machine’s passage is not disclosed eihter. The paper considers interaction between two layers: a layer of shredded logging waste (top layer) and the soil layer of cutting area (bottom layer). Further, we describe a model of rutting that explicitly takes into account the effect caused by differences in physical and mechanical characteristics of the base layers and the change in soil characteristics after machines’ passage. The algorithm calculating the gauge depth after machines’ multiple passages is discribed. The article ends with an example of the rut depth calculation model for 25 machine passages. The study involves cases with variable volumes of pressure exercised by logging machines on the base and thickness of the logging waste layers.
Key words: soil, rutting, cyclic loading, multi-layer ground, logging waste
REFERENCES
1. Ageykin Ya. S. Vezdekhodnye kolesnye i kombinirovannye dvizhiteli: teoriya i raschet [Off-the-road wheeled and combined traction devices: Theory and Calculation]. Moscow, Mashinostroenie Publ., 1972. 184 p.
2. Grigor’ev I. V. Snizhenie otritsatel’nogo vozdeystviya na pochvu kolesnykh trelevochnykh traktorov obosnovaniem rezhimov ikh dvizheniya i tekhnologicheskogo oborudovaniya [Mitigation of negative impacts on the soil caused by wheeled skidders by selecting modes of motion and process equipment]. SPb., LTA Publ., 2006. 236 p.
3. Grigor’ev I. V., Zhukova A. I., Grigor’eva O. I., Ivanov A. V. Sredoshchadyashchie tekhnologii razrabotki lesosek v usloviyakh Severo-Zapadnogo regiona Rossiyskoy Federatsii [Eco efficient technologies of wood harvesting in the North-West region of the Russian Federation]. SPb., LTA Publ., 2008. 176 p.
4. Grigor’ev I. V., Shapiro V. Ya., Rudov S. E., Zhukova A. I. The process of cyclic compaction in the bands adjacent to skid trails model [Model’ protsessa tsiklicheskogo uplotneniya grunta v polosakh, prilegayushchikh k trelevochnomu voloku]. VestnikKrasGAU. № 2. Krasnoyarsk, 2010. P. 8-14.
5. Larin V. V. Metodyprognozirovaniya ipovysheniya opornoy prokhodimosti mnogoosnykh kolesnykh mashin na mestnosti: Diss. ... d-ra tekhn. nauk [Methods for predicting and improving multi-wheeled vehicles locomotion on the ground, Diss. doct. tehn. science]. Moscow, MGTU im. Baumana Publ., 2007. 419 p.
Поступила в редакцию 25.06.2013
