RU2547947C1 - Способ диагностики технического состояния роторных агрегатов - Google Patents

Способ диагностики технического состояния роторных агрегатов Download PDF

Info

Publication number
RU2547947C1
RU2547947C1 RU2014100085/11A RU2014100085A RU2547947C1 RU 2547947 C1 RU2547947 C1 RU 2547947C1 RU 2014100085/11 A RU2014100085/11 A RU 2014100085/11A RU 2014100085 A RU2014100085 A RU 2014100085A RU 2547947 C1 RU2547947 C1 RU 2547947C1
Authority
RU
Russia
Prior art keywords
value
values
diagnostic
factor
critical
Prior art date
Application number
RU2014100085/11A
Other languages
English (en)
Inventor
Владимир Николаевич Костюков
Алексей Владимирович Костюков
Денис Викторович Казарин
Андрей Валерьевич Зайцев
Original Assignee
Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика, надежность машин и комплексная автоматизация"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика, надежность машин и комплексная автоматизация" filed Critical Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика, надежность машин и комплексная автоматизация"
Priority to RU2014100085/11A priority Critical patent/RU2547947C1/ru
Application granted granted Critical
Publication of RU2547947C1 publication Critical patent/RU2547947C1/ru

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

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

Description

Изобретение относится к области диагностики, а именно к способам оценки технического состояния роторных агрегатов, и может быть использовано при оценке состояния подшипниковых узлов, например колесно-моторных блоков (КМБ) подвижного состава железнодорожного транспорта.
Известен способ диагностики подшипников качения и редукторов объектов железнодорожного транспорта, реализованный в комплексе оперативной диагностики ПРОГНОЗ-1М [1, 2], включающий измерение вибрации в контрольных точках на подшипниковых щитах узлов, например КМБ, а также частоты вращения колесной пары. При достижении колесной парой установившегося значения частоты вращения не менее 180 об/мин, с допустимой зоной разброса частоты вращения 20%, производят измерение параметров вибрации в выбранной контрольной точке КМБ, а также измерение сопутствующего фактора, например времени фактической работы объекта и среднего значения температуры окружающей среды за время его фактической работы. Одновременно с измерением, прослушивают сигнал с датчика вибрации через наушники. Техническое состояние объекта оценивают путем сравнения полученных параметров вибрации объекта с заранее заданными критическими и оценочными уровнями, которые корректируют с учетом измеренного сопутствующего фактора и субъективных ощущений, а также путем дополнительного анализа данных о предыдущих ремонтных работах и замечаниях по работе узлов. Формируют отчеты по результатам диагностирования всех узлов в текстовом виде, которые редактируют вручную, с учетом дополнительного анализа и прослушивания сигнала.
Недостатками данного способа являются:
- низкая достоверность диагностирования за счет влияния человеческого фактора на постановку диагноза, связанного с возможностью произвольного субъективного изменения критических значений, разделяющих классы состояний «годен-негоден», и редактирования результатов диагностирования;
- низкая достоверность диагностирования ввиду низкой частоты вращения на уровне 180 об/мин, которая недостаточна для обеспечения нагрузки в зонах поверхностей качения подвешенной колесной пары, на которых образуются эксплуатационные дефекты, кроме того, значения вибрации в принятой допустимой 20% зоне разброса частоты вращения могут значительно отличаться;
- высокая трудоемкость применения способа, связанная с необходимостью дополнительного прослушивания сигнала с датчика вибрации и проведения анализа данных о предыдущих ремонтных работах и замечаниях по работе узлов.
Известен способ диагностики колесно-моторных блоков мотор-вагонного подвижного состава, реализованный в системе [3], включающий одновременную установку датчиков вибрации на диагностируемые узлы, измерение параметров вибрации на заданном скоростном режиме с высокой точностью поддержания стабильной частоты вращения, и оценку технического состояния путем сравнения их значений с критическими.
В качестве наиболее близкого аналога, принятого нами за прототип, является способ оценки технического состояния роторного агрегата [4], на примере центробежного насосного агрегата, по вибрации корпуса, которую измеряют в процессе эксплуатации по совокупности входящих в агрегат элементов: ротора двигателя, исполнительной машины, опорных подшипниковых узлов, соединительной муфты, остова и фундамента, к которому крепится агрегат, в различных частотных полосах, например высокочастотной, среднечастотной и низкочастотной, соответствующих виброускорению, виброскорости и виброперемещению элементов агрегата, определяют эти вибропараметры и другие информативные составляющие сигнала вибрации, используют их в качестве диагностических признаков (ДП), соответствующих совокупности входящих в агрегат элементов, оценку технического состояния производят комплексно по табличной зависимости путем сравнения полученных значений с критическими границами ДП, например аварийными и предупредительными, оповещают персонал о недопустимом состоянии агрегата, при этом критические границы ДП определяют эмпирическим путем. При этом значение вибропараметра ниже предупредительного порога говорит о том, что состояние элемента или агрегата нормальное. Значение вибропараметра между аварийным и предупредительным порогами указывает на критическое состояние. Значение параметра, превышающее аварийный порог, требует немедленной остановки агрегата.
Недостатками данного способа являются:
- недостаточная достоверность диагностирования слабых и зарождающихся дефектов ввиду нестабильного их проявления в сигналах вибрации на низкой частоте вращения (180-210 об/мин), которую выбирают исходя из требований экономии электроэнергии и обеспечения безопасности при диагностировании;
- недостаточная достоверность диагностирования оборудования при единичном диагностировании вследствие отсутствия решения о состоянии объекта при попадании измеренного значения в зону между аварийным и предупредительным порогами.
Кроме того, общим недостатком перечисленных способов является недопустимо высокая ошибка оценки технического состояния объекта, связанная со статистической ошибкой оценки состояния объекта при единичных испытаниях, ограничивающая достоверность диагностирования на уровне менее 90%.
Существующие способы диагностики узлов и агрегатов направлены на разделение их технического состояния на классы - исправное и неисправное. Считается, что приемлемой является ошибка диагностирования на уровне не более 1%, однако при единичных испытаниях получить подобный результат зачастую невозможно, что требует проведения дополнительного диагностирования.
Целью предлагаемого способа является повышение достоверности диагностирования технического состояния роторных агрегатов при одновременном снижении трудоемкости его применения.
Поставленная цель в способе диагностики технического состояния роторных агрегатов достигается тем, что перед измерением параметров вибрации узлов роторных агрегатов, используемых в качестве диагностических признаков, задают величину вероятности ошибки ложной тревоги и пропуска дефекта по каждому диагностическому признаку, для которых устанавливают зону неопределенности принятия решения о техническом состоянии соответствующего узла роторного агрегата, ограничивают ее нижним и верхним критическими значениями, устанавливают критическое значение для другого значения сопутствующего фактора, например для повышенной частоты вращения, определяют значения диагностических признаков при первом значении сопутствующего фактора, сравнивают их с критическими значениями, если значения диагностических признаков меньше соответствующих нижних критических значений, то данные узлы считают исправными, если значения диагностических признаков больше верхних критических значений, то данные узлы считают дефектными, если значения диагностических признаков оказываются между нижним и верхним критическими значениями, их помещают в зону неопределенности, в которой ошибка диагностирования превышает заданную величину, определяют диагностические признаки при другом значении сопутствующего фактора, например, при повышенной частоте вращения, сравнивают их с установленным критическим значением для соответствующего значения сопутствующего фактора, в случае превышения которого узел считают дефектным, тем самым диагностику осуществляют в соответствии с таблицей, при этом ошибка диагностирования не превышает произведения ошибок диагностирования при каждом значении сопутствующего фактора:
Figure 00000001
j, k - значения сопутствующего фактора, например, j=n1 - базовая частота вращения, k=n2 - повышенная частота вращения;
Yj - значения диагностического признака при j-м значении сопутствующего фактора;
Yk - значения диагностического признака при k-м значении сопутствующего фактора;
Figure 00000002
- нижнее критическое значение диагностического признака при j-м значении сопутствующего фактора;
Figure 00000003
- верхнее критическое значение диагностического признака при j-м значении сопутствующего фактора;
Figure 00000004
- критическое значение диагностического признака при k-м значении сопутствующего фактора;
S - состояние узла, где:
S=1 - узел исправен;
S=0 - узел неисправен;
S≠1 & S≠0 - состояние узла не определено.
Анализ отличительных признаков предлагаемого способа диагностики технического состояния роторных агрегатов показал, что:
- установка на каждом диагностическом признаке зоны неопределенности, ограниченной верхним и нижним критическими значениями, позволяет выявить скрытые и зарождающиеся дефекты;
- определение диагностических признаков при других значениях сопутствующего фактора, например при повышенной частоте вращения, в случае попадания в зону неопределенности, обеспечивает снижение вероятности ошибки диагностирования;
- диагностирование на основе предложенной диагностической таблицы обеспечивает высокую достоверность и простоту получения результатов.
Таким образом, предложенная совокупность отличительных признаков, обеспечивающая полученный результат, представляется новой на данном этапе развития науки и техники, превосходит существующий мировой уровень. Изобретение соответствует изобретательскому уровню, поскольку достигаемый результат определяется не только совокупностью отличительных признаков, но и результатом их тесного взаимодействия между собой.
Суть способа поясняется Фиг.1 и 2:
На фиг.1 представлены зависимости значений диагностического признака (ДП) от величины сопутствующего фактора.
На фиг.2 представлен графический расчет вероятности ошибки диагностирования при двух значениях сопутствующего фактора, например частоты вращения.
На Фиг.1 показаны:
1…10 - значения ДП для десяти исправных буксовых узлов, при первом значении сопутствующего фактора n1;
11…20 - значения ДП для десяти дефектных буксовых узлов, при первом значении сопутствующего фактора n1;
21…30 - значения ДП для десяти исправных буксовых узлов, при втором значении сопутствующего фактора n2;
31…40-значения ДП для десяти дефектных буксовых узлов, при втором значении сопутствующего фактора n2;
41 - нижний «хвост» функции распределения ДП
Figure 00000005
для дефектных буксовых узлов, при первом значении сопутствующего фактора n1;
42 - нижний «хвост» обратной функции распределения ДП
Figure 00000006
для исправных буксовых узлов, при первом значении сопутствующего фактора n1;
43 - точка пересечения функции распределения значений ДП для дефектных узлов с обратной функцией распределения для исправных узлов, при первом значении сопутствующего фактора n1;
44 - проекция точки пересечения 43 на ось ДП;
45 - критическое значение ДП
Figure 00000007
при первом значении сопутствующего фактора n1;
48 - нижний «хвост» функции распределения ДП
Figure 00000008
для дефектных буксовых узлов, при втором значении сопутствующего фактора n2;
49 - нижний «хвост» обратной функции распределения ДП
Figure 00000009
для исправных буксовых узлов, при втором значении сопутствующего фактора n2;
50 - точка пересечения функции распределения значений ДП для дефектных узлов с обратной функцией распределения для исправных узлов, при втором значении сопутствующего фактора n2;
51 - проекция точки пересечения 50 на ось ДП;
52 - критическое значение ДП
Figure 00000010
при втором значении сопутствующего фактора n2;
58 - проекция точки пересечения 55 на ось ДП;
59 - нижнее критическое значение
Figure 00000011
при первом значении сопутствующего фактора n1;
60 - проекция точки пересечения 56 на ось ДП;
61 - верхнее критическое значение ДП
Figure 00000012
при первом значении сопутствующего фактора n1;
На Фиг.2 показаны:
41 - нижний «хвост» функции распределения ДП
Figure 00000013
для дефектных буксовых узлов, при первом значении сопутствующего фактора n1;
42 - нижний «хвост» обратной функции распределения ДП
Figure 00000014
для исправных буксовых узлов, при первом значении сопутствующего фактора n1;
43 - точка пересечения функции распределения значений ДП для дефектных узлов с обратной функцией распределения для исправных узлов, при первом значении сопутствующего фактора n1;
44 - проекция точки пересечения 43 на ось ДП;
45 - критическое значение ДП
Figure 00000015
при первом значении сопутствующего фактора n1;
46 - проекция точки пересечения 43 на ось вероятности ошибки;
47 - величина вероятности равных ошибок ложной тревоги и пропуска дефекта q1;
48 - нижний «хвост» функции распределения ДП
Figure 00000016
для дефектных буксовых узлов, при втором значении сопутствующего фактора n2;
49 - нижний «хвост» обратной функции распределения ДП
Figure 00000017
для исправных буксовых узлов, при втором значении сопутствующего фактора n2;
50 - точка пересечения функции распределения значений ДП для дефектных узлов с обратной функцией распределения для исправных узлов, при втором значении сопутствующего фактора n2;
51 - проекция точки пересечения 50 на ось ДП;
52 - критическое значение ДП
Figure 00000018
при втором значении сопутствующего фактора n2;
53 - проекция точки пересечения 50 на ось вероятности ошибки;
54 - величина вероятности равных ошибок ложной тревоги и пропуска дефекта q2;
55 - точка пересечения проекции 53 с «хвостом» 41, при первом значении сопутствующего фактора n1;
56 - точка пересечения проекции 53 с «хвостом» 42, при первом значении сопутствующего фактора n1;
57 - линия равных ошибок, при величине вероятности ошибки q1=q2;
58 - проекция точки пересечения 55 на ось ДП;
59 - нижнее критическое значение
Figure 00000019
при первом значении сопутствующего фактора n1;
60 - проекция точки пересечения 56 на ось ДП;
61 - верхнее критическое значение ДП
Figure 00000020
при первом значении сопутствующего фактора n1;
62 - величина вероятности равных ошибок ложной тревоги и пропуска дефекта
Figure 00000021
на уровне 5%;
63 - проекция величины вероятности ошибки
Figure 00000021
62 на «хвосты» распределений 41 и 42;
64 - точка пересечения проекции 63 с «хвостом» 41, при первом значении сопутствующего фактора n1;
65 - точка пересечения проекции 63 с «хвостом» 42, при первом значении сопутствующего фактора n1;
66 - линия равных ошибок ложной тревоги и пропуска дефекта, при величине вероятности ошибки
Figure 00000021
;
67 - проекция точки пересечения 64 на ось ДП;
68 - нижнее критическое значение
Figure 00000022
для величины вероятности равных ошибок ложной тревоги и пропуска дефекта
Figure 00000021
на уровне 5% при первом значении сопутствующего фактора n1;
69 - проекция точки пересечения 65 на ось ДП;
70 - верхнее критическое значение
Figure 00000023
для величины вероятности равных ошибок ложной тревоги и пропуска дефекта
Figure 00000021
на уровне 5% при первом значении сопутствующего фактора n1;
71 - линия неравных ошибок, когда ошибка пропуска дефекта 54 меньше ошибки ложной тревоги 62, соединяющая точки 55 и 65;
72 - линия неравных ошибок, когда ошибка пропуска дефекта 62 больше ошибки ложной тревоги 54, соединяющая точки 64 и 56.
В основе изобретения лежит экспериментально установленный факт более резкого увеличения значений вибропараметров, например виброускорения, с ростом частоты вращения для дефектных узлов по сравнению с исправными.
Реализуемость способа показана при вибродиагностике подшипниковых узлов роторных агрегатов. В качестве роторного агрегата рассмотрены буксовые узлы КМБ. Исследовали десять буксовых узлов КМБ с исправными подшипниками и десять буксовых узлов КМБ с дефектными подшипниками. Измерение ДП, например виброускорения, проводили для двух значений сопутствующего фактора, например частоты вращения колесной пары. В декартовой системе координат строят зависимость ДП от величины сопутствующего фактора, где по оси абсцисс откладывают значения сопутствующего фактора, а по оси ординат - значения ДП (Фиг.1).
При первой величине сопутствующего фактора, например частоты вращения колесной пары n1=210 об/мин, для десяти исправных буксовых узлов получены значения ДП 1…10, а для десяти дефектных буксовых узлов получены значения ДП 11…20 (Фиг.1). Аналогичным образом получены значения ДП при второй величине сопутствующего фактора, например частоты вращения колесной пары n2=420 об/мин, для исправных буксовых узлов 21…30 и для дефектных буксовых узлов 31…40.
Строят график зависимости величины вероятности ошибки диагностики от критического значения ДП. Для первого значения сопутствующего фактора пересечение функции распределения ДП
Figure 00000024
41 для дефектных узлов с обратной функцией распределения ДП
Figure 00000025
42 для исправных узлов дает точку пересечения 43, проекция 44 на ось ДП соответствует критическому значению ДП
Figure 00000026
45, проекция 46 на ось вероятности ошибки дает величину вероятности ошибки q1 47 (Фиг.2).
Для второго значения сопутствующего фактора пересечение функции распределения ДП
Figure 00000027
48 для дефектных узлов с обратной функцией распределения ДП
Figure 00000028
49 для исправных узлов дает точку пересечения 50, проекция 51 на ось ДП соответствует критическому значению ДП
Figure 00000029
52, проекция 53 на ось вероятности ошибки дает величину вероятности ошибки q2 54 (Фиг.2).
Задают величину вероятности ошибки ложной тревоги и пропуска дефекта и устанавливают зону неопределенности принятия решения о техническом состоянии узла.
Для простоты оценки технического состояния агрегата, ошибки ложной тревоги и пропуска дефекта выбираются равными.
Например, для диагностирования при первом значении сопутствующего фактора задают величину вероятности ошибки q1 равной вероятности ошибки q2 54 при втором значении сопутствующего фактора (Фиг.2). Проекция вероятности ошибки 54 на нижние «хвосты» 41 и 42 дает соответственно точки пересечения 55 и 56, соединив которые получают линию равных ошибок 57, ее проекция на ось значений ДП дает величину зоны неопределенности при равной ошибке ложной отбраковки и пропуска дефекта. Проекция 58 точки 55 на ось значений ДП дает нижнее критическое значение
Figure 00000030
59. Проекция 60 точки 56 на ось значений ДП дает верхнее критическое значение
Figure 00000031
61.
Например, задают величину вероятности ошибки
Figure 00000021
62 на уровне 5%, проекция вероятности ошибки 63 на нижние «хвосты» 41 и 42 дает соответственно точки пересечения 64 и 65, соединив которые получают линию равных ошибок 66, ее проекция на ось значений ДП дает величину зоны неопределенности при равной ошибке ложной отбраковки и пропуска дефекта. Проекция 67 точки 64 на ось значений ДП дает нижнее критическое значение
Figure 00000032
68. Проекция 69 точки 65 на ось значений ДП дает верхнее критическое значение
Figure 00000033
70 (Фиг.2).
В случае попадания измеренного значения ДП в зону неопределенности проводят дополнительное испытание при другом значении сопутствующего фактора.
Использование нижнего критического значения ДП
Figure 00000030
при меньшей вероятности пропуска дефекта, например на уровне 3%, дает большее количество узлов, которые требуют дополнительного диагностирования, и напротив, использование нижнего критического значения
Figure 00000030
при большей вероятности пропуска дефекта, например на уровне 5%, дает меньшее количество узлов которые требуют дополнительного диагностирования.
Использование верхнего критического значения ДП
Figure 00000034
при меньшей вероятности ложной тревоги, например на уровне 3%, дает большее количество узлов, которые требуют дополнительного диагностирования, и напротив, использование верхнего критического значения ДП
Figure 00000034
при большей вероятности ложной тревоги, например на уровне 5%, дает меньшее количество узлов, которые требуют дополнительного диагностирования.
В случае неравных ошибок диагностирования, например низкой ошибки пропуска дефекта и высокой ошибки ложной тревоги, зона неопределенности соответствует линии неравных ошибок 71, а в случае высокой ошибки пропуска дефекта и низкой ошибки ложной тревоги - линии неравных ошибок 72 (Фиг.2).
Результирующая ошибка диагностирования q представляет собой произведение ошибок первого q1 и второго q2 диагностирования q=q1·q2.
Таким образом, взаимодействие использования зоны неопределенности при первом значении сопутствующего фактора и дополнительного испытания при другом значении сопутствующего фактора позволяет повысить достоверность диагностирования технического состояния роторных агрегатов при одновременном снижении трудоемкости его применения.
В качестве примера исполнения (реализации) способа рассмотрен процесс диагностирования буксовых узлов КМБ с помощью системы [3] с введенными в нее установленными критическими значениями ДП.
Для диагностирования букс колесно-моторных блоков электропоезда при первом значении сопутствующего фактора, например частоты вращения колесной пары n1=210 об/мин, величину вероятности ошибки ложной тревоги qЛТ и пропуска дефекта qПД задают равными друг другу, т.е. qЛТ=qПД=q1. Величину вероятности ошибки q1 задают, например, равной величине вероятности ошибки q2 54, при другом значении сопутствующего фактора, соответствующем, например, частоте вращения колесной пары n2=420 об/мин, т.е. q1=q2=3,2% (Фиг.2). В соответствии с выбранной величиной вероятности ошибки получают зону неопределенности принятия решения о техническом состоянии буксы, ограниченную нижним
Figure 00000035
59 и верхним
Figure 00000036
61 критическими значениями, а для другого значения сопутствующего фактора, например частоты вращения колесной пары n2=420 об/мин, устанавливают критическое значение
Figure 00000037
52.
При диагностировании двадцати буксовых узлов при первом значении сопутствующего фактора, например частоты вращения колесной пары n1=210 об/мин, получили значения ДП (виброускорения) 1…20 (Фиг.1). При сравнении полученных величин с критическими значениями выявили:
- значения ДП
Figure 00000038
(1,47…3,7 м/с2) для первого - девятого буксовых узлов 1…9 меньше соответствующего нижнего критического значения
Figure 00000039
, соответственно техническое состояние данных узлов распознают как исправное (Фиг.1) при первом испытании (n1=210 об/мин);
- значения ДП
Figure 00000040
(5,2…9,25 м/с2) для двенадцатого - двадцатого буксовых узлов 12…20 больше соответствующего верхнего критического значения
Figure 00000041
, соответственно техническое состояние данных узлов распознают как неисправное (Фиг.1) при первом испытании (n1=210 об/мин);
- значения ДП
Figure 00000042
и
Figure 00000043
для десятого и одиннадцатого буксовых узлов 10 и 11 находятся между нижним
Figure 00000044
и верхним
Figure 00000045
критическими значениями, т.е. помещаются в зоне неопределенности, и проводят дополнительное испытание при другом значении сопутствующего фактора, например частоты вращения колесной пары n2=420 об/мин. При этом получают соответствующие значения ДП (виброускорения) 30 и 31 (Фиг.1). При сравнении полученных величин с соответствующим критическим значением
Figure 00000046
52 выявили:
- значение ДП
Figure 00000047
для десятого буксового узла 30 меньше соответствующего критического значения
Figure 00000048
, соответственно техническое состояние данного узла распознают как исправное (Фиг.1);
- значение ДП
Figure 00000049
для одиннадцатого буксового узла 31 больше соответствующего критического значения
Figure 00000048
, соответственно техническое состояние данного узла распознают как неисправное (Фиг.1).
Кроме того, при первом значении сопутствующего фактора величины вероятности ошибки ложной тревоги и пропуска дефекта могут задавать не равными друг другу qЛТ≠qПД, при этом существует два варианта задания вероятности ошибки ложной тревоги и пропуска дефекта: первый вариант - вероятность ошибки ложной тревоги меньше вероятности ошибки пропуска дефекта, второй - вероятность ошибки ложной тревоги больше вероятности ошибки пропуска дефекта.
Вариант 1. Величину вероятности ошибки ложной тревоги задают равной qЛТ=3,2% 54, а величину вероятности ошибки пропуска дефекта задают равной qПД=5% 62 (Фиг.2). В соответствии с выбранными величинами вероятности ошибки получают зону неопределенности принятия решения о техническом состоянии буксы, ограниченную нижним
Figure 00000050
59 и верхним
Figure 00000051
62 критическими значениями, а для другого значения сопутствующего фактора, например частоты вращения колесной пары n2=420 об/мин, устанавливают критическое значение
Figure 00000048
52.
При диагностировании двадцати буксовых узлов при первом значении сопутствующего фактора, например частоты вращения колесной пары n1=210 об/мин, получили значения ДП (виброускорения) 1…20 (Фиг.1). При сравнении полученных величин с критическими значениями выявили:
- значения ДП
Figure 00000038
(1,47…3,7 м/с2) для первого - девятого буксовых узлов 1…9 меньше соответствующего нижнего критического значения
Figure 00000039
, соответственно техническое состояние данных узлов распознают как исправное (Фиг.1) при первом испытании (n1=210 об/мин);
- значения ДП
Figure 00000052
,
Figure 00000040
(4,7; 5,2…9,25 м/с2) для десятого, двенадцатого…двадцатого буксовых узлов 10; 12…20 больше соответствующего верхнего критического значения
Figure 00000051
, соответственно техническое состояние данных узлов распознают как неисправное (Фиг.1) при первом испытании (n1=210 об/мин);
- значение ДП
Figure 00000043
для одиннадцатого буксового узла 11 находится между нижним
Figure 00000050
и верхним
Figure 00000051
критическими значениями, т.е. помещаются в зоне неопределенности, и проводят дополнительное испытание при другом значении сопутствующего фактора, например частоты вращения колесной пары n2=420 об/мин. При этом получают соответствующее значение ДП (виброускорения) 31 (Фиг.1). При сравнении полученной величины с соответствующим критическим значением
Figure 00000048
52 выявили:
- значение ДП
Figure 00000049
для одиннадцатого буксового узла 31 больше соответствующего критического значения
Figure 00000048
, соответственно техническое состояние данного узла распознают как неисправное (Фиг.1).
Вариант 2. Величину вероятности ошибки ложной тревоги задают равной qЛТ=5% 62, а величину вероятности ошибки пропуска дефекта задают равной qПД=3,2% 54 (Фиг.2). В соответствии с выбранными величинами вероятности ошибки получают зону неопределенности принятия решения о техническом состоянии буксы, ограниченную нижним
Figure 00000053
68 и верхним
Figure 00000054
61 критическими значениями, а для другого значения сопутствующего фактора, например частоты вращения колесной пары n2=420 об/мин, устанавливают критическое значение
Figure 00000048
52.
При диагностировании двадцати буксовых узлов при первом значении сопутствующего фактора, например частоты вращения колесной пары n1=210 об/мин, получили значения ДП (виброускорения) 1…20 (Фиг.1). При сравнении полученных величин с критическими значениями выявили:
- значения ДП
Figure 00000038
(1,47…3,7 м/с2) для первого - девятого буксовых узлов 1…9 меньше соответствующего нижнего критического значения
Figure 00000055
, соответственно техническое состояние данных узлов распознают как исправное (Фиг.1) при первом испытании (n1=210 об/мин);
- значения ДП
Figure 00000040
(5,2…9,25 м/с2) для двенадцатого - двадцатого буксовых узлов 12…20 больше соответствующего верхнего критического значения
Figure 00000056
, соответственно техническое состояние данных узлов распознают как неисправное (Фиг.1) при первом испытании (n1=210 об/мин);
- значения ПД
Figure 00000057
и
Figure 00000058
для десятого и одиннадцатого буксовых узлов 10 и 11 находятся между нижним
Figure 00000055
и верхним
Figure 00000056
критическими значениями, т.е. помещаются в зоне неопределенности, и проводят дополнительное испытание при другом значении сопутствующего фактора, например частоты вращения колесной пары n2=420 об/мин. При этом получают соответствующие значения ДП (виброускорения) 30 и 31 (Фиг.1). При сравнении полученных величин с соответствующим критическим значением
Figure 00000048
52 выявили:
- значение ДП
Figure 00000059
для десятого буксового узла 30 меньше соответствующего критического значения
Figure 00000048
, соответственно техническое состояние данного узла распознают как исправное (Фиг.1);
- значение ДП для одиннадцатого буксового узла 31
Figure 00000060
больше соответствующего критического значения
Figure 00000048
, соответственно техническое состояние данного узла распознают как неисправное (Фиг.1).
Таким образом, использование зоны неопределенности при оценке технического состояния, а также двух независимых испытаний, дающих ошибку диагностирования равной произведению ошибок первого и второго испытания, обеспечивает достижение поставленной цели - повышение достоверности диагностирования технического состояния роторных агрегатов при одновременном снижении трудоемкости его применения.
Предложенный способ диагностики реализован в системах КОМПАКС-ЭКСПРЕСС и успешно прошел испытания в ряде мотор-вагонных депо ОАО «РЖД».
Литература
1. Руководство по эксплуатации на объектах железнодорожного транспорта 12ДК.318558.005 РЭ. «Комплекс оперативной диагностики ПРОГНОЗ-1М». - с.7-19, 68-77. http://vibraomsk.ru/prod_one.php?id=6&pic=0.
2. Пат. РФ №2284021 С2, G01M 13/00; Заявлено от 04.11.2004.
3. Пат. РФ №2378633 C1, G01M 17/08; Заявлено от 26.09.2008.
4. Пат. РФ 2068553, G01M 15/00; Заявлено от 27.10.1996.

Claims (1)

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

    где: j, k - значения сопутствующего фактора, например, j=n1 - базовая частота вращения, k=n2 - повышенная частота вращения;
    Yj - значения диагностического признака, при j-ом значении сопутствующего фактора;
    Yk - значения диагностического признака, при k-м значении сопутствующего фактора;
    Figure 00000062
    - нижнее критическое значение диагностического признака, при j-м значении сопутствующего фактора;
    Figure 00000063
    - верхнее критическое значение диагностического признака, при j-м значении сопутствующего фактора;
    Figure 00000064
    - критическое значение диагностического признака, при k-м значении сопутствующего фактора;
    S - состояние узла, где:
    S=1 - узел исправен;
    S=0 - узел неисправен;
    S≠1 & S≠0 - состояние узла не определено.
RU2014100085/11A 2014-03-03 2014-03-03 Способ диагностики технического состояния роторных агрегатов RU2547947C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2014100085/11A RU2547947C1 (ru) 2014-03-03 2014-03-03 Способ диагностики технического состояния роторных агрегатов

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2014100085/11A RU2547947C1 (ru) 2014-03-03 2014-03-03 Способ диагностики технического состояния роторных агрегатов

Publications (1)

Publication Number Publication Date
RU2547947C1 true RU2547947C1 (ru) 2015-04-10

Family

ID=53296551

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2014100085/11A RU2547947C1 (ru) 2014-03-03 2014-03-03 Способ диагностики технического состояния роторных агрегатов

Country Status (1)

Country Link
RU (1) RU2547947C1 (ru)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2646207C1 (ru) * 2016-11-23 2018-03-01 Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика", надежность машин и комплексная автоматизация" Способ диагностического мониторинга роторных механизмов
RU2658233C1 (ru) * 2017-09-12 2018-06-19 Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика, надежность машин и комплексная автоматизация" Способ диагностики технического состояния роторных механизмов

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2068553C1 (ru) * 1994-08-29 1996-10-27 Костюков Владимир Николаевич Способ оценки технического состояния центробежного насосного агрегата по вибрации корпуса
RU2371695C1 (ru) * 2008-04-28 2009-10-27 Общество С Ограниченной Ответственностью "Диамех 2000" Система вибрационного контроля
RU96427U1 (ru) * 2010-05-18 2010-07-27 Общество с ограниченной ответственностью "ТМК Инновация" Система мониторинга роторных агрегатов

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2068553C1 (ru) * 1994-08-29 1996-10-27 Костюков Владимир Николаевич Способ оценки технического состояния центробежного насосного агрегата по вибрации корпуса
RU2371695C1 (ru) * 2008-04-28 2009-10-27 Общество С Ограниченной Ответственностью "Диамех 2000" Система вибрационного контроля
RU96427U1 (ru) * 2010-05-18 2010-07-27 Общество с ограниченной ответственностью "ТМК Инновация" Система мониторинга роторных агрегатов

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2646207C1 (ru) * 2016-11-23 2018-03-01 Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика", надежность машин и комплексная автоматизация" Способ диагностического мониторинга роторных механизмов
RU2658233C1 (ru) * 2017-09-12 2018-06-19 Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика, надежность машин и комплексная автоматизация" Способ диагностики технического состояния роторных механизмов

Similar Documents

Publication Publication Date Title
KR101829645B1 (ko) 철도 차량 모니터링 장치 및 이를 이용한 모니터링 방법
US11105712B2 (en) Integrated vibe/ODM fusion and trending analysis for prognostic health management of engine bearing
CN104990709B (zh) 用于检测机车轴承故障的方法
EP2436574B1 (en) State monitoring apparatus and state monitoring method of railway car, and railway car
CN107202907B (zh) 检测车辆安全装置的传感器故障的方法
KR101874472B1 (ko) 진동신호의 주파수 에너지를 이용한 회전체 고장 예측 시스템 및 방법
JP6882397B2 (ja) 異音等検出システム、装置、方法及びプログラム
DE112017004678T5 (de) Überwachungsverfahren und Überwachungsvorrichtung für die Restlebensdauer eines Lagers
KR101539896B1 (ko) 유도전동기 오류 진단 방법
CN106922171B (zh) 用于确定结构状态的方法、***及具有其的目标单元
Melnik et al. Rail vehicle suspension condition monitoring-approach and implementation
JP2012058171A (ja) 移動体異常検知システムおよび移動体
RU2547947C1 (ru) Способ диагностики технического состояния роторных агрегатов
CN115790804B (zh) 动态汽车衡状态监测方法、介质、设备及装置
CN112799380A (zh) 适用于无人驾驶列车的辅助***自检***及方法
RU2551447C1 (ru) Способ вибрационной диагностики технического состояния подшипниковой опоры ротора двухвального газотурбинного двигателя
Nowakowski et al. Diagnostics of the drive shaft bearing based on vibrations in the high-frequency range as a part of the vehicle's self-diagnostic system
Ma et al. A prediction method based on stepwise regression analysis for train axle temperature
RU2658233C1 (ru) Способ диагностики технического состояния роторных механизмов
DE102013105397B4 (de) Zustandsüberwachung eines Schienenfahrzeugs
JP2005283227A (ja) 異音検査方法およびその装置
CN114179858B (zh) 基于车轮健康状态的车轮镟修方法及装置
CN115931399A (zh) 一种实时在线检测高速列车轮多边形故障的方法
CN112330615B (zh) 一种转动部件高强度螺栓状态监测方法及***
Kostrzewski Analysis of selected vibroacoustic signals recorded on EMU vehicle running on chosen routes under supervised operating conditions