Файл: Багин, Б. П. Основы статистической динамики одноковшовых экскаваторов обзор.pdf

Спектральной плотности, описываемой выражением (57),. соответствует корреляционная функция вида:

ризуют частотные свойства нагрузок. При увеличении а кор­ реляционная функция убывает быстро, а колебания нагрузок становятся более резкими и беспорядочными. Параметр р указывает на наличие в спектре процесса Poia(t) круговых частот, близких к частоте р, но со случайными амплитудами и фазами, что и определяет наличие пика в спектральной плотности SPв (со) в зоне частоты со = 5,3 сек~К

Рис. 16. Спектральная плотность высокочастотной составляю­ щей Р01в для забоев «А» и «Б»

Используя полученные выражения корреляционных функ­ ций (55) и (58), определим корреляционные функции слу­

чайных центрированных процессов c(t) и K(t). Принимая во внимание формулу (51), получим

Если принять толщину стружки c(t) в период установив­ шегося копания в хорошо взорванном грунте (сыпучая сре­

да) примерно постоянной (сср), тогда из формулы (52)

40

Коэффициенты вариации толщины стружки Л'„.с и удельного сопротивления копанию Кв-к определяются из соотношений:

Значения удельного сопротивления грунта копанию Ки

средней толщины стружки сср, их среднеквадратические от­ клонения (ок, огс) и коэффициенты вариации для двух забоев «А» и «Б», определенные по приведенным выше зависимо­ стям, представлены в табл. 2.

гут быть использованы в качестве исходных для определения вероятностных характеристик силы сопротивления грунта копанию Pm (t) для других экскаваторов при работе в ука­ занных грунтах, но отличающихся своими динамическими характеристиками.

Определим дисперсию составляющей Лив(0> приведенной

При уср=47° имеем Df „= 4 5 tc2 для забоя «А» и й^в= |100гс2

— для забоя «Б».

41

р7,

примера соотношение АЧХ системы и спектральной плотно­

с дисперсией DFb. Таким образом, при проектировании экска­ ватора необходимо располагать сведениями о характере спектральной плотности внешнего возмущения Poi(t) и дина­ мическими характеристиками системы.

ВЕРОЯТНОСТНЫЕ ХАРАКТЕРИСТИКИ НАГРУЗКИ В МЕХАНИЗМЕ ПОДЪЕМА

Полученные выше вероятностные характеристики силы сопротивления грунта копанию позволяют перейти к опреде­ лению вероятностных характеристик (прогнозированию) на­ грузок в проектируемой машине.

Задача вероятностного определения (или прогнозирова­ ния) нагрузки заключается в определении статистических характеристик нагрузки в механизмах подъема и напора экскаватора.

Для проведения проверочных расчетов деталей механиз­ ма подъема на прочность и выносливость с учетом действи­ тельного случайного процесса нагружения X (t) требуются такие статистические характеристики нагрузки как матема­ тическое ожидание mx (t), среднеквадратическое отклонение нагрузки ах (или коэффициент вариации) и среднее число

выбросов в единицу времени на среднем уровне «*[6].

При копании во взорванном скальном забое с коэффи­ циентом разрыхления /(V)>1,4 (сыпучая среда) среднестати­ стические траектории копания представляют плавные поло­ гие кривые. При этом высота черпания, необходимая для наполнения ковша,'не превышает, как правило, 60% высоты

42

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

(63)

Случайные составляющие ЛиЛО и Ли» (0 являются центрированными, стационарными независимыми процесса­

ми, т. е. их математические ожидания равны нулю. Используя полученные выше результаты, можно записать выражения для корреляционных функций процессов Я,11м(/) и Я01а(£);

Спектральные плотности процессов Р0,„(£) и P„lBU), определенные как преобразование Фурье от корреляционных функций, определяются выражениями:

(66)

Механизм подъема в процессе установившегося копания можно рассматривать линеаризованной стационарной дина­ мической системой. Исчерпывающими динамическими харак­ теристиками линейной системы являются частотная характе­ ристика 'W(Jm) или временная характеристика — импульсная переходная функция g(r). Прогнозирование нагрузки сводит­ ся к определению статистических характеристик «выхода» X(t) линейной системы с динамической характеристикой W(j(o) или g(x) при действии на ее «вход» внешнего случай­ ного воздействия Ли (0 с заданными статистическими харак­ теристиками. Нагрузка в механизме подъема представляется статистической моделью

X (t) = тх {t) + (/х (,t) 4- Zx (t).

Низкочастотная Ux (t) и высокочастотная Zx (l) случайные составляющие являются практически независимыми стацио­ нарными процессами, следовательно, дисперсия процесса X(t) определится выражением Dx = Ни + Dz.

43

Среднее число выбросов в единицу времени на среднем уровне для стационарного процесса определится по формуле

где 3'x=VDx— среднеквадратическое отклонение нагрузки; ох— среднеквадратическое отклонение скорости

изменения нагрузки X(t).

т. е. для определения дисперсии скорости изменения случай­ ного процесса необходимо, чтобы корреляционная функция нагрузки Кх (т) имела вторую производную в нуле. Более

удобным является определение дисперсии а2, если известна

х

спектральная плотность процесса X(t)

Таким образом, для определения среднего числа выбросов в единицу времени необходимо знать корреляционные функ­ ции или спектральные плотности случайных составляющих нагрузки Ux (t) и Zx {t).

Связь между «входом» и «выходом» линейной системы определяется уравнением свертки [17]

Импульсная переходная функция g(r) является реакцией системы на единичный импульс (6-функцию). Между им­ пульсной переходной функцией g(r) и передаточной функ­ цией системы W(р) существует взаимно однозначное соот­ ветствие, определяемое преобразованием Лапласа [17].

Вследствие существенного демпфирования в системе вре­ мя переходного процесса тп< 1 сек. Учитывая, что в данном случае рассматривается процесс установившегося копания,

44

Таким образом, низкочастотные составляющие внешней нагрузки проходят в динамическую систему практически без изменения. К такому же выводу можно прийти при рассмот­ рении АЧХ системы в частотном диапазоне со —O-f-З сек~1. При этом выражение для дисперсии низкочастотной состав­ ляющей нагрузки можно записать в виде

Корреляционная функция Лс/(т) при т = 0 имеет вторую производную, поэтому можно сразу определить дисперсию скорости изменения низкочастотной случайной составляющей

Ux (t):

*ср

Корреляционная функция высокочастотной составляющей нагрузки Zx (t) связана с корреляционной функцией высоко­ частотной случайной составляющей внешнего воздействия /301в(0 соотношением [17]:

Kz (7) = - г т - J J g (о) g Ы КРв + {-п - --2)\ -d-n d-2. (76)

S in - у О О

Вычисление двойного интеграла является весьма трудо­

емким, поэтому для определения дисперсий Dz и Dz лучше пользоваться соотношениями между спектральными плот­ ностями:

Интегралы в этих выражениях могут быть вычислены

•обычными методами графического интегрирования, если из­ вестны графики кривых спектральной плотности S.рв(ы) и ам­ плитудно-частотной характеристики (/.- (м) = | У7 (/и.)|-).

Имеется также и аналитический метод вычисления Dz и Dz , основанный на предположении, что как спектральные плотности, так и частотная характеристика представляют со­ бой дробно-рациональные функции от со [17]. Этот метод позволяет получить в явном виде связи между величиной дисперсий Dz и Dj,, параметрами спектральной плотности

ав и р и динамическими параметрами системы. Для возмож­ ности использования этого метода необходимо выражение для спектральной плотности SPB(co) представить в виде абсо-

45