ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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)подсчета относительной погрешности решения системы ли
нейных алгебраических уравнений для различных отклонений ЬА и 6Ь с целью оценки фактического влияния погрешностей ис ходных данных на погрешность результата.
Программа |
выполнена в коде «Минск-2», |
занимает ячейки |
с 100 по 500. Кроме того, для вычисления числа |
обусловленности |
|
симметричных |
матриц сделана более короткая программа, ра |
ботающая с той же БСП. Она занимает ячейки с 100 по 270. С помощью этих двух программ была проведена оценка точно сти решения систем нормальных уравнений, полученных при балансировке натурных многоопорных роторов энергетических турбоагрегатов. В приведенных ниже таблицах представлены полученные экспериментально комплексные значения динамиче
ских коэффициентов влияния атп, |
являющихся элементами |
матриц. |
|
В таблицах приведены данные атп по валопроводам трех турбоагрегатов, состоящих из турбины типа К-300-420 Ленин градского металлического завода с генератором типа ТВВ-320-2 (табл. 1), турбины того же типа Харьковского турбинного завода с генератором ТГВ-300 (табл. 2) и турбины типа К-ЮО-90 с ге нератором типа ТВ-100-2 (табл. 3).