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

книги из ГПНТБ / Теория и практика балансировочной техники

..pdf
Скачиваний:
64
Добавлен:
27.10.2023
Размер:
31.98 Mб
Скачать

л

 

при

£ K <G„(co„);

 

 

З

 

(ь>в)

 

 

 

 

Т2 :

 

при £ „ > G „ K ) .

 

 

 

 

 

 

 

 

На рис. 6 изображена

блок-схема балансировочной машины,

построенная на основании

 

результатов,

полученных

для

опти­

 

 

мальных измерителей

пара­

ХУД |

 

метров

сигнала от дисбаланса

|

в присутствии помех.

 

 

 

 

Особенностями схемы

явля­

J S 3 J

 

ются:

 

бесконтактный

«генера­

г

 

тор»

(преобразователь)

опор­

 

ного

сигнала, выполненный на

 

 

 

 

базе

импульсно-фазовой

авто-

 

 

 

 

5

 

 

/4

 

 

 

 

 

 

15 13

Ад

12

I

10

 

J

Рис. 6. Блок-схема балансировочной машины:

1, 2 — датчики дисбаланса; 3 — ротор; 4 — фотоэлектри­ ческий датчик опорного импульса; 5 — импульсно-фазовая автоподстройка частоты; 6 — усилитель-ограничитель; 7 — импульсно-фазовыб детектор; 8 — запоминающее устрой­ ство; 9 — фильтр нижних частот; 10 — управляющий эле­ мент; // — подстраиваемый генератор; 12 — фазовраща­ тель; 13 — фазовый детектор; 14 — блок усилителей, на­ стройки и эталонирования; 15 — измеритель дисбаланса

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

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

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

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

создаваемых

в работающей турбомашине, может,

например,

оцениваться

дисперсией 0 ^ = 0,025(2,5 4-3°). Это

позволяет

рекомендовать приведенную схему при проектировании и совер­

шенствовании

измерительных

устройств

 

балансировочных

машин.

 

 

 

 

 

ЛИТЕРАТУРА

 

 

 

 

 

1. Воронин Г.

И., Федосеев Н. М.,

Блохин Ю. И.

Высокоскоростные

совмещенные опоры. «Подшипниковая промышленность»,

1964, №

2.

2. Блохин Ю. И., Панфилов Е. А.,

Федосеев Н. М.

Высокоскоростные

совмещенные опоры. Под ред. Г. И. Воронина. М.,

НИИНавтосельхозмаш,

1966.

 

 

 

 

 

3. Шахгильдян

В. В., Ляховкин А. А. Фазовая автоподстройка частоты.

М., изд-во «Связь»,

1966.

 

 

 

 

4. Пугачев В.

С. Теория случайных

функций. М.,

Физматгиз,

1962.

5.Пестряков В. Б. Фазовые радиотехнические системы. М., изд-во «Со­ ветское радио», 1968.

6.Самсаев Ю. А., Трунаев Ю. В., Феоктистов П. Н. Автоматическое

устройство для получения

опорных импульсов. «Механизация и автоматиза­

ция производства», 1971, №

6.

А. С. ГОЛЬ ДИН

ИСПОЛЬЗОВАНИЕ ЭЦВМ ПРИ УРАВНОВЕШИВАНИИ ТУРБОАГРЕГАТОВ

Увеличение мощности в энергомашиностроении влечет за со­

бой рост материальных потерь при простоях и пусках

агрегатов

на холостом

ходу,

например, при балансировочных

работах.

В этих условиях

целью

уравновешивания является

не просто

достижение

минимальных

амплитуд вибраций в контролируе­

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

ЭЦВМ для балансировочных расчетов.

Последнее

приводит

к необходимости создания специального

алгоритма

обработки

4*

51

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

При внедрении ЭЦВМ в практику уравновешивания необхо­ димо решить две группы задач:

а) разработка алгоритма и программы балансировочных рас­ четов, а также методики уравновешивания;

б) создание оперативной, доступной и надежной системы обеспечения балансировочных работ вычислительными сред­ ствами. Рассмотрим вариант решения этих задач, реализованный в производственных условиях.

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

часть исходных данных для расчета.

Балансировочные

грузы

устанавливают в доступных плоскостях

валопровода,

располо­

женных в разных местах вдоль его оси, и варьируют

по

массе

и углу установки.

 

 

 

Выбор обоснованного варианта установки балансировочных грузов для каждого балансировочного цикла составляет содер­ жание алгоритма балансировочных расчетов. Мы не будем рас­

сматривать

приемы измерений вибрации при

уравновешивании

в условиях

электростанций. Под вибрацией

ниже понимается

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

вует,

то подразумевается

или

модуль комплексной

величины,

или действительная величина.

 

 

 

 

точках і

Если

на валопровод,

имеющий

в контролируемых

вибрации

А І

(і = 1, 2, 3,...,/),

в

различных

балансировочных

плоскостях k

установить

балансировочные

грузы Pk

(k = 1, 2,

З,

К, где / ( < / ) , то для новых значений

вибрации

 

 

 

 

б ; = 1 ; + У я д . А ,

 

 

(1)

 

 

 

 

 

k

 

 

 

 

здесь aih •— балансировочные чувствительности

точек I

к плоско­

стям

k.

 

 

 

 

 

 

 

 

Необходимо определить такие значения

Pk,

при которых до­

стигаются минимальные значения є*. Идеальным было бы, когда ЄЇ = 0, однако при J < К это недостижимо. Таким образом, не­

обходимо сформулировать функцию качества

Ф, минимизация

которой определит оптимальные значения Ph.

Весьма

распро­

страненным является применение метода наименьших

квадра­

тов [2], при этом

 

 

ф = 2 8 ? -

 

(2)

Существуют и другие формулировки функции качества:

Ф = > ? е ? ,

г д е п > 1 ,

(3)

і

 

 

 

И

 

і ф т.

 

Ф = е т , где е.т

> є, при

(4)

Очевидно, что при п - > оо выражение

(3) равносильно

соот­

ношению (4), таким образом наиболее общей является форму­ лировка функции (3).

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

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

(4).

Однако она имеет

слабости. Если происходит снижение гт

за

счет увеличения

2є?_

то это всегда связано с повышением є* для

і Ф т, при этом

но­

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

Sef= min, то это означает, что использованных балансировоч-

І

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

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

І

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

ф = М 2

в? = 2 (в? + Яв,) = min,

(5)

 

 

і

 

 

 

где Dsi — дисперсия є,, являющаяся

функцией

дисперсий

исход­

ных данных.

 

 

 

 

 

Минимизация

суммы

квадратов

амплитуд,

сложенная с их

дисперсиями, при

расчетах уравновешивающих

грузов

исполь­

зуется достаточно часто. При этом предполагается, что значения a%h известны, а их дисперсии назначаются произвольно.

Допустим, что: а) существует нормальный закон распреде­ ления погрешностей; б) имеет место круговое распределение

погрешностей; в) центром распределения является аффикс вво­ димой в расчет комплексной величины; г) погрешностью уста­ новки грузов Ри можно пренебречь.

Принимая во внимание эти допущения, имеем

Dlt=DAi + ^{p\D^^

= DAt + ^ \ { P \ x + Ply)Daik\.

(6)

 

k

 

k

 

Условие минимума Ф:

 

 

 

дФ

дФ

= 0

(fe = 1, 2, 3,. . ., Л').

(7)

дРкх

дРку

 

 

 

Представив комплексные величины в алгебраической форме, произведя соответствующие подстановки и преобразования, из выражений (1), (5), (6) и (7) получим систему линейных урав­ нений с 2/С неизвестными, т. е.

{[aY[a}

+ [D]}{P]

+ [aY[A] = 0,

(8)

где [а] — матрица балансировочных

чувствительностей;

 

[а]т — транспонированная матрица [а];

 

[D]— диагональная матрица дисперсий;

 

[А] — вектор исходных значений

вибрации;

 

[Р] — вектор уравновешивающих

грузов;

 

fllli

вПу

-й\Кх—О-ХКУ

 

Ч\\уй\\х-

-й\Куа\Кх

 

[а] =

 

 

 

 

 

Я/1* О-Лу

• -dJKx — CLjKy

 

.O-llyO-JU-

CtJKxflJKy

 

[D] =

 

Do

 

 

 

 

 

 

 

 

 

 

 

DK D,

 

 

 

= У

Dalk;

 

 

A l x

 

 

~P\x~

 

 

 

 

 

Ply

 

[A] =

 

;

[P]

=

 

 

A.ix

 

 

PKx

 

-Pky -

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

 

[P\ = {[ay[a] + [D])-l[ay[A].

(9)

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

Для получения необходимых сведений при использовании К грузов (систем грузов) необходимо провести К + 1 балансиро­ вочных пусков. Если считать, что каждый последующий пуск

производится после

установки

очередного

пробного груза Pknp,

то

 

 

 

 

 

 

A i k ~ A i ( k - \ ) .

 

u i k

~

 

,

(10)

 

 

Pk

np

 

 

 

 

 

 

Па

Dj<*

+

Dj«*-»)

 

uaik

 

~2

 

 

 

 

Pknp

 

 

где k — номер пуска

(пуск без установки грузов имеет номер 0).

Порядок сравнения пусков, который предусматривается вы­

ражением (10), является

не единственным

и в общем случае не

оптимальным.

 

 

 

 

 

При наличии погрешностей имеет также значение, из замеров какого пуска образуется матрица [А]. Таким образом, если мат­ рицу [А] образовать из замеров пуска 0, а матрицы [а] и [D], пользуясь выражением (10), то решение, полученное из соотно­ шения (9), по критерию (5) не будет оптимальным.

Наиболее естественным путем получения оптимального ре­ шения является рассмотрение всех возможных вариантов обра­ зования этих матриц и выбор того, который обеспечивает mini­ mum minimorum для функции Ф. Поскольку такой путь преду­ сматривает весьма громоздкие вычисления, принята более эко­ номная схема расчета, которая приводит к практически равно­ ценным результатам и является более приемлемой в условиях оперативного использования ЭЦВМ. Матрица [А] образуется из замеров вибрации пуска, для которого ~ZAi — наименьшая. Мат-

i

рицы [а] и [D] образуются путем сравнения пуска, имеющего

наименьшее значение 2 Daih, со всеми остальными.

і

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

1. Определение номера пуска Ь, для которого сумма квадра­ тов вибраций минимальна: £Л,т> = min.

2. Определение номера пуска d, для которого сумма диспер­

сий замеров минимальна: 2

DAici

=

min.

 

 

 

 

 

 

і

 

 

 

 

3.

Образование матрицы [А] из замеров пуска Ь.

4.

Образование

матрицы

[а] путем

замены aih соответствую­

щими разностями

вибраций:

 

 

 

 

для

(k — l)

<d

aik

= Ai(h-l)

A~id\

 

 

для

(k • 1)

= d перейти к следующему

k;

для

(k—l)

> d

aih

= Aik

Aid.

 

 

 

5.

Образование

матрицы [D] путем такой постановки:

для

 

(k — 1)

< d

2

Daik

=

2 [ Щ < № _ „

+ DAid];

 

 

(k 1)

 

і

 

 

і

 

 

 

для

 

— d перейти к следующему

k;

для

 

(k—l)

>d

2 Daik

=

2 (DAik

+

DAid).

іі

6.Нахождение матрицы [P] из выражения (9). В отличие от

уравновешивающих

грузов

члены этой

матрицы обозначены Ps.

 

7. Определение

уравновешивающих

грузов

Ph, которые необ­

ходимо

установить на

ротор, предварительно

сняв все пробные

грузы

Pknp:

 

 

 

 

 

 

при

k-^b

и

k^Cd

Pfe

=

 

 

 

при

k>

b я

k^Cd

Pk

=

 

 

 

npnk^bnk>d

 

Pk

= |

1 +yps

\Pknp,

 

при

k>

b и k> d

Pk

=

 

 

 

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

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

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

балансировки,

складывается из, двух-трех замеров

вибрации

в одной точке

при каждом пуске. Такое количество

замеров,

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

Предложена определенная последовательность образования исходной информации: при каждом пуске производятся два за­

мера

вибрации

во всех

контролируемых точках

А' и А"

через

определенный

промежуток

времени. В

качестве

исходного

при­

нимается

значение А, полученное при

втором замере, т. е. А =

= А",

а

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

разности

двух

замеров

оА

= \А'Г — А'\.

Как известно, DA — а\_

Разность значений А' и А" обычно выходит за пределы воз­

можной

погрешности измерительной аппаратуры

и объясняется

в основном погрешностью объекта измерения, т. е. переменным вибрационным состоянием балансируемого агрегата. Можно считать, что помимо дисбаланса, вибрация зависит от неконтро­ лируемых тепловых деформаций, а также от точности установ­ ки балансировочного числа оборотов. Тепловое состояние агре­ гата меняется во времени, асимптотически приближаясь к уста­ новившемуся. По разным соображениям достичь установивше­ гося теплового состояния агрегата при балансировке не удается. Таким образом, разные балансировочные пуски, строго говоря, несопоставимы. Эта несопоставимость непосредственно обнару­ жится, если агрегат пустить дважды в разное время при одной установке балансировочных грузов. Практически такой способ обнаружения погрешности является неприемлемым. Изменение состояния агрегата между замерами А' и А" моделируют воз­ можное несоответствие состояний агрегата пр_и разных пусках.

Интервал времени между замерами А' и А" определен из рассмотрения изменений вибрации во времени для различных агрегатов и выбран так, что в течение него происходит 20—40% всех изменений вибрации. Это отражает возможное несоответ­ ствие режимов при разных пусках. На основании изложенного рекомендовано первый замер производить через 10 мин после установки балансировочной скорости вращения, второй — через 20 мин после первого. Такая рекомендация касается мощных турбоагрегатов, она будет уточняться в процессе более углуб­ ленного изучения описанных явлений.

Вполне понятно, что в разности А' и А" находит одновремен­ но отражение возможная неточность установки балансировочно­ го числа оборотов, а также погрешность измерений.

При всех обстоятельствах значение \А" — А'\ является слу­ чайной величиной, в связи с чем следует оговорить, что она фи­ гурирует в расчетах не как среднеквадратичное значение по­ грешности замера, а вместо него. Эта подстановка оправдана тем, что на расчет влияет не каждое в отдельности значение DA{,

а их сумма (см. образование матрицы [D]). Интересно отметить, что такая подстановка будет неправомерной, если минимизиро­ вать ожидаемое значение em-

Рассмотрим, как влияет учет дисперсий на результаты расче­

та.

Формально учет дисперсий

выражается в сложении

матри­

цы

[а]т [а] с диагональной матрицей [D],

имеющей положитель­

ные члены. Это способствует улучшению обусловленности

матри­

цы [а]т [а], что особенно важно

в случаях, когда последняя плохо

 

мкм

±

^Ps

 

 

го

Ро

 

 

 

 

1,0

 

 

15

 

0,8

 

 

 

 

0,6

10

 

 

 

• \

0,<>

 

 

 

 

0,1

S)

Рис. 1. Кривые зависимостей массы грузов и расчетных остаточных амплитуд вибраций от дисперсии:

 

при Е = / ( — — 1 д л я

пяти

наибольших

по

м о д у л ю

из

д в е н а д -

цати

введенных

в расчет

А •:

А\

•= 28,

Л 2

=

26,

А3

=

36,

А,

=

= 28,

Аъ

= 50;

б

при

— -

/ ( —

) . Р

-

Л, д л я

D =

0

и

D

— D

 

 

д л я

описанного способа назначения диоперсий

 

 

 

обусловлена

и исключает

таким образом

получение

заведомо

абсурдных

результатов.

 

 

 

 

 

 

 

 

 

 

 

Кривые зависимостей массы грузов и расчётных остаточных амплитуд вибрации от дисперсии приведены на рис. 1. При этом можно отметить:

а) относительные изменения веса грузов превышают относи­ тельные изменения амплитуд вибраций;

б) относительные изменения весов грузов существенно раз­ личаются между собой;

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

Графики отражают определенные закономерности, которые выявляются при изучении алгоритма расчета:

а) введение дисперсий существенно снижает значение веса малоэффективных грузов и практически исключает так называе-

мые нечувствительные сочетания грузов, в то же время мало изменяются эффективные грузы;

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

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

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

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

Программа освоена и эксплуатируется несколькими вычисли­ тельными центрами Министерства энергетики и электрификации

СССР на машинах М-220 БЭСМ-ЗМ и БЭСМ-4.

Соседние файлы в папке книги из ГПНТБ