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

Прогнозирование технического состояния топливорегулирующей аппаратуры газотурбинных двигателей в процессе длительных испытаний Текст научной статьи по специальности «Математика»

CC BY
236
37
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ТОПЛИВОРЕГУЛИРУЮЩАЯ АППАРАТУРА / ГАЗОТУРБИННЫЙ ДВИГАТЕЛЬ / ПОЛУНАТУРНЫЙ ИСПЫТАТЕЛЬНЫЙ МОДЕЛИРУЮЩИЙ СТЕНД / ПРОГНОЗИРОВАНИЕ ТЕХНИЧЕСКОГО СОСТОЯНИЯ / UNIVERSAL CONTROL FUEL PUMP / GAS TURBINE ENGINE / HALF-NATURAL MODEL TEST STAND / TECHNOLOGIC FORECASTING

Аннотация научной статьи по математике, автор научной работы — Шендалева Елена Владимировна

В статье рассматривается проблема прогнозирования технического состояния топливорегулирующей аппаратуры газотурбинных двигателей. Основное содержание исследования составляет анализ применимости метода экстремального прогнозирования для обеспечения гарантированного прогноза технического состояния топливорегулирующей аппаратуры. Метод экстремального прогнозирования позволяет выполнить прогноз на основании результатов измерений, полученных в ходе длительных испытаний на полунатурном моделирующем стенде и последующих плановых испытаний. Прогноз выполняется без анализа статистического распределения результатов измерений. Результаты исследования могут быть использованы для прогнозирования состояния технических систем, состоящих из элементов различной физической природы.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Extreme forecasting for universal control fuel pump and gas turbine engine in service obtaining for a long time trials

The article has considered the problem of technical state forecasting for universal control fuel pump and gas turbine engine in service. The subject matter of this article is the analysis of extreme forecasting method application for providing the assured forecasting of fuel pump technical state. This method makes it possible to forecast technical state with using of measurements results, obtaining for a long time trials of universal control fuel pump on half-natural model test stand and during subsequent planned tests. The forecasting has been performed without statistic distribution analysis of measurements results. The results of investigation may be used for forecasting of different technical system states.

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

sistema [Combined vibration protection system]. Patent RF, no 2015123195/05, 2016.

Бурьян Юрий Андреевич (Омск, Россия) -доктор технических наук, профессор, заведующий кафедрой Основы теории механики и автоматического управления ФГБОУ ВО ОмГТУ (644050, г. Омск, пр. Мира, 11, e-mail: bu-rian7@mail.ru).

Сорокин Владимир Николаевич (Омск, Россия) - доктор технических наук, доцент, профессор кафедры Основы теории механики и автоматического управления ФГБОУ ВО ОмГТУ (644050, г. Омск, пр. Мира, 11, e-mail: sorokin.vn@mail.ru).

Захаренков Николай Владиленович (Омск, Россия) - кандидат технических наук, доцент кафедры Машиноведение ФГБОУ ВО ОмГТУ (644050, г. Омск, пр. Мира, 11, e-mail: znick-bar@mail.ru).

Зелов Александр Федорович (Омск, Россия) -научный сотрудник лаборатории 120 научно-производственное объединение «Прогресс»

(644018, г. Омск, ул. 5-я Кордная, 4, e-mail: aleks.zelov@gmail.com).

Yuriy A. Buriyan (Omsk, Russian Federation) -Doctor of Technical Sciences, Professor, Fundamentals of the Theory of Mechanics and Automatic Control department, Omsk State Technical University (644050, Mira, 11 prospect, Omsk, Russian Federation, e-mail: burian7@mail.ru).

Vladimir N. Sorokin (Omsk, Russian Federation) -Doctor of Technical Sciences, Ass. Professor, Fundamentals of the Theory of Mechanics and Automatic Control department, Omsk State Technical University (644050, Mira, 11 prospect, Omsk, Russian Federation, e-mail sorokin.vn@mail.ru).

Nikolay V. Zakharenkov (Omsk, Russian Federation) - Ph. D in Technical Sciences, Science of Machines department, Omsk State Technical University (644050, Mira, 11 prospect, Omsk, Russian Federation, e-mail znickbar@mail.ru).

Aleksandr F. Zelov (Omsk, Russian Federation) -research engineer, NPP "Progress" (644018, 5th Kordnaya, 4, Omsk, Russian Federation, e-mail aleks.zelov@gmail. com).

УДК 621.45.018.2

ПРОГНОЗИРОВАНИЕ ТЕХНИЧЕСКОГО СОСТОЯНИЯ ТОПЛИВОРЕГУЛИРУЮЩЕЙ АППАРАТУРЫ ГАЗОТУРБИННЫХ ДВИГАТЕЛЕЙ В ПРОЦЕССЕ ДЛИТЕЛЬНЫХ ИСПЫТАНИЙ

Е.В. Шендалева

Омский государственный технический университет, Россия, г. Омск

Аннотация. В статье рассматривается проблема прогнозирования технического состояния топливорегулирующей аппаратуры газотурбинных двигателей. Основное содержание исследования составляет анализ применимости метода экстремального прогнозирования для обеспечения гарантированного прогноза технического состояния топливорегулирующей аппаратуры. Метод экстремального прогнозирования позволяет выполнить прогноз на основании результатов измерений, полученных в ходе длительных испытаний на полунатурном моделирующем стенде и последующих плановых испытаний. Прогноз выполняется без анализа статистического распределения результатов измерений. Результаты исследования могут быть использованы для прогнозирования состояния технических систем, состоящих из элементов различной физической природы.

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

Введение

Газотурбинные технологии находят всё более широкое применение в различных отраслях промышленности - прежде всего в авиации, судостроении, на железнодорожном

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

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

Испытания топливорегулирующей аппаратуры

Критерием достижения заданных технических характеристик согласованной работы газотурбинного двигателя (ГТД) и его топливорегулирующей аппаратуры (ТРА) могут служить результаты длительных испытаний ТРА, например, всережимных топливных регуляторов или дозаторов топлива, на полунатурном испытательном стенде. Данные испытания позволяют контролировать статические и динамические характеристики совместной работы ГТД и ТРА, тем самым, обеспечивая прогнозируемую безопасность эксплуатации ГТД и воздушных судов.

На рис. 1 представлена структурная схема полунатурного испытательного стенда для испытания всережимного топливного регулятора или дозатора топлива - гидромеханической части системы автоматического управления (ГМЧ САУ) [1]. Стенд содержит модель ГТД по контурам регулирования частоты вращения и давления за компрессором, реализованную программно, и преобразователи электрических сигналов с выхода модели ГТД в физические параметры, передаваемые в топливный регулятор. В соответствии с заданным режимом топливный регулятор дозирует топливо, сливаемое затем в закольцованную топливную систему. Сигнал измеренного расхода топлива передаётся в модель ГТД.

При полунатурных стендовых испытаниях топливных регуляторов (дозаторов топлива) совместно моделью ГТД в контуре обратной связи возникают вопросы его оптимального регулирования (отладки) в соответствии с заданными техническими условиями. Отладка топливного регулятора (дозатора топлива) осуществляется путём изменения положения регулировочных элементов (регулировочных винтов), установленных на его внешней поверхности, а также с помощью изменения параметрических характеристик электронной или электронно-гидравлической САУ ГТД. Фиксация положения регулировочных элементов выполняется контрящими элемента-

ми. Отладка топливного регулятора выполняется как вручную, так и в автоматическом режиме с использованием механизмов перемещения регулировочных и контрящих элементов. При этом в качестве последовательности выполняемых операций могут быть использованы типовые операции отладки, отладка на основе регрессионной модели, эталонной либо оптимизированной модели всережимного топливного регулятора (дозатора топлива) [2].

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

Одним из наиболее действенных способов обеспечения безаварийной эксплуатации ТРА ГТД является применение прогнозирующего контроля его надёжности и состояния [3]. По результатам длительных испытаний предсказание технического состояния ТРА обеспечивает её последующее рациональное использование и даёт возможность перейти к обслуживанию ТРА по состоянию. В основу прогнозирования состояния положен тот факт, что сами изменения состояния ТРА вызываются плавными изменениями её параметров, которые могут быть зарегистрированы средствами измерения. Основными же причинами дрейфа параметров ТРА, как показывает опыт эксплуатации, контроля и ремонта, являются процессы старения и износа. Именно эти процессы, идущие на молекулярном уровне, могут приводить к ситуациям, когда значение того или иного параметра выйдет за пределы допуска. Наблюдения за процессами дрейфа параметров в ходе длительных испытаний и последующая обработка результатов этих наблюдений при решении задачи прогноза обеспечивают сохранение работоспособности ТРА, позволяют своевременно принимать меры по предупреждению её повреждений и отказов во время эксплуатации.

Рис. 1. Полунатурный стенд для испытания ТРА ГТД I - полунатурная модель ГТД; II - гидромеханическая часть САУ; III - исполнительные механизмы перемещения регулировочных и контрящих элементов; IV - система управления исполнительными механизмами РЭ и КЭ - регулировочный и контрящий элемент, М - двигатели постоянного тока, Ф и фк - углы поворота регулировочных и контрящих элементов, ип, ир - сигналы управления электроприводом и электропневмопреобразователем, йт - расход топлива

Прогнозирование состояния ТРА. Для

практической реализации прогнозирующего контроля необходимо оценивать техническое состояние ТРА, иначе говоря, использовать модель ТРА в некотором пространстве функционирования. Поскольку ТРА достаточно полно характеризуется определённой совокупностью свойств (параметров), то набор

(вектор) этих параметров у = {у^ }"=о , гДе п -

число параметров, описывающих свойства ТРА, можно использовать для описания технического состояния. Указанный набор состоит из выходных и определяющих координат ТРА, при этом у еУ, где У - множество, в котором у принимает свои значения.

В течение длительных испытаний ? е Т ( Т- период эксплуатации) и последующей эксплуатации состояние ТРА под воздействием дестабилизирующих факторов меняется

случайным образом. Отсюда следует, что поведение ТРА на множестве Т описывается случайным процессом у , ? е Т , у еТ , где Т - множество элементарных событий. На вероятностном пространстве (Т, Г, Р), где Г - алгебра подмножества множества Т, Р -вероятностная мера, случайную векторную

функцию у(;)={у[(?)}П=о можно принять за искомую модель ТРА. Следует отметить, что длина набора {yi(¿УЩ=о может быть единичной, когда состояние ТРА характеризуется только одним параметром у, например, заданным расходом топлива вт. Вопрос о выборе совокупности координат, с достаточной полнотой отражающей свойства ТРА, решается в каждом случае индивидуально, исходя из конструкции, эксплуатационных характе-

ристик, назначения, специфики работы и других показателей ТРА.

В рамках принятой модели ТРА задача прогнозирования дрейфа параметров является задачей прогнозирования случайного процесса у(?). Суть указанной задачи состоит в предсказании траектории движения у(?) на множестве Т. Информационной основой для её решения служат некоторые априорные сведения, а также данные контрольных измерений у(?) при длительных испытаниях ТРА. Совокупность априорных сведений чаще всего ограничивается описанием некоторых структур, пригодных в качестве модели у(?). Параметрические характеристики таких структур редко бывают известны. Для объектов, относящихся к ТРА, наиболее типичным является представление у(?) как

у()_ А -Ф(г),

(1)

где А = ||а|| - матрица случайных коэффициентов размером (п+1)х(т+1), / = 0, ..., п,

} = 0, ..., т, т - размер вектора (?})т_0 случайных коэффициентов непрерывно дифференцируемых детерминированных функций

времени ф(?)_{ф^(^)}т_0 ■ Описание у(?) в

виде (1) можно рассматривать как разложение случайного процесса по детерминированному базису. Справедливость такого представления следует из физической природы процессов старения, износа, происходящих в ТРА, и подтверждается данными имеющейся статистики. Информационная полнота этих данных зависит от многих факторов и, в частности, от таких характеристик как достоверность и форма осуществления контроля.

Влияние окружающей среды, несовершенство средств измерения, недостаточная квалификация обслуживающего персонала приводят к тому, что значение результата измерения отличается от истинного на некоторую случайную величину - ошибку измерения, то есть г (?)_ у(?)+е(?), где г (?) -результат измерения у(?). Мерой точности (прецизионности) служат статистические характеристики ), при этом обычно полагают, что ошибки измерения независимы и распределены по нормальному закону [4]. В слу-

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

р(е(?) < А(?))> Р0, ? е Т. Предельная ошибка А(?) служит заданным показателем точности.

По форме выполнения различают непрерывные и дискретные по времени измерения. На практике возможности проведения непрерывных измерений у(?) ограничены вследствие большой трудоёмкости и дороговизны. Чаще встречаются дискретные по времени измерения как более простые и дешёвые. Результатом таких измерений является последовательность г (?), ? е Ти е Т , где Ти -

период испытаний.

При этом информационная полнота данных измерений ограничена следующим:

- измерение производится со случайной ошибкой ), для оценки влияния которой

м ожно использовать дисперсию )], либо предельную ошибку измерений;

- возможно только дискретное по времени измерение, в результате которого определяется последовательность г(?), ? е Ти е Т .

Установленный уровень информационной полноты характеризует условия, из которых целесообразно исходить при решении задачи прогнозирования у(?).

Алгоритм решения задачи прогнозирования изменения параметров ТРА можно построить на основе существующих статистических методов (метода наименьших квадратов, максимального правдоподобия и т.п.). Однако для практического применения статистических методов необходимо задаться некоторыми вероятностными характеристиками ошибок измерений и модели случайного процесса у(?). В действительности фактические значения этих характеристик не совпадают с прогнозируемыми, что может привести к снижению точности получаемых результатов по сравнению с теоретическими оценками. На точность решения рассматриваемой задачи влияют и погрешности модельных зависимостей [5]. В результате измеренные данные содержат и ошибки модели. Эти ошибки нельзя рассматривать как случайные и поэтому с ними нельзя оперировать как со случайными величинами. Решение задачи про-

гноза в таких условиях с помощью статистических методов может привести к неоправданно оптимистическим оценкам у(г). Если фактическое состояние ТРА окажется хуже предсказанного, то конечная цель прогнозирования - обеспечение эксплуатационной работоспособности ТРА не может быть полностью достигнута [6]. Значительно меньшую опасность представляет получение пессимистических (гарантированных) оценок у(г). Использование результатов гарантированного прогноза при определении рационального режима эксплуатации ТРА обеспечивает своевременное принятие мер по предупреждению её повреждений и отказов в условиях ограниченности исходных данных.

Применительно к модели у(г) вида (1) в этом случае известен только детерминированный базис {у у(г)]\_о , а для ошибок измерения е(г) задана лишь область их возможных значений Е (значения предельных ошибок Аг), то есть г( е Е, г е Ти е Т . Для решения задачи прогнозирования может быть применён метод экстремального (гарантированного) оценивания [7], пригодный для использования в условиях ограничения исходных данных. Применение данного метода наиболее полно обеспечивает использование исходной информации при гарантированной достоверности и точности результатов прогноза. При разработке прогноза отсутствует необходимость в использовании стохастических свойств ТРА, определяемых как собственными параметрами, так и параметрами ГТД и подключённых к ТРА агрегатов.

Суть экстремального оценивания состоит в определении гарантированных пределов изменения оцениваемой величины у(г) при

г е Ти е Т , то есть на получение интервальной оценки у(г). Интервальный прогноз у(г) также может быть выполнен и с помощью статистических методов, однако его достоверность невозможно гарантировать без рассмотрения гипотез о стохастических свойствах у(г) и ошибок измерений в(г) [8].

Рассмотрим теоретические предпосылки использования метода экстремального оценивания в ситуации, когда техническое состояние ТРА характеризуется одним параметром у(^, например, расходом топлива вт. Тогда область возможных значений ошибки измерения можно представить как [Ди, Д^], где Ди, Д2Г - предельные границы для при

т

г е Ти е Т , а у() _ 2 ау (г )фу (г), где

у _0

{ау (г )]т_0 - случайные коэффициенты;

{

фj(tц. о - непрерывные детерминированные функции. По результатам измерений случайного процесса y(t) на интервале времени Ти с Т можно записать систему неравенств

z,-au < ^ < z, +A2i, t е Ти с Т , (2)

где zt - результат измерения y(t) в момент t е Ти. С учетом заданной модели y(t) неравенства (2) принимают вид

z - Alt < A • ф < Z + A2t,

(3)

где ? _ {г(} А1г _ {А1г}, А2г _ {А2г}, Г = ^ ^ ..., к А = {а](г1 )}у_0 ; ф _ ||фу (ц)||, I = 0, ..., к, г1 е Ти е Т, ф - матрица значений функций {фу(гг)}т_о в точках измерения y(t), к -

общее число измерений функции {фу(^)]т_0

в моменты времени tl.

Полученная система неравенств (3) определяет область возможных значений коэффициентов {ау(гг)}т_0, построенную по наблюдениям г, в пространстве Рт+1. Данная область образована пересечением двух полупространств размерностью (к + 1), каждое из которых описано частью неравенства (3),

L = LiH L2,

Lx = {а е Rw+!: z-Дх < А • ф} Lx = Z^/,

1=0

Lll = {А е R-+1: ( А,ф1)> z(tl)-A1(f/)} фг = {фj(tl)j=0;

L2 = {А е Rm+1: z + A2 > А•ф} L2 = [>21,

l=0

L2l = {А е R-+1: (а,ф1 )< z(tl) +A2(tl)} ф1 = {фj(tl)У

Область I. представляет собой выпуклый многогранник. В пространстве Рт+1 область I выделяет все наборы коэффициентов

{а у(?)}И_о, с которыми реализации у(4) могут быть измерены в моменты времени 4,. Тогда в

области L с R"

могут быть найдены гаран-

z

тированные относительно наблюдений пределы изменения {a.j(t)jm_0. Поиск указанных пределов можно осуществить, решая 2(m+1) задач линейного программирования

1) aj = max, j = 0, ..., m,

2) aj = min, j = 0, ..., m

при ограничениях области L, то есть координаты всех точек этой области должны удовлетворять неравенствам |aj| < d, где d - некоторое положительное число; k > m; |zt| < да, t/ е Ти, где (k+1) - число измерений в моменты времени t.

Данный способ приемлем для нахождения векторов P = {pj jm=0 и Q = {qj }m=0 , где

Pj = minL , qj = maxL , j = 0, ..., m, если P

aj aj

и Q используют для определения гарантированных пределов изменения y(t) или определения всех p(t) < y(t) < q(t), t е Tu с T . Каждой точке области L соответствует определенный вектор коэффициентов

А = {a j 1

=0'

при

этом функции

Р]Ф](?) и У]Ф] (г) можно Рас"

у_0 ]_0

сматривать как пределы изменений у(4), ? е Ти е Т . При этом выполняется аппроксимация выпуклого многогранника I в виде гиперпараллелепипедов Ц и 12 (рис. 2).

Рис. 2. Аппроксимация выпуклого многогранника L параллелепипедом Li

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Представление ограничений у($,

т

г е Т > Ти в виде р(груфу(г) и

у=1

т

д(г)= Е дуфу (г) содержит не только реали-

у=1

зации у(^), для которых справедливо неравенство (2), но и такие, для которых оно выполняется частично. Высокая трудоемкость и сложность вычисления р(г) и д(г) ограничивает возможность их использования для решения задачи прогноза у(4). Более точно пределы изменения у(4), г е Т > Ти можно описать согласно свойствам детерминированного базиса {фу(г)}и=о модели (1), путем построения экстремальных полиномов Карлина [7]. Эти полиномы обладает полезным с точки зрения гарантированного прогноза свойством, а именно: на интервале времени Ти е Т они удовлетворяют неравенствам (2), а на интервале времени Т > Ти являются предельными для тех реализаций у(4), которые также удовлетворяют всем неравенствам (2).

Класс функций {фу(г)}^, используемых

на практике для аппроксимации случайных процессов эксплуатационных изменений параметров ТРА, достаточно ограничен. Обычно {фу(г)]т_о - непрерывно дифференцируемые функции, образующие систему Че-бышева [9]. Реальные процессы дрейфа параметров ТРА отличаются большой инерционностью, при стабильных условиях в них отсутствуют резкие изменения [3]. Среди функций рассматриваемого вида наиболее часто для описания у(4) используются степенные или экспоненциальные, т.е. наборы

типа

{ ' т=о - к' }т=о.

и \е 1 ] у=0, где а - показатель

распределения.

Принадлежность {фу(г)]т=0 к классу непрерывно дифференцируемых функций можно использовать для нахождения пределов

изменения у(4), г е Т > Ти . По результатам измерения у(4) при г е Ти е Т могут быть найдены экстремальные полиномы у(4)- и у(4)+, для которых гарантированно выполня-

ется у(4)- < у(4) < у(4)+ (2), г е Ти е Т . При

этом область [у(4)-; у($+] является наименьшей из допустимых, содержащих все возможные результаты наблюдений г реализаций у(4).

Определение у($- и у(4)+ возможно с помощью линейного программирования

1) Еа,- (г*) у (г * )= тах,

2) Еа, (г *) у (г * )= тт

(4)

при ограничениях (2), где 4 - любая фиксированная точка из Т > Ти (возможность использования произвольной точки из Т > Ти вытекает из экстремальности у(4)- и у(£)+ на

Т > Ти).

Применение экстремальных полиномов у(£)- и у(4)+ для прогнозирования у(4) по сравнению с использованием функций р(4) и обеспечивает следующие преимущества:

- трудоёмкость вычисления у(4)- и у(£)+ в т раз меньше (т - размер детерминированного базиса модели у($);

- область [у(4)_, у(4)+], г е Т > Ти содержит только возможные относительно результатов наблюдения г реализации у(4).

Алгоритм прогнозирования изменения параметров ТРА на основе решения (4) может быть представлен в следующем виде:

1) на интервале времени Ти е Т производится к контрольных измерений у(4), причём к > т;

2) данные измерения используются для построения у(£)- и у(4)+ путём решения задачи (4).

В ходе реализации алгоритма определяют пределы изменения у(£), г е Ти е Т (экстремальные полиномы у(£)- и у(4)+) с помощью линейного программирования (4).

Рассмотрим пример определения гарантированного технического состояния ТРА, найденного с помощью приведённого выше теоретического обобщения.

Результаты изменения выходного параметра ТРА (расхода топлива вТ [л/час]) при заданном угле поворота рычага управления

двигателем (аРУД) и частоте вращения ротора ГТД (п) можно аппроксимировать линейной зависимостью вт = вт0 + где вт0, вт1 -

случайные параметры. Результаты измерения вт(^ показали, что в моменты времени = 100 час. и = 200 час. бт(41) = 330 ± в, бт(42) = 331 ± в, где |в| < 3,2 - погрешность измерения. Требуется определить техническое состояние ТРА в момент времени = 6000 час. и после измерения расхода топлива в этот момент времени дать прогноз о поведении вт при = 12000 час.

Ввиду малости интервала А4 = 42 - 41, по которому осуществляют прогноз на длительный период времени, и для улучшения сходимости алгоритма прогнозирования в качестве аргумента функции фу (?) используем /д

4. Решив исходную систему уравнений 1) для положительного и отрицательного значения ошибки в получим 2), 3)

1)

[330 = Gt 0 + Gt1 • lg t1, I331 = Gt 0 + GT 1 • lgt2;

2)

3)

[326,8 = G' 0 + G' 1 • lgt1, [334,2 = G' 0 + G' 1 • lg 12;

[333,2 = G'T' 0 + Gj1 • lg t1, I 327,8 = G'0 + G„ • lgt

l2

и соответствующие экстремальные значения

От:

вт+ = 277,6 + 24,58-/д 4; вт- = 369,1 - 17,94 1д

Согласно данным уравнениям величина бт(43) будет находиться в диапазоне [301,3; 370,5], бт(44) - в диапазоне [295,9; 377,9].

Результат последующего измерения бт(43) составил 315 л/час., что соответствовало первоначальному прогнозу. По результатам измерения бт(43) составлен новый прогноз для расхода топлива бт(44):

[326,8 = GT 0 + GT1 • lgt1, [318,2 = GT 0 + GT 1 • lgt3;

[333,2 = Gf0 + Gj4 • lg t1, [311,8 = GT 0 + GTr lgt3

и получены соответствующие экстремальные значения GT: GT+ = 336,5 - 4,84-/gt; GT- = 357,3 - 12,04-/gt. По результатам данного прогноза величина GT(t4) будет находиться в диапазоне [308,2; 316,7]. Последующий результат измерения расхода топлива составил GT(t4) = 310 ± в, то есть результат измерения соответствовал предсказанному. Используя полученный результат можно продолжить процесс прогнозирования, для этого составляют новые системы уравнений:

[326,8 = GT 0 + GT1 • lgt1, [313,2 = GT 0 + GT1 • lgtA;

[333,2 = GT 0 + GT 1 • lgt1, [306,8 = GT 0 + GT1 • lgt4

и определяют соответствующие экстремальные значения GT: GT+ = 339,9 - 6,54-/gt; GT- = 358,6 - 12,70-/gt. На основании приведённых расчётов можно отметить, что использование результатов дополнительных измерений повышает точность процедуры прогнозирования. Результаты построения расчетов экстремальных полиномов приведены на рис. 3.

Приведенный метод гарантированного прогноза удовлетворяет следующим требованиям:

- выполнение условий несмещённости и сходимости прогноза;

- однозначность результата прогноза;

- минимальная сложность при наличии малого количества измерений.

Рис. 3. Результаты построения экстремальных полиномов

Выводы

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

2. Для обеспечения несмещённости и сходимости прогноза при реализации метода в качестве аргумента экстремальных функций предпочтительным является использование логарифмической функции времени.

3. Метод прогноза имеет оптимальную сложность при минимуме полученной измерительной информации.

4. Однопараметрическое гарантированное прогнозирование является эффективным инструментом обработки результатов измерения, полученных в ходе длительных испытаний ТРА ГТД и последующих плановых испытаний, для обеспечения надёжной эксплуатации системы «ТРА - ГТД».

Библиографический список

1. Шендалева, Е. В. Технология регулирования топливной аппаратуры систем автоматического управления газотурбинных двигателей с использованием моделирующих стендов / Е. В. Шендалева, В. В. Жильцов, В. Ю. Тэттэр // Сборка в машиностроении, приборостроении. - 2005. - № 7. -С. 15 - 21.

2. Жильцов, В. В. Новый метод настройки дозаторов топлива систем автоматического управления газотурбинных установок / В. В. Жильцов, Е. В.

Шендалева // Омский научный вестник. - 2001. -Вып. 17. - С. 95 - 97.

3. Дружинин, Г. В. Надёжность автоматизированных систем. - 3-е изд., перераб. и доп. / Г. В. Дружинин. - М. : Энергия, 1977. - 535 с.

4. Боровиков, В. П. Прогнозирование в системе STATISTICA® в среде Windows. Основы теории и интенсивная практика на компьютере : учеб. пособие / В. П. Боровиков, Г. И. Ивченко - М. : Финансы и статистика, 2006. - 384 с.

5. Методы оптимизации испытаний и моделирования систем управления газотурбинными двигателями / В. Т. Дедеш, В. М. Герман, В. Г. Авгу-стинович и др. ; под общ. ред. В. Т. Дедеша. - М. : Машиностроение, 1990. - 160 с.

6. Прогрессивные технологии моделирования, оптимизации и интеллектуальной автоматизации этапов жизненного цикла авиационных двигателей : монография / А. В. Богуслаев, Ал. А. Олей-ник, Ан. А. Олейник, Д. В. Павленко, С. А. Субботин ; под ред. Д. В. Павленко, С. А. Субботина. -Запорожье : ОАО «Мотор Сич», 2009. - 468 с.

7. Карлин, С. Математические методы в теории игр, программировании и экономике / С. Карлин ; пер. с англ. Н. А. Бодина, Л. И. Горькова и др. - М. : Мир, 1964. - 840 с.

8. Худсон, Д. Статистика для физиков: Лекции по теории вероятностей и элементарной статистике. - 2-е доп. изд. / Д. Худсон ; пер. с англ. В.Ф.Грушина ; под ред. Е. М. Лейкина. - М. : Мир, 1970. - 296 с.

9. Карлин, С. Чебышевские системы и их применение в анализе и статистике / С. Карлин, В. Стадден ; пер. с англ. под ред. С.М. Ермакова. -М. : Наука, 1976. - 568 с.

EXTREME FORECASTING FOR UNIVERSAL

CONTROL FUEL PUMP AND GAS TURBINE ENGINE IN SERVICE OBTAINING FOR A LONG TIME TRIALS

Elena V. Shendaleva

The summary: The article has considered the problem of technical state forecasting for universal control fuel pump and gas turbine engine in service. The subject matter of this article is the analysis of extreme forecasting method application for providing the assured forecasting of fuel pump technical state. This method makes it possible to forecast technical state with using of measurements results, obtaining for a long time trials of universal control fuel pump on half-natural model test stand and during subsequent planned tests. The forecasting has been performed without statistic distribution analysis of measurements results. The results of investigation may be used for forecasting of different technical system states.

Keywords: universal control fuel pump, gas turbine engine, half-natural model test stand, technologic forecasting.

References

1. Shendaleva E. V., Zhiltsov V. V., Tetter V. Y. The fuel apparatus adjustment technology for automatic control systems of gas turbine engines with model stands using. Assembling in mechanical engineering and instrument making, 2005, no 7, pp. 15-21.

2. Zhiltsov V. V., Shendaleva E. V. The new method of fuel dose pump adjustment for automatic control system of gas turbine equipment. Omsk nauchny vestnik, 2001, vol. 17, no. 1, pp. 95-97.

3. Druzhinin G. V. Nadjozhnost avtomatizirovan-nykh sistem [Reliability of computer-aided systems]. Moskow, Energia, 1977. 535 p.

4. Borovikov V. P. Prognozirovanie v sisteme STATISTICA® v srede Windows: Osnovi teorii i intensivnaja praktika na kompjutere [Forecasting at

STATISTICA® system in Windows medium. Fundamentals of theory and intensive practice on computer]. Moskow, Finansy i statistika, 2006. 384 p.

5. Dedehs V. T., German V. M., Avgustovich V. G., Arhipov G. N., Bereznjakov S. V., Pipekin V. I., Rakitin M. M., Ryzhov I. D., Sakhautdinov V. N., Smolko V. V., Tchervonjuk V. V. Metody optimizatsii ispytanij I modelirovanija system upravlenija gazotur-binnymi dvigateljami [The trial and simulation optimization methods of gas turbine engine control systems]. Moskow, Mashinostroenie, 1990. 160 p.

6. Boguslaev A. V., Olejnik Al A., Olejnik An A., Pavlenko D. V., Subbotin S. A. Progressivnye tekhno-logii modelirovanija, optimizatsii i intellektualnoj avto-matizatsii etapov zhiznennogo tsikla aviatsionnikh dvigatelej [The simulation, optimization and intellectual automation progressive technologies at aeroengine life cycle stage]. Zaporozhje, 2009. 468 p.

7. Karlin Samuel. Mathematical methods and theory in games, programming, and economics. Stanford University, Pergamon Press, London - Paris, 1959. 840 p.

8. Hudson Derek J. Statistics: Lectures on Elementary Statistics and Probability. Geneva, 1964. 296 p.

9. Karlin Samuel, Studden William J. Tchebycheff systems: with applications in analysis and statistics. Interscience publishers A. Division of John Wi-ley&Sons, New York, London, Sydney, 1969. 568 p.

Шендалева Елена Владимировна (Омск, Россия) - кандидат технических наук, доцент (Россия), доцент кафедры «Нефтегазовое дело» ФГБОУ ВО «Омский государственный технический университет» (644050, г. Омск, пр. Мира, 11, e-mail: shendalevaev@yandex.ru).

Elena V. Shendaleva (Omsk, Russian Federation) - Ph. D. in Technical Sciences, Ass. Professor, Department of Oil&Gas Busyness, Omsk State Technical University (644050, Mira prospect, 11, Omsk, Russian Federation, e-mail: shendalevaev@yandex.ru).

i Надоели баннеры? Вы всегда можете отключить рекламу.