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

ГОСУДАРС ГВЕННЫЙ НАУЧНЫЙ ЦЕНТР РОССИЙСКОЙ ФЕДЕРАЦИИ

ЦЕНТРАЛЬНЫЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ АВТОМОБИЛЬНЫЙ И АВТОМОТОРНЫЙ ИНСТИТУТ «НАМИ»

На правах рукописи УДК 629.113.001.4

РГБ ОД - 3 ш ш

РЫКОВ СЕРГЕЙ ПЕТРОВИЧ

РАЗРАБОТКА МЕТОДОВ ОЦЕНКИ ПОГЛОЩАЮЩЕЙ И СГЛАЖИВАЮЩЕЙ СПОСОБНОСТИ ПНЕВМАТИЧЕСКИХ ШИН ПРИ РАСЧЕТАХ КОЛЕБАНИЙ АВТОМОБИЛЯ

Специальность 05.05.03 - «Колесные и гусеничные машины»

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

МОСКВА 2000

Работа выполнена на кафедре «Автомобильный транспорт» Братского государственного технического университета

НАУЧНЫЙ РУКОВОДИТЕЛЬ:

доктор технических наук, профессор,заслуженный деятель науки и техники РФ Н.Н.Яценко

ОФИЦИАЛЬНЫЕ ОППОНЕНТЫ: доктор технических наук,

профессор И.В.Балабин

кандидат технических наук, профессор В.В.Таболин

ВЕДУЩЕЕ ПРЕДПРИЯТИЕ - АМО ЗИЛ

Защита диссертации состоится М&& 2000 г. в 14.00 часов нг

заседании диссертационного совета Д 161.01.01 в Государственном научно\ центре Российской Федерации - Центральном научно-исследовательском авто мобильном и автомоторном институте (ГНЦ РФ НАМИ) по адресу: 125438, г Москва, ул.Автомоторная, 2, НАМИ.

С диссертацией можно ознакомиться в научно-технической библиотеке НАМИ.

Отзыв на автореферат (2 экземляра, заверенные печатью) просьба высылать на имя ученого секретаря диссертационного совета по указанному адресу.

Автореферат разослан л* ар/па_2000 г.

Ученый секретарь специализированного совета кандидат технических наук, старший научный сотрудник

Оод-ош.бб-оыцо -V

А.Г.Зубакин

Оы ЧМбоКбО

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

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

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

Объекты исследования: двенадцать шин радиальной и диагональной конструкции легковых автомобилей, микроавтобусов и малотоннажных грузовых автомобилей с посадочным диаметром 14... 16 дюймов отечественного и зарубежного производства, а также грузовой двухосный автомобиль ЗИЛ -4331 и легковой автомобиль высшего класса ЗИЛ — 4102.

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

Научная новизна работы заключается в следующем:

1.Разработана математическая модель движения и теоретический метод определения амплитудно-частотных и спектральных характеристик колебаний динамической системы, эквивалентной подвеске автомобиля, по специальным дорогам автополигона НИЦИАМТа с учетом нелинейной эллиптической модели поглощающей способности шины.

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

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

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

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

б.Экспериментально доказано, что параметры эллиптической модели поглощающей способности шины не зависят от частоты ее радиальной деформации при испытании на невращающемся колесе в диапазоне частот 1...5 Гц, но при более высоком значении параметра пропорциональности модели.

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

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

9.Разработана математическая модель эффекта переменного сглаживания шиной неровностей дороги.

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

11 .Разработана математическая модель движения и численный метод (совместно с программным обеспечением) определения амплитудно-частотных и спектральных характеристик колебаний динамической системы, эквивалентной подвеске автомобиля, и ее звеньев по специальным дорогам автополигона НИ-ЦИАМТа с учетом новой модели эффекта переменного сглаживания шины.

12.Установлены новые закономерности изменения АЧХ звена переменного сглаживания, которые лишены недостатков как модели постоянного сглаживания, так и модели «огибающей», и определяются характеристиками микропрофиля дороги и скоростью движения автомобиля.

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

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

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

Практическая ценность работы.

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

Разработанный шинный стенд и специальное оборудование используются при организации учебного процесса по теории и рабочим процессам автомобиля и курсу НИРС на кафедре автомобильного транспорта БрГТУ.

Реализация работы. Разработанная математическая модель эффекта переменного сглаживания шиной неровностей дороги и ее программное обеспечение внедрены в программную систему FRUND ОГК САПР автозавода AMO ЗИЛ. Акт внедрения прилагается к диссертационной работе.

Апробация работы. Диссертационная работа рассмотрена и одобрена на кафедре «Автомобильный транспорт» Братского государственного технического университета. Основные положения работы докладывались на 7...20-й (1976 ,'.- 1999 гг.) научно-технических конференциях БрГТУ (БрИИ), на Третьей Всесоюзной научной конференции по инерционно-импульсным механизмам, приводам и устройствам (г. Челябинск, 1982 г.), на Первом (г. Москва, 1984 г.), Третьем (г. Москва, 1988 г.) и Четвертом (г. Москва, 1990 г.) Всесоюзных научно-технических совещаниях «Динамика и прочность автомобиля», на Восьмом Международном симпозиуме «Проблемы шин и резинокордных композитов. Дорога, шина, автомобиль» (НИИШП, г. Москва, 1997 г.).

Публикации. По теме диссертаций опубликовано 16 печатных работ, по' лучено одно авторское свидетельство.

Структура и объем работы. Диссертация состоит из введения, анализа состояния проблемы, целей и задач исследования, четырех глав, заключения и выводов, списка использованной литературы и приложений. Работа в целом содержит 318 страниц, 130 рисунков, 12 таблиц, список литературы из 151 наименования.

СОДЕРЖАНИЕ РАБОТЫ

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

Проблеме плавности хода и подрессоривания автомобилей посвящены фундаментальные исследования отечественных ученых И.Г.Пархиловского, Я.М.Певзнера, Р.В.Ротенберга, АЛ.Силаева, А.А.Хачатурова, Н.Н.Яценко и др., в которых разработана современная теория колебаний автомобиля.

В работах И.В.Балабина, Я.М.Горелика, Г.ГГридасова, А.Д. Дербаремди-кера, В.П.Жигарева, А.С.Кольцова, Г.М.Косолапова, Ю.В.Пирковского, А.Е. Плетнева, О.К.Прутчикова, В.М.Семенова, И.С.Степанова, И.Н.Успенского, В.С.Фалькевича, В.С.Щуплякова и др. решены многие важные вопросы действительной оценки рабочих процессов автомобиля в условиях его колебаний от воздействия неровной дороги, реальной нагруженности агрегатов и систем, расчета элементов подвески, ускоренных испытаний автомобильной техники.

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

хода автомобилей, привлекало внимание многих ученых как в России, так и за рубежом. В работах В.Л.Бидермана, Б.Л.Бухина, И.Ф.Бочарова, В.И.Кнороза, К.С. Колесникова, А.С.Литвинова, И.Г.Пархиловского, Я.М.Певзнера, В.А. Петрушова, Р.В.Ротенберга, Я.Е.Фаробина, А.К.Фрумкина, А.А.Хачатурова, Е.А.Чудакова, Н.Н.Яценко, A.Chiesa, M. Julien, E.Marguard, M. Michke, W. Hahn, L.Oberto и др. дан обширный материал по теоретическим и экспериментальным исследованиям процессов, происходящих в пневматической шине при взаимодействии ее с неровной поверхностью дороги и колебаниях радиальной деформации.

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

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

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

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

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

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

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

Для реализации поставленной цели необходимо решить следующие задачи:

1.Обосновать принципы и построить математическую модель колебаний динамической системы, эквивалентной передней (задней) подвеске автомобиля, с учетом нелинейной эллиптической модели H.H. Яценко поглощающей способности шины.

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

2.Обосновать и построить математическую модель корреляционных функций микропрофилей основных испытательных дорог автополигона НИ-ЦИАМТа по экспериментально полученному массиву ординат микропрофилей этих дорог; рассчитать на ЭВМ численные значения коэффициентов модели.

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

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

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

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

Разработать программное обеспечение для персональных ЭВМ численного решения математической модели двухмассовой колебательной системы с учетом эффекта переменного сглаживания и моделирования случайного воздействия, эквивалентного воздействию микропрофиля специальных испытательных дорог автополигона НИЦИАМТа.

• Расчитать и провести графический анализ амплитудно-частотных характеристик звена переменного сглаживания колебательной системы, эквивалент-

ной передней подвеске легкового автомобиля ЗИЛ-4102, при движении с различной скоростью и по различным дорогам автополигона.

5.Разработать программное обеспечение для ввода предлагаемой модели переменного сглаживания шины в программную систему отдела главного конструктора (ОГК) САПР автозавода AMO ЗИЛ для обеспечения формирования дорожного воздействия и получения спектральных характеристик плавности хода и нагруженности проектируемого автомобиля с учетом эффекта переменного сглаживания.

Рассчитать спектральные плотности реакций динамических систем, эквивалентных грузовому автомобилю ЗИЛ-4331 и легковому автомобилю ЗИЛ-4102, в различных точках несущей системы при воздействии микропрофилей дорог автополигона НИЦИАМТа с учетом эффекта переменного сглаживания шин.

Провести ходовые испытания автомобилей ЗИЛ-4331 и ЗИЛ-4102 на специальных; дорогах автополигона НИЦИАМТа, обеспечив запись колебаний в различных точках несущей системы, мостов, элементов подвески и виброзащиты, а также их обработку с целью получения спектральных характеристик плавности хода.

Сравнить результаты ходовых испытаний автомобилей с результатами расчетов по программной системе автозавода AMO ЗИЛ.

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

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

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

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

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

матизированных стендах дает основание аппроксимировать их уравнением эллипса с большим эксцентриситетом:

Р = (1)

где ^ - текущее и амплитудное значения силы неупругого сопротивления в шине;

К, Ь2а - текущее и амплитудное значения радиальной деформации ши-

ны.

а)

20 40 60 80 100 120 б)

Рис. 1. Образцы экспериментальных характеристик упругости шин 245/70 НЯ 16 И-241 «Гранит»: а - в режиме квазистатического нагружения; б - в режиме динамического нагружения.

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

(2)

где коэффициент Нш и показатель степени и являются параметрами модели и характеризуют степень проявления поглощающей способности шин. Дифференциальные уравнения движения масс системы имеют вид: х + И, + 2 кх + &гх =

у

2 к

•»£» (Лг) + <у --х

М

П'

■х = -д,

(3)

где х=г-£ - прогиб упругого элемента подвески; прогиб шины;

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

коэффициент неподрессоренной массы; Яш =НШ /ш - приведенный показатель неупругого сопротивления в шине.

s¡o н

Q -jl jrд

ZI

/6 о/» a=JS,/í h¡ .59 a

Рис. 2. Экспериментальные соотношения между и Ита и их аппроксимация для шины 8.40-15 Я-245.

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

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

Предполагая нелинейность системы (3) в малом, решение ее ищем в виде гармонических функций при входном воздействии гармонического вида q =q0-■sin pt. При этом условии удается получить аналитические выражения частотных характеристик колебательной системы относительно прогиба рессоры и деформации шины в виде трансцендентных уравнений:

Я =

х_

Ч о

= р

Чо

4к2 р10 + (А + а>2В)2 A'+VkpS+BR^1)2

О)

{а + 2kpR Ш/С')2 + [b ~ (П:2 - P2)RmKX Г '

2кр 5 А - (А + со 2 В)(2кр 5 + Яш}г?а-1В)

А(А + со 2 В) + 2кр 5 (2кр

+ Н

п-1

В)

п-\

с» = аг*Я <°2\Ь-№2-р2)Яшк?]-(а +/^„С'ЖХ * [ь - м(П2 - р2 )ЛШЛГ К+с>Ча + ^2крКшкп-х)'

где Л, В, а, б - однотипные составляющие выражений (4) и (5), обособленные для сокращения записи уравнений.

Результаты вычислений амплитудно-частотных характеристик по уравнениям (4) для передней подвески грузового автомобиля, представленные на рис.4 АЧХ по прогибу рессоры, показали, что ввод эллиптической модели по. глощающей способности шины в расчет приводит к заметному снижению амплитуд колебаний, особенно в области резонансных частот. Была доказана также возможность использования только одной кривой из семейства АЧХ для наиболее представительной неровности высотой 2д0 =10 мм; а из двух параметров модели влияние на изменение частотных характеристик коэффициента Нш более значительное, чем показателя степени п.

и* 5мм Чо 10мм 4-Х ММ

\

ч

/

10 20 30 40 50 р,с'

а) б) в)

Рис. 4. Амплитудно-частотные характеристики передней подвески грузового автомобиля по прогабу рессоры: а - учет эллиптической модели; б - влияние коэффициента Нш ; в - влияние высоты неровностей Яо-

Для получения аналитических выражений спектральных плотностей микропрофилей дорог и их воздействия на колебательную систему, были аппроксимированы экспериментальные корреляционные функции трех специальных дорог автополигона НИЦИАМТа выраженями вида:

p(Al) = Axe~a^] + А2е~аЛА1] cos/?, [Д/] + A2e'a'[Al] cos /?2 [дl] (6)

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

Испытательная Дисперсия ор- Параметры аналитической аппроксимации нормированной корреляционной функции

дорога динат

- микропрофиля, D,mm Ai А2 Аз аь -1 M а2, м-1 а3. м-1 Рь м-1 M

Булыжная ровного замо- 258 0,25 0,35 0,4 0,79 0,05 0,29 0,449 0,076

щения Булыжная профилированная 570 0,89 0,37 -0,26 0,59 0,05 0,38 0,766 0,924

"Бельгийская"

мостовая 850 0,75 0,2 0,05 0,8 0,09 0,03 0,247 0,244

Используя преобразование Фурье уравнений (6), полученные выражения частотных характеристик (4) и основное соотношение статистической динамики линейных систем автоматического регулирования - уравнение Винера-Хинчина, были получены аналитические выражения и проведен численный расчет спектральных характеристик колебаний передней подвески грузового автомобиля по прогибу рессоры при движении по специальным дорогам автополигона НИЦИАМТа.

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

ЭДг ммЧ

ш

мы'с

wo

Iii

р

V V л,

V

J 41

О 10 20 30 40 10 60

а)

10 10 30 40 10 60

В)

1 10 10 30 40 10 10 б)

Рис.5. Результаты определения спектральной плотности прогиба рессоры передней подвески автомобиля при испытаниях на "бельгийской" мостовой со скоростью 40 км/ч: а - спектральная плотность воздействия микропрофиля; б - квадрат амплитудно-частотной характеристики подвески; в- спектральная

плотность прогиба рессоры; 1,2 - расчетные спектры без учета и с учетом поглощающей способности шин; 3 - экспериментальный спектр.

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

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

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

Математическое описание условий разработки метода позволило вывести выражение коэффициента эквивалентного вязкого трения вида:

1ш\ zc т р

где Ншип- параметры эллиптической модели; hzcm - статический прогиб шины, соответствующий уровню радиальной нагрузки на колесо и внутреннему давлению воздуха в шине.

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

Во второй главе обоснована и разработана модель эффекта переменного сглаживания шины; разработан математический метод учета новой модели при расчетах колебаний эквивалентных систем автомобиля; разработан автономный программный комплекс спектрального анализа колебаний двухмассовых систем, эквивалентных подвеске автомобиля, с учетом новой модели; проведен численный расчет и графический анализ АЧХ звена переменного сглаживания; разработано программное обеспечение ввода эффекта переменного сглаживания шины в систему FRUND AMO ЗИЛ; выполнены численные расчеты энергетических спектров виброускорений грузового и легкового автомобилей ЗИЛ.

Математическое описание новой модели переменного сглаживания шины базируется на следующих допущениях: 1 - пневматическая шина обладает

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

Тогда основное уравнение модели переменного сглаживания шины может быть выражено оператором скользящего среднего с переменным интервалом осреднения, т.е.

где а (1) - переменная длина пятна контакта; q (1) - функция исходного микропрофиля; qcr (1) - функция сглаженного микропрофиля.

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

Математический метод учета модели переменного сглаживания шины при расчетах колебаний разработан на примере двухмассовой динамической системы, эквивалентной передней (задней) подвеске автомобиля.

Интегро-дифференциальные уравнения движения масс и функционирования звеньев системы имеют вид

| /+а(/)/2

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

Мг + т) (г - #")+ 2ср (г - £) = О,

(4 ~ 4 ) + 2 с (4 ~ Я )~П (г - #)- 2с (г - £) = О, ь 'ш 4 сг ' 4 сг р р

I + а {1)1 2

а{1) =

В общем случае система уравнений (9) не линейна и не стационарна, и не может быть решена методами операционного исчисления. Поэтому был разработан программный комплекс КМРБ для численного решения уравнений системы и спектрального анализа колебаний подвески. Блок-схема программного комплекса приведена на рис.7.

Рис.7

Программный комплекс разрабатывался на следующих принципах: во-первых, вычисления АЧХ звеньев и колебательной системы в целом строились на методах спектрального анализа линейных САР; во-вторых, вычислялись не спектральные плотности входных воздействий и реакций системы, а только спектры амплитуд с использованием алгоритма быстрого преобразования Фурье; в-третьих, программы конструировались так, чтобы можно было выделить

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

При помощи программного комплекса были рассчитаны амплитудно-частотные характеристики звена переменного сглаживания для колебательной системы, эквивалентной передней подвеске легкового автомобиля ЗИЛ — 4102, при движении по двум специальным дорогам автополигона (рис.8).

Анализ графиков показал, что новая модель переменного сглаживания шиной неровностей дороги проявляется в существенном изменении АЧХ звена переменного сглаживания, характер которой определяется параметром осреднения а(1) и зависит от скорости движения автомобиля и дисперсии ординат микропрофиля. Изменение АЧХ звена переменного сглаживания происходит от одного предельного состояния при va = 0, когда она описывается по закону характеристики звена постоянного сглаживания (кривая 1), к другому состоянию, близкому к отрыву шин от опорной поверхности дороги, когда она описывается уравнением огибающей.

ч

\ 2 Л Ч Г~

и-з, Г' /

/

02- Or- lo 1« «0 13 loo U0

роЗ*

а) б)

Рис.8. АЧХ звена переменного сглаживания: а - дорога с асфальтобетонным покрытием; б - дорога с ровным булыжным покрытием; 1- АЧХ звена постоянного сглаживания; 2,3,4 - АЧХ звена переменного сглаживания при различных скоростях движения (У4 > У3 > Уг).

Реализация новой модели эффекта сглаживания шиной неровностей дороги в программной системе FRUND (формирование и решение уравнений нелинейной динамики) автозавода AMO ЗИЛ позволило рассчитать оценочные показатели и характеристики плавности хода (средние квадратические отклонения и спектральные плотности) грузового автомобиля ЗИЛ - 4331 и сравнить их со спектрами по результатам натурных испытаний (рис.9).

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

(К та

Е_1_

I ¡Г I_I I ¡У^-

а) б) Рйс.9. Спектральные плотности вертикальных ускорений на левом лонжероне рамы при движении автомобиля ЗИЛ - 4331 по ровной булыжной дороге со скоростью 30 км/ч с грузом 6т: а - над передним мостом; б - над задним мостом; — эксперимент; - расчет

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

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

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

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

Механизм вынужденных колебаний реализует режим резонансных колебаний динамической системы шинного стенда на испытуемой шине при помо-

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

Механизм передвижения платформы реализует свободное качение колеса для исследования взаимодействия пневматической шины с неровностями гармонического профиля за счет продольного возвратно-поступательного движения платформы с установленными на ней неровностями.

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

В состав измерительно-информационной системы стенда входят пять комплексов: комплекс средств преобразования (датчики, усилители и измерительные преобразователи), комплекс средств измерения и регистрации (промышленные образцы измерительных и регистрирующих приборов аналогового и цифрового типов), комплекс средств автоматизации и управления (САР "Цикл" и САР "Частота"), комплекс средств питания (источники питания регулируемого напряжения постоянного и переменного тока) и комплекс средств коммутации (пульт связи цифропечатающего устройства, пульт коммутации датчиков). В процессе реализации ИИС стенда были разработаны, изготовлены и внедрены: одно- и двухканальный тензоусилители постоянного тока, измеритель частоты вращения вала виброустановки, измеритель температуры покрышки шины, трехканальный регулируемый источник питания постоянного тока, САР циклом "нагрузка-разгрузка" и САР частоты вращения роликового механизма. Это позволило обеспечить необходимые режимы нагружения испытуемой шины, точность задания и регулирования в необходимых диапазонах их параметров.

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

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

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

шин 7-ми типоразмеров на шинном стенде; представлены задачи, программа, оборудование, оценка, погрешности и результаты ходовых испытаний грузового и легкового автомобилей ЗИЛ на специальных дорогах автополигона НИЦИ-АМТа.

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

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

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

Влияние степени износа протектора, материала и числа слоев нитей корда и других конструктивных особенностей шины определяется изменением только коэффициента Нш, показатель степени п остается неизменным. Численные значения показателя п для всех испытанных шин находятся в пределах 0,43 < п <0,60 , т.е. между предельными случаями, когда п=0 (модель «сухого» трения) и п=1 (модель "вязкого" трения).

Численные значения параметров эллиптической модели, коэффициента радиальной жесткости и их статистическая оценка для испытанных шин приведены в табл. 2.

Таблица 2.

Значения параметров поглощающих и упругих свойств шин в режиме квазистатического нагружения

Размеры и модель Ш1ШЫ Параметры эксплуатационного состояния Исследуемые параметры

Pzh> кН Pwh, МПа Нц], Н/мм" кН/мп ОНш, кН/мп п cz, кН/м СТс. кН/м

1 2 3 4 5 6 7 8 9

245/70 HR16 И - 241 9,0 0,24 78.6 2,32 0,04 0,49 0,02 332 9,8

245/70 HR16 И-241 (износ 57%) 9,0 0,24 58.6 1,73 0,03 0,49 0,02 320 9,2

LR78-15 «SUPER ARA MID RADIAL» 9,0 0,24 58.3 1,14 0,02 0,43 0,02 263 7,8

LR70-15 «GOLDEN LIFSAVER» 9,0 0,24 47.8 1,23 0,02 0,47 0,02 262 7,7

185 R 15 Я-288 6,35 0,32 49.5 3,12 0,05 0,6 0,03 . 341 10,2

7.00-15 И-89 6,15 0,22 64.2 3,07 0,05 0,56 0,02 266 7,9

8.40-15 Я-245 5,95 0,20 88.7 5,22 0,08 0,59 0,03 338 10,1

7.35-14 ИД-195 4,35 0,17 72.5 2,0 0,03 0,48 0,02 216 6,4

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

Таблица 3.

Параметры оценки поглощающей способности шин.

Параметр оценки■ Размеры испытуемых шин

73-14 ИД-195 7.00-15 И-89 8.40-55 Я-245 1851*15 Я-288 ЬЯ70- 1Д78-15 БАЯ 245/70 НИ 16 И-241 24570 НЯ16 И-241 и» юс 57%

1 2 3 4 5 6 7 8 9

Ьет, мм 23,2 26,6 19,2 22,2 38,5 38,1 31,1 31,5

Дал, Дж 23,9 33,7 21,6 22,2 32,2 33,4 30,8 42,2

Пхар, Дж 25,2 35,6 26,8 23,1 33,4 34,6 31,8 43,7

кл, кН/м 19,1 19,4 28,9 17,9 9,4 10,2 18,3 13,5

Примечание: Пэл - площадь аппроксимирующего эллипса;

Пхар-плошадь экспериментальной характеристики упругости.

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

Рис.10. Траектории оси колеса при обкатывании шиной гармонической неровности: 1 - профиль неровности; 2 - экспериментальная траектория; 3 -расчетная траектория.

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

Основная задача ходовых испытаний грузового ЗИЛ-4331 и легкового ЗИЛ-4102 автомобилей на специальных дорогах автополигона НИЦИАМТа состояла в получении реализаций ускорений и относительных перемещений в различных точках автомобилей с дальнейшей приборной обработкой их для получения спектральных характеристик колебаний (рис.11).

¿п

+-

,А г—V— 4

! / \ т

/ / 1 и

3 1—

1— ь

1

Л1?

III!

■ , а

а)

б)

Рис.11 .Спектральные плотности вертикальных виброускорений легкового автомобиля ЗИЛ-4102 при движении по ровной булыжной дороге со скоростью 30 км/ч: а - на правом переднем сидении; б - под сидением на кузове

Сопоставление расчетных и экспериментальных спектральных характеристик колебаний легкового автомобиля как по уровню, так и частотам проявления резонансных зон, показало достаточно хорошую их сходимость, особенно в низкочастотной области. Менее хорошее совпадение энергетических спектров наблюдается в области частот 22,4...45 Гц, где проявляются колебания несущей системы автомобиля; это, видимо, связано с недостаточной проработкой на момент счета параметров динамической модели автомобиля. Погрешность несовпадения кривых спектральных плотностей вертикальных виброускорений составляет 10.. .20%, что вполне приемлемо для статистических вычислений.

А

£

ЗАКЛЮЧЕНИЕ И ОСНОВНЫЕ ВЫВОДЫ

Не смотря на сложность и трудоемкость поставленных в данной диссертационной работе задач, все они были успешно решены. А именно:

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

2.Получены аналитические выражения энергетических спектров микропрофиля трех специальных дорог автополигона НИЦИАМТа и их воздействия путем аппроксимации экспериментальных корреляционных функций и применения преобразования Фурье, а также энергетических спектров реакций эквивалентной подвеске системы с использованием уравнения Винера-Хинчина.

3.На примере передней подвески грузового автомобиля ЗИЛ-131 и с использованием разработанной методики учета эллиптической модели были рассчитаны спектральные характеристики колебаний автомобиля по прогибу рессоры при движении по специальным дорогам автополигона; дан сравнительный анализ расчетных и экспериментальных спектров.

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

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

6.Разработан алгоритм и программный комплекс RMPS для малых ЭВМ численного решения интегро-дифференциальных уравнений движения двух-массовых колебательных систем, эквивалентных подвеске автомобиля, с учетом новой модели переменного сглаживания шины.

7.На примере передней подвески автомобиля ЗИЛ-4102 выполнены расчеты и построены амплитудно-частотные характеристики звена переменного сглаживания; дан графический анализ характеристик.

8.Разработано программное обеспечение для ввода новой модели сглаживающей способности шины в программную систему FRUND автозавода AMO ЗИЛ; выполнены расчеты спектральных характеристик колебаний пространственных моделей, эквивалентных грузовому ЗИЛ-4331 и легковому ЗИЛ-4102 автомобилям, с учетом разработанной модели.

9.Получены в результате ходовых испытаний грузового ЗИЛ-4331 и легкового ЗИЛ-4102 автомобилей на специальных дорогах автополигона НИЦИ-АМТа экспериментальные значения энергетических спектров виброускорений и относительных перемещений в различных точках автомобилей; дан графический анализ расчетных и экспериментальных спектральных характеристик испытанных автомобилей.

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

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

На основании проведенных теоретических и экспериментальных исследований можно сделать следующие выводы:

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

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

3.Новая модель переменного сглаживания шиной неровностей дороги проявляется в существенном изменении амплитудно-частотной характеристики звена сглаживания, характер которой определяется параметром осреднения -переменной длиной пятна контакта и зависит от скорости движения автомобиля и дисперсии координат микропрофиля. Изменения АЧХ звена переменного сглаживания происходит от одного предельного состояния при Va =0, когда она описывается по закону характеристики звена постоянного сглаживания, к другому состоянию близкому к отрыву шин от опорной поверхности дороги, когда она описывается уравнением огибающей.

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

5.Разработанный программный комплекс RMPS позволяет проводить численное решение интегро-дифференциальных уравнений модели переменного сглаживания и спектральный анализ систем подрессоривания, эквивалентных подвеске автомобилей, в автономном режиме, т.е. обеспечивая моделирование микропрофилей дорог автополигона НИЦИАМТа, расчет и графическое построение частотных и спектральных характеристик звеньев и системы в целом.

6.Разработанное программное обеспечение для системы FRUND ОГК САПР автозавода AMO ЗИЛ позволило ввести новую модель переменного сглаживания шиной неровностей дороги в процесс формирования возмущаю-

щего воздействия микропрофиля и моделирования динамики движения и колебаний сложных систем, эквивалентных "объемному" автомобилю с упругими рамами и кабинами (кузовами).

7.Проведенные ходовые испытания грузового и легкового автомобилей ЗИЛ на специальных дорогах автополигона НИЦИАМТа позволили получить оценочные параметры и характеристики колебаний в виде диаграмм энергетических спектров и таблиц средних квадратических отклонений вертикальных ускорений и относительных перемещений в различных точках автомобилей.

8.Сопоставление расчетных и экспериментальных характеристик колебаний автомобилей на специальных дорогах автополигона НИЦИАМТа показало достаточно хорошую сходимость энергетических спектров вертикальных ускорений для грузового (ЗИЛ - 4331) и легкового (ЗИЛ - 4102) автомобилей, как по уровню, так и по частотам проявления резонансных зон, особенно в низкочастотной области колебаний; с погрешностью в пределах 10...20%.

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

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

Основное содержание диссертации опубликовано в следующих работах:

1.А.С. 1515077 (СССР) Устройство для измерения сил на колесе транспортного средства / Братский индустр. ин-т; Авт изобр. С.П.Рыков и ВГ.Мелентьев.- Заяв. 04.01.88. № 4382096.

m

27

2.Карцов C.K., Плетнев А.Е., Рыков С.П., Раввин А.Г., Ямпольский Д.А. Учет эффекта переменного сглаживания пневматической шины при формировании случайного воздействия дорог // IV Всесоюзное научно-техническое совещание "Динамика и прочность автомобиля": Тезисы докл.- Москва, 1990.-С. 106-107.

3.Мелентьев В.Г., Рыков С.П., Яценко H.H., Степанов И.С. Изменение динамики подрессоривания автомобиля при переменном сглаживании воздействия дороги // III Всесоюзное научно-техническое совещание "Динамика и прочность автомобиля": Тезисы докл. -Москва, 1988.- С. 169-170.

4.Рыков С.П. Измерение сил при стендовых испытаниях автомобильных шин / Ред. ж. "Автомобильная промышленность".- М., 1990.-10 е.: ил.-Библиогор. 12 назв.- Рус.- Деп. в ЦНИИТЭИ автопроме, № 2022- ап.

5.Рыков С.П., Степанов П.И., Яценко H.H.. Некоторые результаты экспериментальных исследований поглощающей способности ряда шин легковых автомобилей // I Всесоюзное научно-техническое совещание "Динамика и прочность автомобиля": Тезисы докл. - Москва, 1984,- С. 93-94.

6.Рыков С.П., Яценко H.H. К вопросу о сглаживающей способности эластичной шины автомобильного колеса / Братский индустр. ин-т.- Братск, 1989.-5с.: ил.- Библиогр. 4 назв.- Рус.- Деп. в ЦНИИТЭ автопроме, № 1877- ап.

7.Рыков С.П., Ламаков C.B. Стенд для исследований колесных движителей лесных машин // XIII Научно-техническая конференция БрИИ: Тезисы докл.- Братск, 1992,- С. 58.

8.Рыков С.П., Савранский H.H. Тензометрическая площадка для исследований упругих, демпфирующих и сцепных свойств автомобильных шин // XIV Научно-техническая конференция БрИИ: Тезисы докл.- Братск, 1993.- С. 73.

9.Рыков С.П., Коновницын И.П., Пак Д.Л. Результаты стендовых испытаний пневматических шин легковых автомобилей высшего класса // XVIII Научно-техническая конференция БрИИ: Тезисы докл.- Братск, 1998.- С. 88-89.

Ю.Рыков С.П. Тензометрический датчик перемещений и низкочастотных колебаний // XIX Научно-техническая конференция БрИИ: Тезисы докл,-Братск, 1998.- С. 268-269.

П.Рыков С.П., Яценко H.H. Экспериментальные исследования поглощающей способности пневматических шин / Труды Братского госуд. индустр. ин-та: Материалы XX Научно-технической конференции. - Братск, 1999.-Т.2,-С. 131-134.

12.Рыков С.П., Ермаков С.М. Инженерный метод учета поглощающей способности пневматической шины в системах подрессоривания автомобиля/ Труды Братского госуд. индустр. ин-та: Материалы XX Научно-техническая конференции.-Братск, 1999.-Т.2,-С. 134-136.

13.Яценко H.H., Капанадзе Т.Н., Рыков С.П. Колебания подвески с учетом поглощающей способности шин // Автомобильная промышленность 1977.-№7,-С. 15-18.

14.Яценко H.H., Капанадзе Т.Н., Рыков С.П. и др. Колебания подвески с учетом поглощающей способности шин при случайном возмущении // Автомобильная промышленность 1979.- №1.- С. 16-19.

15-Яценко H.H., Рыков С.П., Карцов С.К., Плетнев А.Е., Раввин А.Г. Новая модель сглаживающей способности шин: Расчет колебаний автомобиля // Автомобильная промышленность.- 1992.-№11,- С. 18-21.

16Лценко H.H., Рыков С.П. Моделирование поглощающей способности пневматической шины в системах подрессоривания автомобиля / Труды Братского госуд. индустр. ин-та: Материалы XX Научно-технической конференции.-Братск, 1999.-Т.2.- С. 128-131.

Рыков Сергей Петрович

РАЗРАБОТКА МЕТОДОВ ОЦЕНКИ ПОГЛОЩАЮЩЕЙ И СГЛАЖИВАЮЩЕЙ СПОСОБНОСТИ ПНЕВМАТИЧЕСКИХ ШИН ПРИ РАСЧЕТАХ КОЛЕБАНИЙ АВТОМОБИЛЯ

Специальность: 05.05.03 - Колесные и гусеничные машины

Подписано в печать 8.02.2000. Формат 60 х 84 '/|6 Печать трафаретная Уч.-изд. л. 1,75. Усл. печ. л. 1,6. Тираж 100 экз. Заказ 35.

Отпечатано в издательстве БрГТУ 665709, Братск, ул. Макаренко, 40

ВВЕДЕНИЕ, АНАЛИЗ СОСТОЯНИЯ ПРОБЛЕМЫ, ЦЕЛИ И ЗАДАЧИ

ИССЛЕДОВАНИЯ.

ГЛАВА 1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ КОЛЕБАНИЙ АВТОМОБИЛЯ

С УЧЕТОМ ГОГЛОЩАЮЩЕЙ СПОСОБНОСТИ ШИН.

1.1. Эллиптическая модель поглощающей способности пневматических шин.

1.2. Эквивалентная колебательная система и дифференциальные уравнения движения ее масс с учетом предлагаемой модели поглощающей способности шин.

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

1.4. Построение модели случайного возмущающего воздействия микропрофиля дороги.

1.5. Расчет и анализ колебаний автомобиля с учетом поглощающей способности шин на неровных дорогах.

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

ГЛАВА 2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ КОЛЕБАНИЙ АВТОМОБИЛЯ

С УЧЕТОМ СГЛАЖИВАЮЩЕЙ СПОСОБНОСТИ ШИН.

2.1. Основные допущения для составления модели сглаживающей способности пневматических шин. -.-.

2.2. Эквивалентная колебательная система и дифференциальные уравнения движения ее масс с учетом предлагаемой модели сглаживающей способности шин.

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

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

ГЛАВА 3. РАЗРАБОТКА, ИЗГОТОВЛЕНИЕ И МЕТОДИКА

ИСПОЛЬЗОВАНИЯ НОВОГО ОБОРУДОВАНИЯ ДЛЯ ЭКСПЕРИМЕНТАЛЬНЫХ ИССЛЕДОВАНИЙ ПНЕВМАТИЧЕСКИХ ШИН КАК ЭЛЕМЕНТОВ СИСТЕМ ПОДРЕССОРИВ АНИЯ АВТОМОБИЛЕЙ.

3.1. Универсальный шинный стенд. Назначение и основные требования для проектирования.

3.2. Состав стенда и конструкция его элементов.

3.3. Теоретические основы разработки и конструкции датчиков шинного стенда.

3.4. Измерительно-информационная система шинного стенда.

3.5. Методика статических испытаний шин на стенде. Оценка погрешности экспериментальных данных.

3.6. Теоретические основы и методика испытания шин в динамическом режиме нагружения на стенде.

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

ГЛАВА 4. РЕЗУЛЬТАТЫ ЭКСПЕРИМЕНТАЛЬНЫХ ИССЛЕДОВАНИЙ ПАРАМЕТРОВ ЭКВИВАЛЕНТНЫХ КОЛЕБАТЕЛЬНЫХ СИСТЕМ И КОЛЕБАНИЙ АВТОМОБИЛЯ С УЧЕТОМ ПОГЛОЩАЮЩЕЙ И СГЛАЖИВАЮЩЕЙ СПОСОБНОСТИ

4.1. Испытания и оценка упругих и поглощающих свойств пневматических шин на стенде.

4.2. Испытания и оценка сглаживающих свойств пневматических шин на стенде.

4.3. Исследования колебаний при испытаниях автомобилей на специальных дорогах.

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

Среди основных эксплуатационных свойств автомобиля плавность хода, выражаемая характеристиками колебательных процессов его масс, занимает особое место. Это объясняется существенным влиянием колебаний кузова и колес, возникающих при движении по неровностям дороги, почти на все эксплуатационные качества автомобиля [98]. Вместе с тем, характерной особенностью общей дорожной сети России является большая протяженность дорог с переходным покрытием и грунтовых. В этих условиях особенно велики прямые и косвенные потери на автомобильном транспорте за счет недостаточной защищенности автомобиля от динамических воздействий неровной дороги. Фактические данные показывают, например, увеличение расхода топлива почти в два раза, увеличение себестоимости перевозок 1 т.км более чем в два раза по грунтовым и переходным дорогам по сравнению с перевозками по дорогам с усовершенствованным покрытием.

Проблеме плавности хода и подрессоривания автомобилей посвящены фундаментальные исследования отечественных ученых И.Г.Пархиловского, Я.М.Певзнера, Р.В.Ротенберга, А.А.Силаева, А.А.Хачатурова, Н.Н.Яценко и др., в которых разработана современная теория колебаний автомобиля. Эффективное приложение этой теории развивалось затем как в направлении совершенствования конструирования подвески, так и в новом направлении совершенствования расчетов и оценки нагруженности основных агрегатов трансмиссии и несущей системы в реальных условиях движения, а также построения теории и методов форсированных полигонных испытаний.

В работах И.В.Балабина, Я.М.Горелика, Г.Г.Гридасова, А.Д.Дербаремди-кера, В.П.Жигарева, А.С.Кольцова, Г.М.Косолапова, Ю.В.Пирковского, А.Е.Плетнева, О.К.Прутчикова, В.М.Семенова, И.Н.Успенского, В.П.Шалды-кина и др. решены многие важные вопросы действительной оценки рабочих процессов автомобиля в условиях его колебаний от воздействия неровной дороги, реальной нагруженности агрегатов и систем, расчета элементов подвески, ускоренных испытаний автомобильной техники. Благодаря этому теория автомобиля как научная дисциплина была значительно продвинута вперед. Установление тесной связи переменных воздействий неровной дороги и возникающих колебаний масс автомобиля значительно приблизило методы проектирования систем подрессоривания и методы испытаний новых моделей автомобилей к условиям реальной эксплуатации.

В настоящее время проектные расчеты подвески строятся исходя из рассмотрения колебаний автомобиля как составной части системы дорога - шина -автомобиль - водитель (ДШАВ). Новая технология форсированных полигонных испытаний создается на основе оценки динамических нагрузок основных базовых агрегатов автомобиля при движении по специальным неровным дорогам в процессе его колебаний.

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

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

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

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

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

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

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

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

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

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

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

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

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

Перспективность исследований вопросов внешней механики пневматической шины и, в том числе, ее преобразующих свойств, влияющих на плавность хода автомобиля, привлекало внимание многих ученых как в России, так и за рубежом. В работах В.Л.Бидермана, Б.Л.Бухина, Н.Ф.Бочарова, В.И.Кнороза, К.С.Колесникова, A.C.Литвинова, И.Г.Пархиловского, Я.М.Певзнера, В.А.Пет-рушова, Р.В.Ротенберга, Я.Е.Фаробина, А.К.Фрумкина, А.А.Хачатурова, Е.А.Чудакова, Н.Н.Яценко, A.Chiesa, M.Julien, E.Marguard, M.Michke, W.Hahn, L.Oberto и др. дан обширный материал по теоретическим и экспериментальным исследованиям процессов, происходящих в пневматической шине при взаимодействии ее с опорной поверхностью дороги и колебаниях радиальной деформации.

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

Одной из первых публикаций, посвященных исследованию поглощающей (демпфирующей) способности пневматических шин, была статья К.С.Колесникова [54], в которой приведены методика и результаты экспериментальных исследований внутренних потерь в автомобильной шине 7.50-17 при ее радиальном деформировании. Эксперименты проводились на стенде рычажного типа (рис.В.1), схема которого неоднократно использовалась для создания шинных стендов в последующих исследованиях шин другими учеными.

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

Рис. В. 1. Схемы стенда для исследований упругих и демпфирующих свойств шин: а - в режиме квазистатического нагружения; б - в режиме свободных колебаний.

В качестве основного метода оценки демпфирующих свойств шины в работе [54] был принят метод обработки экспериментальных кривых затухающих колебаний в предположении об убывании размахов колебаний по закону геометрической прогрессии, что означает изменение силы неупругого сопротивления в шине по закону вязкого трения, т.е. пропорционально скорости вертикальных колебаний шины. Теоретическое обоснование методик обработки кривых затухающих колебаний для различных законов изменения силы неупругого сопротивления шины приведено в монографии К.С.Колесникова [55]. В случае сопротивления шины, пропорционального скорости деформации шины, оно оценивается коэффициентом сопротивления, величина которого для приведенной на рис.В.1, б схемы стенда определяется по выражению У кв =

2-1 Лп 5 Ь2-Т в.1) где кв - коэффициент сопротивления шины; - момент инерции стенда относительно оси качания;

1пд - логарифмический декремент затухания; Ъ - расстояние от оси колеса до оси качания; Т - период колебаний стенда на шине.-В результате обработки экспериментальных кривых затухающих колебаний в предположении о вязкой природе сил неупругого сопротивления в шине были получены значения коэффициента сопротивления в зависимости от частоты собственных колебаний стенда и внутреннего давления ^ воздуха в шине (рис.В.2).

70 N кол/сек

Рис.В.2. Зависимость коэффициента сопротивления шины 7.50-17 от частоты собственных колебаний стенда (р- давление воздуха в шине, Сгк - нагрузка на колесо)

Анализ кривых показал, что:

1) коэффициент сопротивления шины с ростом частоты колебаний стенда уменьшается, но в меньшей степени, чем увеличивается частота колебаний;

2) коэффициент сопротивления шины с ростом внутреннего давления в шине уменьшается примерно в одинаковой степени в диапазоне изменения частот от 2 до 8 Гц;

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

4) сила неупругого сопротивления, определенная по упругой характеристике шины, не равна силе сопротивления, определенной по кривым затухающих колебаний;

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

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

Так в работе [3], выполненной под руководством Н.Ф.Бочарова, приводятся результаты экспериментальных исследований демпфирующих свойств пневмокатка модели И-245 при испытании его на автомобиле ЗИЛ-1Э2С методом сбрасывания с заблокированной подвеской. Было принято, что сила трения в пневмокатке пропорциональна первой степени скорости его радиальной деформации, а огибающие кривых затухающих колебаний автомобиля на испытуемой шине описывались уравнением экспоненты

А = А0-е~п\ (В.2) где А - текущее значение амплитуды вертикальных колебаний;

А о - высота сбрасывания; п - коэффициент затухания; время.

Обработка экспериментальных кривых по методике, аналогичной методике [55], позволила определить логарифмический декремент затухания и коэффициент апериодичности колебаний в зависимости от внутреннего давления воздуха в пневмокатке.

В работе [66], выполненной в Тольяттинском политехническом институте, исследовались упругие и демпфирующие характеристики шин легковых автомобилей на стенде, схема которого аналогична схеме стенда [54], но оснащенного более совершенным оборудованием: динамометрической площадкой, реохордным датчиком перемещений, светолучевым осциллографом Н-700. Исходным экспериментальным материалом являлись кривые затухающих колебаний масс стенда на испытуемой шине при изменении внутреннего давления воздуха, нагрузки и скорости радиальной деформации для шин 165-81113 РМу, 165-1ЗР ИЯ-170, 165-13 И-151.

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

2-1-у 1-(1Я2+У2) ги ~ 2 гг ' С™ ~ 2 гр2 ' а -Т а -Т где: / - суммарный момент инерции колеса, грузов и балки стенда относительно оси качания; V - логарифмический декремент затухания; а - расстояние от оси колеса до оси качания стенда; Т- период колебаний стенда.

Результаты обработки кривых затухающих колебаний приведены на рис. В.З иВ.4.

Сщ.кгс/си

ZOO

160

ПО

30

-

У

X / А х-—- —- ~— .J ——•--" -X .—•— 2 -X— 3

-

Сш, пес/см

150 гоо

250 500 й„,пгс гоо

160 по

80

-X-X —г —X— 3

У / \ г : \ \ \ \ лЧ V

V 1—

15 го

25 WjC' а) б)

Рис. В.З. Зависимость коэффициента нормальной жесткости: а - от нагрузки и внутреннего давления воздуха; б - от частоты собственных колебаний; 1 - шина 165-81113; 2 - шина 165-13Р;

3 - шина 165-13 (новая); 4 - шина 165-13 (износ 100%); 2 2 р^-Х^Ъ кгс/см ;.рм> = 2,0 кгс/см . кг с-с/си

V ОМ Ц06 пп,

100 150 WO 250 &к,пгс ' 0 1? 20 25 СО, с'7 а) ö)

Рис. В. 4. Зависимость параметров демпфирования шин: а - коэффициента апериодичности от нагрузки; б - коэффициента сопротивления от частоты собственных колебаний стенда (обозначение шин те же, что и на рис. В.З). Анализ графиков показал, что:

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

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

В работе [64], выполненной под руководством А.С.Литвинова, приведены результаты комплексных исследований большегрузных шин: 18.00-24 И-170, 320-508 ИНВ-12Б, 1500 х 600 х 635 В-77, в т.ч. упругих и демпфирующих свойств при действии радиальной нагрузки. Схема стенда и методика обработ

X— 3 -Tür- ¡f о On N ОО NP к >д

• »Ni ^ X» Л^х О ки экспериментальных кривых затухающих колебаний аналогичны рассмотренным в [54].

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

Схожие закономерности коэффициентов сопротивления и жесткости от частоты свободных колебаний при описании сил неупругого сопротивления в испытуемых шинах вязкой моделью трения были получены и в других исследованиях, например, [75, 123, 124, 150].

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

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

Например, в публикации В.К.Вахламова [21] даны результаты исследования динамического трения в шинах автомобиля Москвич-408, которые проводились на специальном стенде платформенного типа, реализующего динамическое приложение радиальной нагрузки к невращающемуся колесу. Измерение сил, действующих на шину, производилось при помощи тензометрической консоли, а деформации - реохордным датчиком перемещений; запись изменения сил и деформаций во времени осуществлялась при помощи светолучевого осциллографа Н-700 на фоточувствительную ленту. Осциллограммы изменения усилий и деформаций при варьировании частоты нагружения шин от 60 до 600 кол/мин. обрабатывались путем измерения ординат кривых в точках, соответствующих 1/4 и 3/4 периода колебаний. По мнению автора публикации половина разности значений ординат в этих точках (а, б на рис.В.5) в масштабе сил соответствует максимальному значению динамической силы трения в шине.

-Н 0.1 сек п*60нал/мин п-Школ/мин и

Рис.В.5. Осциллограммы изменения усилий при динамической деформации шин автомобиля Москвич-408 (величина деформации Н = 20 мм)

Анализ результатов обработки указанным образом осциллограмм показал, что сила неупругого сопротивления в шине возрастает в-2-.З раза с ростом частоты приложения нагрузки от 60 до 600 кол/мин. Делается вывод о том, что трение в шине при ее работе по своей природе относится к межмолекулярному трению, которое изменяется пропорционально скорости колебаний.

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

В ряде публикаций делается попытка смоделировать демпфирующие свойства пневматической шины при ее вертикальных колебаниях исходя из основных аксиом реологии, полагающих, что всем телам присущи свойства фундаментальных тел: Гука, Ньютона и Сен-Венана, выраженные в разной степени [27]. Так, в работе В.Н.Белковского [9] рассматривается обобщенная модель Кельвина-Фойгта, включающая параллельно соединенные упругий с? и вязкий к элементы, отражающие соответственно высокоэластичные и пластические свойства материалов шины, а также последовательно соединенный с ними упругий элемент сотражающий упругие свойства сжатого в шине воздуха. Сила сопротивления деформированию шины согласно данной трехэлементной модели при вынужденных гармонических колебаниях равна Г

Рш к0-Т

1+ х + с + к0Т2со2 Л 1 + Т2со2 ОС)

В.4) где

V - ■ - — у СО - собственная частота колебаний колеса; Со, Т- коэффициенты, определяемые по выражениям

С = с\'с2 с\+с2

К С 1 т = к сх+с2 сх+с2

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

Все ранее обобщенные работы моделировали упругодемпфирующие свойства пневматических шин в рамках, так называемой, внешней механики, т.е. рассматривали шину как нечто целое, наделенное комплексом выходных характеристик, определяющих ее взаимодействие с другими элементами системы автомобиля. При этом методы исследования базируются на принципах аналитической механики систем с неголономными связями. В то время, как внутренняя механика пневматических шин базируется на принципах механики деформируемого твердого тела и ставит своей задачей определение прочности и долговечности шин в различных режимах работы. Полного решения задачи определения жесткости и демпфирования колеблющейся и одновременно катящейся шины в рамках внутренней механики шин даже при существенных упрощениях до настоящего времени нет. Однако, некоторые модели внутренней механики шин успешно используются в сочетании с моделями внешней механики для построения расчетов упругих и демпфирующих свойств шин в режиме качения колеса. Так, например, в работе [57], выполненной в МВТУ под руководством Н.Ф.Бочарова, для получения зависимости коэффициента демпфирования от скорости качения колеса используется модель шины типа "кольцо на упругом основании". Полная деформация шины при колебаниях и качении колеса разбивается на две составляющие: изгиб беговой части и радиальное сжатие элементов протектора, а коэффициент демпфирования равен сумме

К= кс+крки, (В.5) где кс, ки - составляющие демпфирования при сжатии и изгибе; кф - скоростной коэффициент изгибной деформации, выражение которого выводится из рассмотрения кинематики модели "кольцо".

Основываясь на экспериментальных данных, приводимых в отечественных и зарубежных публикациях [126, 148], было исследовано изменение коэффициентов к9 и кш от скорости качения колеса в диапазоне до 85 км/ч, которые показали, что с ростом скорости качения коэффициент к^ быстро падает до значения 0,1 .0,2 и далее асимптотически приближается к нулю, следовательно, коэффициент кш при этом асимптотически стремится к кс, т.е. к постоянному значению, определяемому сжатием элементов протектора шины. Существенным недостатком предлагаемой модели демпфирования в шине является априорное признание вязкой природы трения для обеих составляющих деформации шины.

В статье Б.И.Морозова [75], посвященной исследованию рабочего процесса пневматической шины, на основании анализа экспериментальных данных по свободным колебаниям шин делается ряд основополагающих заключений, а именно

- рассеяние энергии за счет работы внутреннего трения в резинокордной оболочке шины следует представлять независящим от скорости ее радиальной деформации;

- площади характеристик упругости (петель гистерезиса) с хорошим приближением пропорциональны квадрату амплитуд радиальной деформации шин;

- площади петель гистерезиса, снятые при фиксированной амплитуде деформации и различных частотах нагружения радиальной силой (до 4 Гц), практически не отличаются друг от друга.

Свободные колебания масс стенда на испытуемой шине описывались дифференциальным уравнением вида т'Л = ~СшЛст - к ■ (4 - Лст ) • (0,5signЛ - 0,5) + вК, (В.6) где т - масса колебательной системы; сш - нормальная жесткость шины; к - коэффициент трения;

Л-н > Кт ~ начальная и статическая (равновесная) деформации; 0К - вес колеса.

Сила неупругого сопротивления в шине строится на модели сухого (но не Кулонова) трения и имеет вид

Хр - -к • \Я\ • з1§пХ, (В.7) где А, Я - деформация и скорость деформации шины, причем А=1н-Хст, 1 при X > 0 и signX = -1 при X < 0.

В работе А.Е.Межова [72] оценка поглощающих свойств шины для катящегося колеса строится на модели Кулонова трения, когда шина моделируется при помощи системы пружин и демпферов сухого трения, имитирующих работу элементарных участков шины и равномерно распределенных по ее периметру. Составляя выражения элементарных сил упругости и демпфирования с последующим их интегрированием по длине пятна контакта, были получены выражения силы упругости и демпфирования для всей шины, т.е.

Е = А-къ'\ Р = св.8)

V 0 г а где А - константа для данной шины, зависящая от жесткости шины и радиуса колеса: к - радиальная деформация шины; к - коэффициент демпфирования; Уа - скорость качения колеса;

Г^, Г^ - динамический радиус колеса и скорость его изменения.

На основании анализа полученных выражений в работе делаются следующие выводы:

1) сила упругости изменяется от деформации шины по нелинейному закону;

2 кгд

2) коэффициент демпфирования К, =-пропорционален площади а пятна контакта и обратно пропорционален скорости качения колеса;

3) сила неупругого сопротивления пропорциональна вертикальной скорости оси колеса и, следовательно, правомерна гипотеза о вязком трении в шине.

Несмотря на то, что полученное в работе [72] выражение силы демпфирования достаточно хорошо согласуется с данными других исследований, например, [57, 125] ценность результатов значительно снижается в силу целого ряда допущений и ограничений, в том числе и на диапазон изменения скорости качения колеса.

Значительное число исследований, особенно за рубежом, посвящено экспериментальному определению поглощающей способности шин с учетом и вращения колеса, и его вертикальных колебаний. Обширные исследования шин проведены в работах W.-D. Hahn [140, 141], в которых определяется взаимосвязь между параметрами упруго-вязкой модели шины и условиями ее работы. Эксперименты проводились с шинами легковых автомобилей радиальной и диагональной конструкции при различной нагрузке на шину, давлении воздуха в шине, скорости качения колеса и амплитуды вертикальных колебаний. Общие закономерности зависимостей коэффициента нормальной жесткости и коэффициента демпфирования от скорости качения и частоты колебаний колеса иллюстрируются рис.В.6 [18].

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

Исследования проводились при помощи гидропульсатора и барабанного стенда с заданием частоты колебаний в диапазоне 0 . 20 Гц и скорости качения до 100 км/ч. Механическая модель установки приведена на рис.В.7 [96].

Динамическая жесткость шины определялась следующей зависимостью: 1 cd=c0-av- ad—j==, . (В.9)

-Jl + d2 где С о - жесткость невращающейся шины при отсутствии вертикальных колебаний (С0 — Fö - статическая нагрузка на шину;

S - нормальный прогиб шины;

V^KM/Ч

5 10 20 50 Q,рад/с

Рис.В.6. Зависимости отношения динамической жесткости к статической: а - при V= const; б - при Q = const; в - при Ql V= const; г - коэффициента сопротивления при вертикальных колебаниях в зависимости от скорости качения V и частоты колебаний Q (шина 6.40-13, pw = 0,17 МПа, Р2 = 4кН). sit) м, т s(t);F*(t) 1 т — nt) а) б)

Рис. В.7. Механическая модель установки для испытания шин: а - функциональная схема установки; б - диаграммы изменения возбуждающей силы и деформации; 1 - изменение силы 2 - изменение деформации - динамическая жесткость шины; К - коэффициент вязкого трения; т - масса колеса; Г- гидропульсатор; Мз - месдоза. ау - коэффициент изменения жесткости шины за счет скорости качения колеса (ау—су/с0)\ - жесткость вращающейся шины при отсутствии вертикальных колебаний;

С1 ^ - абсолютный коэффициент повышения жесткости шины; с! - коэффициент потерь (с1=tgS=сок/с; д - сдвиг фаз;

СО - круговая частота изменения нагрузки. Коэффициент демпфирования определялся зависимостью вида а к — с0 • ау а

В. 10) лА+б/2'

Схожие исследования упругих и демпфирующих свойств автомобильных шин представлены в книге Дж.Вонга да]. Шина моделировалась упруго-вязкой моделью. Показатели жесткости и демпфирования для катящейся шины определялись методом сбрасывания. Результаты обработки кривых затухающих колебаний для транспортных и внедорожных шин показали схожие с ранее проводимыми исследованиями закономерности, т.е. жесткость при колебаниях шины выше, чем жесткость при статическом обжатии и возрастает как при росте нагрузки на шину, так и внутреннего давления воздуха в шине. Коэффициент демпфирования с ростом нагрузки на шину возрастает, а с ростом внутреннего давления в шине уменьшается. Показатели жесткости и демпфирования катящейся шины обычно определяются путем измерения реакции шины на известные гармонические возмущения. Реакции регистрируют на ступице колеса, а возмущения воздействуют на протектор шины. Другим методом определения динамической жесткости и демпфирования шины является измерение ее резонансной частоты при качении по барабану или ленте. На рис.В.8 приведены значения динамической жесткости различных типов шин легковых автомобилей, полученные по второму методу. Показано, что жесткость шины снижается, как только она начинает катиться. Однако, при скорости более 20 км/ч это снижение становится менее заметным. Были предприняты попытки определить взаимосвязь между статической и динамической жесткостью шины. Однако, никаких общих выводов не было сделано. Выражения (В.9) и (В. 10) носят чисто эмпирический характер и достаточно приближенно отражают связь между С0 и С^, к и С о.

Демпфирование пневматической шины [151] осуществляется главным образом вследствие гистерезиса ее материала и, в общем случае, имеет место не только Кулоново демпфирование колебаний, но и вязкостное демпфирование, а также сочетание этих двух типов гашения колебаний. Эквивалентный коэффициент вязкостного демпфирования можно получить на основании динамических испытаний. Его значение изменяется в зависимости от конструктивного исполнения шины, а также от условий эксплуатации. Установлено, что демпфирование пневматических шин, изготовленных из синтетических резиновых смесей, значительно меньше демпфирования в подвеске.

Между тем, хорошо известно [81, 89], что основные источники потерь энергии при механических колебаниях можно разделить на следующие три группы: рассеяние энергии в материале циклически деформируемого упругого элемента (пружины) за счет его несовершенной упругости или так называемые потери энергии в материале гистерезисного типа; рассеяние энергии на трение в узлах сочленения колебательной системы, т.е. так называемое конструкционное рассеяние энергии, и, наконец, рассеяние энергии за счет потерь в окружающую среду как подвижную, так и неподвижную или аэрогидродинамические потери энергии. Известно также, что гистерезисные потери нелинейным образом зависят от амплитуды деформации материала и представляются в виде своего рода петель гистерезиса, образующиеся при циклическом обжатии упругого элемента с данной амплитудой независимо от причин, вызывающих рассеяние энергии. При этом площадью петель гистерезиса в координатах напряжение - деформация будет характеризоваться величина рассеяния энергии, т.е. величина энергии, на которую за цикл колебаний будет уменьшаться энергия единичного циклически деформируемого объема "пружины" колебательной системы. Величину полной теряемой колебательной системой энергии за цикл легко получить путем интегрирования по всему объему циклически деформируемого материала "пружины", зная напряженное состояние последней. кН/м

300

250

200 ъ.

-с Со

150

100

50

-д

Диагональная шина

Радиальная шина

Специальная рпдиппьная 'шина

Давление в шинр 2Ч5кОа 1ч7кПа-• о

20 40

Скорость

60 км/я

Рис.В.8. Зависимость динамической жесткости шины от скорости качения (шина 165-400; диаметр барабана 2,5 м: амплитуда.неровностей 5-10 мм)

Такой подход к расчету колебаний механических систем с учетом рассеяния энергии в колебательной системе, базирующийся на оценке потерь энергии площадью петли гистерезиса, известен давно [33, 45, 55, 61]. Однако, теоретическая разработка такого подхода с целью получения аналитического описания силы неупругого сопротивления, исходя из уравнения петли гистерезиса, сдерживалась нелинейностью выражения диссипативной силы, а, следовательно, трудностями введения этой силы в дифференциальные уравнения движения механических систем и их решения. При появлении и широком внедрении в практику инженерной деятельности персональных ЭВМ интерес к данной проблеме возрос. Появились исследования, в которых была разработана и успешно апробирована нелинейная модель учета гистерезисных потерь в материале (первоначально в сталях) при расчете механических колебаний. В исследованиях Г.С.Писаренко [88] эта модель была обобщена на потери энергии в колебательной системе любого типа, т.е. на гистерезисное, аэрогидродинамическое и конструкционное рассеяние энергии при колебаниях.

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

Аналитическое выражение предлагаемой модели поглощения энергии колебаний, отражающее нелинейную зависимость между напряжением о и относительной деформацией материала "пружины" £ с амплитудой для восходящего и нисходящего движения, приводящего за цикл к образованию петли гисгде Е - модуль упругости материала;

27 - декремент колебаний, получаемый на основании экспериментальных данных, в общем случае, никак не связанный с природой рассеянйя энергии. Следовательно, ) + 8(с) + 8{К) + ., (В. 12) где: ё(£с) - декремент, зависящий от амплитуды деформаций; - декремент, зависящий от скорости деформаций; д(с) - декремент, зависящий от жесткости системы; б (К) - декремент, зависящий от сухого трения.

В обобщенном виде, т.е. в координатах сила Р - деформация/^ выражение (В.11) имеет вид з „ г /2Т

Р = са +2/ 8 Л

В.12) V где с - жесткость колебательной системы; /а - амплитуда деформации.

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

В ряде исследований [33, 81] рассеивание энергии в материале при циклическом деформировании упругого-элемента также представлялось площадью петли гистерезиса. При этом для описания закона изменения силы неупругого сопротивления в деформируемом элементе петля гистерезиса аппроксимировалась уравнением эллипса. Функциональная зависимость такой закономерности имеет вид

В. 13) о где е, е0 - текущее и максимальное значения деформации; коэффициенты модели сопротивления, зависящие от свойств материала упругого элемента.

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

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

Первые систематические исследования поглощающей способности автомобильных шин с позиции математического описания экспериментально полученных характеристик упругости как петель гистерезиса были выполнены под руководством Н.Н.Яценко [130, 131, 132]. В этих исследованиях тщательно обоснованы предпосылки представления экспериментальной характеристики упругости шины эллипсом с большим эксцентриситетом, а также взаимосвязь размахов силы неупругого сопротивления с радиальной деформацией шины; получены основные математические соотношения для демпфирующей силы; разработана методика экспериментального определения параметров модели поглощающей способности шины, которая может быть названа нелинейной эллиптической моделью поглощающей способности пневматической шины. Первые положительные результаты эллиптической аппроксимации упругих характеристик шин, полученные на очень ограниченном экспериментальном материале, стали хорошим фундаментом последующих теоретических и экспериментальных исследований поглощающей способности-шин.

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

Исследования обтекающей способности автомобильного колеса как жесткого, так и эластичного представляют интерес для определения нагружен-ности и плавности хода транспортных машин и тракторов, эксплуатируемых по бездорожью, и в настоящее время практически завершены усилиями как отечественных (А.К.Бируля, А.К.Фрумкин, В.А.Щетина, В.М.Гольдштейн и др.), так и зарубежных (К.К1орре1, Н.Морре!!, МЛиНеп, ¡.РаиЬеп и др.) ученых.

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

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

Первые исследования взаимодействия эластичной шины с короткими неровностями дороги посвящены определению жесткостных свойств шин при обжатии на неровностях различной конфигурации и траектории оси колеса при проезде этих неровностей [98]. Полученные при испытаниях шин на неровностях характеристики (рис.В.Ю, а, б) дают только косвенное представление о сглаживающей способности шин, поскольку отражают определенные закономерности перемещения оси колеса и нормальной жесткости шин от характера и степени изменения профиля неровностей. а - статическая характеристика радиальной жесткости шины при обжатии на неровности; 1 - шина 8.25-20 и треугольная неровность; 2, 3 и 4 - шина 200-20 и неровности треугольная, полуцилиндровая, полусферная радиуса К; б - траектории оси колеса с шиной 200-20 при проезде неровности в виде полусферы (Я - 6 см); 5 - жесткое колесо; 6, 7, 8 - скорости качения упругого колеса Уа= 6, 14, 33 км/ч.

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

Так, в работе Э.И.Толстопятенко [118] делается попытка оценить влияние сглаживающего эффекта пневматической шины (в статье эффект назван нивелирующим, что не совсем корректно) на амплитудно-частотную характеристику колебательной системы, эквивалентной передней (задней) части землеройно-транспортной машины. Делаются допущения о том, что в пятне контакта шины отрыв протектора от хмикропрофиля дороги не происходит, скорость качения колеса постоянна, жесткость элементов шины по периметру колеса неизменна, длина пятна контакта постоянна. В одномассовой колебательной системе шина представляется бесконечным множеством эквивалентных пружин, заменяющих ее опорную часть - протектор (рис.В.11).

Рис. В. 11. Схема эквивалентной одномассовой колебательной системы: /г- - расстояние до /-го элемента; / - длина отпечатка шины.

Дифференциальное уравнение движения такой одномассовой системы имеет вид 00 . 00

- - 4 ) • Мг + - ) • С,. = О, (В.14)

1 /=1 где т - часть массы машины, приведенная к центру колеса; ш' ш > " вертикальные перемещение, скорость и ускорение центра колеса; у, у - вертикальные смещение и скорость смещения профиля дороги;

С/ - коэффициент нормальной жесткости эквивалентной пружины;

1 - коэффициент неупругого сопротивления эквивалентной пружины.

00 00

Учитывая, что £ = > 2 сг = сш > где Рш ~ коэффициенты 1 1 жесткости и неупругого сопротивления шины, а также принимая, прямоугольный закон распределения упругих и демпфирующих элементов в зоне контакта и постоянство их характеристик по длине контакта, была получена амплитудно-частотная характеристика в виде

1тт„ ч. (27. ю-Л 2

Ж(7й))\£ш= —г-БШ— -со\1--; ^-тДВ.15) у

CD-l

2¥ У где V - скорость продольного качения колеса; со - условная собственная частота системы. Правую часть зависимости (В. 15) можно записать как произведение двух составляющих, т.е.

К(У ,а>\-\Я(1(о% (в. 16) где - амплитудно-частотная характеристика колебательной системы без учета сглаживающего эффекта шины; ЩУМ - функция (оператор) сглаживания, равная ' , IV . (О • I \K{V,o))\ = ——sin

В. 17) со • / IV

На рис.В.12 представлены графики функций при скорости качения колеса 2 м/с.

МО с'2 т то то

IffMl

600

- Г 1 1

JRMI

КМ У

Рис.В.12. Амплитудно-частотная "\К(ъ-,о)\ характеристика колебательной системы без учета сглаживания и функция сглаживания

Анализ полученных результатов позволил автору сделать следующие выводы:

1) приведенная (рис.В.11) модель шины позволяет вывести уравнения колебаний системы, эквивалентной передней (задней) части землеройно-тран-спортной машины, с учетом сглаживающей способности шины;

2) амплитудно-частотная характеристика такой системы зависит от скорости движения в силу наличия сглаживающей способности шины;

3) сглаживающая способность оказывает существенное влияние на колебания машин при движении на малых скоростях (до 20 км/ч).

Похожее выражение функции сглаживания (В. 17) получено в работе М.Д.Агеева, названное передаточной функцией нивелирования колеса, которая выведена из рассмотрения процесса качения колеса по синусоидальным неровностям и имеет вид sin0L/2

0K(ß)= , (В. 18) eij 2 где 0 - волновая частота;

1К - длина отпечатка шины, которая также принята постоянной в процессе качения колеса по неровной дороге. Несколько другой подход к оценке сглаживающей способности пневматической шины при взаимодействии с неровной дорогой случайного микропрофиля предложен в работах Н.Н.Яценко [36, 129]. При тех же допущениях о безотрывном качении колеса с постоянной скоростью и неизменной длине пятна контакта шины с дорогой шина моделируется тонкой, идеально эластичной оболочкой, плотно облегающей каждую неровность в зоне контакта, с характеристиками упругости и демпфирования сосредоточенного элемента (рис.В.13).

Тогда для плоского сечения такой модели шины вертикальное воздействие дороги на колесо выражается средним значением функции q(l) на длине отпечатка шины а, т.е. 1+а/2

4сг( 0 = -- \q(l)dl, (В.19) а /-а/2 где qcz(l) - функция сглаженного шиной микропрофиля.

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

Sq (А, а) = -уТ'С1' cos fcO'Sq №> (R2°)

Л a где Sq(X) - спектральная плотность исходного микропрофиля;

Sq сг(к., а) - спектральная плотность сглаженного микропрофиля;

Я - волновая частота.

В качестве примера на рис.В.14 изображен эффект сглаживания в среднем при расчете спектральной плотности микропрофиля специальной дороги автополигона для двух шин разной размерности [127]

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

Sq(\),CM3M/paä Л ч

V 4^7 S. 2 J /

К

8 Ю 12 П

16 18 Л,'/м

Рис.В. 14. Спектральные плотности микропрофиля специальной дороги из крупного булыжника: 1 - без учета сглаживающей способности шин; 2 — сглаженного шиной 12.00-18; 3 - сглаженного шиной 15.00-20.

Основным недостатком рассмотренных моделей постоянного сглаживания шины, построенных на осреднении либо микропрофиля и его воздействия [118], либо характеристик упругости и демпфирования [65] по длине отпечатка, является то, что введение этих моделей в расчет колебаний автомобиля приводит к появлению нулей на графиках спектральных плотностей реакций эквивалентных колебательных систем, не подтверждаемые экспериментально. Поэтому в ряде исследований сделана попытка учесть динамику изменения пятна контакта шины и, в частности, его длины при движении автомобиля по неровным дорогам. Так в работе Ю.В. Степанова и др. [112] рассматривается одно-массовая колебательная система, в которой пневматическая шина имитируется набором упругих и демпфирующих элементов. Кроме общепринятых в таких исследованиях допущений также считается, что связь между упругими элементами, лежащими в различных вертикальных сечениях колеса, отсутствует, а величины коэффициентов удельной жесткости и удельного сопротивления постоянны в каждом сечении колеса. Динамическая модель колеса приведена на рис.В.15.

Рис.В. 15. Динамическая модель колеса: ^-абсолютное вертикальное перемещение оси колеса; д*- вертикальное перемещение оси колеса, обусловленное профилем неровностей дороги; с, к- удельные коэффициенты жесткости и сопротивления элемента шины.

Учитывая допущение о пологости профиля неровностей, допускающего аппроксимацию его тремя членами ряда Тейлора, и постоянство площади упру-годеформированной части шины при данной нагрузке, авторами получено уравнение вертикальных перемещений оси колеса с учетом сглаживающего эффекта шины при переменной длине контакта в виде д^{х) = д{х) + д'2{х)

1 Я

-д" (х) 9

32 Я

1-3/1 -д"(х)Я), (в.21) где д{х), д'{х\ д"{х) - функция микропрофиля и ее производные по пути движения;

Я - радиус колеса; площадь упруго деформированной части шины. На основании (В.21) записывается дифференциальное уравнение вертикальных колебаний массы т, которое, в общем случае, получается нелинейным и нестационарным, что значительно затрудняет его решение в операторной форме и делает невозможным оценку статистических характеристик реакций системы по передаточной функции. Результаты решения дифференциального уравнения движения одномассовой колебательной системы на аналоговой вычислительной машине при некоторых упрощениях приведены на рис.В.16.

16 см 2 еглажи&ания у \ со сгп тиканием \ ч 2

20

8 раЗ/с 36

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

В другой работе [26], выполненной на кафедре "Автомобили" БПИ Ю.Ю.Беленьким и др., в качестве сглаженного микропрофиля дсг принимается траектория оси колеса при достаточно медленном обкатывании по неровностям исходного микропрофиля, которая определяется в каждой точке равенством между вертикальной составляющей реакции опорной поверхности на шину Яг и приложенной к колесу статической нагрузкой Сглаженный микропрофиль принимается за кинематическое возмущение и вводится в расчет колебаний автомобиля в предположении точечного контакта шины с дорогой.

Шина представляется в виде бесконечного множества однотипных элементарных радиально расположенных пружин (рис.В.17).

С целью экономии времени счета процедура вычислений сглаженного микропрофиля осуществлялась отдельно от расчетов колебаний автомобиля и хранилась в памяти ЭВМ самостоятельным массивом. Для проверки предложенной методики сглаживающего эффекта шины были рассчитаны среденквад-ратические ускорения на раме автомобиля МАЗ-5549 при движении по булыжной дороге ровного замощения автополигона НИЦИАМТа. Результаты расчетов (рис.В.18) показали, что среднеквадратические ускорения на раме над передней осью в диапазоне частот 0,9.22,4 Гц с учетом сглаживающей способности шины в среднем на 30 % меньше, чем при точечном контакте. Особенно сильно сглаживающая способность проявляется на высоких частотах, где уменьшение ускорений доходит до 1,5.3 раз. Косвенное сравнение результатов расчетов с экспериментальными, проведенными на автополигоне НИЦИ

АМТа на самосвалах МАЗ, показало их хорошую сходимость. у

Рис.В.17. Модель шины и процесс взаимодействия ее с микропрофилем дороги: Здеф ~ площадь упругодеформированной части шины; Во - центр тяжести упругодеформированной части шины; г - свободный радиус колеса.

L КС1

075

0,50

0,25

ЧсС 6

10 Гц

Рис.В.18. Среднеквадратические ускорения в третьоктавных полосах частот на раме автомобиля над передней осью при движении по булыжнику со скоростью Уа=\6,1 м/с:-точечный контакт;

------сглаженный профиль.

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

В работе В.С.Васильева и др. [20] подход к учету эффекта сглаживания шиной неровностей дороги основан на осреднении упругих характеристик шины по длине отпечатка, принимаемого постоянным при движении машины. Частотная характеристика сглаживающего преобразования имеет вид sin Q - ln

В.22)

H(D) = ni о где Q=2k/Ih - дорожная (волновая) частота;

1Н - длина неровности; 21н - длина отпечатка шины с дорогой. Чтобы учесть динамику изменения длины пятна контакта шины, авторы предложили сгладить амплитудно-частотную характеристику эффекта (В.22) огибающей вида

Нсг(П) = а. п+п.

В.23) где &в=1,1/10 ~ волновая частота, соответствующая длине волны Графики амплитудно-частотных характеристик (В.22) и (В.23) приведены на рис.В. 19.

Нш(Не) 1,0

Рис.В. 19. Амплитудно-частотные характеристики сглаживающего эффекта шины:

1- АЧХ постоянного сглаживания;

2- огибающая.

В монографии [35] коллектива авторов под руководством А.А.Хачату-рова на основании модели сглаживающего эффекта шины вида (В.23) были проведены исследования колебаний автомобиля ЗИЛ-130 на специальных дорогах автополигона НИЦИАМТа. Исследования показали, что целесообразность учета сглаживающей способности шины зависит от передаточной функции динамической системы, эквивалентной реальному автомобилю. Кроме того, при малых скоростях Уа сглаживающая способность шины влияет на параметры колебаний (ускорения масс, деформации шин и др.); при больших Уа сглаживающая способность шины практически не отражается на колебаниях автомобиля, т.к. он реагирует на неровности, длина которых больше 1в=2тт/Ов.

Скорость, до которой следует учитывать сглаживающую способность шины,

Г0 = О,9о)м10, (В.24) где СОм - максимальная частота в спектре возмущения, на которую реагирует система.

На рис.В.20 показано влияние сглаживающего эффекта шины на вертикальные ускорения точки кузова над задней осью порожнего автомобиля ЗИЛ-130.

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

6г. м/сг

10 р / А о

Л г

•

У'» i

Рис.В.20. Влияние сглаживающей способности шины на ускорения передних подрессоренных масс автомобиля ЗИЛ-130: 1 - разбитая булыжная дорога; 2 - булыжная дорога в удовлетворительном состоянии; 3 - асфальтированное шоссе; -с учетом сглаживания;

------без учета сглаживания.

25 50 уа,км/ч

Еще один прием преодоления "нулей" в методике учета постоянного сглаживания шины приведен в стандарте предприятия AMO ЗИЛ "Воздействия дорожные. Метод формирования массивов микропрофилей дорог". Учет сглаживающей способности шин осуществляется путем предварительного преобразования исходного микропрофиля дороги функцией сглаживания, моделируемой одномассовой колебательной системой с заданными параметрами (механическим фильтром). Схема такого фильтра приведена на рис.В.21.

Передаточная функция фильтра имеет вид 0 s +2y/ú)0s + Сф где COq — Сф / ТПф - собственная частота колебаний;

В.25) у/ — кф / 2^Сф • ТПф - коэффициент апериодичности.

Для реализации необходимого сглаживания задаются значениями параметров: Шф=\ кг; у/=1; Ьс=0,5 м, где Ьс - длина волны дорожной неровности, соответствующая собственной частоте системы.

Тогда

17тУ" л тл

Уп = а

L. Ф m

Ф' V

L.

1 бтт-V а J кф = 2у/л]сф •тф

Пример сглаживания микропрофиля цементобетонной дороги рассмотренным фильтром приведен на рис.В.22.

Рис.В.21. Схема механического фильтра для моделирования функции сглаживания: д^) - ордината воздействия исходного микропрофиля; ЦфО) ~ ордината сглаженного воздействия. а) б)

Рис.В.22. Амплитудно-частотная характеристика механического фильтра и спектральные плотности цементобетонной дороги: а — АЧХ фильтра; б - спектральные плотности микропрофиля; 1 - исходный микропрофиль; 2 - сглаженный микропрофиль.

Последние два предложения по учету сглаживающей способности шин преодолевают парадокс нулей модели постоянного сглаживания либо переходом от частотной характеристики к уравнению ее огибающей [20], либо путем искусственного ввода механического фильтра с частотной характеристикой, близкой к закону огибающей характеристики эффекта постоянного сглаживания (стандарт ЗИЛа). Хотя оба предложения улучшают выходные параметры сглаженного воздействия, но могут быть приняты чисто условно, т.к. опираются на недоказательные соображения.

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

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

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

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

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

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

Для реализации поставленной цели необходимо решить следующие задачи:

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

2. Обосновать и построить математическую модель корреляционных функций микропрофилей основных испытательных дорог автополигона НИЦИАМТа по экспериментально полученному массиву ординат микропрофилей этих дорог; рассчитать на ЭВМ численные значения коэффициентов модели.

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

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

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

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

Разработать программное обеспечение для персональных ЭВМ численного решения математической модели двухмассовой колебательной системы с учетом эффекта переменного сглаживания и моделирования случайного воздействия, эквивалентного воздействию микропрофиля специальных испытательных дорог автополигона НИЦИАМТа.

Рассчитать и провести графический анализ амплитудно-частотных характеристик звена переменного сглаживания колебательной системы, эквивалентной передней подвеске легкового автомобиля ЗИЛ-4102, при движении с различной скоростью и по различным дорогам автополигона.

5.Разработать программное обеспечение для ввода предлагаемой модели переменного сглаживания шин в программную систему отдела главного конструктора (ОГК) САПР ЗИЛа для обеспечения формирования дорожного воздействия и получения спектральных характеристик плавности хода и нагруженности проектируемого автомобиля с учетом эффекта сглаживания.

Рассчитать спектральные плотности реакций динамических систем, эквивалентных грузовому автомобилю ЗИЛ-4331 и легковому автомобилю ЗИЛ-4102, в различных точках несущей системы при воздействии микропрофилей дорог автополигона НИЦИАМТа с учетом эффекта переменного сглаживания шин.

35

Провести ходовые испытания автомобилей ЗИЛ-4331 и ЗИЛ-4102 на специальных дорогах автополигона НИЦИАМТа, обеспечив запись колебаний в различных точках несущей системы, мостов, элементов подвески и виброзащиты, а также их обработку с целью получения спектральных характеристик плавности хода. сравнить результаты ходовых испытания автомобилей с результатами расчетов по программной системе ОГК САПР ЗИЛ.

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

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

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