Влияние динамических процессов, формирующихся в рабочих режимах, на силовые параметры ножевого исполнительного органа геохода Текст научной статьи по специальности «Механика и машиностроение»

© В.В. Аксенов, В.Ю. Садовец, В.Ю. Бегляков, 2009

УДК 622.232

В.В. Аксенов, В.Ю. Садовец, В.Ю. Бегляков

ВЛИЯНИЕ ДИНАМИЧЕСКИХ ПРОЦЕССОВ, ФОРМИРУЮЩИХСЯ В РАБОЧИХ РЕЖИМАХ, НА СИЛОВЫЕ ПАРАМЕТРЫ НОЖЕВОГО ИСПОЛНИТЕЛЬНОГО ОРГАНА ГЕОХОДА

Исследовано влияние динамических эффектов в приводе на силовые параметры исполнительного органа геохода. Определен коэффициент запаса. Ключевые слова: геоход, диагностика, динамический эффект, ножевой исполнительный орган.

ш ж ривод геохода содержит к силовых домкратов, каждый Л. Л. из которых обеспечивает поворот агрегата на угол

уравновешенность привода и наличие кинематических особенностей (для краткости именуемых дефектами) приводит к формированию в спектре колебаний с высшими частотами тО и пО . Для определенности считаем m < п; m, п = 2,3,....

Диагностика особенностей режимов и распознавание моделей нагружения агрегата возможны, если дополнительные колебания имеют устойчивый гармонический характер, а следовательно, их изображения по Лапласу содержат простые полюса 5 = ±шО и 5 = ±шО . Как показано в [1], существует конечное число моделей систем, передаточные функции которых удовлетворяют этому условию. При квазистатическом воздействии силового домкрата

V) = V • 1(/), где 1(/) - ступенчатая функция Хэвисайда, подобных моделей существует всего 11. Таким образом, проблема моделирования сводится к перебору различимых особенностей или видов состояния систем привода и выбору в качестве расчетного худшего из вариантов.

Сравнение удобно производить с ординарной системой, имеющей 2 степени свободы, передаточная функция которой

создавая среднюю угловую скорость

) = (2^)4

1

1

(1)

2 'Г'2 л 2 2 2 'Г'2 л 2 2

s Т + 4^ т s Т + 4^ п При параллельном соединении элементарных звеньев (по определению, данному в [1]) получаем

1 1 2s

---------------+----------------_ ““5—2-----------------------------------7 (2)

sT - / • 2жт sT + / • 2жт s Т + (2пт)

1 1 2s

или

sT -1 • 2пп sT +1 • 2пп s2Т2 + (2пп)2 На этой основе вводятся следующие модели систем привода

4 ^ 1 1

(3)

У2(э) = (2я) • 2^Т•

У3(У) = (2п)2 • ^Т)

^2Т2 + (2пш)2 s2Т2 + (2жп)2 1 1

(4)

(5)

s2T2 + (2пт)2 s2T2 + (2пп)2

Простое параллельное соединение подсистем приводит к новым вариантам в зависимости от знака связи

УЛ) = (2ж)2

1

1

(6)

s 2T2 + (2кт)2 s2 + (2пп)2

Вводя дополнительно такие же соединения звеньев в подсистемах, получаем две новые группы моделей

1 2sT

У5Лз) = (2л)

У6/ С*) = (2^)2

^ 2Т2 + (2пш)2 s 2Т2 + (2пп)2 2^Т 1

У7/ (^) = (2^)2 • 2^Т

^ 2Т2 + (2пш)2 s 2Т2 + (2пп)2

11

s Т + (2жш) s + (2пп)

(7)

(8) (9)

Таким образом, наряду с ординарным решением (рис.1)

11

V

Ух(г) = 0

22

тп

1 —^—- (п2 cos тО! - т2 cos пО!)

п2 - т2

(10)

т = 2, п = 3

Рис. 1. Изменения окружной скорости на выходе основной динамической системы

Получаем для сравнения совокупность моделей, проиллюстрированную наиболее интересными примерами на рис. 2 при гармониках 2О и 3О.

Изменение амплитудных значений окружной скорости в зависимости от состояния приводов достигает 3-х порядков, что доказывает чрезвычайную важность рассмотрения всей совокупности альтернативных моделей. Понятно, что прочностной расчет должен гарантировать работу привода при модели №3 и, кроме того, зна-

чительный интерес представляют особенности состояния при моделях №6 и №7.

Напротив, при практическом анализе ряд моделей можно отфильтровать. Однако, чтобы избежать пропуска важных

V2(t)

Vз(t)

vлt)

)

Рис. 2. Примеры возможных вариантов динамических процессов при наличии

особенностей (дефектов) в приводе геохода

Рис. 3. Примеры вариантов динамических процессов, формирующихся на оси геохода

эффектов окончательное заключение делается по моделям проекций окружной скорости на координатные оси машины

Ух ^) = V )cos УУ ^) = V ^іп Qt (11)

Здесь следует обратить внимание на то, что, по сути, происходит амплитудная модуляция, приводящая к изменению состава спектра - теперь он содержит гармоники (т ± 1)О и (п ± 1)О .

Это обстоятельство оказывается существенное влияние на трактовку результатов при различных способах измерения колебаний машин.

В качестве примера приведем уравнения координатных составляющих скорости колебаний при модели №3:

cos(m - 1)О! + cos(т + 1)О! -

- cos(n - 1)О! - cos(n + 1)О!

Vз.x (!) = 8^0

3. Х\Ч 2 2

п - т

Ъ, (!) = 8^0

3!\Ч 2 2

п -т

(12)

(13)

sin(m - 1)О! - sin(m + 1)О! -- sin(n - 1)О! + sin(n + 1)О!

Для того, чтобы подчеркнуть эффект изменения частотного состава иллюстрации на рис. 3 даны при т = 3; п = 5 .

Представленный материал позволяет проиллюстрировать многообразие форм колебаний координатных составляющих скорости на оси машины. При наибольшей опасности модели №3 с позиций диагностики большой интерес могут представлять и другие особенности состояния привода.

Одно из основных утверждений заключается в том, что, используя полную группу моделей динамических систем, мы гарантируем представительность выбранных для последующих прочностных расчетов типов динамических процессов. Кроме того, разработанный способ позволяет заранее получить правила проведения и трактовки результатов вибродиагностики без поиска аналогий для оригинальной по конструкции машины.

В развитие задачи моделирования кинематики (координатных составляющих скорости), в которой определена полная группа и представлено многообразие уравнений колебаний для выбора типажа при прочностных расчетах, следует рассмотреть возможные варианты видов движения, которые могут возникать на выходе привода геохода, содержащего дефекты. Эти результаты должны использоваться для описания траекторий движения резцов, отклоняющихся от круговой траектории под действием циклически включающихся силовых домкратов.

Модели кинематики позволяют получить уравнения переносного движения инструмента, отстоящего на радиус R от оси выработки.

Хп () = Ух (Г¥п (Г) = V (Г)R (14)

Можно привести примеры существенно различающихся типов переносных колебаний, возникающих при характерных видах состояния (рис. 4).

№2

т = 3, п = 4

№3

№5

№10

2

0

-2

50

0

-50

10

о

-10

20

о

-20

/\ % / ЧУ С' л / V/ /'Ч. ' \А7 У

0 0.5 т = 3, п = 4 1.5

^ /'\ ЧУ X /Х_,-\ / ~ - /\ Ч./ “

0 0.5 т = 3, п = 4 1.5

Рч /'\ * \У * / \ 'у--, '*'? % / чв/ \У С ч

0 0.5 т = 3, п = 4 1.5

0.5

1.5

Рис. 4. Иллюстрации возможных типов переносных колебаний рабочего органа

Необходимо подчеркнуть, что:

1) модели формирования колебаний в приводе являются между собой линейно независимыми;

2) на графиках видны очевидные амплитудные отличия;

3) даже визуально (например, по положению экстремумов) отмечается разнообразие форм колебаний.

Таким образом, обосновывается постановка задачи о необходимости перебора всех возможных вариантов состояний для получения надежных результатов расчетов.

Для моделирования относительных колебаний машины в плоскости, перпендикулярной к оси выработки, используем простое определение

т = 2, п = 3

1 Ч \ 1 1 I 1 \

0 0.5 1 1.5 т = 2, п = 4

—-\ ^

ч .х / / \ --Л

,%>Г^ :

0 0.5 1 1.5 т = 3, п = 4

| I I

О 0.5 1 1.5

Рис. 5. Примеры наиболее опасных видов колебаний по модели №3

X(Г) = |Ух (3^3; Y(Г) = |у7 (3^3 (15)

0 0 Как и ранее, наиболее опасный вид движения возникает в случае состояния и модели №3. Уравнения колебаний по координатным осям можно получить достаточно просто

4ж¥Т Хъ«) = 4^°Т

3У"' 2 2

п - т

1 . 2п(т - Ш 1

---------------1-------sin

т -1

1 . 2п(п -1)7

- sin

п -1

Т

п +1

Т т +1 2п(п +1У

Т

(16)

1

4яУ Т І 1

гг(і) = 4^°Т 1 1

п2 - т2 І т -1

1 - cos

2л(т -1)? Т

+ -

т +1

1 - cos

2л(т +1)?

Т

п -1

1 - cos

2п(п -1)? Т

п +1

1 - cos

2п(п +1)? Т

(17)

В результате можно выделить наиболее опасные виды колебаний (рис. 5), разнообразные как по амплитуде, так и по форме. Теперь, если образовать из этих координатных составляющих

Рис. 6. Траектории оси машины при относительных колебаниях для модели №3 состояния привода

соответствующую суперпозицию, то можно получить модели траекторий, по которым движется ось машины относительно оси выработки (рис. 6) под действием циклически работающих силовых домкратов.

Главную особенность траекторий составляют особые точки, в

которых имеют место разрывы производной ск / или, иначе, рез/ сх

1

1

1

кие изменения направления движения резцов. Понятно, что реально проявляется угловая суперпозиция переносного и относительного движений, и при больших радиусах особые точки проявляются гораздо слабее.

Можно отметить, что уравнение, описывающее колебания при модели №6, отличается от рассмотренного смешанным характером учета высших гармоник тО и пО, а также присутствием в явном виде гармоники О (рис. 7)

т = 2, п = 3

т = 2, п = 4

Рис. 7. Примеры колебаний машины по модели №6

Рис. 8. Траектории оси машины при относительных колебаниях для моделей №6 и №10 состояния привода

2VT

‘) = 0

1

т I т -1

1 - cos

2ж(т -1)1

Т

+ -

1

т +1

1 - cos

2ж(т +1)1

Т

+

кт

4жп

. . 2nt 1 . 2л (п -1^ 1 . 2л(п +1^

2 эт--------------------------------------------sln-sm

Т п -1

Т

п+1

Т

2КТ Y6(t) = К

т

1 . 2п(т - Ш 1 . 2ж(т + 1У

-эт—1----- -------эт-

т-1

+

КТ

4жп2

2nt Л 21 1 - соэ------1 +

Т

1

п+1

Т ) п -1 2к(п + 1)t

1 - соэ

т+1 Т

2к(п - 1)t

Т

1 - соэ

Т

+

+

(18)

(19)

1

Здесь выделяется совпадение форм колебаний и их независимость от параметра п . Таким образом, получаются два типа одинаковых траекторий (рис. 8).

Наконец, модель формирования №7 приводит к обратному по отношению к №3 типу уравнений движения

2VT

х7(;) = 0

+-

т

п -1

тп

1 - cos

1 - cos

2л(т - 1)і

т -1 2л (п -1)7

Т

Т

+

+

т

п+1

1 - cos

т +1 2л (п +1)7

Т

1 - cos

2л (т +1)7

Т

+

(20)

2УТ ^(7) = ^

тп

п

т -1

біп

2л(т -1)7

п

+ -

т п -1

біп

2л(п -1)7

Т

т

Т т +1 2л(п +1)7

2л(т +1)7

біп—-------— +

Т

п +1

-біп-

Т

(21)

В результате можно построить траектории иного типа (рис. 9), в которых проявляется сильная зависимость формы от параметров высших гармоник. Кроме того, в конечных решениях также можно отметить линейную независимость вариантов между собой.

Обозначив координаты особых точек через х, у, можно использовать простую пропорцию для оценки изменения шага резания на инструменте

п

п

Рис. 9. Траектории оси машины при относительных колебаниях для моделей №7 и №11 состояния привода

ктп^х

....Ах + у

А*—: Р , (22)

2кR

а, следовательно, и приращения усилия резания.

Итак, использование принципов диагностики по отношению к системам привода геоходов позволяет промоделировать значимо различающиеся режимы резания и нагружения рабочего органа.

Геометрические неточности в установке и разброс характеристик домкратов поворота (особенно при их значительном числе) могут привести к тому, что почти ударные воздействия будут перераспределяться в широкой зоне работы ножа исполнительного органа, обеспечивая эффективное разрушение в отдельных, сосредоточенных зонах.

При возникновении динамических эффектов возникают пульсации на траектории движения геохода. Стоит отметить, пульсации на траектории возникают и без динамических эффектов.

В то же время возникающие пульсации могут привести к ухудшению работы ножей исполнительного органа и увеличению необходимых усилий резания, вследствие отклонения углов резания от номинального значения.

Рассмотрим две точки ножа исполнительного органа геохода, которые отстоят друг от друга на расстоянии АД , что соответствует величине смещения траектории движения, вызываемые динамическими погрешностями привода [2].

Угол между плоскостью перпендикулярной оси вращения геохода и передней гранью ножа для первой точки равен

^1 = А + §о, (23)

для второй точки

/2 = А2 + 8о, (24)

причем /1 > /2 (рис. 10).

При возникновении колебаний первая точка ножа исполнительного органа (ИО) геохода смещается на величину АR и попадает на траекторию движения второй точки ножа. Тогда изменение угла резания будет определяться зависимостью:

А§ = §1 — §^ = /1 — /2 . (25)

Подставив в зависимость (25) выражения (23) и (24) , получим:

А§ = |А — А\ = АА . (26)

Возникновение динамических эффектов в приводе геохода приводит не только к изменению траектории движения ножей, но и изменению геометрических параметров резания ножевым исполнительным органом, а, следовательно, и силовых параметров резания. Учет возникающих пульсаций представляется сложным и требует большого объема дополнительных исследований.

Для учета влияния погрешностей, вызванных приводом геохода, на силовые параметры резания введем определяемый коэффициент запаса. При приближении к центру вращения геохода угол наклона винтовой лопасти к оси вращения будет стремиться к 900.

Введем ограничения на угол / < 800 , тогда А/ =/ — §0 .

а

Рис. 10. Геометрические параметры сечений ножа:

циями

а - номинальное положение ножа; б - положение, вызванное пульса-

При угле резания §о = 250 допустимый угол подъема винтовой лопасти, при котором будет происходить эффективное разрушение горного массива, будет равен А/ = 550 .

Радиус, соответствующий углу А/ = 550 , будет определяться исходя из выражения

Го0=/2^Р/, ’ (27)

где Ъв - шаг винтовой лопасти, мм.

Точки ножа исполнительного органа геохода, находящиеся на

окружности радиусом г00 , будут смещаться от номинальной траектории движения геохода на некоторую величину АД . Радиус траектории движения точки ножа, учитывающий возникающие

пульсации, определяется выражением Г = Г00 + АД .

Каждому из выше представленных радиусов соответствуют свои углы подъема винтовой лопасти Г ^ А0 и 0 ^ А00 . Тогда изменение геометрии резания можно записать как

АА = \@00 — А)| или А§ = §00 — §0|.

Согласно методике расчета сил резания по опытным коэффициентам, каждому углу резания соответствует свой коэффициент <Р0, учитывающий свойства разрушаемого горного массива и угол резания. При возникновении погрешности угол резания определится выражением §00 = §0 + АА , и этому углу будет соответствовать свой коэффициент р00 , учитывающий влияния угла резания.

Тогда коэффициент, учитывающий влияние динамических погрешностей привода на силовые параметры ножевого исполнительного органа, назовем его коэффициентом запаса, будет равен

£ = р > 1. (28)

Р0

Динамические эффекты, возникающие в приводе, вызывают изменение траектории движения геохода. Отклонение от заданной

траектории приводит к изменению геометрических параметров резания, что безусловно влияет на силовые параметры резания. Влияние динамических эффектов в приводе геохода на силовые параметры исполнительного органа, учитывается коэффициентом запаса £.

------------------------------------------ СПИСОК ЛИТЕРАТУРЫ

1. Л0г0в А.Б., Замараев Р.Ю. Математические модели диагностики уникальных объектов / Новосибирск, - Издательство СО РАН, 1999. - 228 с.

2. Г0рбун0в В.Ф., Аксен0в В.В., Сад0вец В.Ю. Экспертная оценка влияния особенностей нового класса горнопроходческой техники на методику расчета его параметров. // Вестник КузГТУ - 2004. - № 6.1.- С. 43 - 45. ШХЗ

V. V. Aksenov, V.J. Sadovets, V.J. Begljakov

INFLUENCE OF THE DYNAMIC PROCESSES FORMED IN OPERATING CONDITIONS, ON POWER PARAMETERS OF KNIFE OPERATING BODY OF THE GEOCOURSE

Influence of drive dynamic effects to the executive device force parameters of geohod explored. The coefficient of reserve determinated.

Key words: geocourse, diagnostics, dynamic effect, knife operating body.

— Коротко об авторах ----------------------------------------------------

Аксен0в В.В. - доктор технических наук, ведущий научный сотрудник, Институт угля и углехимии Сибирского отделения РАН (ИУУ СО РАН), [email protected]

Сад0вец В.Ю. - кандидат технических наук, доцент, Государственное образовательное учреждение высшего профессионального образования «Кузбасский государственный технический университет» (ГОУ ВПО КузГТУ), е-так [email protected]

Бегляк0в В.Ю. - старший преподаватель, Юргинский технический институт (филиал) Томского политехнического университета, е-таП: [email protected]_________________________________________