CC BY
103
Научный журнал КубГАУ, №101(07), 2014 года
1
УДК 631.372:001.4 UDC 631.372:001.4
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ MATHEMATICAL MODELING OF MTU
ТЯГОВОЙ НАГРУЗКИ МТА TENSILE LOAD
Кравченко Владимир Алексеевич Kravchenko Vladimir Alekseevich
д.т.н., профессор Dr.Sci.Tech., professor
Азово-Черноморский инженерный институт Azov-Black Sea engineering institute Donskoy State
ФГБОУ ВПО «Донской государственный аграрный Agrarian University, Zernograd, Russia университет» г. Зерноград, Россия
Дурягина Вероника Владимировна Duryagina Veronika Vladimirovna
Г амолина Ирина Эдуардовна Gamolina Irina Eduardovna
к.т.н., доцент Cand.Tech.Sci.
Южный федеральный университет, Таганрог, Southern Federal University, Taganrog, Russia
Россия
В статье дан обзор результатов математического The article is devoted to mathematical modeling of моделирования случайного процесса тяговой random process that occurred in tensile load of serial
нагрузки серийного и опытного МТА. Обсуждается and experienced MTU. The efficiency of EDM эффективность применения УДМ в трансмиссии application in transmission is discussed
трактора
Ключевые слова: МАТЕМАТИЧЕСКОЕ Keywords: MATHEMATICAL MODELING,
МОДЕЛИРОВАНИЕ, МАШИННО-ТРАКТОРНЫЙ MACHINE-TRACTOR UNIT, TRACTIVE АГРЕГАТ, ТЯГОВОЕ СОПРОТИВЛЕНИЕ, RESISTANCE, ELASTICALLY DAMPING
УПРУГОДЕМПФИРУЮЩИЙ МЕХАНИЗМ, MECHANISM, APPROXIMATION, CORRELATION
АППРОКСИМАЦИЯ, КОРРЕЛЯЦИОННАЯ FUNCTION, RANDOM PROCESS
ФУНКЦИЯ, СЛУЧАЙНЫЙ ПРОЦЕСС
При выполнении машинно-тракторным агрегатом (МТА)
энергоемких сельскохозяйственных работ на него действует постоянно изменяющаяся нагрузка, колебания которой достаточно велики и достигают 30-40 % от крюкового усилия трактора.
Возникающие при эксплуатации МТА случайные колебания являются следствием усиления воздействия на МТА неровностей поля и неравномерности взаимодействия рабочих орудий с обрабатываемым материалом: интенсифицируются вертикальные колебания, усиливающие колебания нагрузки на крюке. Эти колебания воздействуют на работу двигателя, приводят к колебаниям остова, поскольку все ходовые системы современных тракторов имеют упругую подвеску и вызывают в деталях и
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
2
узлах переменные по величине напряжения [1, 3, 4]. Колебания тяговой нагрузки и силы сопротивления движению вызывают буксование движителей. Буксование влияет и на эксплуатационно-технологические показатели работы МТА [5].
Одним из способов решения проблемы является установка упругих виброгасителей в элементы и звенья машинно-тракторного агрегата.
Различные научные исследования [1, 2] показывают, что введение упругой связи в отдельные механизмы трактора позволяет снизить динамичность протекающих процессов. Эффективность применения упругодемпфирующих механизмов (УДМ) в трансмиссии объясняется тем, что содержащиеся в них упругие деформирующиеся элементы поглощают и рассеивают энергию колебаний тяговой нагрузки, снижая динамическую нагруженность привода, буксование движителей и т.д.
Представляется актуальным исследование влияния тягового сопротивления на работу МТА с УДМ в трансмиссии. Используемый УДМ [2] в трансмиссии МТА на базе МТЗ-80 является «частотным фильтром». Он ослабляет колебания, действующие на двигатель от внешних воздействий в определенном диапазоне частот, снижая наиболее выраженные гармоники внешних колебаний нагрузки на различных режимах работы трактора. Основные параметры УДМ можно варьировать, добиваясь снижения наиболее выраженных гармоник внешних колебаний нагрузки.
Для количественной оценки сопротивления почвы, которое преодолевается рабочими органами МТА при деформации и разрушении почвенного пласта, используется оценка выходного энергетического показателя работы - оценка тягового сопротивления МТА [3, 4].
Тяговое сопротивление (нагрузка) сельскохозяйственных машин зависит от множества факторов, которые условно можно разделить на три
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
3
группы: конструктивные, почвенно-климатические и эксплуатационные [1, 3, 4].
Интенсивность колебаний нагрузки зависит от [3, 4]: типа машины-орудия (плуга, сеялки и др.) и ее характеристик (количества корпусов, ширины захвата и др.); технологического режима (глубины обработки, угла развала и др.); количества машин-орудий в агрегате; скорости движения агрегата; типа фона и его состояния; зональных условий (удельного сопротивления, глыбистости почвы, влажности и т.д.).
Движение рабочего органа МТА в необработанной почве сопровождается явлениями сжатия и разрушения пласта в различных направлениях, которые обуславливают беспорядочный характер изменения тягового сопротивления [5]. Тяговое сопротивление, вследствие
неоднородности плотности почвы, имеет случайный характер.
Ряд авторов рассматривает тяговую нагрузку МТА как гармоническую функцию с варьирующейся в некоторых пределах частотой. Другие авторы, исследуя тяговое усилие, предлагают использовать некую функциональную зависимость, представляемую в виде пяти основных составляющих колебаний, различающихся между собой амплитудой и частотой [4].
В [1, 11] при исследовании тяговой нагрузки МТА с УДМ в трансмиссии трактора в качестве функции, характеризующей колебания
тягового сопротивления применяется с дискретным набором частот / из диапазона, соответствующего определенной сельскохозяйственной
операции.
Рассмотрим двумерный случай.
Уравнения движения агрегата и внешней нагрузки, колебания остова трактора в продольно-вертикальной плоскости описываются системой уравнений [1, 4]:
, (1)
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
4
ша^с - cr(zi ~ (7l)+ _ 4i.) + crC^i ~ (?з) + ^гС^з ~ Qi), (2)
1сФс = [Cr(A _ Qi) +■ “ <7i)]a — [cr(^i “ Qi) + Or(^= “ Qs)]^ +
[cr(Zj - q7) +
(3)
U +■ Тдз +
гдета - масса агрегата;
%с - отклонение по вертикали центра масс трактора от равновесного положения;
^i, 2= - отклонения по вертикали от положения равновесия центров переднего и заднего моста соответственно;
<7i, Ч* - реализация возмущающего воздействия со стороны микропрофиля на передний и задний мосты трактора;
а , ^ - расстояния от центра тяжести до переднего и заднего
мостов;
?с - тяговое сопротивление на рабочих органах;
- расстояние от линии действия силы тяги до центра масс
трактора;
Ьсг - расстояние от центра колеса до центра масс трактора;
(Р° - угол поворота остова трактора относительно поперечной горизонтальной оси, проходящей через центр масс трактора;
1с - момент инерции трактора относительно центра масс;
ГШ-усилие в контакте отпечатка шины;
^ -продольная деформация шины.
Силу тягового сопротивленияРс(Ппредставим в виде [1,4]:
где - равновесное значение силы сопротивления;
- экспериментальный коэффициент нарастания тягового сопротивления;
http://ej.kubagro.ru/2014/07/pdf/24.pdf
[>j
Научный журнал КубГАУ, №101(07), 2014 года
5
- скорость МТА;
- функция, характеризующая колебания тягового сопротивления и носящая случайный характер.
Основными этапами математического моделирования тяговой нагрузки являются: получение корреляционной функции
экспериментальным путем; нормирование экспериментальной корреляционной функции; аппроксимация корреляционной функции; получение случайного процесса тяговой нагрузки по известной корреляционной функции.
Для описания тягового сопротивления воспользуемся статистическим подходом, учитывающим сложность реальных зависимостей, влияние внешних факторов, приближенность эмпирических измерений. В связи с этим колебания тягового сопротивления будем описывать вероятностными методами, пользуясь результатами натурных экспериментов. Для описания случайных процессов используются следующие статистические функции: среднее значение квадрата
случайного процесса; плотность распределения; автокорреляционная функция; спектральная плотность.
Среднее значение квадрата дает элементарное представление об интенсивности процесса. Плотность распределения характеризует распределение вероятностей процесса в фиксированных точках. Автокорреляционная функция (АКФ) и спектральная плотность дают аналогичную информацию о процессе во временной и частотной областях соответственно [7]. Пусть колебания тяговой нагрузки происходят случайным образом только в вертикальной продольной плоскости. подчинена нормальному закону распределения и является случайной функцией, удовлетворяющей условиям стационарности. А.Я. Хинчин доказал теорему о том, что стационарные случайные процессы являются эргодическими [7].
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
6
Допущение об эргодичности стационарного случайного процесса позволяет предсказать все его статистические характеристики по одной реализации из ансамбля процесса с помощью усреднения.
Графики экспериментальной зависимости автокорреляционной функции для тяговой нагрузки с УДМ и без УДМ приведены на рисунке 1.
1 - серийный агрегат; 2 - агрегат с УДМ в трансмиссии
Рисунок 1. Автокорреляционная функция тяговой нагрузки (экспериментальные зависимости)
При проведении натурного эксперимента получаемая автокорреляционная функция имеет погрешности разного вида, поэтому при моделировании исследуемой функции целесообразно применять не интерполяционные, а аппроксимирующие подходы, основанные на интегральной близости. Имея результаты натурных испытаний в виде автокорреляционной функции, стоит задача аппроксимации с целью моделирования случайного процесса тяговой нагрузки. Существуют различные подходы к аппроксимации, в частности [7]: функциями заданного вида, разложением корреляционной функции в ряд по ортогональным функциям или по дисперсиям производных. При этом в качестве критерия приближения чаще всего используют минимум квадратичной погрешности аппроксимации, минимум экспоненциальновзвешенной аппроксимации и критерий равномерного приближения. Наиболее приемлемыми доступным является метод наименьших квадратов с оценкой аналитического выражения для автокорреляционной функции и
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
7
ее последующая аппроксимация с минимумом квадратичной погрешности.
Как показывают исследования [3, 7], для характеристики случайных процессов обычно применяют затухающие экспоненциальные и экспоненциально-косинусные автокорреляционные функции.
Для аппроксимации нормированной АКФ использовалось аналитическое выражение вида:
НЮ = -ъЛ^'0^, (4)
где^1, а7, Р - коэффициенты, оцененные с помощью метода
наименьших квадратов в пакете прикладных программ MathCad.
Наилучшие результаты (по «невязке») для операции «пахота» приведены ниже.
Таблица 1- Коэффициенты аппроксимации нормированной АКФ тяговой нагрузки серийного и опытного МТА
Наличие УДМ Л, А* сц «S
- 0,208 0,792 1,391 1,211 6,583
+ 0,725 0,275 4,597 0,291 0,319
Следующий этап исследования заключается в моделировании случайного процесса тяговой нагрузки МТА по заданной корреляционной функции.
Наиболее часто применяемым способом для получения дискретных значений случайного процесса в моменты времени, связанные с шагом дискретизации, является построение характеристик формирующего фильтра, описываемых уравнением Виннера-Хинчина [7]. Данное уравнение связывает спектральные плотности входного и выходного процессов линейной динамической системы с помощью частотной передаточной функции. Для построения случайного процесса
используются следующие подходы:
решение системы однородных дифференциальных уравнений [8] формирующего фильтра;
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
8
использование интеграла свертки, заменяемого суммой при малом шаге дискретизации [6];
применение стохастического разностного уравнения, полученного на основе использования теории непрерывных дробей [9].
Недостатками первых двух методов являются отсутствие взаимно однозначного соответствия между непрерывной и дискретной моделями, теоретически необоснованный выбор шага дискретизации [9].
Наиболее целесообразным, достаточно эффективным и простым в реализации является алгоритм моделирования стационарного случайного процесса на основе непрерывных дробей.
На рисунке 2 приведена схема, поясняющая суть метода:
КФ - корреляционная функция; НПФ - непрерывная передаточная функция; ИМ - идентификационная матрица В. Висковатова;
ДПФ - дробно-передаточная функция; СП - случайный процесс
Рисунок 2. Алгоритм моделирования на основе непрерывных дробей
В рассматриваемом случае нормированной корреляционной функции вида(4) соответствует НПФ формирующего объекта:
+ uLs +ип)
Ч-1?-, s Ч-
+ VA
■* 1
(5)
где s - ico А щ (t — Од ) ^j (J — 0,1,2 ) _ действительные функции, зависящие от Аг ^А2 а^ а2 р .
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
9
Непрерывная передаточная функция имеет два нуля:
sL3=7i±0i^ и три полюса:= 7з t ,san =-1?!, где Yi, Уз, &l, Qi -действительные числа. С помощью обратного преобразования Лапласа перейдем к весовой функции формирующего объекта
ь(У) = jiLe-v^ + ti7ey*tcos97t +цэ еУ** sin971 (6)
где th, - действительные коэффициенты.
Шаг дискретизации At выбираем из условия SP-идентифицируемости [9] для получения взаимно однозначного соответствия между непрерывной и дискретной моделями:
Мтахв2\сп . (7)
Так как
то идентифицирующая матрица (модифицированный
алгоритм В. Висковатова) имеет вид:
1 0 0 0 0 0
Ь(0) h(2At) h(nAf)
сгДС) c^fnAf)
■ ■■
«шОО am(nAO
...
(8)
где
(9)
ат_Мп +
«т-:(Ю
«т-1®
а.л (пдГ] = <5(лдГД сеп(?Ш) = hoшд т = 1,2,3,... ? п = ОД,2,3,... Ограничимся вычислениями элементов матрицы при появлении первой «нулевой» строки.
На основе идентифицирующей матрицы сформируем ДПФ:
G(Z ) =
h(Q)
а (0)-z1 1 + ^2 (0)-Z_1
(10)
1+
1+...
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
10
После преобразования для нашего случая ДПФ примет вид:
(11)
Дискретная модель случайного процесса в форме стохастического разностного уравнения имеет вид:
«Белый шум» применяется в тех случаях, когда интервал корреляции шума много меньше всех существенных постоянных времени системы, на которую воздействует шум, и спектральную плотность внешних воздействий можно приближенно считать постоянной. Использование белого шума позволяет упростить математические выкладки. В работе [10] был применен «белый шум» для моделирования тяговой нагрузки МТА.
Имея смоделированный случайный процесс, построим модельную корреляционную функцию ^ууС^)[11]для контроля полученных результатов.
На рисунках 3 и 4 представлены результаты математического моделирования: реализации случайных процессов (а)и соответствующие им автокорреляционные функции (б) тяговой нагрузки серийного и опытного МТА.
(12)
Здесь - искомый случайный процесс, _ «белый шум», с
корреляционной функцией
математические ожидания „ , дисперсии D
через коэффициенты модели ai и &i.
ласно [9] ^у связаны
х и
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
11
<1-------------Г--------------1-------------1-------------1-------------Р-------------1-------------Г
J_______________I_____________j I I___________________I______________L
О 2i -1 6 ft № 12 |Л
:К
а) Реализация случайного процессатяговой нагрузки серийного МТА
1 - Аппроксимация; 2 - Результат численного моделирования
б) - Нормированная АКФ тяговой нагрузки серийного МТА
Рисунок 3 -Результаты математического моделирования тяговой нагрузки
серийного МТА
т-----------1----------1-----------1----------1-----------1----------г
а) Реализация случайного процесса тяговой нагрузки опытного МТА
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
12
б) - Нормированная АКФ тяговой нагрузки опытного МТА.
Рисунок 4 -Результаты математического моделирования тяговой нагрузки
опытного МТА
Анализ результатов математического моделирования показал, что используемый метод дает хорошую согласованность модельной и аппроксимированной АКФ. Таким образом, полученные дискретные модели случайного процесса тяговой нагрузки серийного и опытного МТА применимы для исследования движения агрегата, влияния воздействия внешней нагрузки, эффективности применения УДМ в трансмиссии трактора на эксплуатационные показатели.
Полученные данные могут быть использованы при построении модели функционирования МТА.
Список литературы
1. Кравченко, В.А. Повышение динамических и эксплуатационных показателей сельскохозяйственных машинно-тракторных агрегатов: монография / В.А. Кравченко. - Зерноград: ФГОУ ВПО АЧГАА, 2010. - 224 с.
2. Кравченко В.А., Сенкевич А.А., Сенкевич С.Е., Дурягина В.В., Гончаров Д.А.Патент на изобретение RU 2398147 C1 // Устройство для снижения жесткости трансмиссии машинно-тракторного агрегата //Бюл. №24, 2010.
3. А.Б. Лурье Статистическая динамика сельскохозяйственных агрегатов. - 2-е изд., перераб. и доп.- М.: Колос, 1981. - 382 с.
4. Кутьков Г.М. Тяговая динамика тракторов // Г.М. Кутьков. - М.: Машиностроение, 1980. - 216 с.
5. Кравченко В.А., Дурягина В.В., Гамолина И.Э. Уменьшение техногенного воздействия на почву ходовой системой машинно-тракторного агрегата // Труды
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
13
VIII Международной научно-практической конференции «Дни науки - 2012», Чехия, г. Прага, март 2012, с. 11-13.
6. Рыков С.П. Моделирование случайного микропрофиля автомобильных дорог// С.П. Рыков, Р.С. Бекирова, В.С. Коваль [электронный ресурс] http://brstu.ru/static/unit/iournal system/docs/number8/33-37.pdf
7. Прохоров С.А. Аппроксимативный анализ случайных процессов //Самарский государственный аэрокосмический университет, 2001. - 329 с.
8. Белецкий А.В. Моделирование профиля дорожного основания в задаче анализа динамики трансмиссии колесной машины/ А.В. Белецкий // Электронное издание СДМ - Строительные Дорожные Машины и Техника [электронный ресурс]. - Москва:http://sdm.str-t.ru/insertflles/5.pdf
9. Карташов В.Я., Новосельцева М.А. Цифровое моделирование стационарных случайных процессов с заданной корреляционной функцией на основе непрерывных дробей / В .Я. Карташов, М.А. Новосельцева // УБС, 31 (2010), 4991 с.
10. Степанов В.Е.Влияние колебаний МТА на его энергетические показатели (на примере пахотного агрегата К-701+ПТК-9-35) / В.Е.Степанов: Дис... к. техн. наук. - М., 1985. - 158 с.
11. Кравченко, В.А. Исследование влияния упругодемпфирующего механизма на показатели работы посевного машинно-тракторного агрегата / В.А. Кравченко, А.А. Сенькевич, С.Е. Сенькевич // Механизация и электрификация сельского хозяйства, 2007. - № 11, 21 - 23с.
References
1. Kravchenko, V.A. Povyshenie dinamicheskih i jekspluatacionnyh pokazatelej sel'skohozjajstvennyh mashinno-traktornyh agregatov: monografija / V.A. Kravchenko. - Zernograd: FGOU VPO AChGAA, 2010. - 224 s.
2. Kravchenko V.A., Senkevich A.A., Senkevich S.E., Durjagina V.V., Goncharov D.A.Patent na izobretenie RU 2398147 C1 // Ustrojstvo dlja snizhenija zhestkosti transmissii mashinno-traktornogo agregata //Bjul. №24, 2010.
3. A.B. Lur'e Statisticheskaja dinamika sel'skohozjajstvennyh agregatov. - 2-e izd., pererab. i dop.- M.: Kolos, 1981. - 382 s.
4. Kut'kov G.M. Tjagovaja dinamika traktorov // G.M. Kut'kov. - M.: Mashinostroenie, 1980. - 216 s.
5. Kravchenko V.A., Durjagina V.V., Gamolina I.Je. Umen'shenie tehnogennogo vozdejstvija na pochvu hodovoj sistemoj mashinno-traktornogo agregata // Trudy VIII Mezhdunarodnoj nauchno-prakticheskoj konferencii «Dni nauki - 2012», Chehija, g. Praga, mart 2012, s. 11-13.
6. Rykov S.P. Modelirovanie sluchajnogo mikroprofilja avtomobil'nyh dorog// S.P. Rykov, R.S. Bekirova, V.S. Koval' [jelektronnyj resurs]
http://brstu.ru/ static/unit/j ournal_system/docs/number8/3 3 -37.pdf
7. Prohorov S.A. Approksimativnyj analiz sluchajnyh processov //Samarskij gosudarstvennyj ajerokosmicheskij universitet, 2001. - 329 s.
8. Beleckij A.V. Modelirovanie profilja dorozhnogo osnovanija v zadache analiza dinamiki transmissii kolesnoj mashiny/ A.V. Beleckij // Jelektronnoe izdanie SDM -Ctroitel'nye Dorozhnye Mashiny i Tehnika [jelektronnyj resurs]. - Moskva:http://sdm.str-t.ru/insertfiles/5.pdf
http://ej.kubagro.ru/2014/07/pdf/24.pdf
Научный журнал КубГАУ, №101(07), 2014 года
14
9. Kartashov V.Ja., Novosel'ceva M.A. Cifrovoe modelirovanie stacionarnyh sluchajnyh processov s zadannoj korreljacionnoj funkciej na osnove nepreryvnyh drobej / V.Ja. Kartashov, M.A. Novosel'ceva // UBS, 31 (2010), 49-91 s.
10. Stepanov V.E.Vlijanie kolebanij MTA na ego jenergeticheskie pokazateli (na primere pahotnogo agregata K-701+PTK-9-35) / V.E.Stepanov: Dis... k. tehn. nauk. -M., 1985. - 158 s.
11. Kravchenko, V.A. Issledovanie vlijanija uprugodempfirujushhego mehanizma na pokazateli raboty posevnogo mashinno-traktornogo agregata / V.A. Kravchenko, A.A. Sen'kevich, S.E. Sen'kevich // Mehanizacija i jelektrifikacija sel'skogo hozjajstva, 2007. -№ 11, 21 - 23s.
http://ej.kubagro.ru/2014/07/pdf/24.pdf
