Файл: Теория и практика балансировочной техники..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 283

Скачиваний: 5

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

балансиро­

сечениґі

Номера

вочных

В а р и а н ты у р а в н о в е ш и в а ю щ и х систем Qn в

наград

* *

 

Q,

®эксп

/

 

 

 

 

 

 

 

1,36/260

0,7/85

2

0,57 007

0,14'21

 

0,77/100

0,77/110

0

66 102

0,34/154

3

1,01/36

0,47.348

 

0,18/234

0,18/284

0,55/333

0,3/300

4

2,66 '231

4,33/244

1,65 221

 

 

 

 

1,56/208

 

5

0,68 027

 

 

 

1,02/014

 

 

0,57'134

 

в

 

 

1,59/235

1,78/227

 

1,96 226

2,19/214

1,9/230

Номера точек

1

2

3

4

5

6

7

8

9

10

 

 

^

,

 

.ост

в мкм

 

 

 

 

 

 

Остаточные

вибрации

Ат

 

 

Н а ч а л ь н а я

 

 

 

 

 

 

 

 

А

вибрация

Аост

 

.ост

.ост

 

 

.ост

 

в

мкміград

 

АТ

 

А°ст

.ост

m

 

 

 

л 3

Л 4

 

Л 6

эксп

 

 

17

26

18

6

0

 

6

12

14

 

10/240

26

30

27

15

7

 

13

9

20

 

33/50

15

20

16

28

23

 

17

8

14

 

40/60

21

17

12

10

26

 

8

8

10

 

41/240

11

56

6

12

11

 

15

2

10

 

96 50

12

11

15

10

9

 

12

10

18

 

15/60

21

24

12

15

17

 

11

2

24

 

17/320

7

о

14

38

32

 

30

12

17

 

29/160

39

42

21

16

29

 

12

14

22

 

39/170

28

28

18

21

56

 

27

5

18

 

92170

*

Оптииіальньїй

вариані - вида

Q я '

* *

О п т и \іальньїй

вариант вида

Q

Л2'

 

 

 

 

Выводы

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

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

Сравнение результатов расчета и эксперимента дало хоро­ шее совпадение.

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


Е. А. ГАЛЬПЕРИН, С. И. МИКУНИС, Б. О. MAP ДЕР

ОЦЕНКА ТОЧНОСТИ РЕШЕНИЯ ЗАДАЧИ УРАВНОВЕШИВАНИЯ МНОГООПОРНЫХ РОТОРОВ С ПРИМЕНЕНИЕМ ЭЦВМ

Обработка на ЭЦВМ информации, получаемой при баланси­ ровке однотипных агрегатов, требует решения систем линейных алгебраических уравнений, в которых число уравнений значи­ тельно больше числа неизвестных. Как правило, такие системы несовместны и не имеют точного решения. Приближенное реше­ ние по методу наименьших квадратов сводится к решению си­ стемы линейных уравнений с квадратной матрицей [1], [2]. Одна­ ко в процессе решения возникают трудности, связанные с воз­ можностью плохой обусловленности матрицы системы нормаль­ ных уравнений. Число обусловленности дает оценку того, на­ сколько относительная погрешность результата превосходит погрешность исходной информации. Если число обусловленности велико, то небольшая ошибка в исходных данных приводит к значительным ошибкам в решении. Поэтому оценка обусловлен­ ности матрицы дает существенную характеристику качества ре­ шения.

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

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

предназначен­

ных для решения задачи многоплоскостной

динамической

балансировки гибких валопроводов.

 

Рассмотрим

подробнее проблему обусловленности для си­

стем линейных

алгебраических уравнений с

невырожденной

матрицей

 

 

 

Ax = b,

(1)

где х — n-мерный вектор неизвестных величин;

bзаданный /г-вектор;

А— заданная п X n-матрица, определитель которой \А \ Ф

Пусть

\ = 1. ЭТОГО

всегда

можно

добиться изменением

масштаба одной из неизвестных величин.

 

 

 

Исходные

данные

— элементы

атп, Ьт

— являются

резуль­

татом предыдущих

вычислений или

непосредственных

измере­

ний, или могут быть параметрами

некоторой

приближенной

модели.

Они

всегда

содержат

определенную

погрешность.

Эта погрешность

служит

одним

из

источников

погрешности

решения.

Рассмотрим случай идеального вычислителя, когда машина получает решение х = А~ХЬ совершенно точно. Тем самым выде­ ляется та часть погрешности решения, которая не зависит от


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

Таким образом, A, b, х в выражении (1) означают номиналь­ ные (заданные или вычисляемые) величины. Отклонения 6Л, 6Ь, 6х принципиально ненаблюдаемы, поэтому истинные значе­

ния А + ЬА, b +

6b, х + бх неизвестны и не могут быть

вычисле­

ны.

Через

IUH обозначим

любую подходящую

норму

вектора,

IIЛ || — соответствующая норма матрицы [2].

 

 

 

 

 

 

Записывая соотношение (1) для истинных значений

и

рас­

крывая

скобки, имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ах + АЬх + ЬАх + ЬАЪх^Ь

+ ЬЬ.

 

 

 

 

(2)

 

Вычитая

соотношение

(1)

из выражения (2), получим урав­

нение для

отклонений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 х = А~хЬАх

+ A~xbb

Л^'бЛбд;.

 

 

 

(3)

Переходя

в

выражении

(3)

к нормам

и учитывая

свойства

нормы [2], имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1

1

1 "

 

11

11

11

"

11

п л ||

 

 

1

 

" | | Л | | - | | * | |

 

 

 

 

 

 

+

1И~'II-II6^11-II6*II-

" 4 Н * Н .

 

 

 

(4)

 

 

 

 

Т

"

 

" "

" "

"

| | Л | | . | | * | 1

 

 

 

 

v '

Так

как

\\Ах\\

^

||Л|| • ||х||, то, заменяя в

знаменателе

второго

слагаемого ЦЛЦ-IWI величиной \\Ах\\

— \\Ь\\,

лишь усилим

нера­

венство

(4). Обозначим

\\А~Ц\ • ||Л|| — с(А)

и,

разделив

(4)

на

IUII, окончательно

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

II

 

II < с (

Л

) {

I I м

II

|

Н^Н |

I I м

II Н^Н

\

 

/кч

 

 

II

II

'

' I

II

Л II

 

 

и U II

 

II

Л II

II .. II

)

'

 

'

Если

 

погрешность

имеется

лишь

в

правых

частях

системы

(1), но не в матрице Л,

то, полагая

в

неравенстве

(5)

бЛ = О,

получим неравенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ж

< с

(

л ) Ж .

 

 

 

 

 

 

 

( 6 )

Величина

с (Л)

= ЦЛ|| • ||Л_ 1 ||

называется

числом

обусловлен­

ности матрицы Л и характеризует степень

искажения единичной

окружности после преобразования у = Ах.

 

В эвклидовой

норме

(||х|| д л и н а вектора

х)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С(А)

=

lAmaxAmin.

 

 

 

 

 

 

 

где Атах и АПІІП

максимальное

и минимальное

собственные

чи­

сла симметричной

матрицы Л • Л'.

 

 

 

 

 

 

 

 

 

152


Если ||6Л|| # 0 и величина

p = i _ c W ) l M I L > 0 >

\\а\\

то, разрешая неравенство (5) относительно ||бх]|/||х||, получим гарантированную оценку

||бхЦ

 

с(А)

/

| [ М | |

\\ЬЬ\\

\

 

(7)

11*11

^

Р

V

И Л ||

Ц&ІІ

у

 

 

 

При больших значениях с (А)

будем иметь р ^

0; в этом

слу­

чае оценка (7) несправедлива

и приходится ограничиться

при­

ближенной оценкой, получаемой из неравенства

(5) отбрасыва­

нием члена с(А)-Ц&4Ц

• ||6*||/||Л|| •

второго

порядка

ма­

лости.

число обусловленности с (А)

 

 

Таким образом,

показывает, во

сколько раз может возрасти относительная погрешность резуль­ тата по сравнению с относительной погрешностью исходных данных в случае идеального вычислителя. Расчеты на ЦВМ по­ казали, кроме того, что при больших с (А) вычисление обратной матрицы Л - 1 по существующим стандартным программам дает неудовлетворительный результат [проверено при с(А) = 10 000 для матриц 12-го порядка]. Поэтому при решении систем линей­ ных уравнений необходимо учитывать возможность плохой обу­ словленности матрицы системы.

Программа, разработанная для ЭЦВМ «Минск-2», выполня­ ет следующие операции:

1)подсчета числа обусловленности матрицы А;

2)подсчета относительной погрешности решения системы ли­

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

Программа

выполнена в коде «Минск-2»,

занимает ячейки

с 100 по 500. Кроме того, для вычисления числа

обусловленности

симметричных

матриц сделана более короткая программа, ра­

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

ских коэффициентов влияния атп,

являющихся элементами

матриц.

 

В таблицах приведены данные атп по валопроводам трех турбоагрегатов, состоящих из турбины типа К-300-420 Ленин­ градского металлического завода с генератором типа ТВВ-320-2 (табл. 1), турбины того же типа Харьковского турбинного завода с генератором ТГВ-300 (табл. 2) и турбины типа К-ЮО-90 с ге­ нератором типа ТВ-100-2 (табл. 3).