Режим работы привода ленточного конвейера при случайном грузопотоке Текст научной статьи по специальности «Математика»

© В.П. Дьяченко, 2013

УДК 621.867.2 В.П. Дьяченко

РЕЖИМ РАБОТЫ ПРИВОДА ЛЕНТОЧНОГО КОНВЕЙЕРА ПРИ СЛУЧАЙНОМ ГРУЗОПОТОКЕ

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

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

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

Возможность учесть реальную динамику величины грузопотока дает использование для его описания случайных процессов типа «кенгуру» (КП) [3], которые являются обобщением случайного телеграфного процесса. В работе [2] дана подробная характеристика этих процессов применительно к описанию забойных грузопотоков. Приведем здесь лишь некоторые основные их свойства, отличные от свойств гауссовского экспоненциально коррелированного процесса, которым обычно описывают случайные грузопотоки ленточных конвейеров горных предприятий:

— интенсивность распределения времени пребывания величины грузопотока а на различных уровнях а; функционально связана с величиной этих уровней V =^а), а не является постоянной величиной;

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

— грузопотоки имеют дополнительную характеристику я(а)=-у(а)р(а)Аф — плотность вероятности распределения уровней (состояний) грузопотока (где vсp — средняя интенсивность переходов с одного уровня на другой);

— экспериментальному определению подлежат любые два из трех указанных выше распределений: я(а), v(а)/ и р(а).

При математических выкладках удобнее пользоваться вместо v(a) обратной величиной T(a)=1/v(a) — средним временем пребывания процесса (величины грузопотока) на уровне a.

Рассмотрим сначала случай, когда параметр v(a) не зависит от уровня грузопотока a, т.е. когда v = const . Корреляционная функция загрузки конвейера Q для этого случая определялась в работе [1], но, поскольку целью было вычисление только дисперсии, т.е. величины KQ (т) при т = 0 , то корреляционная функция определена в этой работе только для интервала 0 < т < T k , где Tk — время движения груза по конвейеру.

В действительности, как показано ниже, форма корреляционной функции KQ (т) при |т| > Tk имеет другой вид, нежели при |т| < Tk .

Загрузка конвейера является стационарным случайным процессом и описывается следующим выражением:

Tk

Q = Ja(t )dt

о

Примем корреляционную функцию грузопотока в виде:

Ka(T) = Da- е-*,

где Da — дисперсия величины грузопотока,

2 2

от ^ от ^

Da = J (a - Ma) - exp (- - v |т|) p(a)da = J (a - Ma ) p(a)da. о о

При т=0

2

от ^

Da = J (a- Ma) p(a)da .

о

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

от2

Ma = J a p(a)da .

о

Согласно работе [3], корреляционная функция загрузки Q, как интеграла от a(t) , может быть определена по формуле:

Kq (т) = 2 J (т - 0) - Ka (0 - т) d0 = 2TJ (т - 0) - Da - exp (-v -10 - т|) d0 .

о о

Если Tk > т , т.е. при о < т < Tk знаком абсолютной величины в показателе экспоненты можно пренебречь. Тогда получим вычисленное в работе [1] значение корреляционной функции KQ (т) и дисперсии Q:

KQ (т)= Dai

v2

e-v(Tk-т) + e-v(Tk+т) - 2е-

= Кд (°) =

2Т 2 2е

vTk Л

V а

.2 + " ,,2

При х = Тк получаем:

Ко (Тк) = ^К1 + е-** " 2^Тк ] = " е ^Тк) .

Если х > Тк , получаем другое выражение для корреляционной функции:

Ко (х) = Оа--1^[е^(х-Тк)+ е^ -2е-

При х = Тк это выражение совпадает с полученным выше. Совпадает и производные по х обоих выражений, равные при х = Тк :

±Ко (х) = -Оа(1 -е-- ) <0.

-v(х+Tk)

вх

Вторые производные также совпадают и равны при х = Тк :

Ко (х) = Оа(1 - е-* )2 > 0.

Следовательно, при х = Тк кривая Ко (х) является вогнутой. При х = 0 первая производная Ко (х) равна нулю.

Таким образом, кривая Ко (х) имеет вид, показанный на рис.1. Точка перегиба А определяется координатой хд, которую можно определить из уравнения:

е" = V2е^к -1.

Кд(т)

0

тА т='Гк

Рис. 1. Корреляционная функция загрузки конвейера

т

Зд(ю)

_а_

К

г„

г„

Рис. 2. Спектральная плотность загрузки конвейера

Среднее время корреляции процесса О (t) получается интегрированием нормированной корреляционной функции на отрезках 0 < т < Тк и Тк < т < <»

Т,.

и составляет т = —

К.ср 2

Обобщенная формула для корреляционной функции Кя (т) имеет вид:

Кд (т) = Оа ^(е-*-Тк1 + е-^Тк - 2^) + Оа 2I(Тк - т),

где I (Тк - т) — единичная функция Хевисайда, равная нулю при Тк < т и единице при Тк > т.

Спектральная плотность загрузки конвейера 0(1) определяется как произведение

^ (®) = М05)- |Ф (^®)|^

где Ф (/ш) — преобразование Фурье от переходной функции блока интегрирования по времени; / = \/-Т — мнимая единица.

Функция (ш) может быть определена из выражения _2

м®)=- "

'а\~ / 2,2

К ш + V

Передаточная функция этого блока (преобразование Лапласа) равна W (р) = р (Т - е-ТкР ).

Подставляя в это выражение р = /ш , получаем:

2

0

Ф (j-) = —(1 - e-J-Tk ),

j- v '

\Ф (j-)f = -1(1 - e-aTk )-l(1 - e-aTk )=-,(ch-Tk )=-, (1 - cos ®Tk ).

1 1 j- v ' -j- v ' - -

Отсюда

SQ(-) = §ста —(—-г •(1 - cos aTk).

л - (- + V )

Общий вид этой функции представлен на рис.2. Она является периодиче-

„ , 2л

ской функцией частоты ю с периодом —.

Tk

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

аналогии с предыдущими выкладками получаем:

2

1 (e( + e-vTk - 2e-VH) + -v2 v ' v

w "1 О

KQ (x) = J(a-MJ •[ -I (e-^^l + e"Tk - 2e-vlxl) + 2 I (Tk-x) ] p(a)da,

2

2 V

^ (®) = I (а - Ма )----П ' (1 - ^ ® Тк ) Р(а)ва .

0 ТС ю (со + V )

Функция распределения величины загрузки конвейера, как сумма большого числа малых одинаково распределенных величин, согласно работе [5], близка к нормальному распределению с математическим ожиданием Мд=МаТк, прямо пропорциональным времени Тк, и с дисперсией, определенной выше.

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

1. Шахмейстер Ё.Г., Дмитриев В.Г. Вероятностные методы расчета транспортирующих машин. — М.: Машиностроение, 1983. — 256 с.

2. Дьяченко В.П. Моделирование работы конвейерных транспортных линий/ Горный информационно-аналитический бюллетень, 2006, №12. — М.: МГГУ, с. 195—196.

3. Brissaud A., Frisch U. Solving linear stochastic differential equations. — J. of Math. Phys., 1974, № 5, р. 524-534.

4. Свешников A.A. Прикладные методы теории случайных функций. — М.: Наука, 1968. — 463 с.

5. Байхельт Ф., Франкен П. Надежность и техническое обслуживание. Математический подход: Пер. с нем. — М.: Радио и связь, 1988. — 391 с. ir.ua

КОРОТКО ОБ АВТОРЕ -

Дьяченко Вячеслав Петрович — кандидат технических наук, профессор, kafgmt@msmu.ru, Московский государственный горный университет.