НАЗЕМНЫЕ ТРАНСПОРТНЫЕ СИСТЕМЫ
УДК 629.11.012
Р.И. Жирный1, В.А. Горелов2, Г.О. Котиев2
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРЯМОЛИНЕЙНОГО ДВИЖЕНИЯ АВТОПОЕЗДА С ЖЕСТКОЙ И ГИБКОЙ СВЯЗЬЮ МЕЖДУ ЗВЕНЬЯМИ
ООО «Газсеверинвест»1, Московский государственный технический университет им. Н.Э. Баумана
Представлена математическая модель прямолинейного движения автопоезда по деформируемому грунту. Описаны модели взаимодействия между звеньями для случаев упругой и гибкой связи.
Ключевые слова: автопоезд, проходимость, сцепное устройство, буксование, тягово-энергетическая характеристика, математическая модель, имитационное моделирование, трансмиссия.
Большинство разрабатываемых месторождений нефти и газа находятся в труднодоступных районах Крайнего Севера, в условиях вечной мерзлоты и бездорожья, в связи с чем, доступ к ним весьма затруднён. Эти месторождения, в том числе и на п-ове Ямал, требуют обслуживания и доставки тяжёлых крупногабаритных и неделимых грузов. Большие перспективы для решения подобных транспортных задач имеет направление, связанное с применением многоосных большегрузных автопоездов.
В рамках совместных исследований кафедры «Колёсные машины» МГТУ им. Н.Э. Баумана и ООО «Газсеверинвест» была создана транспортная система (рис. 1), прошедшая предварительные испытания на п-ове Ямал в апреле 2009 года.
Доказано, что в зависимости от характера связи между тягачом и прицепом автопоезд может иметь различные характеристики по проходимости. В связи с этим, необходимо на стадии проектирования решить задачу по прогнозированию характеристик прямолинейного движения автопоезда по бездорожью при буксировании прицепа движущимся тягачом с жесткой связью, буксировании прицепа движущимся тягачом с гибкой связью (односторонней на тросе), буксировании прицепа неподвижным тягачом с использованием лебёдки и анкер-но-упорного устройства. Данная задача решается с помощью имитационного математического моделирования.
Вопросам взаимодействия эластичного колёсного движителя и деформируемого опорного основания посвящено достаточно много трудов как отечественных, так и зарубежных авторов. Прежде всего эти работы направлены на описание процесса взаимодействия колеса с грунтом в функции конструктивных параметров движителя и физико-механических свойств опорного основания. Следует отметить, что теория взаимодействия эластичного колеса с деформируемым опорным основанием в данном представлении весьма развита и позволяет описывать такие явления, как колееобразование, бульдозерный и экскавационный эффекты, уплотнение грунта и т.д. При таком подходе целью исследований является выбор оптимальных параметров конструкции колеса, как правило, по критериям максимального тягового усилия, минимизации затрат на движение и вредного воздействия на почву.
© Жирный Р.И., Горелов В.А., Котиев Г.О., 2010.
Ьук< и|ккиа пришчи ил жестком ЛЫШЧ' с мгпохьхтлнмсм ВТС (> Ь VП К\Ч Ткмммиь-Теммк* > химейммй гигач). Режим '1[»|\н'ин<*'( н пр«( лнижгннн по лмкиормглм
Совчкч ки< 6\к< иропкл при«Н" 1.1 ВТС "Птмрник м и(х юшчк" и ВТС "11илнркнк-Т« \ник" мл жестком ЛМШЛ4*
('(•ЖИМ ПОММ«.'Н»|'Т1 Н II [Н1 М.1 УЧ.Ц IK.1V \Upitip4l.l С ковыип-нным ( |М1р<>ГИИЛ| НЬЫ'М .ШИЖСНИИ ЧИ'рсМЧ'ТМ М'И к ы
ийОЛОК.'ЧНЫ* уч.И
Рис. 1. Транспортные системы
Однако когда движитель уже создан и необходимо исследовать его возможности в составе колёсной машины, например, в случае оценки средней скорости движения на заданном маршруте, с целью решения задачи о рациональном распределении мощности по колёсам для оценки нагруженности элементов трансмиссии представляется целесообразным использовать характеристики тягово-сцепных свойств (рис. 3) и потерь энергии (рис. 2) для различных условий взаимодействия колеса с опорным основанием, полученные при экспериментальных исследованиях в грунтовом канале или на полигоне (рис. 4 - рис. 7).
Характеристиками процесса прямолинейного качения колеса в различных случаях движения являются удельные потери энергии (потери энергии при качении на единицу пройденного колесом пути при единичной вертикальной нагрузке), удельная свободная тяга (продольная сила, приложенная к оси катящегося колеса, при единичной вертикальной
нагрузке на его ось), коэффициент буксования (юза) = к ''к°-— = 1 - —; гк= [1],
®к Гко Гко ®к
где соА, - угловая скорость вращения колеса; гко - радиус колеса в свободном режиме; гк -радиус качения колеса; Ухк - скорость центра масс колеса в продольном направлении.
Результаты испытаний представляются в виде тягово-энергетических /и. = /(ср ) и тягово-сцепных характеристик ср = /(8Б ). Примеры экспериментальных характеристик, виды грунтового канала и экспериментальная установка представлены на рис. 2 - рис. 7 [1, 3].
Рис. 2. Удельные потери энергии при качении модификаций колеса
,11*
Мгов & ¿•Я7-М9«
ъ ,/Л/
и ШЛ • тич • ' . 1Г 1
а ¿У
* } « « * г
№
Рис. 3. Тягово-сцепные характеристики модификаций колеса
Рис. 4. Цилиндрическая модификация колеса
Рис. 5. Образец в грунтовом канале
а б
Рис. 6. Общий вид экспериментальной установки для исследования взаимодействия с грунтом крупногабаритного пневматического колеса при нагрузке 100-250 кН:
а - исследование ведомого колеса; б - опыты с ведущим колесом
а б
Рис. 7. Экспериментальная лабораторная установка для исследования взаимодействия эластичных оболочек с очень слабыми пылевыми и мягко-текучими грунтами:
а - большой канал с пылевидным грунтом; б - малый канал с мягко-сыпучим грунтом
Методика проведения экспериментальных исследований с целью получения представленных характеристик опирается на уравнение энергетического баланса и подробно описывается в [2]. Построение математической модели движения автопоезда рассмотрим на примере трехосного тягача с возможностью реализации индивидуального, дифференциального и блокированного привода колёс и трехосного прицепа в ведомом режиме. Расчётная схема автопоезда, представленная на рис. 8, принята с учётом характера решаемой задачи и описывается соответствующими уравнениями движения для тягача (1) и прицепа (2).
Рис. 8. Расчетная схема движения автопоезда
Для тягача
Для прицепа
где I = 1, ...3; тТ
тт ■ Ух = 2^РХг чп^зш(а) -Р^
г=1
4 А* =МКг -О-^ОСС, +Фг)&2гГКо
тКгаХКг = ,
[КР
1=4
(1)
(2)
тТ
4 А» =(1-^г)(ф/
МКгаХКг = РХг ~ ~ ™юё 8Ш(а) ,
Р - масса тягача и прицепа соответственно; тт - масса колеса; !к
момент инерции колеса; / - номер колеса; Ух , Ухш - продольное ускорение центра масс тягача и прицепа соответственно; со^ - угловое ускорение ьго колеса; аХк - продольное ускорение центра масс ¿-го колеса; g - ускорение свободного падения; а - угол наклона опорной поверхности.
Связь тягача и прицепа реализуется через силу тяги на крюке Р№ . В случае моделирования упругого или жесткого взаимодействия звеньев сила на крюке определяется по зависимости
^КР = (*ст - ^спр)С + (^ХГ - , (3)
где С - коэффициент жесткости сцепного устройства в продольном направлении; В - коэффициент демпфирования сцепного устройства в продольном направлении; ХСт - коорди-
ната центра масс тягача по оси Х; ХСт - координата центра масс прицепа по оси Х; УХТ -скорость центра масс тягача; УХПР - скорость центра масс прицепа.
Таким образом, изменяя коэффициенты жесткости и демпфирования узла сочленения, можно задавать различную связь между тягачом и прицепом.
Для случая соединения звеньев через гибкую связь характер изменения силы тяги на крюке Ркр в зависимости от растяжения троса А/ = ХСт - ХСш, определяется по зависимости (4) с учетом начальных координат положения центров масс звеньев:
РЛР=(А/-4)С,при Д/<0, Ркр =0, (4)
где Ь0 - длина троса в недсформированном состоянии; А/ = (ХСт -ХСпр -¿а).
Сила, действующая на корпус тягача со стороны колеса по оси X
Рх, = (^К, - ¿К, )СПОдв + (^Л-К, - Кхс )5подв , (5)
где Сподв - коэффициент жесткости подвески в продольном направлении; Рподв - коэффициент демпфирования подвески в продольном направлении; ХК{ - расстояние от центра масс до оси колеса по оси Х; ЬК{ - расстояние от центра масс до точки крепления подвески по оси Х; ¥Хс - скорость центра масс автомобиля.
Сила, действующая на корпус прицепа со стороны колеса по оси X
Рх, = (¿к, " )СПОдв + (Ухе - Ухк, )5подв- (6)
Пусть удельная свободная сила тяги соответствует рис. 9 и определяется по зависимости
Л (7)
Ф =«g/?(^)-9max -(1-е 0 )
где ф - максимальное значение коэффициента силы сцепления колеса с опорной поверхностью; - константа.
Удельные потери энергии при качении колеса (рис. 10)
108(1--(8)
^'Шах (1 ^ )
где /м,о - коэффициент удельных энергетических потерь при отсутствии буксования в свободном режиме качения; К - коэффициент пропорциональности.
Рис. 9. Зависимость ф(5Б) Рис. 10. Зависимость f =/(ф)
Вертикальная реакция в пятне контакта колес тягача с опорной поверхностью определяется следующей зависимостью:
RZi=pZi+mKigz os(a) (9)
PZi определяются из решения системы уравнений (9):
рг\ + //2 + ргз = тс8 соз(а), Р2]ЬК1 + Р22ЬК2 + ХМош (10)
~ ^КЗ ) + Р22 (4ГЗ ~ ) + Р2Ъ (4г1 ~~ 4:2 ) = ^ ' где ^Мотн =-{тс -ах -Нс +Ркр -Нкр +Мк(1) +Мк{2) +^-(3)) " сУмма крутящих моментов
относительно проекции центра масс на линию, соединяющую оси колес; Ис - высота центра масс; И^ - высота приложения силы тяги на крюке Ркр ; М гЛ1- крутящий момент, подводимый к соответствующей оси автомобиля.
При определении вертикальных реакций под колёсами прицепа используется такой же подход, как и для тягача, с учетом особенностей, отраженных в расчетной схеме.
Динамика трансмиссии рассмотрена на примере тягача с блокированной раздачей мощности по колёсам каждого борта и описывается системой дифференциальных уравнений с учетом равенства угловых скоростей колес.
Разработанная математическая модель реализована в среде имитационного математического моделирования динамических систем БтиНпк программного комплекса МайаЬ. Схема включает три основные блока: Tyagach, в котором реализуется движение автомобиля-тягача (рис. 12), РтШер, в котором реализуется движение прицепа, Stsepka, в котором реализуется связь между звеньями автопоезда через силу на крюке. Блок-схема модели представлена на рис. 11.
Рис. 11. Блок-схема математической модели
Выводы
Проведенные предварительные теоретические исследования динамики автопоезда при различных упругих и демпфирующих характеристиках сцепного устройства позволили сделать следующие выводы. В том случае, когда параметры сцепного устройства С и В принимаются равными небольшим величинам, тягач «уезжает» от прицепа, скорость которого в определенный момент времени становится равной нулю. В противоположном случае, т.е. когда коэффициенты упругости и демпфирования достаточно велики, тягач оказывается неспособным выполнять транспортировку прицепа. Скорость обеих частей автопоезда в определенный момент равна нулю. При определенных значениях величин С и В оказалось возможным выполнить транспортировку тягачом прицепа при тех же массовых, геометрических и энергетических характеристиках. При начале движения за счёт имеющейся упругой связи тягач сначала разгоняется без приложения силы на крюке, а далее полученной инерции хватает для того, чтобы потянуть за собой прицеп.
Рис. 12. Схема блока для расчета динамики автомобиля-тягача
На основании полученных результатов был сделан вывод о том, что жесткость связи между звеньями оказывает существенное влияние на динамику автопоезда. Таким образом, сравнительная оценка возможностей автопоезда при различной связи между звеньями должна быть проведена ещё до натурных испытаний объекта. Результатом этих исследований являются рекомендации по выбору оптимальных параметров сцепного устройства для повышения уровня опорной проходимости.
Библиографический список
1. Наумов, В.Н. Математическая модель взаимодействия металлоупругого колеса с уплотняющимся грунтом / В.Н. Наумов, Ю.Л. Рождественский // Труды МВТУ им. Н.Э. Баумана, 1980. № 339. С. 84-111.
2. Серебренный, И.В. Повышение проходимости автомобиля за счёт рационального распределения потоков мощности по колёсам / И.В. Серебренный, Г.О. Котиев // Вестник МГТУ им. Н.Э. Баумана. - М.: Машиностроение, 2008. С. 193-201.
3. Кошарный, Н.Ф. Основы теории рабочего процесса и расчета движителей автомобилей высокой проходимости / Н.Ф. Кошарный: дисс. .. .докт. техн. наук: 05.05.03. - Киев, 1978. - 467 с.
Дата поступления в редакцию 25.05.10.
R.I. Zhirnyu, V.A. Gorelov, G.O. Kotiev
MATHEMATICAL MODEL OF RECTILINEAR DYNAMICS OF THE ARTICULATED LORRY WITH ELASTIC AND FLEXIBLE CONNECTION BETWEEN SECTIONS
In article the mathematical model of rectilinear dynamics of the articulated lorry on a deformable ground is presented. Interaction models between sections for cases of elastic and flexible connection are described.
Key words: articulated lorry, passability, coupling, slip, traction-energy characteristic, mathematical model, simulation modeling, transmission.