CC BY
64
УДК 629.114.2
КОНЦЕПЦИЯ И МЕТОДЫ ОБЕСПЕЧЕНИЯ НАДЕЖНОСТИ МОБИЛЬНОЙ ТЕХНИКИ ПРИ ПРОЕКТИРОВАНИИ
А.А. Абызов, И.Я. Березин
Предложена концепция обеспечения надежности мобильной техники, основанная на компьютерном моделировании процесса эксплуатации. Представлена математическая модель гусеничной машины, ориентированная на решение задач прочностной надежности ходовой системы. Модель описывает динамику подрессоренного корпуса и силовой установки при прямолинейном и криволинейном движениях.
Ключевые слова: гусеничные машины, моделирование процессов эксплуатации, криволинейное движение, динамические нагрузки.
Обеспечение конкурентоспособности изделий отечественного машиностроения непосредственно связано с необходимостью повышения уровня научно-технического сопровождения инновационных проектов. В значительной мере это относится к вопросам обеспечения прочностной надежности машин, которая закладывается в процессе выполнения рабочих проектов. Как показывает практика, в большинстве случаев инженерные расчеты ограничиваются определением максимальных напряжений от действия предельных нагрузок и затем путем их сопоставления с характеристиками статической прочности материалов оцениваются коэффициенты запаса ответственных деталей проектируемой конструкции. Очевидно, такой подход ограничивает возможность предсказания усталостных, износовых, термомеханических и других видов повреждений, которые, как показывает опыт, преимущественно являются причиной отказов тяжелонагружен-ных элементов машин. С другой стороны, следует отметить, что предусмотренные техническими заданиями сертификационные характеристики надежности изделий выявляются лишь путем анализа статистической информации, получаемой при эксплуатации серийной продукции.
Применение современных достижений в областях исследований динамики машин, прочности конструкций и вычислительных технологий позволяет принципиально изменить подход к прогнозированию и управлению надежностью изделий путем широкого и опережающего применения методов математического моделирования задач динамики и прочности машин. Постановка задачи прогнозирования прочностной надежности предусматривает последовательное выполнение следующих этапов:
- описание многофакторного случайного воздействия внешней среды на рабочие органы и динамическую систему проектируемых изделий (ветровая нагрузка на высотные сооружения, воздействие профилей дорог на ходовую часть мобильных машин, воздействие со стороны грунта на рабочие органы дорожно-строительной техники т. п.);
- компьютерное моделирование процессов реальной эксплуатации машин;
- определение статистических характеристик силового и кинематического взаимодействия элементов конструкции;
- конечно-элементное моделирование полей напряженно-деформированного состояния с целью выявления потенциально опасных зон ответственных деталей;
- моделирование процессов формирования усталостных, износовых, термомеханических и других видов отказов с учетом случайного характера нагружения и рассеяния прочностных свойств применяемых материалов;
- определение количественных характеристик надежности, таких как функция вероятностей безотказной работы, средняя наработка до отказа, гамма-процентный ресурс тяжелонагруженных элементов машин и т. п.;
- корректировка исходных данных проекта по динамическим свойствам проектируемых изделий и прочностным характеристикам деталей с целью обеспечения требуемых показателей надежности.
Комплекс названных исследований представляет собой инженерную методику имитационного моделирования испытаний. К числу наиболее значимых преимуществ подхода следует отнести возможность проведения параметрического прогноза и количественной оценки влияния характеристик проектируемой машины на эксплуатационную нагруженность и долговечность ответственных элементов. Его применение также позволит повысить эффективность проведения стендовой отработки тяжелонагруженных узлов и обосновать программу натурных испытаний опытных образцов машин.
Предлагаемое исследование иллюстрирует реализацию метода имитационного моделирования испытаний на примере тяжелонагруженных элементов ходовой системы быстроходных гусеничных машин, предназначенных для эксплуатации в условиях необустроенной местности (рис. 1). Актуальность задачи обусловлена тем, что надежность таких машин в значительной степени определяется долговечностью гусеничного движителя, работающего в весьма тяжелых условиях.
Одним из основных элементов рассматриваемой методики является математическая модель гусеничной машины. Применяемые в настоящее время математические модели шасси гусеничной машины можно разделить на несколько групп.
Модели, предназначенные для исследования плавности хода, достаточно подробно описывают подрессоренный корпус. При этом обычно рассматривается только режим прямолинейного движения, силовая установка моделируется одним инерционным элементом [1-3].
Математические модели, ориентированные на исследование динамики силовой установки, отличаются подробным описанием трансмиссии и двигателя машины, при этом подрессоренный корпус описывается достаточно приближенно [4, 5].
При исследовании процессов криволинейного движения уделяется большое внимание описанию взаимодействия гусеницы и грунта, а также работе механизма поворота; корпус машины при этом обычно рассматривается как тело, совершающее плоское движение [6-9].
Для реализации методики имитационного моделирования испытаний предложена математическая модель, ориентированная на решение задач прочностной надежности ходовой части гусеничной машины. Модель характеризуется следующими отличительными особенностями:
- модель описывает динамику связанной нелинейной системы «грунт - гусеница - подрессоренный корпус - силовая установка - водитель»;
- учитывает существенные нелинейности системы, такие, как односторонние связи в гусеничном обводе и опорных катков с грунтом, нелинейные характеристики двигателя, фрикционных элементов и гидромеханических передач;
- отображает изменение структуры силовой установки при переключении передач и работе механизма поворота;
- учитывает варьирование скорости в широком интервале значений при движении машины по труднопроходимой местности [10];
- для расчета сил сопротивления повороту использовано конечно-элементное моделирование взаимодействия траков гусеницы с грунтом, что позволяет учитывать особенности формы опорной поверхности траков, моделировать процессы погружения трака и нагребания грунта [11, 12].
Предложенную модель можно условно рассматривать как совокупность связанных частей:
- основной подсистемы, описывающей динамику подрессоренного корпуса машины;
- модуля для расчета сил взаимодействия опорной поверхности гусеницы с грунтом при прямолинейном и криволинейном движениях;
- подсистемы, описывающей динамику силовой установки машины;
- технико-эргономической подсистемы, описывающей работу механика-водителя по приведению скорости движения машины в соответствие с дорожной ситуацией.
Схема подсистемы, моделирующей подрессоренный корпус машины, представлена на рис. 2,
3, а. Корпус рассматривается как твердое тело, имеющее шесть степеней свободы. С центром тяжести корпуса о связана подвижная система координат о - х - у - г. Положение на местности рассматривается в неподвижной системе 01 - X - У - ^ В качестве обобщенных координат приняты координаты центра тяжести корпуса в неподвижной системе X, У, Z, углы продольно- и поперечноугловых колебаний 3, Д а также угол поворота относительно вертикальной оси «(рис. 3, а). Корпус имеет массу тк и моменты инерции 1х, 1у, 1г относительно соответствующих осей.
МОДЕЛИ ФОРМИРОВАНИЯ ОТКАЗОВ. ХАРАКТЕРИСТИКИ НАДЕЖНОСТИ
Рис. 1. Блок-схема комплексного моделирования
Рис. 2. Расчетные схемы гусеничного движителя, подрессоренного корпуса и силовой установки
На корпус машины действуют сила веса т^, вертикальные усилия со стороны подвесок Иу (7 = 1...2 - номер борта,у = 1...п - номер подвески), силы натяжения рабочих и свободных ветвей гусениц ТрЬ Тс, а также сила сопротивления прямолинейному движению Я. При движении по криволинейным участкам трассы на опорные катки машины действуют нагрузки Еуу, которые через балансиры подвески также передаются на корпус машины.
Подвески машины моделируются эквивалентными вертикальными упругими и демпфирующими элементами с нелинейными характеристиками. При этом учитываются отрывы опорных катков и пробои подвесок на упорах. При движении без отрыва вертикальная координата оси опорного катка принимается равной высоте микропрофиля трассы в соответствующей точке. В случае отрыва соответствующая координата рассчитывается путем интегрирования уравнения движения массы катка. Наряду с усилиями со стороны подвесок на крайние катки действуют усилия со стороны гусениц: Тс, Трг-. Гусеничные ленты моделируются нелинейными упругими элементами одностороннего действия.
При моделировании криволинейного движения машины важной задачей является адекватное определение сил сопротивления, действующих на опорные поверхности гусениц. Для их расчета в настоящее время используется несколько принципиально разных подходов.
Первая группа методов предполагает расчет момента сопротивления повороту и боковой силы в зависимости от свойств грунта, радиуса и смещения центра поворота по эмпирическим формулам [13, 14]. Опыт исследований показывает, что результаты расчета с использованием такой модели удовлетворительно согласуются с результатами испытаний тихоходных машин.
Другой подход основан на использовании математической теории трения [7, 15]; при этом силы взаимодействия гусеницы с грунтом рассматриваются как нагрузка, распределенная по опорной поверхности. Очевидно, что такой метод позволяет получить адекватные результаты в случае движения машины по жесткой, практически недеформируемой поверхности.
В предлагаемой математической модели гусеничной машины использован принципиально иной подход для описания взаимодействия гусеницы с грунтом [6, 8]. В его основу положено экспериментально обоснованное допущение о том, что взаимодействие с грунтом происходит на так называемых «активных» участках, находящихся в окрестности опорных катков машины. При этом распределенную по длине гусеницы нагрузку заменяют сосредоточенными силами, действующими на активных участках ЕХі- и руі- (і = 1, 2 - номер борта машины,] - номер опорного катка, рис. 3, а). Использование такого подхода позволяет учесть реальное распределение вертикальных нагрузок по длине опорной ветви, совместно рассматривать продольное и поперечное проскальзывание, а также определять поперечные нагрузки, действующие на каждый опорный каток машины.
Р(1) 1
Т т т 4^ і
б)
Рис. 3. Нагрузки, действующие на опорную поверхность (а); нагрузки и перемещения на активном участке (б)
Силы ЕХі- и Еуу формируются как сумма нагрузок , рУ- (I - номер трака), действующих
на траки в активном участке. Распределение нагрузок рХ-, рУ- , а также смещения траков 4)
схематично показаны на рис. 3, б. Смещения Л-) могут быть вычислены путем интегрирования по времени скорости проскальзывания трака относительно грунта за время прохождения активного участка. В предлагаемой модели нагрузки рХ- , ру— рассчитывают в зависимости от каткові)
вой нагрузки и смещения Л- с помощью соотношений, полученных по результатам конечноэлементного расчета процесса взаимодействия трака и грунта [12]. В настоящее время метод конечных элементов широко используется для моделирования взаимодействия с грунтом движителей и рабочих органов машин различного назначения [16, 17]; его применение позволяет учесть характеристики деформируемого грунта и форму опорной поверхности траков, а также адекватно описать процессы их взаимодействия, сопровождающиеся большими перемещениями. Для моделирования применен пакет программ Ь8-БУКЛ [11], предназначенный для моделирования процессов, сопровождающихся большими деформациями и разрушением материала. В расчетах использована модель материала FHWA-SOIL (*МЛТ_147), специально предназначенная для описания грунта. Свойства материала задавались на основе справочных данных.
На первом этапе были проведены тестовые расчеты взаимодействия с грунтом одиночного трака быстроходной машины (БМП-2). Конечноэлементные модели трака и грунта показаны на рис. 4, а; при этом гребни трака, не взаимодействующие с грунтом, не моделировались. В процессе расчетов трак нагружался вертикальной силой N и после погружения сдвигался в плоскости грунта в заданном направлении на заданное расстояние (рис. 4, б). В результате получены зависимости полной сдвигающей нагрузки (р), а также ее продольной (РХ) и поперечной (ру) составляющих от соответствующих перемещений.
Выполнен большой объем расчетных исследований для различных видов грунта, траков с грунтозацепами разной высоты и формы, а также для разнообразных режимов нагружения.
В качестве примера на рис. 5, а приведены зависимости силы Г от смещения А при сдвиге трака в поперечном направлении на плотном грунте (глина, целина), а также на мягком (супесь) для различных значений вертикальной нагрузки N. Эти результаты демонстрируют увеличение Г при увеличении N. Для мягкого грунта (супесь) сдвигающая нагрузка существенно меньше, чем для плотного при той же нормальной нагрузке.
Нормальная
Рис. 4. Конечноэлементная модель трака и фрагмента грунта (а); деформированное состояние грунта
при перемещении трака (б)
На рис. 5, б изображены зависимости продольной составляющей нагрузки Гх от поперечной ¥у при сдвиге трака в различных направлениях (Ы = 10 кН); направление сдвига отмечено пунктирной линией. Эти результаты показывают, что направление силы взаимодействия трака с грунтом не совпадает с направлением перемещения. По мере сдвига трака угол между направлением силы и перемещением изменяется.
На основании анализа результатов расчетов получена аппроксимирующая зависимость, связывающая силу, сдвигающую трак в плоскости дороги ¥, с перемещением А и вертикальной нагрузкой:
Г = (N / N0) ”(ф) Г, (ф) [1 - ехр(-с(ф)А)], где N0 - среднее значение вертикальной нагрузки, действующей на трак; ф - угол между направлением перемещения и продольной осью трака.
,х
а)
И - Л=0,01 м + -Д=0,03 м 0 - Д=0,02 м X - Д=0,04 м
------- годограф вектора силы Б
-------направление смещения Д
б)
Рис. 5. Зависимости силы от перемещения при боковом сдвиге трака (а); зависимости сдвигающей нагрузки
от величины и направления сдвига (б)
Также предложена зависимость, связывающая угол между сдвигающей силой и продольной осью трака у с углом ф:
у = a1 + a 2ф + a39 .
Функции га(ф), Fm^), с(ф) и постоянные коэффициенты a1, a2, a3 получены для разных видов грунта и траков с различными грунтозацепами.
В дальнейшем данные соотношения использованы в математической модели гусеничной машины для расчета сил, действующих на опорную поверхность гусеницы.
Одной из подсистем рассматриваемой математической модели является математическая модель силовой установки. В связи с разнообразием применяемых двигателей, трансмиссий и механизмов поворота, структура этой части модели должна определяться особенностями рассматриваемой машины. Поскольку модель машины ориентирована на исследование низкочастотной динамики, признано целесообразным не включать в ее состав упругие элементы. Однако для адекватного описания переходных режимов (переключение передач, изменение подачи топлива, работа механизма поворота) она должна содержать инерционные элементы, элементы, моделирующие работу фрикционов, гидрообъемных и гидродинамических передач, а также учитывать характеристику двигателя. Модель машины с механической трансмиссией представлена в [18].
Для предлагаемой математической модели получены дифференциальные уравнения движения. В качестве примера приводятся дифференциальные уравнения, описывающие движение основной подсистемы модели (подрессоренного корпуса):
X = x cos a + y sin a, Y = y cos a + x sin a,
m*(x + ui + ay) = X (Tpi T - Fci) -R mK (y -Pz -ax) = ]T £ Fyij, i=i i=i j=1
2 n
m* (z -^x + Py) = X (Z Ni, j - Tci sin Vc +Tpi sin Vp ) - mкg,
i=1 j=1
2 n
Iy^ + a(3(Iz -Ix) = X(XNi,}l] + Tpi(Lbsin Vp -Hb cos Vp + Ho) -
i=1 j=1
-Tci L sin Vc - Hb c0s Vc + Ho ) - FciHo ) + Mcy,
IxP + a6( Iy -Iz ) = 0,5B(Z (N1, j - N2, j ) + sin Vp (Tp2 - Tp1) + sin Vc (Tc2 - Tc1)),
j=1
Iza + pe(Ix-Iy) = 0,5B(Tp1 -Tp2 +Tc2 -Tc1 +Fc2 -Fc1)-£XljFyj ,
i=1 j =1
где Ho, Hb - высота над опорной поверхностью центра тяжести машины и оси ведущего колеса; Lb, Lu - расстояние вдоль продольной оси от центра тяжести машины до осей ведущего и направляющего колес; lj - продольная координата j-го активного участка.
При расчете нагрузок, действующих на активный участок, а также натяжения ветвей гусениц используется значение скорости проскальзывания опорных поверхностей i-го борта машины относительно грунта в продольном направлении vpi(t). Эта величина определяется на каждом шаге интегрирования уравнений движения на основе равновесия гусеничного обвода, в соответствии с которым разность усилий в рабочей (Tpi). и свободной (Tci ) ветвях гусеницы должна равняться
сумме продольных усилий, действующих на активных участках (Fxii) (рис. 6).
Для численного интегрирования уравнений движения разработан специальный пакет программ, с помощью которого выполнена серия тестовых расчетов. Анализ результатов и их сопоставление с экспериментальными данными подтверждает адекватность предложенной математической модели. В качестве примера на рис. 7. представлены траектории движения машины (БМП-2) при повороте с фиксированным радиусом на различном грунте (супесь и плотная глина). Расчеты выполнены для машины со стандартными траками гусениц (1, 2) и с изношенными грунтозаце-пами (3). В случае изношенных грунтозацепов, а также на мягком грунте радиус поворота получился больше, чем на жестком грунте. Это объясняется меньшим сопротивлением заносу и большим проскальзыванием опорной поверхности относительно грунта.
Рис. 6. Нагрузки в гусеничном обводе
Таким образом, применение в составе математической модели гусеничной машины предлагаемой методики расчета сил взаимодействия траков гусеницы с грунтом позволяет получить процессы изменения нагрузок, действующих на элементы ходовой части при движении машины по трассе, а также кинематические параметры криволинейного движения с учетом формы опорной поверхности траков и свойств грунта. Процессы изменения нагрузок в дальнейшем используются для расчетной оценки долговечности элементов движителя. Анализ кинематических параметров поворота может быть использован при оценке управляемости машины в различных условиях.
В математическую модель транспортной машины включена подсистема, имитирующая действия водителя по корректировке скорости движения. Скорость движения по местности задается в виде функции пути, зависящей от условий эксплуатации и динамических характеристик проектируемой машины, с одной стороны, и ряда технических и эргономических ограничений - с другой.
Рис. 7. Траектории поворота машины с фиксированным радиусом: 1 - грунт - супесь; 2 - грунт - плотная глина, стандартный трак; 3 -грунт - супесь, трак с изношенными грунтозацепами
При реализации предлагаемой методики для моделирования движения машины по трассе необходимо задание комплекса факторов внешнего воздействия, к которым, в частности, относятся микро- и макропрофиль трассы, радиусы и протяженности криволинейных участков, характеристики одиночных препятствий, характеристики грунта и др. Для задания этих параметров удобно использовать картографические описания трасс известных испытательных станций и полигонов, которыми располагают отрасли-производители сельскохозяйственной, автомобильной, дорожностроительной и другой инженерной техники.
Результаты расчетных исследований, выполненных с помощью предлагаемой математической модели, а также примеры практического применения предлагаемой методики будут представлены в следующей публикации.
Литература
1. Аврамов, В.П. Динамика гусеничной транспортной машины при установившемся движении по неровностям / В.П. Аврамов, Н.Б. Калейчев. - Харьков: Выща школа: Изд-во при Харьковском ун-те, 1989. - 112 с.
2. Дмитриев, А.А. Теория и расчет нелинейных систем подрессоривания гусеничных машин /
А.А. Дмитриев, В.А. Чобиток, А.В. Тельминов. - М. : Машиностроение, 1976. - 207 с.
3. Котиев, Г.О. Прогнозирование эксплуатационных свойств систем подрессоривания военных машин: автореф. дис. ... д-ра техн. наук /Г. О. Котиев. - М. : МГТУ, 2000. -36 с.
4. Вейц, В.Л. Динамика управляемых машинных агрегатов /В.Л. Вейц, М.З. Коловский, А.Е. Ко-чура. - М. : Наука, 1984. - 352 с.
5. Держанский, В.Б. Прогнозирование динамической нагруженности гидромеханических трансмиссий транспортных машин / В.Б. Держанский, И.А. Тараторкин // Приводная техника. -2007. - № 6 (70). - С. 29-35.
6. Красненьков, В.И. Математическая модель криволинейного движения транспортной гу-
сеничной машины по деформируемому основанию / В.И. Красненьков, С.А. Харитонов, А.В. Шумилин //Известия вузов. Машиностроение. - 1989. - № 11. - С. 94-99.
7. Позин, Б.М. Кинематические соотношения при взаимодействии движителя с грунтом при повороте /Б.М. Позин, И.П. Трояновская //Вестник ЮУрГУ. Серия «Машиностроение». - 2005. -Вып. 7. - № 14 (54). - С. 93-96.
8. Шумилин, А.В. Математическая модель криволинейного движения транспортной гусеничной машины по деформируемому основанию / А.В. Шумилин // Тракторы и сельскохозяйственные машины. - 1993. - № 5. - С. 8-11.
9. Al-Milli, S. Track-terrain modelling and traversability prediction for tracked vehicles on soft terrain /
S. Al-Milli, L.D. Seneviratne, K. Althoefer // Journal of Terramechanics. - 2010. - № 47. - Р. 151-160.
10. Савочкин, В.А. Статистическая динамика транспортных и тяговых гусеничных машин /
В.А. Савочкин, А. А. Дмитриев. - М. : Машиностроение, 1993. - 320 с.
11. LS-DYNA user’s manual. Version 970. - USA, Livermore CA: Livermore Software Technology Corp., 2003. - 1564p.
12. Абызов, А.А. Использование метода конечных элементов для моделирования взаимодействия гусеницы с грунтом при криволинейном движении машины / А.А. Абызов // Труды 15 Всерос. науч.-практ. конф. НПО спецматериалов. Т. 3: Бронетанковая техника и вооружение. - СПб., 2012. - С. 184-190.
13. Забавников, Н.А. Основы теории транспортных гусеничных машин / Н.А. Забавников. -М. : Машиностроение, 1968. - 396 с.
14. Никитин, А.О. К вопросу исследования динамических качеств гусеничной машины при повороте / А.О. Никитин // Труды ВАБТВ. - М., 1960. - № 180. - С. 13-25.
15. Гуськов, В.В. Теория поворота гусеничных машин /В.В. Гуськов, А.Ф. Опейко. - М.: Машиностроение, 1984. -168 с.
16. Asaf, Z. Evaluation of link-track performances using DEM / Z. Asaf, D. Rubinstein, I. Shmule-vich // Journal of Terramechanics. - 2006. - № 43. - Р. 141-161.
17. Hambleton, J.P. Modeling wheel-induced rutting in soils: Indentation / J.P. Hambleton, A. Drescher // Journal of Terramechanics. - 2008. - № 45. - Р. 201-211.
18. Березин, И.Я. Моделирование процесса эксплуатации при имитационных ресурсных испытаниях мобильной техники /И.Я. Березин, А.А. Абызов // Техника и технологии строительства и эксплуатации автомобильных дорог: сб. науч. тр. МАДИ. - М. : МАДИ, 2000. - С. 56-74.
Абызов Алексей Александрович. Кандидат технических наук, доцент, доцент кафедры «Прикладная механика, динамика и прочность машин», Южно-Уральский государственный университет (Челябинск). Тел.: (351) 267-92-61; abyzovaa@susu.ac.ru.
Березин Игорь Яковлевич. Доктор технических наук, профессор, профессор кафедры «Прикладная механика, динамика и прочность машин» Южно-Уральского государственного университета (Челябинск). Тел.: (351) 267-92-61; berezinii@susu.ac.ru.
Bulletin of the South Ural State University Series “Mechanical Engineering Industry” ____________________________________________________________2013, vol. 13, no. 2, pp. 19-27
THE CONCEPT AND METHODS OF RELIABILITY IN THE DESIGN OF MOBILE MACHINES
A.A. Abyzov, South Ural State University, Chelyabinsk, Russian Federation, abyzovaa@susu.ac.ru,
I.J. Berezin, South Ural State University, Chelyabinsk, Russian Federation, berezinii@susu.ac.ru
A concept of securing mobile devices, based on a computer simulation of the operation. A mathematical model of tracked vehicle, focused on problem solving strength reliability undercarriage. The model describes the dynamics of the cushioning body and power plant in rectilinear and curvilinear motion.
Keywords: track machines, process modeling, operation, curvilinear motion, dynamic loads.
Поступила в редакцию 3 июля 2013 г.
