Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

INFLUENCE OF THE ELASTICITY OF THE TRACTION BODY ON THE DYNAMICS OF THE ELECTRIC DRIVE IN TRANSIENT MODES

Shabo K.Ya. 1
1 Technical Istitute (brunch) of «North-Eastern Federal University named after M.K. Ammosov»
The presence of an elastic traction body in the electric drive system on one side significantly worsens the operating mode of the electromechanical installation, and on the other hand it makes a very complicated mathematical model necessary for the optimization of the control system. Describing the mechanical part of the electric drive, it is necessary to consider it as a mass elastic system. Here the greatest difficulties arise in obtaining the analytical dependence of displacements and forces in individual sections of the traction organ. In electric drive systems, in which there are elastic links with a small stiffness coefficient, low-frequency oscillations appear, adversely affecting both the durability of the elements of the system and the reliability of the coupling of the drive drums with the traction organ. Therefore, it is advisable, with the help of the control system, to minimize the oscillation processes as much as possible. To describe the processes in mechanical elastic links with distributed parameters, the simplest and most practical is the use of the transfer functions of elementary two mass links obtained by the method of structural modeling with the help of boundary elastic bonds. The advantage of this approach is proved by comparison with the representation of an elastic link with distributed elasticity and mass of a chain system with a large number of concentrated masses
еlectric drive
elastic connections
a lot of mass management system

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

Как показал опыт многочисленных исследований подобных одноприводных систем, элементарной схемой является схема, представленная на рис. 1, а. Все другие представляют определенное сочетание таких элементарных схем. Для описания динамики таких систем предложен ряд методов.

Рассмотрим некоторые методы и оценим их с точки зрения эффективности их применения для решения задачи оптимизации управления многодвигательным электроприводом.

Метод граничных упругих связей

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

hab1.tif

Рис. 1. Расчетные схемы двухмассовых систем с инертной упругой связью

Тогда усилия, действующие со стороны тягового органа на сосредоточенные массы, будут равны

habo01.wmf (1)

где X1y, X2y – перемещения крайних сечений относительно центра тяжести отрезка тягового органа длиной l1; X0 – перемещение тягового органа как абсолютно твердого тела. Чем больше коэффициент Ск, тем меньше разности X1 – X0 – X1y и X2y + X0 – X2.

Искусственное отделение тягового органа от сосредоточенных масс позволяет рассматривать систему как состоящую из трех взаимосвязанных звеньев с определенными входными и выходными переменными (рис. 2).

Принимая упругие силы F1y и F2y направленными от сосредоточенных масс m1 и m2 звенья 1 и 2 будут описываться уравнениями

habo02.wmf (2)

habo03.wmf (3)

где F = F(t) – динамическая составляющая усилия приложенного в массе m1, F = FДВ – FСТ, FСТ – сила статических сопротивлений движению системы, FДВ – сила создаваемая двигателем. Для описания звена нужно в общем случае рассматривать его как звено с распределенными параметрами (массой упругостью). Наиболее близкой аналогией является упругий стержень со свободными концами, к которым приложены силы F1y
и F2y. В системах такого вида для определения перемещения любого поперечного сечения тягового органа длиной l1 без учета рассеивания энергии необходимо решить волновое уравнение

habo04.wmf (4)

где aK – скорость распространения в продольном направлении упругих деформаций, м/с.

hab2.tif

Рис. 2. Декомпозиция двух массовой упругой системы

Для свободных колебаний в случае одновременного действия сил F1y и F2y изображения перемещений крайних сечений принимают вид [2, 3]:

habo05.wmf (5)

habo06.wmf. (6)

Из равенства (7) следует важный вывод, что если рассматривать колебания относительно перемещения тягового органа как абсолютно твердого тела, описываемого уравнением

habo07.wmf (7)

то можно учитывать только нечетные члены рядов [4, 5].

Тогда уравнения (5) и (6) запишутся так:

habo08.wmf (8)

где habo09.wmf – коэффициент учитывающий моменты инерции движущих звеньев; mK – масса ветви тягового органа длиной l1.

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

habo10.wmf. (9)

Составленная по уравнениям (1), (2), (3), (7) и (9) структурная схема по отношению к динамическим усилиям F1y и F2y, после исключения вспомогательных упругих связей (подставляется Ск = ∞) будет иметь вид, показанный на рис. 3.

Передаточные функции по отношению к F1y и F2y:

habo11.wmf; (10)

habo12.wmf; (11)

где

habo13.wmf; habo14.wmf;

habo15.wmf; habo16.wmf;

μK – коэффициент сопротивления.

Структурная схема определения усилий в переменных состояния приведена на рис. 4.

hab4.tif

Рис. 4. Схема упругой системы в переменных состояния

На этой схеме:

habo17.wmf

habo18.wmf

habo19.wmf.

Переменные состояния системы связаны с величиной Е следующими зависимостями:

habo20.wmf

где XI и X2 – значения переменных на N-ом шаге интегрирования;

XI (0) и X2 (0) – значения этих переменных, полученные на предыдущем шаге.

где habo21.wmf.

hab3.tif

Рис. 3. Структурная схема упругой системы

Укрупненная структурная схема алгоритма определения усилий F1y и F2y приведена на рис. 4.

Ориентировочную оценку переходного процесса можно дать непосредственно по изображениям (10) и (11). Так, кроме частоты колебаний ωK, можно, не прибегая к исследованию процесса, определить начальные и установившиеся значения F1y и F2y.

При p → 0 и t = 0

habo22.wmf

habo23.wmf

При p → 0 и t = 0

habo24.wmf

habo25.wmf

Замещение упругой системы с распределенными параметрами дискретной цепной системой

Как показано в [6], упругую систему с распределенными массой и упругостью можно заменить системой с сосредоточенными массами, соединенными упругими идеальными связями, как показано на рис. 5, а.

hab5a.tif

hab5b.tif

Рис. 5. Расчетная (а) и структурная (б) схемы многомассовой упругой системы

Чем больше число участков N, тем ближе по динамическим свойствам дискретная система системе с распределенными параметрами. Общий коэффициент жесткости упругого тягового органа Сy можно определить по скорости распространения упругой деформации aK:

habo26.wmf. (12)

Разделив участок тягового органа, расположенный между сосредоточенными массами m1 и m2, N частей и приравняв элементарную сосредоточенную массу habo27.wmf, найдем коэффициент жесткости элементарной упругой связи CK из равенства

habo28.wmf откуда CK = N?Cy.

В каждую из масс m1 и m2 включены половины масс mK. Обозначив собственные сосредоточенные массы m10 и m20, будем иметь:

habo29a.wmf и habo29b.wmf.

Следовательно, согласно принятым
выше обозначениям habo30.wmf

Запишем систему уравнений динамики упругой системы (рис. 5, а) при ступенчатом приложении усилия F.

habo31.wmf

habo32.wmf

habo33.wmf

habo34.wmf.

Данным уравнениям соответствует структурная схема (рис. 5, б), которую можно упростить, заменив абсолютные значения сосредоточенных масс относительными величинами (рис. 6). Нормированные коэффициенты a1, a2 и βK равны

habo35.wmf (13)

По полученной нормированной структурной схеме легко составляются уравнения состояния и алгоритм анализа переходных процессов для любого числа участков N. Примем N = 6 и значения параметров системы a 11 = 0,5; a12 = 2; aK = 4000; l1 = 1500; F = 200 кН. Уравнения динамики относительно действительных упругих усилий и фиктивных переменных y1 и dy1/dt = V1:

habo36.wmf (14)

hab6.tif

Рис. 6. Универсальная схема многодвигательной упругой системы

Структурная схема алгоритма определения усилий Fi в функции времени приведена на рис. 7.

Графики F1(t) и F2(t) – рис. 8, а, б. На рисунках для сравнения показаны графики усилий F1y и F2y, полученные по формулам (10) и (11). На рис. 8, в приведены графики F1y(t) и F2y(t) при коэффициенте демпфирования μK, рав-
ном 0,01.

hab7.tif

Рис. 7. Структурная схема алгоритма определения упругих сил в многомассовых системах

hab8.tif

Рис. 8. Графики упругих сил, полученные на ЭВМ

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