close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
Научный журнал КубГАУ, №79(05), 2012 года
1
УДК 629.114.2
UDC 629.114.2
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ
ДЕМПФИРОВАНИЯ РАБОЧЕГО
ПРОЦЕССА ГИДРОПРИВОДА МЕХАНИЗМА
ПОВОРОТА КОЛОННЫ ЛЕСНОГО
МАНИПУЛЯТОРА
OPERATION DAMPING MATHEMATICAL
MODEL OF THE FOREST MANIPULATOR
COLUMNS SLEWING MECHANISM OF
HYDRAULIC ACTUATOR
Долженко Сергей Валерьевич
аспирант
Dolzhenko Sergey Valerjevitch
postgraduate student
Попиков Пётр Иванович
д.т.н., профессор
Воронежская государственная лесотехническая
академия, Воронеж, Россия
Popikov Peter Ivanovich
Dr.Sci.Tech., professor
Voronezh State Forestry Engineering Academy,
Voronezh, Russia
Зубков Алексей Владимирович
аспирант
Международный институт компьютерных
технологий, Воронеж, Россия
Zubkov Aleksey Vladimirovich
postgraduate student
International Institute of Computer Technologies,
Voronezh, Russia
Рыкованова Екатерина Павловна
студентка
Воронежский государственный университет,
Воронеж, Россия
Ryikovanova Ekaterina Pavlovna
student
Voronezh State University, Voronezh, Russia
В статье представлена математическая модель
механизма поворота с дополнительным
демпфером
Mathematical model of the slewing mechanism with
additional damper is given in the article
Ключевые слова: ГИДРОМАНИПУЛЯТОР,
ДЕМПФЕР, ГИДРОУДАР, РАСКАЧИВАНИЕ
ГРУЗА
Keywords: HYDRAULIC MANIPULATOR,
DAMPER, HYDRAULIC SHOCK,
COUNTERBALANCE SWAYING
Управление режимами работы гидроманипуляторов, производящееся
гидрораспределителем, в моменты резкой смены режимов приводит к
раскачиванию
груза
и
гидроударам
в
гидравлической
системе
манипулятора [1]. Раскачивание груза затрудняет его позиционирование и
требует дополнительных затрат времени оператора. Проявляющийся
гидроудар
может
вызвать
разрыв
рукавов
высокого
давления
гидравлической системы или отрыв их от мест подсоединения. Одним из
путей устранения гидроударов и уменьшения раскачивания груза является
использование демпфера, встраиваемого в гидросистему манипулятора,
который позволяет сгладить скачки давления [2].
Нами разработан механизм поворота колонны лесозаготовительного
манипулятора
с
дополнительным
http://ej.kubagro.ru/2012/05/pdf/35.pdf
гидромеханическим
демпфером,
Научный журнал КубГАУ, №79(05), 2012 года
2
представленный на рисунке 1 [3]. На колонне 1 стрелового манипулятора,
установленной на раме 2 лесозаготовительной машины, закреплена
шестерня 3, находящаяся в зацеплении с зубчатой рейкой 4, имеющей на
концах поршни гидроцилиндров 5 и 6. Поршни имеют демпфирующие
устройства 7 и 8. В крышках 9 и 10 гидроцилиндров выполнены гнезда
демпфирующих
устройств 7 и 8 и соединенные каналами линий 11 и 12 с
гидрораспределителем
13.
Однако
эти
демпфируюшие
устройства
срабатывают только в крайних положениях поршней.
Рисунок 1 – Расчетная схема демпфера и гидроцилиндров поворота колонны
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
Поэтому
к
линиям
11
и
3
12
подключен
дополнительный
гидромеханический демпфер 14 через обратные клапаны 15 и 16 дроссели
17 и 18. В корпусе демпфера 14 размещен ступенчатый плунжер 19, в
концевых частях которого выполнены дросселирующие каналы с
калиброванными отверстиями.
Гидромеханический демпфер гасит колебания рабочей жидкости при
остановках колонны в промежуточных положениях. Например, при
повороте колонны вправо в режиме «остановка» гидролинии 11 и 12
становятся запертыми, в гидроцилиндре 7 в полости «Л» возрастает
давление рабочей жидкости за счет инерционных сил, возникающих в
колонне манипулятора, в захвате которого имеется пачка бревен.
Жидкость из полости «Л» по гидролинии 11, через обратный клапан 15
поступает в полость «Л2» демпфера 14, перемещая плунжер 19 вправо.
При этом жидкость из полости «П2» через дроссельное отверстие в
плунжере и регулируемый дроссель 18 поступает в полость «П»
гидроцилиндра 6, при этом всплеск давления рабочей жидкости гасится.
Для
проверки
эффективности
демпфера
и
определения
его
оптимальных параметров была разработана имитационная компьютерная
модель гидроманипулятора с демпфером, встроенным в гидросистему
механизма поворота колонны. В качестве объекта исследования в данной
работе был принят лесной гидроманипулятор ЛВ-184А-06, который
серийно выпускается Майкопским машиностроительным заводом.
В
рамках
происходящие
модели
в
воспроизводятся
механической
и
физические
гидравлической
процессы,
подсистемах
манипулятора, оснащенного дополнительным демпфером. Устройство
описывается системой дифференциальных и алгебраических уравнений.
В модели учитываются три механических процесса: вращательное
движение
колонны
и
стреловой
группы
манипулятора
вокруг
вертикальной оси (рисунок 2), поступательное движение плунжера вдоль
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
4
оси демпфера (рисунок 1), а также раскачивание груза G относительно
точки крепления на стреле манипулятора [2, 4]. Для описания данных
процессов используются уравнения классической динамики [5, 6].
MG
Y
G
φ
LG
MГЦ
O
X
Рисунок 2 – Расчетная схема гидроманипулятора
Поворотная колонна манипулятора представляется в виде абсолютно
твердого
тела,
вращающегося
в
горизонтальной
плоскости
XOY
относительно вертикальной оси О. Угловое положение колонны задается
углом φ, измеряемым от направления OY против хода часовой стрелки.
Для описания поворота колонны используется основное уравнение
динамики вращательного движения:
d 2ϕ
J 2 = ∑ Mi ,
dt
i
(1)
где J – момент инерции колонны относительно оси О;
Mi – моменты сил различной природы.
Описание
движения
плунжера
демпфера
под
воздействием
нескомпенсированных давлений в целом базируется на использовании
второго закона Ньютона:
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
d 2xД
= ∑ Fi ,
dt 2
i
где mД – масса плунжера;
mД
5
(2)
xД – положение плунжера в демпфере (рисунок 1);
Fi – некоторые силы, действующие на плунжер.
В рамках данной модели гидравлическая система манипулятора,
оснащенного демпфером, представляется в виде шести отдельных
полостей, содержащих рабочую жидкость: полостей левого и правого
поворотного гидроцилиндров (обозначены буквами "Л" и "П" на рисунке
1); полостей сброса рабочей жидкости демпфера "Л1" и "П1"; запираемых
полостей демпфера "Л2" и "П2";
При перемещении поршней гидроцилиндров или плунжера демпфера
изменяются объемы Vm соответствующих полостей (m означает индекс
полости). Это приводит к изменению давлений Pm в полостях, причем
данные изменения связаны зависимостью [7]:
dPm
E
=−
,
dVm
Vm
(3)
где E – объемный модуль упругости рабочей жидкости.
Если давления в двух полостях, соединенных между собой,
различаются, начинается перетекание рабочей жидкости, при этом расход
Qij определяется по известной формуле:
Qij = kij sign( Pi − Pj ) Pi − Pj ,
(4)
где i и j – индексы полостей;
kij – коэффициент дросселирования;
sign(x) – функция, возвращающая знак переменной x.
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
6
Эта формула используется как для дросселей (коэффициент
дросселирования
достаточно
велик),
так
и
для
трубопроводов
(коэффициент дросселирования мал).
В модели считается, что все дросселирующие отверстия демпфера
имеют
круглое
сечение,
поэтому
коэффициент
дросселирования
определяется через диаметр отверстия dij по формуле [5]:
kij = µ
πdij2
4
2
ρ ,
(5)
где µ – коэффициент расхода;
ρ – плотность рабочей жидкости.
Возможность трубопроводов упруго расширяться под влиянием
давления в модели непосредственно не учитывается, однако косвенно
учитывается упругостью рабочей жидкости, то есть коэффициентом E.
Основную
моделировании
вычислительную
манипулятора
сложность
при
представляют
компьютерном
дифференциальные
уравнения (1), (2), а также множество уравнений вида (3). Ориентируясь на
использование
вычислительных
возможностей
компьютера,
расчет
организован итерационным образом [8]. В каждом шаге просчитываются
элементарные изменения системы - элементарный поворот колонны
манипулятора,
элементарное
элементарное
перемещение
перемещение
груза.
плунжера
Распишем
ниже,
демпфера,
в
какой
последовательности производится расчет параметров системы на каждой
итерации.
Расчет начинается с того, что по текущему значению угла φ поворота
колонны определяется положение поршня xГ в левом гидроцилиндре:
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
7
x Г = x Г 0 + ϕ ⋅ RК ,
(6)
где xГ0 – положение поршня при φ = 0;
RК – радиус шестерни реечной передачи вращения колонны.
Далее, зная положение xГ поршня левого гидроцилиндра и
положение xД плунжера демпфера, рассчитываем объемы полостей
гидроцилиндра (VЛ, VП) и демпфера (VЛ1, VЛ2, VП1, VП2):
VЛ = x Г
πDГ2
4
;
(7)
V П = (L Г − x Г )
VЛ 1 = x Д
πD Г2
(
;
4
π D 2Д − d П2
)
4
(
VЛ 2 = x Д + LКП − LЦ
(10)
(
VП1 = L Д − x Д − H
(8)
;
(9)
)πd4П
) (
2
;
π D 2Д − d П2
4
)
;
(11)
(
VП 2 = L Д + LКП − x Д − H − LЦ
)πd4П
2
,
(12)
где DГ – внутренний диаметр гидроцилиндра;
LГ – длина рабочей полости гидроцилиндра;
LКП – длина запираемых полостей демпфера;
LЦ – длина концевых частей демпфера;
DД – внутренний диаметр демпфера;
dП – диаметр запираемых полостей демпфера;
LД – длина основной полости демпфера;
H – ширина средней части плунжера.
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
8
Изменение объемов полостей при перемещениях поршней и
плунжера приводит к изменению давлений в полостях. Новые давления PЛ,
PП, PЛ1, PЛ2, PП1, PП2, PАЛ, PАП на k-м шаге интегрирования вычисляются по
формуле (3), переписанной в конечных разностях следующим образом:
Pmk
=
Vmk
k −1
Pm − E
− Vmk −1
,
Vmk
(13)
где индекс m означает полость, в которой вычисляется давление, и
может принимать значения "Л", "П", "Л1", "Л2", "П1", "П2".
После расчета новых давлений в полостях некоторые давления,
будучи довольно малыми на предыдущем шаге, могут на шаге k стать
отрицательными. В этом случае производится их корректировка, то есть,
k
k
если Pm < 0 , то производится присваивание Pm = 0 .
Далее в компьютерном расчете производится учет перетекания
жидкости
из
одной
полости
в
другую
под
влиянием
разности
соответствующих давлений. С учетом гидравлической схемы системы
(рисунок 1) возможны следующие варианты перетекания жидкости:
- перетекание "Л" → "Л1":
VЛ = VЛ − k ЛЛ1 PЛ − PЛ 1 ∆t ;
если PЛ > PЛ 1 , то 
VЛ 1 = VЛ 1 + k ЛЛ1 PЛ − PЛ 1 ∆t ;
- перетекание "Л1" → "Л":
VЛ = VЛ + k ЛЛ1 PЛ 1 − PЛ ∆t ;
если PЛ 1 > PЛ , то 
VЛ 1 = VЛ 1 − k ЛЛ1 PЛ 1 − PЛ ∆t ;
- перетекание "Л" → "Л2":
VЛ = VЛ − k ЛЛ 2 PЛ − PЛ 2 ∆t;
если PЛ > PЛ 2 , то 
VЛ 2 = VЛ 2 + k ЛЛ 2 PЛ − PЛ 2 ∆t;
- перетекание "Л1" → "Л2":
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
9
VЛ 1 = VЛ 1 − k Л 1Л 2 PЛ 1 − PЛ 2 ∆t ;
если PЛ 1 > PЛ 2 , то 
VЛ 2 = VЛ 2 + k Л 1Л 2 PЛ 1 − PЛ 2 ∆t ;
- перетекание "Л2" → "Л1":
VЛ 1 = VЛ 1 + k Л 2 Л 1 PЛ 2 − PЛ 1 ∆t ;
если PЛ 2 > PЛ 1 , то 
VЛ 2 = VЛ 2 − k Л 2 Л 1 PЛ 2 − PЛ 1 ∆t ;
- перетекание "П" → "П1":
VП = VП − k ПП1 PП − PП1 ∆t ;
если PП > PП1 , то 
VП1 = VП1 + k ПП1 PП − PП1 ∆t;
- перетекание "П1" → "П":
VП = VП + k ПП1 PП1 − PП ∆t ;
P
>
P
если П1
П , то V
 П1 = VП1 − k ПП1 PП1 − PП ∆t;
- перетекание "П" → "П2":
VП = VП − k ПП 2 PП − PП 2 ∆t ;
P
>
P
если П
П 2 , то V
 П 2 = VП 2 + k ПП 2 PП − PП 2 ∆t ;
- перетекание "П1" → "П2":
VП1 = VП1 − k П1П 2 PП1 − PП 2 ∆t ;
P
>
P
если П1
П 2 , то V
 П 2 = VП 2 + k П1П 2 PП1 − PП 2 ∆t ;
- перетекание "П2" → "П1":
VП1 = VП1 + k П 2 П1 PП 2 − PП1 ∆t ;
если PП 2 > PП1 , то 
VП 2 = VП 2 − k П 2 П1 PП 2 − PП1 ∆t.
Рассчитанные таким образом объемы перетекающей жидкости затем
корректируются с учетом пропускной способности соответствующих
трубопроводов.
В режимах "поворот колонны влево" и "поворот колонны вправо"
необходимо учитывать не только перетекание жидкости из полости в
полость, но и поступление жидкости от гидронасоса и слив жидкости в
сливную магистраль.
В режиме "поворот колонны влево" предварительно рассчитываются
предполагаемые расходы жидкости:
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
-
поступление
10
"ГН"
→
"Л":
если
PГН > PЛ ,
то
"ГН"
→
"Л1":
если
PГН > PЛ 1 ,
то
"ГН"
→
"Л2":
если
PГН > PЛ 2 ,
то
QГНЛ = k ГНЛ PГН − PЛ ;
-
поступление
QГНЛ 1 = k ГНЛ 1 PГН − PЛ 1 ;
-
поступление
QГНЛ 2 = k ГНЛ 2 PГН − PЛ 2 .
Если суммарный расход QГНЛ + QГНЛ1 + QГНЛ2 превышает
номинальный расход гидронасоса Qном, производится корректировка
расходов:
QГНЛ = QГНЛ
QГНЛ
Qном
;
+ QГНЛ1 + QГНЛ 2
(14)
QГНЛ1 = QГНЛ1
QГНЛ
Qном
;
+ QГНЛ1 + QГНЛ 2
QГНЛ
Qном
.
+ QГНЛ1 + QГНЛ 2
(15)
QГНЛ 2 = QГНЛ 2
(16)
После корректировки вычисляются новые объемы жидкости в
полостях:
- если PГН > PЛ , то VЛ = VЛ + QГНЛ ;
- если PГН > PЛ 1 , то V Л 1 = V Л 1 + Q ГНЛ1 ;
- если PГН > PЛ 2 , то VЛ 2 = VЛ 2 + QГНЛ 2 .
В то время как в "левые" полости гидроцилиндра и демпфера
гидронасос подает жидкость, из "правых" полостей жидкость поступает в
сливную магистраль:
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
11
- слив "П" → "А": если PП > PА , то VП = VП − k ПА PП − PА ∆t ;
- слив "П1" → "А": если PП1 > PА , то VП1 = VП1 − k П1А PП1 − PА ∆t.
Аналогично производятся расчеты в режиме "поворот колонны
вправо".
Предварительно
рассчитываются
предполагаемые
расходы
жидкости:
-
поступление
"ГН"
→
"П":
если
PГН > PП ,
то
QГНП = k ГНП PГН − PП ;
-
поступление
"ГН"
→
"П1":
если
PГН > PП1 ,
то
"ГН"
→
"П2":
если
PГН > PП 2 ,
то
QГНП1 = k ГНП1 PГН − PП1 ;
-
поступление
QГНП 2 = k ГНП 2 PГН − PП 2 .
Если QГНП + QГНП1 + QГНП2 > Qном, то:
QГНП = QГНП
QГНП
Qном
;
+ QГНП1 + QГНП 2
(17)
QГНП1 = QГНП1
QГНП
Qном
;
+ QГНП1 + QГНП 2
(18)
QГНП 2 = QГНП 2
Qном
.
QГНП + QГНП1 + QГНП 2
(19)
После корректировки вычисляются новые объемы жидкости:
- если PГН > PП , то VП = VП + QГНП ;
- если PГН > PП1 , то VП1 = VП1 + QГНП1;
- если PГН > PП 2 , то VП 2 = VП 2 + QГНП 2 .
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
12
Слив жидкости из "левых" полостей учитывается в расчете
следующим образом:
- слив "Л" → "А": если PЛ > PА , то VЛ = VЛ − k ЛА PЛ − PА ∆t ;
- слив "Л1" → "А": если PЛ 1 > PА , то VЛ 1 = VЛ 1 − k Л 1 А PЛ 1 − PА ∆t.
Полученные текущие значения объемов жидкости VЛ, VП, VЛ1, VЛ2,
VП1, VП2, VАЛ, VАП в полостях используются на следующем шаге
интегрирования k + 1 для нового расчета давлений Pm по формулам (4).
По известным давлениям в полостях демпфера рассчитываются
силы, действующие на плунжер, и затем рассчитывается новое положение
и
скорость
плунжера.
Уравнение
движения
плунжера
с
учетом
определенных давлений может быть записано следующим образом:
(
)
0, x Д ≥ l П
d 2 xД
π D Д2 − d П2
1 
πd П2
=
P
+
P
+
c
Л1
Л2
П
dt 2
m Д 
4
4
l П − x Д , x Д < l П
π D2Д − dП2
0, LД − x Д − H ≥ lП
πdП2
(
− PП1
4
)−P
П2

−


dxД 
− cП ⋅ 
−
kД ,

l
−
L
+
x
+
H
,
L
−
x
−
H
<
l
4
dt
П
Д
Д
Д
Д
П


(20)
где kД – коэффициент вязкого трения плунжера при перемещениях в
демпфере;
сП – жесткость внутренней пружины-упора;
lП – свободная длина внутренней пружины-упора.
Для интегрирования уравнения (20) используется численный метод –
модифицированный метод Эйлера:
&x&kД
(
π
− PПk1−1
(D
)
2
2
2

 P k −1 π D Д − d П + P k −1 πd П + c
Л2
П
 Л1
4
4

1
=
mД
2
Д
− d П2
4
)−P
0, x kД−1 ≥ l П
⋅
k −1
k −1
l П − x Д , x Д < l П

−



0, LД − xkД−1 − H ≥ lП
k −1
&
− cП ⋅ 
− x Д k Д  ; (21)
k −1
k −1

4
lП − LД + x Д + H, LД − x Д − H < lП

2
k −1 πd П
П2
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
x kД
=
x kД−1 + x& kД−1∆t +
&x&kД (∆t )2
2
13
;
x& kД = x& kД−1 + &x&kД ∆t ,
(22)
(23)
где k и k–1– текущий и предыдущий шаги интегрирования
соответственно.
По сравнению с базовым методом Эйлера, имеющим первый
порядок погрешности, данный метод имеет третий порядок погрешности
для координаты и второй для скорости, а сама вычислительная схема (22,
23) является эффективной и устойчивой [8].
Основное уравнение вращательного движения для колонны имеет
вид:
d 2ϕ
dϕ
J 2 = M ГЦ − M G − M ТР − kTК
+ MУ + M В ,
dt
dt
(24)
где J – момент инерции колонны;
MГЦ – момент, сообщаемый поворотной колонне гидроцилиндром;
MG – момент со стороны приводимого в движение груза;
MТР – момент сил трения в подшипниках колонны;
kТК – приведенный коэффициент вязкого трения в гидроцилиндрах
поворота;
MУ – момент сил от уклона местности;
МВ – момент сил от ветровой нагрузки.
Для расчета момента инерции колонны, последнюю можно считать
телом, состоящим из двух частей более-менее правильной геометрической
формы: стрелы, которую можно рассматривать как однородный стержень
массой mС и длиной LG, и основания колонны, которое можно упрощенно
представить в виде сплошного цилиндра, массой mОК и радиусом RОК. С
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
14
учетом введенных допущений момент инерции рассчитывается по
известным формулам:
1
1
2
J = mC L2G + mОК RОК
.
3
2
Крутящий
момент
со
стороны
(25)
гидроцилиндров
поворота
рассчитывается по формуле:
M ГЦ = FГЦ ⋅ RК ,
(26)
где сила со стороны гидроцилиндров поворота FГЦ рассчитывается
по формуле:
FГЦ = (PЛ − PП )
πD Г2
4
.
(27)
Момент со стороны груза MG рассчитывается следующим образом:
M G = LG (− FGX ⋅ sin ϕ + FGY ⋅ cosϕ ) ,
(28)
где FGX и FGY – декартовы составляющие силы, оказываемой на
стрелу со стороны груза.
Таким образом, окончательное уравнение вращательного движения
колонны можно записать следующим образом.
d 2ϕ
dt2
=
2


1
 RК (PЛ − PП ) πDГ − LG (− FGX ⋅ sinϕ + FGY ⋅ cosϕ) − MТР − kTК dϕ + MУ + MВ  .

1
1
4
dt
2 


mC L2G + mОКRОК
3
2
(29)
В процессе численного интегрирования используется конечноразностная схема уравнения (29), согласно модифицированному методу
Эйлера:
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
ϕ&&k =
(
)
15
(
)
2


1
k −1
k −1
k −1
 RК PЛk −1 − PПk −1 πDГ − LG FGY
⋅ cosϕ k −1 − FGX
⋅ sinϕ k −1 − MТР
− kTКϕ& k −1 + MУk −1 + M Вk −1  ;

1
1
4
2 


mC L2G + mОК RОК
3
2
(30)
ϕ =ϕ
k
k −1
+ ϕ&
k −1
∆t +
ϕ&&k (∆t )2
2
;
(31)
ϕ& k = ϕ& k −1 + ϕ&&k ∆t .
(32)
Груз в модели упрощенно представляется в виде материальной точки
массой mМ. Груз (точка M) взаимодействует со стрелой (точка G)
посредством
невесомого
вязкоупругого
стержня,
имитирующего
устройство захвата. В соответствии с законом динамики поступательного
движения можно записать векторное уравнение движения груза:
d 2 rM
1
=
2
mM
dt
где


  − cM (GM − LM ) − d M d (GM − LM )  GM + mM gr  ,


dt
 GM


(33)
rM – радиус-вектор точки М в декартовой системе координат
XYZ;
сМ и dM – коэффициенты жесткости и вязкости вязкоупругого
взаимодействия посредством устройства захвата;
GM – вектор, исходящий из точки G и оканчивающийся в точке M;
GM – расстояние между точками G и M;
r
g – вектор ускорения свободного падения.
В последнем уравнении выражение, стоящее во внешних скобках,
представляет собой силу, действующую на груз. Противоположная ей сила
действует со стороны груза на стрелу:
r
d (GM − LM )  GM

FG =  cM (GM − LM ) + d M
− mM g .

(34)
dt

 GM
Составляющие последней силы FGX и FGY участвуют в приведенном
выше уравнении вращательного движения поворотной колонны.
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
16
Таким образом, общая система уравнений, описывающая работу
манипулятора в целом, имеет следующий вид:
(
)
2
2
d 2 x
0, xД ≥ lП
πdП2
1  π DД − dП
 Д=
P
P
c
+
+
⋅
−

Л
Л
П
1
2
4
4
 dt2 mД 
lП − xД , xД < lП

dxД 
π D2Д − dП2
0, LД − xД − H ≥ lП

πdП2
P
P
c
kД ;
−
−
−
⋅
−

П2
П
 П1

4
4
lП − LД + xД + H, LД − xД − H < lП dt


d 2ϕ
2


1
 RК (PЛ − PП ) πDГ − LG (− FGX ⋅ sinϕ + FGY ⋅ cosϕ) − MТР − kTК dϕ + MУ + MВ ;
 2 =

1
1
dt
4
2 
 dt


mC L2G + mОКRОК

3
2
d 2 r
F
 M =− G ;
2
mM
 dt

d(GM− LM )  GM
r
FG = cM (GM− LM ) + dM
− mM g.

dt

 GM

(
)
Для решения системы дифференциальных уравнений, положенной в
основу
модели,
экспериментов
"Программа
с
для
и
для
проведения
моделью
составлена
моделирования
различных
компьютерных
компьютерная
гидроманипулятора,
программа
оснащенного
демпфером гидросистемы" на языке Object Pascal в интегрированной среде
программирования Borland Delphi 7.0. Получено свидетельство об
официальной регистрации программы для ЭВМ № 2009610503. Решение
системы уравнений производится путем численного интегрирования, при
этом
выводятся
временные
зависимости
основных
параметров,
характеризующих работу манипулятора и демпфера: давлений в полостях
поворотных гидроцилиндров PП(t) и PЛ(t), положения плунжера демпфера
xД(t), угла поворота колонны φ(t), тангенциальных и радиальных
колебаний груза ∆lτ(t), ∆lr(t).
Список литературы
1.
Хуако З.А. Исследование влияния раскачивания груза при вращении колонны на
производительность и динамическую нагруженность механизмов манипулятора
// 70 лет кафедре механизации лесного хозяйства и проектирования машин
Воронежской государственной лесотехнической академии: Межвуз. сб. науч. тр.
– Воронеж, 2007. – С. 230.
http://ej.kubagro.ru/2012/05/pdf/35.pdf
Научный журнал КубГАУ, №79(05), 2012 года
2.
3.
4.
5.
6.
7.
8.
17
Бартенев И.М. Гидроманипуляторы и лесное технологическое оборудование /
И.М. Бартенев, З.К. Емтыль, А.П. Татаренко, М.В. Драпалюк, П.И. Попиков,
Л.Д. Бухтояров. – М.: ФЛИНТА: Наука, 2011. – 408 с.
А. с. 1792910 СССР, М. Кл. В 66 С 13/42. Механизм поворота колонны
стрелового манипулятора [Текст] / А.П. Нестеров, П.И. Попиков, В.В. Волынко;
(СССР). - № 4835090/29; заяв. 13.04.90; опубл. 07.02.93, Бюл. № 5. - 5 с.: ил.\
Советов Б.Я., Яковлев С.А. Моделирование систем: учебное пособие – М.:
Высш. шк., 1998. – 319 с.
Моделирование сельскохозяйственных агрегатов и их систем управления: учеб.
для вузов / под ред. А. Б. Лурье. – Л.: Колос. Ленингр. отд-ние, 1979. – 312 с.
Расчет и проектирование строительных и дорожных машин на ЭВМ /Под ред.
Е.Ю. Малиновского. – М.: Машиностроение, 1980. – 216 с.
Элементы гидропривода (Справочник). Изд. 2-е, перераб. и доп. /Е.И. Абрамов,
К.А. Колесниченко, В.Т. Маслов. – Киев: Техника, 1977. – 320 с.
Инженерные расчеты на ЭВМ: Справочное пособие / Под ред. В.А. Троицкого. –
Л.: Машиностроение, 1979. – 288 с.
http://ej.kubagro.ru/2012/05/pdf/35.pdf
1/--страниц
Пожаловаться на содержимое документа