книги / Новые подходы к исследованию и идентификации переходных процессов синхронных машин
..pdfИсследования по влиянию величины шага дискретизации на результаты идентификации в опытах ГП с длительными ПП представлены в табл. 39 (для СМ в режиме холостого хода без гасительного сопротивления).
Анализ результатов исследований в табл. 39 показывает, что укрупнение шага дискретизации в длительных ПП в опытах ГП с разомкнутой ОЯ без гасительного сопротивления наблюдается постепенное нарастание погрешностей приближения модели составляющих и всего ПП к их опытным данным. Это объясняется преобладанием зашумлённых опытных данных в конце затухания ПП, снижающих общую погрешность приближения.
Таблица 39
Результаты исследования влияния шага дискретизации дискретного переходного процесса на его идентификацию в опытах ГП
|
Шаг дискретизации |
|
|
|
|
Результаты иденти- |
|
|
|
Результаты идентифи- |
|
|
|
|
||||||||||
ОЯ, |
2U |
|
, |
|
фикации переходной |
|
|
′ , |
кации сверхпереходной |
|
|
|
, |
|||||||||||
Rг |
|
∞ |
|
|
|
составляющей |
|
|
|
K |
|
составляющей |
|
|
|
K |
|
|||||||
о.е |
|
|
|
|
|
|
о.е. |
|
|
|
|
|
|
о.е. |
||||||||||
t′ |
|
– |
τ′опт , с |
|
2U ′ , |
t′′ |
|
– |
|
|
2U′′ , |
|||||||||||||
|
|
|
|
|
н.эф.опт |
|
|
|
|
|
н.эф.опт |
|
τ′′ , с |
|
|
|
|
|
||||||
|
|
|
|
|
tв′.эф.опт , с |
|
0 |
|
|
|
tв′′.эф.опт , с |
|
0 |
|
|
|
|
|||||||
|
|
|
|
|
|
о.е. |
|
|
|
|
о.е. |
|
|
|
|
|||||||||
|
0,01 |
18,75 |
1,435–10,585 |
8,245528 |
267,803 |
0,40596 |
0,005–0,625 |
0,529636 |
|
1,853 |
0,41711 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
K=993 |
|
|
|
|
|
|
K=1063 |
|||||
ХХ, |
0,02 |
17,87 |
1,325–10,525 |
8,278985 |
268,689 |
0,43998 |
0,005–0,625 |
0,517249 |
|
1,847 |
0,44929 |
|||||||||||||
Rг |
|
|
|
|
|
|
|
|
|
|
|
K=497 |
|
|
|
|
|
|
K=532 |
|||||
0,04 |
19,49 |
|
0,805–9,925 |
8,189238 |
267,560 |
0,46156 |
0,005–0,685 |
0,586783 |
|
1,350 |
0,47001 |
|||||||||||||
нет |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
K=249 |
|
|
|
|
|
|
K=267 |
|||||
|
0,08 |
19,49 |
|
0,805–9,925 |
8,178962 |
267,757 |
0,52745 |
0,005–0,565 |
0,441473 |
|
1,156 |
0,53126 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
K=125 |
|
|
|
|
|
|
K=134 |
4.5. ИССЛЕДОВАНИЕ ВЛИЯНИЯ ШАГОВ ВАРИАЦИИ УСТАНОВИВШЕГОСЯ ЗНАЧЕНИЯ ТОКА (НАПРЯЖЕНИЯ) ЯКОРЯ
И НИЖНЕЙ ГРАНИЦЫ ПЕРЕХОДНОЙ СОСТАВЛЯЮЩЕЙ НА
ИДЕНТИФИКАЦИЮ ПЕРЕХОДНОГО ПРОЦЕССА
Во второй главе монографии в подразд. 2.4.3.7 представлена разработка нового оригинального метода минимизации среднеквадратичной погрешности приближения модели переходной со-
171
ставляющей симметричного тока якоря СМ не только для ПП
вопыте ВКЗ, но и для остальных рассматриваемых опытов в данной монографии. Апробация нового оригинального метода минимизации среднеквадратичной погрешности приближения модели переходной составляющей симметричного тока якоря (напряжения в других опытах) к её опытным данным в исследуемых диапазонах ПП с представлением погрешности в виде поверхности
втрёхмерной системе координат успешно подтверждена на опытных данных для всех опытов в 3-й главе данной монографии.
Вданной главе исследованы возможности разработанного метода в направлении сохранения его точности при изменении
шага вариации установившегося значения тока (напряжения) с одновременным изменением шага вариации нижней границы исследуемого диапазона ПП в различных опытах с целью оценки машинного времени минимизации погрешности приближения вместе с построением её в виде поверхности в трёхмерной системе координат.
Исследования, проводимые по идентификации опыта ГП (при разомкнутой ОЯ без гасительного сопротивления), подтвердили, что при оптимизации установившегося значения напряжения (тока) якоря и нижней границы участка переходной составляющей шаг вариации указанных параметров оказывает существенное влияние на время их идентификации при построении дискретной поверхности среднеквадратичной погрешности приближения модели ПП к её опытным данным.
В табл. 40 и на рис. 54 представлены результаты исследования влияния на время обработки ПП в опытах ГП отдельно: шага установившегося значения (верхний ряд на рис. 54) и нижней границы переходной составляющей (нижний ряд на рис. 54).
Анализ результатов в табл. 40 показывает, что при увеличении шага вариации время идентификации ПП значительно сокращается, с практически небольшим разбросом параметров идентификации переходной составляющей, что подтверждается стабильностью величины среднеквадратичной погрешности приближения.
172
Таблица 40
Результаты исследования влияния шагов вариации установившегося значения тока (напряжения) якоря и нижней границы переходной составляющей на идентификацию ПП в опытах ГП
Варьируемая |
Шаг |
2I∞ , |
τ′ , с |
2I′ , |
. . |
|
|
|
′ , |
. . |
Время |
|||||||||
|
|
|
|
|
|
|||||||||||||||
|
2I |
|
или |
t′ |
вари- |
о.е. |
0 |
о е |
K |
|
|
K |
о е |
иденти- |
||||||
|
|
|
∞ |
|
|
|
н |
ации |
|
|
|
|
|
|
|
|
фикации |
|||
|
2I∞ var , о.е. |
0,05 |
17,15 |
8,340080 |
265,09 |
1051 |
0,347746 |
8 мин |
||||||||||||
|
1 сек |
|||||||||||||||||||
t′ |
|
|
|
|
|
|
|
|
|
|
||||||||||
= |
0,01 c |
= |
onst |
0,5 |
17,0 |
8,354930 |
264,99 |
1051 |
0,348383 |
58 |
сек |
|||||||||
н |
|
|
|
|
с |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
5 |
15 |
8,476090 |
266,591 |
1051 |
0,374827 |
13 сек |
||||||
|
|
|
t′ |
, с |
|
0,01 |
17 |
8,354930 |
264,99 |
1051 |
0,348383 |
2 мин |
||||||||
|
|
|
|
н.var |
|
|
|
10 сек |
||||||||||||
2I∞ = 0,2 о.е. = |
|
|
|
|
|
|
|
|
|
|
||||||||||
0,4 |
17 |
8,354930 |
264,99 |
1051 |
0,348383 |
36 сек |
||||||||||||||
|
|
= сonst |
|
|||||||||||||||||
|
|
|
0,1 |
16,6 |
8,373220 |
265,454 |
1051 |
0,350224 |
19 сек |
Исходные данные при проведении исследований оставались неизменными, соответственно объем элементов оставался постоянным. При данных исследованиях диапазон минимизированного объёма ЭТВ ( nэф = 4) автоматически перемещался по исследуе-
мому диапазону переходной составляющей.
В ходе построения поверхности дискретной погрешности приближения алгоритмически определяются координаты минимизированной погрешности с соответствующими параметрами идентификации переходной составляющей симметричного тока (напряжения).
По рис. 54 видно, что при увеличении шага вариации дискретная поверхность среднеквадратичных погрешностей приближения модели к опыту становится менее лохматой (более гладкой) за счёт разрежения построения погрешности приближения зашумлёнными кривыми статистической функции (40) (см. главу 2).
Аналогичные исследования проведены по результатам стендовых испытаний в опытах ВКЗ, ВН, УВ с учетом их особенностей.
Необходимо отметить заранее известную область построения поверхности среднеквадратичной погрешности приближения фор-
173
мируемой модели к опытным данным, связанную с варьированием нижней границы исследуемого диапазона переходной составляющей tн′ − tв′ при изменении установившегося значения от 0 до не более удвоенного значения установившегося тока (напряжения).
Времясчета: 8 мин 1 сек |
Времясчета: 58 сек |
Времясчета: 13 сек |
Времясчета: 2 мин 10 сек |
Времясчета: 36 сек |
Времясчета: 19 сек |
Рис. 54. Результаты исследования влияния на время обработки ПП в опытах ГП шага установившегося значения (верхний ряд)
инижней границы переходной составляющей (нижний ряд)
Витоге область построения поверхности определяет способы варьирования установившегося значения и нижней границы. Использование при этом мелкого шага гарантирует получение неровной, шероховатой поверхности с неярко выраженным «жело-
бом», изобилующей локальными экстремумами (минимумами и максимумами), поэтому найденный на поверхности минимум среднеквадратичной погрешности будет абсолютным в рассматриваемой области, но с большим временем поиска. Использование укрупненного шага позволяет получить гладкую поверхность с ярко выраженным «желобом» минимальных погрешностей, уменьшает время поиска экстремума, но при очень крупных шагах абсолютный минимум может быть утерян (он может оказаться
174
за пределами применяемого шага). Таким образом, методика позволяет изменять величину шага от крупного (на краю исследуемой области, т.е. на склоне «желоба», где заведомо абсолютный экстремум отсутствует) до мелкого (по мере приближения к дну «желоба»), причем шаг можно уменьшать по мере приближения к абсолютному минимуму. Методика также позволяет варьировать установившееся значение и нижнюю границу с различными независимыми друг от друга шагами, а также позволяет применять другие известные способы изменения шага.
По результатам представленных исследований, связанных с минимизацией среднеквадратичной погрешности приближения
вцелом, напрашивается обобщённое резюме по представленному новому оригинальному методу минимизации указанной погрешности.
Разработанный новый оригинальный метод минимизации среднеквадратичной погрешности приближения модели переходной составляющей симметричного тока (напряжения) якоря с визуальным представлением данной погрешности в виде поверхности в системе трёх координат по опытным зашумлённым данным
влюбых симметричных и несимметричных затухающих (возрастающих) синусоидальных ПП СМ позволяет с наивысшей точностью и высокой достоверностью идентифицировать ПП по основному объёму используемой опытной информации на исследуемом
участке tн′ − tв′ . Это происходит за счёт существования в ПП СМ переходной составляющей, занимающей основной объём ПП от его начала и до его установившегося значения. Высокая точность
идостоверность идентификации переходной составляющей в исследованных в 3-й главе различных ПП МСМ обусловлены использованием в разработанном новом методе операции моделирования переходной составляющей минимизированным объёмом ЭТВ из четырёх штук для идентификации данной составляющей, которые по распределению Пуассона отклоняются от МО с минимальной относительной погрешностью (от долей % до 1–2 %)
ивысокой вероятностью их существования до 0,9980 в зашумлённых ПП МСМ (подразд. 2.4.3.7, глава 2).
175
Глава 5. ВОЗМОЖНОСТИ ВЕРОЯТНОСТНО-СТАТИСТИЧЕСКИХ МЕТОДОВ
Разработанные ВСМ могут быть эффективно использованы
вдвух крупных направлениях:
–для исследования и идентификации ПП единичных МСМ
или для малочисленных партий в полном объёме разработанных теоретических основ нового подхода к исследованиям и идентификации данных процессов в опытах ВКЗ, ВН, ГП, УВ;
–для идентификации ПП СМ меньшей мощности при их серийном изготовлении, используя основы ВСМ в объёме экс- пресс-обработки результатов стендовых испытаний.
Первое направление достаточно подробно рассмотрено в данной монографии на примерах МСМ. Второе направление рассмотрено ниже в варианте экспресс-обработки ПП СМ мощностью 800 кВт в опыте ВКЗ.
Рассмотрим наиболее важные для исследования ПП СМ вопросы, публикации по которым практически отсутствуют в печати.
Определение ударного тока из опытов ВКЗ по стандартам имеет много допущений, например, определение ударного тока через 0,01 с от начала ПП, которое требует дополнительных исследований по уточнению данного требования. В данной главе рассмотрен новый подход к достоверному определению ударного тока якоря СМ. При этом возникает много сопутствующих проблем.
Проблема учёта влияния случайных факторов на результа-
ты обработки опыта ВКЗ СМ независимо от уровня их мощности:
–точное использование установившегося значения тока якоря в опытах ВКЗ;
–начальный сдвиг первой вершины тока якоря;
–сдвиг первых пяти-шести вершин тока якоря при сильно выраженной сверхпереходной составляющей в переходном процессе;
–использование эффективных точечных выборок;
–ряд других факторов, практически сопровождающих процесс обработки опыта ВКЗ, например, изменение напряжения
176
возбуждения, напряжения испытания, частоты напряжения питания и т.д.
Интересен вопрос по определению величины расхождения индуктивных переходных и сверхпереходных сопротивлений для ТГ, так как в теории ПП СМ они допускаются примерно одинаковыми, в связи с чем возникает задача уточнения величины отклонения указанных индуктивных сопротивлений по оси q из опыта ВКЗ.
Влияние изменения шага квантования дискретных ПП СМ на результаты исследования и идентификации данных процессов при испытаниях машин в опытах ВКЗ, ВН, ГП, УВ.
5.1. ИДЕНТИФИКАЦИЯ ППИЗ ОПЫТА ВКЗДЛЯ СМ
МОЩНОСТЬЮ 800 КВТ МЕТОДОМ ЭКСПРЕСС-ОБРАБОТКИ
Источники дополнительно к указанным выше факторам случайного происхождения вызывают существенное отклонение реальных ПП от идеального (расчётного) экспоненциального (или показательного) закона затухания, поэтому задачей метода экс- пресс-обработки является оценка реальной погрешности отклонения от экспоненциальности затухающих ПП при серийных испытаниях СМ.
Методы фильтрации, с отбраковкой существенно отклоняющихся значений тока якоря от экспоненциальности, как это предлагается в отдельных публикациях, в принципе не применимы для СМ с малой длительностью ПП. Использование в экспресс-методе случайного признака в виде kj-й ПВ по (20) позволяет не только создавать с некоторой избыточностью генеральную совокупность случайного признака и выборки из неё (даже в случае коротких переходных процессов), но и эффективно упорядочивать данные совокупности случайного признака с полным сохранением исходной (первичной) информации. Это существенно повышает надёжность и достоверность результатов идентификации ПП в условиях влияния на них случайных факторов. Например, для реальных ПП в опыте ВКЗ СМ указанной мощности удаётся подтверждать нормальный закон распределения отклонений случайного признака
177
от МО. Удаётся качественно улучшать уровень идентификации ПП за счёт ЭТВ, легко выделяемых непосредственно из ограниченных объёмов генеральной совокупности, что придаёт исследованиям инженерный характер. Но это становится возможным только с использованием методов ТВ и МС и методов математического моделирования.
Разработанный на вероятностной основе экспресс-метод обработки и исследования переходного процесса в опыте ВКЗ, пригодный также к исследованиям и обработке аналогичных процессов в таких опытах, как ГП, ВН, УВ и др., по праву отнесён к ВСМ.
В процессе обработки опыта ВКЗ учёт погрешности от влияния случайных факторов является главнейшей возможностью метода экспресс-обработки ПП для повышения достоверности результатов идентификации. Исследованиями установлено, что основным источником появления систематической погрешности в процессе обработки является неточное использование установившегося значения тока якоря в опытах ВКЗ. Эта проблема отмечается в научных публикациях давно, но эффективного её решения пока не найдено. Из теории ПП СМ в опытах ВКЗ известно, что наибольшим временем затухания из всех составляющих тока якоря обладает переходная составляющая [5]. При этом также важен выбор участка переходной составляющей, по которому следует оптимизировать установившееся значение тока якоря.
В главе 2 монографии описан новый оригинальный метод минимизации погрешности приближения модели переходной составляющей в исследуемом диапазоне ПП к его опытным данным, эффективно подтверждённый для МСМ в главе 3.
При экспресс-обработке менее мощных СМ из-за ограниченной длительности протекания ПП предложен упрощенный вариант оптимизации установившегося значения тока якоря в опытах ВКЗ, подтверждённый резким снижением среднеквадратичной погрешности приближения для всего ПП.
Рассмотрим влияние установившегося значения тока якоря на результаты идентификации ПП по результатам стендовых испытаний СМ в опыте ВКЗ двумя методами: ГАМ (по стандарту) и ВСМ.
178
Расчёт ПП синхронного турбодвигателя экспресс-методом
В процессе стендовых испытаний синхронного турбодвигателя СТД-800-2 с линейным напряжением испытания на статоре Uст = 0,75Uн = 4512 В при Uн = 6 кВ и мощностью на 800 кВт из
опыта ВКЗ ПИ получена осциллограмма.
Исходными данными для обработки опыта ВКЗ представлены обмеры осциллограммы, полученные заводом в процессе её обработки в соответствии с ГАМ по отечественным стандартам для одной из трёх фаз.
В табл. 40 в колонках 1,2 приведены с шагом 0,02 с измеренные заводом-изготовителем положительные и отрицательные вершины полного тока якоря с начальным фазовым сдвигом первой вершины T = 0,006 с. Установившееся значение тока якоря, измеренное заводом изготовителем, равно 45,3 А (8,2 мм). На рис. 55, а показан дискретный ПП в виде затухающих вершин с шагом 0,02 с, а на рис. 55, б приведён полный дискретный ПП в виде дискретных элементов с шагом 0,01 с. Элемент на первом шаге ПП рассчитан по формулам (5) в главе 2, а элементы со второго шага до конца процесса рассчитаны по формуле (6). Огибающие представляют симметричный ток якоря в виде элементов по (14), (15), ассиметричный ток (среднее значение) получен с учётом (3) по (10), которые в табл. 41 приведены в колонках 5 и 7. В табл. 41 (колонка 10) в соответствии с подразд. 2.4.1; 2.5; 3.6 и рис. 16 теоретических основ нового подхода обоснованы и заданы границы ПП для разделения и идентификации соответствующих составляющих. С помощью ЭТВ по формуле (34) с условием (33) и по формулам (38), (39) идентифицирована переходная составляющая.
Параметры сверхпереходной составляющей определяются тоже с использованием ЭТВ по формулам (41)–(43), для кото-
рых верхняя граница определяется также по условию (33): iв′′(tв′′) = 0,33iн′′(tн′′) .
Идентификация асимметричной составляющей осуществлена по формулам (44) также с использованием ЭТВ.
179
Таблица 41
Исходные данные и промежуточные вычисления по опытным данным ВСМ
Опытные |
Ток якоря между |
Массивы |
Массив |
Модель |
|
||||||||
составляющих |
периодического |
|
|||||||||||
вершины |
огибающими |
периодического |
апериод. |
тока якоря по |
|
||||||||
тока якоря |
|
по (5), (6) |
тока якоря |
тока |
(15) |
Приме- |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
чание |
|
Шаг |
± I , |
i ,мм |
i , мм |
iо′j , мм |
i′′ , мм |
iаj , мм |
ГАМ |
ВСМ |
|||||
|
|||||||||||||
t , |
mj |
по |
оj |
по |
|
||||||||
с |
j |
мм |
иj |
|
оj |
|
по (41) |
|
|
|
|||
|
|
|
|
|
|
(3), (14) |
|
(3), (10) |
|
|
|
||
1 |
|
2 |
|
3 |
|
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
0,006 |
66,0 |
7,1125 |
73,1120 |
|
27,2886 |
29,4440 |
68,8605 |
73,1112 |
T |
||||
0,016 |
–7,9 |
52,2750 |
60,1750 |
|
17,4959 |
22,1875 |
57,5904 |
60,2356 |
tн′′1 =tа.н |
||||
0,026 |
41,7 |
8,5625 |
50,2625 |
|
10,4650 |
16,5687 |
49,4907 |
51,0930 |
|
||||
0,036 |
–9,1 |
35,3250 |
44,4250 |
|
7,2682 |
13,1125 |
43,4857 |
44,4241 |
|
||||
0,046 |
30,0 |
9,5125 |
39,5125 |
|
4,7758 |
10,2437 |
38,8854 |
39,4124 |
|
||||
0,056 |
–9,8 |
25,9250 |
35,7250 |
|
3,2061 |
8,0625 |
35,2451 |
35,5271 |
tа.в |
||||
0,066 |
22,5 |
9,9250 |
32,4250 |
|
1,9385 |
6,2875 |
32,2768 |
32,4219 |
|
||||
0,076 |
–10,0 |
19,8750 |
29,8750 |
|
1,2511 |
4,9375 |
29,7923 |
29,8691 |
|
||||
0,086 |
17,6 |
10,0375 |
27,6375 |
|
0,7205 |
3,7812 |
27,6669 |
27,7181 |
|
||||
0,096 |
–10,0 |
15,7063 |
25,7063 |
|
0,3536 |
2,8531 |
25,8172 |
25,8681 |
|
||||
0,106 |
14,1 |
9,8875 |
23,9875 |
15,7875 |
0,0683 |
2,1062 |
24,1859 |
24,2508 |
tн′1 |
||||
0,116 |
–9,7 |
12,8375 |
22,5375 |
14,3375 |
|
1,5687 |
22,7331 |
22,8188 |
tн′ 2 |
||||
0,126 |
11,75 |
9,3875 |
21,1375 |
12,9375 |
|
1,1812 |
21,4297 |
21,5387 |
|
||||
0,136 |
–9,1 |
10,8813 |
19,9813 |
11,7813 |
|
0,8906 |
20,2544 |
20,3865 |
|
||||
0,146 |
10,1 |
8,8500 |
18,9500 |
10,7500 |
|
0,6250 |
19,1906 |
19,3439 |
|
||||
0,156 |
–8,6 |
9,4000 |
18,0000 |
9,8000 |
|
0,4000 |
18,2252 |
18,3970 |
|
||||
0,166 |
8,8 |
8,3750 |
17,1750 |
8,9750 |
|
0,2125 |
17,3475 |
17,5348 |
|
||||
0,176 |
–8,1 |
8,3250 |
16,4250 |
8,2250 |
|
0,1125 |
16,5484 |
16,7482 |
|
||||
0,186 |
7,9 |
7,7000 |
15,6000 |
7,4000 |
|
0,1000 |
15,8203 |
16,0297 |
|
||||
0,196 |
–7,4 |
7,5188 |
14,9188 |
6,7188 |
|
0,0594 |
15,1564 |
15,3727 |
|
||||
0,206 |
7,2 |
7,2750 |
14,4750 |
6,2750 |
|
0,0375 |
14,5508 |
14,7715 |
|
||||
0,216 |
–7,1 |
6,9625 |
14,0625 |
5,8625 |
|
0,0687 |
13,9983 |
14,2212 |
|
||||
0,226 |
6,75 |
6,8125 |
13,5625 |
5,3625 |
|
0,0312 |
13,4940 |
13,7173 |
|
||||
0,236 |
–6,6 |
6,5688 |
13,1688 |
4,9688 |
|
0,0152 |
13,0336 |
13,2557 |
tв′1 |
||||
0,246 |
6,4 |
6,5125 |
12,9125 |
4,7125 |
|
0,0562 |
12,6134 |
12,8329 |
tв′2 |
||||
0,256 |
–6,4 |
6,2563 |
12,6563 |
4,4563 |
|
0,0718 |
12,2298 |
12,4455 |
|
||||
0,266 |
6,1 |
6,2250 |
12,3250 |
4,1250 |
|
0,0625 |
11,8795 |
12,0905 |
|
||||
0,276 |
–6,1 |
5,9000 |
12,0000 |
3,8000 |
|
0,1000 |
11,5597 |
11,7653 |
|
||||
0,286 |
5,75 |
6,0875 |
11,8375 |
3,6375 |
|
0,1687 |
11,2677 |
11,4673 |
|
180