УДК 622.232.32
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАБОЧЕГО ПРОЦЕССА ГИДРАВЛИЧЕСКОГО УДАРНОГО МЕХАНИЗМА
Устименко Т.А., канд.техн.наук, Кононенко А.П., докт.техн.наук, Селивра С.А., канд..техн.наук, Яценко А.Ф.,канд..техн.наук, Донецкий национальный технический университет
Источник: Наукові праці Донецького національного технічного університету. Серія «Гірничо-електромеханічна». Випуск 16(142).-Донецьк: ДонНТУ, 2008 с.258-264(фахове видання)
Исследована одна из принципиальных схем гидравлического ударного механизма, составлена математическая модель рабочего процесса, содержащая дифференциальные уравнения движения основных звеньев. Путем сравнения данных, полученных в результате моделирования и натурного эксперимента, доказана адекватность разработанной математической модели .
One of the percussion mechanisms schematic circuit was reseached, mathematical model was created. It includes differentials equations main elements movement. It was proved model adequacy to the full-scale experiment.
Проблема и ее связь с научными или практическими задачами. Разработка эффективного инструмента для механизации трудоемких операций по разрушению массива в строительстве, металлургии, угольной промышленности является актуальной задачей. Особенно перспективным является гидравлический инструмент. Если гидравлические механизмы (в основном импортного производства) находят достаточно широкое применение в качестве навесного оборудования, то разработка ручного инструмента является достаточно перспективным направлением.
Анализ исследований и публикаций. Известны несколько подходов по классификации, выделению и исследованию различных структурных схем гидравлических механизмов ударного действия [1,2]. Исследованию работы гидравлических ударных механизмов посвящены работы, в которых рассматриваются устройства c различной управляемостью камер. Однако в данных работах не рассматривается динамика управляющего звена совместно с движением узла ударника.
Постановка задачи. Гидравлический ударный механизм отбойного молотка, работающий в системе объемного гидропривода представляет собой сложную гидродинамическую систему, которая характеризуется большим числом (более 20) факторов. Для всестороннего исследования ее функционирования, выбора основных конструктивных и рабочих параметров, для выявления рациональных режимов, необходимо проведение большого числа экспериментов. Причем проведение каждого опыта будет сопряжено с изготовлением отдельных узлов и деталей устройства. Наиболее эффективным путем решения является математическое, а точнее, имитационное моделирование рабочего процесса установки. В связи с этим ставится задача математического описания рабочего процесса одной из принципиальных схем гидроударного механизма, с целью разработки имитационного алгоритма и программы. Данный подход позволит исследовать большой объем вариантов и произвести выбор оптимальных параметров системы без проведения трудоемких экспериментов.
Изложение материала и результаты. Рассматриваемая принципиальная схема гидроударного механизма приведена на рис.1. На схеме указаны основные узлы устройства: I – узел ударника с рабочим инструментом 1 и поршнем-бойком 2; II – узел управления с поршнем-клапаном 3, а также гидропневмоаккумулятор (ГПА) гидроударного механизма 4 и элементы гидроэнергетического снабжения: насос объемного принципа действия 6 с предохранительным клапаном 7 и емкостью для рабочей жидкости 8, ГПА 9 и напорной магистралью 5, т.е. основные структурные звенья системы, математические модели которых будут построены. При разработке математической модели были приняты следующие допущения:
во входном сечении напорной магистрали предполагается постоянный расход рабочей жидкости;
массы подвижных элементов (поршня-бойка и поршня-клапана) предполагаются сосредоточенными;
волновые процессы в соединительных каналах не учитываются, гидравлические сопротивления каналов считаются сосредоточенными у соответствующих камер;
изменения значений давления и объема воздуха в ГПА подчиняется адиабатному закону;
утечки жидкости в системе отсутствуют;
волновые процессы, происходящие в рабочем инструменте и бойке не учитываются.
На рис.1 указаны направления расходов, выбранные для математического описания процесса. Основным является узел ударника, поршень-боек которого совершает возвратно-поступательное движение под действием сил давления на его рабочие поверхности. Уравнение движения бойка имеет вид:
(1)
где
маса
поршня-бойка; х
–
координата положения поршня-бойка положение;
- давления жидкости в соответствующих рабочих камерах;
- площади рабочих поверхностей поршня-бойка в соответствующих
камерах;
- сила трения, создаваемая резиновыми уплотнениями (кольцами)
поршня-бойка, рассчитывается по формуле:
.
(2)
В соответствии с классической теорией удара скорость после соударения с инструментом определяется формулой
при
,
где
- коэффициент восстановления скорости, который может находиться в
пределах [3] 0,1 …0,3.
Подвижным элементом узла распределителя является поршень-клапан, уравнение движения которого имеет вид:
(3)
Рис. 1. Расчетная гидравлическая схема ударного механизма
здесь
- масса поршня-клапана; z
– координата поршня;
- давления жидкости в соответствующих рабочих камерах;
- площади рабочих поверхностей;
– сила трения, которая рассчитывается по формуле:
При
достижении поршнем-клапаном одного из крайних положений (,
скорость его гасится до нуля.
Давление в камере D определяется давлением ГПА и в соответствии с допущением описывается уравнением адиабаты:
, (5)
начальные
давление и объем закачки газа в ГПА;
-
текущие значения давления и объема ГПА в рабочем процессе; k=1.4
-коэффициент адиабаты. Текущий объем определяется объемом рабочей
жидкости, поступающим и отдаваемым ГПА:
(6)
-
расход рабочей жидкости, поступающий в ГПА.
Давления в рабочих камерах рассчитываются по формулам:
,
– давления в соответствующих рабочих камерах. Расходы
определяются скоростью движения поршня-бойка:
.
– гидравлические сопротивления каналов, соединяющих
соответствующие камеры.
- переменное гидравлическое сопротивление между седлом и клапаном,
рассчитываемое по формуле [ 4]:
где
-
диаметр клапана,
- расстояние от клапана до седла (величина открытия клапана),
- переменный коэффициент расхода, зависящий от положения клапана,
может быть представлен в соответствии с данными [5] таблично заданной
функцией (см. табл.1), позволяющей применить линейную интерполяцию
между соседними координатами.
Табл.1.
|
0,71 |
0,78 |
0,89 |
0,913 |
0,9 |
0,87 |
0,73 |
0,65 |
|
0, 1 |
0, 2 |
0, 4 |
0, 6 |
0, 8 |
1 |
2 |
3 |
|
0, 6 |
0, 56 |
0, 53 |
0, 515 |
0, 5 |
0, 485 |
0, 459 |
0, 431 |
|
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
Давления в камерах клапана-распределителя:
Здесь
- гидравлическое сопротивление канала ВЕ,
– переменное гидравлическое сопротивление между головкой
клапана и седлом со стороны сбросной магистрали, рассчитываемое по
формуле
,
– сосредоточенное сопротивление сбросной магистрали.
– расход, определяемый скоростью перемещения клапана
распределителя,
- расход, поступающий в сбросную магистраль. Давление в камере С
можно выразить и как
.
Связь расходов в рассматриваемой системе (см.рис.1) выражается следующими уравнениями баланса:
-
расход жидкости, поступающий в молоток
равен расходу, потребляемому молотком
;
-
баланс расходов относительно
можно выразить как
.
Здесь
- расход рабочей жидкости, подаваемый к гидроударному механизму из
напорной магистрали (задается исходя их характеристики питающего
насоса).
Таким образом, для определения состояния системы необходимо решение дифференциальных уравнений движения поршня-бойка (1) и клапана (2) с учетом функциональных коэффициентов, задаваемых алгебраическими и интегральными зависимостями (2,4,5,6). Полученная система уравнений не может быть решена аналитически из-за существенной нелинейности граничных условий, поэтому решение осуществлялось численно путем имитационного моделирования рабочего процесса. С целью проверки адекватности разработанной математической модели было проведено сравнение осциллограмм рабочего процесса, полученных опытным путем на натурном образце гидромолотка и в результате моделирования на компьютере. На рис.2 показаны графики зависимости скорости движения бойка для натурного и машинного экспериментов. Сравнение предударной скорости бойка показало, максимальное отклонение составило 15.6 % , а среднеквадратическое отклонение
Рис.2. К проверке адекватности математической модели
не превышает 4%. Это свидетельствует об адекватности имитационной модели натурной с точностью, достаточной для прикладных исследований.
Направление дальнейших исследований. Полученные результаты дают возможность проведения исследований как ручных, так и более мощных механизмов, и использованием построенной математической модели, реализованной в виде имитационного алгоритма, адекватность которой доказана на основании сравнения результатов моделирования и натурного эксперимента. Разработанная имитационная модель может быть использована для оптимизации разрабатываемых конструкций гидравлических ударных механизмов.
Список источников:
Соколинский Б.В. Машины ударного разрушения (Основы комплексного проектирования).-М.: Машиностроение, 1982.-184с.
2. Тимошенко Г.М., Устименко Т.А., Мартыненко В.Ф.Гидравлический отбойный молоток ГМ-9// Уголь Украины,- 1988.-№12.-с.24-26.
Селивра С.А. Разработка гидравлических устройств ударного действия для разрушения горных пород. Дисс. … канд. техн. наук:05.05.06.-Донецк,1986.-194 с.
Гидравлика, гидромашины и гидроприводы: Учебник для машиностроительных вузов/Т.М. Башта, С.С.Руднев, Б.Б. Некрасов и др. – М.: Машиностроение, 1982.-423с.
Караев М.А. Гидравлика буровых насосов. – М.: Недра, 1983.-208 с.