УДК 656.1
М.М. Жилейкин, А.Ю. Захаров, М.В. Паньшин
ПРОВЕРКА АДЕКВАТНОСТИ И ТОЧНОСТИ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ВЗАИМОДЕЙСТВИЯ ЭЛАСТИЧНОГО КОЛЕСА С ДЕФОРМИРУЕМЫМ ОПОРНЫМ ОСНОВАНИЕМ
Московский государственный технический университет имени Н.Э. Баумана
К настоящему времени известно большое число математических моделей движения колеса по деформируемому опорному основанию, основанных на разных подходах к математическому описанию как оболочки шины, так и самого грунта. При разработке любой математической модели всегда возникают вопросы, связанные с ее адекватностью и точностью, что можно проверить только в случае сравнения вычисляемых в модели параметров с их значениями, полученными в ходе проведения измерений на реальных объектах, то есть экспериментальными методами. Целью работы является проверка точности и адекватности математической модели взаимодействия эластичного колеса с деформируемым опорным основанием, разработанной авторами, экспериментальными методами. Установлено, что относительная ошибка при определении различных параметров взаимодействия шины с деформируемым опорным основанием не превышает 7,4%.
Ключевые слова: деформируемое опорное основание, математическая модель колеса, адекватность и точность математической модели.
Введение
При разработке новых вездеходных автомобилей спрогнозировать их проходимость в многообразных дорожно-грунтовых условиях можно только с помощью методов математического моделирования. Широкое практическое использование этих методов затрудняется малой изученностью процессов взаимодействия эластичного колеса с деформируемым опорным основанием, недостаточным объемом материалов по характеристикам грунтов [1].
К настоящему времени известно большое число математических моделей движения колеса по деформируемому опорному основанию, основанных на разных подходах к математическому описанию как оболочки шины, так и самого грунта [2-8].
При разработке любой математической модели всегда возникают вопросы, связанные с ее адекватностью и точностью, что можно проверить только в случае сравнения вычисляемых в модели параметров с их значениями, полученными в ходе проведения измерений на реальных объектах, то есть экспериментальными методами.
Целью работы является проверка точности и адекватности математической модели взаимодействия эластичного колеса с деформируемым опорным основанием, разработанной авторами [9], экспериментальными методами.
Экспериментальное определение параметров взаимодействия колеса с деформируемым опорным основанием
Целью экспериментального исследования является определение параметров взаимодействия одиночного колёсного движителя с деформируемой опорной поверхностью (форма пятна контакта, глубина колеи, «зарывание» колеса в грунт, деформирование опорной поверхности).
В качестве объекта исследования использовался многоцелевой автомобиль МЗКТ-6001 (рис. 1), технические характеристики которого представлены в табл. 1.
© Жилейкин М.М., Захаров А.Ю., Паньшин М.В., 2018.
Рис. 1. Используемый тягач
Таблица 1
Технические характеристики МЗКТ-6001
Колесная формула 6х6
Снаряженная масса ТС, кг 12800
Полная масса ТС, кг 24000
Допустимые осевые массы, кг 3х8000
Максимальная скорость, км/ч 90
Запас хода по контрольному расходу топлива, км 1000
Емкость топливных баков, л 250+250
Высота ТС, мм 3195
Колея ТС, мм 2550
Двигатель - ЯМЗ-7513:
Максимальная мощность при Нд=1900об/мин, л.с. (кВт) 420 (309)
Максимальный крутящий момент при Нд=1100-1300об/мин, кГс ■ м (Н ■ м ) 185 (1813)
Статический радиус, м (шина - 16.00Я20) 0,67
Для определения нагрузочной характеристики шины экспериментальное исследование проводилось на твердой ровной опорной поверхности - асфальтобетонном дорожном покрытии в сухую погоду на полигоне в Дмитровском филиале МГТУ им. Н.Э. Баумана. При проведении эксперимента регистрировались нагрузка на колесо, свободный радиус шины на твердой ровной опорной поверхности, количество грунтозацепов, их высота, глубина колеи, образовавшаяся в области взаимодействия колёсного движителя с грунтом, расстояние от центра вращения колеса до сдеформированного грунта, а также параметры, характеризующие длину поверхности контакта. Было произведено пятнадцать замеров с дальнейшей обработкой полученных результатов измерений с использованием известных статистических методов.
Замеры проводились на твердой ровной опорной поверхности для левого колеса передней оси автомобиля МЗКТ-6001.
Регистрировались следующие параметры, определяющую нагрузочную характеристику шины (рис. 2):
1) нагрузка на колесо Ск, Н;
2) свободный радиус шины гсв, мм;
3) расстояние от центра вращения колеса до опорной поверхности гкп, мм;
4) количество грунтозацепов пг;
5) высота грунтозацепов К2т, мм.
Замеры нагрузки на колесо были проведены при помощи весоизмерительного устройства - подкладных весов «Исток» 7,5-2.
Результатом измерения являются следующие значения искомых параметров:
1) Ск = 45472,21 Н;
2) т"св = 671,5 мм;
3) 2кп = 615 мм;
4) пг = 128;
5) Л2г = 25 мм.
В процессе измерений производился выезд автомобиля на деформируемую опорную поверхность, где осуществлялось наблюдение за взаимодействием одиночного колесного движителя (левое колесо передней оси автомобиля МЗКТ-6001) с деформируемой опорной поверхностью (рис. 3).
Рис. 3. Взаимодействие колесного движителя с грунтом
При этом необходимо было определить следующие параметры взаимодействия колёсного движителя с деформируемым опорным основанием (рис. 4):
1) длина передней части пятна контакта мм;
2) расстояние по горизонтали от центра колеи до точки контакта профиля колеса с недефор-мированным грунтом Хг, мм;
3) длина задней части пятна контакта Ь2, мм;
4) глубина колеи Лг, мм;
Чтобы определить искомые параметры, производилась расчистка области пятна контакта и их регистрация с помощью измерительных инструментов: линейки измерительной металлической и рулетки измерительной (рис. 5).
% : »' ' \ ' - •. ч V
Рис. 5. Измерение параметров взаимодействия колесного движителя с грунтом
Результаты измерения сведены в табл. 2.
Таблица 2
Параметры взаимодействия колеса с грунтом
№ п/п мм Хг, мм ¿2, мм Лг, мм гк, мм
1 183 180 148 10 603
2 180 170 145 12 606
3 183 181 150 20 610
4 180 178 148 10 605
5 178 173 157 8 605
6 185 183 145 15 605
7 188 184 155 13 603
8 182 175 155 13 605
9 176 170 180 15 608
10 196 190 125 13 607
11 188 180 145 18 609
12 186 180 148 15 605
13 190 180 145 10 610
14 200 190 130 15 610
15 188 180 145 15 608
Обработка результатов экспериментальных исследований
Оценка математического ожидания и дисперсии исследуемых параметров представлена в табл. 3.
Таблица 3
Оценка математического ожидания и дисперсии исследуемых параметров колеса на грунте
Параметр Ьг, мм мм ¿2, мм кг, мм гк, мм
Оценка мат. ожидания, т 185.5 179.6 148.8 13.47 606.6
Оценка дисперсии, Б2 39.45 33.44 18.31 9.58 5.57
Вычисление статистики [10] и проверка первого и максимального членов вариационного ряда, сравнивание полученных значений с табличным критическим значением иа (для уровня значимости q = 0,05 и объема выборки п = 15) по каждому из исследуемых параметров показало, что вся совокупность результатов испытаний по всем исследуемым параметрам, кроме совокупности результатов по параметру длины задней части пятна контакта Ь2, принимается в рассмотрение. Для Ь2 было исключено несколько выбросов из рассмотрения, а оценки математического ожидания, дисперсии и среднеквадратического отклонения были скорректированы для 12 замеров без учета исключённых значений.
Проверка соответствия эмпирической функции распределения теоретическому закону распределения с использованием критерия /2 — Пирсона для группированных наблюдений показала, что измеренные результаты по каждому из исследуемых параметров колеса на грунте подчиняются нормальному закону.
Было установлено, что результаты наблюдений подчиняются нормальному закону, поэтому в качестве наилучшей оценки измеряемых величин принимается оценка математического ожидания выборки [10].
Значения математического ожидания выборки и доверительных интервалов измеряемых параметров представлены в табл. 4.
Таблица 4
Математические ожидания и доверительные интервалы исследуемых параметров колеса на грунте
Параметр Ьг, мм мм ¿2, мм кг, мм гк, мм
Оценка мат. ожидания, т 185.5 179.6 148.8 13.47 606.6
Минимальное зна-
чение доверительного интервала, 181.93 176.29 145.995 11.7 605.25
Максимальное
значение доверительного интервала, ХРв 189.13 182.91 151.67 15.24 607.95
Проверка адекватности математической модели взаимодействия эластичного колеса с деформируемым опорным основанием
Проверка проводилась путем сравнения численных значений параметров взаимодействия колеса с опорным основанием, полученных как экспериментальными методами, так и методами имитационного моделирования движения многоцелевого автомобиля МЗКТ-6001 (рис. 6), по деформируемому грунту, характеристики которого получены в ходе проведенных
ранее замеров. Математическая модель взаимодействия колеса с деформируемым опорным основанием представлена в работе [9], модель движения колесной машины - в работе [11].
В качестве модели деформируемого опорного основания принята модель Бернштейна-Летошнева [1]:
Poz = СгрЛ?, МПа,
где сгр, МПа - коэффициент деформации грунта; ц - показатель плотности грунта.
Ранее экспериментальными методами были получены следующие значения сгр и д для исследуемого опорного основания: д = 0.61211; сгр = 0.19343 МПа
В процессе проведения описанных выше испытаний были получены следующие параметры взаимодействия для переднего левого колеса:
• общая длина пятна контакта + L2 = 334,36 мм;
• расстояние от центра вращения колеса до уплотненного грунта zK = 615 мм;
• глубина колеи Лг = 13,47 мм.
Эти же параметры, полученные при моделировании движения МЗКТ-6001 по ровному деформируемому опорному основанию (рис. 6), представлены на рис. 7а, б, в.
Установившиеся значения параметров взаимодействия с грунтом переднего левого колеса, полученные при моделировании:
• общая длина пятна контакта + L2 = 357 мм;
• расстояние от центра вращения колеса до уплотненного грунта zK = 613 мм;
• глубина колеи h = 14,40 мм.
Относительная ошибка при определении общей длины пятна контакта составила 6,8%, при определении расстояние от центра вращения колеса до уплотненного грунта - 0,3%, при определении глубины колеи - 7,4%.
Рис. 6. Моделирование движения МЗКТ-6001 по ровному деформируемому опорному основанию
0.4 0.35 ............................. I
0.25
р
0
0i_I_I_I_I_I_1_
0 0.05 0.1 0.15 0.2 0.25 0.3
tc
а)
в)
Рис. 7. Зависимости от времени параметров взаимодействия левого переднего колеса МЗКТ-6001 с ровным деформируемым опорным основанием:
а - общая длина пятна контакта; б - глубина колеи; в - расстояние от центра вращения колеса до уплотненного грунта
Заключение
Проведена проверка адекватности математической модели взаимодействия эластичного колеса с деформируемым опорным основанием путем сравнения численных значений параметров взаимодействия колеса с опорным основанием, полученных как экспериментальными методами, так и методами имитационного моделирования движения многоцелевого автомобиля МЗКТ-6001. Установлено, что относительная ошибка при определении различных параметров взаимодействия шины с деформируемым опорным основанием не превышает 7,4%.
Библиографический список
1. Вольская, Н.С. Оценка проходимости колесных машин при движении по неровной грунтовой поверхности / Н.С. Вольская. - М.: МГИУ, 2007. - 215 с.
2. Агейкин, Я.С. Проходимость автомобилей / Я.С. Агейкин, Н.С. Вольская, И.В. Чичекин. - М.: МГИУ, 2010. - 275 с.
3. Бабийчук, А.Э. Методика определения потерь мощности на качение колесного движителя с учетом типа трансмиссии и давления воздуха в шинах машины / А.Э. Бабийчук, А.С. Агейкин, Н.С Вольская // Журнал автомобильных инженеров. - 2013. - № 3. - С. 24-27.
4. Лепешкин, А.В. Математическая модель взаимодействия эластичного колеса с деформируемой опорной поверхностью при установившемся прямолинейном качении / А.В. Лепешкин, С.Е. Петров // Материалы 77-й Международной научно-технической конференции ААИ «Автомобиле- и тракторостроение в России: приоритеты развития и подготовка кадров». - М.: МГТУ «МАМИ», 2012. -С. 141-149.
5. Shoop, S. Finite element modeling of tires on snow / S. Shoop, K. Kestler, R. Haehnel // Tire Science and Technology. - 2006. - Vol. 34, № 1. - P. 2-37.
6. Botero, J. A new mathematical model of the traction force in pneumatic tire snow chain systems / J. Botero, M. Gobbi, G. Mastinu // Associazione Italiana Per L'Analisi Delle Sollecitazioni / Dipartimento di Mec-canica, Politecnico di Milano. 2005. 10 p. - URL: www.aiasonline.org/AIAS2005/Articoli/art084.pdf. (дата обращения 09.09.2018).
7. Белкин, А.Е. Конечно-элементный анализ контакта автомобильной шины с опорной поверхностью на основе оболочечной модели / А.Е. Белкин, Н.Л. Нарская // Вестник МГТУ. Серия Машиностроение. - 2004. - № 3. - С. 14-28.
8. Одинцов, О.А. Разработка метода решения нелинейных контактных задач стационарного качения автомобильной шины: дисс. ... канд. техн. наук / О.А. Одинцова. - М.: МГТУ им. Н.Э. Баумана, 2008.
9. Вольская, Н.С. Математическая модель прямолинейного качения эластичного колеса по неровному деформируемому опорному основанию / Н.С. Вольская, М.М. Жилейкин, А.Ю. Захаров // Известия ВУЗов. Машиностроение. - 2017. - № 10. - С. 25-33.
Ю.Орлов, А.И. Прикладная статистика / А.И. Орлов. - М.: Изд-во «Экзамен», 2006. - 672 с. 11.Жилейкин, М.М. Теоретические основы повышения показателей устойчивости и управляемости колесных машин на базе методов нечеткой логики / М.М. Жилейкин. - М.: Издательство МГТУ им. Н.Э. Баумана, 2016. - 238 с.
Дата поступления в редакцию: 09.10.2018
M.M. Zhileykin, A.J. Zakharov, M.V. Panshin
VERIFICATION OF THE ADEQUACY AND ACCURACY OF THE MATHEMATICAL MODEL OF THE INTERACTION OF THE ELASTIC WHEEL WITH A DEFORMABLE SUPPORT BASE
Moscow state university. N.E. Bauman
Purpose: verification of the accuracy and adequacy of the mathematical model of the elastic wheel interaction with the deformable support base, developed by the authors, by experimental methods.
Design/methodology/approach: verification of the mathematical model was carried out by methods of experimental studies of the real object on the ground.
Findings: it was found that the relative error in determining the various parameters of the interaction of the tire with the deformable support base does not exceed 7.4%.
Research limitations/implications: the carried out researches allow to apply the developed mathematical model of rolling of an elastic wheel on roughnesses of the deformable basic basis for forecasting of basic passability of wheel cars. Originality/value: for the first time, detailed experimental studies of the interaction of an elastic wheel with a deformable support base were carried out.
Key words: a deformable support base, a mathematical model of the wheel, the adequacy and accuracy of the mathematical model.