Диссертация (1024744), страница 38
Текст из файла (страница 38)
П.1.График исходной функции результатов измерений и наилучшей моделипрогнозированияВЫВОД:Идентификация функции преобразования при минимальном СМПНсоответствует ММКМНК. Код структуры модели ν: 1110010. Функцияпрогнозирования принимает вид:__S ( x ) = 1854,037 + 12,6368 ⋅ x − 0,2256 ⋅ x 2 + 8,4445 ⋅ x 5С использованием программы «ММИ-поверка 2.0» был определеннаилучший вид распределения (Гаусса) и была получена свертка возможныхзначений основной погрешности СИ. Верхняя граница свертки: 0,1205 нс.Нижняя граница свертки: -0,25 нс.292П.4. ОЦЕНКА ЭКВИВАЛЕНТНОСТИ ИЗМЕРЕНИЙ ГАУССОВСКОЙМОДЕЛИПусть условиями формирования оптимального плана измерения оценкиэквивалентностиматематическоймоделифункциираспределенияслучайного результата будут условия: x – const – измеряемая величина; Y(x) –однократный случайный результат измерения с функцией распределенияF0(y); F(y; θ ) – математическая модель функции распределения;θ = (θ1...θr )–вектор параметров; [a, b] – интервал возможных значений случайнойвеличины Y(x) [94, 97, 101].Очевидно, что как бы ни была сложна математическая модель функциираспределения F(y; θ ), она не может быть тождественна истинной функциираспределения F0(y), т.
е. равенство F(y; θ ) = F0(y) является невозможнымсобытием и потому не может быть принято в качестве нулевой гипотезы H0.По этой причине сформируем эту гипотезу на малых отклоненияхматематической модели от истинной функции распределения следующимобразом.Строим оптимальный план измерений.1) Определяем количество интервалов разбиения:n = 12 .Тр=0,1; α0, β0=0,1, ξ0, ξ1 = 0,2, r=2.На основе этих данных получим:p(n) = 3 n = 2.289ν = n − (r + 1) = 9p=1= 0.083nПриведенный допуск равен Tε p =ε*p =ρ( n )2,15Tεp =1, 0 = 1, 08,22Tp= 1.2p293() = 1.099()ε p0 := ε p ⋅ 1 − ξ 0ε p1 := ε p⋅ 1 + ξ1=1.648ai ( ⋅) =ν + 2 N *ε 2pi= ai ( ν, N * , ε pi ) , i = 0,1,* 2ν + N ε piν*i ( ⋅) = ai−1 ( ν + N *ε 2pi ) = ν*i ( ν, N * , ε pi ) , i = 0,1.a = 1.915 , ν = 55.15600a = 1.96 , ν = 115.48310Согласноa1 ( ν, N * , ε p1 )[231, 232]a0 ( ν, N , ε p 0 )*⋅q p ,ν1* (⋅),β0q p ,ν*0 (⋅),1−α0= λ ( N * ; ν, ε p 0 , ε p1 , α0 , β0 ) = 1. ,тогдаq = 82.4820q = 83.0911aqОкончательно λ = 1 1 = 1.031 .a q0 0Поскольку λ = 1.031 > 1 , то оптимальный план оценки эквивалентностиГауссовскоймоделифункциираспределенияслучайногорезультатаизмерения имеет следующую структуру( n , Nˆ , uˆ ) = (12; 920; 160.408).102) Реализация многократных измерений и определение попаданий взаданныеинтервалынапримереизмеренийфазохронометрическойинформации на турбоагрегате.№1 2 3 4 5 6 78910 11 12интервалаКоличество22 28 47 47 71 99 123 125 112 88 73 52попаданий2943) Определение экспериментального значения аргумента решающейфункции и принятие решенияˆ − νˆ )( Np,2nqp = ∑k =1okˆNpгде νˆ ok , k = 1, n – количество результатов измерений из совокупностиyj(x), j = 1, Nˆ , попавших в интервал Iyk.(( )(( )))0, если q ≤ uˆ − принимается гипотеза H F y;θˆ ≈ F ( y ) ,p000r (q p ) = ˆ1, если q p > uˆ0 − принимается гипотеза H 1 F y;θ ≈/ F0 ( y ) .q p = 27.156 ≤ 160.408 , следовательно принимается гипотеза H0-модель измерений эквивалентна Гауссовской модели функции распределенияслучайного результата измерения.
Может быть реализована методикаповеркиизмерительныхэксплуатации.фазохронометрическихсистемвпроцессе295П.3. РАСЧЕТ ХАРАКТЕРИСТИК РАБОТЫ СИСТЕМЫРЕГУЛИРОВАНИЯ ПРИ РАЗЛИЧНЫХ ЗНАЧЕНИЯХ ПОСТОЯННЫХВРЕМЕНИПри работе канала регулирования по напряжению имеют местоследующие характеристики системы.Данныеотключенномпосравнениюконтурережимоврегулированияработыприпроизводнойвключенномтокаигенератораприведены в Таблице 21.Таблица 21.Максимальные значения параметров рабочего режима турбоагрегата приотсутствии влияния внешней сети и системы возбужденияРежимОтклВклЗна Измене Измене Измен Изме Измен Изм Измене Изменениениеение нение ение ененниениечение амплит крутятокаие амплит амплит ампли токащегоудыпо оси токапотудыудыудысигмоменуглаугла оси d, q, А возбуглаугланалт накач.удитАа, В качания качания кач.еля, муфты муфтеРНД, ротораротораРГ-РВ, возбуАвозбудугл.генератдителугл.ителя,ора,минминя,угл.угл.кН٠мминмин44220,03,4240,03,5250,04,60,380,171,000,181,551,68±1,4±1,2±70,0±60,0∆t 2 =0,0005 сПри возникновении аварийных режимов, коротких замыканий вгенераторе, шинопроводе или трансформаторе, после внезапного отключениягенератора необходимо быстро уменьшить магнитное поле обмоткивозбуждения генератора.
Эта операция носит название гашение поля иосуществляется специальным автоматом гашения поля (АГП). К устройству296АГП предъявляются два основных, иногда противоречащих друг другу,требования: время гашения поля должно быть возможно меньшим, авозникающее при гашении индуктированное перенапряжение в обмоткеротора не должно превосходить допустимых значений.При различных характеристиках данная система будет оказыватьразличное влияние на работу турбоагрегата и раскачку генератора.В Таблице 22 приведены результаты расчетов для коэффициентов:- постоянная времени канала производной тока ротора Ta_1If=0,3;- постоянная времени канала производной тока ротора Ta1If=0,06;Таблица 22.Максимальные значения параметров рабочего режима турбоагрегата приотсутствии влияния внешней сети и системы возбуждения№Измене ИзмененШагЗначиениеение интегрисигн рования амплиту амплитуды угла ды угла,сала,качания качанияВРНД,роторагенерат угл.
минора,угл.мин123±7,5±7±70,000500,000100,000014,412,411,44.612.322.5Измене Измене Изменениениениеамплит тока по токаоси d, по осиудыq, ААуглакач.роторавозбудителя,угл.мин3,812,114,50,860,500,700,160,255,20Времядемпфирования, с0,60,740,74В Таблице 23 приведены результаты расчетов для коэффициентов:- постоянная времени канала производной тока ротора Ta_1If=0,135;- постоянная времени канала производной тока ротора Ta1If=0,02;297Таблица 23.Максимальные значения параметров рабочего режима турбоагрегата приотсутствии влияния внешней сети и системы возбужденияПоствремениканала№ пр.токаротораTa_1IfTa1IfИзменениеЗначамплитуениеШагды углаамплинтегри качаниясигнрования ротораала,,сгенератВора,угл.мин0,14±70,000050,020,14±70,00010,020,10±70,000010,150,05±100 0,000125,000,06±100 0,00010,300,01±100 0,00010,30123456СогласнопеременныхтехническоймоделиконтураИзменениеИзменен амплитудыиеуглаамплитукач.ды углакачания роторавозбудРНД,угл.
мин ителя,угл.минВремядемпфирования, сМаксимальныйразброспериодавращения ТА,мкс8,239,28,40,30288,89.18,90,224311,422.514,50,742938,191,240,50,265713,128,516,80,09893,68,36,50,0243документациизначениярегулированияпараметровнапряжениястатораиисистемного стабилизатора ограничены диапазоном:- постоянная времени канала производной тока ротора Ta_1If – от 0,0до 0,3;- постоянная времени канала производной тока ротора Ta1If – от 0,01до 0,06.298При выключенном канале регулирования по производной тока ротораразброс значений периода вращения (хронограмма вращения ТА) можетсоставлять от 10 до 100 мкс.Рис.
П.2.Хронограмма вращения – результаты эксперимента (разброс 33,6 мкс)299П.4. ОЦЕНКА И ИЗМЕРЕНИЕ ДЕЙСТВУЮЩИХ НА УЗЕЛНАГРУЗОК В СООТВЕТСТВИИ С ТЕХНИЧЕСКИМ ЗАДАНИЕМАнализ статистики аварий на тепловых электрических станциях выявилопределенные повторяющиеся дефекты при эксплуатации турбоагрегатов.Наиболее частым поломкам подвержены шпильки (болты) муфт междувозбудителями и генераторомнератором и валы возбудителей генераторов.На Рис.
П.3. представлены графики измененияэлектрическиххарактеристик работы турбоагрегата ТВВ-200-2ТВВК-200-130130 при двукратнойфорсировке, If=2,0.Рис. П.3.Измерение статического угла скручивания ротора возбудителя и значениеактивной мощности на ВГ при двукратной форсировке300Рис. П.4.Графики изменения электрических характеристик работы турбоагрегатаТВВ-200-2- К-200-130130 при полуторной форсировке, If=1,5Рис.
П.5.Измерение статического угла скручивания ротора возбудителя и значениеактивной мощности на ВГ при полуторной форсировкеВращающий момент может быть определен через выражение.••J 5 ϕ1 + qв (ϕ 5 − ϕ1 ) = M КР ,((П.1)301где J 5 - момент инерции ротора возбудителя;ϕ1 , ϕ 5 - углы поворота дисков на торце возбудителя и на муфте РВ-РГпри скрученном вале;qв - крутильная жесткость ротора возбудителя;M КР - внешний скручивающий момент.Для определения значения крутящего момента можно принять, чтомеждуизмерительнымиметкамидискачастотавращенияявляетсяпостоянной, тогда••J 5 ϕ1 ≈ 0 .