УДК 62-82.003.13
А. И. Павлов, И. А. Полянин, П. Ю. Лощёнов
ТЕОРЕТИЧЕСКИЕ ИССЛЕДОВАНИЯ ДИНАМИЧЕСКИХ СВОЙСТВ ЭЛЕМЕНТОВ ГИДРОПРИВОДА ЛЕСНЫХ МАШИН ПРИ ДИАГНОСТИРОВАНИИ В ФУНКЦИОНАЛЬНОМ РЕЖИМЕ
Представлено теоретическое исследование динамических свойств элементов гидропривода лесных машин при диагностировании в функциональном режиме, которое позволяет обосновать возможность определения технического состояния элементов гидропривода лесных машин.
Ключевые слова: математическая модель, теоретические исследования, гидроударное воздействие, податливость, гидравлическое сопротивление, дифференциальные уравнения, гидропривод, величина повышения давления, диагностирование.
Введение. Гидроприводы лесных машин являются сложными динамическими системами, подверженными воздействию широкого спектра динамических нагрузок переменного характера и работающими в условиях изменяющихся внешних воздействий, обусловленных разнообразными факторами. При эксплуатации эти нагрузки приводят элементы гидропривода к неисправностям. Для предотвращения этих неисправностей рекомендуется применять техническое диагностирование, с помощью которого определяется их текущее техническое состояние и остаточный ресурс.
В настоящее время получило развитие направление диагностирования гидроприводов как в функциональном режиме, так и в процессе технического обслуживания и ремонта, что позволяет изучить причины потери работоспособности, определить степень удаленности от предельного состояния, установить вид и место возникновения дефектов.
Регулярно проводимое диагностирование способствует своевременному выявлению нарушений в техническом состоянии гидроприводов и их устранению, поэтому тема по обоснованию, разработке методов и средств диагностирования гидроприводов лесных машин является весьма актуальной.
Цель работы - теоретическое исследование динамических свойств элементов гидропривода лесных машин при диагностировании в функциональном режиме.
Решаемые задачи:
• разработка математической модели гидропривода лесных машин для исследования его динамических свойств в функциональном режиме;
• исследование упругих свойств элементов гидроприводов лесных машин в функциональном режиме с помощью гидроударного воздействия;
• определение диагностического параметра при гидроударном воздействии на элементы гидропривода.
Для построения математической модели гидропривода используется известный унифицированный метод математического описания механизмов при помощи гидромеханических цепей, позволяющий исследовать влияние гидроударного воздействия на элементы гидропривода.
© Павлов А. И., Полянин И. А., Лощёнов П. Ю., 2012.
В качестве расчётной примем структурную гидромеханическую схему гидропривода стрелы манипулятора марки ОМТЛ-70-02, который содержит в своей структуре элементы с нелинейными характеристиками, оказывающие существенное влияние на динамические процессы, происходящие в гидроприводе.
При составлении гидромеханической цепи предполагается, что каждый конструктивный элемент гидропривода представляется несколькими элементами цепи. Кроме того, каждый элемент цепи представляет лишь одно физическое явление - гидравлическое сопротивление, сжимаемость, источник расхода, источник давления, трение, инерционность, силу и т.д., при этом каждый элемент цепи описывается только одним членом дифференциального уравнения.
Гидромеханическая цепь гидропривода стрелы манипулятора показана на рис. 1.
NK1
R
тр1
72
R
тр1
72
R
шл1
"72
72
R.
тр2
72
^тр2
72
Дщл2 72 R
Г тр2
1 72
шл2
72
Рис. 1. Гидромеханическая цепь гидропривода стрелы манипулятора марки ОМТЛ-70-02 с запорным клапаном
Влияние параметров состояния элементов гидропривода на диагностический параметр предлагается исследовать в переходном режиме, создавая при этом гидроударные воздействия за счет резкого закрытия и открытия запорного клапана.
При составлении цепи приняты следующие нижеописанные допущения и ограничения:
V
7
8
9
1) гидроцилиндр нагружен силой Рпр, направленной в сторону поршневой полости и приведенной к его штоку, которая принимается фиксированной;
2) масса тпр груза и манипулятора является сосредоточенной и приведена к штоку
гидроцилиндра, которая также принимается фиксированной;
3) рукава высокого давления характеризуются гидравлическими сопротивлениями Rшл и объемными деформациями (податливостями) Сшл, причем сопротивления напорной и сливной магистралей приняты равными, а податливости - различными;
4) полости гидроцилиндра характеризуются активными площадями F1, F2 и податливостями С1, С2;
5) сопротивление трения поршня о гильзу цилиндра выражено через коэффициент нелинейного трения п;
6) перетечки масла из напорной полости в сливную оцениваются гидравлическими сопротивлениями Rпер, утечки - гидравлическими сопротивлениями - Rу .
7) гидрораспределитель представлен нелинейными гидравлическими сопротивлениями N и N2, которые в силу конструктивных особенностей равны N1 = N2 = N ;
8) давление Рн перед гидрораспределителем в исходном состоянии принимается
равным номинальному;
9) запорный клапан и обратный клапан представлены нелинейными гидравлическими сопротивлениями Nкл , Nобкл, соответственно.
Гидромеханические цепи удобно анализировать на основе правила контуров или правила узлов. Так как цепь включает гидроцилиндр с различными активными площадями полостей и нелинейные элементы (золотниковый гидрораспределитель и рукава высокого давления), в этом случае целесообразно для составления системы уравнений воспользоваться правилом узлов.
Правилу узлов соответствуют уравнения
¿QI = о, £ ц = о, (1)
1=1 1=1
то есть алгебраическая сумма расходов либо скоростей (ц), жидкостей, притекающих к узлу и отводимых от узла, равна нулю. Правило узлов для гидравлической части цепи представляет собой форму уравнения расходов, для механической части цепи - условие неразрывности (неразрушения) совместно движущихся элементов конструкции [1].
При анализе методом узлов неизвестными являются разности давлений между выбранным опорным узлом и другими узлами цепи.
Разработанная гидромеханическая цепь содержит одну «длинную» ветвь, которая учитывается как один элемент, и 10 узлов, поэтому для цепи необходимо составить одно уравнение «длинной» ветви и 10 уравнений узлов.
Уравнение «длинной» ветви
Рб ■Fl - т*§ ~п' - Рпр - Fl ■ р?=а (2)
Уравнение узла 1
Л1р - Р1 - Т^л/Р1 - Р2 = 0. (3)
^об.кл
Уравнение узла 2
¿V Р1 " Р2 Р2 - Рз =
Уравнение узла 3
л/2 I-
л/Р2 - Рз -
R
шл1
RL1 + dt
Уравнение узла 4
2 (Рз - Р4 У^^Р^ - ^ ^ = °. (6)
II
Ri 1 + ^U4"3 Rmp1v"4 " mp1 dt
Rmpr N
Уравнение узла 6
5 ^ Л/6 ^пер ^ ° * Л
Уравнение узла 7
Р2 -и+ТТ- (Рб - Р7 - Р 8 - С2 = °.
Кпер Ктр2 Л
Уравнение узла 8
' ЛР8
V2 I-
Л/Ру " Р8 -
Rmp2
R2mp2 + Rl, 8 ^ ^Р2 dt
Уравнение узла 9
2 (Р8 - Р9 )-^л/Р9 - Р10 - Сшл2 = °. (11)
\
Rmp2 + Rшл2 Rшл2 dt
Уравнение узла 1°
(4)
2 (Рз-Р<)-Сш1^ = 0. (5)
Уравнение узла 5
^"л/Р4 - Р5 ^Т^л/Р5 - Рб = °. (7)
^ л/Рз - Рб Рб (Рб - Р7)-C % - F-и = 0. (8)
(9)
(Р8 - Р9 )-Cmp2^ = 0. (10)
л/Р9 - Р10 = 0. (12)
л * 9 10 #
При расчётах принимаем, что характер движения масла в трубопроводах ламинарный, а в золотниковых гидрораспределителях, запорном клапане и различных местных сопротивлениях - турбулентный.
Сопротивления трубопроводов и рукавов высокого давления определяются по формуле [2,4]:
1тр I
Nтр = 8 - Ктр • Рж 2 тРг5 ^ Nшл = 8 ' Кшл ' Рж 2 ^5 , (13)
^ - Лтр П ■ Лшл
где 1тр, Лтр, 1шл и Лшл - длина и внутренний диаметр >го трубопровода и РВД, м; Хтр и Хшл - коэффициент Дарси для трубопроводов и рукавов.
Нелинейное гидравлическое сопротивление золотника гидрораспределителя определяется по формуле [3]:
— = К ■ Ь ■ к, (14)
N
К = и ^, (15)
V У
где коэффициент расхода,и = 0,65-0,72;Ь - длина щели, м; h- ширина щели, м; у - кинематическая вязкость, м2/с.
Гидравлические сопротивления утечек и перетечек из полостей гидроцилиндра будут определяться как потери давления при истечении масла через кольцевую щель (посадочную или уплотнительную) по формуле [3]:
я, = , (16)
ж-и-о
где у - кинематическая вязкость, м2/с; р - плотность масла, кг/м3; I - длина уплотнения, м; и - диаметр уплотнения, м; 0 - осевой зазор, м.
Площадь поршневой полости гидроцилиндра можно определить по формуле:
«Д- • (17)
где и - диаметр цилиндра, м.
Коэффициент податливости поршневой полости гидроцилиндра, заполненной рабочей жидкостью, определяется по формуле [1]:
'1 - х и д2- хЛ
С = ж
8 E4-S4 4-E^
(18)
где х - коэффициент Пуассона, х = 0,35; Д1 - внутренний диаметр поршневой полости гидроцилиндра, м; Ец, Еж - модули упругости материала стенок и жидкости (для минеральных масел Еж = (1,4...1,9)-109Па); Sц - толщина стенки гильзы гидроцилиндра, м; х1 - перемещение штока гидроцилиндра, м.
Площадь штоковой полости гидроцилиндра можно определить по формуле:
= -д2 - (19)
2 4 4
где и - диаметр цилиндра, м; - - диаметр штока, м.
Коэффициент податливости штоковой полости гидроцилиндра определяется по формуле [1]:
С2 =Ж
1 - х D3 + D12-(Hmm - х) ^
----1--
, 8 E4-S4 4-Еж ,
(20)
где х - коэффициент Пуассона, х = 0,35; Д1 - внутренний диаметр штоковой полости гидроцилиндра, м; Ец, Еж - модули упругости материала стенок и жидкости, Па; Sц -толщина стенки гильзы гидроцилиндра, м; х1 - перемещение штока гидроцилиндра, м; Ншт - ход штока гидроцилиндра, м.
Коэффициенты податливости рукавов высокого давления, заполненных маслом, определяются по формуле [1]:
С = — -шл 1шл (21)
^шл ~ . ^ , ( )
4-Е
пр
где Епр - приведенный модуль упругости, учитывающий сжимаемость масла и деформацию рукавов, Епр = (1...1,1) - 108Па.
Влияние параметров состояния элементов гидропривода на диагностический параметр предлагается исследовать в переходном режиме, создавая при этом гидроударные
воздействия за счет резкого закрытия и открытия запорного клапана.
Гидроударное воздействие - это кратковременное, но резкое повышение давления в трубопроводе при внезапном торможении двигавшегося по нему потока жидкости в результате мгновенного закрытия и открытия запорного устройства.
Поскольку гидросистема обладает определенной упругостью, то появляющийся при гидроударном воздействии избыток расхода в результате сжимаемости рабочей жидкости и деформации стенок трубопровода компенсируется упругостью гидросистемы:
а -
i
N
-4P3 - сл
об.кл
Фз dt
- Qi = 0.
(22)
Определение параметров диагностирования при гидроударных воздействиях гидропривода в условиях переходного процесса производится с учетом следующих допущений:
1) участок от узла 3 и до полости ^ напорной гидролинии представляем источником расхода Q1. Расход на сжимаемость масла и деформацию участка от узла 3 до полости ^ напорной гидролинии определяется по уравнению [1]:
dp3
Qi = С
dt
(23)
где Сшл - объемные деформации (податливость) рукавов высокого давления, заполненных маслом, которые определяются по формуле (21);
2) время срабатывания предохранительного клапана при гидроударном воздействии больше времени протекания гидроударного воздействия.
С учетом принятых допущений на рис. 2 приведена упрощенная схема гидромеханической цепи напорной гидролинии гидропривода стрелы манипулятора ОМТЛ-70-02 при гидроударном воздействии. Нумерация узлов сохранена в соответствии с рис. 1.
Величина повышения давления при гидроударном воздействии определяется по формуле Жуковского для полного прямого гидравлического удара, так как параметры гидроударного воздействия идентичны параметрам полного прямого гидравлического удара:
Ар = ри0С, (24)
где р - плотность жидкости, кг/м3; о0 - скорость движения жидкости в гидроприводе до гидроудара, которая определяется из выражения, м/с; Су - скорость распространения ударной волны, м/с.
Скорость распространения ударной волны опре-_| деляется по формуле [3]:
Рис. 2. Гидромеханическая цепь напорной гидролинии гидропривода стрелы манипулятора марки МТЛ-70-02 при гидроударном воздействии
C, =
¿ж Р
1 +
D E.,,
(25)
5 E
тр
где
E
Р
- скорость звука в жидкой среде, м/с; Б - диаметр трубопровода, м; 8 -
толщина стенки трубопровода, м; Еж - модуль упругости жидкости, Па; Етр - модуль упругости трубопровода, Па.
Величину повышения давления жидкости с затухающими колебаниями, возникающую при гидроударном воздействии в напорной гидролинии подъема стрелы манипулятора, построим в программе MathCAD по формуле:
x(t) = p0 + Ap • e~Елин (t} • cos с (t), (26)
где х - динамическая составляющая знакопеременной нагрузки; p0 - давление жидкости в напорной гидролинии при подъеме стрелы, Н/м2; Ap - величина повышения давления при гидроударном воздействии, Н/м2; Елин - модуль упругости материала напорной гидролинии подъема стрелы, Па; с = - круговая частота колебаний, с-1; T -
время одного цикла, включающего повышение и понижение давления (фаза гидроударного воздействия), с.
Считая скорость ударной волны при повышении и понижении давления одинаковой, определим фазу гидроударного воздействия [5]:
Т = *,
(27)
где О - скорость распространения ударной волны, м/с; I - длина трубопровода, при которой возможно гидроударное воздействие, м.
Длина трубопровода, при которой возможно гидроударное воздействие, определяется по формуле:
C ■ t l = C
2
(28)
где - время закрытия запорного клапана, с.
При диагностировании напорной гидролинии с помощью гидроударного воздействия в условиях переходного процесса для наиболее точного определения результатов необходимо учитывать следующие требования:
- подъем стрелы при каждом диагностировании необходимо осуществлять при полном вылете манипулятора без груза;
- подъем стрелы необходимо осуществлять при фиксированном угле между стрелой и рукоятью манипулятора;
x1(t) x2(t)
3x10
7
2.5x10",
2x10
1.5x10
1x10
0 0.02 0.04 0.06 0.08 0.1 t
Рис. 3. Зависимости давления от времени при гидроударном воздействии в напорных гидролиниях в момент подъема стрелы манипулятора без нагрузки: / - время, с; х1(1) - давление в гидролинии, не имеющей наработки, Па; х2ф - давление в гидролинии, имеющей предельную наработку, Па
- гидроударное воздействие необходимо осуществлять до момента пока гидроцилиндр стрелы манипулятора не достигнет своего среднего положения от начала подъема.
Полученные в программе MathCAD зависимости давления от времени при гидроударных воздействиях в напорных гидролиниях, не имеющей наработки и имеющей предельную наработку в момент подъема стрелы манипулятора без нагрузки, показаны на рис. 3.
В процессе динамического нагружения гидролинии происходит ее старение, при этом на рабочих поверхностях образуются микротрещины, нарушаются межмолекулярные связи в резине, что приводит к уменьшению модуля упругости гидролинии до предельного значения, при котором она выходит из строя.
Анализ зависимостей показывает, что изменение упругих свойств элементов гидропривода оказывает значительное влияние на колебательные процессы в гидроприводе. Уменьшение коэффициента податливости гидролинии, то есть потеря упругих свойств материала, приводит к уменьшению величины гидроудара и увеличению длительности колебательного процесса при появлении и развитии неисправностей (трещин, микроразрывов и т.д.) в напорной гидролинии, поэтому в качестве диагностических параметров примем величину повышения давления, возникающую при гидроударном воздействии.
Выводы
1. Разработанная математическая модель гидропривода лесных машин и ее анализ позволили обосновать возможность определения технического состояния элементов гидропривода лесных машин в функциональном режиме.
2. Изменение упругих свойств элементов гидропривода оказывает значительное влияние на колебательные процессы в гидроприводе. Уменьшение коэффициента податливости гидролинии, то есть потеря упругих свойств материала приводит к уменьшению величины гидроудара и увеличению длительности колебательного процесса при появлении и развитии неисправностей в напорной гидролинии.
3. Увеличение гидравлического сопротивления элементов гидропривода за счет появления микротрещин, микроразрывов, увеличения шероховатости внутренних рабочих полостей также приводит к уменьшению величины гидроудара.
Таким образом, выбранный нами диагностический параметр является чувствительным и достаточным при проверке работоспособности гидролинии.
Список литературы
1. Хандрос, А.Х. Динамика и моделирование гидроприводов станков /А.Х. Хандрос, Е.Г. Молчанов-ский. - М.: Машиностроение, 1969. - 156 с.
2. Справочное пособие по гидравлике, гидромашинам и гидроприводам / Я.М. Вильнер, Я.Т. Ковалев, Б.Б. Некрасов, Ю.А. Беленков, Ю.Л. Кирилловский. Под. общ. ред. Б.Б. Некрасова. - Мн.: Выш. шк., 1985. - 382 с.: ил.
3. Павлов, А.И. Повышение надежности гидроприводов лесных машин: дис. ... д-ра техн. наук: 05.21.01: защищена 28.12.04: утверждена 5.05.2005/Павлов Александр Иванович. -Йошкар-Ола, 2004.397 с. - Библиогр.: с 292 - 296.
4. Лебедев, Н.И. Гидропривод машин в лесозаготовительной промышленности / Н.И.Лебедев. - М.: Лесная промышленность, 1978. -304 с.
5. Гамынин, Н.С. Гидравлический привод систем управления /Н.С.Гамынин. - М.: Машиностроение, 1972. - 376 с.
Статья поступила в редакцию 01.11.11.
A. I. Pavlov, I. A. Polyanin, P. Yu. Loshchenov
THEORETICAL RESEARCH OF DYNAMIC PROPERTIES OF ELEMENTS OF HYDRAULIC ACTUATOR OF WOOD CARS WHEN DIAGNOSING AT THE FUNCTIONAL MODE
A theoretical research of dynamic properties of the elements of a hydraulic actuator of wood cars when diagnosing in a functional mode is presented. The diagnosis allows to prove the possibility of definition of a technical condition of the elements of a hydraulic actuator of wood cars.
Key words: mathematical model, theoretical research, hydroshock exposure, pliability, hydraulic resistance, differential equations, hydraulic actuator, pressure boost size, diagnosis.
ПАВЛОВ Александр Иванович - доктор технических наук, профессор, проректор по аккредитации, лицензированию и внешним связям МарГТУ. Область научных интересов -повышение надежности гидроприводов лесных машин. Автор более 50 публикаций.
E-mail: [email protected]
ПОЛЯНИН Игорь Александрович - доктор технических наук, профессор, заведующий кафедрой транспортно-технологических машин МарГТУ. Область научных интересов -гидравлика и гидропривод, лесотранспортные машины. Автор более 60 публикаций.
E-mail:[email protected]
ЛОЩЁНОВ Павел Юрьевич - аспирант кафедры транспортно-технологических машин МарГТУ. Область научных интересов - повышение надежности гидроприводов лесных машин. Автор четырех публикаций.
E-mail: [email protected]