Моделирование пуско-тормозных режимов ленточного конвейера

В более ранних работах автора [2] была разработана модель ленточного конвейера. Эта модель представляет собой комплекс из трех математических моделей:

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

Модель асинхронного привода с частотно-векторной системой управления;

Модель вычисления и стабилизации величины тягового фактора конвейера E .

Расчетная схема конвейера представлена на рис. 1.

Здесь лента представлена четырьмя сосредоточенными массами m , m , m , m , соединенными между собой упруго-вязкими связями. Одна масса G представляет собой натяжное устройство, расположенное в хвостовой части конвейера. Именно натяжное устройство, перемещаясь и выбирая слабину, возникающую из-за растяжения ленты, обеспечивает натяжения в ветвях ленты, необходимые для поддержания требуемой величины тягового фактора. В качестве обобщенных переменных приняты перемещения x1, x2, x3, x4, четырех масс, их скорости x1,x2, x3, x4, а так же ход и скорость перемещения натяжного груза 8,8. В работах [2, 4] была получена модель в виде системы дифференциальных уравнений :

(2 т_г + 2 m_n + m_np ) Х 1 + т_гХ 2 + m_nХ 4 + 2Cx 1 - Cx ₂ - Cx ₄ + 0,5G_n lw + 0,5 G_r lw sgn x 1 +

+ 2t] x — x x2 — x xA = —пр- sgn x - x , 1        2        4     R б          с1

тгxv +4тгx2 + тгx3 -Cxx + 2Cx2 — Cx3 GGJwsgnx2 — xxx +2?7x2 -г/x3 = 0, m2x2 + 2тг + 2m, x3 + m„x4 - Cx2 + (2C + 0,25CK)x3 - (C + 0,25CK)x4 -

— 0,5CK8 — r/x2 Vlx x3 — xx^ + (0,5G2 lw + 0,5Gn lw) sgn x3 = 0, m„x + m„x3 + 4m„x. - Cx. - (2C + 0,25C)x3 + (C + 0,25C )x. + 0,5C8 -п1       п 3         п 4        1                       к 3                     к 4

-rj 5c_x— xx₃ + 2 т]Х_л + GJ w sgn , x ₄ = 0,

G

^^ 8- 0,5C к x з + 0,5C к x 4 + C к5 + GHy + GHy f sgn 5 = 0, g где Мпр - момент двигателя, приведенный к валу приводного барабана, Re - радиус приводного барабана, l - длина ленты, GHy - вес груза натяжного устройства, m2 - масса участка грузовой ветви, m„ - масса участка порожней ветви, mnp - масса привода, w -коэффициент сопротивления движению, f - коэффициент сопротивления движению натяжных грузов, т] -вязкость ленты с грузом, C - жесткость ленты, C - жесткость канатов натяжного устройства, G , G - погонный вес движущихся частей соответственно груженной и порожней ветви.

В матричной записи модель имеет вид:

MX N XX + CX + S sgn X + G = P sgn( Xc - X1)Mnp, где матрица M5x5 - матрица масс, N5x5 - матрица коэффициентов вязкого трения, C5х5 - матрица коэффициентов жесткого трения, S5х5 - матрица коэффициентов сил сопротивления движению. Матрицы P5x1 и G5x1 связаны с внешними воздействиями на ленту со стороны движущего привода и натяжного устройства.

В этой модели, используя каноническое правило О. Коши, выберем в качестве координат состояния перемещения и скорости сосредоточенных масс X = ( x 1 , x ₂,..., x ₁₀) ^T ,:

x ₁     x ₁,

^x1 = ^x 6 ,

x ₂    x ₂ ,

x x ,

x ₃    x ₃,

x ₃    x ₈ ,

x ₄    x ₄ ,

x ₄    x ₉ ,

x ₅      ,

x5 = 5.

Модель движения ленты конвейера в пространстве состояний представлена в виде системы нелинейных дифференциальных уравнений: ~~       ~   ~

X = ^AX*B_xsgn x c- x ₆ М пр + B 2 sgn X+BG Hy ,                    (2)

Здесь первое управляющее воздействие обозначим U = Мпр - это движущий    момент, создаваемый приводом, второе управляющее воздействие U2 = sgn X - силы сопротивления движению сосредоточенных масс ленты, третье управляющее воздействие U3 = G - вес натяжного устройства

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

X = AX +~ sgn Х_с -x₆ U + B ~₂ U₂ + B ~₃ U ₃,

~

Матрица состояния системы A представляет собой блочную матрицу,

включающую в себя матрицы М ¹N и М ^ХС , матрицы управления так

~

~~ блочные: матрица Bx включает в себя М 1P, B2 включает в себя

М⁴ 5 ,

же

B ^~ 3

включает в себя М G :

B ^~ 3

10 1

~

A 10 10

0 5 5

M ¹ C 55

E 5

M ¹ N 55

,

B ^~

0 5 1

M ¹ P 5 1

,

B ^~

0 5 5

05 5

0 5 5

_M 1 _S

05 1

M ¹ G 5 1

•

Компьютерное входящей в пакет

моделирование производилось в системе SIMULINK,

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

5с = Ax+Bu, y Cx Du.

Для реализации управления три внешних воздействия объединены в один вектор U , размерность которого 12 1. Матрица управления B= Bx: B2: B3 , размерность матрицы B 10x12. Матрица D - нулевая. В качестве выходных сигналов рассмотрим скорости движения сосредоточенных масс. Тогда матрица C diag 0 0 0 0 0 1 1 1 1 1

Схема моделирования в ППП Matlab приведена на рис. 2.

Рис. 2. Схема моделирования движения ленты конвейера

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

• •отключить момент нагрузки, поступающий на привод со стороны ленты;

• •обнулить сигнал задания по скорости, подаваемый в систему управления приводом.

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

Рис. 3. Разгон и торможение ленты

Отключение движущего привода происходит на 60 секунде. Из рисунка видно, что аппроксимирующие массы движутся плавно и совместно. Полная остановка конвейера происходит за 10 секунд.

Рассмотрим теперь торможение конвейера с использованием тормозных колодок, причем в трех вариантах: использование тормоза на головном барабане, использование тормоза на хвостовом барабане, совместное включение обоих тормозов.

Для того чтобы моделировать эти режимы, следует преобразовать модель, так как в исходной модели не предусмотрены тормоза. Применение тормоза - это приложение тормозящего момента. Если тормозные колодки располагаются на головном барабане, то изменяется первое уравнение системы (1), в него добавляется еще одно управляющее воздействие:

(2 m. + 2 m„ + m „) X, + т_гx^ + m„x ₄ + 2 Cx j - Cx₂ - Cx₄ + 0,5G_n lw + 0,5G_Z lw sgn X, +

+ 21] .X] - X]^- ip. =-- пр Sgn X -X, - ТОР Sgn(0- X j) , 1       2       4      _Rб          с     1       _Rб               1

где M - тормозной момент.

В матричной модели вводим матрицу P 2 0 0 0 0 T , следовательно, блочная матрица  B~ изменит размерность и примет

~ вид: B

1 10 2

05 2

M ¹ P 5 1 M ¹ P 15 1

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

т_гX₂ + 2 ш_г + 2 m_n X₃ + m_nx₄ - Cx ₂ + (2 C + 0,25 С_к ) x ₃ - ( C + 0,25 С_к ) x₄ -

M TO^ sgn(0- X 3 ), R _б

P 0  0  2  0  0 ^T ,

- 0,5С_к5 -i]x2 + 2i] X₃ -1] X₄ + (0,5О_г lw + 0,5G_n lw ) sgn X₃ = -

В матричной модели вводим матрицу следовательно, блочная матрица B~ изменится аналогичным образом.

Если же тормоза располагаются на обоих барабанах конвейера, то в модель нужно ввести матрицу P

T

а матрица B ^~ будет

иметь размерность 10 3 .

Эти изменения размерности следует учитывать при моделировании в ППП Matlab: в блоке State-Space изменяется размерность матрицы управления B = Bx: B2: B3 10x14. Применение тормозов можно моделировать, используя блок Look-Up-Table, формируя в нем постоянную величину, при условии, что скорость сосредоточенной массы отличается от 0. Величины тормозных моментов можно изменять.

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

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

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

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

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

В этом случае движение сосредоточенных масс имеет колебательный характер. Следовательно, при таком режиме остановки конвейера следует применять управляемое натяжное устройство, которое обеспечивает поддержание тяговой способности привода, иначе лента будет проскальзывать на приводном барабане из-за уменьшения натяжения в сбегающей ветви. Разработанная в [2] модель конвейерной установки легко позволяет это сделать. Кроме того, модель дает возможность управлять тормозными режимами, для чего нужно усложнить моделирующие блоки, реализуя в них технологические рекомендации [1]. Например, применять сначала свободный выбег, а после некоторого снижения скорости движения накладывать механические тормоза.

автоматизация, моделирование, Simulink, ленточный конвейер, алгоритм управления, торможение, тяговый фактор, проскальзывание ленты automation, modeling, Simulink, belt conveyor, control algorithm, braking, traction factor, belt sliding