УДК 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 с.