Технические характеристики виброанализатора СД -21

СОДЕРЖАНИЕ

Введение………………………………………………………………………

1. Задачи технической диагностики…………………………………………..

2. Основные положения технической диагностики…………………………

3. Классификация методов диагностики энергетического оборудования…

3.1 Этапы диагностического обеспечения………………………………..

3.2. Общая методика получения диагностического обеспечения…….

3.3. Прогнозирующее диагностирование ………………………

3.4. Основные показатели средств технической диагностики………

3.5.Показатели эффективности процесса диагностирования………

3.6.Оценка эффективности процесса прогнозирования……………

3.7.Диагностические модели с применением теории подобия………

3.8 Классификация диагностических моделей…………………

3.9. Вероятностные диагностические модели……………………

3.10. Диагностическая модель разладки……………………………

3.11. Модели прогнозирования технического состояния оборудования..

3.12. Структура и задачи основных систем отраслевой системы

диагностирования………………………………………………

3.13. Организация локальной системы диагностирования

котлоагрегатов………………………………………………………………

3.14.Диагностические модели поверхностей нагрева…………………

3.15. Автоматизированная система температурного контроля металла

моноблока……………………………………………………………………

3.16. Автоматизированная система учета расхода металла……………

4. Диагностирование роторного и электротехнического оборудования……

4.1. Методы и средства диагностирования роторных механизмов……

4.2. Методы диагностирования турбин………………………………

4.3. Тепловизионный метод……………………………………………

4.4. Автоматизированные диагностические системы электронного

управляющего оборудования…………………………………………

4.5. Особенности диагностирования устройств релейной защиты и

автоматики…………………………………………………………………..

4.6. Методы проверки датчиков преобразования и аналоговых блоков….

4.7. Организация диагностирования в АСУ ТП…………………………

4.8 Диагностирование водно-химических режимов……………………

4.9. Методы прогнозирующего диагностирования для ВХР.......................

5. Основы технической диагностики……………………………………….

5.1. Основной принцип технической диагностики………………………

5.2. Термины и определения……………………………………………

5.3. Разделы технической диагностики…………………………………

5.4. Основные этапы технической диагностики…………………………

5.5. Функциональная и тестовая диагностика………………………………

5.6. Методология технической диагностики……………………………

5.7. Выбор диагностического сигнала………………………………………

6. Общие понятия об акустическом шуме и вибрации……………………

6.1. Акустический шум……………………………………………………….

6.2. Волновое уравнение……………………………………………………..

6.3. Особенности восприятия акустического шума человеком…………….

6.4. Вибрация…………………………………………………………………..

6.5. Единица измерения шума и вибрации…………………………………..

6.6. Анализ вибродиагностических сигналов………………………………..

6.7. Анализ уровней вибрации………………………………………………..

6.8. Анализ ударных импульсов………………………………………………

6.9. Частотный анализ………………………………………………………….

6.10. Частотный анализ с постоянной абсолютной шириной полосы

пропускания..………………………………………………………………

6.11. Частотный анализ с постоянной относительной шириной

пропускания…………………………………………………………………

6.12. Особенности использования частотного анализа в основных методах

мониторинга и диагностики…………………………………………………

7. Средства измерения и анализа виброакустических сигналов……………..

7.1. Измерительные преобразователи……………………………………

7.2. Измерительные преобразователи шума – микрофоны………………

7.3. Измерительные преобразователи вибрации…………………………

7.4. Измерительные преобразователи вибросмещения………………

7.5. Вихретоковые преобразователи……………………………………

7.6. Измерительные преобразователи виброскорости…………………

7.7. Измерительные преобразователи виброускорения…………………

7.8.Анализаторы…………………………………………………………

7.9. Программное обеспечение для вибромониторинга и диагностики……

7.10. Виброакустические диагностические модели узлов и машин………

7.11. Последовательность процедур при построении диагностической

модели………………………………………………………………………

7.12. Основные частоты вибрации подшипников качения………………

7.13. Определения пороговых значений………………………………

7.14. Концепция системы диагностики конденсационной установки

паровой турбины…………………………………………………………

7.15. Диагностика турбогенераторов………………………………

7.16. Диагностика котлоагрегатов……………………………………..

7.17. Диагностика вспомогательного оборудования………………….

8. Основные направления совершенствования организации ремонтов и

эксплуатации оборудования………………………………………………..

8.1. Этапы ввода системы мониторинга и вибродиагностики

технического состояния оборудования (МиВД ТСО)…………………….

8.2. Заключение……………………………………………………………

Литература…………………………………………………………………………

 

1. Коэффициент готовности:

К_г = [Т_р/(Т_р + Т_н.п. )]·100%,

где Т_р- суммарное время работы в часах.;

Т _н.п. - суммарное время непланового (аварийного) простоя, в часах..

2. Коэффициент аварийности:

К_ав = [Т _н.п. (_ав)/(Т _н.п. (_ав)+Т_р)]·100%; или К_ав=(1-К_г)·100%.

3. Коэффициент оперативной готовности:

К_ог = [(Т_р+Т_рез)/Т _к]·100 %,

где Т_рез- cуммарное время простоя оборудования в резерве по диспетчерскому графику в часах;

Т _к - рассматриваемый календарный период (1 год =8760 ч.)

4. Коэффициент рабочего времени:

К_р = [Т_р/Т_к]·100%.

5. Коэффициент неплановых простоев (недоиспользование календарного времени из-за неплановых простоев):

К _н.п. = [Т _н.п. /Т_к]·100%.

6. Коэффициент плановых простоев:

К_пл = [Т_пл/Т_к]·100%,

где Т_пл – время использования установленной мощности.

7. Коэффициент использования установленной мощности:

К_уст = [Э/N_ном·Т_к]·100%,

где Э-выработанная электрическая энергия;

Nном – номинальная мощность.

8. Наработка на отказ, в часах.:

Т_о = Т_р/n

где n-количество отказов.

9. Параметр потока отказов (определяется в расчете на 1000 или 10000 часов):

W= n/Т_р.

 

Рис. 1. Классификация методов диагностирования энергетического оборудования в различных режимах.

 

Для диагностирования любой объект характеризуют совокупностью параметров Х_i, i=1,…k. Любой объект можно описать математической моделью М_j (j = 1, …m) с различной степенью адекватности. Параметры, которые характеризуют существенные свойства объектов, называют определяющими.

 

Система называется полностью наблюдаемой на интервале t _o, t ₁ , если ее начальное состояние Х (t _o) можно определить по измеренному на этом же интервале вектор Y (t ₁). Система называется полностью управляемой, если из некоторого начального состояния Х (t _o) ее можно перевести в любое другое состояние Y (t ₁) за конечный интервал времени t _o, t ₁ воздействием на нее управления G (t) € G, где G – некоторый заданный класс функций. Эти свойства определяют техническое состояние объектов и различимости и оценить возможности их реализации. Условия работоспособности представляют собой правила, разделяющие конечное множество Х возможных состояний объекта на подмножества работоспособных Х _o и неработоспособных Х_f состояний,

Х = Х _o U Х_f. Признак наличия дефекта представляет собой условия различия каждого состояния в подмножествах Х _o, Х_f. Если обозначить через ΔW_i приращение информации, получаемое при наблюдении за i- том (i = 1, n) параметром объекта или за совокупностью их, а через C _i - стоимость такого наблюдения, то отношение ΔW_i / C _i можно использовать в качестве оценки значимости и реализуемости наблюдения за i- тым параметром или за группой параметров.

 

Рис. 2. Структурная схема прогнозирующего диагностирования.

 

Эффективность всей отраслевой системы по отношению к энергоблоку можно оценить коэффициентом готовности:

К_г=W_о/(W_о+ W_n),

где W_о- мощность, отдаваемая энергоблоком за время 0, t;

W_n-мощность, которая теряется в результате простоя энергоблока.

 

Существуют и другие показатели, например: время обнаружения дефекта t_о, время поиска t_п, время устранения неисправности t_у, а также вероятность обнаружения и устранение за время 0,t и др.. На основе понятия элементарных операций рассматриваются элементарные действия, например, O_j (t), где j=1,k, осуществляемые совокупностью π{π _j1 … π _jr} определённых приборов, которыми работает персонал за определенный интервал времени. В качестве таких действий могут быть:

 

Например, показатель достоверности контроля является функцией от элементарных операций D = f {О₁ (t),О₂ (t), …, О_k (t)}. Аналогично можно записать выражения и для других показателей.

 

Схема принятия решения при диагностировании. Таблица 1

Действительное состояние объекта	Принятое решение	Условные обозначения события	Вероятность сложения событий
Исправен Q	Исправен S Неисправен Q‾	Q S Q Ѕ‾	Р1 Рл
Неисправен Q‾	Исправен S Неисправен Ѕ‾	Q‾ S Q‾ Ѕ‾	Рн Р2

 

Вероятность Р1 и Р2 соответствуют правильному принятию решения при контроле. Вероятность Рл - ложная вероятность, Рн - неправильная,

которая возникает вследствие неточности измерительных приборов.

Вероятность Рл- ошибка первого рода, а пропуск дефекта Рн- ошибка второго рода.

Pл = ò_Go W_c (x) dx x € G_o

Рн = ò _Gi W_c (x) dx x € G_1,

W_c (x)- функция распределения ошибки измерения;

G_o,G₁- области интегрирования, которые определяются выбранным допуском.

В этом случае достоверность диагностирования:

D=1-(Рл+Рн)

 

 

Схема принятия решения при диагностировании с учетом ненадёжности средств диагностирования. Таблица 2

Действительное состояние	Принятое решение	Условное обозначение события	Вероятность сложного события
ОД	CД
исправен Q	Исправны К	Исправен S неисправен S‾	QKS QK S‾	Р1 Рл
неисправен Q‾	Исправен S неисправен S‾	Q‾ KS Q‾ K S‾	Рн Р3
исправен Q	Неисправны К‾	Исправен S неисправен S‾	Q К‾ S Q К‾ S‾	Р1 Рл
Неисправен Q‾	исправен S неисправен S‾	Q‾ К‾ S Q‾ К‾ S‾	Рн Р2

 

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

где Р_д=Р (х _i)× Р(k_l /х _i)/R; Р(х _i) – вероятность возникновения отказа по причине х_i, Р(k_l /х_i) – условная вероятность признака k_l при наличии причины отказа х_i;

R =? Р(х _j)×Р(k_l /х _i), где m - число признаков. Составляющая Р_и определяется точностью используемых приборов для измерения и принятия решения (Р_и = D).

 

 

Пусть надежность оценивается коэффициентом готовности и средними потерями. Тогда при текущем диагностировании

_T^o _n - _T^o _об

К_г = (D_т /T^о _п) ò F_о (t) dt,

^o

где D_т- достоверность текущего диагностирования;

Т^о _п и Т^о _об -период и средняя длительность периодического обслуживания объектов;

F_о (t) -вероятность безотказной работы объекта.

В то же время вероятность готовности объекта:

Р_т(t_к)= Р_ср·D_т·F_оо(t_к),

Р_ср -вероятность нахождения средств в рабочем состоянии;

F_оо(t_к)-функциональное распределение длительности обслуживания объекта.

Рассмотрим средние потери W:

Ẅ= С1·Рл+С2·Рн

С1 и С2 -средние потери при ложном отказе и пропуске отказа.

Таким образом, при прогнозном диагностировании:

К_г=D_п(τ)- функция от τ, где D_п(τ) -показатель достоверности прогнозного диагностирования; τ-шаг прогнозирования, который принимает в описанных условиях значения: τ, 2τ, …, n τ.

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

 

Теоретической базой моделирования является теория подобия. Изучение условий подобия явлений производят на основе уравнения Навье-Стокса

- jр_х/jх + mj²w_х/jz² +rg _х=rw _хjw _х/jх + rjw _х/jt.

(P) (T) (G) (I₁) (I₂)

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

Правая часть уравнения учитывает действия инерционной силы (I₁,I₂). Скорость в каждой точке будет меняться в зависимости от времени

(rjw _х/jt –локальное изменение I₂), но в то же время скорость изменяется и от того, что рассматриваемая точка за это время изменила свое положение (jw _х/jх – конвективное изменение I₁).

 

Теоретической базой моделирования является теория подобия. Изучение условий подобия явлений производят на основе уравнения Навье-Стокса

- jр_х/jх + mj²w_х/jz² +rg _х=rw _хjw _х/jх + rjw _х/jt.

(P) (T) (G) (I₁) (I₂)

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

Правая часть уравнения учитывает действия инерционной силы (I₁,I₂). Скорость в каждой точке будет меняться в зависимости от времени

(rjw _х/jt –локальное изменение I₂), но в то же время скорость изменяется и от того, что рассматриваемая точка за это время изменила свое положение (jw _х/jх – конвективное изменение I₁).

Сила трения возникает вследствие градиента скорости, что связано с переносом из одного слоя в другой некоторого количества движения. Интенсивность этого процесса связана с вязкостью жидкости m.

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

 

Таблица 3

 

Отношения, имеющие обратные значения в таблице 3, отмечены звездочкой в практике во внимание не принимаются. Рассмотрим остальные десять комбинаций: P/T, P/G, T/G, P/I₁, T/I₁,G/I₁, I₂ /I₁, P/I₂, T/I₂,G/I₂.

Из ранее изученных дисциплин известны следующие критерии:

P/I₁ " Eu = r p/ r w² – критерий Эйлера; I₁ /T " Re = wlr /µ - критерий Рейнольдса; I₁ /G " Fr = w² /gl – критерий Фруда; I₁ /I₂ " Ho = w t /l – критерий гомохронности.

Все остальные критерии подобия можно получить из этих четырех:

P/T = P/ I₁ * I₁ / T = Eu Re; P/G = P/ I₁ * I₁ /G = Eu Fr; T/G = I₁ /G: I₁ /T = Fr/ Re; P/I₂ = P/I₁* I₁ /I₂ = Eu Ho; T/I₂ = I₁ /I₂: I₁ /T = Ho/ Re; G/I₂ = I₁ /I₂: I₁ /G = Ho/ Fr.

 

Рассмотрим аналитическое прогнозирование.

Пусть контролируемый параметр y(t) в области 0 - τ_n принимает значения y(τ₀), y(τ₁),…,y(τ_n), которые зафиксированы системой диагностирования. Необходимо по известным значениям y(τ_i) контролируемой функции в прошлом предсказать значения величин y (τ_n+1), …, y(τ_n+m), где τ_n+i > τ_n.

Для вероятностного прогнозирования по известным значениям y(τ_i) определяется вероятность того, что значения функции y(τ) не выйдут за допустимые пределы: Р {|y_i (τ_n+i) - y_n (τ) |} < e.

. В этом случае: y_i (τ_n+i) -значение контролируемой функции; y_n (τ) -требуемое изменение функции при исправном состоянии

Для длительной прочности предложена экспериментальная зависимость для твердых материалов:

t_p = t_o exp [(u_о - υ σ)/RT],

где t_o – постоянная; t_p – время до разрушения; u_о -эффективная энергия активизации в ненагруженном состоянии; υ- структурный параметр; σ- напряжение; Т- абсолютная температура; R- газовая постоянная.

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

t_p = AT² σ^-m1 exp [(u_о - υ σ)/RT], где m₁ – коэффициент.

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

Для прогнозирования используется так же зависимость Ларсена-Миллера и ее различные модификации.

Одна из модификаций имеет следующий вид:

t_p = t_o exp{- 2,3 (C – C₁/T)}, где С и С₁ – постоянные коэффициенты; Т-температура.

 

 

Для сложнонапряженного состояния материалов или конструкций используются выражения, которые получаются экспериментальным путем на основе уравнений полинома:

y(τ)=F₁(τ)·A₁ (τ)+F₂(τ)·A₂(τ)+…, где А- весовые коэффициенты;

F(τ) - составляющая функция.

 

 

Рассмотрим прогнозирование для электромеханического оборудования. Зависимость пробивного напряжения “u” от температуры Т имеет следующий вид: u(T)= u(T₁)exp {- a(Т₁-Т)}. Для проводов и кабелей имеется зависимость типа ℓn τ= A/Т+В, где Т- абсолютная температура; А, В –коэффициенты, зависящие от материала и конструкции кабеля.

 

Существует три вида вероятностного прогнозирования:

1. По всему известному прошлому о реализации процесса Z(t), который является результатом преобразования без помех интересующего процесса W(t)

2. По частично известному прошлому о реализации процесса Z(t), который не искажен помехами, а реализации представлена не всеми своими состояниями

3. По частично известному прошлому о наблюдаемом процессе Z (t), который является результатом преобразования W(t) при наличии помех

Математическая постановка задачи требует точного знания статических характеристик Z(t), W(t) и характеристик их совокупности { Z(t), W(t)}.

 

Моделью, которая характеризует состояние пароперегревателей как объекта диагностики, является зависимость, которая связывает временное сопротивление s^t_b с температурой T, давлением p, конструкцией трубы S, а также сортом металла g, т.е. s^t_b = f (T, P, S, g) (1).

Разрушение пароперегревателя происходит тогда, когда главное напряжение достигает некоторого допустимого (предельного) значения s_b ^доп, т.е. выполняется условие: u_b = s_b ^max - s_b ^доп ³ 0 (2)

Зависимость (1) можно найти экспериментально или рассчитать. Для экспериментального нахождения зависимости (1) используется метод планирования эксперимента.

Для расчета s^t_b необходимо знать текущие параметры технологического процесса, т.е. вычислить радиальное s^t_r и тангенциальное s^t_t напряжения в стенках труб с использованием функции напряжения, которая удовлетворяет неоднородному бигармоническому уравнению.

В такой постановке задача является достаточно сложной для практической реализации.

 

 

Рис. 3. Классификация вращающихся (роторных) механизмов.

 

Вероятность правильного диагноза по Байесу –это условная вероятность наличия состояния х_i после появления признака k_j:

P (x_i /k_j) = P (x_i) P (k_j /x_i) / P (k_j)

 

 

 

 

Рис.4. Тепловизор

 

 

Метод мажорирования состоит в том, что для измерения некоторого параметра “х” используют дополнительные датчики “y” и “z”, а результат измерения получают путем присвоения измеренному параметру значение по значениям большинства датчиков. Например, если измеряется цифровая величина, то решение принимается по правилу R = xy Ú xz Ú yz.

 

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

∆=Ө_ф - Ө_эт (1)

 

Рис. 5. Функциональная схема технического диагностирования.

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

2. Работоспособный объект - объект, который может выполнять возложенные на него функции.

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

4. Дефект - изменение состояния объекта в процессе его изготовления, эксплуатации или ремонта, которое потенциально может привести к уменьшению степени его работоспособности.

5. Неисправность - изменение состояния объекта, приводящее к уменьшению степени его работоспособности.

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

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

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

9. Диагностика(диагностирование) - процесс определения состояния объекта.

10. Диагноз - заключение о состоянии технического объекта.

 

Рис. 6. Структура технической диагностики.

 

 

Тренд на рисунке 7 (а) характеризуется 4-мя этапами:

1 этап Т1- приработка машины,

2 этап Т2- нормальная работы,

3 этап Т3- развитие дефекта,

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

 

Тогда остаточный ресурс в первом приближении может быть оценен как разность времени t _пр,соответствующему предельного значению диагностического симптома ∆_пр и времени t _изм, соответствующего значению диагностического симптома ∆_изм на момент измерения диагностического параметра.

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

Во многих практических случаях тренды могут быть немонотонными. На рис. 7,г представлен тренд, участок I которого характеризует развитие одного дефекта, на участке II наблюдается стабилизация уровня вибрации и на участке III производная изменения уровня вибрации увеличивается в результате появления еще одного дефекта. Тогда достоверный прогноз состояния объекта и оценка остаточного ресурса возможны только на последнем участке развития цепочки дефектов.

 

.

Рис.8. Схема основных операций функциональной и тестовой диагностики.

 

Рис. 9. Процесс формирования технического диагноза.

 

Алгоритмы диагностирования составляются при следующих допущениях:

1. Объект может находиться в конечном множестве состояний S, которое

делится на два подмножества S₁ (бездефектные состояния) и S₂ (состояния с различными дефектами, при которых объект остается работоспособным).

Каждое состояние из подмножества S₂ отличается степенью или запасом работоспособности.

2. Состояние объекта характеризуется совокупностью диагностических

показателей d₁, d₂,…..,d_к, которая представляет собой вектор состояния D:

D = (d₁, d₂,…..,d_к).

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

4. Условия работоспособности задают областью работоспособности исходя из следующих предположений:

- вектор состояний оборудования определен,

- существует номинальный вектор состояний,

- отклонения вектора состояний от номинального допускается только в определенных пределах,

- допустимые отклонения определяют область работоспособности.

Условия работоспособности задают по-разному для случаев использования в качестве диагностического показателя параметров или характеристик.

Если в качестве диагностического показателя используется один параметр,то условия работоспособности задаются неравенствами, которые ограничивают его значение с одной или с двух сторон. Таким образом, объект работоспособен, если все неравенства выполняется:

1. d_i >d_iн, d_i<d_iв или d_iн< d_i<d_iв ,

где d_i, d_iн и d_iв - соответственно, текущее, нижнее допустимое и верхнее допустимое значения диагностического параметра.

2. Каждый из диагностических показателей состояния d_j может определяться по совокупности диагностических параметров d_ji,…, d_jl: d_j = d_ji, …, d_jl.

3. Для каждого диагностического параметра d_i существует номинальное значение d_оi, область допустимых отклонений ∆_оi и предельное отклонение (порог опасного изменения параметра) ∆^пр_i, при превышении которого объект считается неработоспособным и должен быть остановлен.

Объект считается бездефектным, если для каждого параметра выполняется неравенство: |d _i -d _оi |≤∆ _оi, где ∆ _оi - порог допустимого отклонения.

Объект считается неработоспособным, если хотя бы для одного из параметров выполняется неравенство: | d _i -d _оi | >∆^пр_i, где ∆^пр_i - порог опасного изменения параметра.

Во всех случаях объект имеет ограниченную работоспособность.

В качестве диагностических показателей могут быть использоваться не только параметры, но и характеристики объекта, типа: у= f(x), где «х», «у» входная и выходная переменные соответственно.

В этом случае условие работоспособности объекта определяется степенью отклонения r_р (f, φ) текущих характеристик f(х) объекта от номинальной φ(х).

Сложные объекты в целом оцениваются как работоспособные при условии работоспособности каждого узла или структурной единицы.

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

 

Рис. 10. Простейшее гармоническое колебание.

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

- Х_пик - пиковое значение, которое характеризует максимальное значение колебаний, но не отображает его развития во времени;

_т

- Х_скз =√^`1/Т ò _о Х² (t) dt – среднеквадратическая (эффективное) значение, которое представляет собой квадратный корень из усредненных во времени, возведённых в квадрат мгновенных значений. Этот параметр имеет особую важность, т.к. связан с энергией, а следовательно, с разрушающей способностью колебаний.

Для гармонических колебаний:

- Х_скз = Х_пик / √^`2;

_т

- Х_сред = 1/Т ò _о /Х/ dt - среднее абсолютное значение, связанное с развитием механических колебаний во времени;

- Х_размах - размах колебаний (двойная амплитуда);

- К= Х_пик/Х_скз - пик-фактор (амплитудный коэффициент), числовое значение тем больше, чем больше выражен импульсный или случайный характер колебаний. Для синусоидальных колебания: К= √^`2.

В то же время для простейшего гармонического колебания можно записать:

Х(t) = Х_о cos (w_о t + j_о) = Х_о e^{- j (wot +jo)},

где Х_о.- амплитуда; w_о –частота; φ_о - начальная фаза; t-время. Следует отметить, амплитуда, частота и начальная фаза сигнала – три независимых параметра, каждый из которых может содержать диагностическую информацию. В то же время кроме гармонических составляющих в реальном сигнале шума или вибрации оборудования присутствуют, как правило, случайные составляющие, представленные рис. 11.

Рис. 11. Случайный сигнал.

 

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

Рис. 12. Сигнал вибрации электродвигателя.

Например, спектр простейшего гармонического колебания имеет вид, представленный на рис. 13,а (справа).

Рис.13. Временные сигналы и их спектры.

 

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

 

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

 

Рис. 14. Временной сигнал вибрации электрической машины и его спектр.

 

Акустическое поле в любой момент времени и в любой точке пространства полностью определяется зависимостью одной из трёх взаимосвязанных величин х, х^•, Р:

- х - вектор смещения частицы среды из положения покоя;

-х^• = dx/dt – вектор колебательной скорости;

- Р - звуковые давления в среде.

Каждая из этих величин в любой момент времени и в каждой точке пространства удовлетворяет волновому уравнению.

Для звукового давления волновое уравнение имеет следующий вид:

Ѳ Р = 1/С ²_о ¶²Р/¶t²,

где Ѳ = [¶²/¶x² + ¶²/¶y² + ¶²/ ¶z² ] - оператор Лапласа; С_о - скорость распространения волны (скорость звука в воздухе С_о » 330 м/с).

 

Плоская волна, представленная на рис.15 - это продольная волна, которая зависит только от одной декартовой координаты и времени.

 

Рис. 15. Плоская волна.

 

В этом случае: ¶²Р/¶y² = ¶²Р/¶z² = 0. Волновое уравнение продольной волны в этом случае в упругой среде без учета потерь для однородной бесконечной среды имеет вид:

¶²Р(х, t)/ ¶x² = 1/ С_о ¶²Р(х, t)/ /¶t².

Решение этого уравнения имеет вид:

P (х, t)= P_о exp j{w (t-x/с_о)},

где P_о - амплитуда переменной составляющей; wt – фаза временная; w (x/с_о) - фаза пространственная.

В безграничной и идеальной среде скорость распространения волны с_о определяется только параметрами среды: С_о = Ö`k/r,

где к- модуль всестороннего сжатия, Па; r_о - плотность воздуха, кг/м³.

Интенсивность звука связана с давлением: J = P²/r с_о. Здесь r с_{о _}- волновое сопротивление, кг/(м²·с (для воздуха r с_о » 440 кг/(м²·с).

 

Главная особенность слуха человека – логарифмическая зависимость восприятия звука от величины переменного давления в воздухе. Эта особенность позволяет человеку эффективно воспринимать звук в диапазоне изменения величины до 10⁶ … 10⁷ раз.

Вторая особенность - чувствительность органа слуха человека зависит от частоты звука, как представлено на графике рис. 16.

 

Рис. 16. Кривая равной громкости А.

 

 

составляющую на определенной частоте, например, движение камертона, колебание которого представлено на рис. 17,а.

Рис.17. Природа механических колебаний.

 

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

В реальных условиях вибрация является сложным механическим колебанием со многими составляющими на разных частотах, которое можно охарактеризовать вибросмещением «х», виброскоростью х^• = dx/dt и виброускорением х^•• = d²x/dt².

Для гармонических колебаний: x(t)= х_о sin (wt+φ_о). В этом случае виброскорость: x^• (t)= х_о ·w·сos (wt+φ_о)= x_о·w·sin (wt+φ_о+π/2). Виброускорение:

x^•• (t)=- х_о·w²·sin (wt+φ_о)= x_о·w²·sin (wt+φ_о+π).

Обозначим: d= x_о -амплитуда перемещения, n= x_о·w- амплитуда виброскорости

a= x_о·w²- амплитуда виброускорения.

Как следует из графиков рис. 18 виброскорость имеет опережение фазы на 90^о относительно вибросмещения, а виброускорение – на 180 ^о .

Рис. 18. Опережение фазы виброскорости и виброускорения относительно вибросмещения.

 

Виброскорость гармонических колебаний может быть определена путем деления их виброускорения на 2p¦, а вибросмещение – делением виброускорения на 4p²¦². В этом случае: n= a/2p¦ = ò a dt; d= a/4p²¦² = òn dt.

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

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

Рис. 19. Спектры виброперемещения, виброскорости и виброускорения асинхронного электродвигателя.

 

Для шума линейными единицами измерения являются:

- паскаль (Па) – при измерении звукового давления;

- ватт на метр квадратный (Вт/м²) – при измерении интенсивности звука.

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

Уровень звукового давления, выраженный в дБ определяется следующей формулой L_p = 20 lg (P/P_nop), где Р –давление, Па; P_nop = 2·10 ^{– 5} Па – порог слышимости человека на частоте 1 кГц.

Интенсивность звука J (Bт/ м²) связана со звуковым давлением Р:

J=Р²/(ρ·c_о), (1)

где ρ· c_о -волновое сопротивление; ρ – плотность среды, равная для воздуха»1,29 кг/ м³; c_о » 330 м/с – скорость звука в воздухе.

Пороговые значения давления Р и интенсивности J выбраны такими, чтобы для воздуха в нормальных условиях они удовлетворяли соотношению (1), поэтому интенсивность звука и звуковое давление в воздухе, выраженные в дБ, совпадают.

Уровень интенсивности звука L_j (дБ) связан с линейными единицами:

L_j = 10 lg (J / J_nop), (2)

В этом случае: J - интенсивность звука, Вт/ м²; J_nop = 10 ^{– 12} Вт/ м² - пороговое значение интенсивности, полученное из выражения (1).

Обычно волновое сопротивление для воздуха: ρ· c_о»440 кг/(м²·с).

 

Для вибрации линейными единицами измерения согласно ГОСТам ИСО являются:

- мкм (микрометр) –для измерения вибросмещения,

- мм/с – для измерения виброскорости,

- м/с² – для измерения виброускорения

Вибрация может выражаться в относительных единицах – децибелах..

Тогда уровни колебаний соответственно виброперемещения, виброскорости и виброускорения (дБ), определяются следующим образом

L_х = 20 lg (Х / Х_nop); (3)

L_х^• = 20 lg (Х^• / Х^•_nop); (4)

L_х•• = 20 lg (Х^•• / Х^••_nop), (5)

где Х_nop, Х^•_nop, Х^••_nop - пороговые значения.

В промышленности используется ряд значений:

Х_nop» 8·10 ^{– 12} ,м; Х^•_nop» 5·10 ^{– 8},м/с; Х^••_nop» 3·10 ^{– 4}, м/с².

Для виброскорости пороговые значения определяются из соотношения, в котором связано давление в упругой волне со скоростью колебания частиц Х^• при замене давления Р на пороговое значение давления P_nop:

P_nop = ρ· c_о Х^•_nop (6).

Тогда для воздуха пороговое значение виброскорости равно

Х^•_nop» 5·10 ^{– 8},м/с. (7).

Следует отметить, что это пороговое значение не зависит от частоты.

Пороговые значения виброперемещения и виброускорения выбираются таким образом, чтобы Х_nop, Х^•_nop, Х^••_nop, измеряемые в одной точке и выраженные в дБ совпадали на частоте ¦₀=1000 Гц.

В настоящее время в соответствии со Стандартом ISO -1683 рекомендуются следующие пороговые значения механических колебаний:

Х_nop» 10 ^{– 12},м; Х^•_nop» 10 ^{– 9},м/с; Х^•_nop» 10 ^{– 6},м/ с², приведенные не к частоте ¦₀=1000 Гц, а, для получения более простых значений к круговой частоте

w₀ = 1000 рад/с (¦₀ » 159 Гц).

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

 

На рис.20 представлен спектр колебаний (обозначен пунктиром) и общее значение амплитуды или общий уровень (обозначен сплошной линией), измеренный в широкой полосе частот.

 

Рис.20. Общее значение амплитуды или общей уровень колебаний в широкой частотной полосе.

 

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

Рис. 21. Структурная схема прибора для измерения сигналов вибрации и шума.

 

высокочастотных сигналах вибрации исправного подшипника качения и подшипника с дефектом (раковиной на поверхности качения), который поясняется графиками рис.22.

 

Рис. 22. Высокочастотная вибрация исправного (а) и дефектного (б) подшипников качения.

 

 

Метод ударных импульсов реализуется достаточно простым прибором, структурная схема которого представлена на рис.23.

Рис. 23. Структурная схема прибора для измерения пикфактора сигнала вибрации.

 

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

 

На рис.24 представлены спектры вибрации.

 

 

 

Рис. 24. Спектры вибрации, измеренной с различной разрешающей способностью.

 

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

 

 

Используется два типа анализаторов спектра:

1. Простейший прибор (рис.25, а), который снабжен перестраиваемым по частоте полосовым фильтром.

2. Прибор с последовательно переключаемыми фильтром (рис.25,б).

 

Рис. 25. Упрощенные структурные схема анализаторов, осуществляющих последовательный спектральный анализ.

 

 

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

 

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

 

Рис. 26. Структурная схема анализатора, работающего в реальном времени и осуществляющего параллельный спектральный анализ.

 

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

 

Амплитудно-частотная характеристика (АЧХ) идеального полосового фильтра представлена на рис. 27.

Рис. 27. Амплитудно-частотная характеристика идеального полосового фильтра.

 

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

 

Рис.28. Амплитудно-частотные характеристики идеальных полосовых

фильтров с постоянными абсолютными значения ширины полосы пропускания частот.

 

Реальные ПФ имеют АЧХ, которая отлична от характеристики идеального фильтра. АЧХ реального фильтра имеет менее крутые фронты, представленные на рис. 29.

 

Рис. 29. Амплитудно-частотная характеристика реального (а) и идеального (б) фильтров.

 

 

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

 

По второму способу ширина полосы пропускания реального фильтра определяется как интервал частот на уровне спада частотной характеристики фильтра на 3 дБ (рис. 29, б).

Для фильтров с постоянной абсолютной шириной полосы пропускания средняя частоты определяется как сумма верхней и нижней частоты фильтра, деленная на два: f_о = (f_в + f_н)/2, а ширина полосы пропускания равна

Df = f_в - f_н, где f_в и f_н - верхняя и нижняя частоты фильтра, соответственно.

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

Рис. 30. АЧХ фильтров анализатора при постоянной ширине их полос пропускания.

 

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

пропускания имеют среднегеометрическая частоту, равную f_о = Öf_вf_н, и ширину полосы пропускания, определяемую как D f_отн = {[(f_в - f_н) ×100]/ f_о}%.

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

Рис.31. Амплитудно-частотные характеристики идеальных фильтров с постоянными относительными значениями ширины полосы пропускания.

 

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

представленных на рис 32.

Рис. 32. Связь между добротностью Q и относительной шириной частотной полосы rƒ_отн.

 

Это позволяет эффективно применять анализаторы с фильтрами с постоянной относительной шириной полосы пропускания при использовании

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

На практике широко используются два типа фильтра с постоянной относительной шириной полосы пропускания:

- октавный фильтр;

- третьоктавный фильтр.

Октава – это диапазон частот, в котором верхняя частота f_в отличается от нижней частоты f_н в два раза. В этом случае для октавного фильтра справедливо соотношение f_в = 2f_н, при этом для среднегеометрического значения частоты имеем f_о = Öf_вf_н = Ö2 ´f_н.

Абсолютная ширина полосы пропускания октавного фильтра равна

Df = f_в - f_н, а для относительной ширины полосы пропускания имеем

D f_отн = {[(f_в - f_н) ×100]/ f_о}%.= 70,7 %.

Если октавный фильтр имеет среднегеометрическую частоту f_о =10 Гц, то абсолютное значение его полосы пропускания равно ∆f=7 Гц.

Если среднегеометрическая частота октавного фильтра f_о =100 Гц, тогда

∆f=70 Гц.

По рекомендациям Международной электротехнической комиссии (МЭК) установлен ряд. Опорным значением частоты этого ряда является f_о =1000 Гц.

Если октавную полосу разделить на три геометрически равных участка, получим третьоктавную полосу. Для октавного и для третьоктавного фильтров абсолютные значения ширины полосы пропускания увеличивается по мере роста их среднегеометрических частот, оставаясь в процентном отношении постоянными – 23% для третьоктавного фильтра и 70% для октавного фильтров. На рис. 33 приведена «гребенка» фильтров третьоктавного анализатора с логарифмической шкалой частот.

 

Рис. 33. АЧХ смежных третьоктавных фильтров с логарифмической шкалой частот.

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

 

Для определения оптимальной ширины полос пропускания частот можно использовать эмпирическое правило: возможный диапазон изменения частоты вращения оборудования может находиться в пределах ±10%. В этом случае изменение частоты в одну сторону составляет 10% и изменение частоты в другую сторону – тоже 10%. Итого 20%, а ширина полосы пропускания третьоктавного фильтра – 23%, что во многих случаях удовлетворяет требованиям к контролю уровней шума и вибрации оборудования.

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

 

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

 

Рис. 34. Случайный амплитудно-модулированный сигнал.

 

В этом случае глубину модуляции m сигнала можно определить в процентах:

m = {[X_max (t) - X_min (t)]/ [X_max (t) + X_min (t)]}100%,

где X_max, X_min - максимальное и минимальное значения огибающей сигнала соответственно.

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

 

 

Рис. 35. Временные сигналы высокочастотной вибрации подшипников качения и спектры ее огибающей: а- исправный подшипник; б – подшипник с износом поверхности трения; в – подшипник с раковиной на поверхности качения.

 

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

 

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

 

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

К достоинствам метода спектрального анализа огибающей следует отнести:

- возможность локализации дефекта;

- высокая чувствительность;

- высокая достоверность определения вида и величины каждого дефекта;

- возможность постановки диагноза и прогноза по однократным измерениям вибрации.

 

В таблица 4 представлены четыре основных вида нагрузки в подшипниках качения.

Таблица 4.

 

 

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

Рис. 37. Схема микрофона.

 

Тонкая никелевая мембрана, закрепленная параллельно проводящему основанию, образуют совместно с ним конденсатор емкостью C_(x). Последовательно к конденсатору C_(x) включены конденсатор C₍₀₎ и источник постоянного напряжения e. В этом случае переменное U_вых(t), пропорционально величине переменного давления воздуха P(t), снимается с емкости C_о.Определим связь перемещения мембраны х(t) с действующим на нее звуковым давлением P(t) следующим образом.

По закону Гука, сила, действующая на мембрану, равна F= kx(t), где k - жесткость мембраны; x(t) – смещение мембраны под действием звукового давления.

С другой стороны, для этой силы справедливо следующее выражение F= P(t)•S, где S- площадь мембраны. Тогда kx(t) = P(t) •S, отсюда x(t) = P(t) •S/ k. Следовательно, перемещение мембраны x(t) прямо пропорционально величине звукового давления P(t).

Для определения связи выходного напряжения U_вых(t) от перемещения x(t) мембраны рассмотрим эквивалентную схему микрофона, представленную на рис. 38.

Рис. 38. Эквивалентная схема конденсаторного микрофона.

 

Для выходного напряжения в соответствии с эквивалентной схемой микрофона имеем U_вых(t)= e C_(x) / (C_(x) + C_o).

При отсутствие перемещения мембраны емкость конденсатора определяется следующим выражением C_(x=о) = e_оeS/d_о, где e_о = 8,85×10^-12 Ф/м - диэлектрическая проницаемость воздуха; e - относительная диэлектрическая проницаемость; d_о – зазор между проводящим основанием и мембраной в исходном состоянии.

При действии звукового давления емкость конденсатора C_(x) определяется выражением C_(x) = e_оeS/[d_о - x(t)], где x(t) – смещение мембраны под действием звукового давления.

Максимальная чувствительность микрофона достигается при условии

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

U_-вых ~ P(t), т.е. переменная составляющая выходного напряжения микрофона прямо пропорцио