Методика оценки режимов нагружения в элементах привода торфяных фрезерующих агрегатов на стадии проектирования

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

МЕТОДИКА ОЦЕНКИ РЕЖИМОВ НАГРУЖЕНИЯ В ЭЛЕМЕНТАХ ПРИВОДА ТОРФЯНЫХ ФРЕЗЕРУЮЩИХ АГРЕГАТОВ НА СТАДИИ ПРОЕКТИРОВАНИЯ

К.В. Фомин, К.С. Крылов, В.Е. Харламов

Тверской государственный технический университет

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

Для оценки режимов нагружения машин в настоящее время [5] широко используют функцию плотности распределения W(x) или интегральную функцию

. (1)

фрезерующий торфяной нагружение прочностной

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

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

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

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

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

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

На схеме введены следующие обозначения Ii - моменты инерции сосредоточенных масс; Mc(, , t) - моменты упругости валов; Мд и Мп- соответственно, моменты двигателя и сопротивления на рабочем органе.

На основании расчетной схемы с использованием методов [6, 10] могут быть получены дифференциальные уравнения, описывающие динамические процессы в приводе:

, (2)

где i - углы закрутки сосредоточенных масс; Mci - упругие моменты в соединениях; Ii - моменты инерции сосредоточенных масс; Mд, Mi - момент двигателя и моменты сопротивления на рабочих органах.

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

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

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

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

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

, (3)

где M - число плоскостей резания; n - номер импульса нагружения на m-й плоскости резания; Мпт() - функция, описывающая изменения момента нагружения в приводе при взаимодействии одиночного режущего элементе в m-й плоскости резания; пт - момент возникновения n-го импульса нагрузки от режущего элемента в m-й плоскости резания; Рпт - случайные параметры n-го импульса от режущего элемента в m-й плоскости резания.

При взаимодействии с древесными включениями для нагрузки в приводе можно записать

, (4)

где S - число актов взаимодействия режущих элементов с n-м древесным включением; - функция, описывающая переходный процесс в элементах привода при взаимодействии одиночного режущего элемента при s-м акте взаимодействия с n-м древесным включением; - момент возникновения импульса нагрузки при s-м акте взаимодействия режущего элемента с n-м древесным включе-нием; Рns - случайные параметры импульса нагружения при s-м акте взаимодействия режущего элемента с n-м древесным включением.

Использование таких моделей формирования момента нагружения в элементах привода позволяет определить его характеристическую функцию и плотность распределения [12, 13].

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

. (5)

Значения математического ожидания mT и дисперсия DT нагрузок могут быть определены с помощью методики, предложенной в [14].

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

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

фn Tn, (6)

где фn - длительность n-го переходного процесса; Tn - интервал времени между началами n-го и n + 1-го переходного процессов.

Величины ф и T являются случайными. Если известно их совместное распределение W(ф,T), то вероятность выполнения условия (6) определится выражением

. (7)

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

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

Характеристическая функция нагрузок в этом случае определится следующим образом. Пусть заданы число I различных форм древесных включений и вероятность встретить их в залежи Pi, а также плотности распределения размеров древесных включений и их глубины залегания в залежи в пределах глубины фрезерования Wd(y), WH(y). Разобьем распределения на достаточно малые интервалы dm…dm+1, Hl…Hl+1 таким образом, чтобы изменения величин d и H в пределах выделенных диапазонов мало изменяло характер и величину переходного процесса в приводе. Каждому сочетанию размеров и глубин залегания древесных включений dmcp=(dm + dm+1)/2 и Hlcp = (Hl+Hl+1)/2, для всех их форм и условий взаимодействия соответствуют значения нагрузок Mi (dmcp, Hlcp, t), которые могут быть определены численно или аналитически. Выделяя в последовательности импульсов нагружения импульсы одинаковой формы, одинаковой амплитуды и длительности (то есть соответствующие одинаковым условиям взаимодействия рабочего органа с древесными включениями одной формы, глубины залегания и размеров) и учитывая, что интервалы между импульсами в данном случае будут распределены также по экспоненциальному закону, запишем выражение для характеристической функции такого процесса [12]

, (8)

где iml - число переходных процессов в единицу времени, возникающих в приводе при взаимодействии с древесными включениями размерами dmcp, глубиной залегания Hlcp и i-й формы.

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

. (9)

Подставляя (8) в (9)

,(10)

где - число древесных включений, попадающих на рабочий орган в единицу времени.

Значения Pm и Pl равны:

, .

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

При уменьшении интервалов разбиения Wd(y) и WH(y), для характеристической функции момента получим

(11)

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

.

Плотность распределения находится с помощью обратного преобразования Фурье характеристической функции [13]

.

В некоторых случаях, (если функция, описывающая изменение нагрузки на рабочем органе имеет сложный вид) для определения плотности распределения нагрузок целесообразно использовать разложение выражения (10) в ряд вида [12, 13]

, (12)

где - кумулянты k-го порядка, определяемые с помощью выражений

. (13)

Обращение характеристической функции в виде (12) по Фурье [12] позволяет найти плотность вероятности

.

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

ч1 = mM, ч2 = у2M = DM, ,

.

Если начиная с k = 3 разложить экспоненту под интегралом по степени нx, то почленное интегрирование приводит к представлению W0(x) в виде ряда по производным Шn(y), называемым рядом Эджворта [12]

.

Для анализа более удобна другая форма записи ряда Эджворта [13]

,

где и - соответственно, 3 и 4-й полиномы Эрмита.

Характеристическая функция и плотность вероятности суммарного момента нагружения в приводе с учетом фрезерования торфяной залежи и взаимодействия с древесными включениями определится с использованием свойства аддитивности кумулянтов для независимых случайных величин [12, 13]

,

где - кумулянт k-го порядка момента нагружения в приводе при взаимодействии с беспнистой залежью; - кумулянт k-го порядка момента в приводе при взаимодействии с древесными включениями

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

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

Полученная информация о режимах нагружения служит исходным материалом для прочностного расчета элементов привода [19] и выбора его оптимальных параметров.

Библиографический список

1. Солопов С.Г., Горцаколян Л. О., Самсонов Л. Н. Торфяные машины и комплексы. - М.: Недра, 1981. 415 с.

2. Самсонов Л.Н. Фрезерование торфяной залежи. - М.: Недра, 1985.

3. Дунаев Ю.М., Жаворонкова В.А., Каменский Ю.А. Надежность трансмиссии торфодобывающих фрезеров //Технология и комплексная механизация торфяного производства. - Калинин: КПИ, 1982. С. 43-46.

4. Лукьянчиков А.Н., Харламов В. Е. Показатели надежности фрезера МТФ-14. //Технология и комплексная механизация торфяного производства. - Калинин: КПИ, 1982. С. 46-49.

5. Решетов Д.Н., Иванов А.С., Фадеев В.З. Надежность машин. - М.: Высшая школа, 1988. 234 с.

6. Светлицкий В.А. Статистическая механика и теория надежности. - М.: МГТУ им. Баумана, 2002. 504 с.

7. Николаенко Н.А., Ульянов С.В. Статистическая динамика машиностроительных конструкций. - М.: Машиностроение, 1977. 367 с.

8. Болотин В.В. Случайные колебания упругих систем. - М.: Наука, 1979. 335 с.

9. Гусев А.С., Светлицкий В.А. Расчет конструкций при случайных воздействиях. - М.: Машиностроение, 1984. 240 с.

10. Докукин А.В., Красников Ю.Д., Хургин З.Я. Статистическая динамика горных машин. - М.: Машиностроение, 1978. 238 с.

11. Моделирование нагрузок на рабочем органе торфяного фрезерующего агрегата. //Перспективы развития Волжского региона. Материалы Всероссийской заочной конференции. - Тверь: ТГТУ, 2002. Вып. 4. С.124-127.

12. Рытов С.М. Введение в статистическую радиофизику. -М.: Наука, 1976. 496 с.

13. Тихонов В.И. Статистическая радиотехника. - М.: Радио и связь, 1982. 624 с.

14. Методика анализа динамических нагрузок в элементах привода торфяного фрезерующего агрегата. //Вестник Тверского государственного технического университета. - Тверь: ТГТУ, 2002. № 1. С.10-14.

15. Малков Л.М., Пономарчук Ф.С. Размеры древесных включений в обработанном слое отремонтированных площадей. //Технология и комплексная механизация торфяного производства. - Калинин, 1986. С. 14-17.

16. Зиновьев Д.А. Исследование пнистости верховой залежи на участке производства фрезерного торфа. //Торф и его переработка. - Л., 1978. С.26-29.

17. Анализ нагрузок на рабочем органе торфяного фрезерующего агрегата при взаимодействии с древесными включениями. //Торфяная отрасль и повышение эффективности использования энергобиоресурсов: Мат. научно-практической конф. - Тверь: ТГТУ, 2000. С.130-132.

18. Параметры импульсов нагружения на фрезе при взаимодействии с древесными включениями. //Технология и комплексная механизация торфяного производства. -Тверь: ТГТУ, 1997. С.146-150.

19. Серенсен С.В., Когаев В.П. Руководство по расчету на усталость деталей машин. - М.: Машиностроение, 1972. 98 с.

20. Когаев В.П. Расчеты на прочность при напряжениях, переменных во времени. - М.: Машиностроение, 1977. 231 с.

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