Динамика машинного агрегата с упругим валом и квадратической характеристикой исполнительного механизма

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Динамика машинного агрегата с упругим валом и квадратической характеристикой исполнительного механизма

Ст.Н. Бъчваров

В.Д. Златанов

С.Г. Делчева-Атанасова

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

Как было отмечено в [2], переходные процессы (разбег и торможение) мало исследованы, при этом механические системы (модели) обычно рассматриваются отдельно от привода и характера рабочих нагрузок.

Цель исследования

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

Динамическая модель

Динамическая модель машины (рис.1), как и в [2], состоит из асинхронного электродвигателя (М) с постоянным моментом инерции ; зубчатой передачи (G), имеющей постоянный (приведенный к валу электродвигателя) момент инерции и постоянное передаточное отношение ; сцепной муфты () между двигателем и зубчатой передачей с постоянным моментом трения при буксовании в зацепленном положении; постоянной муфтой () между зубчатой передачей (G) и исполнительным механизмом (W) ( она постоянно включена и поэтому не изображена на рис. 1).



Рис. 1. Модель агрегата

В настоящем исследовании принимается, что исполнительный механизм (W) - это механизм ротационого типа с постоянным моментом инерции и квадратической механической характеристикой , где - момент сил сопротивления исполнительного механизма; - постоянная величина, - угловая скорость ротора исполнительного механизма.

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

. (1)

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

, (2)

где постоянные величины и определенны техническими данными электродвигателя.

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

Движение агрегата происходит в два этапа:

1) этап пробуксовки муфты, во время которой угловые скорости валов сцепной муфты стремятся к выравнению;

2) этап совместного раскручивания, на котором оба вала ускоряются совместно до установления рабочей (стационарной) угловой скорости агрегата.

Дифференциальные уравнения движения на этапе пробуксовки муфты

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

(3)

При введении новых переменных и :

, (4)

, (5)

система (3) без первого уравнения записывается в виде

(6)

при этом

, (7)

. (8)

Определение закона движения на этапе пробуксовки муфты

Первое уравнение системы (3) представляем в виде

, (9)

где .

Зададим начальные условия движения:

(10)

Учитывая (10) из уравнении (9), находим

, (11)

. (12)

Решение системы дифференциальных уравнений (6) будем искать методами последовательных приближений и малого параметра.

Рассмотрим сначала первое уравнение системы (6). Исходя из (8), для находим:

. (13)

В первом приближении принимаем

, (14)

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

, (15)

.

Решение этого дифференциального уравнения, учитывая начальные условия (10), находим в форме[6]:

. (16)

Подставляя (14) во второе уравнение системы (6) и учитывая (16), имеем

, (17)

.

Решение дифференциального уравнения (17) с учетом нулевых начальных условий (10) выражается интегралом Дюамеля

, (18)

,.

Таким образом, решение (17) можно представить как сумму двух интегралов

, (19)

. (20)

Интеграл (20) вычисляем следующим способом. Функцию представляем в виде ряда [3]:

. (21)

Ряд (21) сходится, если выполнено условие , откуда

,

где - время пробуксовки муфты.

Сохраняя первые три члена разложения (21), находим

. (22)

Четная функция , представлена рядом (22), удовлетворяет условиям Дирихле и может быть представлена в виде суммы ряда Фурье [3,5]:

, (23)

причем

.

Здесь

а



.

Получаем выражение

(24)

позволяющее найти приближенное значение интеграла (20):

где n - число учтенных членов разложения функции в ряд Фурье.

Таким образом, решение дифференциального уравнения (17), примет вид

, (25)

откуда

, (26)

,

,

, .

В первом приближении для угловых скоростей валов , учитывая (7), (8), (16), (26), имеем

(27)

(28)

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

(29)

причем , - малый положительный параметр, для которого принято

,

,

,

а коэффициенты получены из равенства

.

Уравнение (29) таким образом приведено в квазилинейном виде.

Решение дифференциального уравнения (29) ищем в виде ряда по степеням малого параметра с точностью до членов, содержащих в первой степени включительно:

. (30)

Подставляя соотношение (30) в уравнение (29) и приравнивая коэффициенты при одинаковых (нулевой и первой) степенях малого параметра в правой и левой частях уравнения, получаем

, (31)

(32)

Уравнение (31) запишем в виде

,

.

Решение его аналогично решению уравнения (15):

. (33)

Введем обозначения:

,

и сделаем замену переменной

. (34)

Тогда дифференциальное уравнение (33) примет вид

,

а его решение будет [4]:

, (35)

причем

.

Учитывая начальные условия (10), а также (30), (33), (34), находим во втором приближении

, (36)

Второе уравнение системы (6), учитывая (13), записываем в виде:

, (37)

где - малый положительный параметр, для которого принято , , .

Решение диференциального уравнения (37) ищем в виде ряда по степеням малого параметра с точностью до членов, содержащих в первой степени включительно:

. (38)

Подставляя соотношение (38) в уравнение (37) и приравнивая коэффициенты при одинаковых (нулевой и первой) степенях малого параметра в правой и левой частях уравнения, получаем

, (39)

. (40)

Решение уравнения (39) при начальных условиях (10):

. (41)

Подставляя решение (41) в правую часть уравнения (40), имеем:

. (42)

Решение дифференциального уравнения (42) при нулевых начальных условиях (10) может быть представлено в виде интеграла Дюамеля

, (43)

, .

Решение (43) можно представить как сумму интегралов:

, (44)

причем

,

,

где - коэффициенты разложения в ряд Фурье нечётной функции (21). Интеграл вычислен методом, аналогичным использованному при вычислении интеграла (20), при этом для определения , найдено аналитическое выражение интегралов вида:

(45)

Формулы (38)-(45) определяют и во втором приближении.

Функции во втором приближении определены, а это позволило определить изменение угловых скоростей согласно формулам (7), (8), а также крутящий момент на упругом валу.

Пример

Полученные результаты проиллюстрированы численным примером на основе технической характеристики многофункционального прессавтомата с перенастройкой MPA-1PM, производства ТИКЕ-АД, Пловдив, Болгария. Для автомата: , ,,,. Для асинхронного электродвигателя с короткозамкнутым ротором AM-132-M-6: , и . Коэффициент момента сопротивления , а жесткость вала .

Изменение угловых скоростей и нагрузки зубчатой передачи, определены получеными в работе формулами, показаны на рис. 2-5. Вместе с тем, с целью сравнения, построены и графики этих величин, определенные численным интегрированием системы (3) с помощью многошагового метода Хеминга четвертого порядка. Из приложенных рисунков видно хорошее совпадение результатов, полученных двумя методами.

Рис. 2. Изменение угловых скоростей



Рис. 3. Изменение угловых скоростей

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Рис. 4. Изменение крутящего момента

Список литературы

машинный агрегат муфта редуктор

1. Боголюбов Н.Н., Митропольский Ю.А. Асимптотические методы в теории нелинейных колебаний. М.: Физматгиз, 1963. - 463 с.

2. Бъчваров Ст.Н., Златанов В.Д., Делчева-Атанасова С.Г., Янчев И.Г. Динамика машинного агрегата с упругим валом и линейной характеристикой исполнительного мехенизма. [Электронный ресурс]: Электрон. статья (1 файл KБ) // Теория механизмов и машин, вып. 1(3) том 2, 2004.-СПб: СПбГПУ, 2004.-Загл. с титул. экрана.-Электрон. журн.-Свободный доступ из сети Интернет.-http://tmm.spbstu.ru/journal.

3. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Наука, 1971.-1108 с.

4. Зайцев В.Ф., Полянин А.Д. Справочник по нелинейным дифференциальным уравнениям (приложение в механике, точные решения). М.: Наука, 1993.-462 с.

5. Корн Г., Т. Корн. Справочник по математике для научных работников и инженеров. М.:-Наука, 1973.-832 с.

6. Лойцянский Л.Г., Лурье А.И. Курс теоретической механики, часть ІІ. М.:ОГИЗ, 1948.-580 с.

Размещено на Allbest.ru

...