СИСТЕМНЫЙ АНАЛИЗ, УПРАВЛЕНИЕ И ОБРАБОТКА ИНФОРМАЦИИ
УДК 629.735.45:681.5
С. Г. Микрюков, В. Н. Ефанов, А. А. Зайцева
ИНТЕГРИРОВАННАЯ СИСТЕМА УПРАВЛЕНИЯ МНОГОРЕЖИМНОЙ СИЛОВОЙ УСТАНОВКОЙ ВЕРТОЛЕТА
Рассматривается задача синтеза системы управления силовой установкой вертолета, включающей два турбовальных двигателя со свободными турбинами, которые через общую механическую трансмиссию вращают соосные винты вертолета. Обсуждается проблема согласования динамических характеристик отдельных контуров управления при изменении законов регулирования с целью уравновешивания нагрузки каждого из двух двигателей. Учитывая, что разбаланс мощностей происходит в неизвестные заранее моменты времени, предлагается метод синтеза, обеспечивающий желаемый вид нечеткого ожидания выходных переменных силовой установки. Синтез; силовая установка; вертолет; турбовальный двигатель
ВВЕДЕНИЕ
Постоянно расширяющаяся сфера применения вертолетов требует разработки летательных аппаратов с высокими летными характеристиками и технико-эксплуатационными показателями, обеспечивающими требуемый уровень безопасности полетов. Важнейшая роль в решении этой задачи отводится силовой установке вертолета, которая, как правило, строится по двухдвигательной схеме. При такой схеме два турбовальных двигателя со свободными турбинами нагружены на общую механическую трансмиссию, которая приводит во вращение несущие винты вертолета. На систему управления силовой установкой возлагается задача поддержания постоянной частоты вращения винтов во всех рабочих режимах. Поскольку силовые турбины обоих двигателей механически связаны с главным редуктором, который приводит во вращение несущие винты, то требование постоянной частоты вращения винтов оказывается эквивалентным обеспечению постоянных частот вращения силовых турбин. Однако при работающих двух двигателях распределение нагрузки между ними является трудновыполнимой задачей, поскольку любое минимальное различие между двумя системами управления может стать причиной значительной разницы в мощности, отдаваемой каждым из двигателей. В результате возникает так называемый интегративный эффект, суть которого заключается в том, что объединение эффективных подсистем не обеспечивает эффективности всей системы в целом. Отмеченное снижение эффективности совместной работы локальных подсистем в составе интегрированной системы обусловлено
Контактная информация: 8(347)273-78-36
тем, что стремление отдельных подсистем управления оптимизировать собственные показатели качества функционирования препятствуют достижению того же в других подсистемах.
Чтобы устранить подобное несогласованное взаимодействие локальных систем управления каждым из двигателей в составе силовой установки, в настоящее время используется способ, основанный на введении специальных перекрестных связей между каналами управления с целью компенсации их взаимного влияния. Так, с целью перераспределения нагрузки между двумя двигателями при помощи цифрового тракта устанавливается связь, обеспечивающая согласование их характеристик. При этом недогруженный двигатель определяется по величине частоты вращения ротора турбокомпрессора и электронный регулятор повышает частоту вращения этого двигателя до достижения равновесия.
Однако в силу сложности структурной организации системы, такой подход характеризуется приближенным учетом взаимодействий и предполагает экспериментальную настройку структуры и параметров системы управления при стендовых испытаниях двигателя. Это обстоятельство является причиной затягивания этапа технической доводки систем управления и дополнительного расходования материальных ресурсов. Кроме того, использование локальных критериев согласования, ориентированных на обеспечение совместной работы различных агрегатов силовой установки, оставляет открытым вопрос о схеме компромисса, определяющей характер предпочтения показателей функционирования одних подсистем перед другими, а также о допустимой величине уступок, уста-
навливающих пределы снижения показателей качества функционирования подсистем в составе интегрированной системы по сравнению с их величинами при автономном функционировании подсистем. Очевидно, что зависимость уровня согласования подсистем от величин такого рода уступок не является монотонной. После достижения определенного порогового значения дальнейшее снижение качества выполнения собственных функций какой-либо из локальных подсистем приводит не только к ухудшению функционирования силовой установки в целом, но и может вызвать недопустимые режимы работы последней.
Обобщая сказанное, можно сделать вывод
о том, что функциональная разрозненность локальных критериев согласования при отсутствии объединяющего их единого общесистемного показателя не позволяет достичь полного согласования характеристик взаимодействия всего комплекса взаимосвязанных агрегатов силовой установки и, тем самым, не обеспечивает высокого качества функционирования всей системы в целом.
Следовательно, при построении интегрированных систем управления силовой установкой вертолета возникает задача согласования отдельных подсистем исходя из глобальной цели функционирования всей системы при одновременном поддержании в заданных пределах требуемых значений основных газодинамических параметров. В связи с этим в данной работе предлагается метод синтеза подобных систем с использованием обобщенных показателей, позволяющих судить о способности силовой установки выполнять свои главные функции во всем диапазоне изменения условий полета и режимов работы.
1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ СИЛОВОЙ УСТАНОВКИ ВЕРТОЛЕТА
Рассмотрим математическую модель силовой установки вертолета в составе двух турбо-вальных двигателей со свободными турбинами, нагруженных на общую механическую трансмиссию, которая приводит во вращение несущие винты вертолета. В качестве силовых агрегатов рассмотрим турбовальные двигатели типа АИЫШ 2К компании ТИИБОМБСА.
Формирование математической модели проведем в три этапа. Вначале рассмотрим универсальную поэлементную нелинейную модель двигателя, отражающую физическое взаимодействие отдельных агрегатов и внешних уст-
ройств. Далее осуществим линеаризацию этой модели, чтобы получить математическое описание, пригодное для решения задач синтеза. На третьем этапе дополним полученную модель уравнениями, описывающими управляющую часть системы.
Силовая установка служит для преобразования энергии, заключенной в топливе и воздухе, в механическую энергию вращения винтов. В ее состав входят следующие функциональные элементы:
• турбокомпрессор, включающий воздухозаборник кольцевого типа, одноступенчатый центробежный компрессор, кольцевую проти-воточную камеру сгорания, одноступенчатую осевую турбину;
• силовая турбина, представляющая собой одноступенчатую осевую турбину с соосным валом трансмиссии, заключенным в вал турбокомпрессора;
• редуктор совместно с трансмиссией и основным силовым приводом;
• выпускной диффузор.
Основные физические процессы, сопровождающие работу двигателя, заключаются в следующем.
Впуск воздуха через воздухозаборник. При этом на входе в компрессор устанавливаются следующие параметры воздушного потока: дав* лут*
ление рвх и температура твх .
Сжатие воздуха в центробежном компрессоре. Полагаем, что характеристики компрессора заданы
* * ( /—! \ рк = рк \ПГкпр; ^впр р (1)
кК =кК («гкпр; ^впр}
*
где Рк - степень повышения давления в компрессоре, К - эффективный кпд компрессора, пТкпр - приведенная частота вращения турбокомпрессора, Gвпр - приведенный расход воздуха через компрессор.
Тогда при известной частоте вращения турбокомпрессора можно определить расход воздуха через компрессор и температуру воздуха
77*
к , в результате чего можно рассчитать основной параметр работы компрессора - его мощность Ьк.
к-1
1 +
(< Г
*
Лк
(2)
к -1
Свбк я(г;- т*).
Здесь Я - универсальная газовая постоянная, к - показатель адиабаты, коэффициент 5к учитывает уменьшение работы компрессора из-за перепуска воздуха.
Сгорание топливовоздушной смеси в кольцевой противоточной камере сгорания. Параметры газовоздушного потока на выходе из камеры сгорания характеризуются выражениями:
а = а +
Рг
т *
°кс рк ;
= ОсХ + ОтсртТт + Ог Ии Лк
(3)
ОгСрг
**
Здесь, соответственно, иг, рг, Тг - весовой расход, давление и температура газовоздушного потока, в свою очередь, Gт, Тт, Ни - весовой расход, температура и теплотворная способность топлива, срг, срв, срт - удельные теплоемкости газа, воздуха и топлива, пкс - коэффициент полноты сгорания топлива, окс - коэффициент восстановления давления в камере сгорания, который показывает, на сколько изменилось давление газовоздушного потока на ее выходе по сравнению с давлением воздуха за компрессо-
*
ром рк.
Расширение газа на одноступенчатой осевой турбине. Основными параметрами рабочего
процесса на турбине являются: р*Т - степень
понижения полного давления и кТ - КПД турбины. КПД турбины учитывает следующие виды потерь энергии:
• с выходной скоростью, т. е. потери кинетической энергии газа, неиспользованной в турбине;
• потери на преодоление сопротивления в сопловом аппарате и рабочем венце;
• с утечкой газа через зазоры;
• на преодоление трения в подшипниках и т. д .
Зависимости указанных параметров от режима работы турбины определяют ее характеристики, которые обычно задаются в следующем виде:
* * I * \
ЛТ = ЛТ (пТк пр ; РТ )
Огпр Огпр (пТкпр ; РТ )
(4)
''гпр '“'гпр\* Ткпр? Т,
С учетом приведенных соотношений параметры газового потока на выходе из турбины удовлетворяют следующим выражениям:
Рт
1 *
—рг; Рт
*
1 -Лт
( Л X
кг-1 Л
1-
*
V Рт )
(5)
При этом эффективная мощность турбины, потребляемая компрессором и другими агрегатами двигателя, описывается следующим соотношением:
к
кг -
1 Ог Лг (Тг*- Т/).
(6)
Объединяет все процессы, протекающие в турбокомпрессоре двигателя, уравнение движения ротора турбокомпрессора:
ёпТ
1
Тк
& (р/30)2 п
(7)
Тк
где JTk - полярный момент инерции вращающихся масс ротора турбокомпрессора.
Расширение газа в одноступенчатой силовой турбине происходит по закономерностям, аналогичным процессам в турбине турбокомпрессора:
* 1 * рСТ = рТ ;
л
СТ
ГТ1* _ ГТ1*
ТСТ = ТТ
1 -л,
кг -1 Л
Л к г
(8)
Для вычисления этих параметров необходимо знать характеристики силовой турбины
* * *
ЛСТ = ЛСТ \ПСТ пр ; РСТ Ь ) (пСТ пр ; РСТ )
Ог пр = Ог пр (ПСТ пр ; ' '-СТ )
(9)
Вращение основного силового привода через трансмиссию подчиняется следующему уравнению динамики:
СТ
1
-Ь.
агр
СТ
& (л/30)2
(10)
СТ
1
к
к
*
1
1
*
Здесь Ьст — г аЯг (тТ ТСТ )ЛсТ мех
кг - 1
мощность свободной турбины, где ПСТ мех - механический КПД силовой турбины;
Твинт — СвинтРв«в3инт В - суммарная мощность соосных винт0в, где СвИНТ — СвИНТ (фвинт;1 винт ) -
коэффициент мощности винтов, который зави-
сит от 1винт — -
V
поступи винтов и фв
угла установки лопастей винтов; далее, пвинт -частота вращения винтов, В - диаметр винтов, Уп - скорость полета вертолета; Тагр - мощность, идущая на привод агрегатов; птр - КПД трансмиссии.
После силовой турбины отработанные газы попадают в выпускную трубу выпускного диффузора. В выпускном диффузоре и выпускной трубе продолжается расширение газов, сопровождаемое падением их давления и скорости.
Линеаризация приведенных выше уравнений позволяет получить следующую модель, в которой фигурируют относительные отклонения соответствующих переменных.
Для модуля турбокомпрессора.
ёх„
■ — аихп + а!,2хг + а!,зх1к ;
хьТ — а2,1хаг + а2,2хт* + а2,3Хт* ; ХЬК — а3,\хОв +а3,2хтК* ;
х„. — а41 х„,, + а4 2х . + а4 т.х *;
хТ* — а5,\хТк* + а5,2хаг + а5,3хат ;
хбг — ав,1хов + ав,2хат ; С11)
хкТ — а7,1хлТ + а7,2хптк ; х . — а81хГт + а8 2 хп ;
РТ 8,1 8,2 пТк ’
х * ^ а чх * + а ^х *;
Тк 9,1 Кк 9,2 Рк
хкТ — ^0,1-%, + а10,2хпТк ; хРТ — а11,1хОв +а11,2 хпТк .
Аналогично для модуля силовой турбины справедлива следующая совокупность линеаризованных уравнений.
ёх„
пСТ
-------— а121хп + а122 хт + а123хт ;
, 12,1 п СТ 12,2 ТСТ 12,3 твинт
хТСТ — а13,1хСг +а13,2хТТ>+а13,3хТ|*Т ; (12)
хТ —а14 1хп + а14 2 хю ;
твинт 14,1 пвинт 14,2 Твинт
Исключая из (11) и (12) промежуточные переменные, получаем следующую компактную систему уравнений состояния для объекта управления
хпп — р1,1хпТк +01,2хат;
Тк т (13)
хпСТ — 02,1хпТк +02,2хпСТ +02,3%т .
2. АНАЛИЗ ДИНАМИКИ СИСТЕМЫ ПРИ СЛУЧАЙНОМ ИЗМЕНЕНИИ
СТРУКТУРЫ УПРАВЛЯЮЩЕЙ ЧАСТИ
Штатная система управления, поставляемая с двигателем АКШИ8 2К, обеспечивает выполнение следующих функций: пуск, регулировка частоты оборотов, ручное управление, защита от заброса вращения, средства контроля и технического обслуживания. В состав системы входит один блок цифрового управления на каждый двигатель с перекрестным контролем для системы защиты от превышения частоты вращения. Регулирование частоты вращения силовой турбины позволяет поддерживать практически постоянное значение частоты вращения винтов вертолета во всех условиях эксплуатации. При этом положение рычага общего шага винта, соответствующее потребляемой мощности, определяет базовое значение пСТ, которое задает режим работы турбокомпрессора, позволяющий уравновесить мощность, отдаваемую двигателем, с мощностью, потребляемой вертолетом. Контроллер частоты вращения силовой турбины, который представляет собой пропорционально-интегральное звено, измеряет разницу между базовым и измеренным пСТ и рассчитывает значение расхода топлива. Таким образом, пСТ, а следовательно, и пвинт сохраняют постоянное значение без статического зависания.
В тех случаях, когда происходит разбаланс мощностей, отдаваемых каждым из двигателей, осуществляется автоматическое переключение на регулятор частоты вращения ротора турбокомпрессора. Контроллер пТк также является регулятором пропорционально-интегрального типа. Регулирование частоты вращения ротора турбокомпрессора продолжается до полного распределения нагрузки, после чего осуществляется обратное переключение на регулятор пСТ.
Переход на регулирование частоты вращения турбокомпрессора происходит также в случае необходимости ограничения величины пТк с целью защиты главного редуктора от превы-
хп— а15,1хп
СТ
винт
Тк
шения крутящего момента или двигателя от перегрева.
Таким образом, интегрированная система управления многорежимной силовой установкой вертолета должна обеспечивать устойчивость и качество управления не только на всех возможных режимах работы, но и при разнообразных переходах с одного режима на другой. Решение этой задачи усложняется тем, что в реальных условиях существует дефицит достоверной информации об условиях и ожидаемых моментах переключения, а также о возможных длительностях каждого из режимов. При этом использование априорных оценок о том или ином характере случайного процесса переключения может привести к существенному несоответствию расчетных показателей реальным данным, полученным по результатам эксплуатации. В результате создаются предпосылки для автоколебательного режима переключений с одного режима на другой и, как следствие, для потери устойчивости процесса управления, приводящей к развалу функционирования системы. В этом случае приходится использовать специальные процедуры обработки эмпирических данных, отличные от статистических методов, что позволяет описывать поведение системы с помощью более адекватных характеристик, например, функций принадлежности к возможным состояниям системы [1].
При осуществлении процедуры совмещенного синтеза алгоритмов управления для множества возможных состояний системы необходимо описать поведение всего ансамбля ее выходных реакций, конкретная реализация каждой из которых зависит от моментов возможных переключений, одной обобщенной характеристикой. В качестве аналога подобной характеристики может быть использована введенная в [2], так называемая, «средняя реакция», которая получается посредством осреднения вектора выходных реакций системы управления по всем возможным моментам изменения структуры управляющей части. Специфика использования нечетких характеристик, определяющих принадлежность системы к тому или иному режиму, диктует необходимость применения обобщенной характеристики аналогичной природы. В качестве подобной характеристики предлагается использовать нечеткое ожидание выходной реакций системы [3].
Учитывая вышеизложенное, сформулируем задачу построения математической модели сис-
темы с учетом случайного изменения ее структуры.
Для этого дополним модель объекта управления (13) уравнениями, описывающими управляющую часть системы, которая в нашем случае будет содержать:
• контроллер частоты вращения силовой турбины пропорционально-интегрального типа
ёх,
= КI ГёХ^ + х„
(14)
• дополнительное корректирующее звено в цепи обратной связи изодромного типа
ёхос ёхп
*2“ЇГ + Хос =*1—7
-+ Х
Пет
(15)
Здесь хс = хп - хос - сигнал ошибки
є п СТзад ос
управления.
Объединяя (13) с (14), (15), получаем модель системы для режима управления частотой вращения силовой турбины
Р {О )хпст Ь )= во О Кстзад (' ) . (16)
Здесь Р0(О) и в0(О) - дифференциальные операторы вида
ё о ё о ё о ё
Ро (о) = ТГ + аз тг + а2 ~ГГ + а1 ~Т + "°
dt4 3 dt3 ^ dt2 1 dt '
& (D) = £ + ь; £ + ь» А + ь».
dt dt dt
Аналогично при управлении частотой вращения ротора турбокомпрессора управляющая часть будет содержать:
• контроллер частоты вращения ротора турбокомпрессора пропорционально-интегрального типа:
^=4^4 (17)
• корректирующее звено в цепи обратной связи изодромного типа
ёхпс
*4—5 + хос =*3-
хЕ = х - хос.
е пТкзад ос
ёх„
- + х
пТк :
(18)
Учитывая, что структура и параметры управляющей части системы управления частотой вращения ротора турбокомпрессора отличаются от системы управления частотой вращения силовой турбины, то ее модель будет иметь следующий вид:
Р1Ф)хпгт к)= 01 (Dкзад (г), (19)
Тк
&1 )= "ГГ + Ь2 ТТ + Ь1 ~Т + Ьо .
г>!т^\ Т 1 Т 1 Т 1 Т 1
где Р1 Ф)^-7Т + аз-7Г + а2~Г2 + «1^7 + ао =
Тх Тх Тх Тх
т_3
Тх3 ^ dt1 1 Тх
Если изменение режима управления произошло в некоторый случайный момент времени 01, то поведение системы при t > 01 будет описываться следующей системой уравнений:
Р1 ^ )Хпст (Х, 01 )= 01 Ф)хппзад (Х ) (20)
с начальными условиями, определяющими совпадение состояний системы в момент изменения режима управления:
х,„
(х, 01)
= Х„
t=0■
(X)
t=0■
(хпСТ (t, 01 )) =(хпСТ (t))
t=0з
(хпсг (t, 01 )) 0 =(хпсг (t))
X—01
(хпст (t, 01 ))”' = (хпст (х ))"
X— 01
Х—01
X— 01
Такой режим управления сохраняется в течение всего времени, необходимого для восстановления баланса мощностей двигателей силовой установки. Поскольку длительность этого процесса заранее неизвестна, то в случайный момент времени 02 происходит обратное переключение на исходный режим управления
Ро (D^пст (X, 01,02) = 00 (DКтзад (X) . (21)
Начальные условия для этой системы уравнений также должны обеспечивать совпадение состояний в момент обратного переключения:
хпт (t, ^ 02 )
— х„
Х—9л
(х, 01)
Х—9л
(хпст (t, 01, 02 )) — (хпст (t, 01 ))
X —02
(хпст (t, 01, 02 )) — (хпст (t, 01 ))
Х—02
(хпст (t, ^ 02 ))м' — (хпст (х, 01 ))м
Х—02
Х—02
3. СИНТЕЗ ИНТЕГРИРОВАННОМ СИСТЕМЫ УПРАВЛЕНИЯ НА ОСНОВЕ ПРИНЦИПА НЕЧЕТКОЙ СТАБИЛИЗАЦИИ
Для оценки поведения системы в условиях дефицита достоверной информации об интенсивностях возможных переключений режимов
управления воспользуемся нечетким ожиданием ее выходных реакций. Для формализации описания неопределенных моментов переключения режимов воспользуемся функциями принадлежности синтезируемой системы к описанным состояниям:
^ (01) — 1 ехр(-101),
^2(01, 01)—т ехр(-т(02- 01)).
С учетом сделанных замечаний нечеткое ожидание выходных реакций системы, описывающее последовательность изменения режимов работы, приобретает следующий вид:
Е{хпст (х)}—[ хпт (хX ехР(-101 )Т01 +
X
¥ X
[|!ехр(-1(02- 01 ))|1 ехР(- 101 )хпст ^ 01 )Т01Т
+
+
+
[ 1 ехр(-103 )т93 [ [ 1т ехр(-т(02 - 01 ))х
X ехр(- 1(01- 02 ))хпст (t, 01, 0
^иГи2-
Для сохранения значений управляемых переменных в заданных пределах, потребуем, чтобы нечеткое ожидание лежало в области, ограниченной заданными функциями х(1)(к) и х('1>(к), т. е. выполнялось следующее неравенство
х (1>(к) < е{спст (Х)}< х<2) (к). (22)
При решении системы неравенств (22), завершающей процедуру синтеза, могут быть использованы алгоритмы эвристических методов поиска решения, например, метод Хука-Дживса.
4. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
Проведенный анализ особенностей функционирования силовой установки вертолета при регулировании частот вращения силовых турбин, а также при переходе на регулирование частот вращения роторов турбокомпрессоров с целью устранения разбаланса мощностей, отдаваемых каждым из двигателей, позволил выявить основные закономерности взаимодействия локальных систем управления каждым из двигателей в составе силовой установки, что определило методологию разрабатываемой процедуры синтеза интегрированного алгоритма управления.
Предложен метод синтеза интегрированной системы управления многорежимной силовой установкой вертолета, в рамках которой обеспе-
ст
2
0 0
чивается согласование динамических характеристик отдельных контуров управления при изменении законов регулирования с целью уравновешивания нагрузки каждого из двух двигателей.
Учитывая, что разбаланс мощностей происходит в неизвестные заранее моменты времени, разработана аналитическая процедура формализации нечеткой информации в виде нечеткого ожидания выходных переменных силовой установки.
Предложен критерий синтеза, обеспечивающий желаемое поведение всего ансамбля выходных реакций, описывающих траекторию движения системы при переходе на режим выравнивания мощностей турбовальных двигателей с последующим возвратом на исходный автоматический режим регулирования частот вращения силовых турбин.
СПИСОК ЛИТЕРАТУРЫ
1. Нечеткие множества в моделях управления и искусственного интеллекта / А. В. Аверкин [и др.]; под ред. Д. А. Поспелова. М.: Наука, 1986. 312 с.
2. Скляревич А. Н. Линейные системы с возможными нарушениями. М.: Наука, 1975. 352 с.
3. Многоуровневое управление динамическими объектами / В. И. Васильев [и др.]. М.: Наука, 1987. 3о9 с.
ОБ АВТОРАХ
Микрюков Сергей Геннадьевич, управляющий директор ОАО «Кумертауское авиационное производственное предприятие». Окончил Ворошиловское высшее военное авиационное училище штурманов (1977), Военно-воздушную академию имени
Ю. А. Гагарина (1985), Военную академию Генерального штаба ВС РФ (1993). Иссл. в обл. создания авиационных комплексов на базе вертолетов фирмы «Камов».
Ефанов Владимир Николаевич, проф. каф. электроники и биомедицинских технологий. Дипл. инженер электр. техники (УАИ, 1973). Д-р техн. наук по упр-ю в техн. системах (УГАТУ, 1995). Иссл. в обл. создания интеллектуализированных комплексов бортового оборудования.
Зайцева Алина Александровна, ст. преп. каф. информационно-измерительной техники. Дипл. инженер по авиационным приборам и измерительновычислительные комплексам (УГАТУ, 2006). Иссл. в обл. создания интеллектуализированных комплексов бортового оборудования.