CC BY
118
МЕТАЛЛУРГИЯ. МЕТАЛЛООБРАБОТКА. МАШИНОСТРОЕНИЕ
УДК 629.3.001.4
ПРОГНОЗИРОВАНИЕ ПАРАМЕТРИЧЕСКОЙ НАДЕЖНОСТИ ГИДРОУСИЛИТЕЛЯ РУЛЕВОГО УПРАВЛЕНИЯ МОБИЛЬНОЙ МАШИНЫ
Доктора техн. наук, профессора МРОЧЕК Ж. А., ТАРАСИК В. П., асп. МРОЧЕК Т. В.
Белорусский национальный технический университет, Белорусско-Российский университет
Одно из основных требований, предъявляемых к гидроусилителю рулевого управления, -надежность, характеризующая способность гидроусилителя выполнять заданные функции в процессе эксплуатации. Поэтому важной задачей является прогнозирование изменения технического состояния гидроусилителя при эксплуатации, т. е. прогнозирование его надежности. Решение этой задачи позволит сократить затраты времени и средств на длительные и трудоемкие испытания, более рационально использовать потенциальную долговечность гидроусилителя за счет обоснованного построения
системы его технического обслуживания и ремонта.
По данным [1, с. 60], около 60 % отказов гидропривода являются параметрическими, 40 % -функциональными. Таким образом, важное значение имеет прогнозирование параметрической надежности, т. е. надежности, связанной не с отказом функционирования изделия, а с изменением его выходных параметров во времени.
С учетом рекомендаций [1-4] авторами разработан алгоритм прогнозирования параметрической надежности гидроусилителя рулевого управления (рис. 1).
Рис. 1. Алгоритм прогнозирования параметрической надежности гидроусилителя рулевого управления
В базе данных компьютера содержится классификация факторов, влияющих на надежность гидроусилителя, т. е. классификация внешних воздействий и возможных режимов и условий работы усилителя, а также приведено описание физических процессов деградации, приводящих к отказу. Данная информация учитывается при математическом моделировании процессов функционирования гидроусилителя.
Математическая модель гидроусилителя позволяет проводить исследования процессов в гидравлической системе усилителя. С помощью этой модели возможно определение областей состояния для выходных параметров гидроусилителя и его элементов, оптимальных эксплуатационных параметров и предельных режимов эксплуатации на основе планирования многофакторного эксперимента на математической модели гидроусилителя и построения уравнений регрессий. При получении математического описания процессов в гидроусилителе на основе структурно-матричного метода [5] разработана динамическая модель, определены параметры элементов модели и составлена система уравнений математической модели. Математическая модель гидроусилителя рулевого управления включает в себя компонентные и топологические уравнения и имеет следующий вид:
= /(рш, ру1, рщ); фуг ¡Л = / ( а );
Рдг = г (а);
=/ (а),
где рвг - потенциалы внешних воздействий; руг -то же упругих элементов системы; рдг - то же диссипативных элементов; хг - перемещения механических элементов.
Условия проведения испытаний рулевых механизмов совместно с гидроусилителями определяются стандартом СЭВ 1629-79 «Грузовые автомобили и автобусы. Рулевые механизмы с гидравлическим усилителем. Технические требования. Методы стендовых испытаний». Проведенными авторами исследованиями установлено, что требований данного стандарта
иногда оказывается недостаточно для того, чтобы получить всеобъемлющую информацию об объекте испытаний и выявить неисправность. В связи с этим предложена методика испытаний гидроусилителей [6], содержащая ряд дополнений к указанному стандарту. Разработанная методика включает в себя методику общего испытания гидроусилителя (испытание рулевого механизма с гидроусилителем в сборе) и методику поэлементного испытания (определение гидравлических характеристик основных устройств гидроусилителя). Методика общего и поэлементного испытаний гидроусилителя и его устройств на стенде может применяться при испытаниях вновь проектируемых усилителей и диагностировании их в процессе эксплуатации.
В Белорусско-Российском университете разработаны принципиальная гидравлическая схема и конструкция стенда для испытаний гидроусилителей, реализующая предложенную методику. Конструкция стенда защищена патентом Республики Беларусь [7]. На рис. 2 представлены фотографии стендов для общего и поэлементного испытаний гидроусилителя. Стенды установлены в МГКУП «Горэлектротранс-порт» (г. Могилев) и УКПП «Бобруйское троллейбусное управление» и используются для испытаний всех типов гидроусилителей, устанавливаемых на троллейбусы ЗиУ-682 и АКСМ-101.
Стенд позволяет проводить испытания гидроусилителей, выполненных по различным компоновочным схемам, и получать информацию обо всех контролируемых выходных параметрах усилителя для статистической оценки технического состояния усилителя (что необходимо для принятия решения о диагнозе, т. е. отнесения по значениям диагностических параметров состояния испытуемого объекта к одному из ряда возможных его состояний). Для статистической оценки технического состояния гидроусилителя используется метод Байеса либо метод статистических решений [1, с. 227; 8, с. 58]. Метод статистических решений применяется в тех случаях, когда априорных статистических данных по отказам и неисправностям диагностируемой гидросистемы недостаточно и метод Байеса становится неэффективным.
Рис. 2. Стенды для общего и поэлементного испытаний гидроусилителей рулевого управления: а - стенд в МГКУП «Горэлектротранспорт» (г. Могилев); б - стенд в УКПП «Бобруйское троллейбусное управление»
Структурно-следственная диагностическая модель усилителя содержит структурные и диагностические параметры и позволяет по сопутствующим процессам установить характер и причины неисправности. При выборе диагностических параметров для оценки технического состояния может использоваться теория неопределенностей или энтропии [8-12]. При этом анализируются потери информации, вызванные неполнотой контроля параметров.
Для прогнозирования надежности объектов применяют [2, с. 392, 13] методы эвристического прогнозирования (экспертной оценки), методы прогнозирования по статистическим моделям и комбинированные методы, основанные на совместном применении двух первых методов с последующим сравнением результатов. В разработанном алгоритме используются методы второй группы.
Прогнозирование изменения технического состояния проводится на основе метода прогнозирования по среднему статистическому изменению параметра состояния нескольких однотипных изделий или метода прогнозирования по реализации параметра одного конкретного изделия [1, 8]. Для гидроусилителя рулевого управления авторами использован первый метод, поскольку гидроусилитель обслуживается по состоянию. При этом на основании анализа имеющихся статистических данных строится регрессионная модель изменения диагностического параметра. Регрессионная модель базируется на статистических данных об изме-
нении диагностических параметров в функции наработки. Для ее получения строится усредненная кривая изменения диагностического параметра у = f(t) в зависимости от наработки t, выполняется ее аппроксимация путем нахождения параметров кривой с использованием метода наименьших квадратов.
Выбор аппроксимирующей функции существенно влияет на трудоемкость и точность прогнозирования. Аппроксимирующая функция должна быть достаточно универсальной, относительно простой и при этом учитывать физические процессы, приводящие к изменению параметра состояния, и случайный характер скорости такого изменения, а также отражать интегральные закономерности отклонения параметра в зависимости от наработки или срока службы гидроусилителя [8, с. 225].
В [8, 11] указывается, что для агрегатов гидроприводов машин в качестве аппроксимирующей функции целесообразно использовать степенную функцию. В этом случае изменение диагностического параметра Ду определяется как
Ау (t ) = BJßl + АуЭ1 (t ) + АП у,, (1)
где Bc1 - коэффициент, характеризующий скорость изменения 1-го диагностического параметра в зависимости от наработки t; ß, - показатель степени, характеризующий интенсивность изменения параметра во всем диапазоне наработки; Ay31(t) - случайное отклонение диагно-
стического параметра под влиянием внешних эксплуатационных факторов; ЛПу - значение диагностического параметра, численно равное номинальному значению параметра уг0 при t = 0.
Вариация диагностического Луэг(0 по отношению к его номинальному значению под влиянием внешних эксплуатационных факторов обычно значительно меньше, чем его вариация в процессе приработки, поэтому величиной Луэг(0 можно пренебречь.
Значения параметров £сг- и вг определяют экспериментально по результатам обработки статистических данных. При в > 1 или 0 < в < 1 имеет место непрерывная строго монотонно возрастающая или убывающая скорость изменения диагностического параметра.
На основе выражения (1) получена регрессионная модель изменения диагностического параметра устройства гидроусилителя в зависимости от наработки t в виде случайной функции
Остаточный ресурс устройства гидроусилителя равен
У^ ) = Уо ± Vв',
(2)
где у^) - математическое ожидание изменения диагностического параметра при наработке V, у0 - то же диагностического параметра при t = 0.
Оценка качества предсказания, обеспечиваемого полученной моделью, осуществляется по критерию Фишера Б и коэффициенту детерминации В2 [5].
tост tпр
г ' -
Ук - Уо
V ±ВС у
(3)
где - наработка, соответствующая предельному значению диагностического параметра; у'к - значение диагностического параметра в момент контроля
Гидроусилитель рулевого управления включает следующие основные компоненты: гидронасос; регулятор расхода, состоящий из клапана разности давлений, дросселя и предохранительного клапана; гидрораспределитель; гидроцилиндр; фильтр и гидробак.
На рис. 3 представлены регрессионные модели изменения диагностических параметров регулятора расхода и распределителя в функции наработки (тыс. км), построенные по (2). Для построения регрессионных моделей на стенде проводились экспериментальные исследования гидроусилителей, устанавливаемых на троллейбусы АКСМ-101, в ходе которых определялись экспериментальные значения диагностических параметров - расхода на выходе регулятора расхода и расхода внутренних утечек в распределителе. При этом фиксировалась наработка, которая для регуляторов расхода составила от 53 до 281 тыс. км, а для распределителей - от 91 до 382 тыс. км. Объем выборок в обоих случаях составлял 25 устройств. В качестве начальных точек на графиках приняты номинальные значения расходов, соответствующие технической документации.
8,9 г
брр, л/мин 8,7 8,6 8,5 8,4 8,3 8,2 8.1
V
\ ><
N \ V
> X
ж $ ч
\ X
X . X
X
О |Ц|. 2/у:.г, тыс. км
0,6
бВр, л/мин 0,4 0,3 0,2 0,1 0
■ ;<; 1'Ь тыс. км
Рис. 3. Регрессионные модели изменения расхода на выходе регулятора Qрр и расхода внутренних утечек в распределителе Qвр
В табл. 1 приведены значения коэффициентов построенных регрессионных моделей, а также сведения об оценке качества предсказания, обеспечиваемого данными моделями. Так как для регулятора и распределителя полученные значения критерия Фишера Б больше табличного значения критерия Фишера Бт, то регрессионные модели адекватно описывают результаты экспериментов. Значения коэффициента детерминации Я2, определяющего долю рассеяния экспериментальных значений функции отклика, учитываемую регрессионной зависимостью, превышают 0,75, поэтому регрессионные модели работоспособны. По разработанным регрессионным моделям можно проводить прогнозирование изменения технического состояния и остаточного ресурса регулятора расхода и распределителя.
Исчерпывающей характеристикой наработки компонентов гидроусилителя до отказа как случайной величины является вероятностное распределение или функция распределения наработки (модель отказов). Проблема оценки показателей надежности компонентов гидроусилителя сводится к оценке параметров этого распределения.
Строго вероятностные модели отказов (на основе логарифмического, экспоненциального и других законов распределений) [14], устанавливающие связь вероятности появления отказа с наработкой, не связаны с физическими про-
цессами, приводящими к формированию распределения отказов, т. е. не имеют физического обоснования. Вероятностно-физические модели устанавливают связь между вероятностью появления отказа с характеристиками физических процессов деградации, приводящих к отказу. Наибольшими возможностями для расчета надежности систем обладают диффузионное монотонное распределение (ВЫ-
распределение) и диффузионное немонотонное распределение (В^-распределение). Так как физический процесс деградации гидроусилителя описывается случайным процессом с монотонными (необратимыми) реализациями, то распределение отказов гидроусилителя аппроксимируется диффузионным монотонным ВЫ-распределением. Необратимый характер имеют процессы разрушения при усталости, механическом износе, коррозии и старении, присущие механическим объектам.
Плотность вероятности диффузионного монотонного распределения (ВЫ-распределения) определяется по формуле
IЦ) =
? + р
2vt^i2крt
ехр
(' - Р )2
2v V
(4)
где ^ - случайная продолжительность (наработка объекта до отказа, наработка между отказами, ресурс и т. д.); р - параметр масштаба распределения; V - то же формы распределения.
Таблица 1
Параметры регрессионных моделей диагностических параметров
Компонент гидроусилителя Регулятор расхода Распределитель
Наименование диагностического параметра Расход на выходе регулятора Qрр Расход внутренних утечек Qвр
Коэффициенты регрессионной модели диагностического параметра Уо 8,86 0,05
Вс -4,53759 • 10-4 6,38032 • 10-10
в 0,58523 1,58963
Выборочный коэффициент корреляции г -0,918 0,989
Остаточная дисперсия ^т 0,00702 0,00227
Критерий Фишера расчетный Б 4,79663 12,16923
табличный 2,01 2,01
Коэффициент детерминации Я2 0,79986 0,92111
Параметры БМ-распределения определены с использованием моментных оценок.
Моментная оценка параметра масштаба диффузионного распределения
у 1
р = 3
45 - (52 + 3б)Г
(5)
где 5 - выборочное среднее; Б - выборочная дисперсия.
Моментная оценка параметра формы диффузионного распределения
V =
2
№
+ 3Б - 5
45 -л/52
+ 3Б
1/2
(6)
Модель отказов (функция распределения отказов) имеет вид
¥ (0 = БМ (^ р; V) = Ф
t - р
(7)
где Ф(*) - функция нормированного нормального распределения.
Модель надежности (вероятность безотказной работы) имеет вид
Щ) = Ф
р -1
(8)
Модель отказов и модель надежности связаны соотношением ¥ ^) + ) = 1.
Функция нормированного нормального распределения Ф(*) определяется как Ф^*) =
t*
= — [ е 0'5(t ) Л*, где для функции распределил
лени^" отказов ^ =
и t - р
я отказов t* = —
vрt
, а для вероятности
безотказной работы ^ =
р -1
В табл. 2 приведены значения характеристик БМ-распределения отказов регулятора расхода и распределителя.
На рис. 4 показаны графики плотностей вероятности, функции распределения отказов и вероятности безотказной работы регулятора расхода и распределителя в функции наработки t (тыс. км). Как видно из графиков, у регулятора расхода наибольшая плотность вероятности отказов наблюдается при наработке 85 тыс. км, а у регулятора расхода - при 142 тыс. км.
Таблица 2
Характеристики ОМ-распределения отказов регулятора расхода и распределителя
Наименование характеристики БМ-распределения Расчетная формула Регулятор расхода Распределитель
Математическое ожидание м (о = +у ^ 136,63319 • 103 215,06831 • 103
Дисперсия ( 5 2 ^ Б(г) = р V 1 + — 4 V У 6,02186 • 109 1,26994 • 1010
Коэффициент вариации к (Г) = 2 + V2 0,56795 0,52398
Коэффициент асимметрии 4у(б + 11У2 ) ^ = (4 + 5У2 )?/2 1,60315 1,49268
Коэффициент эксцесса Е 6у2 (40 + 93у2 ) Ек = (4 + 5У2 ) 4,15781 3,61907
I ■
Р(Г), Я(г) 0,6 -0,4 ■ 0,2 -О -
\
Р(Г) .
I' I .1 Л .1.. г, тыс. км
■V : V:.. I-:-. г, тыс. км
в
г
Рис. 4. Графики характеристик ВЫ-распределения: а, б - плотность и вероятности отказов регулятора расхода и распределителя соответственно; в, г - модель отказов Б(г) и модель надежности Я(г) регулятора расхода и распределителя соответственно
В Ы В О Д
Таким образом, модели и технические средства, разработанные в соответствии с описанным алгоритмом прогнозирования параметрической надежности, позволяют определять показатели надежности, в том числе остаточный ресурс работы компонентов гидроусилителя рулевого управления. Данный алгоритм апробирован на примере гидроусилителя рулевого управления троллейбусов АКСМ-101.
Л И Т Е Р А Т У Р А
1. Богдан, Н. В. Гидропневмоавтоматика и гидропривод мобильных машин. Эксплуатация и надежность гидро-и пневмосистем: учеб. пособие / Н. В. Богдан, П. Н. Киш-кевич, В. С. Шевченко; под общ. ред. Н. В. Богдана. -Минск: Ураджай, 2001. - 396 с.
2. Проников, А. С. Параметрическая надежность машин / А. С. Проников. - М.: МГТУ им. Н. Э. Баумана, 2002. - 560 с.
3. Гринчар, Н. Г. Методы и средства повышения эксплуатационной надежности гидроприводов дорожных и строительных машин: автореф. дис. ... д-ра техн. наук: 05.05.04 / Н. Г. Гринчар; Московский государственный
университет путей сообщения (МИИТ). - Новочеркасск, 2007. - 35 с.
4. Халамонский, М. Г. Диагностирование гидравлических систем самоходных рабочих машин: автореф. дис. ... д-ра техн. наук: 05.22.10 / М. Г. Халамонский; Белорусская государственная политехническая академия. - Минск, 2000. - 31 с.
5. Тарасик, В. П. Математическое моделирование технических систем: учеб. для вузов / В. П. Тарасик. -Минск: ДизайнПРО, 2004. - 640 с.
6. Тарасик, В. П. Стендовые испытания гидроусилителей рулевых управлений / В. П. Тарасик, Т. В. Мрочек // Материалы, оборудование и ресурсосберегающие технологии: материалы междунар. науч.-техн. конф. - Могилев: Белорус.-Рос. ун-т, 2006. - С. 157.
7. Стенд для общего или поэлементного диагностирования гидроусилителей рулевого управления; пат. 10413 РБ, С1 МПК (2006) О 01М 17/007, В 62Б 5/06 // В. П. Тарасик, В. И. Мрочек, Т. В. Мрочек; заявитель и патентообладатель Бел.-Рос. ун-т. - № а 20060763; заявл. 21.07.2006; опубл. 02.28.2008 // Бюл. - № 1. - 14 с.
8. Техническая диагностика гидравлических приводов / Т. В. Алексеева [и др.]; под общ. ред. Т. М. Башты. -М.: Машиностроение, 1989. - 264 с.
9. Биргер, И. А. Техническая диагностика / И. А. Бир-гер. - М.: Машиностроение, 1978. - 240 с.
10. Сырицын, Т. А. Эксплуатация и надежность гидро- и пневмоприводов: учеб. / Т. А. Сырицын. - М.: Машиностроение, 1990. - 248 с.
11. Михлин, В. М. Прогнозирование технического состояния машин / В. М. Михлин. - М.: Колос, 1976. -288 с.
12. Сухарев, Э. А. Теория эксплуатационной надежности машин. Лекционный курс: учеб. пособие / Э. А. Сухарев. - Рiвне: Удавг, 1997. - 162 с.
13. Надежность в технике. Расчет надежности. Основные положения: ГОСТ 27.301-95. - Минск: Госстандарт Республики Беларусь, 1996. - 16 с.
14. Надежность в технике. Модели отказов. Основные положения: ГОСТ 27.00-97. - Минск: Госстандарт Республики Беларусь, 2005. - 45 с.
Поступила 30.05.08
УДК 621.785.048.7
ПОВЫШЕНИЕ ИЗНОСОСТОЙКОСТИ И РЕСУРСА РАБОТЫ
МНОГОДИСКОВЫХ ТОРМОЗНЫХ МЕХАНИЗМОВ ТРАКТОРА «БЕЛАРУС» ЛАЗЕРНЫМ ТЕРМОУПРОЧНЕНИЕМ БЫСТРОИЗНАШИВАЮЩИХСЯ ДЕТАЛЕЙ
Канд. техн. наук КОБЯКОВ О. С., докт. техн. наук, проф. СПИРИДОНОВ Н. В., инж. БАРКУН А. А.
Белорусский национальный технический университет, Минский тракторный завод
Условия взаимодействия фрикционных механизмов и узлов трения, тормоза трактора «Беларус» модели 1221М, работающих в масляной смазочной среде, характеризуются большим разнообразием одновременно действующих эксплуатационных параметров, включая скоростные режимы относительного скольжения, удельные силовые нагрузки, влияние условий смазки на режимы трения, температурные режимы пар трения.
Одной из проблем является лазерное термоупрочнение внутренних контактных поверхностей корпусных деталей с целью повышения износостойкости и надежности работы тормозных механизмов трактора «Беларус».
Тормозная многодисковая система трактора модели 1221М (рис. 1) отличается повышенной надежностью, долговечностью и отвечает ряду современных требований и международным стандартам качества. Однако, с точки зрения равноресурсности, ряд деталей тормозных механизмов нуждается в повышении ресурса по прочностным и износостойким характеристикам. К числу таких деталей относятся корпусные детали тормоза, в которых повышенному
изнашиванию в процессе эксплуатации подвержены упоры промежуточных дисков тормозных механизмов.
Рис. 1. Сборочный чертеж тормозной системы трактора модели 1221М. Проекция А-А: 1 - корпус тормозной системы; 2 - нажимные диски; 3 - разжимные диски; 4 -промежуточные диски; 5 - шары разжимных дисков; 6 -шарнирный механизм; 7 - возвратные пружины
