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

книги / Новые подходы к исследованию и идентификации переходных процессов синхронных машин

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
17.76 Mб
Скачать

Таблица 20

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

 

Диапазон

 

 

Результаты оптимизации

исследований ПП,

величин 2U

,

tпо опытным

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

 

 

 

 

 

н

 

 

 

 

 

данным с помощью

значение напряже-

моделирования переходной

ния якоря, нижняя

 

составляющей в узлах

граница, перед

 

дискретизации в диапазоне

оптимизацией

 

 

исследования ПП

 

*

2U

 

,

t

 

 

2U

 

,

 

 

,

 

 

 

 

 

 

t

н.эф.опт

 

 

∞. опт

 

 

t, с

 

 

tв.эф.опт , с

 

 

 

 

 

 

K

н

в

о.е.

 

 

 

о.е

 

о.е.

 

 

 

 

 

 

0,505–

225,086

0,315–

 

225,69

 

0,394524

12,305

 

 

 

8,845

 

 

 

 

 

K=1181

Результаты проверки оптимизации величин 2U, tнна модели

переходной

составляющей

 

 

 

 

 

 

 

 

t

2U

 

,

 

 

,

 

 

 

н.эф.опт

 

∞. опт

 

 

tв.эф.опт , с

 

 

 

 

K

о.е

 

о.е.

1,005–

225,69

 

0,000014

9,605

 

 

 

K=1181

3.7.3.8. Идентификация сверхпереходной составляющей симметричного напряжения обмотки якоря

Идентификация сверхпереходной составляющей ПП СМ осуществляется в полном соответствии с п. 2.5 главы 2. Результаты идентификации параметров сверхпереходной составляющей ПП представлены в табл. 21.

Таблица 21

Результаты идентификации параметров сверхпереходной составляющей

t′′

t′′

, с

τ′′

 

, с по (36)

2U′′

,

. .

по

(37)

н.эф

в.эф

 

эф

 

0.эф

 

о е

 

0,005–0,055

 

0,034848

 

1,7555

 

3.7.3.9. Критериальная оценка точности идентификации переходного процесса

Исследования опыта ВН показывают, что наилучшей оценкой подтверждения точности и достоверности идентификации ПП

141

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

 

 

 

 

 

 

 

 

1

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K =

 

[(uоk )мод (uоk )оп]2

,

 

(56)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

где (uоk )мод

– математическая модель дискретно заданного полно-

,

(u

 

)

мод =

2U ′′

 

е

(t′′

)/τ′′

+

2U

 

e

(t)

 

 

получаемая на

оk

 

н.эф

эф

эф

н.эф

опт

 

эф ,

го ПП

 

 

 

0эф

 

 

 

 

0

 

 

 

 

 

 

базе ЭТВ с учётом оптимизации установившегося значения напряжения якоря (2U)опт и нижней границы (tн.эф )опт в диапазоне

tн′′ − tвПП; (uоk )оп – дискретные элементы ПП, полученные по опытным данным (3) в диапазоне tн′′ − tвПП, преобразованного из возрастающего ПП в убывающий. Для учёта первого элемента ПП в формуле модели расчёт составляющих производится для времени T , а опытное значение рассчитывают по выражениям (3). В результате идентификации ПП в опыте ВН погрешность приближения модели симметричного напряжения якоря СМ к опытным данным ПП в диапазоне tн′′ − tв, определяемая по формуле

(56), составила 1231 = 0,500751.

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

142

развития метода с использованием при модернизации и проектирования ЭМ по результатам конструкторских или исследовательских испытаний.

3.8. ОПЫТ ГПИ УВ

Опыты ГП и УВ вызывают нарушения электромагнитного состояния СМ, в результате которых в обмотках якоря и индуктора в зависимости от состояния ОЯ в ней возникают синусоидальные затухающие или возрастающие ПП, состоящие из симметричного тока или напряжения якоря с переходной и сверхпереходной составляющей. Методика разделения данных составляющих с целью их идентификации полностью соответствует описанной для ПП СМ в опытах ВКЗ. ПП в обмотке индуктора в случае затухающей или возрастающей кривой по экспоненциальному закону после их идентификации дают электромагнитные ПВ, практически совпадающие с ПВ, получаемые при идентификации ПП из ОЯ [7].

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

В разделе рассмотрен новый подход к исследованию и идентификации ПП в указанных опытах в полном соответствии с алгоритмом, рассмотренным во 2-й главе для опыта ВКЗ, и распространенной методикой исследования с идентификацией ПП в опыте ВН на опыт УВ.

143

3.8.1. Предварительная обработка первичной информации

Объектом исследования и идентификации ПП в опыте ГП была избрана та же МСМ на 50 МВт. Исходные опытные данные для апробации нового подхода исследования и идентификации ПП получены в процессе её стендовых испытаний с разомкнутой ОЯ (режим холостого хода) без гасительного сопротивления Rг с использованием ЦЗО для представления дискретного ПП

(рис. 7), (табл. 1, прилож. 3).

После предварительной и последующей обработки информации ПП представлен дискретным в виде выделенных вершин с шагом дискретизации 0,02 с (рис. 41) [87] (табл. 2, прилож. 3).

Рис. 41. Выделенные вершины дискретного ПП

На рис. 42 ПП представлен в виде дискретных элементов между огибающими с шагом дискретизации 0,01 с в соответствии с подразд. 2.3.1–2.3.4 (табл. 2, прилож. 3) или в виде экспоненциально затухающих элементов одной полярности (рис. 43), (см. табл. 2, прилож. 3).

Рис. 42. Огибающие в виде дискретных элементов ПП

144

Рис. 43. Представление дискретного ПП в виде экспоненциально затухающих элементов

3.8.2. Задание границ в дискретном переходном процессе

Границы, связанные с обоснованием диапазона ПП для исследования переходной составляющей, затем сверхпереходной, заданы в полном соответствии с подразд. 2.4.1 и рис. 16 (глава 2).

Результаты выбора границ представлены в табл. 22. Таблица 22

Границы диапазона исследований, установившиеся значения напряжения якоря

 

 

 

 

 

 

 

ОЯ

Rг

T ,

tн′′ = T ,

tн′ − tв* ,

2u( T ) ( 2i( T ) ),

2U( 2I),

с

с

с

о.е.

о.е.

 

 

ХХ

нет

0,005

0,005

0,705–10,625

288,2266

8,6468

 

есть

0,005

0,005

0,505–4,495

291,0629

8,7319

КЗ

нет

0,005

0,005

0,405–3,255

0,9013

0,0270

 

есть

0,005

0,005

0,155–3,025

0,9325

0,0280

Примечание: ХХ – режим холостого хода, при этом ОЯ разомкнута, ЦЗО записывает напряжение ОЯ; КЗ – режим короткого замыкания, ОЯ замкнута, ЦЗО записывает ток ОЯ; установившееся значение в колонке 7 согласно [7] взято в размере 3 % от первоначального.

В подразд. 2.4.2 (глава 2) монографии представлено обоснование базового случайного признака для синусоидально затухающих (возрастающих) по экспоненциальному закону ПП в виде полученного выражения (20), которое позволяет оценить качество исходных опытных статистических данных стендовых испытаний и осуществить высокоточную идентификацию ПП МСМ.

145

3.8.3. Проведение статистических исследований

Исследования качества опытных данных проводят в порядке, рассмотренном в опыте ВКЗ и представленном в подразд. 2.4.3, и начинают с формирования ГС случайного признака, выборки из нее, образования вариационных рядов с их свойствами и проверки гипотезы предполагаемого нормального закона отклонения случайного признака от МО исследуемых рядов [66, 67].

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

Объем N генеральной совокупности случайного признака согласно (20) выборками объемом n = 2 с учётом количества элементов K = 993 в соответствии с (21), (22) равен 492 528 и конструируется комбинаторным методом по формуле (23) путем перебора возможных комбинаций каждого элемента со всеми оставшимися элементами в исследуемом диапазоне tн′ − tв* .

По выражению (24) вычислены МО распределения ГС слу-

чайного признака τ′ = 4,854533 с (первый центральный момент)

0

и дисперсия распределения σ02 = 7469 291,26 с2 (второй цен-

тральный момент).

В результате получаем ПВ и нереальную дисперсию, величина которой отражает сильное отклонение от экспоненциального закона случайного признака в данной ГС, что свидетельствует о влиянии различных случайных факторов на опытные данные убывающего ПП данной СМ при её стендовых испытаниях. Для выяснения причин такого отклонения следует перейти от ГС случайного признака к его вариационному ряду. Ввиду огромного размера вариационного ряда случайного признака он представлен фрагментарно в виде его начальной части, середины и конца (табл. 3, прилож. 3). Анализ полученного ряда позволил в опыте ГП при разомкнутой ОЯ без гасительного сопротивления выявить 792 отрицательных и 3 нереальных значения случайного

146

признака, которые образованы по формуле (20). Размах вариационного ряда (27) даёт огромную величину, составляющую несколько сот тысяч (табл. 23). Всё это, вместе взятое, требует предварительного приведения вариационного ряда к удобному для дальнейшего использования виду. Отброс нереальных значений случайного признака по концам вариационного ряда позволил снизить его объём до 490 363 значений случайного признака

 

с МО

τ

 

= 8,889402, дисперсией σ2

= 0,319670 и размахом, пре-

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вышающим

ожидаемое

МО

примерно

в

два

раза (табл.

24),

 

(табл. 4, прилож. 3).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 23

 

 

 

 

 

 

 

 

 

 

Полная ГС без нереальных значений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

, с

 

2

 

2

 

(

)

 

 

,

 

(

) ,

 

 

 

 

 

 

 

ОЯ

Rг

 

 

 

 

N по

 

σ0

, с

 

 

 

τ

kj

 

макс

 

 

τ

kj

 

мин

 

 

 

R, с

 

 

υ, %

 

K

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(21–23)

по (24)

 

по

(24)

 

 

 

 

с

 

 

 

 

 

с

 

 

 

по (27)

 

по (28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ХХ

 

нет

993

492525

4,854533

7469291,26

17576,91

–1914573,731932150,64

56297,91

 

 

есть

400

79799

3,425606

 

27,37

 

 

524,91

 

 

–961,21

1486,13

 

152,74

 

КЗ

 

нет

286

40754

0,760615

 

12,72

 

 

442,40

 

 

–232,31

 

 

674,70

 

468,89

 

 

есть

288

20989

–0,078222

 

184,24

 

1709,08

 

–484,51

2193,59

 

–17352,58

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 24

 

 

 

 

 

 

 

 

 

 

Нормализованный вариационный ряд ГС

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N по

τ, с

 

2

, с

2

 

 

 

)

макс

,

 

)

мин

,

R, с

 

υ, %

 

ОЯ

 

Rг

 

 

 

K

0

 

 

σ0

 

 

 

 

 

kj

 

 

 

kj

 

 

 

 

 

 

 

(21–23)

по (24)

по

(24)

 

 

 

 

 

с

 

 

 

с

 

 

 

по (27)

по (28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ХХ

нет

 

993

490363

8,889402

0,319670

 

 

16,900178

0,869737

16,030441

6,36

 

 

 

 

есть

 

400

79416

3,416943

0,070953

 

 

6,083109

 

0,762030

5,321079

7,80

 

КЗ

 

нет

 

286

40425

0,733899

0,007666

 

 

1,394839

 

0,070496

1,324343

11,93

 

 

 

 

есть

 

288

13382

0,284441

0,004200

 

 

0,564459

 

0,004487

0,559972

22,78

С использованием программы (прилож. 7, программа № 3) на рис. 44 представлена гистограмма распределения нормализованной ГС случайного признака с целью визуального обнаружения предполагаемого нормального закона с количеством интервалов

k= 11. На рис. 45, 46, 47 представлены аналогичные гистограммы

сколичеством интервалов k = 11 при других состояниях ОЯ.

147

Рис. 44. Гистограмма при разомкнутой ОЯ без гасительного сопротивления

Рис. 45. Гистограмма при разомкнутой ОЯ с гасительным сопротивлением

Рис. 46. Гистограмма при замкнутой ОЯ без гасительного сопротивления

148

Рис. 47. Гистограмма при замкнутой ОЯ с гасительным сопротивлением

Гистограммы на рис. 44–47 близки к нормальному распределению случайного признака, определенного по формуле (20).

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

Гипотеза предполагаемого нормального закона распределения базового случайного признака подтверждается расчётами широко используемого в практике ТВ и МС критерия согласия Пирсона, называемого иначе «хи-квадрат». В табл. 3, прилож. 4 по данным полученного нормализованного вариационного ряда ГС случайного признака представлены результаты расчёта критерия согласия Пирсона в соответствии с рассмотренной методикой в подразд. 2.4.3.3. С целью сокращения времени вычислений количество интервалов k было снижено до пяти, а для снижения объёма ГС случайного признака шаг элементов ПП был увеличен с 0,01 до 0,1 с. Расчёт критерия χ2 произведён для объёма ГС вариационного ряда случайного признака, равного N = 4931, МО, равного 8,895642 с и со среднеквадратичным значением дисперсии, равным σ = 0,420744 с. С учётом количества интервалов разбиения, параметров нормального распределения число степеней свободы оказалось равным двум. При уровне значимости, принимаемым обычно равным 0,05, рассчитанное значение критерия согласия Пирсона χ2 = 0,4955388<6 оказалось меньше таблично заданно-

149

го [128], следовательно, гипотеза подтверждается. По результатам расчетов (табл. 3, прилож. 4) на рис. 48 приведена гистограмма распределения случайного признака для вариационного ряда ГС.

3200

3000

2800

2600

2400

2200

2000

1800

1600

1400

1200

1000

800

600

400

200

0

-∞– 7,5200

7,5200 – 8,7284

8,7284 – 9,0449

9,0449 – 9,9463

9,9463 – +∞

Рис. 48. Гистограмма, построенная по результатам расчетов критерия согласия Пирсона

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

3.8.3.3. Оценка уровня зашумлённости и степени отклонения случайного признака

Для оценки уровня зашумленности опытных данных в исследуемом диапазоне переходной составляющей формируется вариационный ряд по выборке из генерального ряда с определением размаха ряда по их концам и коэффициента вариации (табл. 3.5, прилож. 3). Объем выборки случайного признака на исследуемом участке tн′ − tв* , равном nв = 116 250, получен в опыте ГП с разомкнутой ОЯ без гасительного сопротивления аналогично ГС по выражению (19) и K = 125 (из 993) элементам от нижней границы tнперебором по два с оставшимися элементами до верхней

150

Соседние файлы в папке книги