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

Раздел 1. Наземные транспортные средства, энергетические установки и двигатели, ной трансмиссии за счет комбинированного способа регулирования гидромашин. М., «Вестник машиностроения», № 2, 2006, с. 27-32.

2. Лепешкин A.B., Михайлин A.A., Шейпак A.A. Гидравлика и гидропневмопривод. Учебник. Часть 2. Гидравлические машины и гидропневмопривод. Под редакцией A.A. Шей-пака. - М.: МГИУ, 2003. - 352 с.

3. Лепешкин A.B. Математическая модель многоприводной колесной машины в общем случае ее движения. Сборник избранных докладов 49-ой Международной научно-технической конференции Ассоциации автомобильных инженеров (ААИ) России «Приоритеты развития отечественного автотракторостроения и подготовки инженерных и научных кадров» 4-го Международного научного симпозиума, посвященного 140-летию Московского государственного технического университета «МАМИ». Книга 3, М., МГТУ «МАМИ», 2005г., с. 138-158. ISBN 5-94099-036-3.

4. Лепешкин A.B. Математическая модель движения автомобиля «Гидроход-49061» с гидрообъемной трансмиссией в условиях недеформируемой опорной поверхности. Свидетельство о регистрации электронного ресурса ОФАП № 16157 от 06.09.2010.М., ВНТИЦ, 2010, № 50201001473.

5. Курмаев Р.Х., Лепешкин A.B. Повышение точности математической модели движения колесной машины на основании использования результатов ее испытаний. Известия МГТУ «МАМИ». Научный рецензируемый журнал. - М., МГТУ «МАМИ», № 1 (7), 2009, с. 46-56.

6. Городецкий К.И. Механический кпд объемных гидромашин. - М., Вестник машиностроения, 1977, №7.-с. 11-13.

7. Лепешкин A.B. Математическая модель объемного гидропривода вращательного движения. Свидетельство о регистрации электронного ресурса ОФАП № 11446 от 29.08.2008. М., ВНТИЦ, 2008, № 50200801936.

8. Лепешкин A.B. Определение коэффициентов потерь в формулах оценки объемного кпд роторной гидромашины. Свидетельство о регистрации электронного ресурса ОФАП № 11198 от 30.07.2008. М., ВНТИЦ, 2008, № 50200801681.

9. Лепешкин A.B. Определение коэффициентов потерь в формулах оценки механического кпд роторной гидромашины. Свидетельство о регистрации электронного ресурса ОФАП № 11446 от 30.07.2008. М., ВНТИЦ, 2008, № 50200801682.

10.Лепешкин A.B. Методика разработки СААУ трансмиссий многоприводных колесных машин. Монография. Издательство «LAP LAMBERT Academic Publishing GmbH & Co. KG» (Германия). 2011. 102 с. ISBN 978-3-8454-3161-1.

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

прямолинейном качении

к.т.н. проф. Лепешкин A.B., Петров С.Е.

МГТУ «МАМИ» (495) 223-05-23, доб. 1426, /ер a iiiaiiii.ru

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

Ключевые слова: взаимодействие эластичного колеса с недеформируемой опорной поверхностью, математическая модель при установившемся прялю-

Раздел 1. Наземные транспортные средства, энергетические установки и двигатели. линейном качении

В данной статье приведены результаты научных исследований, полученные при выполнении Гранта № 14.В37.21.0290 с финансовой поддержкой государства в лице Ми-нобрнауки России.

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

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

Исходя из опыта создания подобных математических моделей [3], для решения поставленной задачи необходимо, чтобы в ней по известной вертикальной нагрузке Р2 на ось колеса, продольной скорости V х движения его оси и величине продольного скольжения Л' г в пятне его контакта с опорной поверхностью определялись продольная составляющая Рх усилия, передаваемого колесом на остов машины, и момент Мк, возникающий на оси этого колеса.

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

Из рисунка 1 видно, что колесо вращается с угловой скоростью С0К, а его ось имеет продольную скорость Vх. В данной математической модели взаимосвязь между этими параметрами определяется на основании рекомендаций, приведенных в работе [2], и имеет вид:

Кт=сок-гк-

где: Гк - радиус качения эластичного колеса без учета скольжения, для его определения используется формула, предложенная Е. А. Чудаковым:

гк = гк0-1к-Мк,

- коэффициент тангенциальной эластичности шины, определяется для данной шины в соответствии с рекомендациями, приведенными в работе [2];

Гк0 - радиус качения эластичного колеса в ведомом режиме, определяется для данной шины в соответствии с рекомендациями работы [2].

Следует отметить, что формулы (1) и (2), а также рекомендации по определению и

Гк0, приведенные в работе [2], получены независимо от качества опорной поверхности, по

которой катится колесо. То есть они справедливы при качении эластичного колеса как по не-деформируемой, так и по деформируемой опорным поверхностям.

Из рисунка 1 также видно, что колесо через ось воспринимает продольную Рх и нормальную Р2 к опорной поверхности составляющие усилия от остова машины, а также поступающий от привода момент М к .

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

их контакта.

К

Рисунок 1 - Расчетная схема взаимодействия эластичного колеса с грунтом

В соответствии с рекомендациями, высказанными в работе [1], в данной модели пятно контакта шины с грунтом условно рассматривается состоящим из трех областей (рисунок 1):

• первая область - на виде сбоку основная часть этой области представлена дугой АВ в передней части контакта. На видах сверху (в зависимости от соотношения между величиной нормальной деформации грунта Ит и высотой грунтозацепа /?г, протектора шины) эта область обозначена цифрой I. Здесь же обозначено: ¿пр - ширина протектора шины, Ъ -

ширина образующейся колеи;

• вторая область - на виде сбоку эта область представлена горизонталью ВС. На видах

сверху эта область обозначена цифрой //; • третья область - образуется при наличии контакта боковой поверхности шины с грунтом в случае, когда величина нормальной деформации грунта hr превышает высоту грунто-зацепа hri шины. На соответствующем виде сверху эта область обозначена цифрой III. Кроме этого данная область показана на виде спереди колеса.

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

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

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

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

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

Таким образом, учитывая то, что часть нормальной нагрузки, приходящейся на ось колеса со стороны остова машины, реализуется в первой и третьей областях пятна контакта в виде силы R'z, нормальная деформация шины hz в модели вычисляется по величине реакции Rz, определяющейся из уравнения равновесия колеса в вертикальном направлении по формуле:

Отсюда, величина hz, в соответствии с рекомендациями [2], определяется по формуле:

где: Кг - коэффициент, постоянный для данной шины, определяется по данным [2]; рм, - давление воздуха в шине.

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

Rz = Pz - R'z.

К =-,

1 + Ю-pw

Раздел 1. Наземные транспортные средства, энергетические установки и двигатели, ная работа равна суммарной реакции деформируемого грунта в направлении нормали к опорной поверхности. Учитывая принятые допущения о том, что эта сила формируется в первой и третьей областях пятна контакта (рисунок 1), и рекомендации [1] Р'2 в принятом случае определяется по формуле:

К=кг-Ьщ>-\р{Ь)-М + {\-кг)-Ъщ>- | р(И)-М +

2

р

3

где: к 2 ~ коэффициент насыщенности протектора шины;

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

нормальной деформации грунта /?(/?), из уравнений (3), (4) и (5) определяются величины

нормальных деформаций шины Н2 и грунта Ит в данных условиях качения колеса.

Далее рассмотрим работу эластичного колеса во второй области пятна контакта (поверхность ВС на рисунке 1).

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

/ш=/о+/м,

где: /о ~~ коэффициент сопротивления качению эластичного колеса в ведомом режиме. На основании рекомендаций, приведенных в работе [2], равен:

/о--:->

1 + Рм.

ОС у и Р у - коэффициенты, определяющиеся для данной шины в соответствии с рекомендациями, приведенными в работе [2]; /ш - дополнительная составляющая коэффициента сопротивления качению эластичного колеса, нагруженного крутящим моментом. По данным [2] определяется по формуле:

К  V  V

'кО к

Как отмечалось выше, в пределах второй области пятна контакта (поверхность ВС) за счет взаимодействия протектора шины (грунтозацепов) с грунтом формируется также продольная составляющая Ях усилия, которая, как известно [1, 2], в общем случае определяется по формуле:

=

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

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

• если Аг >/?гз, то Ях =5Я -\_к2 -т^ +

2 • /

• если Аг < Л13, то Ях = к2 -т^ +ИТ -Ьир----тср,

1 г

здесь: Б ц - площадь второй области пятна контакта (поверхность ВС на рисунке 1);

Т^ - напряжения трения в контакте наружней поверхности грунтозацепов с опорной поверхностью, определяется в соответствии с рекомендациями [1];

ТСр - напряжения сопротивления срезу грунта грунтозацепами во впадине между ними, определяется в соответствии с рекомендациями [1];

1\ ~ половина длины пятна контакта ВС., определяющаяся из формулы (17);

12 - шаг грунтозацепов шины.

Следует отметить, что результатом взаимодействия эластичного колеса с грунтом в первой и третьей областях пятна контакта, кроме силы 11'2, является также продольная сила, обозначенная на рисунке 1 Я'х. По сути своей эта сила создает дополнительное сопротивление, возникающее при качении эластичного колеса по деформируемому основанию. В соответствии с рекомендациями [1] сила Я'х создается составляющей усилия от нормальной деформации грунта (образование колеи), которая характеризуется коэффициентом /г сопротивления качению эластичного колеса по грунту, обусловленным потерями на его деформацию, и силой от нагребания грунта перед колесом (так называемое бульдозерное действие [1]). Для оценки этого дополнительного сопротивления в соответствии с рекомендациями [1] используется специальный коэффициент _/б сопротивления качению колеса по деформируемой поверхности от нагребания грунта перед колесом.

Отсюда в математической модели принято: Я'х =

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

^ =

В соответствии с рекомендациями, высказанными в работе [1], принято, что линия действия вектора равнодействующей сил Я'2 и Я'х проходит через ось вращения колеса, образуя некоторый угол Рт с опорной поверхностью. Это допущение, по утверждению [1], может быть принято потому, что поверхность контакта колеса с грунтом в первой области (дуга АВ) по форме близка к боковой поверхности кругового цилиндра с осью, совпадающей с осью колеса. Следствием этого является то, что эта равнодействующая не создает дополнительного момента вокруг оси колеса.

Отсюда следует, что Я'2 = Я'х  .

Для определения величины угла Р т при моделировании делается допущение о том, что линия действия вектора равнодействующей образует с прямой ОВ центральный угол, равный

а /

у^ , где ОС х ~ центральный угол, определяющийся величиной дуги АВ. Это допущение

основано на том, что в пределах первой области пятна контакта (дуга АВ) эпюра распределения относительной величины нормального к опорной поверхности усилия по форме близка к треугольной (см. рисунок 1).

Из геометрических соображений угол ОС х определяется по формуле:

ос т = arccos

Ггс-(И2+Ит)л v

V с У

-а,

где: Гс ~ свободный радиус эластичного колеса, определяется в соответствии с рекомендациями [2], а на рисунке 1 равен отрезкам ОА и ОВ\ hz ~ нормальная к опорной поверхности деформация шины; hT ~ нормальная к опорной поверхности деформация грунта;

ОС ~ центральный угол, определяющийся из соотношения:

ос = arcsm

'll + cЛ

Y ' с

¡1 ~ половина длины пятна контакта ВС (вторая область), определяющаяся из формулы:

1л=^{2-гс-Н2-Н1)-

С ~ дополнительный снос вертикальной реакции во второй области при работе эластичного колеса в ведущем режиме качения.

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

с = Хк-М к-ф.т.

Отсюда величина угла Р т, в соответствии с рисунком 1 определяется из формулы:

Р.т =~-

Таким образом, приведенные зависимости (6) - (19) позволяют в данных условиях качения эластичного колеса по деформируемому основанию при известной величине реализуемого на оси колеса моменте Мк оценить продольную составляющую усилия Рх, воспринимаемую остовом машины через ось рассматриваемого колеса.

Момент Мк, реализуемый на оси рассматриваемого ведущего колеса, определяется из уравнения его мощностного баланса [1, 2]. В соответствии с рекомендациями, приведенными в работе [1], при качении эластичного колеса по деформируемому основанию это уравнение имеет вид:

Ик = N. + N. + ЫГш + ЫГт + + N.,

где: = Мк • С0К - мощность, подведенная к колесу от привода;

Мт = Рх  Ух = (ф т — /т — /б ) • Р2  Гк • СОк • (1 — 'х ) - тяговая мощность колеса;

Nц = Рх • V8 = (ф т — /т — Уб ) • Р2 • Гк • СО к • ях - мощность потерь на буксование;

М/ш =  Я2 • Гк • СОк - мощность потерь на качение эластичной шины (обусловлены потерями на нормальную и окружную деформацию шины);

= /, • 1}2  Ух - мощность потерь на деформацию грунта колесом;

А'г, = Уб • Р2 • Ух - мощность потерь от нагребания грунта перед колесом (бульдозерное

действие), оценивается на основании рекомендаций работы [1]; Ыа = /д • Рг • Гк • Юк - мощность потерь, обусловленная липкостью грунта, оценивается на основании рекомендаций работы [ 1 ]. После алгебраических преобразований с учетом приведенных выше формул и рекомендаций, приведенных в работе [2], получим:

М =iVr,

кО

Ф.г + /о

1

Р.

(/г + /б ) • sx + /л

Z У

Анализ формулы (21) показывает, что в ней используются только те параметры, которые участвовали при получении зависимостей (1) - (19), а это значит, что формула (21) замыкает систему уравнений, характеризующих работу эластичного колеса при его взаимодействии с деформируемой опорной поверхностью. То есть приведенная система уравнений (1) - (21) является математическим описанием процессов, эту работу характеризующих.

0,50

Кт [-1

-0,10

-суглинок

-заболоченный луг

Рисунок 2 - Зависимость коэффициента тяги колеса по сцеплению kT = ф т — fT — f6 от

продольного скольжения S х в пятне контакта

На основании этого математического описания была разработана соответствующая математическая модель в виде программного модуля, написанного на языке Fortran в среде Compaq Visual Fortran под Windows.

На рисунках 2 и 3 в качестве примера приведены графики взаимозависимостей параметров, характеризующих взаимодействие эластичного колеса марки 1300*530-533 с тремя вариантами горизонтальных деформируемых оснований (с условными названиями: «суглинок», «песок» и «заболоченный луг»), полученные в результате расчета по разработанной математической модели. Параметры эластичного колеса марки 1300*530-533 взяты из работы [2], а характеристики опорных поверхностей выбраны из данных, приведенных в работе [1]. Условия функционирования задавались следующими параметрами: давление воздуха в шине - pw = 0,3 МПа, вертикальная составляющая нагрузки на ось колеса - Pz = 20 кН, угловая скорость вращения колеса - С0К = 2,24 рад/с, а также диапазоном изменения продольной величины скольжения в пятне контакта колеса с опорной поверхностью - S х =0 ... 0,5.

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

Рисунок 3 - Взаимозависимость между суммарной продольной составляющей Рх усилия, создаваемой двухколесной осью, от подведенного к ней через трансмиссию

суммарного момента Мк

Литература

1. АгейкинЯ.С. Проходимость автомобилей. М., «Машиностроение», 1981, 231с.

2. Петрушов В.А., Московкин В.В., Евграфов Ф.Н. Мощностной баланс автомобиля. М., «Машиностроение», 1984, 160с.

3. Свидетельство об отраслевой регистрации разработки № 11447 от 29.08.2008г. «Математическая модель работы эластичного колеса при его взаимодействии с недеформируемой опорной поверхностью». Авт.: C.B. Бахмутов, A.B. Лепёшкин.

Топливный насос нового поколения

д.т.н. проф. Максимов Ю.В., д.т.н. проф. Соковиков В.К., к.т.н. доц. Бекаев A.A.,

Строков П И.

Университет машиностроения bekaev(a)Jist.ru

Аннотация. На основе электрогидравлического эффекта, также известного как «эффект Юткина», была разработана принципиально новая конструкция электрогидравлического топливного насоса (ЭГ-насос) высокого давления, обеспечивающего превосходящие (по давлению и расходу) рабочие характеристики современных топливных насосов высокого давления двигателей внутреннего сгорания (ДВС). Применение разработанного ЭГ-насоса в ДВС позволит полностью отказаться от использования какого-либо механизированного привода, являющегося неотъемлемой частью современных топливных насосов, а значит, существенно упростить конструкцию и принцип работы топливной системы, уменьшить ее металлоемкость, стоимость и уровень шума.

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

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