ЭНЕРГЕТИЧЕСКАЯ ОЦЕНКА ПРОЦЕССА КАЧЕНИЯ ПНЕВМАТИЧЕСКОГО  КОЛЕСА ПО ДЕФОРМИРУЕМОМУ ГРУНТУ

 

Синицын С.С., Синицына Т.С.  (БГИТА, г. Брянск, РФРФ)

 

They Are Received accounting dependencies, allowing conduct the energy estimation of the interaction of the wheel machines with soil.

 

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

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

Лесотранспортные машины работают главным образом  в лесисто- болотистой местности, где имеет место сопоставимость деформаций эластичного колеса и почво-грунтов. Причем, поверхность контакта весьма далека от цилиндрической  и представляет из себя совокупность  плоской и криволинейной  (в плоскости Z, Y) зон.

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

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

Заменяя действительный путь частицы грунта по циклоиде S его вертикальной составляющей h по оси Z колеса,  принимаем, что реакция  q грунта на элементарную площадку протектора шины будет пропорциональна глубине её погружения в грунт. Сумма произведений этих сил на координаты Х их приложения определяет момент сопротивления  М_f  качению колеса вследствие деформации грунта. Уравнение для  М_f  можно представить в следующем виде

 

 

Рисунок 1 -  Расчетная схема взаимодействия колеса с опорной поверхностью: v -  вектор линейной  скорости;  x,  z  -  декартовы координаты; М_к – крутящий момент, подведенный к колесу от трансмиссии; М_f -  момент сопротивления качению колеса; s – траектория перемещения элементарной частицы грунта из точки 1 в точку 2; q – элементарная реакция грунта; h_Г – глубина колеи;  h_Ш -  деформация шины; R  - свободный радиус колеса; А = R -  h_ш - h_Г.

 

.                                                   (1)

Поскольку  [3], где A  = R – h - ∆ ,  а , то тогда получаем

.                                              (2)

В данных выражениях (1 и 2) Y представляет собой ось координат, перпендикулярную плоскости  рис. 1, а h₀ -  базовую деформацию почво-грунта, численно равную 0,01 м.

Разложив (z-A)^μ в степенной ряд и поставив вместо z его значение из характеристического уравнения поверхности шины, а именно:

, получим 

            (3)

Приняв, что х =R u, а  в = в v,  и   определив   якобиан   преобразования   I=R в >0, для первого интеграла получаем рекуррентное соотношение

Переходя к полярным координатам , , приводим интеграл к следующему виду:

      (4)

После преобразований внутреннего интеграла получаем

а это не что иное, как частный случай интеграла от биноминального дифференциала  [4]

,                                                   (5)

где m = 2, .

Полученный интеграл не подходит ни под один из трёх случаев выражения его через элементарные функции. Тогда, применяя подстановку t = r², преобразуем его к виду

                                           (6)

а  это  есть  частный  случай  интегрального  представления  неполной  бета- функции [4]

                                       (7)

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

F (а, в; с; х) следующим образом [4]:

.                                       (8)

Гипергеометрическая функция, в свою очередь,  может быть выражена через гамма-функцию, а именно [4]:

.                                            (9)

Тогда

.

Взяв внешний интеграл в уравнении (4), получим

.        (10)

Так как F_{(1+ p)} =  р Г_(р),  а Г₍₁₎=1, то приводим правую часть уравнения (10) к следующему виду:

.                                         (11)

Использовав для гамма-функций асимптотическую формулу, получаем ,     .

Тогда

.                      (12)

Решив последующие интегралы уравнения (3), для жесткого колеса окончательно получаем

                                    (13)

Проинтегрировав уравнение (3) в диапазоне от 0 до r₁, определяем потери на деформацию грунта эластичным колесом

                               (14)

С учетом того, что

,  ,  ,

имеем

(15)

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

 .                                          (16)

Сопротивление грунта качению колеса может быть определено если известна работа, затрачиваемая на образование колеи. В свою очередь эта работа пропорциональна перемещению точки поверхности колеса в грунте из положения 1 (рисунок 1) в положение 2 по циклоиде S, которое может быть заменено нормальным перемещением h по оси  z колеса.

С учетом того, что  h_Г = z - A, работа колееобразования за один оборот эластичного колеса описывается следующей зависимостью:

(17)

После интегрирования первого тройного интеграла по  dz и dl получаем

.                     (18)

После замены переменных ,   выражение принимает следующий вид:

.                  (19)

Полученный интеграл  есть частный случай интегрального представления неполной бета-функции В_х (р, q) (уравнение 7), которая  аналогичным  образом (зависимости 8 и 9) может быть выражена через гипергеометрическую функцию F(а, в; с;  х)  и в общем виде через гамма-функцию

 

.

С учетом того , что Г_(р+1) = р Г_(р) ,  Г₍₁₎ = 1, а  получаем

.

Поскольку , а  , то тогда

.

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

.

После преобразований приводим к следующему виду:

.                        (20)

Проведя аналогичные выкладки для второго тройного интеграла, получим

.         (21)

 

С учетом того,  что А= R –h_Г – h_ш, разлагаем выражение в скобках в биноминальный ряд и, ограничившись линейным приближением этого ряда, окончательно получаем

.                                      (22)

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

 

Список использованных источников

1. Наземные тягово-транспортные системы. Энциклопедия. Т 1-2 [Текст] / И.П. Ксеневич, В.А. Гоберман, Л.А. Гоберман. – М.: Машиностроение, 2003. – 743 с.

2. Левин, М.А. Теория  качения деформируемого колеса [Текст] /  М.А. Левин. – М.: Наука, 1989. – 270 с.

3. Синицын, С.С. Математическая модель колеса с эластичной шиной [Текст] / С.С. Синицын // Эксплуатация  лесовозного  подвижного  состава:  Межвуз. сб. науч. тр./ Изд-во УПИ – 1985. - С. 34-41.

4. Янке, Э.Я.Справочник по специальным функциям [Текст] / Э.Я. Янке.  – М.: Наука, 1979. - 832 с.