Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Вибродиагностика машин и механизмов. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Рис. 3.17

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

стандартном отклонении.

Другая задача анализа состоит в определении границ зоны допустимых значений измеряемых параметров. Если количество измерений достаточно велико ( N і 30 ), то пороги определяются довольно просто. Для этого необходимо лишь задать вероятность p того, что (t) ока-

жется выше величины порога

 

п .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ґ

-

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

 

2

 

 

p

 

(t) > n

 

=

 

 

 

 

 

т e

2

 

d .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

Для p { (t) >

n } = 0.05 ,

 

 

 

 

 

 

 

 

 

 

 

n =

+ 2 .

 

 

 

 

Для p { (t) >

n } = 0.005 ,

 

 

 

 

 

 

 

 

 

 

n =

+ 3 .

 

 

 

 

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

случайной величиной и подчиняется распределению 2 .

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

Первое условие: оптимальный выбор интервалов и минимального количества измерений.

Второе условие связано с возможным появлением случайных (ошибочных) измерений.

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

сле ремонта или по мере износа.

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

41

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

3.3. АНАЛИЗ ВИБРАЦИЙ В ЧАСТОТНОЙ ОБЛАСТИ

3.3.1. ФИЛЬТРАЦИЯ СИГНАЛА

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

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

Рассмотрим различные формы фильтров (рис. 3.18). Пунктиром на этом рисунке показаны идеальные фильтры.

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

K

 

K

 

 

Фильтр

 

 

 

Широкополосный

 

нижних

 

 

 

 

 

частот

 

 

 

f

 

f

 

 

 

1000

 

100

300

42

K K

Фильтр пробка

Узкополосный фильтр

 

 

f

 

 

f

50

100

150

50

100

150

Рис. 3.18

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

Следует отметить, что АЧХ фильтра является только вещественной частью коэффициента передачи фильтра, который имеет вид

K( j ) = K( )e j ( ) .

Здесь K ( ) – амплитудно-частотная характеристика;

– фазочастот-

ная характеристика (ФЧХ);

 

– круговая частота.

 

 

Основное требование к ФЧХ – это ее линейная зависимость от

частоты, для сохранения формы сигнала.

 

 

 

 

 

 

В большинстве случаев АЧХ и ФЧХ связаны друг с другом. Пред-

ставим K ( j ) в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K( j ) = eA( )+ j ( ) = eG( j ) ,

 

где

A( ) = ln K ( )

логарифмическое затухание;

G( j ) = A( ) +

+ j

( ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Связь между A(

)

и

(

) определяется преобразованием Гильберта:

 

 

 

 

 

 

2

 

 

Ґ

 

A(

)

 

 

 

 

 

 

 

(

i ) =

 

i

т

 

 

 

d

,

 

 

 

 

 

 

 

 

 

 

2

-

2

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

i

 

 

 

 

 

A(

 

i ) =

2

Ґ

 

 

 

 

 

(

)

d .

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

-

2

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

43

 

 

 

 

 

 

 

 

Многие параметры фильтров обычно нормализуются. Особое место занимает АЧХ, причем нормируется логарифмическая АЧХ.

А

3 кл.

2 кл. 1 кл.

f

Рис. 3.19

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

Из полосовых фильтров применяются широкополосные с постоянной относительной полосой пропускания (третьоктавные, октавные и др.). В этих фильтрах нормируются: неравномерность АЧХ в пределах полосы фильтра, затухание. Иногда нормируется центральная частота и т.д.

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

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

Все эти характеристики используются при эксплуатации фильтров. При разработке фильтров используются еще ряд дополнительных характеристик.

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

44

дискретизации сигнала, которая должна быть в два раза выше fв

верхней частоты сигнала.

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

электронный фильтр для подавления частот выше fв .

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

3.3.2. СПЕКТРАЛЬНЫЙ АНАЛИЗ

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

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

 

a0

Ґ

 

 

a0

Ґ

 

 

x(t) =

+ е An cos(n

1t -

n ) =

+ е (an cos n

1t + bn sin n

1t) ,

 

 

2

n=1

 

2

n=1

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

= arctg

bn

 

где

A =

a2

+ b2

,

n

,

 

 

n

n

n

 

 

an

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

b =

2

2

x(t)sin n tdt,

 

 

 

 

n

 

1

 

T

 

-тT

 

 

 

2

 

a =

2

2

x(t)cos n tdt,

 

 

 

 

 

 

 

 

n

 

 

T

-тT

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

a

 

1

2

 

 

 

 

0

 

=

 

 

т x(t)dt ,

 

 

2

 

T

 

 

 

 

-

T

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

а в комплексном виде как

45

 

 

 

 

 

 

T

 

 

x(t) =

1

Ґ A*e jn

1t dt, A =

2

 

2

 

x(t)e- jn

1t dt .

 

 

 

 

n

 

n

 

 

тT

 

 

 

2 -тҐ

 

 

T

 

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

 

 

 

2

 

 

Для периодического сигнала ряд Фурье заменяется интегралом

Фурье:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

x(t) =

т S* ( )ei t d ,

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ґ

S* ( ) = т x(t)e- j t dt.

Второе выражение называется прямым интегральным преобра-

зованием Фурье, первое – обратным преобразованием.

Сравнение спектров периодического сигнала, имеющего период T, с непрерывным сигналом в той же полосе частот и одинаковой длительности позволяет сопоставить An ( ) и S ( ) следующим образом:

 

 

 

A =

1

S(

 

) ,

 

 

 

n

 

 

 

n

 

 

где

= n - n-1 =

2

.

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

Современные средства спектрального анализа сигналов используют преобразование Фурье в дискретной форме.

Для этого из сигнала x(t) на интервале 0–Т берется N равноотстоящих отрезков. Заменяя Т на количество отрезков N и обозначая x(ti ) как xi , с учетом, что n = 2 n / N , получим прямое и обратное преобразование Фурье следующего вида:

 

1

N -1

- j

2

ni

 

*

 

 

е x(t)e

 

N

 

Sn

=

 

 

, xi

 

 

 

 

 

 

N i=0

 

 

 

 

( N -1)/2

2

 

= е

Sn*e j

N

ni .

-( N -1)/2

 

 

46

Полученная система уравнений решается в матричном виде, и для

ее решения необходимо провести N 2 операций.

Интервал, на котором производится дискретизация сигнала x(t) ,

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

(рис. 3.20).

А

 

8

А

 

8,56

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

n

Рис. 3.20

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

Обычно в задачах мониторинга рекомендуется усреднять не менее четырех спектров. При решении динамических задач усреднение проводится по 10…20 спектрам.

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

ние

 

» fвр / 10 , где fвр – циклическая частота вращения ротора, на-

2

 

 

пример, ~ 5 Гц.

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

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

47

G12* ( ) = S1* ( )S2* ( ) .

Вещественная часть этого взаимного спектра называется синфазным (совпадающим) спектром, мнимая – квадратичным спектром.

Взаимный спектр связан с функцией взаимной корреляции преобразованием Фурье (теорема Винера–Хинчина), где

G12* ( ) = т R12 (

)ei

d ,

R12 ( ) = тG12* ( )ei

d .

Количественная оценка связи двух сигналов определяется нормированным взаимным спектром, называемым функцией когерентно-

сти,

 

 

 

 

G*

(

)

 

 

 

 

122

( ) =

 

 

12

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

S* ( )

 

S

2

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

минимального количества усреднений для нахождения 122 .

Опыты показывают, что необходимо усреднить не менее 30 спектров.

Кроме спектров взаимной мощности в некоторых задачах, связанных с определением форм колебаний, можно использовать фазовые взаимные спектры.

3.3.3. СПЕКТРАЛЬНЫЙ АНАЛИЗ ОГИБАЮЩЕЙ

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

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

48

Формирование огибающей сигнала по времени является наиболее эффективным способом выделения моделирующего компонента. В этом случае частоты огибающей и модулируемого сигнала различны. При этом рассматривают два случая:

1)оба компонента периодические.

2)один из компонентов – стационарный случайный процесс. Кроме спектрального анализа огибающей сигнала по времени, су-

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

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

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

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

x(t) = (t) (t) ,

где (t), (t) – моделирующая и стационарная случайная составляю-

щая сигнала; (t) = 1 + m cos

t , – частота модуляции.

Тогда x(t) = (1 + mcos t)

(t) .

3.3.4. КЕПСТРАЛЬНЫЙ АНАЛИЗ

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

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

49

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

3.4. МНОГОПАРАМЕТРИЧЕСКИЙ АНАЛИЗ

3.4.1. МНОГОМЕРНЫЙ АНАЛИЗ

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

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

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

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

50

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]