CC BY
40
УДК 630*332.2.001.57
UDC 630*332.2.001.57
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ УПРАВЛЕНИЯ ПРОЦЕССОМ ОБРЕЗКИ КРОН ДЕРЕВЬЕВ МАШИНОЙ
MATHEMATICAL MODEL OF MANAGEMENT OF SCRAP PROCESS OF CRONES OF TREES BY MANIPULATOR
МАНИПУЛЯТОРНОГО ТИПА С ДИСКОВОИ WITH A DISK SAW ПИЛОЙ
Попиков Пётр Иванович д.т.н., профессор
Popikov Peter Ivanovich Dr.Sci.Tech., professor
Драпалюк Михаил Валентинович д.т.н., профессор
Drapaluk Michael Valentinovich Dr.Sci.Tech., professor
Попиков Виктор Петрович к.т.н., старший преподаватель
Воронежская государственная лесотехническая академия, Воронеж, Россия
Popikov Viktor Petrovich Cand.Tech.Sci., senior teacher
Voronezh State Academy of Forestry and Technologies, Voronezh, Russia
В статье представлена математическая модель управления взаимодействия рабочего органа обрезчика с ветвями деревьев. Выполнена оптимизация конструктивных параметров машины для обрезки деревьев
In this article, the mathematical model of management of interaction of working body of the cutter with branches of trees is submitted. Optimization of design data of the machine for scrap of trees is executed
Ключевые слова: ОБРЕЗКА ВЕТВЕЙ, МАТЕМАТИЧЕСКАЯ МОДЕЛЬ
Keywords: SCRAP OF BRANCHES, MATHEMATICAL MODEL
В общем объеме лесохозяйственных работ важное место занимает переход к внедрению промышленных методов лесовыращивания на основе лесного селекционного или сортового семеноводства. Лесным кодексом РФ предусматривается создание постоянных лесосеменных участков и лесосеменных плантаций. Одним из основных направлений реализации федеральной целевой программы «Развитие лесного семеноводства на период 2009-2020 годы» является обеспечение выполнения всего комплекса мероприятий по содержанию имеющихся лесосеменных объектов и повышение их урожайности за счет лесоводственных и агротехнических уходов, включающих обрезку кроны, внесение удобрений и др.
Существующие конструкции машин для обрезки ветвей деревьев обеспечивают высоту обрезки не более 4,5 м, в то время как на лесосеменных плантациях и участках требуемая высота обрезки 7-8 м. В настоящее время обрезка ветвей деревьев лесосеменных плантаций производится рабочими, поднятыми гидроподъемниками в крону деревьев,
с помощью ручных ножовок и бензопил. Любые повреждения при обрезке вредны для дерева и способствуют возникновению различных инфекционных болезней. Поэтому при обрезке крон необходимо обеспечивать качественный срез ветвей без задиров коры и расщепов. Задача повышения качества и производительности обрезки крон деревьев вызывает острую необходимость в разработке средств механизации этого сложного технологического процесса.
Из анализа исследований параметров существующих технических устройств для обрезки крон деревьев следует, что наиболее перспективными являются машины манипуляторного типа, рабочим органом которых являются дисковые пилы с гидроприводом. Однако рабочие процессы и параметры технологического оборудования машин для обрезки крон деревьев лесосеменных плантаций недостаточно исследованы, нет четкого представления о том, какими должны быть параметры гидропривода и зубьев дисковых пил для качественной обрезки крон деревьев.
Предлагается новый рабочий орган в виде дисковой пилы с гидроприводом, смонтированный на манипуляторе посредством шарнирного соединения маятникового рычага с удлинителем рукояти, поворот которого осуществляется поворотным гидродвигателем (ротатором) на угол более 180°, что расширяет зону действия манипулятора и его технологические возможности (рисунок 1) (свид-во на пол. модель № 30057). На маятниковом рычаге предусмотрена установка V-образного кронштейна для устранения зажима пилы и удаления срезанных ветвей из кроны дерева. Кроме того, для повышения качества обрезки предложена односторонняя заточка длинной и короткой боковых режущих кромок зубьев дисковой конической пилы и односторонний развод зубьев в сторону отделяемой части ветви (рисунок 1) (патент изобретения № 2374824).
При обрезке ветвей деревьев резание производится длинной боковой режущей кромкой, которая имеет форму ножевого секатора, поэтому происходит как бы силовое резание, а на ветви остается участок с гладким срезом, так как задний угол близок к нулю, при этом деформируется только отделяемая часть ветвей. При ежегодной подрезке однолетних побегов семенных деревьев и кустарников используется эта же пила, но с обратным направлением вращения. Резание в этом случае производится короткой боковой режущей кромкой, которая не дает отклоняться побегам в процессе срезания в сторону. При математическом моделировании устройство для обрезки ветвей рассматривалось как сложная многозвенная система, обладающая большим количеством параметров и погруженная при этом в среду с целой совокупностью параметров, не являющихся частью разрабатываемой системы.
13 12 11 10 9 8 7 6 5
Рисунок 1 - Схема машины для обрезки крон деревьев: 1 - базовая машина;
2 - поворотная колонна; 3 - гидроцилиндр стрелы; 4 - стрела; 5 -
гидроцилиндр рукояти; 6 - рукоять; 7 - гидроцилиндр удлинителя; 8 -
удлинитель рукояти с винтовой парой; 9 - гидродвигатель (ротатор); 10 -маятниковый рычаг; 11 - гидромотор пилы; 12 - дисковая пила; 13 - V-образный упор.
При этом окружающая обрезчик среда и сам обрезчик, в зависимости от комбинации своих внутренних параметров, взаимодействуют друг с другом, оказывая влияние на качественные показатели процесса подрезки корней.
При моделировании все факторы, оказывающие влияние на технологический процесс обрезки ветвей, были разделены на первичные (к которым можно отнести ветви, вершины, сучья т.д.) и вторичные (к которым можно отнести технологические особенности обрезки ветвей, механическое воздействие рабочих органов на ветви, а также конструктивные особенности обрезчика).
Для решения задачи моделирования обрезчика как механической системы был проведен анализ рабочего процесса обрезки ветвей с учётом факторов, оказывающих влияние на технологический процесс обрезки ветвей. Результаты анализа представлены в виде блок-схемы на рисунке 2.
Рисунок 2. Блок-схема динамической системы рабочего процесса обрезчика Обрезчик является сложной динамической системой, состоящей из взаимодействующих между собой подсистем: базовой машины, системы управления, рабочего оборудования, деревьев и перерезаемых ветвей.
Базовая машина представлена подсистемами: двигатель, трансмиссия и ходовое оборудование. Система управления представлена подсистемами:
человек-оператор и устройство управления. Рабочее оборудование представлено подсистемами: гидропривод и рабочий орган. Ветви
представлены подсистемами: реакция ветвей на рабочий орган и физикомеханические свойства ветвей.
При работе устройства для обрезки ветвей человек-оператор (или автоматическое устройство), исходя из условий эффективности рабочего процесса (А1, А2), оказывает влияние (А3) на устройство управления и формирует управляющие воздействия на трансмиссию (А5) базовой машины и на гидропривод (А4) рабочего оборудования.
Подсистема базовой машины, воспринимая воздействия от устройства управления, перераспределяет энергию между двигателем и трансмиссией (А6, А7), оказывая при этом влияние на ходовое оборудование (А8) и человека-оператора (А9). Ходовое оборудование, в свою очередь, воздействует (А11) на устройство управления, на которое одновременно может оказывать влияние и двигатель (А10) базовой машины.
Подсистема рабочего оборудования, воспринимая воздействия от устройства управления, передаёт энергию от гидропривода рабочему органу и в обратном направлении (А13, А14). На гидропривод может оказывать влияние и трансмиссия (А12) базовой машины.
Подсистема обрезаемых ветвей воспринимают воздействия от рабочего органа (А15), оказывая на него ответное влияние (А16). При этом реакции ветвей зависит от их физико-механических свойств (А17).
На основе анализа рабочего процесса подрезки ветвей была составлена математическая модель данного процесса. При моделировании взаимодействия пилы с ветвью согласно методу конечных элементов участок ветви представляется в виде совокупности большого количества элементарных кубов малого размера ё. Ветвь в модели первоначально
представляет собой геометрическую область, имеющую форму цилиндра радиусом Яв и высотой Ьв (рисунок 3).
Каждый зуб пилы в модели представляет собой пятигранник (рисунок 4). Вершинами зуба-пятигранника являются шесть точек, три из которых (точки 1, 2, 3) лежат в нижней плоскости пильного диска, а оставшиеся три
точки 4, 5, 6) - в верхней плоскости пильного диска.
Координаты точек (х1, у1, 21) пятигранника по отношению к системе координат, связанной с пильным диском, выражаются следующим образом:
Рисунок 3 - Область пространства размером Ьх □ □ Ьг, в которой
производится моделирование и первоначальное представление ветви в форме цилиндра.
Рисунок 4 - Пятигранная форма зуба пилы в модели (вид сверху на плоскость диска, зуб не разведен)
х1 = Г2-соб ф2; у1 = Г2-біп ф2; 2І = 0; х2 = Г2-соБ(ф2 + 2л/№уб); у2 = Г2-БІп(ф2 + 2л/№уб); 22 = 0; х3 = (г2 + Дг)-соб(ф2 + Дф); у3 = (г2 + Дг)-БІп(ф2 + Дф); 23 = 0; х4 = хМС; у4 = уМС; 24 = Ь2; х5 = Г2-соб ф2; у5 = Г2-БІП ф2; 25 = Ь2; х6 = Г2-соБ(ф2 + 2л/№уб); у6 = Г2-БІп(ф2 + 2л/№уб);
26 = Ь2, (1)
где, хх - радиус пильного диска; №уб - количество зубьев на диске; Дг = 1,06858-2л/№уб-гх - высота зуба; Дф = 0,38893-2л/№уб - угловое расстояние между основанием и краем зуба; хМС и уМС - координаты точки 4, определяемые методом Монте-Карло в программе, реализующей модель; Ьх - толщина диска.
Составлены и проанализированы дифференциальные уравнения, описывающие расход рабочей жидкости гидропривода и рабочий процесс резания дисковой пилой при повороте маятникового рычага поворотным гидродвигателем:
= -ГГ- (Ч н • Щн - Ч мЩ - ауР - Ч г • 9 п )' к р
пр
за р ёК
п ! т -Г
п "1 т Г _________ К
км
к И +
п под
с И
рез под
Яр (щщ- к щЩ
(2)
где р - давление рабочей жидкости в гидросистеме, Па; 1 - время, с; Кр -коэффициент податливости упругих элементов гидропривода, м5/Н; qн -рабочий объём насоса, м3/об; qм - рабочий объём гидромотора, м3/об; □ н - угловая скорость вращения насоса, с-1; □ - угловая скорость вращения вала гидромотора, с-1; ау - коэффициент утечек,м5/(Нс); qг - рабочий объём ротатора, м3/об; фп - угол поворота маятника, рад.; Дпр -приведенный момент инерции вращающихся масс к валу гидромотора, включающий момент инерции роторной группы гидромотора и момент инерции рабочего органа, кг-м2; пп - полный объёмный КПД гидромотора; ^0 - объёмный КПД гидромотора; кМ - коэффициент, определяющий силу сопротивления при удалении элементарного куба, с; Кр - количество
V
И = 0,12^-под под 2 па
удаленных элементарных кубов ветви; 2 па - подача на режущий
зуб, м; упод - скорость подачи,м/с; кп - удельная сила резания, Н/м2; ^ -коэффициент трения древесины о зуб; 5 - угол резания передней режущей
1
Р
кромки; срез - коэффициент пропорциональности, постоянный для данной обрабатываемой древесины; Яр - радиус резания; кю - коэффициент вязкого сопротивления резанию.
Система дифференциальных уравнений (2) решена методом численного интегрирования - модифицированным методом Эйлера-Коши. Для решения системы дифференциальных уравнений, положенной в основу модели, и для проведения различных компьютерных экспериментов с моделью составлена компьютерная программа на языке Object Pascal в интегрированной среде программирования Borland Delphi
7.0. На экран выводятся графики угловой скорости и давления на гидромоторе пилы.
ю, рад/с 1
t, с
Рисунок 5 - Характер зависимостей: 1 - угловой скорости вращения диска ю от времени 1, 2 - давления на гидромоторе Р от времени 1
Зависимость силы резания Бшах от скорости подачи (упод) имеет приближенно линейный характер на участках от 0,01 до 0,05 м/с (рисунок
6, а). Возрастающий характер зависимости можно объяснить тем, что с увеличением скорости подачи более 0,06 м/с возрастает объем древесины, срезаемый пилой в единицу времени ёУд/Л, что может привести к зажиму пилы в пропиле. Увеличение работы А срезания ветви с увеличением упод (рисунок 6, б) можно объяснить тем, что с увеличением скорости подачи возрастают силы сопротивления резанию.
А, кДж 15 -
10-
10
20
30
р, градусы
40
р, градусы
Rm,max? мкм
500 400300 200100
0,00
симметричныи развод верхний развод
10
20
30
40
р, градусы
Рисунок 6 - Зависимости максимальной силы резания на зубе Fmax (а), работы A срезания ветви (б) и шероховатости среза Rm,max (в) от угла развода зубьев пилы р
Проведена серия компьютерных экспериментов с различными схемами развода зубьев пилы: симметричный развод - зубья разведены через один вверх и вниз на одинаковый угол р от плоскости диска; верхний развод - зубья с нечетными номерами не отклонены от плоскости диска, а зубья с четными - отклонены вверх на угол р. Во всех сериях угол р изменяли от 5° до 40° с шагом 5°. Возрастающий характер зависимостей Fmax^) и А(р) для двух схем развода можно объяснить тем, что с увеличением р возрастает как ширина, пропила, так и интенсивность воздействия зубьев на древесину в области резания (рисунок 6 а, б, в). При одностороннем разводе наблюдается снижение энергоемкости пиления по сравнению с симметричным разводом от 5 до 10 кДж (рисунок 6 б). При верхнем разводе шероховатость Rm, max среза ветви минимальна и составляет около 200 мкм (рисунок 6, в). В то же время при симметричном разводе шероховатость принимает большие значения (от 260 мкм) и возрастает с увеличением угла развода до 500 мкм.
5
0
б
a
в
В процессе обрезки кроны дерева часто приходится удалять ветви, расположенные не перпендикулярно к плоскости пилы. Угол встречи 0 между плоскостью пильного диска и стволом ветви может составлять от 60° до 120°. Установлено, что Бшах имеет параболическую зависимость от угла 0 (рисунок 7, а). Работа А увеличивается с увеличением площади среза, поэтому имеет зависимость, близкую к закону А=А0/соб 0 (рисунок
7, б).
Проведена серия экспериментов, при которой угол резания 5 изменялся от 60° до 120° с шагом 10°.
^шах, Н А, кДж
Рисунок 7 - Зависимости максимальной силы резания на зубе Бшах (а) от скорости упод подачи пильного диска (1) и угла встречи Д9 (2) и работы А срезания ветви (б) от скорости упод подачи пильного диска (1) и угла встречи 9 (2)
Установлено, что Бшах практически не зависит от угла 5 (рисунок 8, а), в то время как работа срезания ветви А уменьшается с увеличением 5 (рисунок 8, б). В интервале от 120 до 160° работа резания практически не изменяется и лежит в пределах 7,5...8,0 кДж. Поэтому предлагается угол резания для длинной боковой режущей кромки принимать в пределах
140...160°. Максимальная сила резания Бшах незначительно изменяется от угла заточки в зуба (рисунок 8, а). В то же время работа резания ветви увеличивается с 7,5 кДж до 9,5 кДж при увеличении угла заточки в от 40°...80° (рисунок 8, б). Поэтому угол заточки зуба принимается равным 30°...40° с минимальными энергозатратами.
Порода деревьев оказывает существенное влияние на возникающие на зубьях пилы силы и энергетические затраты по обрезке ветви. Проведена серия экспериментов, при которых менялась твердость древесины от 4,6 до 50,6 Н/мм2 с шагом 4,6 Н/мм2.
Максимальная сила резания меняется от 57,7 Н для легкообрабатываемой древесины (сосна) до 173,7 Н для труднообрабатываемой древесины (береза). Работа резания ветви изменяется от 6,28 до 14,77 кДж соответственно.
40
60
80
ß, град.
40
60
80
б
Рисунок 8 - Зависимости максимальной силы резания Бшах (а) от угла резания 5 (1) и угла заточки в (2) и работы А срезания ветви (б) от угла резания 5 (1) и угла заточки в (2)
На величину скорости подачи пилы, которая определяет качество среза оказывают влияние следующие величины: подача на зуб - №, число зубьев - 7, диаметр пилы - Б и окружная скорость вращения пилы - Ур. Целевая функция при ограничениях варьируемых факторов имеет вид:
V(U2,Vp,Z,D) ® min,
при ограничениях : 0.001 < Uz < 0.006;20 < Vp < 40;90 < Z < 120; 0.4 < D < 0.65. (3)
В результате проведенной оптимизации выявлено, что оптимальными параметрами резания дисковой пилой являются окружная скорость 30 м/с, подача на зуб 0,012 мм, диаметр пилы 0,5 м, количество зубьев на дисковой пиле 100 шт.
Таким образом, при решении математической модели процесса резания с использованием метода конечных элементов, при решении
ß
а
которой установлено, что при одностороннем разводе зубьев пилы шероховатость среза ветвей составляет около 200 мкм, в то же время при симметричном разводе шероховатость среза находится в пределах
260...500 мкм. Энергоемкость процесса резания при одностороннем разводе зубьев пилы по сравнению с симметричным снижается в 1,4...1,5 раза. Рекомендуются следующие параметры дисковой конической пилы: диаметр пилы 0,5 м, количество зубьев на дисковой пиле 100 шт, подача на зуб 0,012.0,024 мм, скорость резания 26.30 м/с, развод зубьев 0,2.0,4 мм, угол резания для длинной боковой режущей кромки 160°, угол заточки зуба 30°...40°, для короткой боковой режущей кромки угол резания 50°.
Литература
1. Попиков, В. П. Механизация обрезки крон деревьев и кустарников [Текст] / Д. Д. Репринцев, М. В. Драпалюк, В. П. Попиков. // Лесн. хоз-во. - 2006. - № 1. -С. 45.
2. Попиков, В. П. Имитационное моделирование технологического процесса лесной машины с гидроприводом дискового рабочего органа [Текст] / В. П. Попиков, В. Н. Коротких, М. В. Драпалюк // Вестн. КрасГАУ. - 2009. - № 5 -С. 129-132.
3. Попиков, В. П. Моделирование процесса обрезки ветвей деревьев дисковой пилой на лесосеменных плантациях [Текст] / В. П. Попиков, Л. Д. Бухтояров // Вестн. КрасГАУ. - 2009. - № 8 - С. 3-7.
4. Свид-во 30057 РФ, МПК7 А 01 О 3/04. Рабочий орган машины для подрезки крон деревьев [Текст] / В. П. Попиков, Д. Д. Репринцев, П. И. Попиков ; заявитель и патентообладатель ВГЛТА. - № 2002129704/20 ; заявл. 05.11.02 ; опубл. 20.06.03, Бюл. № 17. - 2 с.
5. Пат. 2374824 РФ, МПК7 А010 3/00. Рабочий орган машины для подрезки крон деревьев [Текст] / В. П. Попиков, Л. Д. Бухтояров, М. Д. Драпалюк ; заявитель и патентообладатель ВГЛТА. - № 2008106730/12 ; заявл. 21.02.2008 ; опубл. 10.12.09, Бюл. № 34. - 7 с.
